Re: [gdal-dev] Transforming Raster Points from GeoTIFF with GCPs

2021-12-10 Thread Yann-Sebastien Tremblay-Johnston
Simon, See the Transformer API, in particular https://gdal.org/api/gdal_alg.html?highlight=transformer#_CPPv419GDALTransformerFunc (or gdal.Transformer in Python). Sebastien On Fri, Dec 10, 2021 at 10:24 AM Simon Eves wrote: > So we have a GeoTiFF of Sentinel-1 imagery which `gdalinfo` reports

[gdal-dev] Clarifying AxisMappingStrategy for SpatialReference read from GeoJSON in gdal > 3

2020-04-28 Thread Yann-Sebastien Tremblay-Johnston
Hi, Tested using conda gdal 3.0.4 and gdal bundled with QGIS (3.0.4) >From the documented changes in gdal > 3 with respect to honouring coordinate order i.e. lat/lon order for EPSG:4326, I would expect that when from osgeo import ogr > ds = ogr.Open('test.geojson') > layer = ds.GetLayer() > srs

Re: [gdal-dev] Interpolation used for egm96_15.gtx when using RPC_DEM_APPLY_VDATUM_SHIFT

2020-04-17 Thread Yann-Sebastien Tremblay-Johnston
VRT with the correct compoundCRS. Interpolation of the geoid grid is done by PROJ itself, and this is > bilinear interpolation Awesome, good to know. Is there anyway of configuring this? Thanks for the quick response. Sebastien On Fri, Apr 17, 2020 at 10:09 AM Even Rouault wrote:

[gdal-dev] Interpolation used for egm96_15.gtx when using RPC_DEM_APPLY_VDATUM_SHIFT

2020-04-17 Thread Yann-Sebastien Tremblay-Johnston
Hi, I spent a bit of time learning about vertical datum shifts when performing RPC coordinate transformations using a DEM this week, and I was curious about how the egm96_15.gtx file is interpolated. I was able to find in gdal_rpc.cpp that the required vdatum shift is calculated using an OGR Coo

Re: [gdal-dev] IAU2000 support?

2015-09-24 Thread Yann
Thanks for the info both of you. It is a recurring issue for me (and my work with gdal/GRASS), please contact me when you are up to tricks on that, I will try to find time to give a hand. Cheers, Yann On 24/09/2015 14:26, Hare, Trent wrote: Yann, Yes there has been no movement on pushing

Re: [gdal-dev] IAU2000 support?

2015-09-24 Thread Yann
Thank you Even, at which version this OGC URN support was started? Cheers, Yann On 24/09/2015 11:45, Even Rouault wrote: Yann, I have recovered a old script I made about 1.5 year ago (and about 1+ version of Ubuntu away) to process Chandrayaan-M3 Moon data in IAU2000:30100 and it seems that

[gdal-dev] IAU2000 support?

2015-09-24 Thread Yann
. We had discussions earlier about planetary support in this list, among whom Trent was involved (IIRC). Does anybody have any update about the status change for the newest version of GDAL? For the record: Ubuntu 15.10 (Wily), v1.11.2 Thanks, Yann -- - Mobile: +33 (0)7 83 85 52 34

Re: [gdal-dev] PDS Reader error

2014-12-08 Thread Yann Chemin
Hi Trent, Just wondering, as even ISIS does not import it, do you know any guys actually reading the rings images? What software do they use? It might be worth just clarifying how they actually read them, even just for the sake of pds2isis and sanitizing the gdal driver ... Cheers, Yann On 9

Re: [gdal-dev] Update of EPSG database to v8.4

2014-12-04 Thread Yann Chemin
Hi all, just want to pick up the status of this, anything new/started somewhere I could help with? Trent, anything happened on your side about planetary proj support? I am going to have a nother look about it, once answers/comments reach. Cheers, Yann On 16 May 2014 at 04:46, Etienne Tourigny

Re: [gdal-dev] Motion: Adopt RFC 46 : GDAL/OGR unification

2014-05-19 Thread Yann Chemin
+1, much needed On 20/05/2014, Even Rouault wrote: > Hi, > > I think that the points raised in the discussion have been answered. I've > just > done a minor edit to the RFC text to mention RFC 36, as rightly suggested by > > Ivan. > > So: > > Motion : I move to adopt RFC 46: GDAL/OGR unification

Re: [gdal-dev] Update of EPSG database to v8.4

2014-05-15 Thread Yann Chemin
Hi Even, I would be interested to include it. But I see that the 1Mb uncompressed text is an issue. If it was OK for the additional Mb in the source code, what would be the place to look into gdal code to create the patch needed? Thanks, Yann PS: did not copy other lists as I am not a member

Re: [gdal-dev] Update of EPSG database to v8.4

2014-05-14 Thread Yann Chemin
Hi, is there planetary datum support in this new version (i.e. Moon 2000, or etc.)? Yann On 14/05/2014, Even Rouault wrote: > Hi, > > I've followed the update process of the EPSG SRS database to latest v8.4, > and > just committed the updated files into libgeotiff, GDAL

[gdal-dev] a small GDAL+wxPython RS tool (wxGIPE) in github

2014-01-30 Thread Yann Chemin
/blob/master/screenshot.png svn check out: svn co https://github.com/YannChemin/wxGIPE Yann -- ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] Accessing a 2D array stored in 1D by GDALRasterIO()

2013-12-12 Thread Yann Chemin
Hi, I loaded an image with: GDALRasterIO(hB,GF_Read,0,0,nX,nY,oneDarray,nX,nY,GDT_Float32,0,0); to run some statistics, however, a new requirement is to make stats per land-use, is there a way to get the index of a given pixel in the 1D array from a given set of row & column? Thank you,

Re: [gdal-dev] again .hdf to .tif

2009-03-26 Thread Yann Chemin
Yes Vincent, that rings a bell indeed. Thanks, Yann 2009/3/26 Vincent Schut : > Yann Chemin wrote: >> >> Hi there, >> >> I am using this one right now on Level 2 and 3 MODIS products >> >> --- >> #!/bin/bash

Re: [gdal-dev] again .hdf to .tif

2009-03-26 Thread Yann Chemin
subdataset#*=} > tempfile.txt out=$(sed 's/\(.*\):\(.*\):"\(.*\).A\(.*\)\.\h\(.*\)\.\(.*\)\.\(.*\).hdf":\(.*\):\(.*\)/\3\_\4\_h\5\_\8\_\9\.tif/' tempfile.txt) gdalwarp -of GTiff -t_srs EPSG:4326 ${subdataset#*=} processed/$out done done -

[gdal-dev] problem accessing rows of an image in parallel

2009-03-01 Thread Yann Chemin
: IReadBlock failed at X offset 0, Y offset 157 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 157 On the original standard serial code, all is fine. Thanks for any suggestions, Yann -- Yann Chemin Perso: http://www.freewebs.com/ychemin YiKingDo: http://yikingdo.unblog.fr

[gdal-dev] Re: satellite image processing in Python

2009-02-11 Thread Yann Chemin
(output_filename, result) del output, geot, proj still loading the created file makes a segfault in ReadArray it seems. I am using 1.5.2 for Lenny. Maybe should try 1.6 Yann for i in range (0,6): tmp = gdal.Open(F[i

[gdal-dev] Re: satellite image processing in Python

2009-02-10 Thread Yann Chemin
#SaveArray(result, output_filename, 'GTiff') #Empty memory objects del result, b, out #NDVI NIR = LoadFile('b4_ref.tif') Red = LoadFile('b3_ref.tif') ndvi = (NIR-Red)/1.0*(NIR+Red) ___ 2009/2/11 Yann Chemin : >

[gdal-dev] satellite image processing in Python

2009-02-10 Thread Yann Chemin
b) 2 - A segmentation fault is crashing Python when using LoadFile() on any one of the files generated in the loop. Thank you, Yann - # Landsat 7 ETM + image downlaoded from GLCF #ftp.glcf.umiacs.umd.edu/glcf/Landsat

[gdal-dev] reproject points using raster info

2009-01-09 Thread Yann Chemin
++ examples actually, and have difficulties to find a proper C example somewhere to do that. thanks for any help/suggestion, Yann ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] using GDALCreateCopy and changing only the datatype

2008-11-24 Thread Yann Chemin
not georeferenced. need a bit of reading indeed. Thanks, Yann 2008/11/25 Frank Warmerdam <[EMAIL PROTECTED]>: > Yann Chemin wrote: >> >> Hello, >> >> I would like to use GDALCreateCopy and change the datatype from int32 >> to float32. >> >>

Re: [gdal-dev] using GDALCreateCopy and changing only the datatype

2008-11-24 Thread Yann Chemin
file as int32 so far for allowing change... Yann 2008/11/24 Kolberg Sjur A <[EMAIL PROTECTED]>: > > Hi Yann, > > I had this same problem and found no other solution than to switch from > CreateCopy() to Create(). > > Good luck, > > Sjur K :-) > > >

[gdal-dev] using GDALCreateCopy and changing only the datatype

2008-11-24 Thread Yann Chemin
; hBandOut0 = GDALGetRasterBand(hDatasetOut0,1); thanks, Yann -- Yann Chemin International Rice Research Institute Office: http://www.irri.org/gis Perso: http://www.freewebs.com/ychemin YiKingDo: http://yikingdo.unblog.fr/ ___ gdal-dev mailing list gdal-dev

Re: [gdal-dev] loading many files with C API

2008-11-24 Thread Yann Chemin
ES+1]; float *pafScanline[MAXFILES+1]; int nXSize1, nYSize1; out = argv[1]; //Loading input files n_imgs = argc - 2; for (i=0;i: > Yann, Frank, > > Another thing to considerer is that any process you might want to do on that > file

[gdal-dev] loading many files with C API

2008-11-23 Thread Yann Chemin
Hello, I need to load about 400 files with the C API to work on each pixel temporal dimension. I would like to make a look to load them but fail so far... GDALDatasetH hDataset[MAXFILES]; GDALAllRegister(); GDALDriverH hDriver[MAXFILES]; GDALRasterBandH hBand[MAXFILES]; float *pafScanline[MAXFIL