[Gfoss] cml2shp.py

[Luca Casagrande]

Ciao Bud,

complimenti davvero per il lavoro.
Potremo svilupparlo come plugin per Qgis in modo da avere lo shape già 
caricato nella mappa..

Se ti va posso darti una mano.

Luca

Bud P. Bruegger ha scritto:
> Ciao a tutti,
>
> Ho finalmente messo mano su un file CXF e uno CML, tutte e due
> disponibile ugualmente dal portale per comuni del catasto.  Ho capito
> dal file (e non l'avevo capito dalla spec) che l'attributo e' gia'
> logicamente associato alla geometria.   Cosi il parsing e' molto piu'
> facile che pensavo.  
>
> Allora mi ho messo a scrivere un piccolo parser per CML che e' xml e
> cosi piu' facile che quello di CXF.  
>
> La versione attuale prende solo Bordi (incl. particelle e fabbricati).
> Fin ora lo mette in una rappresentazione in memoria.  Il prossimo passo
> sara' di scrivere su un shape file usando la shapelib di python.  Spero
> di mandare una versione che include questo piu' tardi oggi.
>
> Commenti molto benvenuti
>
> saluti
> -b
>
>
>
> #######################################################################
> ## =======
> ## cml2shp
> ## =======
> ##
> ## Converts a subset of content from a file in
> ## Cadastral Markup Language (CML) [1] to a
> ## shape file
> ##
> ## [1] CML spec:
> ## http://www.agenziaterritorio.it/servizi/comunieistituzioni/..
> ## ..fornitura_dati_catastali/specifica%20tecnica%20CML.doc
> ##
> ## Copyright:  Comune di Grosseto
> ## Author:     Bud P. Bruegger <bud at comune.grosseto.it>
> ## License:    GPL (any version)
> #######################################################################
>
> import operator
>
> import shapelib
> import elementtree.ElementTree as ET
>
> version = "0.1 14/12/2007"
>
> #fileName = "testdata/E202_000100.CMF"
> fileName = "testdata/vertisola.CMF"
>
> class Feature(object):
>     "a simple in memory rep of a feature"
>     def __init__(self, bordoType, mappaID, particellaID, valenza,
> esterconf, geom): self.bordoType = bordoType
>         self.mappaID = mappaID
>         self.particellaID = particellaID
>         self.valenza = valenza
>         self.esterconf = esterconf
>         self.geom = geom
>         #geometry is represented by a list of rings,
>         #where every ring is a list of x/y tuples
>         #the first ring is the main one,
>         #the others are islands
>
> #-------------------------------------
> def parseFabbriPart(codbo):
>     "parses a codbo of a fabbricato or particella"
>     bordoType = "PARTICELLA"
>     particellaID = codbo
>     if codbo[-1] == '+':
>         bordoType = "FABBRICATO"
>         particellaID = codbo[:-1]
>     return bordoType, particellaID
>         
> def codboParse(codbo, mapName):
>     "parses a codbo"
>     #returns: bordoType, particellaID
>     codbo = codbo.strip()
>     mappaID = mapName
>     if codbo == mapName:
>         bordoType = "MAPPA"
>         particellaID = ""
>     elif codbo == "STRADA":
>         bordoType = "STRADA"
>         particellaID = ""
>     elif codbo == "ACQUA":
>         bordoType = "ACQUA"
>         particellaID = ""
>     else:
>         bordoType, particellaID = parseFabbriPart(codbo)
>     return bordoType, particellaID
>
> def vertConsume(nVerts, vertList):
>     return vertList[:nVerts], vertList[nVerts:]
>
> def coordParse(nVert, vertIsolaList, coordStr):
>     noIsolaVerts = reduce(operator.add, vertIsolaList, 0)
>     vertexList = [subStr.split(',') for subStr in coordStr.split()]
>     geom = []
>     mainRing, remainingVerts = vertConsume((nVert - noIsolaVerts),
> vertexList) geom.append(mainRing)
>     for noIslandVerts in vertIsolaList:
>         islandRing, remainingVerts = vertConsume(noIslandVerts,
> remainingVerts) geom.append(islandRing)
>     return geom
>
>
> #-- main parsing ---------------------
> mapTree = ET.parse(fileName)
> mapRoot = mapTree.getroot()
> mapID = mapRoot.find("INFOMAPPA").get("nome")
> featureList = []
>
> #-- process BORDO elements ------------------
> for border in mapRoot.getiterator("BORDO"):
>     bordoType, particellaID = codboParse(border.get("codbo"), mapID)
>     valenza = border.get("valenza")
>     esterconf = border.get("esterconf")
>     gbordo = border.find("GBORDO")
>     #nIsole = gbordo.get("n.isole")
>     nVert = int(gbordo.get("n.vert"))
>     vertIsolaList = map(int, [i.text for i in gbordo.getiterator
> ("VERTISOLA")]) coordStr = gbordo.find("COORD").text
>     geom = coordParse(nVert, vertIsolaList, coordStr)
>     featureList.append(Feature(bordoType, mapID, particellaID,
>                                valenza, esterconf, geom))
>
>
>