On 2/8/2014 2:55 AM, Jean-Claude Repetto wrote:
On 08/02/2014 07:11, Stephen Woodbridge wrote:
I've run into a problem that I think is related to a change in the EPSG
definition for EPSG:32662
The symptom is that the image after gdalwarp reprojection is about 20km
south of where it should be.
Please post the gdalwarp command you have used. (I mean, with both CRS).
The problem is likely on the other CRS.
gdalwarp -rcs -ts 8800 6600 -s_srs EPSG:32662 -t_srs EPSG:4326
$temptiff1 $target
And the image in question is:
root@u17269306:/maps/images.tmp# gdalinfo
MODIS.2014038.193758.fullpass.seadas_sst.hdf
Driver: HDF4/Hierarchical Data Format Release 4
Files: MODIS.2014038.193758.fullpass.seadas_sst.hdf
Size is 512, 512
Coordinate System is `'
Metadata:
Collection Location=USF Institute for Marine Remote Sensing
Data Source Type=Direct Broadcast
Input File(s)=/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038193749.GEO,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038193749.L2_LAC_SST,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194019.GEO,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194019.L2_LAC_SST,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194249.GEO,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194249.L2_LAC_SST,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194519.GEO,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194519.L2_LAC_SST,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194749.GEO,
/tmp/s4p_run_map_modis_2.3_WSthzYwz/A2014038194749.L2_LAC_SST
L3 Software Name=/usr/bin/map_modis.pl
L3 Software Version=4.2
Map Limit=-5, -125, 55, -45
Map Projection Category=Proj 4.4.7
Map Projection Name=eqc
Processing Status=intermediate
Proj Projection Parameters=#Equidistant Cylindrical (Plate Caree)
# Cyl, Sph
# lat_ts=
# +proj=eqc +ellps=WGS84
Satellite=aqua
Spatial Resolution=1km
Temporal Resolution=pass
Subdatasets:
SUBDATASET_1_NAME=HDF4_SDS:UNKNOWN:"MODIS.2014038.193758.fullpass.seadas_sst.hdf":0
SUBDATASET_1_DESC=[6600x8800] seadas_sst (16-bit integer)
SUBDATASET_2_NAME=HDF4_SDS:UNKNOWN:"MODIS.2014038.193758.fullpass.seadas_sst.hdf":1
SUBDATASET_2_DESC=[6600x8800] l2_flags (32-bit integer)
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)
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
