Re: [gdal-dev] Issue with Custom Transformation Pipeline in GDALReprojectImage

2025-01-09 Thread Even Rouault via gdal-dev
Anyway a proper pipeline I use for years is: +proj=pipeline +step +inv +proj=krovak +lat_0=49.5 +lon_0=24.8 +alpha=30.2881397527778 +k=0. +x_0=0 +y_0=0 +ellps=bessel +step +inv +proj=hgridshift +grids=Slovakia_JTSK03_to_JTSK.gsb +step +proj=krovak +lat_0=49.5 +lon_0=24.8

Re: [gdal-dev] Issue with Custom Transformation Pipeline in GDALReprojectImage

2025-01-09 Thread Martin Ždila via gdal-dev
On Thu, Jan 9, 2025 at 4:53 PM Even Rouault wrote: > > > Regarding the pipeline, it could be a bad copy-paste or garbled by > > ChatGPT 😀. > > I wouldn't trust ChatGPT to come with geodesic-relevant pipelines... If > you aren't really sure about what you are doing, you'd better let PROJ > infer a

Re: [gdal-dev] Issue with Custom Transformation Pipeline in GDALReprojectImage

2025-01-09 Thread Even Rouault via gdal-dev
Regarding the pipeline, it could be a bad copy-paste or garbled by ChatGPT 😀. I wouldn't trust ChatGPT to come with geodesic-relevant pipelines... If you aren't really sure about what you are doing, you'd better let PROJ infer a reasonable pipeline. I'm also questioning the use of the gr

Re: [gdal-dev] Issue with Custom Transformation Pipeline in GDALReprojectImage

2025-01-09 Thread Martin Ždila via gdal-dev
Thank you Even, I went more low-level and successfully used GDALChunkAndWarpImage . I also posted answer on SO . Regarding the pipeline, it

Re: [gdal-dev] Issue with Custom Transformation Pipeline in GDALReprojectImage

2025-01-09 Thread Even Rouault via gdal-dev
Martin, GDALReprojectImage() expects that a correct geotransform is already set on the target dataset. You may use GDALSuggestedWarpOutput2() to get a guessed geotransform Otherwise for full automated way, you may also using the C function GDALWarp(). By the way, your below pipeline looks