You can try converting each layer to data frame instead of all the layers
at once time. You can do this using a loop.
r <- terra::rast(file.nc)
df<- NULL
for( i in 1:nlyr(r)){
r1 <- r[[i]]
r1.df <- as.data.frame(r1,xy =T)
df <- rbind(df, r1.df)
}
Best regards
Vlad Alex
gb_1km_daily_20150101-20150131.nc",
subds = "yourvarname")
or with raster library:
r1 <- raster::brick("chess-met_precip_gb_1km_daily_20150101-20150131.nc")
same,
r1.sub <- raster::brick("chess-met_precip_gb_1km_daily_20150101-20150131.nc",
varname = "your
NuRD/file6583fb9afa9/stars_0.4-4.tar.gz’
had non-zero exit status
>
Best regards,
Vlad Alexandru AMIHAESEI
[[alternative HTML version deleted]]
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Hi Advait,
if you want to keep all time properties, the package timetools may
help you to do so.
This package defines the class TimeIntervalDataFrame which can be a
good start to look at.
Best,
Vlad
-
Vladislav Navel
SCAN
www.scan-datamining.com
Le 13/03/2012 11:03, Advait Godbole
))
sum.val - sum (tapply(x[x0], grps[x0], min))
return (c(nb.grps, sum.val))
} else {
return(c(0, 0))
}
}
calc (cntneg, theFun)
Best,
Vlad
---
Vladislav Navel
SCAN
www.scan-datamining.com
Le 08/03/2012 11:46, sajid pareeth a écrit :
Dear all
I
=FALSE)
spPolys - SpatialPolygons(xx)
polyData - data.frame (test=1:length(spPolys)) # actually the datayou
want to use
dfReturn - SpatialPolygonsDataFrame (spPolys, polyData, match.ID=FALSE)
Best,
Vlad
---
Vladislav Navel
SCAN
www.scan-datamining.com
Le 07/03/2012 04:22, Benjamin Stewart
(sorry about that), so if the
documentation is not understandable, ask me.
I hope this can help.
Vlad
Vladislav Navel
Ingénieur d'études
Atmo Poitou-Charentes
05 46 44 83 88
Le 14/10/2011 07:24, Mathieu Rajerison a écrit :
Hi,
And will contiguous raster cells have the same value