問題タブ [shapely]

Windows上のPythonでOpenCV 3アルファを使用しています。グラブカットを使用した画像セグメンテーションにつながるバックグラウンド減算法があります。そのため、前景と背景の可能性のある情報を提供するMOG検出器があります。たとえば、ここに現在の画像があります（視覚化のために与えられた四角形）。

そして、これがMOG検出器からの出力です。

この情報を cv2.grabcut にフィードしたいと考えています。画像全体をセグメント化する必要がなく、既知のオブジェクトの周囲の領域を指定して、前景と背景の可能性が高いものを渡す方が高速 (?) になることを願っています。ブロブは、境界が xmin、ymin、xmax、ymax である形の良いポリゴンとして格納されます。

これは常に私に黒いイメージを与えます。すべての背景。

cv2.GC_INIT_WITH_MASK で初期化すると、正しく動作します (赤い四角は無視してください)。ただし、四角形の境界のかなり外側に推定された前景が含まれることがあるため、四角形は確実に無視されます (この場合は表示されません)。

rectを間違って保存していますか? x,y,w,h じゃないの？rect を指定すると実際に高速になりますか、それとも画像をトリミングする必要がありますか?

形の良いポリゴンの外側のポイントを削除したい 10 x 10 のグリッドがあります。

結果の図は次のとおりです。

そして、最終結果を次のようにしたい:

配列を 100 x 2 のグリッド ポイントの配列に変更できることはわかっています。

ポイントがポリゴン内にあるかどうかを簡単に確認できます。

しかし、この情報を取得して元のグリッドを変更するのに問題があります

これら 2 つのオブジェクトはLineStringfromを使用して作成されていますshapely.geometry

を使用して交点を見つけることができますが、これらの 2 つの線の交点が自動的に行われている のLine1_string.intersection(line2_string).coords線分も定義したいと思います。Line_string1

これまでのところ、この回答 [リンク] ( Point の両側の LineString の頂点を取得する)に出くわしましたが、定義された分割関数から何も出力されません。

何が間違っているのか、または他の提案についてのアイデアはありますか?

で簡単に視覚化できる (緯度、経度) ポイントのセットがありますmatplotlib.pyplot.scatter(list_of_longitudes, list_of_latitudes)。しかし、たまたまシェープファイルにある道路網の上にこれらのポイントをプロットしたいと思います。2 つのファイルを Python と一緒に 1 つのマップ/図に単純にプロットするにはどうすればよいですか?

特にShapelyのようなものを探しています