Více

Spojovací body umístěné na opačné straně mezinárodního poledníku pomocí Pythonu GDAL / OGR?

Spojovací body umístěné na opačné straně mezinárodního poledníku pomocí Pythonu GDAL / OGR?


Vzhledem k seznamu souřadnic musím vytvořit liniový tvar spojující každý bod pomocí skriptu Python GDAL / OGR. Zde popsaný postup Převést XY bodů na čáru? funguje dobře téměř ve všech případech, kromě případů, kdy body překračují poledník změny data. Když k tomu dojde, čára spojí dva body křížením celého světa, což je zjevná chyba.

Jak mohu vynutit skript, aby spojil 2 body na nejkratší vzdálenost?

Zde je několik příkladů souřadnic:

lon, lat 172,5, 50,0 175,0, 43,4 -175,0, 35,0 -173,1, 31 -172,1, 28,0 -171,3, 26,5 -173,9, 21,8

Na nejkratší vzdálenost se možná budete chtít podívat do výpočtu velkých kruhů.

Zde je jeden příklad provedený v QGIS: velké kruhy v QGIS a export do 3857 webové mapy

Místní CRS se na datovém řádku nerozbije.

Pokud je vaše oblast zájmu kolem datové čáry, můžete použít místní projekci příčného mercatoru:

+ proj = tmerc + lat_0 = 50 + lon_0 = 172,5 + x_0 = 0 + y_0 = 0 + a = 6371000 + b = 6371000 + jednotky = m + no_defs

nebo se řiďte mými radami na QGIS display world country shape files centered on the Pacific ocean using Robinson, Miller Cylindrical or other projection

Chcete-li své řádky zobrazit v „normální“ projekci, musíte svoji geometrii čáry zhutnit a oříznout ji na datové čáře (podobně jako u řádků u reprojektovaných sp objektů s projekcí Mollweide)


Podívejte se na video: CSV File to Point Shapefile - Python GDALOGR