使用GDAL GetGeoTransform()读取TIFF范围时如何保留高精度小数位？

不管是GDAL还是rasterio，返回的地理坐标/边界值本身都是双精度浮点数（double类型），不存在精度丢失的情况——你看到的“丢失最后两位小数”只是默认打印输出时的截断显示，而非数值本身精度受损。以下是具体的处理方案：

使用GDAL的解决方案 #

1. 直接获取原始的地理变换参数，这些参数本身携带完整精度
2. 需要输出或转换为字符串时，显式指定小数位数来完整展示

from osgeo import gdal

# 打开TIFF文件
ds = gdal.Open("your_tiff_file.tif")
if not ds:
    raise Exception("无法打开TIFF文件")

# 获取地理变换参数（双精度浮点数，无精度丢失）
geotransform = ds.GetGeoTransform()

# 计算完整的边界范围
min_x = geotransform[0]
max_y = geotransform[3]
max_x = min_x + geotransform[1] * ds.RasterXSize
min_y = max_y + geotransform[5] * ds.RasterYSize

# 显式格式化输出，保留10位小数（可根据需求调整位数）
print(f"完整边界范围：")
print(f"最小X: {min_x:.10f}, 最小Y: {min_y:.10f}")
print(f"最大X: {max_x:.10f}, 最大Y: {max_y:.10f}")

# 后续计算直接使用原始变量即可，它们保持完整精度

使用rasterio的解决方案 #

rasterio的.bounds返回的BoundingBox对象内部存储的也是双精度值，同样只需格式化输出即可展示完整小数位：

import rasterio

with rasterio.open("your_tiff_file.tif") as src:
    # 获取边界对象，内部数值为双精度，无精度丢失
    bounds = src.bounds

    # 显式格式化输出完整精度
    print(f"完整边界范围：")
    print(f"最小X: {bounds.left:.10f}, 最小Y: {bounds.bottom:.10f}")
    print(f"最大X: {bounds.right:.10f}, 最大Y: {bounds.top:.10f}")

    # 直接调用属性获取原始值用于后续计算，精度不受影响
    left = bounds.left

关键说明 #

QGIS能显示完整小数位，是因为它在界面中对坐标值做了全精度的格式化渲染；而GDAL、rasterio的默认打印行为是简化显示，并不代表返回的数值丢失了精度。只要你不主动截断，直接使用返回的原始浮点变量，就能保留所有精度。

内容的提问来源于stack exchange，提问作者cris