J'essaie de remapper de Geostationary à Lambert conforme en utilisant gdalwarp. Mes données d'entrée sont en netcdf, et sont en coordonnées géographiques (degrés) et je voudrais sortir les données remappées vers netcdf. J'ai créé un fichier vrt correspondant pour les données d'entrée netcdf. Gdalwarp affichera le fichier netcdf, mais les données de sortie sont toutes des zéros et je reçois l'erreur suivante:
Creating output file that is 5120P x 5120L.
Processing input file netcdf.vrt.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.
Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping.
0...10...20...30...40...50...60...70...80...90...100 - done.J'ai tenté la commande suivante:
/usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 +lon_0=-75 +x_0=-0.151844 +y_0=0.151844 +a=6378140 +b=6356754.99999591 +units=degrees +no_defs" -t_srs "+proj=lcc +ellps=clrk66 +a=6378137 +b=6378137 +e=0.0818191910435 +lat_0=24.9999 +lon_0=-95 +lat_1=24.9999 +lat_ts=25.0001 +units=meters +no_defs" -te -1952976.3246 -828316.5944 3248431.6754 4373091.4056 -of netCDF -geoloc -overwrite -r bilinear -ts 5120 5120 netcdf.vrt out.ncEst-ce que gdalwarp peut repapper des coordonnées géographiques en projeté? Ou dois-je d'abord traduire géographique par projeté? En outre, gdalwarp peut-il lire les informations de projection directement à partir de netcdf ou avez-vous BESOIN d'écrire d'abord en .vrt?
Voici ce que gdalinfo sort du fichier d'entrée: (c'est un fichier GOES 13 de CLASS)
Size is 512, 512
Coordinate System is `'
Metadata:
NC_GLOBAL#Conventions=CF-1.4
NC_GLOBAL#Satellite Sensor=G-13 IMG
NC_GLOBAL#Source=McIDAS Area File
Subdatasets:
SUBDATASET_1_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":auditTrail
SUBDATASET_1_DESC=[3x80] auditTrail (8-bit character)
SUBDATASET_2_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":data
SUBDATASET_2_DESC=[1x665x2036] data (32-bit floating-point)
SUBDATASET_3_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lat
SUBDATASET_3_DESC=[665x2036] lat (32-bit floating-point)
SUBDATASET_4_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lon
SUBDATASET_4_DESC=[665x2036] lon (32-bit floating-point)
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 512.0)
Upper Right ( 512.0, 0.0)
Lower Right ( 512.0, 512.0)
Center ( 256.0, 256.0)Et des informations supplémentaires sur la variable de données:
Driver: netCDF/Network Common Data Format
Files: goes13.2013.100.174518.BAND_04.nc
Size is 2036, 665
Coordinate System is `'
Metadata:
data#coordinates=lon lat
data#long_name=0-255 Brightness Temperature
data#type=VISR
NC_GLOBAL#Conventions=CF-1.4
NC_GLOBAL#Satellite Sensor=G-13 IMG
NC_GLOBAL#Source=McIDAS Area File
NETCDF_DIM_EXTRA={time}
NETCDF_DIM_time_DEF={1,4}
NETCDF_DIM_time_VALUES=1365615900
time#long_name=seconds since 1970-1-1 0:0:0
time#units=seconds since 1970-1-1 0:0:0
Geolocation:
LINE_OFFSET=0
LINE_STEP=1
PIXEL_OFFSET=0
PIXEL_STEP=1
SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]
X_BAND=1
X_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lon
Y_BAND=1
Y_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lat
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 665.0)
Upper Right ( 2036.0, 0.0)
Lower Right ( 2036.0, 665.0)
Center ( 1018.0, 332.5)
Band 1 Block=2036x1 Type=Float32, ColorInterp=Undefined
NoData Value=9.96920996838686905e+36
Metadata:
coordinates=lon lat
long_name=0-255 Brightness Temperature
NETCDF_DIM_time=1365615900
NETCDF_VARNAME=data
type=VISRToute aide sera fortement appréciée!
1
La projection géographique n'utilisera pas de degrés; essayez les compteurs. Où obtenez-vous les valeurs + x_0 / + y_0? Basé sur gdalinfo, je ne suis pas sûr que le raster en entrée soit géoréférencé. Dans les srs cibles, vous avez + a = + b qui est une sphère, mais aussi set + e. Cependant, + ellps est pour un ellipsoïde entièrement différent. Les différentes valeurs + lat semblent également toutes étranges. lat_ts est lat de l'échelle réelle, donc un parallèle standard comme lat_1.
—
mkennedy
Je vous remercie. Je vais essayer d'utiliser des compteurs. J'obtiens les x_0 et y_0 (échelle et décalages) de la définition GOES bien que ce ne soient pas des entrées requises pour + proj = geos, donc je peux essayer de les supprimer. Et merci d'avoir souligné l'ajout d'ellipsoïdes + e Les définitions lat pour les t_srs sont pour la définition AWIPS de lambert (une taille de sortie spécifique). J'ajouterai ce que gdal info recrache pour la variable de données spécifique au message de question (trop long pour un commentaire)
—
Katie J
La définition AWIPS à laquelle je fais référence est décrite sur cette page: nws.noaa.gov/noaaport/html/icdtb48_2.html (la première est le Lambert que j'essaie de remapper )
—
Katie J
Hmmm, donc il dit lat / lon WGS84, mais les coordonnées de coin signalées me concernent parce que ce ne sont que des valeurs de cellules brutes. Le LCC est un cas tangent - un seul parallèle standard / latitude d'origine sont tous à 25N. Je n'ai travaillé avec aucune de ces données, donc je vais juste par les informations de métadonnées.
—
mkennedy
L'image n'est pas géoréférencée, mais une source srs est fournie. Quelques questions: * Pouvez-vous exécuter avec CPL_DEBUG = GDAL_netCDF? Donc CPL_DEBUG = GDAL_netCDF / usr / bin / gdalwarp ... Je soupçonne qu'il peut y avoir un problème avec les tableaux de géolocalisation. * Pouvez-vous rendre vos données disponibles?