Dear All,
I want to subset a large spatial grid (2000*2000) into many smaller
grids for further processing in my simulation models. Could anyone do it
before and any tips will be highly appreciated?
Library(rgdal)
Map - readGDAL(yong.img) # one image from erdras imagine
Map$va - Map$band1
Hi everyone:
with the package spatstat I can do the kernel density estimation in the
function density.ppp. but How do I use R for Adaptive kernel density
estimation.
Thanks in advance
Gordon
2008-5-20
___
R-sig-Geo mailing list
Thanks Edzer and Hengl. It works.
But R is really shaking when large grids are processed. We wish there is
improvement of sp in the future like GIS handling grids.
Regards
Yong
From: Edzer Pebesma [mailto:[EMAIL PROTECTED]
Sent: Tue 20/05/2008 6:11 PM
To:
Yong, I would have loved to see your R shaking.
This argument has been heard more often here. When processing massive
grids, please have a look at the function getRasterData in package
rgdal. It allows you to _read_ only portions of a grid.
--
Edzer
Yong Li wrote:
Thanks Edzer and Hengl. It
On Tue, 20 May 2008, Edzer Pebesma wrote:
Yong, I would have loved to see your R shaking.
Unless you do seismology, of course?
This argument has been heard more often here. When processing massive grids,
please have a look at the function getRasterData in package rgdal. It allows
you to
Hi,
I am trying to display a geotiff using displayDataset command from rgdal .
which i can do pretty successfully but i cannot display axes for my map i
tried to use axes = TRUE and i got the error that axes is matched by multiple
arguments and if i try to set xaxs and yaxs as equal
On Tue, 20 May 2008, Monica Pisica wrote:
Hi,
I am trying to display a geotiff using displayDataset command from rgdal
. which i can do pretty successfully but i cannot display axes for
my map i tried to use axes = TRUE and i got the error that axes is
matched by multiple
If it is of any help, you might want to take a look at this guide:
Creating Maps for Publication using R Graphics
http://www.nceas.ucsb.edu/scicomp/GISSeminar/UseCases/MapProdWithRGraphics/OneMapProdWithRGraphics.h
tml
The second example demonstrates how to visualize a remote sensing band using
Hi,
Thanks a lot. This really does the trick. I need to figure now how i am
displaying the color bar ... so i am off to figure out image.plot ;-)
Monica Date: Tue, 20 May 2008 16:38:06 +0200 From: [EMAIL PROTECTED] To:
[EMAIL PROTECTED] CC: r-sig-geo@stat.math.ethz.ch Subject: Re:
Hi,
I am impressed by the second example . I will have to look at it closely
I am preparing a ppt for my office in which i want to show some of the
options R has in dealing with spatial data, map graphics and such and your
example is just perfect. The audience is almost virgin
Hi,
I'm working on Windows with 2GB RAM and am trying to run the errorsarlm
with 32 predictor variables and 10,071 records. I created my listw object
using the dnearneigh function with an upper distance limit of 3000m. I used
the following code with method=Matrix to deal with memory
On Tue, 20 May 2008, Sylvia Churcher wrote:
Hello all,
Being a new and cautious user of the rgdal library I am having trouble
confirming that my results are correct.
transforming vectors with long and lat values to easting and northing in zone
16, datum=NAD83
I'll obviously omit the entire
Thanks for your response. Here's what I have now.
On May 20 2008, Roger Bivand wrote:
A bit more context would help. The output of sessionInfo() at least (R and
package versions may matter), plus the preceding code, including the
definition of rams_dak.dnear3k.listw - is it
On Tue, 20 May 2008, [EMAIL PROTECTED] wrote:
Thanks for your response. Here's what I have now.
On May 20 2008, Roger Bivand wrote:
A bit more context would help. The output of sessionInfo() at least (R and
package versions may matter), plus the preceding code, including the
definition of
Dear list,
I tried to fit the model variogram to the experimental semivariogram. I used
the following code:
A-read.table(file.choose(new=FALSE), sep = ,skip=8,header=FALSE)
names(A) - c(x, y, B1,B2,B3,B4)
library(sp)library(gstat)coordinates(A) =
Hi,
I've used R for a few years but I'm a novice regarding the topics
discussed in this list. I would like to know if there is a package or
script that can easily access and manipulate the NARR (North American
Regional Reanalysis) database. I have a model that takes these data as
input at
16 matches
Mail list logo