Question

this is probably simple to those who know (I hope!)

I have an Oracle spatial database with a geometry column containing a list of node points, in northing/easting format (if it’s relevent!)

I need to select those objects which fall within a given radius of a given point.

Northings and Eastings are 1 meter apart which makes it a bit easier.

Ideally this should include objects which cross the area even if their node points fall outside it.

Is this an easy-ish query? Maybe using SDO_WITHIN_DISTANCE?

The table looks like this:

MyTable
ID NUMBER
NAME VARCHAR2(20)
DESC VARCHAR2(50)
GEOM SDO_GEOMETRY

Thanks for any help!

Answer 1

You can do this one of two ways. First, as you mentioned, SDO_WITHIN_DISTANCE is a valid approach.

select 
    *
from center_point a
inner join target_points b
    on a.id = 1
    and sdo_within_distance( b.shape, a.shape, 'distance = 10' ) = 'TRUE'
;

In this case, the distance is in linear units defined by a’s spatial reference. Oracle treats the coordinates as Cartesian so you will need to make sure you have a linear coordinate system before using this operator (as opposed to angular lat/lon units). Since you are working with northings/eastings, I think you’ll be okay as long as the points you are comparing against are in the same spatial reference.

This approach uses an inner-loop to solve the query so not very efficient if you have a lot of points to compare against. Also, Oracle Spatial is VERY picky about the order of operands in the SDO functions so you might need to play around with parameter order to find the sweetspot. If your query runs for a long period, try switching the first and second parameter of your sdo operator. You can also play with the order of the ‘from’ and ‘inner join’ tables using the /*+ ORDERED */ hind after SELECT.

Another approach is to buffer the geometry and compare against the buffer.

select 
    *
from center_point a
inner join target_points b
    on a.id = 1
    and sdo_relate( b.shape, sdo_buffer(a.shape, 0.05 ), 'mask=anyinteract' ) = 'TRUE'
;

Keep in mind that whatever is in the second parameter of the SDO_RELATE (called the window) will not have a spatial index if you transform it like we are here with the buffer.

If you plan on doing this with several points, it is recommended to build a table where all of the source points are buffered. Then create a spatial index against the buffered areas and compare that to your target points.

For example:

create table point_bufs unrecoverable as
select sdo_buffer (a.shape, b.diminfo, 1.35)
from centerpoint a, user_sdo_geom_metadata b
where table_name='CENTERPOINT'
  and column_name='SHAPE';

select
    a.gif,
    b.gid 
from target_points a, 
     point_bufs b
where sdo_relate(a.shape, b.shape, 'mask=anyinteract querytype=join') = 'TRUE'
;

NOTE: When intersecting points with polygons, you always want to polygon to be in the window position of the sdo_relate (which is the second parameter). This will ensure your spatial index is used correctly.