<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:10pt;color:#111111;"><div>Dear GDAL'ers,<br><br>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?<br><br>>> Here is the gdalinfo -nogcp for the original image:<br>***********<br>Driver: L1B/NOAA Polar Orbiter Level 1b Data Set<br>Files: NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC<br>Size is 2048, 191<br>Coordinate System is `'<br>Metadata:<br> DATASET_NAME=NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC<br> SATELLITE=NOAA-12(D)<br> DATA_TYPE=AVHRR LAC<br> REVOLUTION=17013<br> SOURCE=Fairbanks, Alaska,
USA (formerly Gilmore Creek)<br> PROCESSING_CENTER=NOAA/NESDIS - Suitland, Maryland, USA<br> START=year: 1994, day: 235, millisecond: 58454903<br> STOP=year: 1994, day: 235, millisecond: 58494236<br> LOCATION=Descending<br>Corner Coordinates:<br>Upper Left ( 0.0, 0.0)<br>Lower Left ( 0.0, 191.0)<br>Upper Right ( 2048.0, 0.0)<br>Lower Right ( 2048.0, 191.0)<br>Center ( 1024.0, 95.5)<br>Band 1 Block=2048x1 Type=UInt16, ColorInterp=Undefined<br> Description = AVHRR Channel 1: 0.58 micrometers -- 0.68 micrometers<br>Band 2 Block=2048x1 Type=UInt16, ColorInterp=Undefined<br> Description = AVHRR Channel 2: 0.725 micrometers -- 1.10 micrometers<br>Band 3 Block=2048x1 Type=UInt16, ColorInterp=Undefined<br> Description = AVHRR Channel 3: 3.55 micrometers
-- 3.93 micrometers<br>Band 4 Block=2048x1 Type=UInt16, ColorInterp=Undefined<br> Description = AVHRR Channel 4: 10.3 micrometers -- 11.3 micrometers<br>Band 5 Block=2048x1 Type=UInt16, ColorInterp=Undefined<br> Description = AVHRR Channel 5: 11.5 micrometers -- 12.5 micrometers<br>*************<br><br>>> A sample of GCPs returned by gdalinfo with no flags<br>*************<br>GCP[ 0]: Id=, Info=<br> (25.5,0.5) -> (-142.140625,54.984375,0)<br>GCP[ 1]: Id=, Info=<br> (225.5,0.5) -> (-134.1015625,54.8984375,0)<br>GCP[ 2]: Id=, Info=<br> (105.5,0.5) -> (-138.1640625,55.0078125,0)<br>GCP[ 3]: Id=, Info=<br> (625.5,0.5) -> (-126.4375,54.3125,0)<br>GCP[ 4]: Id=,
Info=<br> (185.5,0.5) -> (-135.28125,54.9453125,0)<br>GCP[ 5]: Id=, Info=<br> (225.5,0.5) -> (-134.1015625,54.8984375,0)<br>GCP[ 6]: Id=, Info=<br> (1225.5,0.5) -> (-119.21875,53.28125,0)<br>GCP[ 7]: Id=, Info=<br> (305.5,0.5) -> (-132.0703125,54.7890625,0)<br>GCP[ 8]: Id=, Info=<br> (345.5,0.5) -> (-131.1875,54.734375,0)<br>GCP[ 9]: Id=, Info=<br> (385.5,0.5) -> (-130.375,54.671875,0)<br>GCP[ 10]: Id=, Info=<br> (2025.5,0.5) -> (-103.5546875,49.1796875,0)<br>**************<br><br><br>>> Issuing gdalwarp -tps
NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC sample4.tif<br>**************<br> There is a problem to invert the interpolation matrix<br> There is a problem to invert the interpolation matrix<br>Creating output file that is 2042P x 420L.<br>Processing input file NSS.LHRR.ND.D94235.S1612.E1622.B1701313.GC.<br> There is a problem to invert the interpolation matrix<br> There is a problem to invert the interpolation matrix<br>0...10...20...30...40...50...60...70...80...90...100 - done.<br>***************<br><br>>> gdalinfo on sampl4.tif returns<br>***************<br>Driver: GTiff/GeoTIFF<br>Files: sample4.tif<br>Size is 2042, 420<br>Coordinate System is:<br>GEOGCS["WGS 72",<br> DATUM["WGS_1972",<br> SPHEROID["WGS 72",6378135,298.2600000000056,<br>
AUTHORITY["EPSG","7043"]],<br> AUTHORITY["EPSG","6322"]],<br> PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4322"]]<br>Origin = (40192486126073056.000000000000000,2696296989708398592.000000000000000)<br>Pixel Size = (6420214358024886.000000000000000,-6420214358024886.000000000000000)<br>Metadata:<br> AREA_OR_POINT=Area<br>Image Structure Metadata:<br> INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left (40192486126073056.000,2696296989708398592.000) (-2147483648d-2147483648'144692950053862998016.00"E,-2147483648d-2147483648'9706669162950232637440.00"N)<br>Lower Left (40192486126073056.000,-193040662053528.000) (-2147483648d-2147483648'144692950053862998016.00"E,-2147483648d-2147483648'694946383392700800.00"S)<br>Upper Right (13150270205212891136.000,2696296989708398592.000)
(-2147483648d-2147483648'47340972738766404845568.00"E,-2147483648d-2147483648'9706669162950232637440.00"N)<br>Lower Right (13150270205212891136.000,-193040662053528.000) (-2147483648d-2147483648'47340972738766404845568.00"E,-2147483648d-2147483648'694946383392700800.00"S)<br>Center (6595231345669481472.000,1348051974523172608.000) (-2147483648d-2147483648'23742832844410135248896.00"E,-2147483648d-2147483648'4852987108283421556736.00"N)<br>Band 1 Block=2042x1 Type=UInt16, ColorInterp=Gray<br>Band 2 Block=2042x1 Type=UInt16, ColorInterp=Undefined<br>Band 3 Block=2042x1 Type=UInt16, ColorInterp=Undefined<br>Band 4 Block=2042x1 Type=UInt16, ColorInterp=Undefined<br>Band 5 Block=2042x1 Type=UInt16, ColorInterp=Undefined<br>******************<br><br>Thanks,<br>Yehya<br></div></div><br>
<p class="MsoNormal"> </p>
<tbody><tr>
<td style="padding: 0.75pt;">
<div class="MsoNormal" style="text-align: center;" align="center"><font face="Times New Roman" size="3"><span style="font-size: 12pt;">
<hr align="center" size="1" width="100%">
</span></font></div>
<p class="MsoNormal"><font face="Times New Roman" size="3"><span style="font-size: 12pt;"><img id="_x0000_i1026" src="http://us.i1.yimg.com/us.yimg.com/i/ca/iotg_search.jpg" align="absbottom" border="0" height="25" hspace="4" width="25"><a href="http://ca.toolbar.yahoo.com/" target="_new"><b><span style="font-weight: bold;" lang="NO-BOK">Yahoo!
Canada Toolbar :</span></b><span lang="NO-BOK"> Search from anywhere on
the web and bookmark your favourite sites. Download it now! </span></a>
</span></font><span lang="NO-BOK"><o:p></o:p></span></p></body></html>