<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7654.12">
<TITLE>Modify PCS Units</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P><FONT SIZE=2 FACE="Arial">Good afternoon all:</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">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:</FONT></P>
<P><FONT SIZE=2 FACE="Arial">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:</FONT></P>
<P><FONT SIZE=2 FACE="Arial">INSERT INTO public.spatial_ref_sys (srid, auth_name, auth_srid, srtext, proj4text)</FONT>
<BR><FONT SIZE=2 FACE="Arial">VALUES</FONT>
<BR><FONT SIZE=2 FACE="Arial">(</FONT>
<BR><FONT SIZE=2 FACE="Arial">104307, </FONT>
<BR><FONT SIZE=2 FACE="Arial">'NONE', </FONT>
<BR><FONT SIZE=2 FACE="Arial">104307, </FONT>
<BR> <FONT SIZE=2 FACE="Arial">'PROJCS["NAD83 / Florida GDL Albers",</FONT>
<BR> <FONT SIZE=2 FACE="Arial">GEOGCS["GCS_North_American_1983",</FONT>
<BR> <FONT SIZE=2 FACE="Arial">DATUM["D_North_American_1983", SPHEROID["GRS_1980",6378137,298.257222101]],</FONT>
<BR> <FONT SIZE=2 FACE="Arial">PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],</FONT>
<BR> <FONT SIZE=2 FACE="Arial">PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["False_Easting",248.548],</FONT>
<BR> <FONT SIZE=2 FACE="Arial">PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-84],</FONT>
<BR> <FONT SIZE=2 FACE="Arial">PARAMETER["Standard_Parallel_1",24],PARAMETER["Standard_Parallel_2",31.5],</FONT>
<BR> <FONT SIZE=2 FACE="Arial">PARAMETER["Latitude_Of_Origin",24],UNIT["Mile_US",1609.347218694438]]',</FONT>
<BR><FONT SIZE=2 FACE="Arial">'+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'</FONT></P>
<P><FONT SIZE=2 FACE="Arial">);</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">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.</FONT></P>
<P><FONT SIZE=2 FACE="Arial">Cheers,</FONT>
<BR><FONT SIZE=2 FACE="Arial">Rob</FONT>
</P>
<P><FONT FACE="Eras Demi ITC">Robert E. Clift</FONT>
<BR><FONT SIZE=2 FACE="Eras Medium ITC">Florida Department of Health</FONT>
<BR><FONT SIZE=2 FACE="Eras Medium ITC">Division of Health Access & Tobacco</FONT>
<BR><FONT SIZE=2 FACE="Eras Medium ITC">4025 Esplanade Way, Room 330B</FONT>
<BR><FONT SIZE=2 FACE="Eras Medium ITC">Tallahassee, FL 32311</FONT>
<BR><FONT SIZE=2 FACE="Eras Medium ITC">Phone: 850.245.4444 ext 2471</FONT>
<BR><FONT SIZE=2 FACE="Eras Medium ITC">Fax: 850.414.6470</FONT>
</P>
</BODY>
</HTML>