問題タブ [shapely]

私は Shapely に慣れていません (しかし熱心です)。最近、ちょっとした道のりを発見しました。

フィオナ経由で読み込んでいるポリゴン シェープファイルがあります。このシェープファイルには、ポリゴンとマルチポリゴンの両方のアイテムが含まれており、その中のすべての座標 (つまり、外部と内部の両方) の各フィーチャに対して配列を作成する必要があります。特に、2 つのポリゴン アイテムには内部リングがあります (そしてそれらは有効です)。

ポリゴン/マルチポリゴンの外部座標へのアクセスに問題はないようですが、内部座標については何も取得していません。

ここで新しいアプローチを取る必要がありますか (つまり、LinearRings)...?

私は GIS を行うために形を整えて作業していますが、ジオメトリが非常にギザギザで複雑であるため、各郵便番号のジオメトリをメモリにロードするときにメモリ エラーが発生します。

形状をあまり歪ませずに、境界点の数をできるだけ少なくして、メモリ内の形状を小さくしたいと思います。凸包を使用することは、境界から多くのポイントを単に捨てることができるように、1 つの潜在的な答えのように思えます。この問題を解決する何かが既にそこにあるかどうか疑問に思っています。

ジオメトリ コレクションとしてジオメトリ タイプのレイヤーを作成できません。したがって、シェープファイルに出力できません。以下にコードを添付しました。dstLayer がヌルです。これは、geom タイプが多角形または複数線ストリングまたはその他の場合には発生しません。

どんな助けでも大歓迎です。

ありがとう

fiona のインストールに問題があります。

fiona の公式ドキュメント ( https://github.com/sgillies/fiona ) と以下の関連スレッドをフォローしました。

Python cpl_error.h に Fiona をインストールする際の問題: そのようなファイルまたはディレクトリはありません

pip install shapely に問題はなく、brew install gdal を使用しました。コマンド pip install fiona を使用すると、いくつかのエラーが発生します

などの gdal-config に関する警告と同様に"Failed to get options via gdal-config"

問題の原因となっている gdal-config の場所に問題があるのではないかと心配しています。brew install を使用する前に、QGIS で使用するために kyngchaos で GDAL を以前にインストールしていた可能性があると思います。

どんな洞察も大歓迎です。fiona のインストールで問題が発生したのは私だけではないと思いますが、Google 検索だけではまだ解決策を見つけることができませんでした。

Kartograph.py を使用して SVG マップをレンダリングしようとしています。TypeError がスローされます。Pythonコードは次のとおりです。

これがスローされるエラーです

Python 3.4 と shapely 1.3.2 を使用して、長い/緯度の座標ペアのリストから Polygon オブジェクトを作成し、それらを解析するために既知のテキスト文字列に変換しています。このような Polygon は次のようになります。

shapely は射影を処理せず、デカルト空間のすべてのジオメトリ オブジェクトを実装するため、次のようにそのポリゴンで area メソッドを呼び出します。

その多角形の面積を平方度の単位で教えてくれます。平方メートルなどの平面単位で面積を取得するには、別の投影法 (どれ?) を使用してポリゴンの座標を変換する必要があると思います。

pyproj ライブラリがこれを行う方法を提供する必要があることを何度か読みました。pyprojを使用して、形の整ったPolygonオブジェクト全体を別の投影に変換してから面積を計算する方法はありますか?

私は自分のポリゴンを使って他のことを行います (あなたが今考えていることではありません)。特定の場合にのみ、面積を計算する必要があります。

これまでのところ、次の例しか見つかりませんでした: http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/

これは、各 Polygon オブジェクトを外側のリングに分割し、存在する場合は内側のリングに分割し、座標を取得し、座標の各ペアを別の投影に変換して Polygon オブジェクトを再構築し、その面積を計算することを意味します (とにかく、単位は何ですか?)。これは解決策のように見えますが、あまり実用的ではありません。

より良いアイデアはありますか？

緯度/経度座標の長いリストを、それらが属する米国の州 (または郡) に変換したいと思います。状態のジオメトリがある場合、考えられる解決策の 1 つは、すべての状態に対して各ポイントをチェックすることです。

おそらくO（1）で、これを行うためのより最適化された方法はありますか？

1 つの形状に結合しようとしている 7 つの重なり合う円と楕円のクラスターがありますが、実行するcascaded_union()とエラーが発生します。

これまでに書いたものは次のとおりです。

私の目標は、右の図のようなものを取得することです: コードの残りの部分は、不規則な形状の境界内にランダムな分布のポイントがいくつあるかを決定します。「null値」コメントが返されたときに、エラーが何を参照しているのか混乱しています。個々の形状の重なりを正しく考慮していませんか? 私が検索したものから、すでにcascaded_union形状の配列の入力を取得していますが、何らかの理由でこの場合は機能していません。