問題タブ [geopy]

これは、 Python での地理空間分析へのフォローアップの質問です。

質問はhttps://gis.stackexchange.com/questions/84114/shapely-unable-to-tell-if-polygon-contains-pointに似ていますが、緯度を逆にすると問題は解決したようですが、そうではありません私を助けて。

Uber のデータは緯度/経度で与えられますが、これは非常に簡単です。geopandas で使用すると、逆ルックアップで住所が得られます。

ただし、問題は、シェープファイルで緯度/経度を検索することです。Uber のデータは次のようになります。

形状ファイルのポリゴン境界は次のようになります

したがって、 polygon/polygon.contains(point) のポイントが機能していません。データを見ると、緯度経度は形の整ったファイルと比較して非常に小さく、ある単位を別の単位に変換する必要があるかどうかわかりません。まったく異なるメートル法に見えます:) 以下はコードです:

Uber データと SFO シェープファイルへのリンクはこちらhttp://hortonworks.com/blog/magellan-geospatial-analytics-in-spark/#comment-606532

geopy ライブラリを使用して場所を操作します。https://pypi.python.org/pypi/geopy

初めてインストールして使用しましたが、問題なく動作しました。次の数日、同じコードを実行したところ、例外が発生しました。

また、ライブラリを再インストールして、インタープリターのパスを再設定しようとしました。Windows 7、Eclipse Mars.1 Release (4.5.1) を pydev で使用しています。Python のバージョンは 2.7.10 です。

私はこのpythonスクリプトで供給しているこのCSVファイルを持っています

CSVファイルの抜粋

しかし、元のcsvは次のようになります

これには 15k 以上の行があり、お気づきのように、いくつかの場所には同じエントリが複数回含まれています。私のワークフローは、sqlite データベースからこの大きな csv にデータを取得することでした。その後、二重のエントリをフィルタリングしたので、すべての場所/州を一度だけ取得しました。したがって、ジオコーダーのクエリを減らすことができました (制限がまったくないと考えるのはとても単純でした)。私の問題は、緯度/経度情報を使用した単純な「すべてのエントリの検索と置換」がまだ困難で時間のかかる作業であるため、再度まとめることができない 2 つのデータセットがあることです。また、エントリ番号が一致しないためです。私には2つのアイデアがあり、どちらもPythonを初めて使用し、このデータに終日取り組んでいるため、そこにたどり着く方法がわかりません。

1. (都道府県/エリア)を(都道府県/エリア、緯度、経度)で検索・置換する処理を自動化

2. 事件ごとに、同じ日に同じ場所で亡くなった人が何人かいます。技術的にそれらを合計して、データ全体が小さくなるようにすることができます（どこでこれを行う必要がありますか。csvを使用したPythonのSQLiteまたはPythonで）。したがって、このデータをジオコーダーに照会して、ファイル全体を再び取得できます。州と地域をジオコーダーに問い合わせると、次のようになります。

   /li>

あなたのアイデアをありがとう

編集1

編集3

変更を試みましたが、機能しません

これらのエラーで

私のCSVには、個人、死亡日、死亡地域に関する情報が含まれています。また、geopy を使用してこの領域を特定しています。この部分は機能していると思いますが、実際には個人を知りたくありませんが、1 日あたりのインシデントごとおよびエリアごとの合計ボディカウントが必要です。そのため、死亡者数、地理位置情報、日付、および発生した地域/州に関する情報を含む CSV が得られます。これが私の CSV のスニペットです。

これが私のコードです

機能していないために何が表示されないのですか?

ターミナルにこのメッセージが表示されます

geopy を使用して、経度と緯度に基づいて 2 つの住所間の距離を計算すると、個々のデータ ペアで正常に機能します。しかし、より多くのデータがある場合、常に次のエラーが表示されます。

ファイル "/Library/Python/2.7/site-packages/geopy/geocoders/osm.py"、193 行目、ジオコードの self._call_geocoder(url, timeout=timeout), exact_one ファイル "/Library/Python/2.7/site- packages/geopy/geocoders/base.py", line 171, in _call_geocoder raise GeocoderServiceError(message) geopy.exc.GeocoderServiceError: urlopen エラー [Errno 65] ホストへのルートがありません

どうすればこの問題を回避できるか知っていますか?

私のコードは単純です: (このためのデータ入力には多くのデータのペアがあります)

geopyこのオブジェクトの出力をファイルに書き込もうとしていcsvますが、各文字が別の列に配置され、緯度と経度が別の行に出力されます。どうすれば修正できますか？

この関数を別の時間に実行して、新しいアドレスを次の行に出力できるようにしたいと考えています。データを上書きせずに保存します。csvこれは Pythonでの書き込みで実行できますか?

github からダウンロードしたコードを使用して、geopy を使用して住所をジオコーディングしようとしています。しかし、コードのエラーは、str obj がアイテムの割り当てをサポートしていないということgeocode_addresses address_dict["error"] = ""ですか? なぜそのエラーが発生するのですか。以下はコードです。

データセットの例