GeoModels (version 2.2.1)

GeoNeighIndex: Spatial or spatiotemporal near neighbour indices.

Description

The function returns the indices associated with a given spatial (temporal) neighbour and/or distance

Usage

GeoNeighIndex(coordx,coordy=NULL,coordz=NULL,
coordt=NULL,coordx_dyn=NULL,distance="Eucl",neighb=4,maxdist=NULL,
maxtime=1,radius=1,bivariate=FALSE)

Value

Returns a list containing the following components:

colidx

First vector of indices

rowidx

Second vector of indices

lags

Vector of spatial distances

lagt

Vector of temporal distances

Arguments

coordx

A numeric (\(d \times 2\))-matrix or (\(d \times 3\))-matrix Coordinates on a sphere for a fixed radius radius are passed in lon/lat format expressed in decimal degrees.

coordy

A numeric vector giving 1-dimension of spatial coordinates; Optional argument, the default is NULL.

coordz

A numeric vector giving 1-dimension of spatial coordinates; Optional argument, the default is NULL.

coordt

A numeric vector assigning 1-dimension of temporal coordinates. Optional argument, the default is NULL then a spatial random field is expected.

coordx_dyn

A list of \(m\) numeric (\(d_t \times 2\))-matrices containing dynamical (in time) spatial coordinates. Optional argument, the default is NULL

distance

String; the name of the spatial distance. The default is Eucl, the euclidean distance. See the Section Details of GeoFit.

neighb

Numeric; an optional (vector of) positive integer indicating the order of neighborhood. See the Section Details for more information.

maxdist

A numeric value denoting the spatial distance Details.

maxtime

A numeric value denoting the temporal distance Details.

radius

Numeric; a value indicating the radius of the sphere when using the great circle distance. Default value is 1.

bivariate

Logical; if FALSE (the default) the data are interpreted as univariate spatial or spatial-temporal realisations. Otherwise they are intrepreted as a a realization from a bivariate field.

Author

Moreno Bevilacqua, ,https://sites.google.com/view/moreno-bevilacqua/home, Víctor Morales Oñate, , https://sites.google.com/site/moralesonatevictor/, Christian", Caamaño-Carrillo, ,https://www.researchgate.net/profile/Christian-Caamano

Details

The function returns the spatial or spatiotemporal indices of the pairs tha are neighboords of a certain order and/or with a certain fixed distance

Examples

require(GeoModels)
NN = 400
coords = cbind(runif(NN),runif(NN))
scale=0.5/3
param = list(mean=0,sill=1,nugget=0,scale=0.5/3,smooth=0.5)
corrmodel = "Matern"; 

param = list(mean=0,sill=1,nugget=0,scale=scale,smooth=0.5)
set.seed(951)
data = GeoSim(coordx = coords,corrmodel = corrmodel,
                  model = "Gaussian",param = param)$data

sel=GeoNeighIndex(coordx=coords,neighb=5)        

data1=data[sel$colidx]; data2=data[sel$rowidx]
## plotting pairs  that are neighboord of order 5
plot(data1,data2,xlab="",ylab="",main="h-scatterplot, neighb=5")