Top 5 Recent Articles
ARTICLES CATEGORIES
- Algorithms (20)
- All (400)
- Biography (1)
- Blog (45)
- Business Requirements (1)
- Commentary (1)
- Customers (2)
- Data Models (1)
- Education (2)
- GeoRaptor (13)
- GPS (1)
- Image Processing (2)
- Import Export (8)
- Licensing (2)
- LiDAR (1)
- Linear Referencing (4)
- Manifold GIS (3)
- Mapping (1)
- MySQL Spatial (7)
- Networking and Routing (including Optimization) (4)
- Open Source (18)
- Oracle Spatial and Locator (193)
- Partitioning (1)
- PostGIS (34)
- Projections (1)
- Published Articles (1)
- qGIS (1)
- Recommendations (1)
- Services (1)
- Software Change Log (1)
- Source Code (35)
- Space Curves (9)
- Spatial Database Functions (108)
- Spatial DB comparison (1)
- Spatial XML Processing (11)
- SQL Server Spatial (92)
- Standards (3)
- Stored Procedure (15)
- Tessellation or Gridding (10)
- Tools (2)
- Topological Relationships (1)
- Training (2)
- Triangulation (2)
Free Union, Intersection, Xor and Difference Functions for Oracle Locator – Part 2
In this blog we will execute the four overlay types on two simple rectangular polygons.
In addition, I have not tested “real world” polygon data using multiple projections. I am just showing what is possible.
Here is the SQL that I use to create two test polygon geometries.
SELECT 'GEOM1' AS gtype, sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(1,1,10,10)) AS geom FROM dual UNION ALL SELECT 'GEOM2' AS gtype, sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(5,5,15,15)) AS geom FROM dual
Here are what they look like.
Here is the SQL that I used to combine the above test polygon geometries.
SELECT 'UNION' AS gtype, SC4O.ST_Union(sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(1,1,10,10)), sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(5,5,15,15)), 1) AS rGeom FROM dual UNION ALL SELECT 'INTERSECTION' AS gtype, SC4O.ST_Intersection(sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(1,1,10,10)), sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(5,5,15,15)), 1) AS rGeom FROM dual UNION ALL SELECT 'XOR' AS gtype, SC4O.ST_XOr(sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(1,1,10,10)), sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(5,5,15,15)), 1) AS rGeom FROM dual UNION ALL SELECT 'DIFFERENCE' AS gtype, SC4O.ST_Difference(sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(1,1,10,10)), sdo_geometry(2003,82469,NULL,sdo_elem_info_array(1,1003,3),sdo_ordinate_array(5,5,15,15)), 1) AS rGeom FROM dual;
The following image shows the result of executing four Overlay types on two simple rectangular polygons visually using GeoRaptor’s neat Thumbnail capability.
NOTE 1: JTS doesn’t currently support geodetic computations. If you wish to use this package with geodetic data then you should:
- Transform (sdo_cs.transform) your data to a suitable projected coordinate system, eg UTM;
- Call the relevant JTS function;
- Transform the result back to you original geodetic SRID.
Documentation
- MySQL Spatial General Function Documentation
- Oracle LRS Object Documentation
- Oracle Spatial Exporter Package Documentation
- Oracle Spatial Object Function Documentation
- Oracle Spatial Object Function Documentation (Multi Page Version)
- PostGIS pl/pgSQL Function Documentation
- SC4O Oracle Java Topology Suite (Stored Procedures) Package Documentation
- SQL Server Spatial General TSQL Function Documentation
- SQL Server Spatial LRS TSQL Function Documentation