How to work at the native resolution?

[liuyadong]

Hi, if not specify an view in raster_cube, the default overview has a coarse resolution, which maybe not what we want. So what’s the easiest way to work on the native resolution? Suppose I want to stitch files without any other operations.

[appelmar]

Hi, the idea of the default view was simply to provide a quick look at the data. I think that in most cases it does make sense to let users explicitly tell "I want to work at full resolution", because this may involve long computation times. The only way to do this at the moment is to provide a custom data cube view, and setting dx, dy, dt, and srs accordingly.

For image stitching (spatial mosaics?), gdalcubes calls gdalwarp to reproject, rescale, resample, and crop images according to the grid defined in the data cube view. It can only work if

• all images have the same SRS (no reprojection needed)
• all bands have the same spatial resolution (no rescaling needed),
• pixels of all images must be aligned with pixels of the target grid (but not overlap parts of pixels) (no resampling needed),
• the extent of the data cube view covers the extent of all images (no cropping needed)

What dataset are you working with?

[liuyadong]

It does make sense to have a default low-resolution view for a quick try, but there is no straightforward way to work at the full resolution. I think adding an option would make it much easier. I'm trying to stitch MODIS tiles and export as an tif file using this package (not MRT or HEG) 😄