v0.2.53..v0.2.54 changeset DualHighwaySplitter.h

diff --git a/hoot-core/src/main/cpp/hoot/core/algorithms/splitter/DualHighwaySplitter.h b/hoot-core/src/main/cpp/hoot/core/algorithms/splitter/DualHighwaySplitter.h
new file mode 100644
index 0000000..83da11d
--- /dev/null
+++ b/hoot-core/src/main/cpp/hoot/core/algorithms/splitter/DualHighwaySplitter.h
@@ -0,0 +1,133 @@
+/*
+ * This file is part of Hootenanny.
+ *
+ * Hootenanny is free software: you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation, either version 3 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ * --------------------------------------------------------------------
+ *
+ * The following copyright notices are generated automatically. If you
+ * have a new notice to add, please use the format:
+ * " * @copyright Copyright ..."
+ * This will properly maintain the copyright information. DigitalGlobe
+ * copyrights will be updated automatically.
+ *
+ * @copyright Copyright (C) 2015, 2017, 2018, 2019, 2020 DigitalGlobe (http://www.digitalglobe.com/)
+ */
+
+#ifndef DUALHIGHWAYSPLITTER_H
+#define DUALHIGHWAYSPLITTER_H
+
+// GEOS
+#include <geos/geom/Coordinate.h>
+
+// Hoot
+#include <hoot/core/util/Units.h>
+#include <hoot/core/ops/OsmMapOperation.h>
+#include <hoot/core/criterion/HighwayCriterion.h>
+
+#include <unordered_set>
+
+namespace hoot
+{
+
+class OsmMap;
+class Way;
+
+/**
+ * Splits all "divided=yes" highways into two oneway streets. All intersections should be split
+ * before this method is called. Failure to do so will result in undefined behavior. A number of
+ * assumptions must be made to do this including assumptions about the direction of travel on
+ * roads (right or left hand drivers).
+ */
+class DualHighwaySplitter : public OsmMapOperation
+{
+public:
+
+  static std::string className() { return "hoot::DualHighwaySplitter"; }
+
+  static int logWarnCount;
+
+  typedef enum DrivingSide
+  {
+    Left,
+    Right
+  } DrivingSide;
+
+  DualHighwaySplitter();
+  DualHighwaySplitter(const std::shared_ptr<const OsmMap>& map, DrivingSide drivingSide,
+                      Meters splitSize);
+
+  std::shared_ptr<OsmMap> splitAll();
+
+  static std::shared_ptr<OsmMap> splitAll(const std::shared_ptr<const OsmMap>& map,
+                                          DrivingSide drivingSide, Meters defaultSplitSize);
+
+  virtual void apply(std::shared_ptr<OsmMap>& map);
+
+  virtual QString getInitStatusMessage() const
+  { return "Splitting divided highways into two one way streets..."; }
+
+  virtual QString getCompletedStatusMessage() const
+  { return "Split " + QString::number(_numAffected) + " divided highways"; }
+
+  virtual QString getDescription() const
+  { return "Splits all 'divided=yes' highways into two one way streets"; }
+
+  /**
+   * @see FilteredByCriteria
+   *
+   * This isn't actually using HighwayCriterion in the filtering, but for the purposes of reducing
+   * unnecessary conflate ops we don't need to run it unless we're running road conflation.
+   */
+  virtual QStringList getCriteria() const
+  { return QStringList(QString::fromStdString(HighwayCriterion::className())); }
+
+  virtual std::string getClassName() const { return className(); }
+
+private:
+
+  Meters _defaultSplitSize;
+  DrivingSide _drivingSide;
+  std::shared_ptr<const OsmMap> _map;
+  std::shared_ptr<OsmMap> _result;
+  std::unordered_set<long> _nodes;
+
+  // temporary variables for convenience
+  std::shared_ptr<Way> _left, _right;
+  std::shared_ptr<const Way> _working;
+  Meters _splitSize;
+
+  void _addConnector(long nodeId);
+  std::shared_ptr<Way> _createOneWay(const std::shared_ptr<const Way>& w, Meters bufferSize,
+                                     bool left);
+  void _createStub(const std::shared_ptr<Way>& dividedWay, long centerNodeId, long edgeNodeId);
+  double _dotProduct(const geos::geom::Coordinate& c1, const geos::geom::Coordinate& c2) const;
+  void _fixLanes(const std::shared_ptr<Way>& w);
+
+  /**
+   * Returns the node id of the nearest node to nid on w
+   */
+  long _nearestNode(long nid, const std::shared_ptr<const Way>& w);
+
+  geos::geom::Coordinate _normalizedVector(long nid1, long nid2);
+  bool _onRight(long intersectionId, const std::shared_ptr<Way>& inbound, long leftNn,
+                long rightNn);
+  void _reconnectEnd(long centerNodeId, const std::shared_ptr<Way>& edge);
+  void _splitIntersectingWays(long nid);
+  void _splitWay(long wid);
+};
+
+}
+
+#endif // DUALHIGHWAYSPLITTER_H