Even, Frank,
thanks for your answers. I have run into a problem with using
gdalwarp with geoloc transform.
In ticket #1895 there is an hdf file for which the HDF4 driver creates
GEOLOCATION metadata.
Trying gdalwarp on it gives an error, and the result transform is
incorrect. Am I doing something wrong, or is this a bug?
command is : gdalwarp -geoloc -t_srs EPSG:4326
HDF4_SDS:UNKNOWN:"A2006005182000.L2_LAC_SST.x.hdf":13 tmp1.tif
more info below:
$ gdalinfo HDF4_SDS:UNKNOWN:"A2006005182000.L2_LAC_SST.x.hdf":13
Driver: HDF4Image/HDF4 Dataset
Files: A2006005182000.L2_LAC_SST.x.hdf
Size is 389, 602
Coordinate System is `'
Metadata:
[...]
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=HDF4_SDS:UNKNOWN:"A2006005182000.L2_LAC_SST.x.hdf":11
Y_BAND=1
Y_DATASET=HDF4_SDS:UNKNOWN:"A2006005182000.L2_LAC_SST.x.hdf":12
Corner Coordinates:
[...]
$ gdalwarp -geoloc -t_srs EPSG:4326
HDF4_SDS:UNKNOWN:"A2006005182000.L2_LAC_SST.x.hdf":13 tmp1.tif
Creating output file that is 31119P x 5864L.
Processing input file HDF4_SDS:UNKNOWN:A2006005182000.L2_LAC_SST.x.hdf:13.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.
0...10...20...30...40...50...60...70...80...90...100 - done.
$ gdalinfo tmp1.tif
Driver: GTiff/GeoTIFF
Files: tmp1.tif
Size is 31119, 5864
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-90.198287963867188,90.300000000000011)
Pixel Size = (0.015398926739995,-0.015398926739995)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -90.1982880, 90.3000000) ( 90d11'53.84"W, 90d18' 0.00"N)
Lower Left ( -90.1982880, 0.0006936) ( 90d11'53.84"W, 0d 0' 2.50"N)
Upper Right ( 389.001, 90.300) (Invalid angle, 90d18' 0.00"N)
Lower Right ( 389.001, 0.001) (Invalid angle, 0d 0' 2.50"N)
Center ( 149.4013126, 45.1503468) (149d24' 4.73"E, 45d 9' 1.25"N)
Regards,
Etienne
On Tue, Feb 7, 2012 at 7:05 PM, Frank Warmerdam<[email protected]> wrote:
On Tue, Feb 7, 2012 at 7:40 AM, Etienne Tourigny
<[email protected]> wrote:
Hi all,
I would like to have information on the status of geolocation array
support in GDAL.
The relevant RFC4 (http://trac.osgeo.org/gdal/wiki/rfc4_geolocate) is
still "under development", and the information there is that
GDALCreateGeoLocTransformer "is currently partially implemented".
Etienne,
I think the geoloc transformer received some fixes from someone
other than myself and now is working reasonable well. It would
be desirable to freshen up the RFC and get it passed but it isn't
particularly a priority of mine so I'm not likely to drive it in the
very near future.
Best regards,
--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, [email protected]
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Software Developer
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
