[postgis-users] Modify PCS Units

[Robert_Clift at doh.state.fl.us]

Good afternoon all:

I'm working on a project in which I feed parameters to our intern for
input into R. She's using the spatstat package to do some point pattern
analysis. I'll enable PL/R and incorporate her code once she's got the
method for returning distance and density plots neatly nailed to the
floor. Here's my issue:

The plots she's returned from the sample data I've provided are measured
in units of meters or square meters. I thought it would be simple to
customize our standard projection (Florida Albers - EPSG 3086) to return
point pairs in miles instead of meters. I've been striking out and hope
someone can point out what I'm doing wrong. I've tried many variations
on:

INSERT INTO public.spatial_ref_sys (srid, auth_name, auth_srid, srtext,
proj4text)
VALUES
(
104307, 
'NONE', 
104307, 
	'PROJCS["NAD83 / Florida GDL Albers",
	GEOGCS["GCS_North_American_1983",
		DATUM["D_North_American_1983",
SPHEROID["GRS_1980",6378137,298.257222101]],
	PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],
	
PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["False_Easting",248.548]
,
	PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-84],
	
PARAMETER["Standard_Parallel_1",24],PARAMETER["Standard_Parallel_2",31.5
],
	
PARAMETER["Latitude_Of_Origin",24],UNIT["Mile_US",1609.347218694438]]',
'+proj=aea +lat_1=24 +lat_2=31.5 +lat_0=24 +lon_0=-84 +x_0=248.548
+y_0=0 +a=3963.19 +rf=298.2533112582781456953642384106 +units=mi'
);

I converted the False_Easting value to miles and added the '+a' and
'+rf' switches as a last-ditch fix to no avail. Thanks in advance for
your advice.

Cheers,
Rob

Robert E. Clift
Florida Department of Health
Division of Health Access & Tobacco
4025 Esplanade Way, Room 330B
Tallahassee, FL 32311
Phone: 850.245.4444 ext 2471
Fax: 850.414.6470

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20100803/cc51dc8f/attachment.html>