所用輸入數據為遙感影像數據,可以是單波段也可以是多波段數據。
運行環境:windows10 pycharm python3.7.7 GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。
完整代碼如下:
采用批處理方式,輸入和輸出為文件路徑。只需要修改輸入輸出文件路徑即可計算所需結果。
"""
此段代碼將遙感圖像轉歸一化處理
"""
import numpy as np
from osgeo import gdal
import os
#將遙感影像歸一化處理 寫成函數
def img_normalized(inputpath,output_filepath):
gdal.UseExceptions()
ds = gdal.Open(inputpath)
band01 = ds.GetRasterBand(1)
im_width, im_height=band01.XSize, band01.YSize
dim_z = ds.RasterCount #圖像通道數
# 獲取輸出文件完整路徑
output_image = input_image_filepath.split("/")[-1]
print("output_image:", output_image)
output_image = output_filepath + output_image
output_image