Answers:
总有一种“蛮力”方法：
拍摄具有已知坐标系的图层，该图层应与未知图层重叠。
现在对未知层的投影进行一些有根据的猜测。（UTM，Plate Carree等）。将您已知的坐标系图层投影到每个投影中，直到找到与未知图层尽可能匹配的图层为止。
我通过添加使用他们的技术的答案来促进Mapperz和Brandon Copeland的评论。
该技术利用了ArcMap可以将数据转换为内存中新坐标系这一事实。没有已知coordsys的数据无法转换，因此只能显示。通过将已知数据投影到内存中，我们可以快速查看其在各种可能的坐标系中的坐标。
Esri有两个很棒的链接，对此进行了详细说明：
在上一份工作中，我有一层名为“ FSU_Geol.shp”的地质数据（多边形）。我的老板把它给我，并要求我找出一些事情。首先，客户端将他交给了这个shapefile，并且没有.prj文件，所以他希望我自己弄清楚。他还希望我弄清楚什么是地质学。我可以继续前进，但是随它去吧……。我用谷歌搜索了，最后到了这里。
命名约定中的首字母缩写词代表“前苏联的地表地质”，而Google搜索将我“直接”引到了来历（USGS）。我可能需要了解的有关shapefile的所有信息都位于我点击的顶部链接。我并不是说Google可以找到任何东西，但我刚大学毕业，只是“暗中摸索”，然后看看我得到的反馈！
在另一种情况下，办公室中的某人已使用批处理功能下载了许多shapefile。我现在不记得软件包的名称，但是缺少.prj文件。我只是进入ArcCatalog中的元数据，并且实际上在那里找到了答案。同样，不要引用我，但是我认为那里有一个指向源的URL。
从那一天开始，如果元数据中没有任何线索，我会从Google搜索开始！
我真的很喜欢这个网络工具：http : //projfinder.com/。检查您的文件并查找一些坐标。放大到地球上大约相同的位置并进行猜测。
不是一种工具（我不知道有一种工具可以让您做到这一点），但请查看@mkennedy的回复，我将如何将此点转换为WKID 4326？。她解释了如何得出正确的空间参考。 SpatialReference.org和耐心将成为您的朋友。
此外，Esri提供了有关如何猜测坐标系的指南（尽管如果您对数据有更多了解，我更喜欢mkennedy的方法）。
对于这个问题，我找到的最佳答案是非技术性的：找出数据的来源。机构和组织倾向于与其使用预测相一致。知道它来自您所在州的DOT吗？查看他们的其余数据，看看它能告诉您什么。不知道它来自哪里？有根据的猜测很可能使您走上正确的道路。
至少这使得用蛮力解决问题变得更加可行！
此时严重过时，但是Werner Flacke和Birgit Klaus 于2007年在ArcScripts上发布了Find Projection。不幸的是，我认为源代码不存在。它基于VBA，因此仅可在ArcGIS Desktop v9.2和9.3中使用。它的确包含两个shapefile，它们具有EPSG大地测量参数数据集所关注的区域，可用于缩小可能性。
Blue Marble Geographic Calculator和Geographic Transformer具有坐标系统恢复工具。
尝试使用GDAL随附的ogrinfo程序。
因此，如果没有.prj文件，ogrinfo可能不会为您提供投影信息，但是它仍然是帮助您调查潜在投影列表的有用工具。
例如：
Geometry: Polygon
Feature Count: 269
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000)
该响应向我表明shapefile正在使用基于米而非度的参考系统。
假设您知道数据的大致来源，那么现在可以得到的投影比较短。
其他最佳做法可能是：
最好从一些已知的位置开始，而不是从头开始尝试蛮力。
要添加到本文中的其他响应中，我将添加以下内容：
提供有关NAD1927的特定检查
如果坐标以十进制度为单位，例如在-180和+180之间，在-90和+90之间，则标识用于数据的地理坐标系（基准）。在9.2版之前，默认情况下，ArcMap将GCS_Assumed_Geographic_1分配给数据。这会将数据放置在NAD_1927基准面上
此外
如果数据在美国，并且显示小数点左边的坐标是6、7或8位数字的程度，则该数据可能会投影到州平面或UTM坐标系。
该页面提供了一些确定投影的检查，这对协助识别未知坐标系非常重要。
以下概念是理解ArcGIS中地图投影用法的基础。请注意，尽管预测的主题非常广泛，但本文仅涉及几个重要主题。
坐标系，也称为地图投影，是对空间数据的任意指定。它们的目的是为地球表面上特定位置或区域的交流提供通用基础。处理地图投影时，最关键的问题是知道投影是什么，并具有与数据集关联的正确坐标系信息。
在设计第一个地图投影时，错误地认为地球是平坦的。后来对假设进行了修改，并且假设地球是一个完美的球体。在18世纪，人们开始意识到地球并不是完美的圆形。这是制图球体概念的开始。
为了更准确地表示地球表面的位置，地图制作者研究了地球的形状（大地），并提出了椭球的概念。然后设计了地理坐标系（GCS），其中包括一个基准，度量单位和一个本初子午线。基准将椭球体连接到地球表面的特定部分。最近的基准面被设计为非常适合整个地球的表面。
北美最常用的基准是：
•使用Clarke 1866椭球的NAD 1927（北美基准1927）
•使用GRS 1980椭球的NAD 1983（北美基准1983）
•使用WGS 1984椭球的WGS 1984（世界大地测量1984）较新的椭球体是根据卫星测量结果开发的，比Clarke在1866年开发的椭球体更精确。术语“地理坐标系”和“基准”可互换使用，但是如上所述，GCS包括基准，椭球体，度量单位和本初子午线。
- 即使数据坐标使用相同的地图投影和参数，数据坐标也会根据这些坐标所基于的基准和椭球体而变化。
例如，下面的地理坐标是使用3种不同基准的华盛顿贝灵汉市内单个点的：
代码：DATUM X坐标Y坐标NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299
- 良好的数据管理原则是从提供数据的数据源中获取投影参数。不要对数据的投影进行有根据的猜测，因为会导致不正确的GIS数据库。必要的参数如下：
•投影
•度量单位
•ZONE（对于UTM）
•FIPS区域（对于国家平面）
•基准根据投影，可能需要其他参数。例如，Albers和Lambert投影需要以下参数：
•第一标准平行度，以度，分钟和秒（DMS）为单位
•第二标准平行度（DMS）
•中央子午线（DMS）
•投影原点的纬度（DMS）
•虚假的东距和度量单位
•虚假的北距和度量单位
• X位移和度量单位
•Y位移和度量单位
看这个网站：
http://www.egger-gis.at/shapefile-projectionfinder/
我是这个免费工具的开发商。也许此工具可以帮助您找到并定义shapefile的正确投影。它基于Aaron Racicot 的http://projfinder.com/的思想。
您也可以尝试以下特殊解决方案：
奥地利（德国）：https：//www.data.gv.at/anwendungen/checkaustrianprojection/
您可以通过Fiona做到。
import fiona as f
a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
print(a.crs)
您将得到答案
{'init': 'epsg:4326'}