hi all,
i have one question ***does it matter if the KNN spatial weights matrix is
asymmetric or symmetric for linear spatial autoregressive models? *** i am
getting slightly different findings depending on the weights matrix
specification for the Lagrange multiplier diagnostic for spatial depende
R-sig-geo group,
I have found similar problems with the BBMM package in R. I ran my datasets
through brownian bridge movement model as implemented in Animal Space Use,
http://www.cnr.uidaho.edu/population_ecology/animal_space_use.htm and the R
package adehabitat,
http://cran.r-project.org/web/p
see this link -
http://r-sig-geo.2731867.n2.nabble.com/potential-BUG-in-BBMM-td5289846.html
--
View this message in context:
http://r-sig-geo.2731867.n2.nabble.com/Re-potential-BUG-in-BBMM-tp5840462p5844166.html
Sent from the R-sig-geo mailing list archive at Nabble.com.
___
Christian,
Could be quicker if you make a matrix out of the array (or, if
possible, avoid the array in the first place). See below
library(raster)
# create an array
m <- matrix(1:100, ncol=10, byrow=TRUE)
a <- array(list(m, m, m, m, m))
# empty RasterBrick
x <- brick(ncol=10, nrow=10)
# make m
Dear all,
I needed to convert an array into a RasterBrick object. The output did
not need to be written to memory, as assessed by canProcessInMemory().
The only way to put the data into the RasterBrick that worked for me
was layer by layer:
for (l in 1:nlayers(x)) x <-setValues(x,values=a
On Thu, 16 Dec 2010, Agustin Lobo wrote:
Roger,
Your new rgdal version from rgdal_0.6-32.tar.gz works fine once I've
removed gdal1.8dev (using gdal1.7.3 now):
overwrite_layer=T works and the files go to the directory indicated by dsn
I've now released rgdal_0.6-33 to CRAN, so we ought to be r
Hi Robert
Looking over your code, Wouldn't it be much faster
to use Rgooglemaps or something where the map is already rasterized?
The newer versions of R have the lowlevel raster function for plotting
images.
Also, even in base graphics you can use alpha blending so for your
network lines
lines(rt
Browsing through your script, another issue I noted is that you
interpolate in long/lat, but the software assumes Euclidian distances
make sense as a distance measure in this space. They don't.
Wbr,
On 12/16/2010 03:52 PM, Tarek Benamara wrote:
> Hallo,
>
> I have used rgdal to generate a GeoTIF
Hey Tarek,
I guess you have this issue because the g object has two bands
(the estimated mean and variance or sd of your kriging, someone will
correct me). So your writeGDAL() should be
writeGDAL(g[,,"band1"], "ordinary_kriging_1.tif",
drivername="GTif
Hallo,
I have used rgdal to generate a GeoTIFF for my krige Interpolation but the
output image is a 2 Band image und without RGB Color, my Question is haw to
give the RGB Color to my GeoTFF and how to increase resolution with rgdal.
And this is my R-File
## load some libraries first:
librar
On Fri, Dec 10, 2010 at 1:46 PM, Massimo Di Stefano
wrote:
> the R packages installed on the Osgeo-Live are :
>
> https://svn.osgeo.org/osgeo/livedvd/gisvm/trunk/app-conf/R/installRpackages.r
>
> packagelist is what actually is installed
>
> while packagelist2 is commented out (more complete list
On Thu, 16 Dec 2010, Roger Bivand wrote:
On Thu, 16 Dec 2010, Wang, Kevin (SYD) wrote:
Hi all,
I've managed to draw a map using the following code:
lga1 <- readShapePoly("D:/Documents and Settings/kevinwang/My
Documents/SMS/ASGC/LGA10aAust", IDvar = "LGA_CODE10")
plot(lga1[lga
Thanks, Roger,
Got it to work.
Got another question, is it possible to "enlarge" the map either using the plot
function, or another function?
Cheers
Kevin
sent from my hTC Desire handset
- Reply message -
From: "Roger Bivand"
Date: Thu, Dec 16, 2010 9:22 pm
Subject: [R-sig-Geo]
On Thu, 16 Dec 2010, Wang, Kevin (SYD) wrote:
Hi all,
I've managed to draw a map using the following code:
lga1 <- readShapePoly("D:/Documents and Settings/kevinwang/My
Documents/SMS/ASGC/LGA10aAust", IDvar = "LGA_CODE10")
plot(lga1[lga1$STATE_CODE == 7,], border="blue")
-- Forwarded message --
From: Agustin Lobo
Date: 2010/12/16
Subject: Re: Problem with writing parameters in overlay()
To: "Robert J. Hijmans" , roger.biv...@nhh.no,
alobolis...@gmail.com
Robert,
The weird problem at writing raster layers in some cases is also
solved once I remov
Roger,
Your new rgdal version from rgdal_0.6-32.tar.gz works fine once I've
removed gdal1.8dev (using gdal1.7.3 now):
overwrite_layer=T works and the files go to the directory indicated by dsn
Thanks
Agus
> setwd ("/media/FREECOM_HDD/GRAVI")
> require(rgdal)
Loading required package: rgdal
Geosp
On 12/15/2010 09:22 PM, Barry Rowlingson wrote:
...
Ah okay - feel free to ask on the r-sig-geo mailing list. It should
be trivial for Linux and Windows systems, I think there are some
issues with Mac OS, but the boys on r-sig-geo will be helpful assuming
you've had a go yourself first accordin
On 16/12/10 09:06:33, dorina.lazar wrote:
> How to create a shapefile including my data about some socio-demographic
> indicators (7), for about 30 countries around the world (to be read with
> readShapePoly)?
Hi Dorina,
I'm not quite sure what you have and what you want.
I assume that you'd li
Hi Hadley,
That would great if rasters could be printed more efficiently. Could you
send me an e-mail, or post on the r-sig-geo list when you are ready to
receive any contributed code?
Paul
On 12/15/2010 10:20 PM, Hadley Wickham wrote:
This is something that I hope the student would explore
19 matches
Mail list logo