光栅化。合并删除图像

2024-05-14 01:24:58 发布

您现在位置:Python中文网/ 问答频道 /正文

我正在尝试合并陆地卫星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

这似乎更好,但也不正确。我有点困惑到底发生了什么-如何才能正确地合并所有图像?任何帮助都可以,我对拉斯特里奥是新手


Tags: tonamefromimportmodefilesopenmerge
1条回答
网友
1楼 · 发布于 2024-05-14 01:24:58

万一有人遇到这个问题,就会有很多问题

基本上,陆地卫星7 B1与陆地卫星8 B1不同,您需要将陆地卫星7 B1与陆地卫星8 B2合并

  • 陆地卫星8号传感器的比例与陆地卫星7号不同。将陆地卫星8重新缩放到0到255之间

相关问题 更多 >