在未知时识别Shapefile的坐标系？

#1 楼

总有一种“蛮力”方法：


采用已知坐标系的图层应与未知图层重叠。
现在对未知层的投影进行一些有根据的猜测。 （UTM，Plate Carree等）。将您已知的坐标系图层投影到每个投影中，直到找到与未知图层尽可能匹配的图层为止。

#2 楼

我通过添加使用他们的技术的答案来促进Mapperz和Brandon Copeland的评论。


如果具有ArcGIS Desktop，则将一些覆盖相同区域的参考数据添加到ArcMap。此参考数据必须具有有效的坐标系（投影）定义。
将数据添加到未知坐标系中。
研究在此区域中可能使用的坐标系。检查的地方之一是EPSG注册表。
尽可能设置数据框的坐标系。如果带有已知坐标的数据与不带有已知坐标的数据对齐，则您可能已经确定了坐标系。

该技术利用了以下事实： ArcMap可以将数据转换为内存中的新坐标系。没有已知coordsys的数据无法转换，因此只能显示。通过将已知数据投影到内存中，我们可以快速查看它在各种可能的坐标系中具有什么坐标。

#3 楼

由于这个问题永远不会变老，所以我建立了一个执行蛮力方法的站点。如果将压缩的shp + shx拖到地图上，它将在PostGIS中可用的每个坐标系中进行映射。假设您知道“正确”的外观，则可以缩放到该区域，然后单击多边形以从epsg.io获取.prj文件。

#4 楼

在上一份工作中，我有一层名为“ FSU_Geol.shp”的地质数据（多边形）。我的老板把它给我，并要求我找出一些事情。首先，客户端将他交给了这个shapefile，但没有.prj文件，因此他希望我自己弄清楚。他还希望我弄清楚什么是地质学。我可以不停地前进，但是随它去吧……。我用谷歌搜索，然后我在这里结束。

命名约定的首字母缩写表示“前苏联”，而Google搜索将我“直接”带到了消息来源（USGS）。我可能需要了解的有关shapefile的所有信息都位于我点击的顶部链接。我并不是说Google可以找到任何东西，但我刚从大学毕业，就只是“暗中摸索”，然后看看我得到的反馈！

在另一个例如，办公室里有人使用批处理功能下载了许多shapefile。我现在不记得软件包的名称，但是缺少.prj文件。我只是进入ArcCatalog中的元数据，并且实际上在那里找到了答案。同样，不要引用我，但我认为那里有一个指向源的URL。

从那天起，如果元数据中没有任何线索，我就从Google搜索开始“！”

#5 楼

Esri有两个重要的链接，对此进行了详细介绍：


使用ArcMap识别未知的投影坐标系
识别空间参考，投影或数据坐标系

#6 楼

我真的很喜欢这个网络工具：http://projfinder.com/。
从未知文件中获取一些X，Y坐标，将其缩放到地球上大约相同的位置，粘贴坐标，然后猜测。

#7 楼

询问制作人。

如果您知道您所在地理区域的通用CRS，可以尝试其中的一些。但是问比较好。

#8 楼

不是一种工具（我不知道有一种工具可以让您做到这一点），但请查看@mkennedy的回复，我将如何将此点转换为WKID 4326？她解释了如何得出正确的空间参考。 SpatialReference.org和耐心将是您的朋友。

此外，Esri还提供了有关如何猜测坐标系的指南（尽管如果您对数据有更多了解，我更喜欢mkennedy的方法）。

#9 楼

对于这个问题，我找到的最佳答案是非技术性的：找出数据的来源。机构和组织倾向于与其使用预测相一致。知道它来自您的州DOT吗？查看他们的其余数据，看看它能告诉您什么。不知道它来自哪里？有根据的猜测很可能会使您走上正确的道路。

至少可以使用蛮力解决问题更可行！

#10 楼

此时严重过时，但是Werner Flacke和Birgit Klaus在2007年发布了关于ArcScripts的“查找投影”。不幸的是，我认为没有源代码。它基于VBA，因此仅可在ArcGIS Desktop v9.2和9.3中使用。
它确实包含两个shapefile，它们具有EPSG大地参数数据集中的感兴趣区域，可用于缩小可能性。

Blue Marble Geographic Calculator和Geographic Transformer具有坐标系统恢复工具。

#11 楼

请尝试使用GDAL随附的ogrinfo程序。

请参阅如何使用OGR访问Shapefile元数据？

因此，如果没有.prj文件，ogrinfo可能不会提供投影信息，但它仍然是帮助您调查潜在预测列表的有用工具。例如，

Geometry: Polygon
Feature Count: 269 
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000) 

此回答对我而言很有意义

假设您知道数据的大致来源，因此现在有了较短的可能投影列表。

其他最佳做法可能是：


首先问给您数据的人，即使他们不确切知道应该列出哪些可能的预测，可以首先查看。
看看是否有任何关联的元数据。

最好从已知的位置开始，而不是从头开始尝试蛮力。

#12 楼

要添加到此帖子的其他响应中，我将添加以下内容：

如何：使用ArcMap识别未知的投影坐标系

对NAD1​​927进行特定检查


如果坐标以十进制度为单位，例如在-180和+180之间，在-90和+90之间，则标识用于数据的地理坐标系（基准）。在9.2版之前，默认情况下，ArcMap将GCS_Assumed_Geographic_1分配给数据。这会将数据放置在NAD_1927基准上


此外


如果数据在美国，并显示坐标的范围小数点左边的是6、7或8位数字，数据可能投影到了国家平面或UTM坐标系。


常见问题解答：投影基础知识：GIS专业人士需要什么认识

本页提供了一些确定投影的检查，这对协助识别未知坐标系非常重要。


以下概念对于了解ArcGIS中地图投影的使用。请注意，尽管投影的主题非常广泛，但是本文仅涉及几个重要主题。


坐标系（也称为地图投影）是任意的空间数据的名称。它们的目的是为有关地球表面上特定位置或区域的交流提供通用基础。处理地图投影时最关键的问题是知道投影是什么，并具有与数据集相关的正确坐标系信息。
在设计第一个地图投影时，错误地认为地球是平坦的。后来对假设进行了修改，并且假设地球是一个完美的球体。在18世纪，人们开始意识到地球并不是完美的圆形。这是制图球体概念的开始。
为了更准确地表示地球表面的位置，地图制作者研究了地球的形状（大地构造），并创建了球体概念。然后设计了地理坐标系（GCS），其中包括一个基准，度量单位和一个本初子午线。基准将椭球体连接到地球表面的特定部分。最新的基准面被设计为适合整个地球的表面。
北美最常用的基准面是：

•使用Clarke 1866椭球体的NAD 1927（北美基准面1927）
•使用GRS 1980椭球体的NAD 1983（北美基准面1983）
•使用WGS 1984椭球体的WGS 1984（世界大地测量局1984）

根据卫星测量结果开发了新的椭球体，比克拉克（Clarke）在1866年开发的精确度更高。
术语“地理坐标系”和“基准”可互换使用，但是如上所述，GCS包括基准，球体，度量单位和本初子午线。


即使使用相同的地图投影和参数，数据的坐标也会根据这些坐标所基于的基准和椭球体而变化。

例如，下面的地理坐标是针对位于美国华盛顿州贝灵厄姆市内的单个点ng 3个不同的基准：

代码：
DATUM X坐标Y坐标
NAD_1927 -122.466903686523 48.7440490722656
NAD_1983 -122.46818353793 48.7438798543649
WGS_1984 -122.46818353793 48.7438798534299


良好的数据管理原则是从提供数据的数据源中获取投影参数。不要对数据的投影进行有根据的猜测，因为会导致不正确的GIS数据库。必要的参数如下：

•投影
•度量单位
•区域（对于UTM）
•FIPS区域（对于国家平面）
•基准

可能需要其他参数，具体取决于投影。例如，Albers和Lambert投影需要以下参数：

•第一标准平行度，以度，分和秒（DMS）
•第二标准平行度（DMS）
•中央子午线（DMS）
•投影原点的纬度（DMS）
•虚假的东移和度量单位
•虚假的北移和度量单位
•X偏移和度量单位测量
•Y位移和测量单位

#13 楼

请访问以下网站：

http://www.egger-gis.at/shapefile-projectionfinder/

我是这个免费工具的开发者。也许此工具可以帮助您找到并定义shapefile的正确投影。它基于Aaron Racicot的http://projfinder.com/。

您还可以尝试以下特殊解决方案：

奥地利（德语）：https ：//www.data.gv.at/anwendungen/checkaustrianprojection/

澳大利亚：https://maegger.github.io/map_australia.html

UTM-区域： https://maegger.github.io/map_utm.html

#14 楼

在我处理项目（DEM）时遇到了同样的问题。 DEM之一是未知坐标，无法投影，因此我利用ArcMap中的地理配准工具进行了“地理配准”。地理配准后，您将获得prj.file。

#15 楼

您可以区分地理坐标和制图坐标，
第一个纬度/经度以度为单位
第二个北纬度/以度为单位
如果文件为形状，则可以转储几何的坐标使用shapelib-1.3.0，然后查看值...

#16 楼

您可以通过Fiona做到。

    import fiona as f
    a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
    print(a.crs)

您会得到

    {'init': 'epsg:4326'}

的答案