Več

Določitev azimutne enako oddaljene projekcije vremenske radarske karte Environment Canada


Okolje Kanada koristno zagotavlja nedavne zajeme vseh državnih vremenskih radarskih postaj na naslovu http://dd.meteo.gc.ca/radar/PRECIPET/GIF/ v intervalih po deset minut. To so iste slike GIF kot tiste, objavljene na njihovi vremenski spletni strani (npr. Http://weather.gc.ca/radar/index_e.html?id=xsm) in vstavljene spodaj.

To vprašanje koristno opisuje, kako prenesti trenutne vremenske radarske podatke v formatih KML iz okolja Kanada. Vendar moram avtomatizirati obdelavo globokega arhiva slik GIF, posnetih v intervalih po deset minut.

Tudi jaz razumem, da so surovi radarski podatki v prostorskih oblikah na voljo pri podjetju Environment Canada po cenah za povrnitev stroškov, ki bodo imele veliko višjo tematsko in prostorsko ločljivost, kot je mogoče od teh rasterjev GIF. Vendar taka natančnost za to aplikacijo ni pomembna.

Nato je izziv opredeliti koordinatni sistem teh GIF-ov v znani projekciji. Lokacija radarja v središču zemljevida ("+") je znana in objavljena (51.0246, -113.3981 za radarsko postajo na zgornji sliki). Na podlagi poizvedb v Kanadi za okolje sem ugotovil tudi, da so te radarske slike v azimutni enako oddaljeni projekciji.

Kako bi nekdo oblikoval niz PROJ4 ali v ArcGIS definiral projekcijo, ki bi opisala te slike GIF v njihovem nepreoblikovanem stanju, in tako izkoristil ta brezplačni vir.


ArcGIS podpira krogle ali elipsoide (sferoide) z azimutno enako oddaljeno projekcijo.

Tu je vzorec niza prj (wkt) za "World_Azimuthal_Equidistant", ki temelji na WGS 1984.

PROJCS ["World_Azimuthal_Equidistant", GEOGCS ["GCS_WGS_1984", DATUM ["D_WGS_1984", SPHEROID ["WGS_1984", 6378137.0,298.257223563]], PRIMEM ["Greenwich", 0.0], UNIT [1997] PROGRAM ["Azimuthal_Equidistant"], PARAMETER ["False_Easting", 0.0], PARAMETER ["False_Northing", 0.0], PARAMETER ["Central_Meridian", 0.0], PARAMETER ["Latitude_Of_Origin", 0.0], UNIT ["Meter", 1.0] ]

Če želite podati kroglo, nastavite drugi številski parameter v ključni besedi SPHEROID na 0 (nič). Kot je AndreJ dodal v komentarjih, bi spremenili osrednji poldnevnik in vrednosti parametrov zemljepisne širine, da se ujemajo z določenim radarskim mestom.

Če vsi gifi predstavljajo isto območje v azimutni enako oddaljenosti, bodo koordinate vogalov (ki jih boste morali georeferencirati), enake, ne glede na koordinate središčne točke.


Določiti morate projekcijo na kroglo, ker se zdi, da proj podpira le sferične formule te projekcije:

+ proj = aeqd + R = 6371000 + lat_0 = 51.0246 + lon_0 = -113.3981

Za nekaj primerov ploskov glejte Manipuliranje azimutnih ekvivalentnih projekcij v QGIS


Poglej si posnetek: FS 19 mapa Agrospol 1dil naživo (Oktober 2021).