<div dir="ltr"><div>Brent<br><br></div>Great, informative answer.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 22, 2015 at 3:14 PM, Brent Wood <span dir="ltr"><<a href="mailto:pcreso@pcreso.com" target="_blank">pcreso@pcreso.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:verdana,helvetica,sans-serif;font-size:13px"><div dir="ltr"><span>Hi Aaron,</span></div><div dir="ltr"><br></div><div dir="ltr">Hopefully this (simplistic) description helps.<br></div><div dir="ltr"><br><span></span></div><div dir="ltr"><span>Setting a SRID value for a feature tells Postgis what coordinate reference system (CRS) the coordinates are in. This includes the unit, which can be any linear unit of measurement, such as degrees (measured at the surface of the earth), meters, feet, etc. Obviously getting this right is important to calculate distances correctly.<br></span></div><div dir="ltr"><span><br></span></div><div dir="ltr"><span>Different CRS's generally apply to different parts of the globe, and include the projection parameters. Note that every projection applies some distortion - trying to map (project) a part of a spherical surface to a flat plane has this problem. There are three main types of distortion - angular (changes shapes), area and distance. Normally to measure distance using a projected CRS, you'd use an equidistant projection (which minimises distance distortions) centered near the location you are trying to measure.</span></div><div dir="ltr"><span><br></span></div><div dir="ltr"><span>An alternative approach is to measure it on a spheroid, in 3D space, instead of on a projected 2D space. This is basically what a Postgis geography allows you to do. But the coordinate units in a geography are degree coordinates, so you need to specify that your coordinates are unprojected lon/lat values when you use them in a geography. The SRID (Spatial Reference ID) for such a CRS is 4326. <br></span></div><div dir="ltr"><span><br></span></div><div dir="ltr"><span>In your case, try:</span></div><div dir="ltr">SELECT * FROM users<br>WHERE ST_DWithin (users.location::geography, st_setsrid(st_makepoint (146.0,<br>138.19), 4326)::geography, 100);</div><div><br></div><div dir="ltr">This assumes your location is a geometry with a SRID of 4326 (ie: the coordinate values are unprojected lon/lat degrees). It then converts this geometry to a geography datatype, which it tests against a point geometry in SRID 4326, which is also converted to a geography for the test, to see if it is within 100m. So this SQL tests geography against geography datatype.</div><div dir="ltr"><br></div><div dir="ltr">If your location feature is not SRID 4326, you'll need to reproject (transform) it to 4326 to for this to work:</div><div dir="ltr"><br></div><div dir="ltr">SELECT * FROM users<br>WHERE ST_DWithin (ST_Transform(users.location,4326)::geography, st_setsrid(st_makepoint (146.0,<br>138.19), 4326)::geography, 100);</div><div dir="ltr"><br></div><div dir="ltr">(I haven't tested it, but I think this should work)</div><div dir="ltr"><br></div><div dir="ltr">Cheers</div><div dir="ltr"><br></div><div dir="ltr">Bent Wood<br></div><div><br>  </div><div style="font-family:verdana,helvetica,sans-serif;font-size:13px"> <div style="font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"> <div dir="ltr"> <hr size="1">  <font face="Arial"> <b><span style="font-weight:bold">From:</span></b> Aaron Lewis <<a href="mailto:the.warl0ck.1989@gmail.com" target="_blank">the.warl0ck.1989@gmail.com</a>><br> <b><span style="font-weight:bold">To:</span></b> <a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a> <br> <b><span style="font-weight:bold">Sent:</span></b> Monday, March 23, 2015 12:09 AM<br> <b><span style="font-weight:bold">Subject:</span></b> [postgis-users] Need help on basic concepts,        I just need really simple calculations<br> </font> </div> <div><br>Hi,<br><br>I've been searching online for days. Trying to understand why SRID is<br>required. So I picked some random value.<br><br>Now I'm need to retrieve POINTs within a circle range, e.g a circle at<br>(146.0, 138.19) with radius of 100 meters:<br><br>SELECT * FROM users<br>WHERE ST_DWithin (users.location, st_setsrid(st_makepoint (146.0,<br>138.19), 2600), 100);<br><br>It's very simple, but the result seems wrong. I have a record contains<br>a POINT(55 43) that matches this query.<br><br>Anyone know what's wrong with it?<span class="HOEnZb"><font color="#888888"><br><br>-- <br>Best Regards,<br>Aaron Lewis - PGP: 0x13714D33 - <a href="http://pgp.mit.edu/" target="_blank">http://pgp.mit.edu/</a><br>Finger Print:   9F67 391B B770 8FF6 99DC  D92D 87F6 2602 1371 4D33<br>_______________________________________________<br>postgis-users mailing list<br><a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br><a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users</a><br><br><br></font></span></div> </div> </div>  </div></div><br>_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org">postgis-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users</a><br></blockquote></div><br></div>