遥感图像增强教学内容.pptVIP

第三节 遥感图像增强;、遥感图像光学（RGB）原理 三色光合成加、减法;;3、彩色合成 原始图像色系：;;一、遥感图像数字增强意义;二、图像增强处理技术可以分类;（一）空间域处理 g(x,y)=EH[f(x,y)] g(x,y)处理后图像， f(x,y)为点,某面域（邻域） 1、基于像素处理 逼近点 2、基于模板处理 对小的子图像（模板）;空域变换增强： g(x,y)=EH[f(x,y)] g(x,y) 值由(x,y)点的f(x,y)值决定， EH为灰度变换，则 t= EH(s) 具体方法： ;1、直接对每个像素进行灰度调整操作 a.增强对比度 ;EH如图;;;;增强对比度代码：;b图像求反;X1=imread(2zong.jpg); figure,imshow(X1) f1=200;%f1和 g1分别为f，g的最大值 g1=256; figure,plot([0,f1],[g1,0]) axis tight,xlabel(f),ylabel(g) title(intensity transformation) k=g1/f1; [m,n]=size(X1); X2=double(X1); for i=1:m %转换矩阵的每个元素 for j=1:n f=X2(i,j); g(i,j)=0; if (f=0)(f=f1) g(i,j)=g1-k*f; else g(i,j)=0; end end end figure,imshow(mat2gray(g)) ;c动态范围压缩 ;X1=imread(2zong.jpg); figure,imshow(X1) %绘制变换曲线， c=255/log(256); x=0:1:255; y=c*log(1+x); figure,plot(x,y) axis tight,xlabel(f),ylabel(g) title(intensity transformation) %转换矩阵中的每个元素 [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n g(i,j)=c*log(X2(i,j)+1); end end figure,imshow(mat2gray(g));d二值化处理：;EH如图;EH如图; figure,plot([0,f1,f1,255],[0,0,255,255]) axis tight,xlabel(f),ylabel(g) title(intensity transformation) % [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n f=X2(i,j); g(i,j)=0; if (f=0)(f=f1) g(i,j)=0; elseif (f=f1)(f=255) g(i,j)=255; end end end figure,imshow(mat2gray(g));E直方图均衡化模式;;I = imread(tire.tif); J = histeq(I); imshow(I) figure, imshow(J) figure,imhist(I,64) figure,imhist(J,64);F图像间的代数运算模式;%用原图像与压缩后的图像相减 I = imread(1zong.jpg); I = double(I)/255; %??算离散变换矩阵，返回结果为双精度型 T = dctmtx(8); %实现图像的显示块操作 B = blkproc(I,[8 8],P1*x*P2,T,T); mask =[1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; B