Suite

Cale de reprojection des données manquantes en stéréographie


J'essaie de reprojeter une image de projection équirectangulaire en stéréographie autour du pôle nord de mars à l'aide de gdalwarp et je ne suis pas sûr de ce qui se passe dans mon image de sortie.

L'image d'entrée est une image rectangulaire d'aspect normal sans zones noires. En zoomant sur la région de coin manquante, je peux voir que tous les pixels de l'image d'entrée ne sont pas présents dans l'image de sortie.

Ma commande gdalwarp :

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

sortie gdalinfo sur l'image d'entrée :

Pilote : GTiff/GeoTIFF
Fichiers : projection_space_basemap_top.tif
La taille est 11520, 1035
Le système de coordonnées est :

PROJCS["sans nom", GEOGCS["ellipse sans nom", DATUM["inconnu", SPHEROID["sans nom",3396200,169.8100000000004]], PRIMEM["Greenwich",0], UNIT["degré",0.0174532925199433]], PROJECTION["Équirectangulaire"], PARAMÈTRE["latitude_d'origine",0], PARAMÈTRE["central_meridian",0], PARAMÈTRE["standard_parallel_1",0], PARAMÈTRE["false_easting",0], PARAMÈTRE["false_northing",0 ], UNITÉ["mètre",1, AUTORITÉ["EPSG","9001"]]]

La chaîne PROJ.4 est :
'+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=3396200 +b=3376200 +units=m +no_defs '
Origine = (-10669477.099999999627471,5334738,5999999999627471)
Taille de pixel = (1852.339774305555466,-1852.339791666666542)
Métadonnées :
AREA_OR_POINT=Zone
TIFFTAG_RESOLUTIONUNIT=2 (pixels/pouce)
TIFFTAG_XRESOLUTION=1024
TIFFTAG_YRESOLUTION=1024
Métadonnées de la structure de l'image :
INTERLEAVE=PIXEL
Coordonnées du coin :
En haut à gauche (-10669477.100, 5334738.600) (179d59'59.99"E, 90d 0' 0.01"N)
En bas à gauche (-10669477.100, 3417566.916) (179d59'59.99"E, 57d39'22.50"N)
En haut à droite (10669477.100, 5334738.600) (179d59'59.99"W, 90d 0' 0.01"N)
Inférieur droit (10669477.100, 3417566.916) (179d59'59.99"W, 57d39'22.50"N)
Centre (0,000, 4376152,758) ( 0d 0' 0,01"E, 73d49'41,26"N)

Une idée de ce qui ne va pas ici ?


J'ai pu obtenir la sortie correcte en ajoutant les options suivantes à ma commande gdalwarp :

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

Merci user30184 pour les avoir signalés !