view flys-backend/contrib/shpimporter/shpimporter.py @ 2861:b0132e1b9719

Added further shape importers and added the option to reproject shapes during the import process. flys-backend/trunk@4342 c6561f87-3c4e-4783-a992-168aeb5c3f6f
author Ingo Weinzierl <ingo.weinzierl@intevation.de>
date Thu, 03 May 2012 14:19:52 +0000
parents bd9e76e0b55d
children 998f1a7dcfde
line wrap: on
line source
import ogr

import utils
from uesg  import UESG
from axis  import Axis
from km    import KM
from lines import Line
from fixpoints import Fixpoint
from buildings import Building
from crosssectiontracks import CrosssectionTrack
from floodplains import Floodplain
from boundaries import HydrBoundary, HydrBoundaryPoly

DBCONN='OCI:user/pass@host'
PATH='/path/to/Gewaesser/Elbe'
RIVER_ID=the_river_id
DEST_SRS=31467


def getImporters():
    return [
        Axis(DBCONN, RIVER_ID, DEST_SRS),
        KM(DBCONN, RIVER_ID, DEST_SRS),
        CrosssectionTrack(DBCONN, RIVER_ID, DEST_SRS),
        Line(DBCONN, RIVER_ID, DEST_SRS),
        Fixpoint(DBCONN, RIVER_ID, DEST_SRS),
        Building(DBCONN, RIVER_ID, DEST_SRS),
        Floodplain(DBCONN, RIVER_ID, DEST_SRS),
        HydrBoundary(DBCONN, RIVER_ID, DEST_SRS),
        HydrBoundaryPoly(DBCONN, RIVER_ID, DEST_SRS),
        UESG(DBCONN, RIVER_ID, DEST_SRS)
        ]


if __name__ == '__main__':
    importers = getImporters()

    types = {}

    for importer in importers:
        shapes = utils.findShapefiles(importer.getPath(PATH))
        print "Found %i Shapefiles" % len(shapes)

        for shpTuple in shapes:
            geomType = importer.walkOverShapes(shpTuple)
            try:
                num = types[geomType]
                types[geomType] = num+1
            except:
                types[geomType] = 1

    for key in types:
        print "%i x geometry type %s" % (types[key], key)

http://dive4elements.wald.intevation.org