Re: [gdal-dev] Help with overlapping images in tiles generated by gdal2tiles.py

2024-12-04 Thread Rodriguez, Jose M CTR (USA) via gdal-dev
Hello, Yes my images (single band greyscale) come from the Clementine mission and are processed first with a) four corner points as gcps are inserted in the TIF file and b) warped so that the image's upper raster is aligned west to east: gdal_translate -q -of GTiff -gcp 1 1 $UL_LON $UL_LAT

Re: [gdal-dev] Help with overlapping images in tiles generated by gdal2tiles.py

2024-12-04 Thread Laurențiu Nicola via gdal-dev
Hi Jose, I'm not sure if this is happening in your case, but I've seen something similar with rotated rasters that used 0 or a negative value for the outside pixels. If you're doing interpolation or reprojection, those values will be used in the computation and affect neighboring pixels. But if

[gdal-dev] Help with overlapping images in tiles generated by gdal2tiles.py

2024-12-04 Thread Rodriguez, Jose M CTR (USA) via gdal-dev
Hello, I have been using GDAL to process a large collection of images of the Moon with the ultimate goal of generating a tile set for use in CesiumJS. I use gdalbuildvrt to build of virtual mosaic of the images in question and then use gdal2tiles.py to generate the tile set. All these steps wor

Re: [gdal-dev] Help using gdal_grid when src_datasource is a NetCDF File

2024-08-22 Thread Mike Taves via gdal-dev
On Fri, 23 Aug 2024 at 04:03, Ewan Tarrant via gdal-dev wrote: > Thanks in advance for your help, To read point data from a netCDF file, refer to the vector driver: https://gdal.org/drivers/vector/netcdf.html The data would need to conform to CF-1.8 for Point data: https://cfconventions.org/Data

[gdal-dev] Help using gdal_grid when src_datasource is a NetCDF File

2024-08-22 Thread Ewan Tarrant via gdal-dev
Hello gdal community! I would like to use gdal_grid (https://gdal.org/programs/gdal_grid.html) to rasterize scattered data stored in a NetCDF file. My scattered data consists of a single scalar value, z, defined at an Easting and Northing coordinate pair. I also have an EPSG code. I am having

Re: [gdal-dev] help

2024-01-25 Thread Abel Pau via gdal-dev
au via gdal-dev Enviado el: dijous, 25 de gener de 2024 15:58 Para: Marián Hadri mailto:marian.ha...@rib-software.com>> CC: gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> Asunto: Re: [gdal-dev] help Hi,. I have in my windows environment variables (not in any other place) this one: P

Re: [gdal-dev] help

2024-01-25 Thread Marián Hadri via gdal-dev
dal-dev mailto:gdal-dev-boun...@lists.osgeo.org>> En nombre de Abel Pau via gdal-dev Enviado el: dijous, 25 de gener de 2024 15:58 Para: Marián Hadri mailto:marian.ha...@rib-software.com>> CC: gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> Asunto: Re: [gdal-dev] hel

Re: [gdal-dev] help

2024-01-25 Thread Abel Pau via gdal-dev
with this one: cmake --build D:\GitHub-repository\GDAL\build I hope being helpful De: gdal-dev En nombre de Abel Pau via gdal-dev Enviado el: dijous, 25 de gener de 2024 15:58 Para: Marián Hadri CC: gdal-dev@lists.osgeo.org Asunto: Re: [gdal-dev] help Hi,. I have in my windows environment

Re: [gdal-dev] help

2024-01-25 Thread Abel Pau via gdal-dev
thing :) De: gdal-dev En nombre de Marián Hadri via gdal-dev Enviado el: dijous, 25 de gener de 2024 14:58 Para: gdal-dev@lists.osgeo.org Asunto: [gdal-dev] help Hi, I'm trying to build the GDAL dll and lib files, so GDAL can be used for linker and includes in an other C++ project. We

[gdal-dev] help

2024-01-25 Thread Marián Hadri via gdal-dev
Hi, I'm trying to build the GDAL dll and lib files, so GDAL can be used for linker and includes in an other C++ project. We're needing to upgrade from GDAL 3.2.2 to 3.8.3. I try to create project file for VS 2019. In GDAL 3.2.2 was for that purpose used "generate_vcxproj.bat", which was in root

Re: [gdal-dev] help please with close/reopen and GDALSubdatasetInfo

2023-12-03 Thread Michael Sumner via gdal-dev
Awesome, thank you! I was a bit lost - learnt a lot from exploring this though. Cheers, Mike On Mon, Dec 4, 2023 at 12:11 AM Even Rouault wrote: > Hi Michael, > > It was most missing a "break;" statement when the subdataset of interest > has been found to exit the loop. Without the break, th

Re: [gdal-dev] help please with close/reopen and GDALSubdatasetInfo

2023-12-03 Thread Even Rouault via gdal-dev
Hi Michael, It was most missing a "break;" statement when the subdataset of interest has been found to exit the loop. Without the break, the loop will continue to read papszSubdatasets[] which has been invalidated by the poSrcDS->ReleaseRef(). Ah, memory unsafe languages :-). You also had a f

[gdal-dev] help please with close/reopen and GDALSubdatasetInfo

2023-12-02 Thread Michael Sumner via gdal-dev
May I please ask for assistance with this code? I'm trying to close a dataset with subdatasets and then reopen if the requested subdataset is found. It segfaults for all but the last subdataset name ... and I'm at a loss for what I'm doing wrong. Thank you. Code below, and this gist documents i

Re: [gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-11 Thread Catania, Luke A ERDC-RDE-GRL-VA CIV via gdal-dev
he hatches not converting properly and missing symbology for the geometry. -Original Message- From: gdal-dev On Behalf Of Scott via gdal-dev Sent: Wednesday, October 11, 2023 2:07 PM To: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] help needed with pygdaltools library, documentati

Re: [gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-11 Thread Scott via gdal-dev
rdc.dren.mil>>; Scott <mailto:pub...@postholer.com>>; gdal-dev@lists.osgeo.org <mailto:gdal-dev@lists.osgeo.org> Subject: Re: [gdal-dev] help needed with pygdaltools library, documentation sparse Hi, The actual error seems to be that you try to write some geometry that is not

Re: [gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-11 Thread Catania, Luke A ERDC-RDE-GRL-VA CIV via gdal-dev
...@maanmittauslaitos.fi>> Sent: Wednesday, October 11, 2023 2:12 AM To: Catania, Luke A ERDC-RDE-GRL-VA CIV mailto:luke.a.cata...@erdc.dren.mil>>; Scott mailto:pub...@postholer.com>>; gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> Subject: Re: [gdal-dev] help needed

Re: [gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-10 Thread Rahkonen Jukka via gdal-dev
.59 Vastaanottaja: Scott ; gdal-dev@lists.osgeo.org Aihe: Re: [gdal-dev] help needed with pygdaltools library, documentation sparse I have not even gotten to see that anything is being ignored because I have not yet added any options. My first attempt was with adding skipfailures, but I can

Re: [gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-10 Thread Catania, Luke A ERDC-RDE-GRL-VA CIV via gdal-dev
pylint tells me it is not callable -Original Message- From: gdal-dev On Behalf Of Scott via gdal-dev Sent: Wednesday, October 11, 2023 1:45 AM To: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] help needed with pygdaltools library, documentation sparse -nln is being ignored because the

Re: [gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-10 Thread Scott via gdal-dev
-nln is being ignored because the data layer you're -updating -appending to exists. You can't rename a layer if it has a name. Remove -nln option Next you are trying to write a collection as a linestring. Use -explodecollections (if exists in your install) and/or use -nlt PROMOTE_TO_MULTI Tr

[gdal-dev] help needed with pygdaltools library, documentation sparse

2023-10-10 Thread Catania, Luke A ERDC-RDE-GRL-VA CIV via gdal-dev
I am running simple code to convert dxf to gpkg and I get the error below. -am an trying to figure out how to add a bunch of options ( -f GPKG -append -update --config DXF_ENCODING CP1252 --config DXF_INCLUDE_RAW_CODE_VALUES TRUE --config DXF_TRANSLATE_ESCAPE_SEQUENCES FALSE --config DXF_MERGE_

Re: [gdal-dev] help with (docker) build workflow

2023-02-19 Thread Michael Sumner
Ah, awesome - thanks Even and Daniel, both very helpful. Cheers, Mike On Mon, Feb 20, 2023 at 5:36 AM Even Rouault wrote: > Michael, > > for your use case, you are only interested in the "builder" image, not the > "runner" one. > > So just build it with something like: docker build . --target

Re: [gdal-dev] help with (docker) build workflow

2023-02-19 Thread Even Rouault
Michael, for your use case, you are only interested in the "builder" image, not the "runner" one. So just build it with something like: docker build . --target builder -t whatever_tag_name_you_want Even Le 19/02/2023 à 19:08, Michael Sumner a écrit : I would like to build GDAL in an image

Re: [gdal-dev] help with (docker) build workflow

2023-02-19 Thread Daniel Evans
Hi Michael, The idea of a two-image build like this is explicitly to create a Runner without all the compile dependencies, reducing the amount of baggage that users of the Docker image end up having. That's why it doesn't copy over all the dependencies for compilation, as you note! One of the down

[gdal-dev] help with (docker) build workflow

2023-02-19 Thread Michael Sumner
I would like to build GDAL in an image on the basis of a dockerfile used for CI, I'm confused by the layer process that copies the build outputs for the final steps https://github.com/OSGeo/gdal/blob/master/docker/ubuntu-small/Dockerfile#L231 I want to be able to run that image and then do my own

[gdal-dev] Help on DWG data extraction with GDAL/ODA

2021-12-03 Thread Maxime Colmant
Hello, I'm currently playing with GDAL and ODA to extract the data from a DWG file. I would like to extract the block attributes and I did not find any solution. I tried to use the DWG_INLINE_BLOCKS but nothing linked to these attributes are extracted. I'm not sure if it is supported right now. I

Re: [gdal-dev] Help with OGDI vector driver usage

2021-12-02 Thread Carl Godkin
Hi Even, > try: > > ogrinfo gltp:/vrf/d:/users/carl/vvod That works! I don't know why I didn't try that variation. Thank you so much! As I said, I'm using OSGeo4W to test this out. But for my real use case, I will need to build OGDI into GDAL here. The instructions for the OGDI library loo

Re: [gdal-dev] Help with OGDI vector driver usage

2021-12-02 Thread Even Rouault
Le 02/12/2021 à 05:18, Carl Godkin a écrit : Hi, I've got some data that is in VPF (Vector Product Format) format and it appears that the OGDI vector should be able to consume it. I have GDAL 3.4.0 (the OSGeo4W build) and cannot figure out how to even run "ogrinfo" on the data set. From th

[gdal-dev] Help with OGDI vector driver usage

2021-12-01 Thread Carl Godkin
Hi, I've got some data that is in VPF (Vector Product Format) format and it appears that the OGDI vector should be able to consume it. I have GDAL 3.4.0 (the OSGeo4W build) and cannot figure out how to even run "ogrinfo" on the data set. >From this page (https://gdal.org/drivers/vector/ogdi.html)

Re: [gdal-dev] help with GDALWarpApp - resample options

2021-05-14 Thread Michael Sumner
Oh, I see - you just append it all on. Thanks! C++ is a real struggle for me, but this is going to be excellent. Best, Mike On Sat, 15 May 2021, 07:18 Even Rouault, wrote: > Michael, > > > > how do I set the equivalent of eResampleAlg for the GDALWarp *app* ? > > > > > https://gdal.org/api/gda

Re: [gdal-dev] help with GDALWarpApp - resample options

2021-05-14 Thread Even Rouault
Michael, how do I set the equivalent of eResampleAlg for the GDALWarp *app* ? https://gdal.org/api/gdal_utils.html#gdal__utils_8h_1ab2094d5332364d63e0851dcb12398469 I think you can't?  GDALWarpAppOption

[gdal-dev] help with GDALWarpApp - resample options

2021-05-14 Thread Michael Sumner
Hello, how do I set the equivalent of eResampleAlg for the GDALWarp *app* ? https://gdal.org/api/gdal_utils.html#gdal__utils_8h_1ab2094d5332364d63e0851dcb12398469 I think you can't? GDALWarpAppOptions is for the -wo command line equivalents, not the -r (right?) I can see how to use the impleme

[gdal-dev] help request about WarpRegionToBuffer

2020-09-23 Thread Maud Ferrato
Dear GDAL community I'm currently working with gdal api C/C++ and I'm facing an issue with gdal warp region to buffer functionality (WarpRegionToBuffer). When my destination dataset is not strictly contained in the frame of my source dataset, the area where there should be no data values is fil

Re: [gdal-dev] Help request for data alignment issue

2020-06-17 Thread jratike80
Hi, You do not give enough details so I could try to follow your steps. It is odd that gdal_edit does not give the desired result. But have you tried to use "GEOTIFF_KEYS_FLAVOR"="ESRI_PE" directly as a gdalwarp creation option? -Jukka Rahkonen- Jessica Evans wrote > Hi > > Wondering if anyone

[gdal-dev] Help request for data alignment issue

2020-06-17 Thread Jessica Evans
Hi Wondering if anyone can help. I am currently trying to use osgeo gdal.Warp to resample and reproject a GeoTiff (from EPSG:26916 NAD83 / UTM zone 16N to EPSG:102008 North America Albers), but also to get the data to align with other data that has been processed using ArcGIS. I can get the d

[gdal-dev] Help building GDAL with ODA for the first time

2020-05-06 Thread Jorge Gustavo Rocha
Hi, I would like to compile GDAL with DWGv8 support in Ubuntu. I already have an ODA account and I've downloaded the Kernel and Drawings SDK (and an activation header): Kernel_lnxX64_5.3dll_21.3 Drawings_lnxX64_5.3dll_21.3 I'm trying to compile GDAL with: CPPFLAGS='-I/home/jgr/dev/gdal/ODA/Kerne

Re: [gdal-dev] Help with finding residual error for a given GCP.

2020-04-22 Thread Hays Barrett
That works great! Thanks Evan! From: "Even Rouault" To: gdal-dev@lists.osgeo.org Cc: "Hays Barrett" Sent: Tuesday, April 21, 2020 7:18:04 AM Subject: Re: [gdal-dev] Help with finding residual error for a given GCP. Hays, GCPsToGeoTransform() is not used

Re: [gdal-dev] Help with finding residual error for a given GCP.

2020-04-21 Thread Even Rouault
Hays, GCPsToGeoTransform() is not used directly by gdalwarp (although it should likely be equivalent to running -order 1 polynomial), but by other code (OZI .map exporting, .tab exporting) You rather want to use the gdal.Transformer() API which is tested in https://github.com/OSGeo/gdal/blob/

[gdal-dev] Help with finding residual error for a given GCP.

2020-04-20 Thread Hays Barrett
Hello all, I want to get the error for each GCP before warping so I can calculate the RMSE before georeferencing. With a list of GCPs, is it possible to get the eventual map location of each pixel coordinate to calculate the transformation error for each GCP? I have tried to use GCPsToGeoTrans

Re: [gdal-dev] Help with CPL_DEBUG

2020-02-17 Thread Mateusz Loskot
On Mon, 17 Feb 2020 at 18:34, Joaquim Manuel Freire Luís wrote: > On Mon, 17 Feb 2020 at 18:05, Joaquim Manuel Freire Luís > wrote: > >> > >> Hi Even, > > > FYI, you have e-mailed the public gdal-dev mailing list. > > Yes, I know but since Even responds all the times ... I took the liberty I ho

Re: [gdal-dev] Help with CPL_DEBUG

2020-02-17 Thread Joaquim Manuel Freire Luís
On Mon, 17 Feb 2020 at 18:05, Joaquim Manuel Freire Luís wrote: >> >> Hi Even, > FYI, you have e-mailed the public gdal-dev mailing list. Yes, I know but since Even responds all the times ... I took the liberty I also read that page, and used "CPL_DEBUG=ON" in my command, but no message was

Re: [gdal-dev] Help with CPL_DEBUG

2020-02-17 Thread Mateusz Loskot
On Mon, 17 Feb 2020 at 18:05, Joaquim Manuel Freire Luís wrote: > > Hi Even, FYI, you have e-mailed the public gdal-dev mailing list. > I have no idea why that might happen and tried to run it with CPL_DEBUG as in > gmt grdgdal lixo.csv -R0/10/0/10 -Gjunk.nc -I0.05 -F"-a nearest -co > CPL_DEBUG

[gdal-dev] Help with CPL_DEBUG

2020-02-17 Thread Joaquim Manuel Freire Luís
Hi Even, We are having a strange failure in one CI test of GMT that calls GDAL. The thing happens at this line https://github.com/GenericMappingTools/gmt/blob/master/src/gmt_gdal_librarified.c#L355 and since bUsageError comes != 0 the error message at https://github.com/GenericMappingTools/gmt/b

Re: [gdal-dev] Help with band pixel functions

2019-07-15 Thread Adriana Parra
Hello Steven, Thank you for your reply. I should have mentioned that my ultimate goal is to mosaic more than two images, so I don't know if the workflow that you suggest would work in that case. I did use your suggestion of masking values, and changed the VRT python function to this: import numpy

Re: [gdal-dev] Help user when /vsicurl/ reads GetFileSize=0

2019-06-25 Thread Even Rouault
On mardi 25 juin 2019 07:02:06 CEST Rahkonen Jukka (MML) wrote: > Hi, > > Have a look at > https://gis.stackexchange.com/questions/326773/importing-shapefile-from-web > site-if-download-link-is-compressed-file/326778 > > I tried what happens with command > > ogrinfo > /vsizip/vsicurl/https://www

[gdal-dev] Help user when /vsicurl/ reads GetFileSize=0

2019-06-25 Thread Rahkonen Jukka (MML)
Hi, Have a look at https://gis.stackexchange.com/questions/326773/importing-shapefile-from-website-if-download-link-is-compressed-file/326778 I tried what happens with command ogrinfo /vsizip/vsicurl/https://www2.census.gov/geo/tiger/GENZ2010/gz_2010_36_620_l2_500k.zip --debug on HTTP: libcu

Re: [gdal-dev] Help with band pixel functions

2019-06-24 Thread Stephen Woodbridge
Hi Adriana, Here are a few ideas: import numpy as np ima = gdal.Open(filea, 0) banda = ima.GetRasterBand(1) dataa = banda.ReadAsArray() print type(dataa), dataa.shape, isinstance(dataa, np.ma.MaskedArray) imb = gdal.Open(fileb, 0) bandb = imb.GetRasterBand(1) datab = bandb.ReadAsArray() prin

[gdal-dev] Help with band pixel functions

2019-06-24 Thread Adriana Parra
Hello, I have been trying to mosaic two images that partially overlap using band pixel functions in Python. My images contain unsigned integers (uint8) and the no data value is 255. I considered using the nanmean function so that I would get average values in overlapping areas of the two images an

[gdal-dev] Help with GPM Data

2019-06-14 Thread Rajesvari Parasa
Hello, This is the first time that I am dealing with Global Precipitation Measurement (GPM) datasets. I am downloading these from GES DISC. From the site, only NetCDF and HDF5 format files are available for download. When I open these in QGIS, I see that these are wrongly oriented. As if lat, long

Re: [gdal-dev] Help with GDAL/OGR Python API and OCI

2019-03-25 Thread Alejandro Saucedo
Yes, it was that. I replaced the library it was linked to with the one from my GDAL I built. Thank you very much! Regards. On Mon, Mar 25, 2019 at 5:56 PM Even Rouault wrote: > On lundi 25 mars 2019 17:26:17 CET Alejandro Saucedo wrote: > > Hello, > > Please, your help with the following probl

Re: [gdal-dev] Help with GDAL/OGR Python API and OCI

2019-03-25 Thread Even Rouault
On lundi 25 mars 2019 17:26:17 CET Alejandro Saucedo wrote: > Hello, > Please, your help with the following problem. > I am not able to see the Oracle OCI driver from python using the GDAL/OGR > bindings. > The following code prints always None: > from osgeo import ogr > print(ogr.GetDriver

[gdal-dev] Help with GDAL/OGR Python API and OCI

2019-03-25 Thread Alejandro Saucedo
Hello, Please, your help with the following problem. I am not able to see the Oracle OCI driver from python using the GDAL/OGR bindings. The following code prints always None: from osgeo import ogr print(ogr.GetDriverByName( 'OCI' )) Also, the driver is never listed here: drivers = [og

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-22 Thread Alex HighViz
From: Alex Highviz Sent: 22 June 2018 18:27:14 To: even.roua...@spatialys.com Subject: Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark On 14/06/2018 21:22, Even Rouault wrote: > > [ ] > - 1000x1000 is somewhat small

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-20 Thread Even Rouault
> I ended up with the following timings: > > GDAL-Only using ReadBlock(...) : 7.7 seconds > GDAL-Only using RasterIO(...): 30.1 seconds > GDAL-Only using GetLockedBlockRef(...) : 8.9 seconds > Pronto Raster using GetLockedBlockRef(...

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-20 Thread Alex HighViz
On 19/06/2018 18:41, Mateusz Loskot wrote: > On 19 June 2018 at 11:50, Alex HighViz wrote: >> On 19/06/2018 10:28, Mateusz Loskot wrote: >>> On 19 June 2018 at 11:22, Alex HighViz wrote: Regrettably I overstated the performance in my previous post, due to a bug that is now solved. With

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-19 Thread Mateusz Loskot
On 19 June 2018 at 11:50, Alex HighViz wrote: > On 19/06/2018 10:28, Mateusz Loskot wrote: >> On 19 June 2018 at 11:22, Alex HighViz wrote: >>> Regrettably I overstated the performance in my previous post, due to a bug >>> that is now solved. With the solved bug, Pronto is about 50% slower than

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-19 Thread Alex HighViz
On 19/06/2018 10:28, Mateusz Loskot wrote: > On 19 June 2018 at 11:22, Alex HighViz wrote: >> Regrettably I overstated the performance in my previous post, due to a bug >> that is now solved. With the solved bug, Pronto is about 50% slower than >> GDAL directly. I believe this is still good co

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-19 Thread Mateusz Loskot
On 19 June 2018 at 11:22, Alex HighViz wrote: > > Regrettably I overstated the performance in my previous post, due to a bug > that is now solved. With the solved bug, Pronto is about 50% slower than GDAL > directly. I believe this is still good considering the benefits, but not as > glamorous

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-19 Thread Alex HighViz
> Even Rouault wrote > A few thoughts: [ ] > 1000x1000 is somewhat small. Perhaps benchmark on larger rasters I am still working on your other suggestions, but thought it would be fair to report on my progress. Regrettably I overstated the performance in my previous post, due to a bug that is

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Alex HighViz
Great, will try and let you know. From: Mateusz Loskot Sent: 15 June 2018 00:27:34 To: Alex HighViz Cc: Even Rouault; gdal-dev@lists.osgeo.org; Hagen-Zanker AH Dr (Civil & Env. Eng.) Subject: Re: [gdal-dev] Help with idiomatic GDAL solution for raster alg

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Mateusz Loskot
On 15 June 2018 at 01:22, Alex HighViz wrote: > I need to look into sse2/ave2. My feeling is that will be easier > and more effective for the reference case than for Pronto, because pronto is > geared to express all operations at the pixel level. You may try with auto vectorization first (could b

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Alex HighViz
pronto is geared to express all operations at the pixel level. From: Even Rouault Sent: 14 June 2018 21:22:36 To: gdal-dev@lists.osgeo.org Cc: Alex HighViz; Hagen-Zanker AH Dr (Civil & Env. Eng.) Subject: Re: [gdal-dev] Help with idiomatic GDAL solution for ra

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Mateusz Loskot
On 14 June 2018 at 18:47, Alex HighViz wrote: > Would you by any chance be able (and willing) to conjure up a simple example > of "OUT = 3 * A + B * C" using Boost.GIL? Quick update: out of the box, GIL provides related algorithm that accepts one or two sources: transform_pixels Here is example

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Mateusz Loskot
On 14 June 2018 at 18:47, Alex HighViz wrote: > I remember when looking at it at the time that I felt that using Boost.GIL > would give me more difficulty than reward. It is not easy indeed. Things had not been improved for long time, esp. documenation. It's changing and improving though. > I w

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Even Rouault
Hi, Interesting work. A few thoughts: - the benchmark_3_rasters_reference() could receive a few comments for non- advanced users to mention that it is correct only if all the rasters use Byte datatype, and all have the same block size. - It would be interesting to see a variant of this referen

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Alex HighViz
6:34 To: Alex HighViz Cc: gdal-dev@lists.osgeo.org; Hagen-Zanker AH Dr (Civil & Env. Eng.) Subject: Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark On 14 June 2018 at 14:17, Alex HighViz wrote: > > I few times I have posted to the list trying to promote the idea of

Re: [gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Mateusz Loskot
On 14 June 2018 at 14:17, Alex HighViz wrote: > > I few times I have posted to the list trying to promote the idea of > providing iterators over pixels in a raster band , and more generally to > make raster data accessible using (future) standard conforming ranges. In case you are not aware, Boos

[gdal-dev] Help with idiomatic GDAL solution for raster algebra benchmark

2018-06-14 Thread Alex HighViz
Dear all, I few times I have posted to the list trying to promote the idea of providing iterators over pixels in a raster band , and more generally to make raster data accessible using (future) standard conforming ranges. It would make implementing algorithms on raster data a lot more intuitive

Re: [gdal-dev] help printing the standard out from gdal_grid in python

2018-04-12 Thread Sean Gillies
Hi Nicolas, check_output() is higher level and easier to use. I recommend it: https://docs.python.org/2.7/library/subprocess.html#subprocess.check_output. There's a note in there about capturing stderr as well. If gdal_grid is running successfully, there may not be anything in stdout. I think I s

[gdal-dev] help printing the standard out from gdal_grid in python

2018-04-12 Thread Nicolas Cadieux
Hi, I am not an experienced programmer.  I am calling gdal_grid from a python script.  I am using p.imap_ unordered to multiprocess the script.  The idea is the following: def gdal_Grid(fileList):     gdal_Grid_Cmd = 'gdal_grid -l '     shell = subprocess.Popen(r'c:\

Re: [gdal-dev] [help][gdalwarp]Use of gdalwarp with data overlapping antimeridian

2018-03-13 Thread LAVOYE Christophe
18 15:41 À : LAVOYE Christophe Cc : gdal-dev@lists.osgeo.org Objet : Re: [gdal-dev] [help][gdalwarp]Use of gdalwarp with data overlapping antimeridian Hi Christophe, Have a look at the following: https://lists.osgeo.org/pipermail/gdal-dev/2016-January/043358.html https://gis.stackexchange.com

Re: [gdal-dev] [help][gdalwarp]Use of gdalwarp with data overlapping antimeridian

2018-03-08 Thread Even Rouault
On jeudi 8 mars 2018 14:40:44 CET Robert Coup wrote: > Hi Christophe, > > Have a look at the following: > > https://lists.osgeo.org/pipermail/gdal-dev/2016-January/043358.html > https://gis.stackexchange.com/questions/37790/how-to-> > reproject-raster-from-0-360-to-180-180-with-cutting-180-merid

Re: [gdal-dev] [help][gdalwarp]Use of gdalwarp with data overlapping antimeridian

2018-03-08 Thread Robert Coup
I’m new on this mail list do I have to “Reply All”? or just reply to > you? > > > > Best regards, > > > > > > [image: cid:image003.png@01D1FD34.3B7C8CC0] > > *Christophe LAVOYE - **Analyste* > * AKKA Technologies *Informatique & Systèmes > 7 boulevard Henri Ziegle

Re: [gdal-dev] [help][gdalwarp]Use of gdalwarp with data overlapping antimeridian

2018-03-07 Thread Robert Coup
Hi Christophe, What do you *expect or want* to happen with your images? Are your 4326 image coordinates in -180° to +180° range (ie. the right longitude is < the left); or do they go past +180°? 3857 is a projected coordinate system, so doesn't have a builtin concept of wrapping. Some clients im

[gdal-dev] [help][gdalwarp]Use of gdalwarp with data overlapping antimeridian

2018-03-01 Thread LAVOYE Christophe
Hello, I'm using GDAL library in a C# project. To reproject some geotiff files i'm calling gdalwarp tool as external process with this command: gdalwarp.exe -s_srs EPSG:4326 -t_srs EPSG:3857 -dstalpha input.tiff output.tiff All works fine except for data overlapping antimeridian such as the

[gdal-dev] Help: No coordinate conversion using OGRSpatialReference/OGRCoordinateTransformation in version 2.2.3

2018-01-15 Thread Davide De Marchi
Dear GDAL list, I'm having problems using the standard way to convert coordinate points from a source to a target SRS. After updating to version 2.2.3, code that always worked started to leave the input coordinate unchanged, without giving any error or warning.This is the code I use: bool Co

Re: [gdal-dev] help

2018-01-13 Thread Rasika Chinchwade
I have used STD dataset of INSAT-3D https://www.mosdac.gov.in/catalog?field_satellitenode_nid=69 you can download the datset from www.mosdac.gov.in HEM daily ,LST product. product of INSAT3D .This all are in disc shape. I am trying to convert the .h5 dataset to GeoTiff format using gdal. This produ

Re: [gdal-dev] help

2018-01-13 Thread Even Rouault
On samedi 13 janvier 2018 14:02:50 CET Rasika Chinchwade wrote: > hello all, > Why I am getting this error while performing same operation on same dataset > given on link Providing an actual link to a dataset could be helpful -- Spatialys - Geospatial professional services http://www.spatialys

[gdal-dev] help

2018-01-13 Thread Rasika Chinchwade
hello all, Why I am getting this error while performing same operation on same dataset given on link https://gis.stackexchange.com/questions/253884/projecting-insat3d-l1b-data-from-hdf5-to-geotiff ERROR 1: tolerance condition error ERROR 1: tolerance condition error ERROR 1: tolerance condition e

[gdal-dev] Help with getting gdal python module working

2017-05-15 Thread Candy, Liam
Hi, I've managed to get GDAL 2.1.3 to build with ECW 5.3 support on my 64-bit system. The last hurdle I'm struggling with is getting the python module to install and work. The lines in my nmake.opt file relating to python are; !IFNDEF PYDIR PYDIR = "... \AppData\Local\Programs\Python\Python36"

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-28 Thread Grégory Bataille
Hey, just to let you know, I found the issue. So I think it has to do with the fact that I keep some handles opened on the original TIFF In the parent thread. Even setting the variable to None does not change anything, but I guess it's because the garbage collector is asynchronous. So if I do my i

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-21 Thread Grégory Bataille
hum... debugging in the C layers, that'll be interesting... Ok, I'll continue to dig in, thanks --- Gregory Bataille On Tue, Mar 21, 2017 at 8:44 AM, Even Rouault wrote: > On mardi 21 mars 2017 05:51:18 CET Grégory Bataille wrote: > > > Hey Even, > > > > > > reaching out again, because I just

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-21 Thread Even Rouault
On mardi 21 mars 2017 05:51:18 CET Grégory Bataille wrote: > Hey Even, > > reaching out again, because I just don't know... > > so here is gdal2tiles_parallel strategy (as far as I understand it) > - in each process, open the input file, compute the autowarped vrt (in the > case of the dataset I

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-20 Thread Grégory Bataille
Hey Even, reaching out again, because I just don't know... so here is gdal2tiles_parallel strategy (as far as I understand it) - in each process, open the input file, compute the autowarped vrt (in the case of the dataset I use at least) loop over tiles to generate and work only depending on the

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-16 Thread Grégory Bataille
no luck. I tried this config. I also tried (with this config) to do a copy of the vrt file (with shell utility, not gdal) in each thread before opening it (still goes to the same TIFF file though), but no luck. But at least now that you confirm that it's related to concurrent reads, I have somethi

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-16 Thread Even Rouault
On jeudi 16 mars 2017 18:45:03 CET Grégory Bataille wrote: > Does that mean a different vrt file or simply reopen the file to create a > gdal object in each process. Because I'm reopening the vrt file in each > thread Ah, that must be issue of http://gdal.org/gdal_vrttut.html#gdal_vrttut_mt Try d

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-16 Thread Grégory Bataille
Does that mean a different vrt file or simply reopen the file to create a gdal object in each process. Because I'm reopening the vrt file in each thread On Thu, 16 Mar 2017 at 17:23, Even Rouault wrote: > On jeudi 16 mars 2017 17:16:20 CET Grégory Bataille wrote: > > > Hello all, > > > > > > Re

Re: [gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-16 Thread Even Rouault
On jeudi 16 mars 2017 17:16:20 CET Grégory Bataille wrote: > Hello all, > > Reaching out to the community caused I have failed for the past few days. > > *short version*, I'm trying to multithread the gdal2tiles utility, and I'm > getting this > > Generating Base Tiles: > ERROR 1: LZWDecode:Wron

[gdal-dev] Help requested: Concurrent read of a GeoTiff

2017-03-16 Thread Grégory Bataille
Hello all, Reaching out to the community caused I have failed for the past few days. *short version*, I'm trying to multithread the gdal2tiles utility, and I'm getting this Generating Base Tiles: ERROR 1: LZWDecode:Wrong length of decoded string: data probably corrupted at scanline 256 ERROR 1:

Re: [gdal-dev] Help working with hdf5 images

2017-02-19 Thread Stephen Woodbridge
On 2/19/2017 10:15 AM, Stephen Woodbridge wrote: Hi all, I'm trying to create GTiff files from IMARS Modis hdf5 images. I did something like this about 3-4 years ago when they were using hdf4 datasets but things have changed since then. So I have two questions: 1. How do I extract the sst data

[gdal-dev] Help working with hdf5 images

2017-02-19 Thread Stephen Woodbridge
Hi all, I'm trying to create GTiff files from IMARS Modis hdf5 images. I did something like this about 3-4 years ago when they were using hdf4 datasets but things have changed since then. So I have two questions: 1. How do I extract the sst data (subdataset 3) into a GTiff and create a mask

[gdal-dev] help

2016-10-17 Thread jp internet
___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] Help regarding DXF format interpretation

2016-08-04 Thread Avyav Kumar
Hi all, I am a student who is participating under the GSoC program, and I need the help of the community! https://trac.osgeo.org/gdal/ticket/6246 contains an issue that I aim to resolve. However, reading the DXF format is a difficulty (for me). I don't understand how the indexing corresponding to

Re: [gdal-dev] help about the GDAL efficiency

2015-05-20 Thread Even Rouault
Xia Lang, > > > // create out file > eDT = GDALDataType(iBand1->GetRasterDataType()); > GDALDriver *hDriver = (GDALDriver *)GDALGetDriverByName( "GTiff" ); > CPLAssert( hDriver != NULL ); > hDstDS = (GDALDataset *)GDALCreate( hDriver, "G:\\data\\TM\\out.tif", > nPixels, nLines, GDALGetRaste

[gdal-dev] help about the GDAL efficiency

2015-05-19 Thread xialang2012
Hi, I used GDAL to reproject the TM images, and the original projection of the TM images was longitude/latitude WGS84, the final projection was UTM 51. The codes I used is as following: void doCore() { GDALDataType eDT; GDALDataset *hDstDS; GDALDataset *hSrcDS; GDALAllRegister(); // c

[gdal-dev] Help formatting 2 million images for geoserver WMS via GDAL

2014-05-21 Thread Peter Lenou
I am trying to use the GDAL command line utilities to format 2 million images (geotiff) of variable size in a form that can be uploaded to geoserver to be served via WMS. I am having problems determining what is the best method to format the data i.e. do I pyramid the whole set using gdal_retile.p

Re: [gdal-dev] Help request to make Python bindings to build again with MacOS X 10.9

2014-05-03 Thread Even Rouault
> Maybe there is something weird in the Travis system. They must have changed something to fix it, since builds now pass. > > - > William Kyngesburye > http://www.kyngchaos.com/ > > "Time is an illusion - lunchtime doubly so." > > - Ford Prefect -- Geospatial professional services htt

Re: [gdal-dev] Help request to make Python bindings to build again with MacOS X 10.9

2014-04-20 Thread William Kyngesburye
On Apr 20, 2014, at 4:39 AM, Even Rouault wrote: >> Well, you could try leaving out the os.uname test. clang is clang no matter >> what version of OS X, and it wouldn't hurt to include the option when it's >> already default. Also, testing the OS version doesn't necessarily say what >> clang v

Re: [gdal-dev] Help request to make Python bindings to build again with MacOS X 10.9

2014-04-20 Thread Kurt Schwehr
http://en.wikipedia.org/wiki/Xcode#Composition says that clang started at xcode 4.2, for which the version numbering is decoupled from the darwin version. There are two tables at the bottom that only sort of tell the whole story... http://en.wikipedia.org/wiki/Xcode#Versions_compare_table e.g. o

Re: [gdal-dev] Help request to make Python bindings to build again with MacOS X 10.9

2014-04-20 Thread Even Rouault
> Well, you could try leaving out the os.uname test. clang is clang no matter > what version of OS X, and it wouldn't hurt to include the option when it's > already default. Also, testing the OS version doesn't necessarily say what > clang version is installed. Could you try leaving out test os.

Re: [gdal-dev] Help request to make Python bindings to build again with MacOS X 10.9

2014-04-19 Thread William Kyngesburye
On Apr 19, 2014, at 7:44 PM, Even Rouault wrote: > Selon William Kyngesburye : > >> On Apr 19, 2014, at 2:56 PM, Even Rouault >> wrote: >> >>> Le samedi 19 avril 2014 21:19:53, William Kyngesburye a écrit : It looks like a bit of a coordination and packaging snafu with Apple. Apple b

  1   2   3   >