要获取地理学家角的坐标, 请按照以下步骤操作:
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