問題タブ [geopandas]

データのインポートに GeoDataFrame を使用しています。しかし、以下の問題があります。実際、この関数は一部のシェープファイルではうまく機能しますが、一部の特定のシェープファイルではうまく機能しません。

geopandas を使用したデータの読み取りに関する問題がありますが、この投稿によると、エラーは fiona が原因のようです。 fiona を 1.1.6 バージョンに再インストールしようとしましたが、それでも機能しません。

src として fiona.open('taz.shp') を指定して fiona をインポートします: for feature in src: print feature

しかし、他の同様のスクリプトを書くと、うまくいくようです

2 つのバラバラなグループ / ポリゴンの「島」があるとします (隣接していない 2 つの郡の国勢調査区を考えてみてください)。私のデータは次のようになります。

各島内のポリゴンが、そのグループを表す ID (任意の可能性があります) を取得するようにします。たとえば、左下の 3 つのポリゴンは IslandID = 1 を持つことができ、右上の 3 つのポリゴンは IslandID=2 を持つことができます。

これを行う方法を開発しましたが、それが最善/最も効率的な方法かどうか疑問に思っています。私は次のことを行います：

1) マルチポリゴン単項ユニオン内のポリゴンに等しいジオメトリを持つ GeoDataFrame を作成します。これにより、「島」ごとに 1 つずつ、合計 2 つのポリゴンが得られます。

2) 各グループの ID を作成します。

3) 島を元のポリゴンに空間的に結合し、各ポリゴンが適切な島 ID を持つようにします。

これは実際にこれを行うための最良/最も効率的な方法ですか?

この質問は、私の前の 2 つの質問に関連しています。

Python での地理空間分析

Python 2.6 に geopandas をインストールする

システム全体をアップグレードすることで、ジオパンダが機能するようになりました。フィオナを機能させることができません。以下は、試したときに表示されるエラーです

ラインストリングの GeoDataFrames に次のものがあり、そのうちの 1 つは道路を表し、もう 1 つは等高線を表すとします。

これらをプロットすると、次のようになります。

3 つの等高線と 2 つの道路があります。各道路に沿った各ポイントの標高を見つけたい。基本的に、道路と等高線 (これにより 12 ポイントが得られるはずです) を交差させ、両方のジオデータフレーム (道路名と標高) の属性を保持したいと考えています。

2 つのジオデータフレームの和集合の交点を使用して、12 のポイントを生成できます。

しかし、これら 12 のポイントそれぞれの道路名と標高を取得するにはどうすればよいでしょうか? 空間結合は期待どおりに動作せず、4 つのポイントのみを返します (定義上はそのように作成されているため、12 個すべてがライン ファイルと交差する必要があります)。

これを行う方法について何か提案はありますか?

編集:問題は交点の作成方法に関係しているようです。道路と等高線をごく少量だけバッファすると、交差点は期待どおりに機能します。下記参照：

上記は目的の出力ですが、線の交点から作成された点と交差するように線をバッファリングする必要がある理由はわかりません。

geopandas を使用して何かをプロットしようとしており、bokeh を使用して同様のマッピングを実行しようとしていましたが、同じエラーが発生し、グーグルで何時間もかけて修正できなかったようです。

virtualenvでipythonノートブック（python 2.7.10）を実行しているMac OS X 10.10.5で実行しています。システムに ActiveTLC 8.6.4 をインストールしました。私は困惑しています

2 つのジオデータフレームをマージしようとしています (各ポイントがどのポリゴンにあるかを確認したい)。

次のコードでは、最初に警告 (" CRS does not match!") が表示され、次にエラー (" RTreeError: Coordinates must not have minimums more than maximums") が表示されます。

何が間違っているのですか？CRSは座標系ですか? もしそうなら、なぜそれらは同じようにロードされないのですか?

再現用データへのリンク: Shapefile、 GPS ポイント