Re: [R-sig-Geo] Reference for join-count model

2009-08-25 Thread Roger Bivand
Catmogs are online at: http://qmrg.org.uk/catmog/ (the series published by GeoBooks). Lots of useful stuff, a bit "legacy", but often that is what is needed. There are few join count references since the mid 1980's, so some time ago, quite undeservedly, I think, they are useful for categorical

Re: [R-sig-Geo] Kriging

2009-08-25 Thread Thomas . Adams
Maybe something like this: http://www.ars.usda.gov/sp2UserFiles/ad_hoc/1200SpatialWorkshop/01VinyardOverview.pdf - Original Message - From: milton ruser Date: Tuesday, August 25, 2009 7:50 pm Subject: Re: [R-sig-Geo] Kriging > Hi Kabeli, > > I never saw Brian Vinyard slides 40! :

Re: [R-sig-Geo] Matern

2009-08-25 Thread TANKISO THEJANE
Dear R users   Greetings from a novice. I would be grateful if you can help me complete the plot of a matern correlation function Matern Correlation function with   kappa =1,1.5,2,2.5,3 , phi = 1  The problem is that i can manage to get the result when t> 0 using the following code:   curve(ma

Re: [R-sig-Geo] Kriging

2009-08-25 Thread milton ruser
Hi Kabeli, I never saw Brian Vinyard slides 40! :-) It is accessible on a internet site? bests milton On Tue, Aug 25, 2009 at 6:56 PM, KABELI MEFANE wrote: > Hi all > > Please help and correct me, to predict Z(S0) at (0.5,0.5) given Z(S1) = 3 > at (0,0), Z(S2) = 5 at(0,1), Z(S3) = 6(1,0) and

[R-sig-Geo] Kriging

2009-08-25 Thread KABELI MEFANE
Hi all   Please help and correct me, to predict Z(S0) at (0.5,0.5) given Z(S1) = 3 at (0,0), Z(S2) = 5 at(0,1), Z(S3) = 6(1,0) and Z(S4) = 4(1,1). Let γ (h)=h^2  h<1    = 1 h>=1 using geo.   I did this: coords<-matrix(c(0,0,1,1,0,1,0,1), nrow=4, ncol=2) > data<-c(3,5,6,4) > coordata<-dat

Re: [R-sig-Geo] Generic mapping tools

2009-08-25 Thread Robert J. Hijmans
Ogbos, You could have a look at pointsToRaster in the raster package to count the number of points (events?) per raster cell. to install raster: install.packages("raster", repos="http://R-Forge.R-project.org";) Robert On Tue, Aug 25, 2009 at 10:56 AM, ogbos okike wrote: > Hello to everybody, >

[R-sig-Geo] Generic mapping tools

2009-08-25 Thread ogbos okike
Hello to everybody, I have a huge geographic data labeled latitude and longitude (global). I am attempting to use Generic mapping tools to plot the distribution of the events. I wish to obtain number of events per region on the map. Using factor function, I have summed up the number of events on

Re: [R-sig-Geo] Overlay on raster loses points when merging with data frame

2009-08-25 Thread Mehdi Khan
Ok it seems like the limits match. Therefore values inside this range are somehow being deleted. Is there any way I can determine which points are being deleted? On Tue, Aug 25, 2009 at 4:21 AM, milton ruser wrote: > Hi there, > > (1) try make a summary of X and Y and see if the minimum/maximum

Re: [R-sig-Geo] GDAL Error - .grd not recognised as a supported file

2009-08-25 Thread Darren Norris
Many thanks all for your tips and suggestions, Yes, I'm using the standard windows binary and unfortunately am rather clumsy when comes to PCs so will refrain from compiling / rebuilding myself . The spatial functionality developed in R is helping me a great deal and at this stage the options

Re: [R-sig-Geo] GDAL Error - .grd not recognised as a supported file

2009-08-25 Thread Michael Sumner
Hello, (I was trying to send this last night, so it's been covered by Barry and Roger, but I'm now sending it as a test of this email address subscription). Unfortunately you won't be able to read NetCDF with the CRAN version of Windows rgdal. If you build the package yourself against FWTools, or

[R-sig-Geo] transparency in raster images

2009-08-25 Thread Barry Rowlingson
I've just been hacking image.SpatialGridDataFrame to support transparency, either via specifying an alpha band or a constant value of alpha: image(foo,red=1,green=2,blue=3,alpha="band4",add=TRUE) image(foo,red=1,green=2,blue=3,alpha=0.5,add=TRUE) I don't think what I've done works properly with N

Re: [R-sig-Geo] Overlay on raster loses points when merging with data frame

2009-08-25 Thread milton ruser
Hi there, (1) try make a summary of X and Y and see if the minimum/maximum limits are included on your minimum/maximum limit of your grid; (2) try to see a frequency table of X/Y combination, to have sure if some of your X/Y combinations have frequencies>1. good luck, milton On Mon, Aug 24, 200

[R-sig-Geo] converting from images or spatialgriddataframes to kasc

2009-08-25 Thread Els Ducheyne
Dear list I have been trying to convert input data from either images or spatialgriddataframes to kasc for enfa modelling However I keep getting a strange error that the cell size is not the same in the different dimensions - see below. Is there anyone that can enlighten because I am rea

Re: [R-sig-Geo] GDAL Error - .grd not recognised as a supported file

2009-08-25 Thread Roger Bivand
I think Darren is using the standard Windows binary rgdal package, which does not include any GDAL drivers that depend on external DLLs, so does not include netCDF. Unless he wants to rebuild against the FWTools or OSGEO4W DLLs (as described in the file readme.WINDOWS in the package), another op

Re: [R-sig-Geo] Reference for join-count model

2009-08-25 Thread Marco Helbich
Dear Michael, this one comes to my mind: @BOOK{Goodchild1986, title = {Spatial Autocorrelation}, publisher = {Geo Books, Norwich}, year = {1986}, author = {Goodchild, M}, timestamp = {2007.05.14} } hope it helps marco > -- > > Message: 10 > Date: Tue, 25

[R-sig-Geo] Reference for join-count model

2009-08-25 Thread Haenlein, Michael
Dear all, I'm using the joincount.test function in the spdep package to test for the presence of spatial autocorrelation in a nominal variable. Given that I'm doing this analysis for a paper I'm currently working on, I would need a reference that explains the logic behind the join-count model. T