Rasterio.merge删除图像

Question

我正在尝试合并陆地卫星8波段。我的数据集有几个相邻的B3波段

import rasterio 
from rasterio.plot import show
from rasterio.merge import merge
from rasterio.warp import reproject

import os 
from matplotlib import pyplot


files_to_merge = []
for file in os.listdir():
    if file.endswith(".TIF") and ("GM" not in file):
        files_to_merge.append(rasterio.open(file))

test, out_trans = merge(files_to_merge)
show(test)

当我在上面的代码之上时，我得到：

merged image

要合并的文件列表具有以下对象：

 <open DatasetReader name='LE07_L1TP_017015_20210503_20210529_02_T1_B3.TIF' mode='r'>,
 <open DatasetReader name='LC08_L1TP_016015_20210504_20210517_02_T1_B3.TIF' mode='r'>,
 <open DatasetReader name='LE07_L1TP_017016_20210503_20210529_02_T1_B3.TIF' mode='r'>,
 <open DatasetReader name='LC08_L1TP_016014_20210504_20210517_02_T1_B3.TIF' mode='r'>,
 <open DatasetReader name='LE07_L1TP_017014_20210503_20210529_02_T1_B3.TIF' mode='r'>,

我有点困惑，为什么其他的光栅图像没有显示在合并的图像中。当我删除第一个波段时，我可以重新运行merge命令：

test, out_trans = merge([files_to_merge[1],files_to_merge[2], files_to_merge[3], files_to_merge[4]])
show(test,cmap='terrain')

merged image 2

这似乎更好，但也不正确。我有点困惑到底是怎么回事--我怎么才能正确地合并所有的图像呢？如果有任何帮助，我是第一次接触rasterio。

Answer 1

如果有人遇到这个问题，这里有一堆错误的东西

• 这不是你合并陆地卫星7号和8号图像的方式。如果你做了最少的研究(我没有)，你会发现这个网站：https://www.usgs.gov/faqs/what-are-band-designations-landsat-satellites?qt-news_science_products=0#qt-news_science_products

基本上，陆地卫星7 B1与陆地卫星8 B1不同，您希望将陆地卫星7 B1与陆地卫星8 B2合并。

• 陆地卫星8传感器的比例与陆地卫星7不同。重新调整陆地卫星8的比例以在0到255

之间移动