Re: [gdal-dev] Extract Extent or Bounding Box for a subset of shapes for a shapefile

2011-04-22 Thread Marius Jigmond
Luca, You can also loop through the features with relatively simple Python code to get the extents. See sample below (adapt at will): def findPoints(geometry, results): for i in range(geometry.GetPointCount()): x,y,z = geometry.GetPoint(i) if results['north'] == None or results['north']

Re: [gdal-dev] Polygonize from C# ?

2011-04-22 Thread Luiz Motta
Jay, I made the plugin for QGIS, where, calculate the boundary of image with pixel value upper to zero, i believe is what you want, "footprint shape". The plugin read all images of directories and create one vector layer with all boundary of images. For calculate, the background of image need h

Re: [gdal-dev] Polygonize from C# ?

2011-04-22 Thread Elijah Robison
Hi Jay, you might want to review the following thread (Get extent from rasters). Based on a cursory read of his advice, I believe Trent offers a solution that might accommodate your task: http://lists.maptools.org/pipermail/fwtools/2010-October/001886.html /Elijah Jay Jennings wrote: Corr

RE: [gdal-dev] Polygonize from C# ?

2011-04-22 Thread Jay Jennings
Correct, we just need the footprint of the "good" (non-NODATA) area, which will often approximate a simple 4-point polygon with no holes. -Original Message- From: Elijah Robison [mailto:eli...@villagis.com] Sent: Friday, April 22, 2011 2:07 PM To: Jay Jennings Cc: gdal-dev@lists.osgeo.or

Re: [gdal-dev] Polygonize from C# ?

2011-04-22 Thread Elijah Robison
Jay, do you basically want a footprint shape of the useful data in the Sat images (say, without the collar), or is the problem more complex, where perhaps you have a mosaic of several images and, like swiss cheese, there are no-data holes within the image? It sounds to me like you just want a

RE: [gdal-dev] Polygonize from C# ?

2011-04-22 Thread Jay Jennings
I see, thanks. But since I have satellite images, this part of the doc is alarming: "...For non-thematic raster data (such as satellite images) the result will essentially be one small polygon per pixel, and memory and output layer sizes will be substantial. The algorithm is primarily intended

Re: [gdal-dev] gdal2tiles.py issue

2011-04-22 Thread Jerl Simpson
Thank you thank you thank you, GREATLY appreciated! Typed with thumbs... On Apr 22, 2011 12:25 PM, "Klokan Petr Pridal" wrote: > Hi Jerl, > > please read this issue report: > > http://code.google.com/p/maptiler/issues/detail?id=6 > > It gives you also explanation and also workaround for this pro

Re: [gdal-dev] Polygonize from C# ?

2011-04-22 Thread Even Rouault
Selon Jay Jennings : Generally, refering to the C/C++ doc is a good way of finding the info : See http://gdal.org/gdal__alg_8h.html#3f522a9035d3512b5d414fb4752671b1 > Hi list, > Can anyone point to any guidance or (preferably) examples on using > GDAL.Polygonize() with C# bindings ? I see what

[gdal-dev] Polygonize from C# ?

2011-04-22 Thread Jay Jennings
Hi list, Can anyone point to any guidance or (preferably) examples on using GDAL.Polygonize() with C# bindings ? I see what the interface is, but there are some unexplained parameters, e.g. "string[] options". My eventual goal is to open a raster and produce a WKT POLYGON string enclosing the

Re: [gdal-dev] Extract Extent or Bounding Box for a subset of shapes for a shapefile

2011-04-22 Thread Luca Sigfrido Percich
Thank you Chaitanya. I think that this makes great sense when dealing with a summary of the whole layer, especially because a lot of formats store the extent in the layer metadata/header. But if I need to dump all the features, or to get a dump or summary of a filtered set of features (in which c

Re: [gdal-dev] Extract Extent or Bounding Box for a subset of shapes for a shapefile

2011-04-22 Thread Chaitanya kumar CH
Sig, ogrinfo doesn't force the recomputation of extent. It is considered as a costly operation. ogrinfo just reports the stored extent. Whereas in ogr2ogr, creating a new file automatically computes the extents. On Fri, Apr 22, 2011 at 7:44 PM, Luca Sigfrido Percich wrote: > > Hi Matthew, > > Yo

Re: [gdal-dev] Extract Extent or Bounding Box for a subset of shapes for a shapefile

2011-04-22 Thread Luca Sigfrido Percich
Hi Matthew, Your question is interesting. ogrinfo -where seemed to be the answer, but... I tried ogrinfo with the -where "sql where" option with both MapInfo and PostGIS layers. In your case it should be: ogrinfo -ro -so -al -where "$WHERECLAUSE" "whole_shapefile.shp" The returned extent is al

[gdal-dev] Extract Extent or Bounding Box for a subset of shapes for a shapefile

2011-04-22 Thread Matthew Pettis
Hi All, I think this is the right list for the question below, as I cannot find a different list that matches my question better. I'm using ogr2ogr and ogrinfo (v1.4), and I want to extract the extent from the result of a subsetting where argument. I have a solution here that works but I'm wonde

Re: [gdal-dev] cropping a tif with lat/lon

2011-04-22 Thread Chaitanya kumar CH
Dan, Can you tell what your error message was? On Thu, Apr 21, 2011 at 1:41 AM, spalt wrote: > Hi there! > > Please forgive this terribly newbie question. Been searching all day and I > just don't get it. > > I have a geotiff VFR image from the FAA with the following coordinates as > reported