Add a +to_meter parameter to your proj4 definition, that might nail it. I'm not sure specifying +units actually does any math.
P On Tue, Aug 3, 2010 at 10:48 AM, <[email protected]> wrote: > 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 > > _______________________________________________ postgis-users mailing list [email protected] http://postgis.refractions.net/mailman/listinfo/postgis-users
