[gdal-dev] FWTools 2.3.0 does find State Plane projection information

Charles Karney ckarney at sarnoff.com
Mon May 18 16:04:23 EDT 2009


FWTools 2.2.8 and 2.3.0 seem not to know about some state plane
coordinate systems, while gdal under Linux and Cygwin report the correct
projection data.

For example, I have some imagery (a GeoTiff) for Florida downloaded from
the USGS site.  The cygwin version of listgeo reports:

  Geotiff_Information:
     Version: 1
     Key_Revision: 1.0
     Tagged_Information:
	ModelTiepointTag (2,3):
	   0                0                0                
	   419093.218       185211.481       0                
	ModelPixelScaleTag (1,3):
	   0.30480061       0.30480061       0                
	End_Of_Tags.
     Keyed_Information:
	GTModelTypeGeoKey (Short,1): ModelTypeProjected
	GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
	GTCitationGeoKey (Ascii,271): "IMAGINE GeoTIFF Support\n
	 Projection Name = NAD_1983_StatePlane_Florida_North_FIPS_0903_Feet\n
	 Units = us_survey_feet\n
	 GeoTIFF Units = meters"
	ProjectedCSTypeGeoKey (Short,1): PCS_NAD83_Florida_North
	PCSCitationGeoKey (Ascii,189): "IMAGINE GeoTIFF Support\n
	 State Plane Zone -903\n
	 NAD = 83"
	End_Of_Keys.
     End_Of_Geotiff.

  PCS = 26960 (NAD83 / Florida North)
  Projection = 10933 (SPCS83 Florida North zone (meters))
  Projection Method: CT_LambertConfConic_2SP
     ProjFalseOriginLatGeoKey: 29.000000 ( 29d 0' 0.00"N)
     ProjFalseOriginLongGeoKey: -84.500000 ( 84d30' 0.00"W)
     ProjStdParallel1GeoKey: 30.750000 ( 30d45' 0.00"N)
     ProjStdParallel2GeoKey: 29.583333 ( 29d35' 0.00"N)
     ProjFalseEastingGeoKey: 600000.000000 m
     ProjFalseNorthingGeoKey: 0.000000 m
  GCS: 4269/NAD83
  Datum: 6269/North American Datum 1983
  Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
  Prime Meridian: 8901/Greenwich (0.000000/  0d 0' 0.00"E)
  Projection Linear Units: 9001/metre (1.000000m)

  Corner Coordinates:
  Upper Left    (  419093.218,  185211.481)  ( 86d23'15.45"W, 30d39'26.37"N)
  Lower Left    (  419093.218,  183070.562)  ( 86d23'14.12"W, 30d38'16.86"N)
  Upper Right   (  420296.571,  185211.481)  ( 86d22'30.26"W, 30d39'27.02"N)
  Lower Right   (  420296.571,  183070.562)  ( 86d22'28.94"W, 30d38'17.50"N)
  Center        (  419694.894,  184141.021)  ( 86d22'52.19"W, 30d38'51.94"N)

The version of lisgeo in FWTools2.3.0 reports the same information up to
the PCS line.  Then it has

  PCS = 26960 (name unknown)
  Projection = 10933 ()

and the corner points are listed without the lat/long equivalents.

I see that there's a stateplane.csv in the FWTools distribution.  But it
appears not to be read?  (I can provide the GeoTiff, if necessary.)

Thanks.

-- 
Charles Karney <ckarney at sarnoff.com>
Sarnoff Corporation, Princeton, NJ 08543-5300

Tel: +1 609 734 2312
Fax: +1 609 734 2662


More information about the gdal-dev mailing list