ç´æ¥è¿è¡è¿ä¸ªæ
clear all
%ä¸ï¼å¾åçé¢å¤çï¼è¯»å
¥å½©è²å¾åå°å
¶ç°åº¦å
PS=imread('s7.jpg'); %读å
¥JPG彩è²å¾åæ件
figure(1);subplot(2,2,1);imshow(PS);title('åå¾åç°åº¦å¾');
%äºï¼ç»å¶ç´æ¹å¾
[m,n]=size(PS); %æµéå¾å尺寸åæ°
GP=zeros(1,256); %é¢å建åæ¾ç°åº¦åºç°æ¦ççåé
for k=0:255
GP(k+1)=length(find(PS==k))/(m*n); %计ç®æ¯çº§ç°åº¦åºç°çæ¦çï¼å°å
¶åå
¥GPä¸ç¸åºä½ç½®
end
figure(1);subplot(2,2,2);bar(0:255,GP,'g') %ç»å¶ç´æ¹å¾
title('åå¾åç´æ¹å¾')
xlabel('ç°åº¦å¼')
ylabel('åºç°æ¦ç')
%ä¸ï¼ç´æ¹å¾åè¡¡å
S1=zeros(1,256);
for i=1:256
for j=1:i
S1(i)=GP(j)+S1(i); %计ç®Sk
end
end
S2=round((S1*256)+0.5); %å°Skå½å°ç¸è¿çº§çç°åº¦
for i=1:256
GPeq(i)=sum(GP(find(S2==i))); %计ç®ç°ææ¯ä¸ªç°åº¦çº§åºç°çæ¦ç
end
figure(1);subplot(2,2,4);bar(0:255,GPeq,'b') %æ¾ç¤ºåè¡¡ååçç´æ¹å¾
title('åè¡¡ååçç´æ¹å¾')
xlabel('ç°åº¦å¼')
ylabel('åºç°æ¦ç')
%åï¼å¾ååè¡¡å
PA=PS;
for i=0:255
PA(find(PS==i))=S2(i+1); %å°å个åç´ å½ä¸ååçç°åº¦å¼èµç»è¿ä¸ªåç´
end
figure(1);subplot(2,2,3);imshow(PA) %æ¾ç¤ºåè¡¡ååçå¾å
title('åè¡¡ååå¾å')
imwrite(PA,'PicEqual.bmp');
温馨提示:答案为网友推荐,仅供参考