遥感数字图的滤波处理.docxVIP

实验三：遥感数字图像的滤波处理 实验目的 本次实验目的有两个：掌握ERDAS软件的滤波功能的使用和操作；掌握用编程实现图像滤波的过程。 实验要求 1. 要求利用ERDAS软件，使用成都市Landsat真彩色合成影像进行均值滤波、中值滤波处理，给出实验步骤和结果，及相应的描述。 2. 要求采用编程(Matlab或其他平台)实现对成都市Landsat真彩色合成图像进行均值滤波和中值滤波，给出实验步骤和结果，及相应的描述。给出相应的代码文件和工程文件。 实验内容 1. 基于ERDAS软件的遥感数字图像滤波，包括均值滤波和中值滤波。 2. 基于编程的遥感数字图像滤波，包括均值滤波和中值滤波。 注：编程处理数据时，由于童鞋们的电脑硬件不是太优秀，可以从成都市影像中裁剪出合适的大小来做。 实验过程 （一）基于ERDAS软件的遥感数字图像滤波 1、均值滤波 方法步骤： （1）在ERDAS菜单上选择Main/Imagin Interpreter/Spatial Enhancement/Focal Analysis命令，打开Focal Analysis对话框。 （2）在对话框中打开待处理的图像文件，设置输出文件路径并命名。 （3）Subset Dfintion为处理范围，这里不用输入数值，采用默认的数值。数据输出类型为Unsigned 8 bit。 （4）定义处理函数，由于这里是均值滤波，所以这里选择Mean。选中Ignore Zero in States复选框，表示在输出数据统计时忽略0值。 （5）单击OK按钮，执行聚焦分析操作，结果如下图。 2、中值滤波 这里的操作方法与均值滤波基本相同，只是在定义处理函数时，由于这里是中值滤波，所以这里选择Median。 执行中值滤波操作的结果如下图。 （二）基于编程的遥感数字图像滤波 1、中值滤波 代码： clc,clear,close all; img = imread(C:\Users\wqy\Desktop\456.jpg);%读取原图像 gray = rgb2gray(img); salt = imnoise(gray,salt pepper,0.08); %加入椒盐躁声 subplot(221) imshow(salt); title(salt); sort=2; [ height, width ] = size(gray); x1 = double(salt); x2 = x1; name = {salt 3x3 mid 5x5 mid 7x7 mid}; n1 = [3 5 7]; for n = n1; for i = 1: height-n+1 for j = 1:width-n+1 mb = x1( i:(i+n-1), j:(j+n-1) ); mb = mb(:); mm = median(mb); x2( i+(n-1)/2, j+(n-1)/2 ) = mm; end end img = uint8(x2); subplot(2,2,sort) imshow(img); title(name{sort}); sort = sort+1; end 结果： 2、均值滤波 clc,clear,close all; img = imread(C:\Users\wqy\Desktop\456.jpg);%读取原图像 gray = rgb2gray(img); [height,width] = size(gray); %输入图像是height x width的,且hightn,widthn x1 = double(gray); x2 = x1; name={3x3 Laplace,5x5 Laplace,7x7 Laplace,9x9 Laplace} sort=1; n1=[3 5 7 9] for n = n1; a(1:n,1:n) = 1; %n×n模板,元素全是1 for i = 1:height-n+1 for j = 1:width-n+1 c = x1(i:i+(n-1),j:j+(n-1)).*a; %取出x1中从(i,j)开始的n行n列元素与模板相乘 s = sum(sum(c)); %求c矩阵中各元素之和 x2(i+(n-1)/2,j+(n-1)/2) = s/(n*n); %将与模板运算后的各元素的均值赋给模板中心位置的元素