Viac

Klin chýbajúcich údajov sa replikuje do stereografických snímok

Klin chýbajúcich údajov sa replikuje do stereografických snímok


Snažím sa pomocou gdalwarpu reprojektovať ekvirektangulárny projekčný obraz na stereografický okolo severného pólu Marsu a nie som si istý, čo sa deje v mojom výstupnom obrázku.

Vstupný obrázok je normálne vyzerajúci obdĺžnikový obrázok bez čiernych oblastí. Pri priblížení k chýbajúcej klinovej oblasti vidím, že nie všetky pixely zo vstupného obrázku sú prítomné vo výstupnom obrázku.

Môj príkaz gdalwarp:

gdalwarp map_upper_crop.tif map_north_pole.tif -t_srs " +proj = stere +lon_0 = 0 +lat_0 = 90 +a = 3396200 +b = 3376200 jednotiek = m"

výstup gdalinfo na vstupnom obrázku:

Ovládač: GTiff/GeoTIFF
Súbory: projection_space_basemap_top.tif
Veľkosť je 11520, 1035
Súradnicový systém je:

PROJCS ["nemenovaný", GEOGCS ["nemenovaný elipsa", DATUM ["neznámy", SPHEROID ["nemenovaný", 3396200,169,8100000000004]], PRIMEM ["Greenwich", 0], UNIT ["stupeň", 0,0174532925199433]], PROJECTION ["Equirectangular"], PARAMETER ["latitude_of_origin", 0], PARAMETER ["central_meridian", 0], PARAMETER ["standard_parallel_1", 0], PARAMETER ["false_easting", 0], PARAMETER ["false_northing", 0 ], UNIT ["meter", 1, AUTHORITY ["EPSG", "9001"]]]

Reťazec PROJ.4 je:
„ +proj = eqc +lat_ts = 0 +lat_0 = 0 +lon_0 = 0 +x_0 = 0 +y_0 = 0 +a = 3396200 +b = 3376200 +jednotiek = m +no_defs '
Pôvod = (-10669477.099999999627471,5334738,599999999627471)
Veľkosť pixelov = (1852,339774305555466, -1852,339791666666542)
Metadáta:
AREA_OR_POINT = Oblasť
TIFFTAG_RESOLUTIONUNIT = 2 (pixely/palec)
TIFFTAG_XRESOLUTION = 1024
TIFFTAG_YRESOLUTION = 1024
Metadáta štruktúry obrázku:
INTERLEAVE = PIXEL
Súradnice rohu:
Vľavo hore (-10669477.100, 5334738.600) (179d59'59,99 "E, 90d 0 '0,01" N)
Vľavo dole (-10669477.100, 3417566.916) (179d59'59,99 "E, 57d39'22,50" N)
Vpravo hore (10669477.100, 5334738.600) (179d59'59,99 "W, 90d 0 '0,01" N)
Vpravo dole (10669477.100, 3417566.916) (179d59'59,99 "W, 57d39'22,50" N)
Stred (0,000, 4376152,758) (0d 0 '0,01 "E, 73d49'41,26" N)

Máte predstavu, čo sa tu deje?


Podarilo sa mi získať správny výstup pridaním nasledujúcich možností do príkazu gdalwarp:

-wo sample_grid = yes, -wo sample_steps = 1000, -wo source_extra = 2 -wo CUTLINE_ALL_TOUCHED = true -r cubic

Ďakujeme používateľovi30184 za upozornenie!