Advertisement
Guest User

Untitled

a guest
Jul 27th, 2017
75
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.38 KB | None | 0 0
  1. ogr2ogr --config GDAL_DATA "C:OSGeo4W64sharegdal" -f "ESRI Shapefile" C:kmfeltshapeutm322017%%a_p.shp PG:"host=%Server_host% user=%Schema_user% dbname=%Database% password=sl_nab" -sql "SELECT * FROM kmfelt_o2_p where flatenr=%%a AND komid < 1800" -t_srs EPSG:%SRID% -geomfield geo
  2.  
  3. call python %prog%5_deletefields_p.py C:kmfeltshapeutm322017%%a_p.shp
  4.  
  5. ...
  6. # Set POINT-shapefile
  7. fn = r""+sys.argv[1]+"" # Read filename from input
  8. print('filnavn er ' +fn)
  9.  
  10. #get the driver
  11. driver = osgeo.ogr.GetDriverByName('ESRI Shapefile') # will select the driver for our shp-file
  12.  
  13. # open the data source
  14. dataSource = driver.Open(fn, 1) # 0 means read-only. 1 means writeable (read/write)
  15.  
  16. # get layer
  17. layer = dataSource.GetLayer()
  18. geotype = layer.GetGeomType()
  19. print(geotype)
  20.  
  21. ...
  22. if layer.GetFeatureCount() == 0:
  23. km = ogr.FieldDefn('KMNR', ogr.OFTInteger)
  24. km.SetWidth(4) #4 integer width
  25. print "Legger til egenskapen KMNR som siste kolonne"
  26. layer.CreateField(km)
  27. # After adding fields, get the layer definition
  28. featureDefn = layer.GetLayerDefn()
  29. # Create a feature from scratch
  30. feature = ogr.Feature(featureDefn)
  31. # Populate the attribute table with a default (dummy) value
  32. print "Setter KMNR lik 9999"
  33. feature.SetField("KMNR", 9999)
  34. # Finally, feature must be written to the layer
  35. layer.CreateFeature(feature)
  36. # Clean up memory
  37. feature.Destroy()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement