/****************************************************************************** * $Id: ili2reader.cpp 10646 2007-01-18 02:38:10Z warmerdam $ * * Project:  Interlis 2 Reader * Purpose:  Implementation of ILI2Reader class. * Author:   Markus Schnider, Sourcepole AG * ****************************************************************************** * Copyright (c) 2004, Pirmin Kalberer, Sourcepole AG * * Permission is hereby granted, free of charge, to any person obtaining a * copy of this software and associated documentation files (the "Software"), * to deal in the Software without restriction, including without limitation * the rights to use, copy, modify, merge, publish, distribute, sublicense, * and/or sell copies of the Software, and to permit persons to whom the * Software is furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included * in all copies or substantial portions of the Software. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER * DEALINGS IN THE SOFTWARE. ****************************************************************************/#include "ogr_ili2.h"#include "cpl_conv.h"#include "cpl_string.h"#include "ilihelper.h"#include "iomhelper.h"#include "ili2reader.h"#include "ili2readerp.h"using namespace std;CPL_CVSID("$Id: ili2reader.cpp 10646 2007-01-18 02:38:10Z warmerdam $");//// constants//static const char *ILI2_TID = "TID";static const char *ILI2_REF = "REF";static const int ILI2_STRING_TYPE = 0;static const int ILI2_COORD_TYPE = 1;static const int ILI2_ARC_TYPE = 2;static const int ILI2_POLYLINE_TYPE = 4;static const int ILI2_BOUNDARY_TYPE = 8;static const int ILI2_AREA_TYPE = 16; // also SURFACEstatic const int ILI2_GEOMCOLL_TYPE = 32;static const char *ILI2_COORD = "COORD";static const char *ILI2_ARC = "ARC";static const char *ILI2_POLYLINE = "POLYLINE";static const char *ILI2_BOUNDARY = "BOUNDARY";static const char *ILI2_AREA = "AREA";static const char *ILI2_SURFACE = "SURFACE";// // helper functions// int cmpStr(string s1, string s2) {    string::const_iterator p1 = s1.begin();  string::const_iterator p2 = s2.begin();  while (p1 != s1.end() && p2 != s2.end()) {    if (toupper(*p1) != toupper(*p2))      return (toupper(*p1) < toupper(*p2)) ? -1 : 1;    ++p1;    ++p2;  }  return (s2.size() == s1.size()) ? 0 :         (s1.size() < s2.size()) ? -1 : 1;}string ltrim(string tmpstr) {  unsigned int i = 0;  while (i < tmpstr.length() && (tmpstr[i] == ' ' || tmpstr[i] == '\t' || tmpstr[i] == '\r' || tmpstr[i] == '\n')) ++i;  return i > 0 ? tmpstr.substr(i, tmpstr.length()-i) : tmpstr;}string rtrim(string tmpstr) {  unsigned int i = tmpstr.length() - 1;  while (i >= 0 && (tmpstr[i] == ' ' || tmpstr[i] == '\t' || tmpstr[i] == '\r' || tmpstr[i] == '\n')) --i;  return i < tmpstr.length() - 1 ? tmpstr.substr(0, i+1) : tmpstr;}string trim(string tmpstr) {  tmpstr = ltrim(tmpstr);  tmpstr = rtrim(tmpstr);  return tmpstr;}int getGeometryTypeOfElem(DOMElement* elem) {  int type = ILI2_STRING_TYPE;    if (elem && elem->getNodeType() == DOMNode::ELEMENT_NODE) {    if (cmpStr(ILI2_COORD, XMLString::transcode(elem->getTagName())) == 0) {      type = ILI2_COORD_TYPE;    } else if (cmpStr(ILI2_ARC, XMLString::transcode(elem->getTagName())) == 0) {      type = ILI2_ARC_TYPE;    } else if (cmpStr(ILI2_POLYLINE, XMLString::transcode(elem->getTagName())) == 0) {      type = ILI2_POLYLINE_TYPE;    } else if (cmpStr(ILI2_BOUNDARY, XMLString::transcode(elem->getTagName())) == 0) {      type = ILI2_BOUNDARY_TYPE;    } else if (cmpStr(ILI2_AREA, XMLString::transcode(elem->getTagName())) == 0) {      type = ILI2_AREA_TYPE;    } else if (cmpStr(ILI2_SURFACE, XMLString::transcode(elem->getTagName())) == 0) {      type = ILI2_AREA_TYPE;    }  }  return type;}char *getObjValue(DOMElement *elem) {  DOMElement *textElem = (DOMElement *)elem->getFirstChild();    if ((textElem != NULL) && (textElem->getNodeType() == DOMNode::TEXT_NODE))    return CPLStrdup(XMLString::transcode(textElem->getNodeValue()));    return NULL;}char *getREFValue(DOMElement *elem) {    return CPLStrdup(XMLString::transcode(elem->getAttribute(XMLString::transcode(ILI2_REF))));}OGRPoint *getPoint(DOMElement *elem) {  // elem -> COORD (or ARC)  OGRPoint *pt = new OGRPoint();    DOMElement *coordElem = (DOMElement *)elem->getFirstChild();  while (coordElem != NULL) {    if (cmpStr("C1", XMLString::transcode(coordElem->getTagName())) == 0)      pt->setX(atof(getObjValue(coordElem)));    else if (cmpStr("C2", XMLString::transcode(coordElem->getTagName())) == 0)      pt->setY(atof(getObjValue(coordElem)));    else if (cmpStr("C3", XMLString::transcode(coordElem->getTagName())) == 0)      pt->setZ(atof(getObjValue(coordElem)));    coordElem = (DOMElement *)coordElem->getNextSibling();  }  pt->flattenTo2D();  return pt;}OGRLineString *ILI2Reader::getArc(DOMElement *elem) {  // elem -> ARC  OGRLineString *ls = new OGRLineString();  // previous point -> start point  OGRPoint *ptStart = getPoint((DOMElement *)elem->getPreviousSibling()); // COORD or ARC  // end point  OGRPoint *ptEnd = new OGRPoint();  // point on the arc   OGRPoint *ptOnArc = new OGRPoint();  double radius = 0; // radius    DOMElement *arcElem = (DOMElement *)elem->getFirstChild();  while (arcElem != NULL) {    if (cmpStr("C1", XMLString::transcode(arcElem->getTagName())) == 0)      ptEnd->setX(atof(getObjValue(arcElem)));    else if (cmpStr("C2", XMLString::transcode(arcElem->getTagName())) == 0)      ptEnd->setY(atof(getObjValue(arcElem)));    else if (cmpStr("C3", XMLString::transcode(arcElem->getTagName())) == 0)      ptEnd->setZ(atof(getObjValue(arcElem)));    else if (cmpStr("A1", XMLString::transcode(arcElem->getTagName())) == 0)      ptOnArc->setX(atof(getObjValue(arcElem)));    else if (cmpStr("A2", XMLString::transcode(arcElem->getTagName())) == 0)      ptOnArc->setY(atof(getObjValue(arcElem)));    else if (cmpStr("A3", XMLString::transcode(arcElem->getTagName())) == 0)      ptOnArc->setZ(atof(getObjValue(arcElem)));    else if (cmpStr("R", XMLString::transcode(arcElem->getTagName())) == 0)      radius = atof(getObjValue(arcElem));    arcElem = (DOMElement *)arcElem->getNextSibling();  }  ptEnd->flattenTo2D();  ptOnArc->flattenTo2D();  interpolateArc(ls, ptStart, ptOnArc, ptEnd, arcIncr);  return ls;}OGRLineString *getLineString(DOMElement *elem) {  // elem -> POLYLINE  OGRLineString *ls = new OGRLineString();    DOMElement *lineElem = (DOMElement *)elem->getFirstChild();  while (lineElem != NULL) {    if (cmpStr(ILI2_COORD, XMLString::transcode(lineElem->getTagName())) == 0)      ls->addPoint(getPoint(lineElem));    else if (cmpStr(ILI2_ARC, XMLString::transcode(lineElem->getTagName())) == 0) {      // end point      OGRPoint *ptEnd = new OGRPoint();      // point on the arc       OGRPoint *ptOnArc = new OGRPoint();      // radius      double radius = 0;            DOMElement *arcElem = (DOMElement *)lineElem->getFirstChild();      while (arcElem != NULL) {                if (cmpStr("C1", XMLString::transcode(arcElem->getTagName())) == 0)          ptEnd->setX(atof(getObjValue(arcElem)));        else if (cmpStr("C2", XMLString::transcode(arcElem->getTagName())) == 0)          ptEnd->setY(atof(getObjValue(arcElem)));        else if (cmpStr("C3", XMLString::transcode(arcElem->getTagName())) == 0)          ptEnd->setZ(atof(getObjValue(arcElem)));        else if (cmpStr("A1", XMLString::transcode(arcElem->getTagName())) == 0)          ptOnArc->setX(atof(getObjValue(arcElem)));        else if (cmpStr("A2", XMLString::transcode(arcElem->getTagName())) == 0)          ptOnArc->setY(atof(getObjValue(arcElem)));        else if (cmpStr("A3", XMLString::transcode(arcElem->getTagName())) == 0)          ptOnArc->setZ(atof(getObjValue(arcElem)));        else if (cmpStr("R", XMLString::transcode(arcElem->getTagName())) == 0)          radius = atof(getObjValue(arcElem));                arcElem = (DOMElement *)arcElem->getNextSibling();      }            ptEnd->flattenTo2D();      ptOnArc->flattenTo2D();      OGRPoint *ptStart = getPoint((DOMElement *)lineElem->getPreviousSibling()); // COORD or ARC      interpolateArc(ls, ptStart, ptOnArc, ptEnd, PI/180);    } /* else { // FIXME StructureValue in Polyline not yet supported    } */            lineElem = (DOMElement *)lineElem->getNextSibling();  }    return ls;}OGRLineString *getBoundary(DOMElement *elem) {    DOMElement *lineElem = (DOMElement *)elem->getFirstChild();  if (lineElem != NULL)    if (cmpStr(ILI2_POLYLINE, XMLString::transcode(lineElem->getTagName())) == 0)      return getLineString(lineElem);    return new OGRLineString;}OGRPolygon *getPolygon(DOMElement *elem) {  OGRPolygon *pg = new OGRPolygon();    DOMElement *boundaryElem = (DOMElement *)elem->getFirstChild(); // outer boundary  while (boundaryElem != NULL) {    if (cmpStr(ILI2_BOUNDARY, XMLString::transcode(boundaryElem->getTagName())) == 0)      pg->addRing((OGRLinearRing *)getBoundary(boundaryElem));            boundaryElem = (DOMElement *)boundaryElem->getNextSibling(); // inner boundaries  }    return pg;}OGRGeometry *ILI2Reader::getGeometry(DOMElement *elem, int type) {  OGRGeometryCollection *gm = new OGRGeometryCollection();  DOMElement *childElem = elem;  while (childElem != NULL) {    switch (type) {      case ILI2_COORD_TYPE :         if (cmpStr(ILI2_COORD, XMLString::transcode(childElem->getTagName())) == 0)          return getPoint(childElem);        break;      case ILI2_ARC_TYPE :        // is it possible here? It have to be a ARC or COORD before (getPreviousSibling)        if (cmpStr(ILI2_ARC, XMLString::transcode(childElem->getTagName())) == 0)          return getArc(childElem);        break;      case ILI2_POLYLINE_TYPE :        if (cmpStr(ILI2_POLYLINE, XMLString::transcode(childElem->getTagName())) == 0)          return getLineString(childElem);        break;      case ILI2_BOUNDARY_TYPE :        if (cmpStr(ILI2_BOUNDARY, XMLString::transcode(childElem->getTagName())) == 0)          return getLineString(childElem);        break;      case ILI2_AREA_TYPE :        if ((cmpStr(ILI2_AREA, XMLString::transcode(childElem->getTagName())) == 0) ||          (cmpStr(ILI2_SURFACE, XMLString::transcode(childElem->getTagName())) == 0))          return getPolygon(childElem);        break;      default :         if (type >= ILI2_GEOMCOLL_TYPE) {          int subType = getGeometryTypeOfElem(childElem); //????          gm->addGeometry(getGeometry(childElem, subType));        }        break;    }        // GEOMCOLL    childElem = (DOMElement *)childElem->getNextSibling();  }    return gm;}const char* ILI2Reader::GetLayerName(IOM_BASKET model, IOM_OBJECT table) {    static char layername[512];    IOM_OBJECT topic = GetAttrObj(model, table, "container");    layername[0] = '\0';    strcat(layername, iom_getattrvalue(GetAttrObj(model, topic, "container"), "name"));    strcat(layername, ".");    strcat(layername, iom_getattrvalue(topic, "name"));    strcat(layername, ".");    strcat(layername, iom_getattrvalue(table, "name"));    return layername;}void ILI2Reader::AddField(OGRLayer* layer, IOM_BASKET model, IOM_OBJECT obj) {  const char* typenam = "Reference";  if (EQUAL(iom_getobjecttag(obj),"iom04.metamodel.LocalAttribute")) typenam = GetTypeName(model, obj);  if (EQUAL(typenam, "iom04.metamodel.SurfaceType")) {  } else if (EQUAL(typenam, "iom04.metamodel.AreaType")) {  } else if (EQUAL(typenam, "iom04.metamodel.PolylineType") ) {  } else if (EQUAL(typenam, "iom04.metamodel.CoordType")) {  } else {    OGRFieldDefn fieldDef(iom_getattrvalue(obj, "name"), OFTString);    layer->GetLayerDefn()->AddFieldDefn(&fieldDef);    CPLDebug( "OGR_ILI", "Field %s: %s", fieldDef.GetNameRef(), typenam);  }}int ILI2Reader::ReadModel(char **modelFilenames) {  IOM_BASKET model;  IOM_ITERATOR modelelei;  IOM_OBJECT modelele;  iom_init();  // set error listener to a iom provided one, that just   // dumps all errors to stderr  iom_seterrlistener(iom_stderrlistener);