In Article previous we have discussed some SQL functions for spatial processing . In this article we will see a function a little more complicated to put at work: the aggregation function for two layers of polygons. This operation is the one called ** Union** in ArcGis. The purpose of the operation is to obtain a layer with all the intersections of the entities of the two original layers, and, also, all non-intersecting entities.

Let’s use the following example. We have a first layer “rectangles”

A second layer “circles”

That overlaps spatially:

What we want to obtain is a layer containing all the polygons of both layers but creating polygons to overlay areas:

We will proceed by following the same method used in the previous article, i.e. by developing an SQL query including several sub queries.

In the first place, we will develop a query that will build the polygons of the overlapping areas (the intersections) of both layers. These polygons appear in yellow in the next picture.

We will resume the request used in the previous article:

*select**circle.id as id1, rectangles.id as*** id2**,

*st_intersection*

*(***circle.geom**

*, rectangles.geom*

*) as*

*the_geom*

*from*

*circle, rectangles*

*where*

*st_intersects*

*(*

*circle.geom, rectangles.geom*

*)*

*group by*

*circle.id, rectangles.id*You will notice that the select is built not just to retrieve intersected geometries, but also the identifiers of the original entities (*circle.id***and** ** rectangles.id)**. Of course, you can recover other attributes, or all by changing this part of the request.

To these polygons, we will add the rest of the areas occupied by the rectangle layer (in yellow in the following picture)

**select** *0 as id1, rectangles.id as***id2**, **st_difference** **(** ** rectangles.geom**,

**(select** **multi** **ST_** **(ST_Union** **(the_geom)) as** **the_geom** **from (****select** **st_intersection** **(**circle.geom*, rectangles.geom***) as** **the_geom** **from** *circle, rectangles***where** **st_intersects** **(** *circle.geom, rectangles.geom***))** **as** **the_geom**)**)** **as** **the_geom** **from rectangles**

What we get to do with this request is to extract all the intersection areas (as in the previous query) and convert them in a single multi-polygon. The request removes the difference of this multi- polygon with the layer of rectangles.

Finally, these polygons, we are going to add the rest of the areas occupied by the layer circles to these polygons (yellow in the following image)

We will proceed likewise the rectangles, by using a multi- polygon of intersection areas and calculating the difference with the layer circles:

**select** *circle.id as id1, 0 as***id2** *,***st_difference (** *circle.geom***, (** **select** **st_multi** **(** **st_union** **(** **the_geom** **)) as** **the_geom** **from (** **select** **st_intersection** **(** *circle.geom, rectangles.geom***) as** **the_geom** **from** *circle, rectangles***where** **st_intersects** **(** *circle.geom, rectangles.geom***))** **as** **the_geom** **)) as** **the_geom** **from** **circle**

The union of these three requests produces the outcome layer we expected :

**select *** **from** **(select circle.id as id1, rectangles.id** **as** **id2, st_intersection (** **circle.geom, rectangles.geom** **) as** **the_geom** **from** **circle, rectangles** **where** **st_intersects** **(** **circle.geom, rectangles.geom** **)** **group by** **circle.id, rectangles.id** **union all** **select 0 as id1, rectangles.id as id2, st_difference (** **rectangles.geom** **, (** **select** **st_multi** **(** **st_union** **(** **the_geom** **)) as** **the_geom** **from (** **select** **st_intersection** **(** **circle.geom, rectangles.geom** **) as** **the_geom** **from** **circle, rectangles** **where** **st_intersects** **(** **circle.geom, rectangles.geom** **))** **as** **the_geom** **)) as** **the_geom** **from rectangles** **union all** **select circle.id as id1, 0 as id2, st_difference (** **circle.geom** **, (** **select** **st_multi** **(** **st_union** **(** **the_geom** **)) as** **the_geom** **from (** **select** **st_intersection** **(** **circle.geom, rectangles.geom** **) as** **the_geom** **from** **circle, rectangles** **where** **st_intersects** **(** **circle.geom, rectangles.geom** **))** **as** **the_geom** **)) as** **the_geom** **from** **circle** **) as T1**

Thank you for this