python gdal 裁剪_用Python实现矢量对栅格数据的裁剪

原标题：用Python实现矢量对栅格数据的裁剪

欢迎关注天善智能 hellobi.com，我们是专注于商业智能BI，大数据，数据分析领域的垂直社区，学习、问答、求职，一站式搞定！

对商业智能BI、大数据分析挖掘、机器学习，python，R等数据领域感兴趣的同学加微信：tstoutiao，邀请你进入数据爱好者交流群，数据爱好者们都在这儿。

矢量对栅格数据的裁剪在GIS软件中是基本功能，本文使用Python来实现该功能。其中，矢量数据是面(Polygon)类型，且矢量数据和栅格数据的坐标系一致。在这里，用到的矢量数据为geojson格式，栅格数据为tif格式。 数据读取 import geopandas as gpd from osgeo import gdal, gdalnumeric ht = gdal.Open(raster

_file_path) data

_df = gpd.read_file(vector

_file_path) poly = data_df.ix[

0][

'geometry']

接下来，将根据poly这一面状要素，对栅格数据中相应范围进行裁剪。 地理坐标转像素坐标

利用GDAL中的GeoTransform参数来实现。

defgeoToPixel(record,GeoTransform):

""" http://www.gdal.org/gdal_tutorial.html GeoTransform[0] /* top left x */ GeoTransform[1] /* w-e pixel resolution */ GeoTra