Krzysztof Dyba

You can also consider: - [openeo](https://cran.r-project.org/web/packages/openeo/index.html) - [RStoolbox](https://cran.r-project.org/web/packages/RStoolbox/index.html) - [CDSE](https://cran.r-project.org/web/packages/CDSE/index.html) - [rsat](https://cran.r-project.org/web/packages/rsat/index.html)

I compared this example using `terra` and it basically works (probably doesn't draw that one missing geometry?). ```r library(terra) f = "https://naciscdn.org/naturalearth/50m/cultural/ne_50m_admin_0_countries.zip" v = vect(paste0("/vsizip/vsicurl/", f), what = "geoms") v...

@Kodiologist, I think the easiest way is to check the [NEWS file](https://github.com/rspatial/terra/blob/master/NEWS.md) (version 1.8-15): > improved estimate of available memory on linux systems https://github.com/rspatial/terra/issues/1506 by Cedric Rossi (or check the...

```r project(r, "EPSG:3857", method = "bilinear", gdal = TRUE, filename = output_terra) #> class : SpatRaster #> dimensions : 12685, 1048, 3 (nrow, ncol, nlyr) #> resolution : 38239.51, 38239.55...

You probably have a mistake in your code somewhere. I tested this simplified version and the result looks ok: ```r library(sf) sampling_points = st_read("site_data_all_gck.shp") weather_stations = st_read("gcr_weather_stations.shp") idx = st_nearest_feature(sampling_points,...

Both Sentinel-1 and Sentinel-2 products are directly supported by GDAL. Maybe the Landsat driver should be there too? 1. https://gdal.org/en/stable/drivers/raster/safe.html 2. https://gdal.org/en/stable/drivers/raster/sentinel2.html

In `crop.jl`: ``` rast = Raster(RasterStack(raster_files; name=band_names, lazy=true)) x = crop(rast; to=ext) extent(x) #> Extent(X = (598515.0, 727485.0), Y = (5.682105e6, 5.780985e6)) ``` xmax and ymax are different, because `rasters.jl`...

In `terra`: ``` x = aggregate(ras, fact = 3, fun = "mean") global(x[[1]], "mean", na.rm = TRUE) #> 9767.289 x[1000,1000,1] #> 10125.44 ``` In `rasters_jl`: ``` x = Rasters.aggregate(mean, raster,...

`extract-points.jl`: Everything is missing. Output should have values. ``` @NamedTuple{geometry::Union{Missing, Tuple{Float64, Float64}}, B1::Union{Missing, UInt16}, B10::Union{Missing, UInt16}, B11::Union{Missing, UInt16}, B2::Union{Missing, UInt16}, B3::Union{Missing, UInt16}, B4::Union{Missing, UInt16}, B5::Union{Missing, UInt16}, B6::Union{Missing, UInt16}, B7::Union{Missing, UInt16},...

In `terra`: ``` x = lapp(ras, fun = ndvi) global(x, mean, na.rm = TRUE) #> 0.285815 ``` In `rasters_jl`: ``` x = get_ndvi(red, nir) mean(skipmissing(x)) #> 0.331073887918443 ```