(2) |
rabdomain
Figure 4 Domain of inputs and outputs of a two-dimensional prediction-error filter. |
rabtent
Figure 5 Placement of tent-like weighting function in the space of filter inputs and outputs. |
We need a weighting function that vanishes where the filter has no outputs. The amplitude of the weighting function is not very important because we have learned how to put signals back together properly for arbitrary weighting functions. We can use any pyramidal or tent-like shape that drops to zero outside the domain of the filter output. The job is done by subroutine tentn(). A new parameter needed by tentn is a, the coordinate of the beginning of the tent. tenttent weights
In applications where triangle weights are needed on the inputs (or where we can work on a patch without having interference with edges), we can get ``triangle tent'' weights from tentn() if we set filter dimensions and lags to unity, as shown in Figure .
windwt90
Figure 6 Window weights from tentn() with nwind=(61,19), center=(31,1), a=(1,1) . |
Triangle weighting functions triangle weighting functions can sum to a constant if the spacing is such that the midpoint of one triangle is at the beginning of the next. I imagined in two dimensions that something similar would happen with shapes like Egyptian pyramids of Cheops, 2-|x-y|+|x+y|. Instead, the equation (1-|x|)(1-|y|) which has the tent-like shape shown in Figure adds up to the constant flat top shown in Figure . (To add interest to Figure , I separated the windows by a little more than the precise matching distance.) In practice we may chose window shapes and overlaps for reasons other than the constancy of the sum of weights, because mkwallwt accounts for that.
wallwt90
Figure 7 (Inverse) wall weights with n1=100, w1=61, k1=2, n2=30, w2=19, k2=2 |
Finally is an example of filtering a plane of uniform constants with an impulse function. The impulse function is surrounded by zeros, so the filter output patches are smaller than the input patches back in Figure . Here in Figure , both axes need more window density.
cinloip90
Figure 8 Filtering in patches with the same parameters as in Figures and . Additionally, the filter parameters are a1=11 a2=5 lag1=6 lag2=1 . Thus, windows are centered on the 1-axis and pushed back out the 2-axis. |