如何从翻转或 tfw 文件中获取投影(epsg 编号)

how to get projection(espg number) from fif or tfw file

我有一个 tfw 文件的 tif 文件

 30.00000000
 0.00000000
 0.00000000
-30.00000000
4558664.95812873
4002879.73390571

文件所有者说他们的数据在wgs84 中并用arcgis 生成tif。最后两个数字太大,不是wgs84。

一些 .tif 文件在文件中嵌入了 GeoTIFF 标签,这些标签指定了栅格的空间参考。要查看 .tif 是否具有内置的空间参考,请导航到 ArcMap 目录 window 中的 .tif,右键单击 .tif,选择属性,然后向下滚动到常规选项卡上的空间参考属性。或者,如果您想在 Python 代码中执行此操作,只需检查空间参考的工厂代码(也称为众所周知的 ID 或 WKID)是否为零:

import arcpy
raster = arcpy.Raster(r"C:\path\to\raster.tif")
wkid = raster.spatialReference.factoryCode
print("spatial reference " + str(wkid) if 0 != wkid else "NO spatial reference!")

请参阅 ArcGIS 的 ArcPy 模块中 Raster and SpatialReference Python 类 的文档。

正式称为 world file 的 .tfw 文件不包含空间参考。它提供地理配准信息,但不是实际的空间参考。您将无法从中获取空间参考。但是,您可以获得一些线索。例如,世界文件中的最后两个数字是左上角像素的 x 和 y 坐标。显然,这些值不是经度和纬度,因此如果世界文件正确,则数据集不太可能在 WGS84 中。