使用MATLAB读取landsat影像文件夹，并将其合成为指定波段的文件，进行显示以及保存（学习记录）

使用场景：

在学习使用matlab对遥感图像的处理后，特此记录，欢迎交流补充。

使用MATLAB读取landsat影像文件夹，并将其合成为指定波段的文件，进行显示以及图片保存。

可用波段组合（landsat8）

除了常见的 RGB 三波段组合外，还有一些其他常用的三波段组合，如下：

5-4-3： 常用于植被研究，绿度指数（NDVI）等
7-6-4： 常用于矿物识别，土地利用/覆盖分类等
6-5-4： 常用于城市区域研究，水体和土地分类等
6-7-2： 常用于植被和岩石识别，土地利用/覆盖分类等
7-5-3： 常用于植被研究，水体和土地分类等
7-6-5： 常用于火山和矿物研究等

不同的三波段组合在不同的研究领域和应用中都有其特定的优势。

---

代码

提示：需要自行更改图像地址

%读取波段
band1 = imread('E:/landsat影像/LC08_L1TP_123039_20210806_20210811_01_T1_B3.TIF');
band2 = imread('E:/landsat影像/LC08_L1TP_123039_20210806_20210811_01_T1_B4.TIF');
band3 = imread('E:/landsat影像/LC08_L1TP_123039_20210806_20210811_01_T1_B5.TIF');

% 将近红、红、绿波段数据组合成RGB图像（假彩色合成）
rgb = cat(3, band3, band2, band1);

% 使用imcrop函数裁剪图像，选择裁剪矩形的位置和大小（裁剪）
cropped_img = imcrop(rgb, [3000, 3000, 500, 500]);

% 对图像进行拉伸处理，增强对比度（增强）
stretched_img = imadjust(cropped_img, [0 1], [0.05 0.95], 1);

% 显示裁剪后的假彩色合成图像
imshow(stretched_img);

%保存整张图像
[A,R] = geotiffread('E:/landsat影像/LC08_L1TP_123039_20210806_20210811_01_T1_B3.TIF') ;
info=geotiffinfo('E:/landsat影像/LC08_L1TP_123039_20210806_20210811_01_T1_B3.TIF');
geotiffwrite('change_image.tif', rgb, R, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);