matlab遥感代码,遥感数字图像处理（内含Matlab)代码

【实例简介】

1 读取BIP 、BIL、 BSQ文件

2 均值滤波 中值滤波

3 边缘信息提取

4 DFT FFT

5 主成分变换

6 缨帽变换

7 图像分类(K—均值分类、最小距离分类、最大似然分类)

8 大气校正 反射率 地表温度的反演

9 Habib教授课程总结

N=6

for 1=I N

img=fopen (D: MATLAB Sample BIP,rb,)

bo=fread(img, i-1,uint8);

b=fread(img, lines*samples,'uinL8, (N-1)

band cov=reshape(b, samples, lines)

band cov 2=band Cov: 0/0x0A

band uint8=uint&(band cov2);

tif - imadjust(band uint8)

mkdir(E: MATLAB, tifbands')

name=[E: MATLAB \tifbands tif, int2str(i), tif

imwrite(tif, name, tif ); %imwrite(A, filename, fmt)

lL='波段

subplot(3, 2, 1), imshow ( tif); title(tilt)

fclose(img);

3.读取BIL的源代码

clear all

clc

lines =400

samples=640

N=6:

for i=l:N

bi-zeros(lines, samples)

for j=l: samples

img=fopen(D: MATLAB\Sample BIL, 'rb)

bb=fread(img, (i-1)*640, uint8)

bo=fread(img, IG-1),uint8);

bandi linej-fread(img, lines, uint8, 1(N*Samples-1)

fclose(img)

i( j)=bar

band uint8=uint(bi)

tif-imadjust(band uint8)

mkdir(D: MATLAB, tifbands')

name=l'D: MATLAB\tifbands tif, int2str(),tit

imwrite til,name, til

tit[波段,int2str()

subplot(3, 2, 1), imshow(tif) title(tilt

end

四、运行结果

波段

淀段4

波段

派段

图1:读取文件的六个波段图

实习2均值/中值滤波、边缘信息提取

、实验目的与原理

各种图像滤波算子可以实现图像的增强,去噪,边缘提取等。图像增强的目

的在于:1.釆用一系列技术改善图像的视觉效果,提高图像的清晰度,2.将图像