Veel

Gdalwarp jätab 18026 pikkuselt suure augu, kui taandub 4326-lt kohandatud AEqD-le

Gdalwarp jätab 18026 pikkuselt suure augu, kui taandub 4326-lt kohandatud AEqD-le


Püüan loodusliku Maa I 1: 50m rastri ümber projekteerida asimutaalseks võrdse kaugusega projektsiooniks, mille keskpunkt on minu kodulinn. Nagu te võib-olla teate, katab lähteraster vahemikus -180 ° kuni 180 ° laiust ja -90 ° kuni 90 °. Siin on käsk, mida kasutan:

gdalwarp -multi -dstnodata 0 -baidine -t_srs "+ proj = aeqd + lat_0 = 48.037 + lon_0 = -122.408 + ellips = sfäär + ühikud = m + no_defs" -r kuup-GTiff -ts 10800 10800 NE1_50M_SR_W.tif natural_earth .tif

See töötab üsna hästi, välja arvatud see, et see jätab 180 ° läände suure koleda tühimiku. Samuti lõikab see osa väljundrastri mõlemast küljest.

Olen katsetanud raadiuse seadistamist+ ellipsidparameeter, ilma muutusteta. WGS84 kasutamine annab väljundina vaid väikese osa Maast.

Siin on, mida ma proovides sain+ R_Aselle asemel+ ellps = kera:

Tulemus WGS84-ga oli paljuski sama. Enne minu kohandatud AEqD-le uuesti proovimist proovisin ka vaheprojektsioonina sfäärilisse Mercatorisse ümberprojekteerida, ilmuvad samad esemed.

Niisiis, kuidas ma tühimiku täidan? Pikslid on sisendrastris kindlasti olemas, nad lihtsalt ei jõua väljundisse.


Kuupäevarea lõhe vabanesin nende võimalustega:

gdal_translate -a_srs "+ proj = longlat + ellps = sfäär" NE1_50M_SR_W.tif ne.tif gdalwarp -dstnodata 0 -bait -wo SOURCE_EXTRA = 1000 -t_srs "+ projekt = aeqd + lat_0 = 48.037 + lon_0 = -122.408 + sfäär + ühikud = m + no_defs "-of GTiff ne.tif neaeqd.tif

TheSOURCE_EXTRAkirjutamisvalik tagab kattuvuse:

SOURCE_EXTRA: See on arv lisapiksleid, mis on antud päringu jaoks lisatud lähteakna ümber ja ümardamisvea eest hoolitsemine on vaikimisi 1. Selle suuruse määramine suurendab lugemist vajavate andmete hulka, kuid võib vältida lähteandmete puudumist.

Ikka saan atolerantsi tingimuse viga, mis võib tulla antipodaalsest punktist, lõpmatuseni ümber lükatud. Külgede jaotustükid jäävad: