[gdal-dev] can't handle grib file

2012-08-22 Thread Margherita Di Leo
Hi all, I have a problem handling some grib files, I get a warning message and I can't understand how to deal with it.. running gdalinfo, it says: Un-handled possible ensemble section center 78 subcenter 255 [... several times, cut] Un-handled possible ensemble section center 78 subcenter 255 Dr

Re: [gdal-dev] [Qgis-developer] Large raster (ecw) identify very long

2012-08-22 Thread Even Rouault
Le mercredi 22 août 2012 20:13:39, Radim Blazek a écrit : > Even, > thanks for exhaustive explanation and testing. > > On Wed, Aug 22, 2012 at 2:37 PM, Even Rouault > > wrote: > >> I found in GDAL ecwdataset.cpp that it is treating single row > > > >> requests in IRasterIO in a special way: >

Re: [gdal-dev] [Qgis-developer] Large raster (ecw) identify very long

2012-08-22 Thread Radim Blazek
Even, thanks for exhaustive explanation and testing. On Wed, Aug 22, 2012 at 2:37 PM, Even Rouault wrote: >> I found in GDAL ecwdataset.cpp that it is treating single row >> requests in IRasterIO in a special way: > > I tried the following Python script that must be representative of how QGIS >

Re: [gdal-dev] problem with CloseRings

2012-08-22 Thread Chaitanya kumar CH
Doug, You could rebuild individual polygons by attempting to close the constituent rings one-by-one while discarding those that give an error with the CloseRings() method. On Wed, Aug 22, 2012 at 10:41 PM, wrote: > > Oliver, > > Would it just be a matter of counting the vertices of each polygon

Re: [gdal-dev] cutting up large TIFF files with gdal_translate

2012-08-22 Thread Even Rouault
Le mercredi 22 août 2012 18:43:12, Frank Warmerdam a écrit : > On Wed, Aug 22, 2012 at 9:38 AM, Zoltan Szecsei wrote: > > Hi, > > I have a large image bounded by: > > > > Upper Left ( 447007.599, 7132849.501) ( 32d28'15.03"E, 25d55'19.12"S) > > Lower Left ( 447007.599, 7127851.511) ( 32d28'1

Re: [gdal-dev] problem with CloseRings

2012-08-22 Thread Doug_Newcomb
Oliver, Would it just be a matter of counting the vertices of each polygon, selecting the polygons with 3 vertices, and comparing the value of the first and last vertex as you are building wkt value string below? Doug >The problem with your rings isn't that they are not closed, it is that

Re: [gdal-dev] problem with CloseRings

2012-08-22 Thread Frank Warmerdam
Oliver, The problem with your rings isn't that they are not closed, it is that they are degenerate. For instance: (559560.19 142169.94,559550.64 142173.43,559560.19 142169.94) This goes from PointA to PointB back to PointA. It is closed, it just is degenerate. It has no area. I'm not aware of

Re: [gdal-dev] cutting up large TIFF files with gdal_translate - Solution

2012-08-22 Thread Zoltan Szecsei
Hi Frank, Thanks for the prompt reply and explanation. Kind regards, Zoltan On 2012/08/22 18:43, Frank Warmerdam wrote: On Wed, Aug 22, 2012 at 9:38 AM, Zoltan Szecsei wrote: Hi, I have a large image bounded by: Upper Left ( 447007.599, 7132849.501) ( 32d28'15.03"E, 25d55'19.12"S) Lower L

Re: [gdal-dev] cutting up large TIFF files with gdal_translate

2012-08-22 Thread Frank Warmerdam
On Wed, Aug 22, 2012 at 9:38 AM, Zoltan Szecsei wrote: > Hi, > I have a large image bounded by: > > Upper Left ( 447007.599, 7132849.501) ( 32d28'15.03"E, 25d55'19.12"S) > Lower Left ( 447007.599, 7127851.511) ( 32d28'14.30"E, 25d58' 1.59"S) > Upper Right ( 451005.901, 7132849.501) ( 32d30'38

[gdal-dev] cutting up large TIFF files with gdal_translate

2012-08-22 Thread Zoltan Szecsei
Hi, I have a large image bounded by: Upper Left ( 447007.599, 7132849.501) ( 32d28'15.03"E, 25d55'19.12"S) Lower Left ( 447007.599, 7127851.511) ( 32d28'14.30"E, 25d58' 1.59"S) Upper Right ( 451005.901, 7132849.501) ( 32d30'38.75"E, 25d55'19.63"S) Lower Right ( 451005.901, 7127851.511) ( 32

[gdal-dev] problem with CloseRings

2012-08-22 Thread Oliver Christen
dear community first of all please accept my apologies if this is not the correct place to ask about such question. I am trying to manipulate some geometric data using the gdal/ogr python lib and Im facing a problem when trying to transform invalid ring into valid ring. I have polygons with uncl

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Smith, Michael ERDC-RDE-CRREL-NH
Jukka, Actually, it doesn't matter since only one value is used in the query, min or max return the same result. Its just using the min function with the group by to aggregate the values into one row. All the other values are null. Mike -- Michael Smith Remote Sensing/GIS Center US Army Corps

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Rahkonen Jukka
Hi, This part of query takes UB min values as maxx and maxy. Wouldn't it be better to select the max values instead? min(case when r=1 then sdo_ub else null end) maxx, min(case when r=2 then sdo_ub else null end) maxy -Jukka- Smith, Michael > > Jukka, > > I'm responsible for the code the c

Re: [gdal-dev] [Qgis-developer] Large raster (ecw) identify very long

2012-08-22 Thread Even Rouault
Selon Radim Blazek : > On Wed, Aug 22, 2012 at 12:29 PM, haubourg > wrote: > > > > Radim Blazek-2 wrote > >> > >> > >> QGIS is using GDALRasterIO() which reads a single pixel on original > >> resolution. AFAIK, ECW is using tiles internally so it should be all > >> very fast. I can imagine 2 prob

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Rahkonen Jukka
Hi, If I do ogrinfo with "--debug on" against Oracle l can see updates for USER_SDO_GEOM_METADATA. Same happens if I do ogr2ogr with output to for example GML. For me is seems not to happen only with loading and updating. I have not an easy access to Oracle logs and I have studied only the OGR

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Stefano Iacovella
2012/8/22 Rahkonen Jukka > Stefano Iacovella wrote: > > > >> I guess that meaning is to check all the metadata rows from the > >> ALL_SDO_GEOM_METADATA view which are referring to the queried table. > >> There can be several rows created be different Oracle users. Next the > >> maximum BBOX is co

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Smith, Michael ERDC-RDE-CRREL-NH
Actually, thinking about it, ALL_SDO_GEOM_METADATA has to be used to get the metadata info otherwise cross schema access would not be possible. Eg I'm connected as user_a but want to read data from user_b. Normal Oracle grants would control whether reading (or other operation) is possible. Mike

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Smith, Michael ERDC-RDE-CRREL-NH
Jukka, I'm responsible for the code the calculates the min/max extents for inserting/updating the USER_SDO_GEOM_METADATA view. This is supposed to only be run via ogr2ogr when loading/updating tables (I didn;t write that part, just the query to get those extents). There shouldn't be any updating g

Re: [gdal-dev] [Qgis-developer] Large raster (ecw) identify very long

2012-08-22 Thread Radim Blazek
On Wed, Aug 22, 2012 at 12:29 PM, haubourg wrote: > > Radim Blazek-2 wrote >> >> >> QGIS is using GDALRasterIO() which reads a single pixel on original >> resolution. AFAIK, ECW is using tiles internally so it should be all >> very fast. I can imagine 2 problems: >> >> - the tiles in ECW file are

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Rahkonen Jukka
Stefano Iacovella wrote: >> I guess that meaning is to check all the metadata rows from the >> ALL_SDO_GEOM_METADATA view which are referring to the queried table. >> There can be several rows created be different Oracle users. Next the >> maximum BBOX is constructed and that is updated into >> U

Re: [gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Stefano Iacovella
2012/8/22 Jukka Rahkonen > Hi, > > I learned to use ogrinfo with --debug on lately and it revealed what > all happens after sending plain ogrinfo through OCI driver. A few > lines interest me. What is the meaning to doing the select from > ALL_SDO_GEOM_METADATA then insert into USER_SDO_GEOM_META

Re: [gdal-dev] [Qgis-developer] Large raster (ecw) identify very long

2012-08-22 Thread Radim Blazek
On Wed, Aug 22, 2012 at 11:40 AM, haubourg wrote: > Hi all, > I'm ready to support developpement now, so let's undig a topic discussed > previously here: > http://hub.qgis.org/issues/4594 http://hub.qgis.org/issues/4594 > > Identify for very large ecw (from 1 to 100 Go - 1 080 000 * 1 090 000 > pi

[gdal-dev] When and why ogrinfo sends update for Oracle?

2012-08-22 Thread Jukka Rahkonen
Hi, I learned to use ogrinfo with --debug on lately and it revealed what all happens after sending plain ogrinfo through OCI driver. A few lines interest me. What is the meaning to doing the select from ALL_SDO_GEOM_METADATA then insert into USER_SDO_GEOM_METADATA as follows? OCI: Prepare(select