SpatialVx (version 0.1-5)

FeatureMatchAnalyzer: Analyze Features of a Verification Set

Description

Analyze matched features of a verification set.

Usage

FeatureMatchAnalyzer(x, y = NULL, matches = NULL, object = NULL,
                     which.comps = c("cent.dist", "angle.diff", "area.ratio", "int.area",
                     "bdelta", "haus", "ph", "mhd", "med", "msd", "fom", "minsep", "bearing"),
                     sizefac = 1, alpha = 0.1, k = 4, p = 2, c = Inf, distfun = "distmapfun", ...)
## S3 method for class 'FeatureMatchAnalyzer':
summary(object, ...)
## S3 method for class 'FeatureMatchAnalyzer':
plot(x, ...)
FeatureComps(Y, X, which.comps=c("cent.dist", "angle.diff", "area.ratio", "int.area",
             "bdelta", "haus", "ph", "mhd", "med", "msd", "fom", "minsep", "bearing"), 
             sizefac=1, alpha=0.1, k=4, p=2, c=Inf, distfun="distmapfun", deg = TRUE,
             aty = "compass", ...)

Arguments

x,y,matches

x, y and matches are list objects with components as output by deltamm or similar function. Only one is used, and it first checks for matches, then y, and finally x

X,Y

list object giving a pixel image as output from solutionset from package spatstat for the verification and forecast fields, resp. These arguments are passed directly to the locperf function.

object

Obligatory argument for use with FeatureSuite function. Not used by FeatureMatchAnalyzer (so far). For the summary method function for FeatureMatchAnalyzer, this is the returned value from the self-sam

which.comps

character vector indicating which properties of the features are to be analyzed.

sizefac

single numeric by which area calculations should be multiplied in order to get the desired units. If unity (default) results are in terms of grid squares.

alpha

numeric value for the FOM measure (see the help file for locperf.

k

numeric indicating which quantile to use if the partial Hausdorff measure is to be used.

p

numeric giving the value of the parameter p for the Baddeley metric.

c

numeric giving the cut-off value for the Baddeley metric.

distfun

character naming a distance functions to use in calculating the various binary image measures. Default is Euclidean distance.

deg, aty

optional arguments to the bearing function.

...

Additional arguments to deltametric from package spatstat. In case of the summary method function, additional optional arguments may be passed, which include silent (logical, should the information be pri

Value

• FeatureMatchAnalyzer returns a list of list objects. The specific components depend on the 'which.comps' argument, and are the same as those returned by FeatureComps. These can be any of the following.
• cent.dist
  numeric giving the centroid (Euclidean) distance.
• angle.diff
  numeric giving the orientation (major axis) angle difference.
• area.ratio
  numeric giving the area ratio, which is always between 0 and 1 because this is defined by Davis et al. (2006) to be the area of the smaller feature divided by that of the larger feature regardless of which field the feature belongs to.
• int.area
  numeric giving the intersection area of the features.
• bdelta
  numeric giving Baddeley's delta metric between the two features.
• haus, ph, mhd, med, msd, fom, minsep
  numeric, see locperf for specific information.
• bearing
  numeric giving the bearing from the forecast object centroid to the observed object centroid.
• The summary method for FeatureMatchAnalyzer invisibly returns a matrix with the same information, but where each matched object is a row and each column is the specific statistic. Or, if optional interest argument is passed, a list with components:
• match.properties
  above mentioned matrix.
• object.interest
  numeric vector giving the total interest for each feature.
• interest.values
  numeric vector giving the interest values used (i.e., as passed through argument interest.

Details

FeatureMatchAnalyzer is designed to be used with FeatureSuite. It is set up to calculate the values discussed in sec. 4 of Davis et al. (2006) for a single verification set (i.e., mean and standard deviation are not computed because it is only a single case). If criteria is 1, then features separated by a distance D < the sum of the sizes of the two features (size of a feature is defined as the square root of its area) are considered a match. If criteria is 2, then a match is made if D < the average of the sizes of the two features. Finally, criteria 3 decides a match as being anything less than a pre-determined constant.

FeatureComps is the primary function called by FeatureMatchAnalyzer, and is designed as a more stand-alone type of function. Several of the measures that can be calculated are simply the binary image measures/metrics available via, e.g., locperf. It calculates comparisons between two matched features (i.e., between the verification and forecast fields).

The summary method function for FeatureMatchAnalyzer allows for passing a function, con, to determine confidence for each interest value. The idea being to set the interest to zero when the particular interest value does not make sense. For example, angle difference makes no sense if both objects are circles. Currently, no functions are included in this package for actually doing this, and so the functionality itself has not been tested.

References

Davis, C. A., Brown, B. G. and Bullock, R. G. (2006) Object-based verification of precipitation forecasts, Part I: Methodology and application to mesoscale rain areas. Mon. Wea. Rev., 134, 1772--1784.

See Also

locperf, FeatureSuite, convthresh, deltamm, deltametric, bearing

Examples

data(ExampleSpatialVxSet)

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

hold <- FeatureSuitePrep("x", "xhat")
look <- convthresh( hold, smoothpar=1.5)
look2 <- centmatch(look, object=hold)
tmp <- FeatureMatchAnalyzer(matches=look2)
summary(tmp)
plot(tmp)

data(pert000)
data(pert004)
hold <- FeatureSuitePrep("pert000", "pert004")
look <- convthresh( hold, smoothpar=10.5)
look2 <- centmatch(look, object=hold)
tmp <- FeatureMatchAnalyzer(matches=look2)
summary(tmp)
plot(tmp)