Ok, thanks for the answer. It must be a problem of disk space.
I'll try what you say.
Best,
Mathieu
2011/11/9 Robert J. Hijmans
> Mathieu/Carsten,
>
> Carsten's suggestion is also good:
>
> elev<-overlay(detm@layers[[1]],detm@layers[[2]],
> fun=function(a,b) {return(a-b)} )
>
> although I wou
Hi
I am plotting series of raster* objects with 'plot'. The raster* objects are
change maps with variable scales.
I want to plot them all with the same scale centred on zero such that it's easy
to see where there is a positive or negative response.
I am using zlim=c(display_min,display_max) to
Hello!
Trying to decide the best way to store a large amount of data in a
function I am writing using the raster package.
The goal is speed. I need to flexibly and efficiently handles cases
where 10 to 5000 copies of a raster are required. The raster that
requires replication would typ
Leni,
If you have raster data in Arc and export as GRID it will create a database
(folder, including a parallel folder called info) that you can treat like a
file and open in R. Alternatively you can export to tif, img (erdas) and
other formats (avoid ascii) that raster can read (most of them via
Mathieu/Carsten,
Carsten's suggestion is also good:
elev<-overlay(detm@layers[[1]],detm@layers[[2]],
fun=function(a,b) {return(a-b)} )
although I would do:
elev<-overlay(detm, fun=function(a,b) {return(b-a) } ,
filename="disk/with/space.tif")
Robert
On Wed, Nov 9, 2011 at 8:46 AM, Carsten Ne
Mathieu,
The error message says that something is going wrong when writing the
results to disk. This smells like bug, but it is almost certainly because
there is no more disk space available in your temp folder/disk. To avoid
using a temp file (used for large files only), you can use calc with a
f
Please list what version of raster you are using. A bug was submitted and
corrected on this issue
sometime ago so I suspect you may have an old version.
sessionInfo() please when you have a question.
On Tue, Nov 1, 2011 at 8:36 AM, Aidin Niamir wrote:
> Dear all
>
> I am experiencing a str
Am 09.11.2011 15:39, schrieb Mathieu Rajerison:
> Hi,
>
>
> I have a RasterStack detm composed of a dem and a dtm.
>
> I wanted to apply the difference between dem and dtm so as to get the
> elevation.
>
> But I got the following message (in French):
> Erreur dans result[, i]<- readBin(raster@fil
Roger Setting the repositories as suggested worked well, for these and one
other package I'd been having trouble with. gLength worked fine once I got them
working. I was only talking about single-ring polygon objects.
Robert The shapefiles were indeed in lat/long, so the fact that
library("
Hi,
I have a RasterStack detm composed of a dem and a dtm.
I wanted to apply the difference between dem and dtm so as to get the
elevation.
But I got the following message (in French):
Erreur dans result[, i] <- readBin(raster@file@con, what = dtype, n =
ncols, :
le nombre d'objets à remplacer
Hi,
I want to with raster data generated from a 3D pointcloud in R with the
package "raster". But I have problems exporting my raster dataset from
ArcMap. By "export data"/GRID I get a file with an *.aux extension. What I
need is *.grd.
Can anybody help me out?
Cheers,
Leni
--
View this mess
I have two shapefile layers with identical prj file, but when I open
the first one
in QGIS Properties/General I see "EPSG:23031 - ED50 / UTM zone 31N"
defined as
+proj=utm +zone=31 +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +units=m +no_defs
(which is what I want)
while when I open the second one it
assuming your geometry column is 'the_geom': while querying
SELECT ... FROM ... WHERE the_geom @ ST_geomfromtext(...);
with ST_geomfromtext() you can set up your "custom bbox".
But of course you can load everything into R and then use subset() if your data are in a data.frame
or use over() with
Hi Anthony,
if I am not mistaken, David already started to translate his filter into
R. At least I know that he wanted to do so. If you really are keen, we
could give it a try. On the other hand the filter was implemented in
www.movebank.org
Theoretically, you can upload your data to movebank,
Hi,
once a Postgis layer X is retrieved, how can I juggle with subsets of X?
eg. I am working with bubble plots on a EU dataset but I want to do test it
only using a 30x30Km in the nearby of Rome and I have no attribute data but
the value I use for the bubbles. Should I use "subset"? Can I set a
These filters are easy to implement and are available in contributed R
packages, but they are not very good. When you say "satellite"do you mean
Argos or GPS?
Argos recently updated their proprietary algorithms to give better
estimates but it it is still a rather black box.To me the issue is about
I would like to remove a second order trend from elevation data before kriging.
Was thinking on using surf.ls to fit the trend but not sure how to pass it onto
the krige argument. Any suggestions? Thanks
[[alternative HTML version deleted]]
__
17 matches
Mail list logo