Re: [gdal-dev] Transformations with same horizontal datum but different geoid outputs the same Z value

2022-10-14 Thread Andrew C Aitchison
On Fri, 14 Oct 2022, Diogo wrote: Hello everyone! I'm merging multiple elevation datasets from different countries and I would like to re-project all into the same coordinates system, 4326+3855 using gdalwarp. I'm using `gdaltransform` to confirm that there is in fact a conversion from the orig

[gdal-dev] Transformations with same horizontal datum but different geoid outputs the same Z value

2022-10-14 Thread Diogo
Hello everyone! I'm merging multiple elevation datasets from different countries and I would like to re-project all into the same coordinates system, 4326+3855 using gdalwarp. I'm using `gdaltransform` to confirm that there is in fact a conversion from the original projection to the aspected one.

Re: [gdal-dev] Create tile index by time?

2022-10-14 Thread Stefan Gofferje
Hi, On 10/14/22 12:52, Thomas Knudsen wrote: Sounds like a STAC related task to me - https://stacspec.org/en after having a quick look, I have to say that I didn't understand anything about that STAC thingy. Not sure how that is supposed to help me. I need to creat

Re: [gdal-dev] Create tile index by time?

2022-10-14 Thread Thomas Knudsen
Sounds like a STAC related task to me - https://stacspec.org/en Den fre. 14. okt. 2022 kl. 09.49 skrev Stefan Gofferje < li...@home.gofferje.net>: > Morning all! > > I have a little hobby project where I pull Sentinel imagery of my town's > area and create various composites. I run that every mor

Re: [gdal-dev] C++ Getting OGRGeometry out of GDALDataset

2022-10-14 Thread Andrew C Aitchison
On Fri, 14 Oct 2022, Johannes Paul wrote: Hello, Yes i would need a polygon with the raster footprint. apps/gdaltindex.cpp implements the gdaltindex application and is written in C++, so should do what you wish. Note that this assumes the raster geolocation is given by the affine transformat

Re: [gdal-dev] ogr2ogr preserve fid behavior

2022-10-14 Thread Alexandre Gacon
Hi I managed to get the expected result by removing the options -preserve_fid and --config PGSQL_OGR_FID id and by using an existing table set correctly with the right column as primary key. Alexandre Le jeu. 13 oct. 2022 à 10:28, Rahkonen Jukka < jukka.rahko...@maanmittauslaitos.fi> a écrit :

Re: [gdal-dev] C++ Getting OGRGeometry out of GDALDataset

2022-10-14 Thread Laurențiu Nicola via gdal-dev
In that case you need to create it manually. There's an example here, using Python, but the API is pretty similar: https://pcjericks.github.io/py-gdalogr-cookbook/geometry.html#create-a-polygon. Remember to close it by duplicating the first vertex. To get the footprint itself, you'll need to lo

Re: [gdal-dev] C++ Getting OGRGeometry out of GDALDataset

2022-10-14 Thread Laurențiu Nicola via gdal-dev
PS: Sorry Johannes for getting your name wrong, my email client has a fun little bug. On Fri, Oct 14, 2022, at 10:55, Laurențiu Nicola wrote: > Hi Mats, > > Vector and raster datasets are pretty much different things in GDAL. And > generally it would be a good idea to check whether e.g. GetLaye

Re: [gdal-dev] C++ Getting OGRGeometry out of GDALDataset

2022-10-14 Thread Johannes Paul
Hello, Yes i would need a polygon with the raster footprint. Thanks Le ven. 14 oct. 2022 à 10:01, Laurențiu Nicola via gdal-dev < gdal-dev@lists.osgeo.org> a écrit : > Hi Mats, > > Vector and raster datasets are pretty much different things in GDAL. And > generally it would be a good idea to che

Re: [gdal-dev] C++ Getting OGRGeometry out of GDALDataset

2022-10-14 Thread Laurențiu Nicola via gdal-dev
Hi Mats, Vector and raster datasets are pretty much different things in GDAL. And generally it would be a good idea to check whether e.g. GetLayer actually returns a pointer or nullptr. What exactly do you mean by "associated OGRPolygon"? Do you want a polygon with the same footprint as the ra

[gdal-dev] C++ Getting OGRGeometry out of GDALDataset

2022-10-14 Thread Johannes Paul
Hello, Is there a way to get the OGRGeometry object out of a GDALDataset object in c++? My GDALDataset is a raster and i would like to get the associated OGRPolygon out of it. I've used: srcDataset->GetLayer(0)->GetNextFeature()->GetGeometryRef()->toPolygon() However i keep getting segmentation fau

[gdal-dev] Create tile index by time?

2022-10-14 Thread Stefan Gofferje
Morning all! I have a little hobby project where I pull Sentinel imagery of my town's area and create various composites. I run that every morning, so I have a bunch of raster images which represent the same geographic area at different times. I recently also started playing with mapserver and