需手动解析GeoKeyDirectoryTag（34735）等标签：先用TIFFGetField读取ushort数组，前4元素为版本等元信息，后续每4元素含KeyID/TagLocation/Count/ValueOffset；TagLocation为0时值在ValueOffset字段，为34736或34737时分别索引GeoDoubleParamsTag或GeoAsciiParamsTag；注意校验count≥4且数组长度被4整除，BigTIFF需新版libtiff支持。
如何用 libtiff 读取 GeoTIFF 的 GeoKeyDirectoryTag 和 GeoDoubleParamsTag
GeoTIFF 的地理信息不存于常规 EXIF 区域，而是通过 TIFF 标签嵌入：核心是 GeoKeyDirectoryTag（34735），它定义坐标系、投影参数和键值映射；辅以 GeoDoubleParamsTag（34736）存浮点参数（如标准纬线、中央经线），GeoAsciiParamsTag（34737）存字符串（如椭球名、单位名）。libtiff 不自动解析这些标签，必须手动提取并按规范解码。
实操要点：
-
TIFFGetField(tif, TIFFTAG_GEOKEYDIRECTORY, &count, &dir)获取GeoKeyDirectoryTag原始 ushort 数组；前4个元素为版本/修订/主键数/扩展键数，后续每4个元素为KeyID, TagLocation, Count, ValueOffset -
TagLocation == 0表示值直接存在第4字段（仅限 short 类型）；TagLocation == 34736指向GeoDoubleParamsTag数组索引；TagLocation == 34737指向GeoAsciiParamsTag字符串偏移 - 务必检查
count >= 4且数组长度能被4整除，否则是损坏或非标准 GeoTIFF
为什么 TIFFTAG_GEOKEYDIRECTORY 读出来全是 0？
常见原因是 TIFF 文件使用了 BigTIFF 格式，而旧版 libtiff（GeoKeyDirectoryTag 支持不完整，或编译时未启用 LIBTIFF_BIGTIFF_SUPPORT。另一个隐蔽原因是图像被压缩（如 LZW、ZIP），某些 libtiff 版本在未调用 TIFFSetDirectory(tif, dirnum) 切换到主图像目录前，无法正确读取自定义标签。
排查步骤：
立即学习“C++免费学习笔记（深入）”；
- 用
tiffinfo your.tif确认是否含GeoKeyDirectoryTag及其长度（非零） - 检查 libtiff 版本：
std::cout ，建议 ≥ 4.5.0 - 确保在
TIFFOpen()后立即调用TIFFSetDirectory(tif, 0)（即使单页图） - 若仍失败，尝试用
TIFFReadCustomDirectory()手动定位并解析 tag offset（需从 IFD 头手动跳转）
如何把 GCS_WGS_84 或 PCS_NAD83_UTM_zone_10N 转成 EPSG 代码？
GeoTIFF 规范中 GTModelTypeGeoKey（2048）、GTRasterTypeGeoKey（2049）、GeographicTypeGeoKey（2049）等键值对应预定义常量，但不直接提供 EPSG。你需要查表映射——最可靠方式是硬编码常用组合，而非依赖运行时网络查询。
关键映射逻辑：
- 若
GTModelTypeGeoKey == ModelTypeProjected且ProjectedCSTypeGeoKey == PCS_NAD83_UTM_zone_10N（= 26910），则 EPSG = 26910 - 若
GTModelTypeGeoKey == ModelTypeGeographic且GeographicTypeGeoKey == GCS_WGS_84（= 4326），则 EPSG = 4326 - 注意：EPSG 代码不是所有 GeoKey 都有对应，例如自定义投影需结合
ProjLinearUnitsGeoKey、ProjStdParallel1GeoKey等重建 WKT
可复用的常量定义见 geotiff.h（来自 libgeotiff），但 libgeotiff 是独立库，不随 libtiff 自动链接——若不想引入新依赖，直接查 libgeotiff 的 geotiff.h 复制 key 定义更轻量。
解析出的坐标系参数为何和 QGIS 显示不一致？
根本原因在于 GeoTIFF 存储的是“模型空间”参数（即栅格像素坐标 → 地理坐标变换所需参数），而非最终地理范围。你拿到的 ModelTransformationTag（33550）或 ModelPixelScaleTag（33551）+ ModelTiepointTag（33552）才是定位依据。很多开发者只读 GeoKey，却忽略这些基础几何标签，导致坐标计算偏移。
必须协同使用的标签：
-
ModelPixelScaleTag（33551）：3 元 double 数组，表示 X/Y/Z 方向每像素代表的真实单位（如度或米） -
ModelTiepointTag（33552）：至少 6 元 double 数组，格式为[i j k x y z]，其中(i,j,k)是像素坐标（通常 i=j=0,k=0），(x,y,z)是对应地理坐标 - 两者结合可构造仿射变换矩阵：
X_geo = tie_x + i * scale_x，Y_geo = tie_y - j * scale_y（注意 Y 常反转，因图像坐标原点在左上）
真正棘手的是多 Tiepoint 或非正交缩放——这时必须用 ModelTransformationTag 的 16 元素 double 矩阵做齐次变换，且要确认矩阵存储顺序（行优先还是列优先，GeoTIFF 规范要求列优先）。
