空间域滤波演示文稿.ppt

反锐化掩模(Unsharp masking?(USM)) 57/66 Matlab 函数:unsharp Demo I= imread(hestain.png); h= fspecial(gaussian,[3 3],0.5); blurred= imfilter(I,h,replicate); alpha=7; sharpened=I+alpha*(I-blurred); figure,imshow([I, sharpened]) 58/66 Matlab 函数:unsharp Demo I= imread(hestain.png); h= fspecial(unsharp); K= imfilter(I,h,replicate); imshowpair(I,K,montage) 58/66 作业2： Detector Striping：Landsat Missions 坏条带的由来：2003年5月31日，Landsat-7ETM+机载扫描行校正器(SLC) 故障，导致此后获取的影像出现了数据条带丢失，严重影响了Landsat ETM遥感影像的使用。此后Landsat 7 ETM SLC-on是指2003.5.31日Landsat 7SLC故障之前的数据产品，Landsat 7 ETMSLC-OFF则是故障之后的数据产品。 66/66 Google关键词： Landsat stripe noise /detector-striping 拓展阅读: 1. Gabor Filter 2. bilateral Filtering /sparis/bf_course/ 3. guided image filtering /eccv10/ 4. domain transform filter http://inf.ufrgs.br/~eslgastal/DomainTransform/ 更多有趣的滤波 :http://jiaya.me/ 60/66 I= double(imread(lady.png)); H = RF(I,30,100) - I + 128; G = imfilter(H,fspecial(gaussian,[3 3],100)); opacity = 50; Dest = (I*(100 - opacity)+(I+2*G-256)*opacity)/100; imshow([uint8(I) uint8(Dest)]); Ref. 1. /Imageshop/ 2. http://inf.ufrgs.br/~eslgastal/DomainTransform/ 61/66 62/66 I= double(imread(lady1.jpg)); H = RF(I,30,100) - I + 128; G = imfilter(H,fspecial(gaussian,[3 3],100)); opacity = 50; Dest = (I*(100 - opacity)+(I+2*G-256)*opacity)/100; imshow([uint8(I) uint8(Dest)]); 63/66 64/66 * * * * * * * Numerical Example 225 225 225 226 226 226 226 226 225 225 255 226 226 226 225 226 226 226 225 226 0 226 226 255 255 226 225 0 226 226 226 226 225 255 0 225 226 226 226 255 255 225 224 226 226 0 225 226 226 225 225 226 255 226 226 228 226 226 225 226 226 226 226 226 0 225 225 226 226 226 226 226 225 225 226 226 226 226 226 226 225 226 226 226 226 226 226 226 226 226 225