<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi:<div><br></div><div>I'm just not having any luck determining if a point is in a polygon.</div><div><br></div><div>1. I add my SRID for <span class="Apple-style-span" style="font-family: Courier; font-size: 13px; ">NAD_1983_UTM_Zone_10N to the spatial ref.</span></div><div>2. Importing the City of Vancouver's boundary (MULTILINESTRING) with -s 1 (SRID=1) and load the SQL file.</div><div>3. Pick a point within the city from google map: 49.284562,-123.097172</div><div>4. Check if the point is in the boundary:</div><div><br></div><div><div>covciv=# select ST_Contains(c.the_geom, ST_GeomFromText('POINT(49.284562 -123.097172)', 1)) from city_boundary c;</div><div> st_contains </div><div>-------------</div><div> f</div><div>(1 row)</div><div><br></div><div>covciv=# select ST_Contains(c.the_geom, ST_GeomFromText('POINT(-123.097172 49.284562)', 1)) from city_boundary c;</div><div> st_contains </div><div>-------------</div><div> f</div><div>(1 row)</div></div><div><br></div><div>I've also tried this with MULTIPOLYGON data in case MULTILINESTRING was the problem.</div><div> </div><div>Any tips would be appreciated.</div><div><br></div><div><br></div><div>The shape file's proj:</div><div><br></div><div><div>PROJCS["NAD_1983_UTM_Zone_10N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",-123],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0],UNIT["Meter",1]]</div></div><div><br></div><div>My proj with SRID = 1;</div><div><br></div><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 13px/normal Courier; ">INSERT INTO spatial_ref_sys VALUES (1,NULL,NULL,'PROJCS[''NAD_1983_UTM_Zone_10N'', GEOGCS[''GCS_North_American_1983'', DATUM[''D_North_American_1983'', SPHEROID[''GRS_1980'',6378137,298.257222101]], PRIMEM[''Greenwich'',0], UNIT[''Degree'',0.0174532925199433]], PROJECTION[''Transverse_Mercator''], PARAMETER[''False_Easting'',500000.0], PARAMETER[''False_Northing'',0.0], PARAMETER[''Central_Meridian'',-123.0], PARAMETER[''Scale_Factor'',0.9996], PARAMETER[''Latitude_of_Origin'',0.0], UNIT[''Meter'',1.0]]');</div></div></body></html>