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

2016-07-15 Thread Even Rouault
Le vendredi 15 juillet 2016 19:32:14, François ALLAIN a écrit : > I am now trying to do this process as an in-memory process. > I coudl figure out how to load the .tif in-memory : > > # Create in-memory file and initialize it with the content > gdal.FileFromMemBuffer('/vsimem/tiffinmem', self.fil

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

2016-07-15 Thread François ALLAIN
I am now trying to do this process as an in-memory process. I coudl figure out how to load the .tif in-memory : # Create in-memory file and initialize it with the content gdal.FileFromMemBuffer('/vsimem/tiffinmem', self.file) # Open the in-memory file ds = gdal.Open('/vsimem/tiffinmem', gdal.GA_

Re: [gdal-dev] SPOT-7 and Worldview 3

2016-07-15 Thread Shawn Gong
Thank you Even. Shawn From: Even Rouault [even.roua...@spatialys.com] Sent: July-15-16 12:02 PM To: gdal-dev@lists.osgeo.org Cc: Shawn Gong Subject: Re: [gdal-dev] SPOT-7 and Worldview 3 Le vendredi 15 juillet 2016 17:29:19, Shawn Gong a écrit : > Hi all

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

2016-07-15 Thread Benjamin Deschamps
Created a ticket here: https://trac.osgeo.org/gdal/ticket/6594#ticket First ticket, so hopefully I followed protocol. Thanks for your help. On Fri, Jul 15, 2016 at 12:07 PM, Even Rouault wrote: > Le vendredi 15 juillet 2016 15:31:09, Benjamin Deschamps a écrit : > > Hi Even, > > > > The origin

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

2016-07-15 Thread Sam Wilson
Thank you Andrew! That did it. Cheers, Sam On 7/14/16, 11:12 PM, "Andrew C Aitchison" wrote: >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

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

2016-07-15 Thread Even Rouault
Le vendredi 15 juillet 2016 15:31:09, Benjamin Deschamps a écrit : > Hi Even, > > The original file is a Landsat-8 image, I can confirm that it has tags set > in the imd metadata domain with gdalinfo -mdd imd > LC80040242014245LGN00_B3.TIF but I'm creating a brand new raster, not > making a copy,

Re: [gdal-dev] SPOT-7 and Worldview 3

2016-07-15 Thread Even Rouault
Le vendredi 15 juillet 2016 17:29:19, Shawn Gong a écrit : > Hi all, > > > > I have downloaded and built gdal 2.1.1 > > Ran some tests on SPOT-7 and Worldview 3 > > > > SPOT-7 tiled dimap: the DIM_SPOT7_*.xml file was loaded fine. But the > georeferencing was wrong. The xml seems to indicat

[gdal-dev] SPOT-7 and Worldview 3

2016-07-15 Thread Shawn Gong
Hi all, I have downloaded and built gdal 2.1.1 Ran some tests on SPOT-7 and Worldview 3 SPOT-7 tiled dimap: the DIM_SPOT7_*.xml file was loaded fine. But the georeferencing was wrong. The xml seems to indicate ITRF2005 (epsg 4896) but it has no projection information in gdal. The individ

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

2016-07-15 Thread Benjamin Deschamps
Hi Even, The original file is a Landsat-8 image, I can confirm that it has tags set in the imd metadata domain with gdalinfo -mdd imd LC80040242014245LGN00_B3.TIF but I'm creating a brand new raster, not making a copy, so I'm not sure why the metadata is being transferred. The following basic exam

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

2016-07-15 Thread François ALLAIN
My bad... I finally got it working with your code : import os import sys from osgeo import gdal, ogr, osr geojson_filename = "my.geojson" driver = ogr.GetDriverByName("GeoJSON") vector_ds = gdal.OpenEx( geojson_filename, gdal.OF_VECTOR ) geotiff_filename = "my.tif" ds = gdal.Open(geotiff_filena

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

2016-07-15 Thread Antonio Falciano
Il 15/07/2016 02:13, Sam Wilson ha scritto: 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: