python gdal 抽取tif一个波段

最新推荐文章于 2023-03-28 21:13:11 发布

邰志敏

最新推荐文章于 2023-03-28 21:13:11 发布

阅读量605

点赞数 1

分类专栏： gis 文章标签： python gdal

本文链接：https://blog.csdn.net/taizhimin/article/details/126864088

版权

gis 专栏收录该内容

10 篇文章 0 订阅

订阅专栏

gdal 抽取一个波段

from osgeo import gdal
import os
basedir = '源文件目录'
for dir in os.listdir(basedir):
    dir = os.path.join(basedir, dir)
    if dir.endswith('.tif'):
        try :
            srcDs = gdal.Open(dir)
            srcBand = srcDs.GetRasterBand(1)
            driver = gdal.GetDriverByName("GTiff");
            dstDs = driver.Create('目的地址'),srcDs.RasterXSize, srcDs.RasterYSize, 1,options=['TILED=YES','COMPRESS=LZW'])
            dstDs.SetProjection(srcDs.GetProjection());
            dstDs.SetGeoTransform(srcDs.GetGeoTransform());
            dstDs.GetRasterBand(1).WriteArray(srcBand.ReadAsArray())
            del dstDs
        except Exception as e :
            print(e)