Trouble with Projection

[Frank Warmerdam]

Jim Shanton wrote:
> I’m trying to plot a property parcel shp file, and it keeps plotting in
> Pennsylvania instead of North Carolina.
...
>   PROJECTION
>
>      proj=lcc
>
>      lat_0=33.75              #Latitude of False Origin
>
>      lon_0=-79.0              #Longitude of Natural Origin (Central
> Meridian)
>
>      lat_1=34.33333333333334  #Latitude of First Standard Parallel
>
>      lat_2=36.16666666666666
>
>      x_0=2000000.002616666    # False Easting
>
>      y_0=0.0                  # False Northing
>
>      ellps=GRS80              # SPHEROID
>
>      a=6378137.0              # from SPHEROID
> ["GRS_1980",6378137,298.257222101]],
>
>      rf=298.257222101         #from SPHEROID
> ["GRS_1980",6378137,298.257222101]],
>
>      datum=NAD83
>
>   END
...
> UNIT["Foot_US",0.3048006096012192]]

Jim,

I think your main problem is the linear units.  The default units in
PROJ.4 is meters, but you need US feet.  Try adding a "units=us-ft"
in your PROJECTION declaration for the layer.  There is also some slight
chance you will need to alter the x_0 parameter by multiplying (or
dividing) by the meter/feet conversion ratio, but I doubt that will
be needed.

Good luck,

--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent