Mercurial > dive4elements > river
view flys-backend/contrib/shpimporter/fixpoints.py @ 5202:646c154477fe
SCHEMA CHANGE: removed Waterlevel and WaterleveDifferences, added lookup table for wst
author | Tom Gottfried <tom@intevation.de> |
---|---|
date | Fri, 08 Mar 2013 12:45:48 +0100 |
parents | edba92337204 |
children | 71a24278322c |
line wrap: on
line source
try: from osgeo import ogr, osr except ImportError: import ogr, osr from importer import Importer import logging logger = logging.getLogger("Fixpoints") fixpoints_no_km_logged=False TABLE_NAME="fixpoints" PATH="Geodaesie/Festpunkte" NAME="Fixpoints" class Fixpoint(Importer): def getPath(self, base): return "%s/%s" % (base, PATH) def getTablename(self): return TABLE_NAME def getName(self): return NAME def isGeometryValid(self, geomType): return geomType == 1 def isShapeRelevant(self, name, path): return True def createNewFeature(self, featureDef, feat, **args): newFeat = ogr.Feature(featureDef) geometry = feat.GetGeometryRef() newFeat.SetGeometry(geometry) newFeat.SetField("name", args['name']) if self.IsFieldSet(feat, "river_id"): newFeat.SetField("river_id", feat.GetField("river_id")) else: newFeat.SetField("river_id", self.river_id) if self.IsFieldSet(feat, "KM"): newFeat.SetField("km", feat.GetFieldAsDouble("KM")) elif self.IsFieldSet(feat, "ELBE_KM"): newFeat.SetField("km", feat.GetFieldAsDouble("ELBE_KM")) else: if not fixpoints_no_km_logged: logger.error("Could not find KM attribute") fixpoints_no_km_logged = True return None if self.IsFieldSet(feat, "X"): newFeat.SetField("x", feat.GetFieldAsDouble("X")) if self.IsFieldSet(feat, "Y"): newFeat.SetField("y", feat.GetFieldAsDouble("Y")) if self.IsFieldSet(feat, "HPGP"): newFeat.SetField("HPGP", feat.GetField("HPGP")) return newFeat