Subject: [ggl] intersection of two vectors of polygons
From: Barend Gehrels (barend)
Date: 2011-04-29 18:20:25
Let me just answer on my own mail:
> It is not only this. A multi_polygon must be seen as a single
> geometry. It is one geometry, with one attribute (or one set of
> attributes). All the individual polygons making up that geometry share
> the same attribute(s).
>
> If you're into (spatial) databases, see it as one row, with a
> geometry-column (the multi-polygon) and one or more other columns (the
> attributes).
While this is the case in general and in most Geometry libraries and
databases, in Boost.Geometry the multi-polygons *can* have an attribute
(e.g. ID) for each polygon it contains, and this is (more or less by
chance, or by the approach) the case in the code I created. This might
be useful sometimes, but it is not required by the OGC specs and not
used in databases.
In the code I created it is useless, the output multi-polygon contains
zero or more polygons of the input polygons, but all with ID zero
(because of the default constructor).
If you want to work it out further, you have to create another model of
a multi-polygon, having two ID's (owner_id and mineral_rights_id), and
assigning these ID's right before or after intersection (comparable to
the mentioned property merge).
Each derived multi-polygon must be registered as a multi-polygon (I
didn't comment on that, but a large piece of the sample code I created
(everything in namespace traits) is necessary for registration of the
polygon_with_id - it might be convenient to have a macro for this in the
future).
Regards, Barend
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/ggl/attachments/20110430/9501437f/attachment.html