[gdal-dev] ogr + python: find polygons in x,y

2010-12-08 Thread Sebastian E. Ovide
Hi All, Using Python + Gdal/OGR is there any efficient way of finding all the polygons that contains the point X,Y ? thanks ? -- Sebastian E. Ovide ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] gdal_translate: GDAL: Should not happen. Cannot find laea_grid.asc

2010-12-08 Thread Hermann Peifer
Hi, This is most likely a minor issue, but it looks to me that the above should not happen message always appears when gdal_translating a .vrt file that refers to an .asc file. There are 2 GDALClose() for the .asc. Hermann $ gdal_translate --debug on laea_grid.vrt laea_grid.tif \ -co

Re: [gdal-dev] ogr + python: find polygons in x,y

2010-12-08 Thread Jose Gomez-Dans
Hi, On 8 December 2010 09:38, Sebastian E. Ovide sebastian.ov...@gmail.com wrote: Using Python + Gdal/OGR is there any efficient way of finding all the polygons that contains the point X,Y ? This is my try. Not the most efficient, but you could always use a spatial index and a spatial filter:

[gdal-dev] Problem using Union - want to join 2 geometry

2010-12-08 Thread Namrata
Hi All, I am trying to use Union function in OGRGeometry. But the unified geometry is giving just a single point. So I am wondering what is going wrong?? Basically I want to join two geometry (can say disjoint) so as to have a single geometry. So the method can compute the close points of those

Re: [gdal-dev] gdal_translate: GDAL: Should not happen. Cannot find laea_grid.asc

2010-12-08 Thread Even Rouault
Herman, yes the warning was a minor annoyance but should not have caused any problem. I've just pushed as fix in SVN so it doesn't appear in that situation. It's linked to the caching mechanism behind VRTs that differ the opening of dataset until real I/O is required. There are some

Re: [gdal-dev] ogr + python: find polygons in x,y

2010-12-08 Thread Sebastian E. Ovide
Hi Jose, thanks for that. I've used SetSpatialFilterRect(x,y,x+0.01,y+0.01) On Wed, Dec 8, 2010 at 12:16 PM, Jose Gomez-Dans jgomezd...@gmail.comwrote: Hi, On 8 December 2010 09:38, Sebastian E. Ovide sebastian.ov...@gmail.com wrote: Using Python + Gdal/OGR is there any efficient way of

Re: [gdal-dev] Re: several No data value

2010-12-08 Thread Frank Warmerdam
On 10-12-06 09:53 AM, ludovic wrote: Frank, thanks for your response : On 3 déc, 22:05, Frank Warmerdamwarmer...@pobox.com wrote: GDAL only supports single nodata values as part of the data model. You can add arbitrary metadata describing other nodata values, but generally GDAL itself and

[gdal-dev] Python geometry.Transform() returns code 6

2010-12-08 Thread Gregor at HostGIS
Hey all. I am having a problem with the geometry.Transform() method in OGR for Python. Short version: geometry.Transform() returns a code 6. A code 6, according to the OGR source is simply OGRERROR II guess the rest is best explained by example: from osgeo import ogr, osr shapefile =

[gdal-dev] Re: Problem using Union - want to join 2 geometry

2010-12-08 Thread Namrata
Probably I think I must use some different method to merge all the polygons. Please help me finding out to do so. Thanks in advance! -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Problem-using-Union-want-to-join-2-geometry-tp5815141p5818134.html Sent from the GDAL -