第六章图像变换的不变性与偏微分方程.ppt

由medk的定义可得?m,|Xm-x|k≥1∕2因为Xm是一个单调减的集合序列，同时具有有限的测度，根据实变函数中的Lebesgue极限收敛定理|Xm-x|k→|Xl-x|k并且|Xl-x|k≥1/2，再由medk的定义，有x∈medkXl即medk(Xl)?∩m﹤lmedk(Xm)■第61页，共97页，星期日，2025年，2月5日其中可测，且其中可测，且则有若第62页，共97页，星期日，2025年，2月5日定义9：图像的加权中值滤波器基于一个结构元素B，其定义为显然，medk(u)是个supinf型的，是一个平移不变的形态学算子。如果在记号上不区分图像中值滤波器和集合中值算子，都记为medk。medk作为一个单调的集合算子，可以用最大值表现公式扩展到一个图像的变换T（定理4）Tu(x)=sup{u(x)|x∈medk(X)}同时满足cl(Tu)=medk(clu)第63页，共97页，星期日，2025年，2月5日由于medk作为一个单调、平移不变的集合算子，所以T是一个单调的、对比不变的变换，根据定理7其中B={B,0∈medk(B)}。所以有：定理13：medk是图像加权中值滤波器的伴随集合算子。6.4.2中值算子的离散算法已知图像u的分辨率为m×n，u[i,j]（i=1,…m,j=1…n）表示在[i,j]处的灰度值，取结构元素为B，t为参数，V[i,j]表示运算的结果：第64页，共97页，星期日，2025年，2月5日Median(u,B,t)NB=结构元素B的像素数目FORx0=1tomFORy0=1tonFOREACH(B中的像素点p)得到p的坐标(px,py)将u[x0+px,y0+py]的灰度值存入数组array中对array进行排序，排在第（NB/2+1）的值med_array令v[x0,y0]=med_arrayENDFOREACHENDEND第65页，共97页，星期日，2025年，2月5日噪声图像及32：32：224水平线中值滤波后图像及32：32：224水平线第66页，共97页，星期日，2025年，2月5日Matlab源码：Image=imread(lena.bmp);subplot(2,3,1),imshow(Image);I=imnoise(Image,saltpepper,0.05);subplot(2,3,2),imshow(I);L_s_I=Level_Set_Line(I,32,32,256);subplot(2,3,3),imshow(uint8(L_s_I));M=medfilt2(I,[3,3]);subplot(2,3,5),imshow(M);L_s=Level_Set_Line(M,32,32,256);subplot(2,3,6),imshow(uint8(L_s));第67页，共97页，星期日，2025年，2月5日6.5欧氏不变的形态学算子6.5.1定义和微分性质记H[a]=T[x1+ax22](0)其中T[x1+ax22]表示T(u)，u=x1+ax22。如果T是单调的，那么H也是单调的。定义Th[x]=h·T[x]=h·H[0]Th[x1+ax22](0)=h·T[x1+h·ax22](0)=h·H[a·h]用D(0,M)表示圆心在0，半径为M的圆。第68页，共97页，星期日，2025年，2月5日定义5：集合算子T′是平移不变的，如果tx(T′(X))=T′(txX)定义图像变换是平移不变的，如果tx(T(u))=T(txu)定理6：（Matheron）令T′是平移不变的单调集合算子，那么存在一个集合族BB={X,0∈T′(X)