GetFeature ( i ) # Create output Feature outFeature = ogr. GetFeatureCount ()): # Get the input Feature inFeature = inLayer. GetLayerDefn () # Add features to the ouput Layer for i in range ( 0, inLayer. CreateField ( fieldDefn ) # Get the output Layer's Feature Definition outLayerDefn = outLayer. GetFieldCount ()): fieldDefn = inLayerDefn. GetLayerDefn () for i in range ( 0, inLayerDefn. wkbPoint ) # Add input Layer Fields to the output Layer inLayerDefn = inLayer. CreateLayer ( "states_centroids", geom_type = ogr. CreateDataSource ( outShapefile ) outLayer = outDataSource. DeleteDataSource ( outShapefile ) # Create the output shapefile outDataSource = outDriver. GetDriverByName ( "ESRI Shapefile" ) # Remove output shapefile if it already exists if os. GetLayer () # Create the output Layer outShapefile = "states_centroids.shp" outDriver = ogr. Open ( inShapefile, 0 ) inLayer = inDataSource. GetDriverByName ( "ESRI Shapefile" ) inDataSource = inDriver. CreateFeature ( feature ) feature = None # Save and close DataSource inDataSource = None outDataSource = Noneįrom osgeo import ogr import os # Get the input Layer inShapefile = "states.shp" inDriver = ogr. CreateField ( idField ) # Create the feature and set values featureDefn = outLayer. wkbPolygon ) # Add an ID field idField = ogr. CreateLayer ( "states_convexhull", geom_type = ogr. ConvexHull () # Save extent to a new Shapefile outShapefile = "states_convexhull.shp" outDriver = ogr. GetGeometryRef ()) # Calculate convex hull convexhull = geomcol. wkbGeometr圜ollection ) for feature in inLayer : geomcol. GetLayer () # Collect all Geometry geomcol = ogr. CreateFeature ( feature ) feature = None # Save and close DataSource inDataSource = None outDataSource = Noneįrom osgeo import ogr import os # Get a Layer inShapefile = "states.shp" inDriver = ogr. CreateLayer ( "states_extent", geom_type = ogr. AddGeometry ( ring ) # Save extent to a new Shapefile outShapefile = "states_extent.shp" outDriver = ogr. GetExtent () # Create a Polygon from the extent tuple ring = ogr. exit ( 'ERROR: can not find layer in service' ) else : passįrom osgeo import ogr import os # Get a Layer's Extent inShapefile = "states.shp" inDriver = ogr.
GetLayerByName ( "largelayer" ) if not layer : sys. feat = None # Get a specific layer layer = wfs_ds. GetNextFeature () while feat is not None : feat = layer. ExportToWkt ()) # iterate over features feat = layer. GetSpatialRef () print 'Layer: %s, Features: %s, SR: %s. exit ( 'ERROR: can not open WFS datasource' ) else : pass # iterate over available layers for i in range ( wfs_ds. Open ( 'WFS:' + url ) if not wfs_ds : sys. SetConfigOption ( 'OGR_WFS_PAGE_SIZE', '10000' ) # Open the webservice url = '' wfs_ds = wfs_drv. SetConfigOption ( 'OGR_WFS_PAGING_ALLOWED', 'YES' ) gdal. Works on WFS 2.0 services and WFS 1.0 and 1.1 with some other services. SetConfigOption ( 'OGR_WFS_LOAD_MULTIPLE_LAYER_DEFN', 'NO' ) # Set config for paging.
GetDriverByName ( 'WFS' ) # Speeds up querying WFS capabilities for services with alot of layers gdal.
How to create a ring around a point in mapinfo 10 driver#
exit ( 'ERROR: cannot find GDAL/OGR modules' ) # Set the driver (optional) wfs_drv = ogr. Import sys try : from osgeo import ogr, osr, gdal except : sys. stderr, ': you must pass at least one argument - the layer name argument' sys. GetPrecision () print fieldName + " - " + fieldType + " " + str ( fieldWidth ) + " " + str ( GetPrecision ) # Close connection conn = None if _name_ = '_main_' : if len ( sys. GetFieldTypeName ( fieldTypeCode ) fieldWidth = lyrDefn. GetLayerDefn () for i in range ( lyrDefn. stderr, ': layer name = " %s " could not be found in database " %s "' % ( lyr_name, databaseName ) sys. GetLayer ( lyr_name ) if lyr is None : print > sys. From osgeo import ogr import sys databaseServer = "" databaseUser = "" databasePW = "" connString = "PG: host= %s dbname= %s user= %s password= %s " % ( databaseServer, databaseName, databaseUser, databasePW ) def GetPGLayerFieldTypes ( lyr_name ): conn = ogr.