[gdal-dev] How to set BoundingBox within a shapefile using OGR
Pythons bindings
Daniele Romagnoli
daniele.romagnoli at geo-solutions.it
Wed Dec 16 12:50:59 EST 2009
Hi list,
I'm creating a shapefile using the OGR python bindings. (from the 1.6.1
installer)
I have some issues in creating the boundingBox for it.
This is a snippet of the code I'm using:
-----------------------------------------------------
from osgeo import ogr
from osgeo import osr
....
....
northings_llc = dd['northings_llc']
northings_urc = dd['northings_urc']
eastings_llc = dd['eastings_llc']
eastings_urc = dd['eastings_urc']
t_srs = osr.SpatialReference()
t_srs.SetFromUserInput('EPSG:32632')
drv = ogr.GetDriverByName('ESRI Shapefile')
ds = drv.CreateDataSource(fileout)
x0, y0, x1, y1 = str(eastings_llc[0][0]), str(northings_llc[0][0]),
str(eastings_urc[0][0]), str(northings_urc[0][0])
layer = ds.CreateLayer(ds.GetName(), geom_type = ogr.wkbPolygon, srs =
t_srs)
geom = ogr.Geometry(type = layer.GetLayerDefn().GetGeomType())
geom.AssignSpatialReference(t_srs)
wkt = 'POLYGON('+x0+','+y0+','+ x0+','+ y1+','+ x1+',' +y1+','+ x1+','
+y0+','+ x0+',' +y0+')'
print(wkt)
geom.importFromWkt(wkt)
...
...
When running, I get this output:
POLYGON(602759.984846,4739727.96408,602759.984846,4739720.3473,602760.45247,4739720.3473,602760.45247,4739727.96408,602759.984846,4739727.96408)
Traceback (most recent call last):
File "prova.py", line 42, in <module>
geom.importFromWkt(wkt)
File "c:\python26\lib\site-packages\osgeo\ogr.py", line 2672, in <lambda>
__getattr__ = lambda self, name: _swig_getattr(self, Geometry, name)
File "c:\python26\lib\site-packages\osgeo\ogr.py", line 34, in
_swig_getattr
raise AttributeError,name
AttributeError: importFromWkt
Indeed, looking at the ogr.py there is no "importFromWkt" in the Geometry
definition.
How can I set the bounding box within the geometry?
Any help would be greatly appreciated.
Best Regards,
Daniele
--
-------------------------------------------------------
Eng. Daniele Romagnoli
Software Engineer
GeoSolutions S.A.S.
Via Carignoni 51
55041 Camaiore (LU)
Italy
phone: +39 0584983027
fax: +39 0584983027
mob: +39 328 0559267
http://www.geo-solutions.it
-------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20091216/baacc4a4/attachment.html
More information about the gdal-dev
mailing list