Re: [gdal-dev] ReadAsArray

2015-04-07 Thread Rutger
oid the inefficiency. I'll try to see if increasing the CACHE_MAX might do this 'on the fly' without me having to worry about the alignment of the block definitions. Regards, Rutger -- View this message in context: http://osgeo-org.1560.x6.nabble.com/gdal-dev-

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Stefan Keller
Hi Yuta Thanks! I'll try this. I'm still wondering which approach is faster: With gdal or rasterio... Cheers, Stefan 2015-04-03 17:46 GMT+02:00 Yuta Sato : > > > On Fri, Apr 3, 2015 at 11:08 PM, Stefan Keller wrote: >> >> Hi Yuta and Even >> >> Excuse me when I creep in that discussion - but I

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Yuta Sato
On Fri, Apr 3, 2015 at 11:08 PM, Stefan Keller wrote: > Hi Yuta and Even > > Excuse me when I creep in that discussion - but I'm acutally having > quite a similar beginners question (I'm rather a vector guy :-) ): > > Yuta: Can you post the code fragment here in order to make benchmarks? > > I h

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Stefan Ziegler
Hi Stefan You propably want to take a look at rasterio: https://github.com/mapbox/rasterio Which is also great to work with python and raster data. Regards Stefan Von meinem iPhone gesendet > Am 03.04.2015 um 16:08 schrieb Stefan Keller : > > Hi Yuta and Even > > Excuse me when I creep in th

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Stefan Keller
Hi Yuta and Even Excuse me when I creep in that discussion - but I'm acutally having quite a similar beginners question (I'm rather a vector guy :-) ): Yuta: Can you post the code fragment here in order to make benchmarks? I have a GeoTIFF (representing heights from SRTM3) in Mercator CRS of the

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Yuta Sato
Thank you very much "Even Rouault" for making me understood. On Fri, Apr 3, 2015 at 8:13 PM, Even Rouault wrote: > Le vendredi 03 avril 2015 12:46:48, Yuta Sato a écrit : > > Dear Even Rouault: > > > > Thank you very much. > > What about setting these parameters "used with .ReadAsArray()", thoug

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Even Rouault
Le vendredi 03 avril 2015 12:46:48, Yuta Sato a écrit : > Dear Even Rouault: > > Thank you very much. > What about setting these parameters "used with .ReadAsArray()", though I > did not know their meanings? > buf_xsize=None, buf_ysize=None, buf_obj=None (Please keep the list CC'ed) buf_xsize an

Re: [gdal-dev] ReadAsArray

2015-04-03 Thread Even Rouault
Le vendredi 03 avril 2015 12:22:00, Yuta Sato a écrit : > Dear Respected GDAL Developers and Users: > > What parameters should I set beforehand in order to accelerate the reading > of a GeoTiff file? > > I am using as follows: > > data = src_dataset.GetRasterBand(1).ReadAsArray(xoff,yoff,xsize,y

[gdal-dev] ReadAsArray

2015-04-03 Thread Yuta Sato
Dear Respected GDAL Developers and Users: What parameters should I set beforehand in order to accelerate the reading of a GeoTiff file? I am using as follows: data = src_dataset.GetRasterBand(1).ReadAsArray(xoff,yoff,xsize,ysize) Does setting the following parameters accelerate? GDAL_CACHEMAX,

Re: [gdal-dev] ReadAsArray: putting data in an already existent array

2013-03-23 Thread Even Rouault
Le vendredi 22 mars 2013 18:47:49, Diego Gnesi Bartolani a écrit : > Hi, > > I'm programming in Python and I need to "scan" a large image with GDAL > reading n rows at the time. I'm using the ReadAsArray method, but every > time I execute it I get a new array. I think this solution consumes a

Re: [gdal-dev] ReadAsArray: putting data in an already existent array

2013-03-23 Thread Jan Heckman
Hi Diego, I haven't checked this out, but my guess is that you should make a (python)list of the arrays if you need to keep them. If you don't need to keep them, you might try to del(ete) myarray before the next ReadAsArray call. Of course there is a lower level function which accepts an address a

[gdal-dev] ReadAsArray: putting data in an already existent array

2013-03-22 Thread Diego Gnesi Bartolani
Hi, I'm programming in Python and I need to "scan" a large image with GDAL reading n rows at the time. I'm using the ReadAsArray method, but every time I execute it I get a new array. I think this solution consumes a lot of resources. Is there a method or option to pass to ReadAsArray an array

Re: [gdal-dev] ReadAsArray with Windows 7

2011-10-21 Thread Tamas Szekeres
2011/10/21 Frank Warmerdam > > If the problem with your script persists when using the OSGeo4W > python environment then I should be able to reproduce it with a > suitable bug report. If it does not happen with OSGeo4W then > the issue likely relates to Tamas's use of SWIG version for the > pytho

Re: [gdal-dev] ReadAsArray with Windows 7

2011-10-21 Thread Tamas Szekeres
2011/10/21 Frank Warmerdam > On Fri, Oct 21, 2011 at 6:34 AM, Discourse Maps > wrote: > > Hi, > > > > I was wondering if the issue mentioned in a past thread has every been > > addressed. > > > http://osgeo-org.1803224.n2.nabble.com/Gdal-Python-Win7-crash-with-ReadAsArray-td4398401.html > > > >

Re: [gdal-dev] ReadAsArray with Windows 7

2011-10-21 Thread Frank Warmerdam
Oct 21, 2011 at 10:36 AM, Discourse Maps wrote: > Sure.  Though I'm not too saavy with the open source programming. > How do I go about acquiring and installing OSGeo4W binaries? > >> Date: Fri, 21 Oct 2011 10:14:39 -0700 >> Subject: Re: [gdal-dev] ReadAsArray with Wind

Re: [gdal-dev] ReadAsArray with Windows 7

2011-10-21 Thread Frank Warmerdam
On Fri, Oct 21, 2011 at 6:34 AM, Discourse Maps wrote: > Hi, > > I was wondering if the issue mentioned in a past thread has every been > addressed. > http://osgeo-org.1803224.n2.nabble.com/Gdal-Python-Win7-crash-with-ReadAsArray-td4398401.html > > I am having much the same issues with GDAL 1.6.1,

[gdal-dev] ReadAsArray with Windows 7

2011-10-21 Thread Discourse Maps
Hi, I was wondering if the issue mentioned in a past thread has every been addressed. http://osgeo-org.1803224.n2.nabble.com/Gdal-Python-Win7-crash-with-ReadAsArray-td4398401.html I am having much the same issues with GDAL 1.6.1, Windows 7, and using 'ReadAsArray'. Python crashes almost ever

Re: [gdal-dev] ReadAsArray misbehaving...

2010-09-06 Thread Riaan van den Dool
Even Thanks for pointing that out. I have now been able to solve that problem by keeping the dataset in scope... Riaan On Fri, Sep 3, 2010 at 7:52 PM, Even Rouault wrote: > Riann, > > hum, I think you've hit a classical problem with the use of the GDAL python > bindings. The root of your proble

Re: [gdal-dev] ReadAsArray misbehaving...

2010-09-03 Thread Even Rouault
Riann, hum, I think you've hit a classical problem with the use of the GDAL python bindings. The root of your problem is likely the following line: reader = RasterReaderWriter(gdal.Open(iFilename).GetRasterBand(1)) Currently the GDAL python bindings do a poor job with life cycle management of

Re: [gdal-dev] ReadAsArray misbehaving...

2010-09-03 Thread Riaan van den Dool
Thanks for the suggestion. I have tried: tileRaster = self.raster.ReadAsArray(xoff=offsetX, yoff=offsetY, win_xsize=width, win_ysize=height) and tileRaster = self.raster.ReadAsArray(xoff=offsetX, yoff=offsetY, buf_xsize=width, buf_ysize=height) but the result is something like: (array([],

[gdal-dev] ReadAsArray misbehaving...

2010-09-03 Thread Riaan van den Dool
I have recently started a new open-source project that will use gdal extensively, named scikits.eartho. The vision is to implement some advanced algorithms and ideas that we develop and work with at the South African Satellite Applications Centre (SAC) for use in Python (SciPy). I am experiencing