Question

我downloaded的圣路易斯市包裹数据，得到的积分信息看起来像(564433.152, 1037664.554)。我认为这可能仅需除以10,000或得到学位即可，但这显然不是直接映射。我做了一些研究，认为这些可能是UTM坐标或其他东西，但是也没有成功。作为参考，圣路易斯的90.5W功率约为38.5N。有任何想法吗？代码：

import shapefile
sf = shapefile.Reader("prcl")
for i in sf.iterShapeRecords():
     parcel = i.record[0]
     points = i.shape.points
     p2 = list(map(lambda x: [y/10000 for y in x], i.shape.points))
     print(parcel,points,p2)

Answer 1

事实证明，比起数据有问题或我使用库不正确的想法，这里的答案更加有趣和深刻。我忘记了使用GDAL的ogr2ogr制作的KML。我想，好吧-显然GDAL知道如何转换这些坐标-但它是如何做到的。我认为查看JSON文件可能会更容易一些，所以我运行了ogr2ogr -f GeoJSON prcl.json prcl.dbf。正如预期的那样，这产生了一个JSON文件-但在前几行中有一个非常有趣的条目：

"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::26796" } }

事实证明，这是一个仅适用于美国一小部分地区的特殊坐标系，您可以read more about it here。不幸的是，我还没有一种很好的编程方式来转换系统，但是很有趣的是知道这些数字是从哪里来的。

与ESRI来源相距甚远

1 个答案: