Chris Toney

Would the return of `ds$fetch(1)` in place of `ds$getNextFeature()` change anything here? Currently, the return value of `ds$getNextFeature()` is basically `ds$fetch(1)` from the current cursor position, but with list column(s)...

Right, `get_pixel_line()` just applies the inverse geotransform and does no bounds checking. Thanks for pointing this out. I imagined it used as a lower level function with caller doing validation...

>should be a method on the class though, ds$get_pixel_line(xy) Yes that makes sense for it to exist in both forms. Done in #341. That PR also allows `xy` to be...

This seems like a worthwhile improvement that could be useful in general, but I haven't thought through all the details yet. >Could we have the destination dataset passed in as...

`GDALAutoCreateWarpedVRT()` might be useful for this, https://gdal.org/api/gdalwarp_cpp.html#_CPPv423GDALAutoCreateWarpedVRT12GDALDatasetHPKcPKc15GDALResampleAlgdPK15GDALWarpOptions Earlier I added the exported method `GDALRaster::setFilename()` as mentioned above, but it has ```cpp if (fname_in == "" && filename != "" &&...

This has been implemented for `warp()` in #469. Please see comments and the code example there. Possibly related is `autoCreateWarpedVRT()` which also returns a dataset object (#471). `GDALRaster::setFilename()` is now...

Unfortunately there is not a way to run `gdal2tiles.py` using `gdalraster`. The GDAL utilities that are exposed in the GDAL C API are generally available (most of "core" executables, with...

Thanks, I didn't mean to imply that porting `gdal2tiles.py` is not feasible in general, but kind of sounded like that. What I should have said is, I don't think it's...

I think I'm following but apology in advance if I just add noise here. It sounds like you do need just a subset of `gdal2tiles.py`, only to generate the tiles....

@wiesehahn, thanks for that description. Sorry, I was slow to understand the whole workflow but I think I'm getting there now. Your hint of possibly using VRT wasn't what I...