基于矢量文件的栅格数据批量裁剪与导出

John H.

已于 2022-09-29 17:18:31 修改

阅读量860

点赞数

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

于 2021-08-10 17:22:13 首次发布

本文链接：https://blog.csdn.net/weixin_42778485/article/details/119574868

版权

本文介绍了如何利用Python结合GDAL库，对矢量文件定义的区域进行栅格数据批量裁剪，并导出处理后的结果。通过示例代码，详细展示了具体的实现步骤。

摘要由CSDN通过智能技术生成

见如下代码。

def clip_raster_by_shpfile_save_into_tif_formate_jpg_formate(input_raster, input_shpfile,
                                                             output_raster, output_jpg):
    '''这个函数是用shpfile裁剪raster，同时保存两个结果到指定文件夹中
    一个是tif格式带坐标系的tif栅格数据，一个是jpg格式无坐标是的jpg图片数据
    input_raster待裁剪的栅格数据完整路径
    input_shpfile用于裁剪栅格的矢量数据完整路径
    output_raster输出栅格数据的完整路径
    output_jpg输出图片数据的完整路径
    -----------------------------------------------------------------------------------
    gdal.Warp参数说明如下：
        outputRasterNamePath裁剪后文件保存路径及文件名
        input_RASTER待裁剪影像源文件
        cutlineDSName用于裁剪栅格数据的矢量文件完整路径
        cropToCutline保证裁剪后影像大小跟矢量文件的图框大小一致
        cropToCutline设置为False时，结果图像大小会跟待裁剪影像大小一样，则会出现大量的空值区域
    '''
    shapefile_data_reader = shapefile.Reader(input_shpfile)
    # 获取该shp文件有多少条记录
    shapefile_data_record = shapefile_data_reader.numRecords
    flag_value = 0
    if shapefile_data_re