python 读取geotiff_从GeoTIFF文件中获取经度和纬度

要获取地理学家角的坐标, 请按照以下步骤操作:

from osgeo import gdal

ds = gdal.Open('path/to/file')

width = ds.RasterXSize

height = ds.RasterYSize

gt = ds.GetGeoTransform()

minx = gt[0]

miny = gt[3] + width*gt[4] + height*gt[5]

maxx = gt[0] + width*gt[1] + height*gt[2]

maxy = gt[3]

但是, 它们可能不是纬度/经度格式。正如贾斯汀指出的那样, 您的Geotiff将存储在某种坐标系中。如果您不知道坐标系是什么, 可以通过运行

gdalinfo

来找出:

gdalinfo ~/somedir/somefile.tif

输出是什么:

Driver: GTiff/GeoTIFF

Size is 512, 512

Coordinate System is:

PROJCS["NAD27 / UTM zone 11N",

GEOGCS["NAD27",

DATUM["North_American_Datum_1927",

SPHEROID["Clarke 1866",6378206.4,294.978698213901]],

PRIMEM["Greenwich",0],

UNIT["degree",0.0174532925199433]],

PROJE