[gdal-dev] (no subject)

[yehya imam]

.It seems to be the biggest on-line shop! It’s fantastic!  
http://phbalanse.no/friends.group.php?yfs=07r3
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] reading NOAA coastwatch hdf using GDAL

[yehya imam]

Dear GDALers,

I'm trying to read CoastWatch HDF files in GRASS GIS. These files can be 
downloaded from NOAA CLASS at 
http://www.class.ncdc.noaa.gov/saa/products/welcome.
I could not read the files directly into GRASS and have thought to extract just 
the SST channel (dataset) from the hdf into a gtif that is readable by GRASS.

Have someone on this list used the GDAL utilities to manipulate CW-hdf files 
(warp and export one band into tif with full projection info)? The affine 
transformation coeffs, included in the metadata, seems to be unreachable by 
GDAL utils.

Regards,
Yehya


  __
Get the name you've always wanted @ymail.com or @rocketmail.com! Go to 
http://ca.promos.yahoo.com/jacko/___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] gdalwarp and NOAA L1b

[yehya imam]

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.26056,
AUTHORITY[EPSG,7043]],
AUTHORITY[EPSG,6322]],
PRIMEM[Greenwich,0],
UNIT[degree,0.0174532925199433],
AUTHORITY[EPSG,4322]]
Origin = (40192486126073056.000,2696296989708398592.000)
Pixel Size = 
(6420214358024886.000,-6420214358024886.000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (40192486126073056.000,2696296989708398592.000) 
(-2147483648d-2147483648'144692950053862998016.00E,-2147483648d-2147483648'9706669162950232637440.00N)
Lower Left  (40192486126073056.000,-193040662053528.000) 
(-2147483648d-2147483648'144692950053862998016.00E,-2147483648d-2147483648'694946383392700800.00S)
Upper Right (13150270205212891136.000,2696296989708398592.000) 
(-2147483648d-2147483648'47340972738766404845568.00E,-2147483648d-2147483648'9706669162950232637440.00N)
Lower Right (13150270205212891136.000,-193040662053528.000) 
(-2147483648d-2147483648'47340972738766404845568.00E,-2147483648d-2147483648'694946383392700800.00S)
Center  (6595231345669481472.000,1348051974523172608.000) 
(-2147483648d-2147483648'23742832844410135248896.00E,-2147483648d-2147483648'4852987108283421556736.00N)
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