Thank you so much.  The areas do not match the qGIS derived area exactly but that appears to be a well known issue.<br>They're certainly close enough for my needs.<br><br>The bonus is that I can now tell my boss that open source software is not the problem!<br>
<br>I can't thank you enough! <br>Charles <br><br><div class="gmail_quote">On Mon, Jun 1, 2009 at 5:22 PM, Paragon Corporation <span dir="ltr"><<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Charles,<br>
<br>
Which SRID did you use to import your data?  Sounds like that is the<br>
problem.  If you are using QGIS to view the data, I am guessing you used the<br>
shape file?  And it was reading the prj file to determine the native<br>
projection, but you overwrote in SPIT 3310 (which is wrong).<br>
<br>
For census data -- you should set SRID  to 4269 in SPIT (NAD 83 long lat)<br>
and then your area transform below should work fine.<br>
<br>
So bring the data in as SRID=4269<br>
Then transform 3310 as you are doing  (national atlas equal area - 2163<br>
might also be suitable and that covers all of US, but probably not as<br>
accurate as California Albers)<br>
<br>
Leo<br>
<div><div></div><div class="h5"><br>
-----Original Message-----<br>
From: <a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a><br>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net">postgis-users-bounces@postgis.refractions.net</a>] On Behalf Of Charles<br>
Blankenship<br>
Sent: Monday, June 01, 2009 8:04 PM<br>
To: <a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
Subject: [postgis-users] trouble getting correct units from "area(transform"<br>
<br>
Disclaimer: I'm very new to GIS.<br>
<br>
I want to calculate the areas of school districts in Nevada (and California<br>
at some point) using census shape files.<br>
<a href="http://www.census.gov/geo/cob/bdy/sn/sn00shp/sn32_d00_shp.zip" target="_blank">http://www.census.gov/geo/cob/bdy/sn/sn00shp/sn32_d00_shp.zip</a><br>
<br>
The data is unprojected so I know I must project it onto something good for<br>
equal areas.  I'm trying NAD83 California Albers (SRID = 3310).<br>
<br>
It seems to me from reading the documentation and numerous other<br>
posts/tutorials, that the following code should return the area of the<br>
districts in meters (or something I can convert to meters).<br>
<br>
select area(transform(the_geom, 3310)) from mytable<br>
<br>
This returns a column of areas which match my unprojected areas exactly!<br>
<br>
I've fiddled with setting the constraint to SRID=3310, importing through the<br>
qGIS SPIT plugin, selecting into a new table, delete and resetting<br>
constraints, and on and on...<br>
<br>
<br>
FWIW, I can load the file into qGIS, set the projection to 3310, and the<br>
identify area feature calculates the "derived area" exactly as it should be.<br>
I can't for the life of me figure out how to even get into the approximately<br>
correct units.  (I've tried reprojecting to quite different SRIDs and can<br>
see the area change slightly, but it still appears to be in native map units<br>
or degrees).<br>
<br>
Thanks for any assistance,<br>
Charles<br>
<br>
</div></div>_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
</blockquote></div><br>