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

2010-12-09 Thread Paolo Corti
# And now the error: a code 6 but no further explanation geom.Transform(coordTrans) Hi Gregor your code is correct. I think it may be related to your dataset/feature geometry: could you kindly provide the wkt of the geometry that is failing (print geom.ExportToWkt()) ? P -- Paolo Corti GIS

[gdal-dev] image-io extension

2010-12-09 Thread Viji V Nair
Hi, I am one of the GDAL co-maintainer for Fedora, I would like to have the image-io extension support, currently it been supported through the following patch from image-io. But I dont think this is a good idea as the modifications are bit heavy.

[gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread Jukka Rahkonen
iomeneandrei aborruso at tin.it writes: Hi all, is there in OGR a command that I can use to convert multipart polygons in singleparts? If it does not matter to use some other tool than OGR, OpenJUMP can explode multifeatures and geometry collections. I made a little test and it correctly

[gdal-dev] importFromEPSG fails in a VS2008 project

2010-12-09 Thread Helm, P.W. (Pim) van den
Hi all, I have downloaded gdal-1.7.3 and build this via VS2008 succesfully. Furthermore I created an evironment variable GDAL_DATA that points to my location of the folder gdal-1.7.3\data. I copied the gdal17.dll towards my release location of my VS2008 application that will use this dll and

Re: [gdal-dev] importFromEPSG fails in a VS2008 project

2010-12-09 Thread Ivan Lucena
Helm, Have you used the Microsoft Process Monitor. That could help. Ivan ---Original Message--- From: Helm, P.W. (Pim) van den pim.vandenh...@tno.nl To: gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.org Subject: [gdal-dev] importFromEPSG fails in a VS2008 project Sent: Dec 09

RE: [gdal-dev] Creating aux.xml files for GTiff using Python

2010-12-09 Thread Jason Roberts
Matt, Since nobody else more knowledgeable replied and I was in the original discussion, I'll take a shot. To review: for some situations and formats, the following code would work: statistics = band.GetStatistics(False, True) # Force calculation of statistics

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

2010-12-09 Thread Gregor at HostGIS
On 12/9/2010 12:34 AM, Paolo Corti wrote: dataset/feature geometry: could you kindly provide the wkt of the geometry that is failing (print geom.ExportToWkt()) ? Sure. It's a zip code area for Arizona, just some SHP I grabbed off National Atlas. POLYGON ((-110.347639068730317

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

2010-12-09 Thread Frank Warmerdam
On 10-12-08 07:40 PM, Gregor at HostGIS wrote: 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 ... sourceSR = osr.SpatialReference() sourceSR.ImportFromProj4('+proj=longlat +ellps=WGS84 +datum=WGS84

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

2010-12-09 Thread Howard Butler
On Dec 9, 2010, at 11:08 AM, Frank Warmerdam wrote: On 10-12-08 07:40 PM, Gregor at HostGIS wrote: 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 ... sourceSR = osr.SpatialReference()

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

2010-12-09 Thread Gregor at HostGIS
IIRC, older versions of the bindings silenced error messages to stdout but didn't throw exceptions by default. Gregor might try issuing 'osr.UseExceptions()' Ah, good advice hobu. I had set ogr.UseExceptions() in my real code, but not osr's version. The resulting error is slightly less

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

2010-12-09 Thread Gregor at HostGIS
On 12/9/2010 9:08 AM, Frank Warmerdam wrote: ERROR 1: latitude or longitude exceeded limits I'm not sure why you did not see the error message I am seeing. No idea here. Maybe an older version of this or that? I believe the problem is that your data (-110 longitude) is no where near the

[gdal-dev] basic gdal_translate question

2010-12-09 Thread Honey, Steve
Hi I'm new to gdal so please forgive me if this is a really dumb question. I'm running a file through the command line utility gdal_translate. Everything works as I expect except the corners of the equirectangular grid are in Northing/Easting. I would like it to output in degrees. If

Re: [gdal-dev] basic gdal_translate question

2010-12-09 Thread Gregor at HostGIS
I would like it to output in degrees. If there an option I can specify on gdal_translate to put it in degrees? Or is there another utility such as gdalwarp which will do the conversion. You said it: gdalwarp is how you would reproject a image to a different SRS, say from UTM to lon-lat. --

Re: [gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread iomeneandrei
Hi Jukka, thank you for your time and for your advice. Rahkonen Jukka wrote: If it does not matter to use some other tool than OGR, OpenJUMP can explode multifeatures and geometry collections. My goal is to make this transformation in Python, and I did not want to write either unnecessary

Re: [gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread Paolo Corti
My goal is to make this transformation in Python, and I did not want to write either unnecessary code or complicated code (for me). If I have a solid and interesting result, I will publish it here. Ciao Andrea this code will do the trick (note that it will create a polygon for each ring, so

Re: [gdal-dev] image-io extension

2010-12-09 Thread Even Rouault
Hi Viji, I maintain the GDAL Java bindings, but the image-io GDAL plugin is another project whose maintener is Daniele I think. He'll correct me if I'm wrong, but I saw his name on http://docs.codehaus.org/display/GEOTOOLS/ImageIO-EXT+GDAL ;-) . I've CC'ed him so he has a chance to comment and

[gdal-dev] GDALDeregister_GTiff()

2010-12-09 Thread David Hildebrand
I'm not sure if this is the right list to post this question or not. I have a MapServer application built with GDAL support but when I try to access a 3-band TIFF image I get the error GDAL: GDALDeregister_GTiff() called.\r showing up in my log file. Does anyone know what the

Re: [gdal-dev] GDALDeregister_GTiff()

2010-12-09 Thread Frank Warmerdam
On 10-12-09 04:06 PM, David Hildebrand wrote: I'm not sure if this is the right list to post this question or not. I have a MapServer application built with GDAL support but when I try to access a 3-band TIFF image I get the error GDAL: GDALDeregister_GTiff() called.\r showing up in my log

[gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread iomeneandrei
Caro Paolo, Paolo Corti wrote: Ciao Andrea this code will do the trick (note that it will create a polygon for each ring, so you will get a polygon for each hole as well) [...] writing in groups like this is always a surprise and a great pleasure. Thank you, thank you, thank you. I have

Re: [gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread Even Rouault
I have modified a little your code - https://gist.github.com/735330 - and I have this message: Polygon added. Traceback (most recent call last): File PaoloIsGreat.py , line 32, in module multipoly2poly(in_lyr, out_lyr) File PaoloIsGreat.py , line 8, in multipoly2poly for geom_part in

RE: [gdal-dev] GDALDeregister_GTiff()

2010-12-09 Thread David Hildebrand
Frank: I turned on CPL_DEBUG in my apache local configuration file and re-started the server. However, I am getting no other messages. Incidentally, I am posting messages to the mapserver-users list with respect to my problem of displaying 3-band images. Thanks for your help.

[gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread iomeneandrei
Hi Even, Even Rouault wrote: Geometry have been made iterable only since GDAL 1.7.0. You need to upgrade your GDAL version or modify a bit the script to make the iteration with the OGR API. instead of : for geom_part in geom: addPolygon(geom_part.ExportToWkb(), out_lyr)

Re: [gdal-dev] Re: Multipart to singlepart

2010-12-09 Thread Sebastian E. Ovide
Hi Andrea, Just a note of warning: if you have a point x,y that is inside the hole (so it doesn't belong the the polygon), and you run some point in polygon query (intersection for example) you will get that it is no part of any polygon but if you convert the ring in two polygons using this