如何利用 GDAL-Python裁剪栅格数据(DEM等)

ARCGIS裁剪数据量比较大的栅格很容易卡死,故用gdal。

可批量裁剪多个shp的组合。

翻阅了一些博客,这么一个简单的功能写得很复杂。提供代码供有需要的各位参考。

from osgeo import gdal
 
input_raster = r"F:\JUNK\dem90.tif"
# or as an alternative if the input is already a gdal raster object you can use that gdal object
input_raster=gdal.Open(input_raster)
input_shape = r"F:\JUNK\gla-shp.shp" # or any other format
output_raster=r'F:\JUNK\test2.tif'   #your output raster file
 
ds = gdal.Warp(output_raster,
              input_raster,
              format = 'GTiff',
              cutlineDSName = input_shape,      # or any other file format
              cutlineWhere="FIELD = 'whatever'",# optionally you can filter your cutline (shapefile) based on attribute values
              dstNodata = -9999)              # select the no data value you like
ds=None     #do other stuff with ds object, it is your cropped dataset. in this case we only close the dataset.

 这是小编裁剪好的DEM高程数据,精度为12.5m。

目前已裁剪为全国分省市县的12.5高程dem。

数据下载链接:数据下载链接   

数据来源地址:地理科学生态网 ,网站地址www.csdn.store

评论 9
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值