Hello,
How package include function quantile.Spatial()?
Thank You
Agnieszka
[[alternative HTML version deleted]]
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
I just did a test using dropLayer on a similar brick of global
temperatures
class : RasterBrick
dimensions : 36, 72, 1332 (nrow, ncol, nlayers)
resolution : 5, 5 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84
values : no
Hello Christian:
Yes indeed, my experience is that Raster package operations on very
large datasets take a long time,
on either Linux or Windows desktops.
As an alternative to the Raster package, I have had good experience
using the GDAL raster image utilities,
embedded within shell or Python
Hi there,
I have a Digital Elevation Map (mnt50.asc) for which I would like to
calculate slope, curvature etc. using rsaga.local.morphometry (from
RSAGA).
I am on Ubuntu. I installed saga in /usr/bin (i.e. I have saga_cmd here)
The path to mnt50.asc is /home/kati/projects/cemagref/maps
I get the
Dear all,
The package 'raster' provides excellent tools to work on geographic
data. Subsetting a RasterBrick object, however, takes ages.
I have the following RasterBrick object:
Temp.CRU.cropped
class: RasterBrick
filename:
dimensions:43, 99, 1200 (
Edzer, I will try also your suggestion. Thank you.
Best wishes,
-Original Message-
From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org]
On Behalf Of Edzer Pebesma
Sent: Friday, July 01, 2011 4:27 PM
To: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] space-time
Pinar / Tom,
another route you might follow is (improvised):
require(spacetime)
prec_data$YearDate = as.Date(paste(prec_data$Year, "1", "1",sep="-"))
x = stConstruct(prec_data, c("X", "Y"), "YearDate")
class(x) # STIDF
x = as(x, "STSDF") # will sort out common locations & times
Bests,
On 07/01/
Thank you Tom! Now index is correct, comprised from 5365 record. I created
STSDF object with this index.
Best regards,
Pinar
-Original Message-
From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org]
On Behalf Of Tom Gottfried
Sent: Friday, July 01, 2011 3:35 PM
T
You could try
index <- cbind(as.integer(as.factor(prec_data$Station)),
as.integer(as.factor(prec_data$Year)))
Tom
Am 01.07.2011 14:02, schrieb Pınar Aslantas Bostan:
Hi Edzer,
Thank you for your mail.
Prec data has location and time information.
prec_data[1:10,]
Station X
Hi Edzer,
Thank you for your mail.
Prec data has location and time information.
> prec_data[1:10,]
Station X Y Prec Year Month Day
117020 693894 4931348 1274.4 1970 1 1
217020 693894 4931348 1134.3 1971 1 1
317020 693894 4931348 940.7 1972 1 1
4
Hi Pinar,
as you say, some data are missing, but the way you create index suggests
none are missing. Does your table "prec" have information for each
observation on which location and at which time it has been taken? Could
you maybe show how prec[1:10,] looks like?
Wbr,
On 07/01/2011 01:12 PM, P
Dear all,
I am working about space-time (ST) kriging. I want to make ST ordinary
kriging with sparse data set. Dependent data is annual precipitation
measured from 257 meteorological station. Observation period starts from
1970 and finishes at 2008, so temporal entity of my study contains 39
ye
12 matches
Mail list logo
