python gdal 按照属性提取矢量，并进行影像裁切

小小毛喵

已于 2024-04-06 08:44:02 修改

阅读量115

点赞数

CC 4.0 BY-SA版权

分类专栏： python 文章标签： gdal python 矢量裁切栅格

于 2018-07-04 17:33:47 首次发布

本文链接：https://blog.csdn.net/qq_37798011/article/details/80916118

＃ -  *  -  coding：utf-8  -  *  - 
“”“ 
/ ******************* ************ 
ClipRaster.py:1 
    。按照属性提取矢量（在上层进行操作）
    2。矢量转为同源数据的栅格掩膜（某些矢量存在交叉点，无法使用gdalwarp进行裁切）
    3.栅格掩膜与源数据进行运行，得到裁切后的影像@version <1.1> 2018-07-03 Wujd：创建。“” 来自osgeo import ogr 来自osgeo import gdal 来自osgeo import gdal_array作为ga import os，math，ogr，osr ＃解决SHAPE文件的属性值










gdal.SetConfigOption('SHAPE_ENCODING', '')
# #解决中文路径乱码问题
# gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "NO")

"""
路径设置
"""
input_vector ='./data/lns_210000_city.shp'
input_raster ='./data/GF1_WFV4_E124.8_N43.3_20180524_L1A0003210738_cal_orth.tiff'
extent_vector = './result/clip/lns_210000_city_42000102_clip.shp'
extent_raster = './result/clip/image_extent.tiff'
clip_raster = './result/clip/image_42000102_clip.tiff'


def copyInLayer(input_vector,extent_vector,filterStr):

    ds = ogr.Open(input_vector)
    driver = ogr.GetDriverByName("ESRI Shapefile")

    if os.access(extent_vector, os.F_OK):
        driver.DeleteDataSource(ex