[gdal-dev] Reprojecting image with lat/lon bands - best practice?

Even Rouault even.rouault at spatialys.com
Fri Oct 18 10:50:42 PDT 2024


Conrad,

Try something like:

gdal_translate input.tif imagery.vrt -b 3

gdalwarp imagery.vrt imagery_warped.tif -geoloc -to X_DATASET=input.tif 
-to X_BAND=2 -to Y_DATASET=input.tif -to Y_BAND=1 -to PIXEL_OFFSET=0 -to 
PIXEL_STEP=1 -to LINE_OFFSET=0 -to LINE_STEP=1 -to SRS=EPSG:4326 -a_srs 
EPSG:4326  -overwrite

Obviously I have most certainly got something wrong in the above, but 
hopefully with a tiny tweaking that should put you on the right track.

Reference: https://gdal.org/en/latest/development/rfc/rfc4_geolocate.html

Even

Le 18/10/2024 à 12:37, Javier Jimenez Shaw via gdal-dev a écrit :
> Is it an actual grid? in the meaning of having constant step size in X 
> and Y.
> In that case the geolocation is just the corner and the x and y sizes. 
> You can convert to a georeference raster, and warp it.
> If it is not the case, you have something more like a 2D pointcloud, 
> or a bunch of poins in a strange vector format.
>
> On Fri, 18 Oct 2024 at 12:20, Conrad Bielski via gdal-dev 
> <gdal-dev at lists.osgeo.org> wrote:
>
>     Hello GDAL-experts,
>
>     normally when I use GDAL for reprojecting imagery, the projection
>     information that I use is the source spatial reference (SRS)
>     associated with the imagery. However, now I have imagery which is
>     lat/lon geographic and I have two separate bands which also carry
>     the pixel geographic information. So the following raster inputs
>     all the same size:
>     1. Band 1 = latitude
>     2. Band 2 = longitude
>     3. Band 3 = imagery
>
>     The question I have is how best to integrate this information into
>     a reprojection workflow?
>
>     I presume that gdalwarp is the best option here, but how can I
>     take advantage of the individual pixel location information
>     (rather than just the extents for example)? I know that I can
>     mosaic into an existing file that I have already created in the
>     target projection. Is this the best way to apply gdalwarp in this
>     context?
>
>     I'm just wondering what is the best way to integrate the lat/lon
>     pixel information into my warping using gdalwarp.
>
>     Thanks in advance for your help,
>     Conrad
>     _______________________________________________
>     gdal-dev mailing list
>     gdal-dev at lists.osgeo.org
>     https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
Butcher of all kinds of standards, open or closed formats. At the end, this is just about bytes.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20241018/eb2e896a/attachment-0001.htm>


More information about the gdal-dev mailing list