《数字图像处理作业(第三次)直方图处理》由会员分享，可在线阅读，更多相关《数字图像处理作业(第三次)直方图处理（7页珍藏版）》请在金锄头文库上搜索。
1、数字图像解决作业（第三次)直方图解决第一题第二题（例3.6)1.问题:1)编写函数Y=Hieq(X）, 对灰度图象X进行直方图均衡化，得到均衡化后的图像Y,假设的数据类型为unt8。2)运用自己编写的直方图均衡化函数，重现例.6中的实验，即对图.16中的图像进行直方图均衡化。2程序:(1) 函数程序:funtioY = MyHst(X)%对灰度图象X进行直方图均衡化,得到均衡化后的图像 m,n=sze(X)； Azos(1,256)； fo =0:255 A(k+)=length(in(X=k)）/(*n）； nd 1=s(1,26）； for 1:56 orj:i S1()=（j)1(i);
2、 %算Sk en end S2=round（S*56）; or i=1:256 B（i)=sum(A（find（S2i）); end fure,ar(:55，,b); ttle（均衡化后的直方图）； lbel（灰度)，yabel(浮现概率) YX; for i=0:255 Y(find(X=)=S2(i+1）; endnd(2) 主程序:=imred(g36（1)(to_eft)f);gue,imsho（X);YMyHte(X);fgure,imow(Y)；3. 实验成果:第三题(例3.12)问题:编写程序，重现例3.12中的实验，产生图.2。2程序：X=irad(F27(a)(tungstn
3、_original)tf);igure,iow（）;Y=histq(X）;fiure,sw(Y)；=4.0;00.4;1=.02；2=0.；,n=si（);X=X;M=(X1(:）/（m*n);g=um(sum（(X-g).2)/(m*n);fri=2:m-1for =2:n-1X2=X1(i-1：i+1,j-1:j+1)；M=s（X2（:）)/；Vs=sm(sum(2Ms).2);if(Msk*Mg&k1*g=sVn,qn X1=duble(X); X2=X1; fr i=1:heig-n1 or j=1:hig-n+1 C=X1(:i+(-1）,j:j(n1); %取出X1中从(,)开始的n行列元素,即模板（nn的） =C(1,:); 是矩阵的第一行 fr 2:n EE，(,:); %将C矩阵变为一种行矩阵 ed mm=mn（E); %m是中值 X(i+(n-1)2,j+(n-1）/2)=mm;%将模板各元素的中值赋给模板中心位置的元素 end end Yunt8(X2);nd(2) 主程序：=imrea(Fig0335(a）(cktboard_ltpp_probpt05).tif);igue,mhow()；n=3;Y=mid_iltr(X，n）；figure，imso(Y);3. 实验成果:被椒盐噪声污染的电路板X射线图像：用33中值滤波器降噪后的图像：
