src/noding/FastNodingValidator.cpp

   1 /**********************************************************************
   2  *
   3  * GEOS - Geometry Engine Open Source
   4  * http://geos.osgeo.org
   5  *
   6  * Copyright (C) 2006 Refractions Research Inc.
   7  *
   8  * This is free software; you can redistribute and/or modify it under
   9  * the terms of the GNU Lesser General Licence as published
  10  * by the Free Software Foundation.
  11  * See the COPYING file for more information.
  12  *
  13  **********************************************************************
  14  *
  15  * Last port: noding/FastNodingValidator.java rev. ??? (JTS-1.8)
  16  *
  17  **********************************************************************/
  18
  19 #include <geos/noding/FastNodingValidator.h>
  20 #include <geos/noding/MCIndexNoder.h> // for checkInteriorIntersections()
  21 #include <geos/noding/SingleInteriorIntersectionFinder.h>
  22 #include <geos/util/TopologyException.h> // for checkValid()
  23 #include <geos/geom/Coordinate.h>
  24 #include <geos/io/WKTWriter.h> // for getErrorMessage()
  25
  26 #include <string>
  27 #include <iostream>
  28
  29 namespace geos {
  30 namespace noding { // geos.noding
  31
  32 /*private*/
  33 void
  34 FastNodingValidator::checkInteriorIntersections()
  35 {
  36         isValidVar = true;
  37         segInt.reset(new SingleInteriorIntersectionFinder(li));
  38         MCIndexNoder noder;
  39         noder.setSegmentIntersector(segInt.get());
  40         noder.computeNodes(&segStrings);
  41         if (segInt->hasIntersection()) {
  42                 isValidVar = false;
  43                 return;
  44         }
  45 }
  46
  47 /*public*/
  48 std::string
  49 FastNodingValidator::getErrorMessage() const
  50 {
  51         using geos::io::WKTWriter;
  52         using geos::geom::Coordinate;
  53
  54         if (isValidVar) return std::string("no intersections found");
  55
  56         //return std::string("found non-noded intersection etc etc");
  57
  58         const std::vector<Coordinate>& intSegs = segInt->getIntersectionSegments();
  59         assert(intSegs.size() == 4);
  60         return "found non-noded intersection between "
  61                 + WKTWriter::toLineString(intSegs[0], intSegs[1])
  62                 + " and "
  63                 + WKTWriter::toLineString(intSegs[2], intSegs[3]);
  64 }
  65
  66 void
  67 FastNodingValidator::checkValid()
  68 {
  69         execute();
  70         if (! isValidVar)
  71         {
  72                 //std::cerr << "Not valid: " << getErrorMessage() << " interior intersection: " << segInt->getInteriorIntersection() << std::endl;
  73                 throw util::TopologyException(getErrorMessage(), segInt->getInteriorIntersection());
  74         }
  75 }
  76
  77 } // namespace geos.noding
  78 } // namespace geos
  79
  80 /**********************************************************************
  81  * $Log$
  82  **********************************************************************/