Changes between Version 33 and Version 34 of BluePrint/GIS/ShapefileLayers

Timestamp:

05/17/13 18:50:57 (12 years ago)

Author:

Fran Boon

Comment:

--

BluePrint/GIS/ShapefileLayers

@@ -43 +43 @@
 
 === Import into native Tables ===
+SHP should be converted to an lxml.etree for import using XSLT.
 
-See below for code to convert SHP to an lxml.etree for import using XSLT
-
-Paul Porthouse [2013-04-27 at the [http://mapaction.org MapAction] Hackathon]:
-
-Our idea is to use a basic web form with a file upload (similar to the Import from !OpenStreetMap) where you can select which Shapefile to use. You should also be able to set basic information including database connection details. This will then upload the shape and use OGR2OGR to convert the Shapefile to geometry and import it into the specified database.
-
-A basic OGR2OGR which should accomplish the import is:
-{{{OGR2OGR -overwrite -f "PostgreSQL" PG:dbname=databaseconnection shapefile.sh}}}
-
-This could be called directly from Python by running it as a sub process.
-
-Alternatively, to use Python bindings, you could perform the import using GDAL which would allow more control over the import directly within Python, but should still handle the import into Postgres.
-
+Working sample script below has been started to be imported into the right place in Sahana here:
+* https://github.com/flavour/eden/blob/master/modules/s3/codecs/shp.py#L228
 
 {{{
@@ -74 +64 @@
     print "Usage: importshape.py filename"
     sys.exit(0)
-    
+
 shapefilename = sys.argv[1]
 
@@ -80 +70 @@
 
 # Create the datasource
-ds = ogr.Open( shapefilename )
+ds = ogr.Open(shapefilename)
 
 # Open the shapefile
@@ -86 +76 @@
     print "Open failed.\n"
     sys.exit(0)
-    
+
 # Get the layer and iterate through the features
 lyr = ds.GetLayer(0)
 
-root = etree.Element("shapefile",name=layername)
+root = etree.Element("shapefile", name=layername)
 
 for feat in lyr:
-        
-    featurenode = etree.SubElement(root,"feature")
-    
+
+    featurenode = etree.SubElement(root, "feature")
+
     feat_defn = lyr.GetLayerDefn()
-                
+
     for i in range(feat_defn.GetFieldCount()):
-                
+
         field_defn = feat_defn.GetFieldDefn(i)
-                        
-        fieldnode = etree.SubElement(featurenode,field_defn.GetName())
-                
+
+        fieldnode = etree.SubElement(featurenode, field_defn.GetName())
+
         if field_defn.GetType() == ogr.OFTInteger:
             fieldnode.text = str(feat.GetFieldAsInteger(i))
@@ -112 +102 @@
         elif field_defn.GetType() == ogr.OFTString:
             FieldString = str(feat.GetFieldAsString(i))
-            fieldnode.text = FieldString.decode(encoding='UTF-8',errors='strict')
-            
-    wktnode = etree.SubElement(featurenode,"wkt")
-        
+            fieldnode.text = FieldString.decode(encoding="UTF-8", errors="strict")
+
+    wktnode = etree.SubElement(featurenode, "wkt")
+
     geom = feat.GetGeometryRef()        
     wktnode.text = geom.ExportToWkt()    
-                
 
 # Test the etree object