[Gdal-dev] new bug in FAST format

[Fu Chen]

Hi!
I test one of my old program. But find GDAL120b has some problems in parsing the FAST format.

#include "gdal_priv.h"
#include "cpl_conv.h"

int main ( int argc, char ** argv)
{
    GDALAllRegister();
	GDALDataset * poDataset= (GDALDataset *) GDALOpen( "h:\\chinagisdata\\landsatdata\\L71118038_03820020111_HPN.FST", GA_ReadOnly );
	const char * sDriverDesc = poDataset->GetProjectionRef();
	printf(sDriverDesc);
}

gdal119 output:

PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.
257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","632
6"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519
9433,AUTHORITY["EPSG","9108"]],AXIS["Lat",NORTH],AXIS["Long",EAST],AUTHORITY["EP
SG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0]
,PARAMETER["central_meridian",123],PARAMETER["scale_factor",1],PARAMETER["false_
easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]

gdal120b output:

PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.
257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","632
6"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519
9433,AUTHORITY["EPSG","9108"]],AXIS["Lat",NORTH],AXIS["Long",EAST],AUTHORITY["EP
SG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0]
,PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_ea
sting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

You can see in central_meridian false_easting are not correctly interpreted in 120b. What's change in the new version on the FAST format? 

any idea?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: L71118038_03820020111_HPN.FST
Type: application/octet-stream
Size: 4608 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/gdal-dev/attachments/20040316/7119dcb4/L71118038_03820020111_HPN.obj