(help request) using ArchGDAL.gdalwarp() for resampling a raster dataset

Specific Domains Geo
,
2

Okay! I think I figured it out…?

ds = ArchGDAL.read(path) do source
    ArchGDAL.gdalwarp([source], ["-t_srs", "$(wkt)",
                               "-te", "$(west_bound)", "$(south_bound)",
                               "$(east_bound)", "$(north_bound)",
                               "-te_srs", "$(wkt)",
                               "-tr", "$(step(x_dim))", "$(step(y_dim))",
                               "-r", "cubic"]) do warped
        band = ArchGDAL.getband(warped, 1)

        ArchGDAL.read(band)
    end
end
2 Likes