python gdal 基于栅格shp文件裁剪geotif图
import numpy as npimport gdalimport operatorfrom functools import reduceimport osimport shapefiledef geo2imagexy(dataset, x, y):‘’’根据GDAL的六 参数模型将给定的投影或地理坐标转为影像图上坐标:param dataset: GDAL地理数据:param x: 投影或地理坐标x:param y: 投影或地理坐标y:return: 影坐标或地理坐标(x,