Questions tagged [stars]
The R package providing classes and methods for reading, manipulating, plotting and writing spatiotemporal data cubes.
24
questions
0
votes
0answers
10 views
st_warp a star_proxy object with multiple files
I have a s stars_proxy object with 12 attributes in files, i want to reproject each files in a new crs and resolution stars_proxy object, i use this command in R :
r <- stars::st_warp(src = s, dest ...
0
votes
0answers
8 views
Difference in CRS when reading with read_stars und read_ncdf
I have a NetCDF file that I would like to load into a stars object.
The two functions st::read_stars and st::read_ncdf both load it in, but the data ends up in different locations when I plot it using ...
1
vote
0answers
27 views
What is the canonical way to perform raster calculations using R stars?
I have a data cube that I am handling in the R stars package. It has 3 dimensions (latitude, longitude and time) and 2 attributes (northwards and eastwards components of wind). A dummy equivalent ...
0
votes
1answer
31 views
abind Error when converting a .tif raster data object to an sf object using st_as_sf() in R
When I try to convert a .tif raster data object to an sf object using st_as_sf() in R, I get the following error:
Error in abind(x, along = length(dim(x)) + 1) : along must be between 0 and 2
Does ...
2
votes
1answer
41 views
Why krige function return all NA's values for the prediction when trying to interpolate a variable?
I am trying to build a map with kriging interpolation for a variable. When running the krige function the return is all NA values and the kriged map plot only one value. I couldn't find out why and I ...
0
votes
0answers
45 views
Making stars::st_rasterize equivalent to raster::raster using meuse.grid
I'm trying to learn my way around stars. Given the meuse.grid data, how can I make a stars object that has the same dimensions, extent, etc. as I can get from going from SpatialPointsDataFrame to ...
1
vote
0answers
34 views
Getting an sf boundary of a raster or stars object
I am trying to figure out the best approach to convert an raster or stars object into its boundary returned as an sf object. I think I have a pretty decent handle on how to do it with raster (shown ...
1
vote
1answer
101 views
Convert large stars object to sf object
My ultimate goal is to convert landcover raster (.tif) objects to an sf object representing the raster's grid and the original values of each cell within each geometry. I have been able to do this for ...
0
votes
0answers
58 views
Extract all multiband pixel values inside the polygons from stars-proxy
I use stars and sf packages of R programming language for satelite image manipulation and data mining (with tidyverse packages).
I try to extract all pixel values from stars proxy object inside the ...
1
vote
1answer
66 views
temporal rolling mean in stars
I need to find temporal rolling means (or other stats - max, min, etc) in stars spatial obects. For stars raster cubes, this would be directly analogous to the calc(rasterBrick, function(x) movingFun(...
1
vote
1answer
139 views
Error while converting raster to ploygon using stars package
I am trying to convert raster to polygon using stars package as suggested in this and this. I have used the following code
library(raster)
library(stars)
library(sf)
library(magrittr)
#read raster
r1&...
1
vote
1answer
136 views
Exporting ASCII file from 'stars' package
I'm playing around with the stars package (0.4-1) in R, but I can't get how to export a grid object in ASCII format, as the traditional writeRaster(r, filename=file.path(dir.list[j], 'test.asc'), ...
2
votes
1answer
112 views
Raster file statistics in R (stars and raster package) don't match gdalinfo
I'm trying to read a gridded binary file in R using read_stars. Comparing the stats from stars package to the stats from gdalinfo, I could see differences in min/max values. Below is the sample code.
...
0
votes
1answer
128 views
Read NetCDF file with read_stars in R
I am trying to work with cordex climatic data (curvelinear and rotated data) in the stars package, and I am having a lot of issues, so decided to simply reproduce online example, and I fail reading in ...
3
votes
1answer
95 views
Generate a stars object from tabular spatial temporal data?
I have tabular data with columns longitude, latitude, time, and intensity, like so:
library(lubridate)
## mock example data
ex <- tibble(longitude = 1:10,
latitude = 10:1,
...
2
votes
1answer
102 views
r stars: calculate ndvi for every time step on stars proxy
I would like to calculate the ndvi on every time step of a 4d data cube (x, y, bands, time) in on an stars proxy object.
Do you have any ideas of how to achieve this on the proxy object?
My current ...
1
vote
0answers
99 views
Copernicus ERA5 crs/projecition for NetCDF
i'm trying to work from netCDF from this link:
https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form
but i can't find the correct EPSG.
In documentation page i ...
3
votes
1answer
496 views
stars package in R. How to add time dimension in a stars object?
I have a raster stack (my_stack) and I converted it into a stars object:
mystar <- st_as_stars(my_stack)
That results in:
stars object with 3 dimensions and 1 attribute
attribute(s), summary of ...
2
votes
0answers
152 views
Stratified ordinary kriging using gstat and stars packages
The gstat package supports stratified ordinary kriging using sp objects as shown at https://github.com/r-spatial/gstat/blob/master/demo/examples.R (the code below).
This year (https://r-spatial.github....
1
vote
2answers
309 views
NetCDF time aggregation using R stars package
I have a monthly temperature netCDF dataset with aprox 20 years and I'd like to calculate a monthly climatology (average of all Jan, Feb, Mar...).
I know I can do this with stackapply from package ...
4
votes
1answer
347 views
Convert XYZ to stars curvilinear grid and plot
I have data in XYZ format (long, lat, elevation) that does not fit on a rectangular grid, but I would like to plot it as grid cells of different sizes.
To illustrate what I mean, the last two plots ...
1
vote
0answers
245 views
Creating map projection for net of geodesic dome?
I am working on a project that involves mapping the night sky's northern hemisphere onto the exterior surface of a small (16 cm radius) geodesic dome. The problem is, I don't know what map projection ...
1
vote
0answers
299 views
R stars package and set dimensions to a NetCDF file
Related: GDAL and reprojecting a NetCDF file
Following above, I would like to reproject a NetCDF ( NOAA GOES-16 file ) using R stars package.
However, I have been unsuccessful at the first step; I ...
2
votes
1answer
1k views
Error from aggregate function from the R package sf
I am currently reading the new (and free!) book by Edzer Pebesma and Roger Bivand (you can find it here).
In chapter four they describe how one can aggregate all the measurement stations within one ...
