Do you mean a MultiLineString?
This piece of code is working for me (I hope without any bug).
if (wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString)
multiLineStringGeometry(poGeometry);
...
multiLineStringGeometry(OGRGeometry* poGeometry)
{
Thanks Even for notifying this.
I see the fix will be included in the next releases 3.7.3 and 3.8.0, both
planned for November 1st (just in a couple of weeks)
https://github.com/OSGeo/gdal/milestones
On Mon, 16 Oct 2023 at 16:42, Even Rouault via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
+1 Javier
On Wed, 11 Oct 2023 at 21:09, Howard Butler via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Yes. Well "up to", of course, but Alessandro is providing us some more
> capacity.
>
> > On Oct 11, 2023, at 1:56 PM, Sean Gillies
> wrote:
> >
> > Hi Howard,
> >
> > To be clear, this is a
Hi Even. Thanks, it sounds good.
However I see a potential problem. I see that you use once "SetCacheMax".
We should not forget about that in the future for sensible tests. The cache
of gdal is usually a percentage of the total memory, that may change among
the environments and time.
On Wed, 11
+1 Javier
On Wed, 1 Nov 2023 at 12:19, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: keskiviikko 1. marraskuuta 2023 13.14
> Vastaanottaja:
+1 Javier
On Tue, 31 Oct 2023, 18:04 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> +1 Even
>
> Le 31/10/2023 à 17:59, Howard Butler via gdal-dev a écrit :
> > Dear PSC,
> >
> > Sorry for the short notice, but I would like to motion that the GDAL
> Sponsorship Program
Wow Even. There's a lot of new things there! Thanks to all the contributors!
On Fri, 3 Nov 2023, 18:53 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
>
> > I've prepared a beta1 of GDAL 3.8.0 to get feedback from earlier testers.
> >
> > Sorry no updated NEWS.md file yet
> ==> now
I do not know that code. But having a look I have the impression that you
have to set the nodata also to papszCreateOptions, so the created dataset
has it.
On Fri, 3 Nov 2023 at 21:36, Fitch, Simeon via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> I'm one of the maintainers at the
Hi
Using a GeoTIFF with a single band and float values can be very useful to
show any distribution over the terrain. One typical example is DSM (digital
surface model), where the value is the elevation on every point.
Is there any standard (or accepted) metadata to say the units of those
values?
Thanks!
On Thu, 26 Oct 2023 at 12:27, Even Rouault via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Javier,
> Le 26/10/2023 à 11:59, Javier Jimenez Shaw via gdal-dev a écrit :
>
> Hi
>
> Using a GeoTIFF with a single band and float values can be very useful to
>
This is working for me (and also in gdal 3.8.0):
gdalwarp 3635_rasters_agreges.jp2 salida.tif -cutline
geometry_extraction.shp -crop_to_cutline -dstnodata 0 -cl
geometry_extraction -overwrite -of GTiff
About the "white" pixels inside the image, it could be that a single band
has a value of 0
Hi Naima
I have been testing with your dataset. To me, using the GDAL in Ubuntu
22.04 (3.4.1) seems to be a problem with the JP2 output format. If you
output as geotiff it works fine.
On Sun, 5 Nov 2023 at 19:43, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
>
>
>
This is producing something nicer, without those "white pixels". It is
adding an alpha band:
gdalwarp 3635_rasters_agreges.jp2 salida.jp2 -cutline
geometry_extraction.shp -crop_to_cutline -cl geometry_extraction -overwrite
-dstalpha
On Mon, 6 Nov 2023 at 14:05, Javier Jimenez Shaw wrote:
> This
I never used that env bastante before, but it seems to be what you
describe. See that it is enabled also in the tests, as Even says in the
pull 2666
I'm curious. Why do you need to do something different on GDAL compilation?
(I assume it is in a exported header)
On Thu, 21 Sept 2023, 15:56 Abel
Are you running pycharm under the conda environment? If you lose any
setting while combining both it can be a problem.
Have you tried setting PROJ_DEBUG=3 ?
It will show some traces from PROJ, like the location of proj.db used. The
point is if you can see the std out.
On Wed, 11 Oct 2023, 15:18
GDAL got 25 years old on 2023-10-17 !
I found it in this post
https://geoobserver.wordpress.com/2023/10/17/25-jahre-gdal-happy-birthday/
.___ ._ ..._ .. . ._. .___ .. __ . _. . __.. ... ._ .__
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
+1 Javier
On Thu, 9 Nov 2023, 11:19 Rahkonen Jukka via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: torstai 9. marraskuuta 2023 12.12
> Vastaanottaja:
Assuming that it does not happen in the main area of Ireland...I guess it
is because it is outside of the area of usage of the datum transformation
https://epsg.org/transformation_1641/TM65-to-WGS-84-2.html
Probably (this is just a guess) gdalwarp is not taking it into
consideration, making a
If I remember correctly, DXF does not support multipolygons with different
elevations.
On Wed, 24 Jan 2024 at 16:15, Abel Pau via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> HI,
>
> I’ve checked that your file is correct but for some reason DXF does not
> get this capacity of variable
Hi Mike
Out of curiosity, are run running it in a virtual machine?
A few year ago I had problems running a program in a virtual machine
(virtualbox, but I read it happens in others) due to a missing SSE
instruction. The solution there was to "enable" the missing instructions in
the virtual
Could you set up your VMs to include those SSE instructions? I think that
keeping VMs that "old" configured is a source of problems using
pre-compiled binaries.
The same way GDAL updates dependencies of compilers and other libraries to
something more modern (but not too modern), those SSE
+1 Javier
On Fri, 9 Feb 2024 at 10:13, Even Rouault via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> Motion:
>
> Adopt GDAL 3.8.4RC1 as 3.8.4 release
>
> Starting with my +1
>
> Even
>
> --
> http://www.spatialys.com
> My software is free, but my time generally not.
>
>
The origin of this question is this ticket in QGIS:
https://github.com/qgis/QGIS/issues/56288
The summary is that with rasters ordered bottom-up, QGIS is doing a Warped
VRT. But Warp us using square pixels... that is distorting too much the
image in some cases (and I would say it is more
Hi
I am calling the python method osgeo.gdal.VectorTranslate to convert a
GeoJSON (with CRS) to shapefile.
To keep consistency among EPSG versions, I have the pipeline of the
conversion, using what in pyproj was "always_xy=True"
gdal.VectorTranslate(
dst_file,
Thank you.
Yes, I got confused with gdal.VectorTranslate and gdal_translate, when I
really meant ogr2ogr. But you realized that already. Great catch.
On Wed, 7 Feb 2024 at 18:09, Even Rouault
wrote:
> Javier,
>
>
> Is there a way to tell gdal.VectorTranslate that transformation is
>
+1 Javier
On Sat, 16 Dec 2023, 15:01 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> I have prepared an advanced-of-time GDAL/OGR 3.8.2 release candidate,
> mostly
> to fix a 3.8.0 regression that has been reported lately
> (https://github.com/OSGeo/gdal/issues/8967),
ent, but no CRS keys.
>
> Even
> Le 02/01/2024 à 13:14, Javier Jimenez Shaw via gdal-dev a écrit :
>
> Hi
>
> Does it make any sense to provide COG via HTTP POST method?
> If I understand correctly, the "magic" of downloading only what is needed
> from the server
+1 Javier
On Tue, 2 Jan 2024, 12:56 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> I have prepared an advanced-of-time GDAL/OGR 3.8.3 release candidate,
> mostly
> to fix a 3.8.0 regression that has been reported lately
> (https://github.com/OSGeo/gdal/issues/8998), that
Hi
Does it make any sense to provide COG via HTTP POST method?
If I understand correctly, the "magic" of downloading only what is needed
from the server is done with HTTP Range requests. But I am not sure if that
works in POST method (it does work in GET method, when enabled ;).
In case it could
+1 Javier
On Thu, 4 Jan 2024, 19:32 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> A third (and hopefully final) release candidate for 3.8.3 with the
> following additional fixes on top of RC2:
> - re-enable multi-threaded ArrowArray interface in GPKG driver now that
>
+1 Javier
On Wed, 3 Jan 2024, 21:46 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> I've issued a 2nd release candidate for 3.8.3 with the following
> additional fixes:
> - ogr2ogr: do not use ArrowArray interface if -clipsrc, -clipdst, -gcp
> or wrapdateline are
This is the discussion I have with my colleagues every time (but with
EPSG). The wanted to use an integer, and I ask them again and again to use
a better notation like "EPSG:{code}". The method that Even mentions,
SetFromUserInput, works with "EPSG:code", with "ESRI:code", with a WKT,
with a
+1 Javier
On Tue, 28 Nov 2023 at 17:29, Howard Butler via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1 Howard
>
> > On Nov 28, 2023, at 9:57 AM, Rahkonen Jukka via gdal-dev <
> gdal-dev@lists.osgeo.org> wrote:
> >
> > +1
> >
> > -Jukka-
> >
> > -Alkuperäinen viesti-
> > Lähettäjä:
> * use TILED=YES by default in GTiff driver
>
> Was that decided? As some people pointed (including me) a big tiled
geotiff without overviews take ages to load in QGIS.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
+1 Javier
On Mon, 27 Nov 2023 at 17:30, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: perjantai 24. marraskuuta 2023 13.06
> Vastaanottaja:
+1 Javier
On Fri, 1 Dec 2023 at 16:32, Even Rouault via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> Motion: adopt RFC 98: Build requirements for GDAL 3.9
>
> https://github.com/OSGeo/gdal/pull/8802
>
> Starting with my +1,
>
> Even
>
> --
> http://www.spatialys.com
> My software is
On Thu, 23 Nov 2023 at 15:01, Even Rouault
wrote:
>
> Le 23/11/2023 à 14:53, Javier Jimenez Shaw via gdal-dev a écrit :
>
> Hi
>
> A colleague of mine is compiling in Mac. I am aware of the problems with
> arrow and libkml:
> https://gdal.org/development/building_f
My colleague (I have only Linux) will give more info.
On Thu, 23 Nov 2023 at 16:25, Even Rouault
wrote:
>
>
>> I thought that it was enough disabling GDAL_USE_EXTERNAL_LIBS (we are not
>> enabling explicitly arrow or kml). Am I wrong?
>>
>> He might have to remove CMakeCache.txt because it
Hi
A colleague of mine is compiling in Mac. I am aware of the problems with
arrow and libkml:
https://gdal.org/development/building_from_source.html#building-on-macos
But we are compiling with GDAL_USE_EXTERNAL_LIBS=OFF
However he has those problems with arrow and libkml, and had to unistall
+1 Javier
On Wed, 22 Nov 2023 at 18:18, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: keskiviikko 22. marraskuuta 2023 14.06
> Vastaanottaja:
and call SetMetadataItem() on it
>
> Even
> Le 04/12/2023 à 16:14, Javier Jimenez Shaw via gdal-dev a écrit :
>
> Hi
>
> I am trying to convert a "normal" tiff file into a COG. For that I am
> using an intermediate vrt file to add the geolocation param
Hi
I am trying to convert a "normal" tiff file into a COG. For that I am using
an intermediate vrt file to add the geolocation parameters, and finally
copy the vrt into a COG.
I am doing it in C++, but it should be similar doing it on the command line
I hope. (In some cases I am using also other
What about NaN?
On Thu, 25 Jan 2024, 18:07 Abel Pau via gdal-dev,
wrote:
> Hi,
>
> there is any value in GDAL for VECTORS that indicates that a concrete
> value of a Z is not known (z nodata value)?
>
> I couldn’t find it anywhere.
>
>
>
> In MiraMon format we use one concrete number documented
+1 Javier
On Wed, 15 Nov 2023 at 21:44, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: keskiviikko 15. marraskuuta 2023 11.52
> Vastaanottaja:
+1 Javier
On Wed, 8 Nov 2023, 10:16 Rahkonen Jukka via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: keskiviikko 8. marraskuuta 2023 10.39
> Vastaanottaja:
Hi Abel,
If you have clear when it started happening, maybe you can try to find the
reason using https://git-scm.com/docs/git-bisect
specially as you are reproducing it in your computer (CI is working fine,
so that should be something related to your configuration).
Replicate it working fine with
o, Windows debug. That was not easy (when I worked in Windows, more
than 10 years ago). It is using different names for libraries in debug and
release.
'ogr_VDV.dir\Debug\ogr_VDV.lib': could you check that is the content of
that directory? I wouldn't be surprised if the library file there
I am not an expert in Windows at all. But the message says that cannot find
"TransformWithErrorCodes...blablah" in ogr2gr.exe. Could it be that this
method is actually in a dll, and the dll is not found? can be a PATH
problem? There are tools to see the exported methods from a dll. Check that
it
I hate that.
I never use "long" exactly for that problem. When the size is important,
uint64_t and friends are very useful, and well defined.
On Fri, 8 Mar 2024 at 10:16, Abel Pau via gdal-dev
wrote:
> FINALLY!
>
> I discovered the problem.
>
>
>
> I was using “unsigned long” as 4 bytes
Hi Even.
In principle the idea sounds good.
How is it parsing the numbers? is it locale agnostic? I think it is not,
because it is using "strtod". That means that if I have my locale in
Spanish, French, German, ... it will expect "," as the decimal separator,
right?
... well, how is GDAL
I don't like that the behaviour of a command line depends on the
configuration of the user (that is usually not aware of). So a command that
works for me doesn't work for you. That is bad.
I don't know in GDAL, but in proj there is the option --bbox, that has
comma separated coordinates. That can
Probably here:
https://github.com/OSGeo/gdal/blob/master/frmts/gtiff/libtiff/tiffvers.h
On Mon, 18 Mar 2024 at 11:28, Luís Moreira de Sousa via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi all,
>
> I am currently facing the error "PREDICTOR=2 is only supported with 64
> bit samples starting
In my company we confirmed that "Windows heap allocation mechanism sucks."
Closing the application after using gtiff driver can take many seconds due
to memory deallocations.
One workaround was to use tcmalloc. I will ask my colleagues more details
next week.
On Thu, 21 Mar 2024, 01:55 Even
+1 Javier
On Tue, 2 Apr 2024 at 14:49, Howard Butler via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> Howard
>
> > On Apr 2, 2024, at 6:09 AM, Rahkonen Jukka via gdal-dev <
> gdal-dev@lists.osgeo.org> wrote:
> >
> > +1
> >
> > -Jukka Rahkonen-
> >
> > -Alkuperäinen viesti-
> >
What I do? When opening the layer I create some memory blocs of
> 250 Mb (due to the format itself) and I use that created memory to manage
> whatever I need. And when closing, I free it.
>
>
>
> While doing that I observed that sometimes I have to use GDAL code that
> doesn’t take it in
://spatialreference.org/ref/epsg/900913/proj4.txt
On Sat, 13 Apr 2024 at 11:12, Sebastiaan Couwenberg via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> On 4/13/24 10:19 AM, Javier Jimenez Shaw via gdal-dev wrote:
> > Bas, are they really equivalent?
> As far as I know they are, whe
If what you need is really EPSG:3857, yes, use it.
However I have seen strange parameters on your projection. The radius of
the sphere is the "average" 3671 km, and you set a false easting and
northing of just 4.4 km. Is that trying to correct the radius of the
sphere? I do not know why you need
aphic 2D EPSG 6422 ESRI 106047
0
See that EPSG ones are deprecated. (surprisingly the ellipsoid EPSG:7048 is
not deprecated, but the datum that uses it is deprecated).
On 4/13/2024 4:19 AM, Javier Jimenez Shaw via gdal-dev wrote:
>
> If what you need is really EPSG:3857, yes, use it.
>
&g
Hi Nathan
Unfortunately there is not ColorInterp for NIR or RedEdge. But you can use
the "Description" of the bands. QGIS will use that name. I do not know how
to do it in gdal_merge (I did in C++).
About computing the NDVI, as far as I know there is no trick to do it in
the "Layer Styling". If
In addition to all the previous, what about gamma correction?
Maybe your RGB images are gamma corrected, while the other bands are not.
In that case you will be comparing apples and oranges (in addition to the
scaling problems described before).
That also means that using the multispectral images
Hi
I am accessing a big file in S3. All works fine (when I learned how to set
the authentication environment variables), but gdalinfo is checking for a
lot of files that do not exist. I know there is only a GeoTIFF - COG (no
side car files)
However gdalinfo is trying to find several side car
https://gdal.org/user/configoptions.html). Typically, I set
> GDAL_DISABLE_READDIR_ON_OPEN=EMPTY_DIR which will disable this.
>
>
>
> Mike
>
>
>
>
>
> --
>
> Michael Smith
>
> Remote Sensing/GIS Center
>
> US Army Corps of Engineers
>
&g
Not knowing any detail about TileDB, I find good this bump in GDAL 3.9.0,
that is already changing many other dependency minimum versions.
On Wed, 24 Apr 2024 at 20:45, Even Rouault via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
> A future TileDB version will remove various deprecated
Javier Jimenez Shaw via gdal-dev, <
> gdal-dev@lists.osgeo.org> wrote:
>
>> Hi
>>
>> I would like to read in QGIS or GDAL an interpolated value in a DSM
>> (well, actually it is a geoid model, but it is the same behaviour). See
>> that I do not
Hi
I would like to read in QGIS or GDAL an interpolated value in a DSM (well,
actually it is a geoid model, but it is the same behaviour). See that I do
not want the pixel value, but the linear interpolation among the neighbour
pixels, assuming that the pixel value is in the center of the pixel.
I know that this is not what you want... but have you considered to use
COG? Openlayers supports it, Leaflet I am not sure (I hope it does, at
least via plugin).
To improve performance, I warped to web mercator with GDAL, so Openlayers
is not wasting time on the reprojection.
On Tue, 30 Apr 2024
Maybe I don't understand your question, but in gdalinfo you have the four
corners as lat-lon
On Fri, 3 May 2024, 21:46 Simon Eves via gdal-dev,
wrote:
> So we are trying to optimize our raster import process, and particularly
> the steps to derive the final WGS84/4326 bounding box for a raster
Now I think I understand what you mean.
One effective way is to convert the four sides of the image. To avoid the
conversion of every pixel, you can start a partition process. Then compare
the difference between the transformed central point and just the mid point
between the extremes. If that
Your file prj.adf is
ProjectionLAMBERT_AZIMUTHAL
Units METERS
ZunitsNO
Xshift0.0
Yshift0.0
Parameters6378137.0 6378137.0
6378137.0 /* radius of the sphere of reference
20 0 0.0 /* longitude of center of projection
5 0 0.0 /* latitude of center of
+1 Javier
On Wed, 8 May 2024 at 13:15, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> +1
>
> -Jukka Rahkonen-
>
> -Alkuperäinen viesti-
> Lähettäjä: gdal-dev Puolesta Even
> Rouault via gdal-dev
> Lähetetty: keskiviikko 8. toukokuuta 2024 14.13
> Vastaanottaja:
asterBandH, double dfPixel, double
>> dfLocation, GDALRIOResampleAlg eInterpolation, double *pdfValue), that
>> could be used by a new mode of gdallocationinfo. The GDALRPCGetDEMHeight()
>> function in alg/gdal_rpc.cpp is a plausible candidate implementation for
>>
Hi
The github actions on a branch I created two days ago is failing in Conda
environment, and I do not understand the reason
https://github.com/jjimenezshaw/gdal/actions/runs/9161579309/job/25186735095#step:8:244
I think it is not related to my code.
Does anybody know what is it about?
Thanks!
On Mon, 20 May 2024 at 20:17, Even Rouault
wrote:
>
> As suggested in the comment of #939, attempting
> https://github.com/rouault/gdal/commit/1000cfaa950a18bd14b0c72109439cad58c789d5
>
> That fixed it. Merged into master
>
> -- http://www.spatialys.com
> My software is free, but my
Hi
I am using vscode to edit (and build and debug) gdal. So far so good.
The compilation works with a simple task that runs a shell (see code below
if you need it). Debug also works fine.
My problem is syntax highlighting. Apparently it does not realize that
"GDAL_COMPILATION" is defined, so
On Wed, 29 May 2024 at 08:59, Rahkonen Jukka via gdal-dev <
gdal-dev@lists.osgeo.org> wrote:
> Hi,
>
>
>
> When you have a RasterBand from the WMTS data source, it is abstracted
> and you can read the raster data just like from any other data source and
> raster band
>
75 matches
Mail list logo
