gdal影像镶嵌并裁剪

陨星落云

已于 2022-04-05 23:03:18 修改

阅读量1.5k

点赞数 4

分类专栏： gdal（python，C++）文章标签： GDAL

于 2022-04-05 23:01:12 首次发布

本文链接：https://blog.csdn.net/qq_28368377/article/details/123979641

版权

gdal（python，C++）专栏收录该内容

38 篇文章

订阅专栏

# -*- coding: utf-8 -*-
"""
Created on Tue Apr  5 20:18:49 2022

@author: 陨星落云
"""
import os 
from osgeo import gdal

def mergeAndClip(img_dir,clip_shp,out_dir):
    files = [os.path.join(img_dir,s)for s in sorted(os.listdir(img_dir))]
    merged_vrt = os.path.join(img_dir,"merged.vrt")
    clip_merged = os.path.join(img_dir,"clip_merged.vrt")
    
    gdal.BuildVRT(merged_vrt,files)
    #translate_options = gdal.TranslateOptions(format = 'GTiff',creationOptions = ['TFW=YES', 'COMPRESS=LZW'])
    
    gdal.Warp(clip_merged,merged_vrt,cutlineDSName=clip_shp,cropToCutline=True,dstNodata=-9999)
    gdal.Translate(out_dir,clip_merged)
    os.remove(merged_vrt)
    #os.remove(clip_merged)

if __name__ == '__main__': 
    img_dir = "D:\\Lake\\2017\\tiles" # 待镶嵌的影像
    clip_shp = "D:\\Lake\\2017\\clip.shp" # 裁剪范围
    out_dir = "D:\\Lake\\2017\\cliped_img.tif" # 裁剪的影像
    mergeAndClip(img_dir,clip_shp,out_dir)