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
_______________________________________________
postgis-users mailing list
[email protected]
http://postgis.refractions.net/mailman/listinfo/postgis-users
