[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