The VARIOGRAM Procedure

Pair Formation

The basic starting point in computing the empirical semivariance is the enumeration of pairs of points for the spatial data. Figure 19 shows the spatial domain D and the set of n measurements z Subscript i, i equals 1 comma ellipsis comma n, that have been sampled at the indicated locations in D. Two data points upper P 1 and upper P 2, with coordinates bold-italic s 1 equals left-parenthesis x 1 comma y 1 right-parenthesis and bold-italic s 2 equals left-parenthesis x 2 comma y 2 right-parenthesis, respectively, are selected for illustration.

A vector, or directed line segment, is drawn between these points. If the length

bar upper P Subscript i Baseline upper P Subscript j bar equals bar bold-italic s 2 minus bold-italic s 1 bar equals left-parenthesis x 2 minus x 1 right-parenthesis squared plus left-parenthesis y 2 minus y 1 right-parenthesis squared

of this vector is smaller than the specified DEPSILON= value, then the pair is excluded from the continuity measure calculations because the two points upper P 1 and upper P 2 are considered to be at zero distance apart (or collocated). Spatial collocation might appear due to different scales in sampling, observations made at the same spatial location at different time instances, and errors in the data sets. PROC VARIOGRAM excludes such pairs from the pairwise distance and semivariance computations because they can cause numeric problems in spatial analysis.

If this pair is not discarded on the basis of collocation, it is then classified—first by orientation of the directed line segment bold-italic s 2 minus bold-italic s 1, and then by its length bar upper P Subscript i Baseline upper P Subscript j bar. For example, it is unlikely for actual data that the distance bar upper P Subscript i Baseline upper P Subscript j bar between any pair of data points upper P Subscript i and upper P Subscript j located at bold-italic s Subscript i and bold-italic s Subscript j, respectively, would exactly satisfy bar upper P Subscript i Baseline upper P Subscript j Baseline bar equals bar bold-italic h bar equals h in the preceding computation of ModifyingAbove gamma With caret Subscript z Baseline left-parenthesis bold-italic h right-parenthesis. A similar argument can be made for the orientation of the segment bold-italic s 2 minus bold-italic s 1. Consequently, the pair upper P 1 upper P 2 is placed into an angle and distance class.

The following subsections give more details about the nature of these classifications. You can also find extensive discussions about the size and the number of classes to consider for the computation of the empirical semivariogram.

Figure 19: Selection of Points upper P 1 and upper P 2 in Spatial Domain D

 Selection of Points P1 and P2 in Spatial Domain


Last updated: December 09, 2022