<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Thanks Regina that is it, I didn't have the right projection to begin with.<BR>
I needed NAD_1983_StatePlane_Nevada_East_FIPS_2701_Feet<BR>
and I was using NAD83 / Nevada East <BR>
<BR>
The projection I needed wasn't in the default spatial_ref_sys? I assumed all of the srid's would be in there, I guess not? I was able to find it at <A href="http://spatialreference.org">http://spatialreference.org</A> and I inserted into spatial_ref_sys. Now it is working as expected.<BR>
<BR>
Thanks again.<BR><BR>> Subject: RE: [postgis-users] Need Help with Transform<BR>> Date: Thu, 11 Sep 2008 15:07:47 -0400<BR>> From: robe.dnd@cityofboston.gov<BR>> To: postgis-users@postgis.refractions.net<BR>> <BR>> RH,<BR>> <BR>> I'm not aware of an easy way to tell if its in range (I think you can<BR>> from the proj4text lon a b blah blah blah args) defined in<BR>> spatial_ref_sys, but I'm not that technically savvy.<BR>> <BR>> If you get something like all the same long lat when you project back<BR>> to 4326 (WGS 84), or you get weird numbers like 135, 90 etc. then its<BR>> probably wrong.<BR>> <BR>> Given that all your below seem to project to below - I would say you are<BR>> wrong<BR>> POINT(-115.583333333333 90)<BR>> <BR>> First of all when I look in spatial ref<BR>> <BR>> SELECT srtext from spatial_ref_sys where srid = 32107 <BR>> <BR>> I get the below<BR>> "PROJCS["NAD83 / Nevada<BR>> East",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS<BR>> 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","<BR>> 6269"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0<BR>> 1745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]],PRO<BR>> JECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",34.75],PAR<BR>> AMETER["central_meridian",-115.5833333333333],PARAMETER["scale_factor",0<BR>> .9999],PARAMETER["false_easting",200000],PARAMETER["false_northing",8000<BR>> 000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32107"]]<BR>> "<BR>> <BR>> so if you know your data is not in Nevada - you guessed wrong about the<BR>> spatial reference system<BR>> <BR>> <BR>> Well even if its in the range doesn't necessarily mean you have the<BR>> right projection. <BR>> <BR>> The trick I normally use when I have no clue what projection a data set<BR>> is in but I have some guesses e.g. I know its a meter projection or the<BR>> state etc) and I know say the long lat of a point that should be in that<BR>> range would be something like below (and that usually 0s down to 1 or 2<BR>> likely candidates). If the projection doesn't exist in spatial_ref_sys<BR>> (or I have the ESRI .prj file) - then I go to<BR>> http://spatialreference.org and do some searches and load up some likely<BR>> culprits in my spatial_ref_sys table and retry the below db or upload<BR>> the prj to spatial ref.<BR>> <BR>> SELECT srid, srtext, proj4text<BR>> FROM spatial_ref_sys <BR>> WHERE proj4text ILIKE '%units=m%' AND srtext ILIKE '%NAD%' <BR>> AND ST_Transform(ST_SetSRID(ST_MakePoint(-71.015,42.372), 4326), srid)<BR>> && <BR>> ST_Expand(ST_SetSRID(ST_MakePoint(239898, 902552), srid),1000)<BR>> <BR>> <BR>> Hope that helps,<BR>> Regina<BR>> ________________________________<BR>> <BR>> From: postgis-users-bounces@postgis.refractions.net<BR>> [mailto:postgis-users-bounces@postgis.refractions.net] On Behalf Of R H<BR>> Sent: Thursday, September 11, 2008 11:42 AM<BR>> To: PostGIS Users Discussion<BR>> Subject: RE: [postgis-users] Need Help with Transform<BR>> <BR>> <BR>> Thanks for the Help Regina, below are some of the points I attempted.<BR>> <BR>> setsrid(MakePoint(834934.2263,26712690.6089, 32107) <BR>> setsrid(MakePoint(834533.210694,26712881.1962, 32107) <BR>> setsrid(MakePoint(834224.086891,26712318.3331, 32107) <BR>> <BR>> Is there a way to check if it is in the range of the projection? <BR>> <BR>> <BR>> <BR>> ________________________________<BR>> <BR>> <BR>> From: lr@pcorp.us<BR>> To: postgis-users@postgis.refractions.net<BR>> Subject: RE: [postgis-users] Need Help with Transform<BR>> Date: Fri, 5 Sep 2008 18:19:58 -0400<BR>> <BR>> <BR>> <BR>> RH,<BR>> <BR>> This kind of thing happens if your input points are not actually in the<BR>> projection you say they are with setsrid.<BR>> <BR>> Give us an example 32107 point you are using. I suspect it doesn't fall<BR>> in the range of a valid 32107 point.<BR>> <BR>> Hope that helps,<BR>> Regina<BR>> <BR>> ________________________________<BR>> <BR>> From: postgis-users-bounces@postgis.refractions.net<BR>> [mailto:postgis-users-bounces@postgis.refractions.net] On Behalf Of R H<BR>> Sent: Friday, September 05, 2008 1:13 PM<BR>> To: postgis-users@postgis.refractions.net<BR>> Subject: [postgis-users] Need Help with Transform<BR>> <BR>> <BR>> <BR>> <BR>> ________________________________<BR>> <BR>> <BR>> Transform is returning to the same value, for different inputs? Below<BR>> is what I have done...<BR>> <BR>> select AddGeometryColumn('data_stage', 'table', 'geo4269', 4269,<BR>> 'POINT', 2)<BR>> select AddGeometryColumn('data_stage', 'table', 'geo32107', 32107,<BR>> 'POINT', 2)<BR>> <BR>> update table<BR>> set geo32107= setsrid(MakePoint(x, y), 32107) <BR>> where ....<BR>> <BR>> update table set geo4269 = transfrom(geo32107, 4269)<BR>> where ....<BR>> <BR>> column geo4269 all have the same vaule. select AsText(geo4269) returns<BR>> 'POINT(-115.583333333333 90)' for all rows. And select AsText(geo32107)<BR>> returns the point as input.<BR>> <BR>> Am I using transform wrong? Or any other problems? Any help many<BR>> thanks...<BR>> <BR>> <BR>> <BR>> <BR>> -----------------------------------------<BR>> The substance of this message, including any attachments, may be<BR>> confidential, legally privileged and/or exempt from disclosure<BR>> pursuant to Massachusetts law. It is intended<BR>> solely for the addressee. If you received this in error, please<BR>> contact the sender and delete the material from any computer.<BR>> <BR>> _______________________________________________<BR>> postgis-users mailing list<BR>> postgis-users@postgis.refractions.net<BR>> http://postgis.refractions.net/mailman/listinfo/postgis-users<BR><BR></body>
</html>