問題タブ [geopandas]

2 つの x アレイ データセットがあります。

• 1 つは呼び出され、変数を含むfooディメンションを持ちます。 x,y,timea,b,c,d
• 2 番目が呼び出されbar、 の時系列分析の結果ですfoo。
  • これには、ディメンションと、 のすべての値がと の間にx,yある単一の変数しかありません。 ee01

で処理されるデータ エントリbarの量をフィルタリングまたは制限するために使用したいと思います。ここでは、 at eachののみを考慮する必要があります。 x,y,timefoox,y,timee > 0.8x.y

私はまだXアレイにかなり慣れていません。私の質問は、2 つの xarray データセットを何らかの形で AND 演算または交差させ、より小さなデータセットを期待するという「セット操作」の考え方が、x-array について推論する現実的な方法であるかどうかです。

これが私がこれまでに持っているものです。

foo = loadDataset()
bar = perform_timeseries_analysis()
filtered_bar = bar > 0.8
#TODO: Use bar to reduce the size of foo

mssql 空間フィールドを geopandas に直接インポートできないようです。通常の mssql テーブルを Pymssql を使用して問題なく pandas にインポートできますが、空間フィールドを形状の良いジオメトリにインポートする方法がわかりません。mssql の OGR ドライバーで処理できることはわかっていますが、これを理解するには SQL の知識がありません。ポイントは mssql フィールドから x 座標と y 座標に変換される可能性があるため、これはラインとポリゴンにとってより大きな問題です。ありがとう！

python / を使用して、ポイントと国の国境の間の距離を計算したいshapely。point.distance(poly) で問題なく動作するはずです。たとえば、ここで示されています。Polygon Shapely の最も近いポイントの座標を検索しますが、使用するgeopandasと次の問題に直面します。 'GeoSeries' object has no attribute '_geom'

データの処理で何が問題になっていますか? 私の国境データセットはhttp://www.gadm.org/からのものです

ジオメトリ フィールドが単一の(x,y) に設定されたgeopandasデータフレームがあります。次に、ジオメトリ フィールドがポリゴンとその他の情報に設定された 2 番目のデータフレームがあります。例えば：APointB

と

のポイントが のポリゴンの内側にある場合にのみ、 AwithBのフィールドに新しい列を追加するにはどうすればよいですか?infoAB

私は次のようなもので終わりたいと思います

（GeoPandasでロードされた）とPandasデータフレーム（プロットしたいデータを含む）のmerge間にあるこのGeoPandasデータフレームがあります。shapefile

与えます：

GeoPandas のドキュメント (ここ: http://geopandas.readthedocs.io/en/stable/mapping.htmlおよびgpd.plot関数のドキュメント) によると、クロロプレス マップをプロットしたい場合は、次のようにする必要があります。

しかし、これにより次のエラーが発生します。

matplotlibで作業するために何をする必要があるのか​​ わかりませんかGeoPandas？

ポイントの列の最初の座標の値 (x 値) に基づいて geopandas データフレームをサブセット化したいと思います。

したがって、上記のデータが与えられた場合、最初の座標が (たとえば) 3 より大きいかどうかを確認したいと思います[False, True, False]。これにより、元のデータフレームをサブセット化できるように の配列が返されます。

ポリゴンを含むシェープファイルのセットが 2 つあります。シェープファイルの 1 つのセットは、私が興味を持っている米国の郡だけであり、これは会社や年によって異なります。シェープファイルのもう 1 つのセットは、企業の事業領域であり、もちろん、これは企業や年によって異なります。毎年、各企業のこれら 2 つの層の交点を取得する必要があります。これまでのところ、関数 overlay(df1, df2, how = 'intersection') は私の目標を達成しました。しかし、企業年ごとに約300秒かかります。私は長い会社のリストと何年にもわたるリストを持っていることを考えると、これを完了するには数日かかるでしょう. このパフォーマンスを向上させる方法はありますか?

同じことを ArcGIS で行うと、300 秒が数秒に短縮されることに気付きました。しかし、私は ArcGIS の新しいユーザーであり、まだ Python に慣れていません。

私は簡単な演習のために geopandas ライブラリをテストしています: マップ上にいくつかのポイントを表示し、その上に大きな円を重ねて差分メソッドでそれらの一部を削除します。

変換が正常に機能することを確認するために、iPython ノートブックを使用してさまざまなレイヤーを表示しています。

だから、これが私の操作の始まりです：

次に、私のポイントを確認するために、次のように書きました。

結果は次のとおりです。

それはとてもいいです。ここで、このレイヤーに大きな半径のポイントを追加したいだけです。私はこれを試しました：

次に、これらのコマンドで十分だと思いました：

しかし、グラフィカルな表示の代わりに、私のノートブックは以下を返します:

そして、これまでのところ何が間違っているのかわかりませんでした...