Mercurial > dive4elements > river
view flys-backend/contrib/shpimporter/km.py @ 4794:a7d080347ac3
MINFO: Allow two methods for outlier test in SQ relation.
* Methods can be switched as option in conf.xml.
* Methods:
- Find outliers via multiples of the standard deviation.
- Grubbs (used in Fix-Analysis)
author | Raimund Renkert <rrenkert@intevation.de> |
---|---|
date | Fri, 11 Jan 2013 13:57:38 +0100 |
parents | 59ca5dab2782 |
children | f1c01fecf194 |
line wrap: on
line source
import ogr from importer import Importer TABLE_NAME="river_axes_km" PATH="Geodaesie/Flussachse+km" NAME="KMS" class KM(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 name == "km" def createNewFeature(self, featureDef, feat, **args): newFeat = ogr.Feature(featureDef) newFeat.SetGeometry(feat.GetGeometryRef()) 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.IsDoubleFieldSet(feat, "km"): newFeat.SetField("km", feat.GetFieldAsDouble("km")) elif self.IsDoubleFieldSet(feat, "KM"): newFeat.SetField("km", feat.GetFieldAsDouble("KM")) else: return None return newFeat