使用仿射變換矩陣,將像素坐標映射到世界坐標.例如,使用

affine包. (還有其他方法可以使用簡單的數學方法.)

from affine import Affine

fname = '/path/to/raster.tif'

以下是獲得仿射變換矩陣T0的兩種方法.例如,使用GDAL / Python：

from osgeo import gdal

ds = gdal.Open(path, gdal.GA_ReadOnly)

T0 = Affine.from_gdal(*ds.GetGeoTransform())

ds = None # close

import rasterio

with rasterio.open(fname, 'r') as r:

T0 = r.affine

GDAL(T0)使用的變換數組的約定是引用像素角.您可能希望改為引用像素中心,因此需要將其翻譯50％：

T1 = T0 * Affine.translation(0.5, 0.5)

現在要從像素坐標轉換為世界坐標,將坐標與矩陣相乘,這可以通過一個簡單的函數來完成：

rc2xy = lambda r, c: (c, r) * T1

現在,獲取第一行,第二列(index [0,1])中柵格的坐標：

print(rc2xy(0, 1))

另外,請注意,如果需要從世界坐標獲取像素坐標,可以使用反轉仿射變換矩陣~T0.

總結

以上是生活随笔為你收集整理的geotif 添加坐标_python – 如何获取geotif中单元格的坐标？的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。