回答:
常に「ブルートフォース」方式があります：
未知のレイヤーとオーバーレイするはずの既知の座標系を持つレイヤーを取得します。
次に、未知のレイヤーがどのような投影になる可能性があるかについて、知識に基づいた推測を行います。（UTM、Plate Carreeなど）。既知の座標系レイヤーを各投影に投影し、未知のレイヤーに可能な限り一致するものを見つけます。
MapperzとBrandon Copelandのコメントを、そのテクニックを使用した回答を追加することで促進しています。
この手法は、ArcMapがデータをメモリ内の新しい座標系に変換できるという事実を利用しています。既知の座標なしのデータは変換できないため、表示されるだけです。既知のデータをメモリに投影することにより、さまざまな可能な座標系でどの座標を持つかをすばやく確認できます。
この質問は決して古くならないので、ブルートフォースメソッドを実行するサイトを構築しました。圧縮されたshp + shxをマップにドラッグすると、PostGISで使用可能なすべての座標系にマップされます。「正しい」がどのように見えるかを知っていると仮定すると、その領域にズームし、ポリゴンをクリックして、epsg.ioから.prjファイルを取得できます。
Esriには、これに関する詳細なリンクが2つあります。
私の最後の仕事で、「FSU_Geol.shp」という名前の地質学データ（ポリゴン）のレイヤーがありました。私の上司は私にそれを与え、いくつかのことを見つけるように頼みました。最初に、彼はクライアントからこのシェープファイルを渡されましたが、.prjファイルがなかったため、彼にそれを理解してもらいたいと考えました。彼はまた、地質学のカテゴリーが何であるかを私に理解してほしかった。私はどんどん行くことができましたが、追いかけましょう....私はそれをグーグルで調べました、そして、私はここに行きました。
命名規則の頭字語は「旧ソビエト連邦の地表地質学」を表しており、Googleの検索により「直接」ソース（USGS）に導かれました。このシェイプファイルについて知っておくべきことはすべて、私がヒットしたトップリンクにありました。Googleが何かを見つけることができると言っているわけではありませんが、大学を卒業したばかりで、「暗闇で撮影」して、得たフィードバックを見てください！
別の例では、オフィスの誰かがバッチ機能を使用して多くのシェープファイルをダウンロードしていました。現時点ではソフトウェアパッケージの名前は覚えていませんが、.prjファイルが見つかりませんでした。ArcCatalogのメタデータにアクセスしただけで、実際に答えが見つかりました。繰り返しになりますが、引用はしませんが、そこにはソースへのURLがあったと思います。
その日以来、メタデータに手がかりがなければ、Google検索を「開始」します。
http://projfinder.com/の Webツールが大好きです。ファイルを確認して、座標を調べてください。地球上のほぼ同じ場所にズームして、推測させます。
ツールではありません（それを可能にするツールは存在しません）が、このポイントをWKID 4326に変換するにはどうすればよいか、@ mkennedyの回答をご覧ください。。彼女は正しい空間参照に到達した方法を説明します。 SpatialReference.orgと忍耐があなたの友達になります。
さらに、Esriは座標系を推測する方法についてのガイドを提供します（ただし、データについてもう少し知っている場合は、mkennedyの方法を好みます）。
この質問に対する私が見つけた最良の答えは技術的ではありません。データがどこから来たのかを見つけてください。代理店と組織は、予測の使用と一貫している傾向があります。それがあなたの州のDOTから来たことを知っていますか？残りのデータを見て、それが何を伝えているかを見てください。どこから来たのか分からない？経験に基づいた推測は、あなたを正しい道に導きそうです。
少なくとも、ブルートフォースで問題に対処することがもう少し実行可能になります！
この時点では非常に時代遅れですが、Werner FlackeとBirgit Klaus は2007年にArcScriptsにFind Projectionを投稿しました。残念ながらソースコードはそこにあるとは思いません。VBAベースなので、ArcGIS Desktop v9.2および場合によっては9.3でのみ使用できます。可能性を狭めるために使用できるEPSG Geodetic Parameter Datasetの関心領域を持つ2つのシェープファイルが含まれています。
Blue Marble Geographic CalculatorとGeographic Transformerには、座標系回復ツールがあります。
GDALの一部として提供されているogrinfoプログラムを試してください。
OGRを使用してシェープファイルメタデータにアクセスする方法をご覧ください。
そのため、ogrinfoは.prjファイルなしでは投影情報を提供しないかもしれませんが、潜在的な投影のリストを調査するのに役立つ便利なツールです。
例えば：
Geometry: Polygon
Feature Count: 269
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000)
この応答は、シェープファイルが度ではなくメートルに基づく参照システムを使用していることを示唆しています。
データがどこから来たかを知っていると仮定すると、可能性のある予測のリストが短くなりました。
その他のベストプラクティスは次のとおりです。
ゼロから力ずくで試行するよりも、既知のある位置から開始する方がおそらく良いでしょう。
この投稿の他の応答に追加するには、次を追加します。
NAD1927に特定のチェックを提供する
座標が経度-180〜+180、緯度-90〜+90などの10進度の場合、データに使用される地理座標系（データム）を特定します。バージョン9.2より前のArcMapは、デフォルトでGCS_Assumed_Geographic_1をデータに割り当てます。これにより、データがNAD_1927データムに配置されます
加えて
データが米国にあり、小数の左側の座標が6、7、または8桁の範囲を示す場合、データはおそらく州平面またはUTM座標系に投影されます。
FAQ：投影の基本：GISプロフェッショナルが知っておくべきこと
このページには、不明な座標系を特定するのに役立つ投影法の決定に関するいくつかのチェックがあります。
次の概念は、ArcGISでの地図投影法の使用を理解するための基本です。ただし、予測のトピックは非常に広範であり、この記事ではいくつかの重要なトピックに触れることしかできないことに注意してください。
マップ投影とも呼ばれる座標系は、空間データの任意の指定です。それらの目的は、地球の表面の特定の場所または領域に関するコミュニケーションの共通基盤を提供することです。マップ投影を扱う際の最も重要な問題は、投影が何であるかを把握し、データセットに関連付けられた正しい座標系情報を取得することです。
最初の地図投影法が考案されたとき、地球は平らであると誤って想定されていました。その後、仮定が修正され、地球は完全な球体であると仮定されました。18世紀に、人々は地球が完全に丸くないことを認識し始めました。これは、地図製作スフェロイドの概念の始まりでした。
地球の表面上の位置をより正確に表すために、地図作成者は地球の形状（測地）を研究し、回転楕円体の概念を作成しました。次に、データム、測定単位、および子午線を含む地理座標系（GCS）が考案されました。データムは、回転楕円体を地球の表面の特定の部分にリンクします。最近のデータムは、地球の表面全体に適合するように設計されています。
北米で最もよく使用されるデータムは次のとおりです。
•Clarke 1866スフェロイドを使用したNAD 1927（北米データム1927）
•GRS 1980スフェロイドを使用したNAD 1983（北米データム1983）
•WGS 1984 スフェロイドを使用したWGS 1984（World Geodetic Survey 1984）新しいスフェロイドは衛星測定から開発され、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データベースが結果となるため、データの投影について経験に基づいた推測を行わないでください。必要なパラメーターは次のとおりです。
•投影
•測定単位
•ゾーン（UTM用）
•FIPSゾーン（州平面用）
•データム投影法に応じて、他のパラメーターが必要になる場合があります。たとえば、AlbersおよびLambertの投影には、次のパラメーターが必要です。
•度、分、秒で第一標準パラレル、（DMS）
•第二の標準的なパラレル（DMS）
中央子午線•（DMS）
•投影原点（DMS）の緯度
偽東距と測定単位•
測定の偽北上し、単位•
• Xシフトと測定単位
•Yシフトと測定単位
このウェブサイトを見てください：
http://www.egger-gis.at/shapefile-projectionfinder/
私はこの無料ツールの開発者です。このツールは、シェープファイルの正しい投影を見つけて定義するのに役立つかもしれません。これは、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
あなたはフィオナでそれを行うことができます。
import fiona as f
a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
print(a.crs)
あなたは答えを得るでしょう
{'init': 'epsg:4326'}