使用GDAL剪切数字高程模型

作者: v-star*위위
发布时间: 2024-12-31 06:03:06 (22天前)
转自：

2 条回复

0#
回复此人
小狮子 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 您可以使用gdal.Warp（）和shapefile作为切割线以更简单的方式实现这一点 </p> <pre> <code> from osgeo import gdal input_raster = "path/to/yourDEM.tif" # or as an alternative if the input is already a gdal raster object you can use that gdal object input_raster=gdal.Open("path/to/yourDEM.tif") input_shape = "path/to/yourShapefile.shp" # or any other format output_raster="path/to/outputDEM.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. </code> </pre> </DIV>

编辑

登录后才能参与评论