Hi, we also have a bilinear sampling function here: lib/gis/sample.c (also bug affected?) -> raster_sample_bilinear() Probably it makes sense to reduce the number of functions and maybe move the r.proj functions to libgis for wider usage (e.g. r.resamp.interp)? Maybe a dumb question... I don't know the internals of these modules. Markus