栅格数据波段处理及合并
发布时间
阅读量：
阅读量
该函数通过运用RGB三波段影像，计算出RGBGDSI植被指数，并将该指数作为一个新增波段，与原始影像进行融合，从而生成包含四个波段的影像数据，以提升影像信息的丰富度。鉴于影像数据规模可能较大，一次性将单个波段加载至内存中存在导致内存溢出的风险。为此，该函数采用了分块读取与写入的方式处理数据，其中block_size参数用于设定分块的尺寸，可根据实际需求进行调整，默认情况下每次读取和写入的数据量为5000行和列。
import sys
import math
import numpy
from osgeo import gdal
from osgeo import gdal_array
def muti_bands_combine(input_raster, output_raster, block_size=5000):
'''input_raster输入影像的完整路径
output_raster输出影像的完整路径'''
in_daopen = gdal.Open(input_raster)
input_band1_inform = in_daopen.GetRasterBand(1)
全部评论 (0)
还没有任何评论哟~