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] Merging GeoTiff files

2015-04-03 Thread Newcomb, Doug
Steve, By merge, do you mean files that cover the same area, or merging adjoining images into a single larger image? Doug On Fri, Apr 3, 2015 at 10:54 AM, Steve Cagle wrote: > I have several GeoTiff files from a (commercial vendor) to merge. I have > used gdal, gdal_merge.py, the merge with QGI

[gdal-dev] Merging GeoTiff files

2015-04-03 Thread Steve Cagle
I have several GeoTiff files from a (commercial vendor) to merge. I have used gdal, gdal_merge.py, the merge with QGIS and get the same results.  The output file is the proper size but only contains the last input .tiff file.  I have struggled with this for a while now, pretty sure it's somethin

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] gdal_translate and PostGIS

2015-04-03 Thread Martin Landa
Hi Even, 2015-04-03 11:37 GMT+02:00 Even Rouault : > As the error message suggests, no, not at the time with the PostGISRaster > driver which basically only supports creation of a new raster through > duplicatin an existing postgisraster table thanks to SQL copy mechanisms. > > You have to use ras

Re: [gdal-dev] gdal_translate and PostGIS

2015-04-03 Thread Even Rouault
Le vendredi 03 avril 2015 11:31:09, Martin Landa a écrit : > Hi, > > is it possible to import raster data into PostGIS using gdal_translate? > > $ gdal_translate -of PostGISRaster dmt.tif PG:dbname=pokusnik > Input file size is 19615, 9 > > fails with > > ERROR 6: PostGISRasterDataset::Crea

[gdal-dev] gdal_translate and PostGIS

2015-04-03 Thread Martin Landa
Hi, is it possible to import raster data into PostGIS using gdal_translate? $ gdal_translate -of PostGISRaster dmt.tif PG:dbname=pokusnik Input file size is 19615, 9 fails with ERROR 6: PostGISRasterDataset::CreateCopy() only works on source datasets that are PostGISRaster Thanks, Martin