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 135 136 137 138 139 140 141 142 143 144 145 146
|
/*
Copyright (c) 2009-2013,
Ken Arroyo Ohori g.a.k.arroyoohori@tudelft.nl
Hugo Ledoux h.ledoux@tudelft.nl
Martijn Meijers b.m.meijers@tudelft.nl
All rights reserved.
This file is part of pprepair: 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.
Licensees holding a valid commercial license may use this file in
accordance with the commercial license agreement provided with
the software.
This file 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.
*/
#ifndef IOWORKER_H
#define IOWORKER_H
#include "FaceInfo.h"
#include "definitions/CGALDefinitions.h"
class IOWorker {
public:
// Construction
IOWorker();
// Main operations
bool addToTriangulation(Triangulation &triangulation, TaggingVector &edgesToTag, const char *file, unsigned int schemaIndex);
bool tagTriangulation(Triangulation &triangulation, TaggingVector &edgesToTag);
bool makeAllHolesValid(Triangulation &triangulation);
bool splitRegions(Triangulation &triangulation, double ratio);
bool repairTrianglesByNumberOfNeighbours(Triangulation &triangulation, bool alsoUniverse);
bool repairTrianglesByAbsoluteMajority(Triangulation &triangulation, bool alsoUniverse);
bool repairTrianglesByLongestBoundary(Triangulation &triangulation, bool alsoUniverse);
bool repairRegionsByLongestBoundary(Triangulation &triangulation, bool alsoUniverse);
bool repairRegionsByRandomNeighbour(Triangulation &triangulation, bool alsoUniverse);
bool repairByPriorityList(Triangulation &triangulation, const char *file);
bool repairEdgeMatching(Triangulation &triangulation, const char *file);
bool matchSchemata(Triangulation &triangulation);
void removeConstraints(Triangulation &triangulation);
void removeVertices(Triangulation &triangulation);
bool reconstructPolygons(Triangulation &triangulation, std::vector<std::pair<PolygonHandle *, Polygon> > &outputPolygons);
bool exportPolygons(std::vector<std::pair<PolygonHandle *, Polygon> > &outputPolygons, const char *file, bool withProvenance);
bool exportTriangulation(Triangulation &t, const char *file, bool withNumberOfTags, bool withFields, bool withProvenance);
// Printing functions
void insertToStream(std::ostream &ostr, OGRFeatureDefn *layerDefinition, unsigned int indentation = 0, int schemaIndex = -1);
void insertToStream(std::ostream &ostr, const OGRFieldType &ft);
void insertToStream(std::ostream &ostr, const OGRwkbGeometryType >);
void insertTriangulationInfo(std::ostream &ostr, const Triangulation &t);
private:
// Data structures
struct FieldDescriptor {
char *file;
int layer;
int field;
FieldDescriptor(char *fl, int l, int fld) {
file = fl;
layer = l;
field = fld;
}
bool operator<(const FieldDescriptor &fe) const {
//std::cout << "COMP {" << std::endl << "\t" << file << "\t" << layer << "\t" << field << std::endl << "\t" << fe.file << "\t" << fe.layer << "\t" << fe.field << std::endl;
if (file < fe.file) return true;
if (file > fe.file) return false;
if (layer < fe.layer) return true;
if (layer > fe.layer) return false;
if (field < fe.field) return true;
return false;
}
};
struct FieldDefinition {
char *name;
OGRFieldType type;
OGRJustification justification;
int width;
int precision;
FieldDefinition(const char *n, OGRFieldType t, OGRJustification j, int w, int p) {
name = new char[strlen(n)+1];
strcpy(name, n);
type = t;
justification = j;
width = w;
precision = p;
}
~FieldDefinition() {
delete name;
}
bool matches(FieldDefinition *f) {
// For the moment, only exact matches.
// if we had reference data, we could match them together (change values to the ones of the reference and return true)
if (strcmp(name, f->name) != 0) return false;
if (type != f->type) return false;
if (justification != f->justification) return false;
if (width != f->width) return false;
if (precision != f->precision) return false;
return true;
}
};
struct FieldComparator {
bool operator() (Field * const &f1, Field * const &f2) const {
return (*f1) < (*f2);
}
};
// What is kept from input
std::vector<char *> fileNames;
std::vector<PolygonHandle *> polygons;
OGRFieldType schemaFieldType;
std::vector<FieldDefinition *> fields;
std::map<FieldDescriptor, unsigned int> fieldEquivalencies;
OGRSpatialReference* spatialReference;
// Internal special tags
PolygonHandle universe;
// Cached values
Triangulation::Face_handle startingSearchFace, startingSearchFaceInRing; // faces that are expected to be close to the next point to be added
// Helper functions
unsigned int removeDuplicateVertices(std::list<Point> &ring);
std::vector<Ring *> splitRing(Ring &ring);
void testRings(std::vector<Ring *> &outerRings, std::vector<Ring *> &innerRings, std::vector<std::vector<Ring> > &classification, long fid);
void copyFields(OGRFeature *ogrfeature, PolygonHandle *handle);
void tagStack(std::stack<Triangulation::Face_handle> &positiveStack, std::stack<Triangulation::Face_handle> &negativeStack, PolygonHandle *positiveHandle, PolygonHandle *negativeHandle);
void tagStack(std::stack<Triangulation::Face_handle> &stack, PolygonHandle *handle);
std::list<Triangulation::Vertex_handle> *getBoundary(Triangulation::Face_handle face, int edge, PolygonHandle *polygon);
void addtoCount(std::map<PolygonHandle *, unsigned int> &count, PolygonHandle *ph);
void addToLength(std::map<PolygonHandle *, double> &lengths, PolygonHandle *ph, double length);
};
#endif
|