国際日付変更線を横切るポリゴン[-180 .. + 180経度]

9

衛星軌道帯用のポリゴンを生成しようとしています。これまでのところ、[lat、long]で各帯のエッジを表す2つの線を生成する方法があります。一部のスワスは国際日付変更線と交差しているため、次のように囲みます。

私はこれを解決することができましたogr2ogr -wrapdateline：

ogr2ogr -wrapdateline  -f "ESRI Shapefile" test.shp orbits.shp

おそらく行を分割します

次に、両方の線の内側にポリゴンを生成できるようにしたいと思います。したがって、たとえば、スワスの1つのエッジが日付ラインと交差する場合、次のように、ポリゴンが反対側に現れたときに塗りつぶします。

タスクを何度も繰り返す必要があるため、自動化された方法が必要です。できればpythonで行を生成します。以下は、これらの行を含む2つのシェープファイルです。wraparound ; 固定日付

— ジェームズ
ソース

追加のアイデアについては、gis.stackexchange.com / questions / 429およびgis.stackexchange.com/questions/18562で関連するスレッドを参照してください。おそらく、gis.stackexchange.com / questions / 17788で提示されたアイデアも参考になるでしょう。しかし、「内部」とはどういう意味かと思います。これらのポリゴンは明確に定義されていないため、少なくとも、（a）各ポリラインのどちら側が「内部」と見なされるか、および（b）方法を示す情報を提供する必要があります。極の近くでそれらを切り離すため。

— whuber

2

スワスのほぼ中心に配置されたカスタムメルカトル図法を作成できます。たとえば、スワス25に使用します。

+proj=merc +lon_0=-140 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

この予測では、帯は日付変更線によって分割されていません。ラインからポリゴンを作成できます。

次に、EPSG：4326で-179.95°Eと179.95°Eの間のカットポリゴンを作成します。

Nr;WKT
1;POLYGON ((-179.95 89, 179.95 89, 179.95 -89, -179.95 -89, -179.95 89))

それをカスタムCRSにも再投影し、スワスポリゴンから減算します。

EPSG：4326に再投影した後、帯は日付線で正しく分割されます。

日付変更線を横切るすべての帯を続行します。

— AndreJ
ソース

1

帯状線の生成プロセスを書き直して、同じ連続した縦長の空間で開始および終了します。つまり、ラインが170°で始まり-170°で終わった場合、プロセスを書き直して、-180,180でラップせずに190°で終わるようにします。

次に、線の間に切れ目のないポリゴンを作成できます。

次に、クリッププロセスを使用してポリゴンを180、-180ラインで分割し、-180、180スペースの外側にあるパーツを、360度適切に加算または減算してシフトします。

特定のプロジェクション/データムで保存する前に、すべて完了してください

— パープル氏
ソース