如何解决如何使用 GDAL VRT 在单个 GeoTiff 中将图像堆叠为光栅带?
我有一堆 .tif 图像,但我想将它们全部堆叠为 1 个 .Tif 图像。我如何使用 gdalVrt 堆叠所有 .Tif 文件??
from osgeo import gdal
from PIL import Image
import numpy as np
from numpy import asarray
import matplotlib.pyplot as plt
import os
os.listdir('../LC08_L1TP_137042_20210301_20210301_01_RT/')
band1=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B1.TIF')
band2=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B2.TIF')
band3=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B3.TIF')
band4=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B4.TIF') # near infra red
band5=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B5.TIF') # near infra red
band6=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B6.TIF')
band7=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B7.TIF')
band8=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B8.TIF')
band9=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B9.TIF')
band10=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B10.TIF')
band11=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_B11.TIF')
band12=gdal.Open('LC08_L1TP_137042_20210301_20210301_01_RT_BQA.TIF')
band1_array = band4.ReadAsArray()
band2_array = band4.ReadAsArray()
band3_array = band4.ReadAsArray()
band4_array = band4.ReadAsArray()
band5_array = band4.ReadAsArray()
band6_array = band4.ReadAsArray()
band7_array = band4.ReadAsArray()
band8_array = band5.ReadAsArray()
band9_array = band4.ReadAsArray()
band10_array = band4.ReadAsArray()
band11_array = band4.ReadAsArray()
band12_array = band4.ReadAsArray()
stack= (band1_array+band2_array+band3_array+band4_array+band5_array+band6_array+band7_array+band8_array+band9_array+band10_array+band11_array+band12_array)
print("The stacked images is \n \n",plt.imshow(stack))
我确定这不是正确的方法..
解决方法
步骤 1.) 使用选项“separate=True”创建虚拟光栅 (VRT) 以将图像堆叠为单独的波段:
[root@ip-10-0-0-218 ~]# df -h
Filesystem Size Used Avail Use% Mounted on
devtmpfs 1.7G 0 1.7G 0% /dev
tmpfs 1.7G 0 1.7G 0% /dev/shm
tmpfs 1.7G 22M 1.7G 2% /run
tmpfs 1.7G 0 1.7G 0% /sys/fs/cgroup
/dev/nvme0n1p2 16G 7.1G 8.2G 47% /
/dev/nvme0n1p1 599M 6.4M 593M 2% /boot/efi
tmpfs 345M 0 345M 0% /run/user/1000
步骤 2.) 将虚拟栅格 (VRT) 转换为 GeoTiff:
from osgeo import gdal
ImageList = ['Band1.tif','Band2.tif','Band3.tif'] # or use sorted(glob.glob('*.tif')) if input images are sortable
VRT = 'OutputImage.vrt'
gdal.BuildVRT(VRT,ImageList,separate=True,callback=gdal.TermProgress_nocb)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。