1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
/**********************************************************************
* $Id: PointBuilder.cpp 1820 2006-09-06 16:54:23Z mloskot $
*
* GEOS - Geometry Engine Open Source
* http://geos.refractions.net
*
* Copyright (C) 2001-2002 Vivid Solutions Inc.
*
* This is free software; you can redistribute and/or modify it under
* the terms of the GNU Lesser General Public Licence as published
* by the Free Software Foundation.
* See the COPYING file for more information.
*
**********************************************************************/
#include <geos/operation/overlay/PointBuilder.h>
#include <geos/operation/overlay/OverlayOp.h>
#include <geos/geomgraph/Node.h>
#include <geos/geomgraph/EdgeEndStar.h>
#include <geos/geomgraph/Label.h>
#include <vector>
#include <cassert>
#ifndef GEOS_DEBUG
#define GEOS_DEBUG 0
#endif
using namespace std;
using namespace geos::geomgraph;
using namespace geos::geom;
namespace geos {
namespace operation { // geos.operation
namespace overlay { // geos.operation.overlay
/*
* @return a list of the Points in the result of the specified
* overlay operation
*/
vector<Point*>*
PointBuilder::build(OverlayOp::OpCode opCode)
{
extractNonCoveredResultNodes(opCode);
return resultPointList;
}
/*
* Determines nodes which are in the result, and creates Point for them.
*
* This method determines nodes which are candidates for the result via their
* labelling and their graph topology.
*
* @param opCode the overlay operation
*/
void
PointBuilder::extractNonCoveredResultNodes(OverlayOp::OpCode opCode)
{
map<Coordinate*,Node*,CoordinateLessThen> &nodeMap =
op->getGraph().getNodeMap()->nodeMap;
map<Coordinate*,Node*,CoordinateLessThen>::iterator it=nodeMap.begin();
for (; it!=nodeMap.end(); ++it)
{
Node *n=it->second;
// filter out nodes which are known to be in the result
if (n->isInResult()) continue;
// if an incident edge is in the result, then
// the node coordinate is included already
if (n->isIncidentEdgeInResult()) continue;
if ( n->getEdges()->getDegree() == 0 ||
opCode == OverlayOp::opINTERSECTION )
{
/**
* For nodes on edges, only INTERSECTION can result
* in edge nodes being included even
* if none of their incident edges are included
*/
Label *label=n->getLabel();
if (OverlayOp::isResultOfOp(label,opCode))
filterCoveredNodeToPoint(n);
}
}
}
void
PointBuilder::filterCoveredNodeToPoint(const Node *n)
{
const Coordinate& coord=n->getCoordinate();
if(!op->isCoveredByLA(coord)) {
Point *pt=geometryFactory->createPoint(coord);
resultPointList->push_back(pt);
}
}
} // namespace geos.operation.overlay
} // namespace geos.operation
} // namespace geos
/**********************************************************************
* $Log$
* Revision 1.21 2006/06/05 15:36:34 strk
* Given OverlayOp funx code enum a name and renamed values to have a lowercase prefix. Drop all of noding headers from installed header set.
*
* Revision 1.20 2006/03/17 13:24:59 strk
* opOverlay.h header splitted. Reduced header inclusions in operation/overlay implementation files. ElevationMatrixFilter code moved from own file to ElevationMatrix.cpp (ideally a class-private).
*
* Revision 1.19 2006/03/02 12:12:01 strk
* Renamed DEBUG macros to GEOS_DEBUG, all wrapped in #ifndef block to allow global override (bug#43)
*
* Revision 1.18 2006/02/19 19:46:49 strk
* Packages <-> namespaces mapping for most GEOS internal code (uncomplete, but working). Dir-level libs for index/ subdirs.
*
* Revision 1.17 2005/11/15 12:14:05 strk
* Reduced heap allocations, made use of references when appropriate,
* small optimizations here and there.
*
* Revision 1.16 2005/06/28 01:07:02 strk
* improved extraction of result points in overlay op
*
* Revision 1.15 2005/06/25 10:20:39 strk
* OverlayOp speedup (JTS port)
*
* Revision 1.14 2005/02/05 05:44:47 strk
* Changed geomgraph nodeMap to use Coordinate pointers as keys, reduces
* lots of other Coordinate copies.
*
**********************************************************************/
|
