Re: [gdal-dev] gdal_calc.py floating point arithmetic

2016-07-14 Thread Andrew C Aitchison
On Fri, 15 Jul 2016, Sam Wilson wrote: Hi, I’m attempting to linearly interpolate two geotiff files using gdal_calc.py. The geotiff files are 3 hours apart, so for example hour 0 and hour 3. I’d like to end up with hourly data, so I need hours 1 and 2. Below are the commands I’m using: Ho

[gdal-dev] gdal_calc.py floating point arithmetic

2016-07-14 Thread Sam Wilson
Hi, I’m attempting to linearly interpolate two geotiff files using gdal_calc.py. The geotiff files are 3 hours apart, so for example hour 0 and hour 3. I’d like to end up with hourly data, so I need hours 1 and 2. Below are the commands I’m using: gdal_calc.py -A hour0.tif -B hour3.tif —outf

Re: [gdal-dev] Fwd: Problem with gdal.Rasterize utility : argument GDALDatasetShadow in wrapper_GDALRasterizeDestName

2016-07-14 Thread François ALLAIN
It did resolve this error, but now I've got this : ERROR 10: Pointer 'hDS' is NULL in 'GDALGetProjectionRef'. Warning 1: The input vector layer has a SRS, but the output raster dataset SRS is unknown. Ensure output raster dataset has the same SRS, otherwise results might be incorrect. ERROR 10: P

Re: [gdal-dev] Fwd: Problem with gdal.Rasterize utility : argument GDALDatasetShadow in wrapper_GDALRasterizeDestName

2016-07-14 Thread Even Rouault
Le jeudi 14 juillet 2016 21:13:24, François ALLAIN a écrit : > Hello, > > I am trying to convert a gdal_rasterize command to a gdal.Rasterize call > with the Python bindings > (Python 2.7.9 64bits, GDAL 2.1.0 compiled with libproj0 libgeos-dev on > Debian 8.3). > > The gdal_rasterize command wor

Re: [gdal-dev] IMD files being created by gdalwarp and driver.Create

2016-07-14 Thread Even Rouault
Benjamin, > > I am seeing .IMD files being created first by driver.Create() with the > GTIFF driver, and then also by gdalwarp (doing a simple reprojection), with > both input and output in GTIFF format. > > I've never seen these files before... Can I prevent them from being > written? If they

[gdal-dev] Fwd: Problem with gdal.Rasterize utility : argument GDALDatasetShadow in wrapper_GDALRasterizeDestName

2016-07-14 Thread François ALLAIN
Hello, I am trying to convert a gdal_rasterize command to a gdal.Rasterize call with the Python bindings (Python 2.7.9 64bits, GDAL 2.1.0 compiled with libproj0 libgeos-dev on Debian 8.3). The gdal_rasterize command works perfectly : gdal_rasterize -i -b 4 -burn 0 my.geojson my.tif With the in

[gdal-dev] IMD files being created by gdalwarp and driver.Create

2016-07-14 Thread Benjamin Deschamps
Hi, I am seeing .IMD files being created first by driver.Create() with the GTIFF driver, and then also by gdalwarp (doing a simple reprojection), with both input and output in GTIFF format. I've never seen these files before... Can I prevent them from being written? I'm using GDAL 2.0.2 with the

Re: [gdal-dev] Get Metadata from VRT file, store 4D Dataset into VRT

2016-07-14 Thread Even Rouault
Hannes, > > 1st Question: How can I get the metadata from each band e.g. the > "SourceFilename"? > I tried the following: > /band1=ds.GetRasterBand(1) > metadata_band1=band1.GetMetadata()/ > Unfortunately, there is not output, although /gdalinfo/ (outside python) > provides all the metadata of th

[gdal-dev] Error building 1.11.1 on MingW64 for 32-bit

2016-07-14 Thread William Kelly Magee
Hello all, Trying to build a 32-bit version of gdal v1.11.1 on Windows10 using Mingw64 gcc 4.9.1 and getting an unresolved error at link time; C:/gcc-4.9.1/bin/../lib/gcc/i686-w64-mingw32/4.9.1/../../../../i686-w64-mingw32/lib/../lib/libiconv.a(localcharset.o):localcharset.c:(.text+0x7):

Re: [gdal-dev] Error in GDALWarp to NWT_GRD

2016-07-14 Thread Even Rouault
Le jeudi 14 juillet 2016 13:07:42, jramm a écrit : > I added the following to the end of the Create method in > frmts/northwood/grddataset.cpp: > > > vsi_l_offset nFileSize = 1024 + nXSize * nYSize * 2; --> beware of the potential int32 overflow in nXSize * nYSize > if (VSIFTruncateL(po

Re: [gdal-dev] Error in GDALWarp to NWT_GRD

2016-07-14 Thread jramm
I added the following to the end of the Create method in frmts/northwood/grddataset.cpp: vsi_l_offset nFileSize = 1024 + nXSize * nYSize * 2; if (VSIFTruncateL(poDS->fp, nFileSize) != 0) { CPLError(CE_Failure, CPLE_FileIO, "Failed to allocat

[gdal-dev] Get Metadata from VRT file, store 4D Dataset into VRT

2016-07-14 Thread hannes.diedr...@gfz-potsdam.de
Hello, I just started to use virtual raster bands with gdal (python) and have two questions (see below). My aim: I like to manage Sentinel 2 data with VRB preferably in an 4D environment (x,y,time, spectral band) This is what I have done so far: I managed to build a vrt file by adding multiple j

Re: [gdal-dev] GDAL/OGR 1.11.5, 2.0.3 and 2.1.1 released

2016-07-14 Thread Homme Zwaagstra
Hi, For Docker users the releases are packaged up as Docker images at . These are based on the GDAL TravisCI test environment and as such are compiled with a broad range of drivers - the following commands list the available raster formats for the thre