图像增强参考程序.docVIP

反色function Inverse_Callback(hObject, eventdata, handles) x=imread(‘*.bmp’) y=255-x; 直方图均衡HisteQ_Callback(hObject, eventdata, handles) x= imread(‘*.bmp’) N=256; [row,col]=size(x); p=zeros(N); for i=1:row*col; p(x(i)+1)=p(x(i)+1)+1; %%%%%%%计算灰度值为i的个数，类似%%%%%p(i+1)=length(find(x==i))/(m*n);? end p=(double(p))/(row*col); for k=2:N p(k)= p(k) + p(k-1) ; end s=fix(256*p); for i=1:row*col; x(i)=s(x(i)+1)-1; end 阈值变换 x= imread(‘*.bmp’) prompt={Threshold(0~255):}; name=Input for Thresholding; numlines=1; defaultanswer={50}; anss=inputdlg(prompt,name,numlines,defaultanswer); th=str2num(anss{1}); [row,col]=size(x); for i=1:row*col; if x(i)th x(i)=256; else x(i)=0; end end 全局线性变换.. f= imread(‘*.bmp’) [row,col]=size(f); f1=double(f); prompt={the min graylevel for output image(0~255): the max graylevel for output image(0~255):}; name=Input for ContrastStretching; numlines=1; defaultanswer={0 255}; anss=inputdlg(prompt,name,numlines,defaultanswer); b=str2num(anss{1}); a=str2num(anss{2}); fminGray=min(min(f1)); fmaxGray=max(max(f1)); g=double(zeros([row,col])); for x=1:row for y=1:col g(x,y)=[(a-b)*(f1(x,y)-fminGray)/(fmaxGray-fminGray)]+b; end end 分段线性变换... f= imread(‘*.bmp’) [row,col]=size(f); f1=double(f); N=256; prompt={the x of the first point: the y of the first point: the x of the second point: the y of the second point: the maximum graylevel for output image:}; name=Input for ContrastStretching; numlines=1; defaultanswer={80 70 170 180 255}; anss=inputdlg(prompt,name,numlines,defaultanswer); x1=str2num(anss{1}); y1=str2num(anss{2}); x2=str2num(anss{3}); y2=str2num(anss{4}); gmaxGray=str2num(anss{5}); fmaxGray=max(max(f1)); for x=1:row for y=1:col if f1(x,y)=x1 g(x,y)=x2*f1(x,y)/x1; else if (f1(x,y)x1)(f1(x,y)=x2) g(x,y)=(y2-y1)*(f1(x,y)-x1)/(x2-x1)+y1; else g(x,y)=(gmaxGray-y2)*(f1(x,y)-x2)/(fmaxGray-x2);