佛山网站优化包年,iis 网站制作,thinkphp建站网址,大同推广型网站建设使用仿射变换矩阵,将像素坐标映射到世界坐标.例如,使用affine包. (还有其他方法可以使用简单的数学方法.)from affine import Affinefname /path/to/raster.tif以下是获得仿射变换矩阵T0的两种方法.例如,使用GDAL / Python#xff1a;from osgeo import gdalds gdal.Open(pa…使用仿射变换矩阵,将像素坐标映射到世界坐标.例如,使用affine包. (还有其他方法可以使用简单的数学方法.)from affine import Affinefname /path/to/raster.tif以下是获得仿射变换矩阵T0的两种方法.例如,使用GDAL / Pythonfrom osgeo import gdalds gdal.Open(path, gdal.GA_ReadOnly)T0 Affine.from_gdal(*ds.GetGeoTransform())ds None # closeimport rasteriowith rasterio.open(fname, r) as r:T0 r.affineGDAL(T0)使用的变换数组的约定是引用像素角.您可能希望改为引用像素中心,因此需要将其翻译50T1 T0 * Affine.translation(0.5, 0.5)现在要从像素坐标转换为世界坐标,将坐标与矩阵相乘,这可以通过一个简单的函数来完成rc2xy lambda r, c: (c, r) * T1现在,获取第一行,第二列(index [0,1])中栅格的坐标print(rc2xy(0, 1))另外,请注意,如果需要从世界坐标获取像素坐标,可以使用反转仿射变换矩阵~T0.
