Becker, Please take a look at RFC#4: http://trac.osgeo.org/gdal/wiki/rfc4_geolocate
Regards, Ivan > -------Original Message------- > From: Becker, Thomas <[email protected]> > Subject: [gdal-dev] reprojection HDF / GCP extraction > Sent: Aug 10 '09 07:00 > > Hi all, > > I have a couple NOAA-15 images as HDF-EOS files which I would like to > reproject to WGS84. The data carry a subgroup called 'Geolocation Fields' holding latitude and longitude for each pixel in separated arrays. Reading the files with gdalinfo leads to the output at the end of this posting. > > Out of the HDF-EOS format I would like to create a correct georeferenced > image (GTIF). One which will not look like a swath scanned image, but a skewed and correctly referenced one. > > Using gdalwarp leads to the message: There is no affine transformation and > no GCPs. > > Since the coordinates for the centre of each pixel are stored in the before > mentioned arrays, is there a way to extract GCPs out of the data. > > Thanks for your help. > > Best regards, > Thomas > > gdalinfo output: > > Driver: HDF4/Hierarchical Data Format Release 4 > Files: amsua15_2000.034_08964_0104_0258_GC.eos > Size is 512, 512 > Coordinate System is `' > Metadata: > HDFEOSVersion=HDFEOS_V2.4 > Subdatasets: > > SUBDATASET_1_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:23800.37 MHz > SUBDATASET_1_DESC=[851x30] 23800.37 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_2_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:31400.42 MHz > SUBDATASET_2_DESC=[851x30] 31400.42 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_3_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:50299.91 MHz > SUBDATASET_3_DESC=[851x30] 50299.91 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_4_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:52799.39 MHz > SUBDATASET_4_DESC=[851x30] 52799.39 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_5_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:53595.41 +- 115 MHz > SUBDATASET_5_DESC=[851x30] 53595.41 +- 115 MHz Orbit 8964 (32-bit > floating-point) > > SUBDATASET_6_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:54399.53 MHz > SUBDATASET_6_DESC=[851x30] 54399.53 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_7_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:54940.64 MHz > SUBDATASET_7_DESC=[851x30] 54940.64 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_8_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:55498.70 MHz > SUBDATASET_8_DESC=[851x30] 55498.70 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_9_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:57290.33 MHz > SUBDATASET_9_DESC=[851x30] 57290.33 MHz Orbit 8964 (32-bit floating-point) > > SUBDATASET_10_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:57290.33 +- 217 MHz > SUBDATASET_10_DESC=[851x30] 57290.33 +- 217 MHz Orbit 8964 (32-bit > floating-point) > > SUBDATASET_11_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:57290.33 +- 322.2 +- 48 MHz > SUBDATASET_11_DESC=[851x30] 57290.33 +- 322.2 +- 48 MHz Orbit 8964 (32-bit > floating-point) > > SUBDATASET_12_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:57290.33 +- 322.2 +- 22 MHz > SUBDATASET_12_DESC=[851x30] 57290.33 +- 322.2 +- 22 MHz Orbit 8964 (32-bit > floating-point) > > SUBDATASET_13_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:57290.33 +- 322.2 +- 10 MHz > SUBDATASET_13_DESC=[851x30] 57290.33 +- 322.2 +- 10 MHz Orbit 8964 (32-bit > floating-point) > > SUBDATASET_14_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:57290.33 +- 322.2 +- 4.5 MHz > SUBDATASET_14_DESC=[851x30] 57290.33 +- 322.2 +- 4.5 MHz Orbit 8964 > (32-bit floating-point) > > SUBDATASET_15_NAME=HDF4_EOS:EOS_SWATH:"amsua15_2000.034_08964_0104_0258_GC.eos":Orbit > 8964:88997.00 MHz > SUBDATASET_15_DESC=[851x30] 88997.00 MHz Orbit 8964 (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) > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev > _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
