Re: [gdal-dev] reprojection HDF / GCP extraction

[Lucena, Ivan]

Becker,

Please take a look at RFC#4: http://trac.osgeo.org/gdal/wiki/rfc4_geolocate

Regards,

Ivan


>  -------Original Message-------
>  From: Becker, Thomas <thob at dmu.dk>
>  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
>  gdal-dev at lists.osgeo.org
>  http://lists.osgeo.org/mailman/listinfo/gdal-dev
>