Krzysztof Dyba

Note that the tools in `terrain()` are originally implemented in [GDAL](https://gdal.org/programs/gdaldem.html), which `{terra}` just wraps. If we want to avoid new dependencies in `{terra}`, that means someone has to write...

Wouldn't it be easier to just use the existing tools in [SAGA](https://saga-gis.sourceforge.io/saga_tool_doc/9.0.3/ta_hydrology_29.html) or [GRASS](https://grass.osgeo.org/grass82/manuals/r.flow.html) from R?

I suspect this problem may be due to processing raster in blocks from disk. If the raster is too large to be loaded into memory, it is processed sequentially in...

> Or is the n argument in mem_info() supposed to be the number of layers in the SpatRaster? `mem_info()` gives the required memory for the whole object -- all layers,...

I think I can reproduce this and this is some bug in `project()`: ```r library("terra") # 1.7.55 gdal(lib = "") #> gdal proj geos #> "3.6.2" "9.2.0" "3.11.2" set.seed(1) n...

@mdsumner, I just wanted to create some reprex. In my example, for identical inputs, shouldn't the result be identical regardless of the number of layers? I would expect that the...

I tried to reproduce this issue using the code below, but it works ok. Maybe this problem is more specific to the dataset? ```r library("terra") terraOptions(memfrac = 0.1) f

FWIW: It seems that it doesn't work in `{stars}` either (I checked on Windows 10 and Ubuntu 22.04): ```r url = "https://www.unidata.ucar.edu/software/netcdf/examples/sresa1b_ncar_ccsm3-example.nc" stars::read_stars(url) #> area, trying to read file: NETCDF:"/vsimem/http_2/sresa1b_ncar_ccsm3-example.nc":area...

BTW: Shouldn't `crds` be `crds[, 2:3]` to omit column with ID? Nevertheless, the crash still occurs on Windows.

@paleolimbot, do you have any tips on how to load a .parquet file with geometry? I once tried to load using `{arrow}`, but I couldn't find a way to convert...