[gdal-dev] gdalwarp and NOAA L1b

yehya imam yehya_imam at yahoo.com
Thu May 14 15:21:11 EDT 2009


Dear GDAL'ers,

I've downloaded AVHRR L1b images from NOAA CLASS. I'm trying to rectify the images via gdalwarp -tps. The output is obviously erroneous as per gdalinfo (also tried to open image in GRASS). I would appreciate your help in resolving this issue. Alternatively, could you please recommend another public tool to perform the task?

>> Here is the gdalinfo -nogcp for the original image:
***********
Driver: L1B/NOAA Polar Orbiter Level 1b Data Set
Files: NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC
Size is 2048, 191
Coordinate System is `'
Metadata:
  DATASET_NAME=NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC
  SATELLITE=NOAA-12(D)
  DATA_TYPE=AVHRR LAC
  REVOLUTION=17013
  SOURCE=Fairbanks, Alaska, USA (formerly Gilmore Creek)
  PROCESSING_CENTER=NOAA/NESDIS - Suitland, Maryland, USA
  START=year: 1994, day: 235, millisecond: 58454903
  STOP=year: 1994, day: 235, millisecond: 58494236
  LOCATION=Descending
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  191.0)
Upper Right ( 2048.0,    0.0)
Lower Right ( 2048.0,  191.0)
Center      ( 1024.0,   95.5)
Band 1 Block=2048x1 Type=UInt16, ColorInterp=Undefined
  Description = AVHRR Channel 1:  0.58  micrometers -- 0.68 micrometers
Band 2 Block=2048x1 Type=UInt16, ColorInterp=Undefined
  Description = AVHRR Channel 2:  0.725 micrometers -- 1.10 micrometers
Band 3 Block=2048x1 Type=UInt16, ColorInterp=Undefined
  Description = AVHRR Channel 3:  3.55  micrometers -- 3.93 micrometers
Band 4 Block=2048x1 Type=UInt16, ColorInterp=Undefined
  Description = AVHRR Channel 4:  10.3  micrometers -- 11.3 micrometers
Band 5 Block=2048x1 Type=UInt16, ColorInterp=Undefined
  Description = AVHRR Channel 5:  11.5  micrometers -- 12.5 micrometers
*************

>> A sample of GCPs returned by gdalinfo with no flags
*************
GCP[  0]: Id=, Info=
          (25.5,0.5) -> (-142.140625,54.984375,0)
GCP[  1]: Id=, Info=
          (225.5,0.5) -> (-134.1015625,54.8984375,0)
GCP[  2]: Id=, Info=
          (105.5,0.5) -> (-138.1640625,55.0078125,0)
GCP[  3]: Id=, Info=
          (625.5,0.5) -> (-126.4375,54.3125,0)
GCP[  4]: Id=, Info=
          (185.5,0.5) -> (-135.28125,54.9453125,0)
GCP[  5]: Id=, Info=
          (225.5,0.5) -> (-134.1015625,54.8984375,0)
GCP[  6]: Id=, Info=
          (1225.5,0.5) -> (-119.21875,53.28125,0)
GCP[  7]: Id=, Info=
          (305.5,0.5) -> (-132.0703125,54.7890625,0)
GCP[  8]: Id=, Info=
          (345.5,0.5) -> (-131.1875,54.734375,0)
GCP[  9]: Id=, Info=
          (385.5,0.5) -> (-130.375,54.671875,0)
GCP[ 10]: Id=, Info=
          (2025.5,0.5) -> (-103.5546875,49.1796875,0)
**************


>> Issuing gdalwarp -tps NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC sample4.tif
**************
There is a problem to invert the interpolation matrix
 There is a problem to invert the interpolation matrix
Creating output file that is 2042P x 420L.
Processing input file NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC.
 There is a problem to invert the interpolation matrix
 There is a problem to invert the interpolation matrix
0...10...20...30...40...50...60...70...80...90...100 - done.
***************

>> gdalinfo on sampl4.tif returns
***************
Driver: GTiff/GeoTIFF
Files: sample4.tif
Size is 2042, 420
Coordinate System is:
GEOGCS["WGS 72",
    DATUM["WGS_1972",
        SPHEROID["WGS 72",6378135,298.2600000000056,
            AUTHORITY["EPSG","7043"]],
        AUTHORITY["EPSG","6322"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4322"]]
Origin = (40192486126073056.000000000000000,2696296989708398592.000000000000000)
Pixel Size = (6420214358024886.000000000000000,-6420214358024886.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (40192486126073056.000,2696296989708398592.000) (-2147483648d-2147483648'144692950053862998016.00"E,-2147483648d-2147483648'9706669162950232637440.00"N)
Lower Left  (40192486126073056.000,-193040662053528.000) (-2147483648d-2147483648'144692950053862998016.00"E,-2147483648d-2147483648'694946383392700800.00"S)
Upper Right (13150270205212891136.000,2696296989708398592.000) (-2147483648d-2147483648'47340972738766404845568.00"E,-2147483648d-2147483648'9706669162950232637440.00"N)
Lower Right (13150270205212891136.000,-193040662053528.000) (-2147483648d-2147483648'47340972738766404845568.00"E,-2147483648d-2147483648'694946383392700800.00"S)
Center      (6595231345669481472.000,1348051974523172608.000) (-2147483648d-2147483648'23742832844410135248896.00"E,-2147483648d-2147483648'4852987108283421556736.00"N)
Band 1 Block=2042x1 Type=UInt16, ColorInterp=Gray
Band 2 Block=2042x1 Type=UInt16, ColorInterp=Undefined
Band 3 Block=2042x1 Type=UInt16, ColorInterp=Undefined
Band 4 Block=2042x1 Type=UInt16, ColorInterp=Undefined
Band 5 Block=2042x1 Type=UInt16, ColorInterp=Undefined
******************

Thanks,
Yehya



      __________________________________________________________________
Get a sneak peak at messages with a handy reading pane with All new Yahoo! Mail: http://ca.promos.yahoo.com/newmail/overview2/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090514/e9d9060e/attachment.html


More information about the gdal-dev mailing list