SpatialVx (version 0.1-7)

saller: Analyze Features in a Field (Image)

Description

Functions to analyze features in a field (image) in conjunction, e.g., with the FeatureSuite function.

Usage

saller(object, x, y = NULL, matches = NULL, d = NULL, time.point = 1, 
    model = 1)
## S3 method for class 'saller':
print(x, ...)
## S3 method for class 'saller':
summary(object, ...)

Arguments

x

saller: x is a list object returned by convthresh or other feature identification function that returns a list with components X.feats, Y.feats (themselves lists with owin class objects defi

object

saller: A list object of class SpatialVx as returned by make.SpatialVx providing information about the original fields.

summary: object the returned by saller.

y

(optional) if not NULL, it is used instead of x, and otherwise has the same stucture thereof. When used with FeatureSuite, one might want to merge objects in one or both fields before calling the analysis function. If so,

matches

Not used by saller function. If matching of features between fields is required, this will be a two-column matrix describing which features from the forecast (first column) match to those of the verificaiton field (second column).

d

(optional) the SAL (saller) method requires division by the longest distance between two border points. If NULL, this is taken to be simply the length of the longest side.

time.point

numeric or character indicating which time point from the SpatialVx verification set to select for analysis.

model

numeric indicating which forecast model to select for the analysis.

...

not used.

Value

• saller returns a list with components:
• A
  numeric giving the amplitude component.
• L
  numeric giving the lcoation component.
• S
  numeric giving the structure component.
• print invisibly returns a named vector with S, A and L.

  summary does not return anything.

Details

saller: Computes S, A, and L of the SAL method introduced by Wernli et al. (2008).

References

Wernli, H., Paulat, M., Hagen, M. and Frei, C. (2008) SAL--A novel quality measure for the verification of quantitative precipitation forecasts. Mon. Wea. Rev., 136, 4470--4487, doi:10.1175/2008MWR2415.1.

See Also

FeatureSuite, centroid.owin, connected, tiles, tess, deltamm, make.SpatialVx

Examples

data(ExampleSpatialVxSet)

x <- ExampleSpatialVxSet$vx
xhat <- ExampleSpatialVxSet$fcst

q <- mean( c(c(x[x>0]),c(xhat[xhat>0])), na.rm=TRUE)

hold <- make.SpatialVx(x, xhat, field.type="contrived", units="none",
    data.name=c("Example", "x", "xhat"))

hold2 <- convthresh(hold, smoothpar=5, thresh=q)
plot(hold2, set.pw=TRUE)

look <- saller(hold, hold2)   
summary(look)

data(pert000)
data(pert004)
data(ICPg240Locs)

hold <- make.SpatialVx(pert000, pert004, loc=ICPg240Locs,
    projection=TRUE, map=TRUE,
    field.type="Precipitation", units="mm/h",
    data.name=c("Perturbed ICP Cases", "pert000", "pert004"))

look <- convthresh(hold, smoothpar=10.5)
summary(look)
plot(look, set.pw=TRUE)

saller(hold, look)