next up previous
Next: Field Dynamics Up: Multilinear Operations Previous: Convolution Connections

Convolution over Transformed Coordinates

In the definitions of correlation and convolution, Eqs. 1 and 2, the expressions s-r and r-s show us that these operations are sensitive to distance and direction in the domains of the fields, that is, they depend on the coordinates over which the fields are defined. For example, if results from by a coordinate transformation, , then the results of convolving with a Gaussian will not be the same as the results of convolving with . The convolution averages over regions that are circular in 's domain, whereas averages over circular regions in 's domain. For example, because of the logmap transformation between the retina and VI, a Gaussian convolution in VI will not have the effect of a Gaussian convolution in retinal coordinates or vice versa. This sensitivity of convolutions and correlations to the coordinate system can be a problem that needs to be solved or a computational resource that can be exploited.

Suppose we have two domains and such that fields over are transformations of fields over ; let be the coordinate transformation (an isomorphism). For example, and might be two brain regions (such as the retina and VI), or one or the other might be an external region (such as physical space around the body). Let and be two fields over and suppose we want to compute the convolution ; for example we might want to do a Gaussian convolution in retinal space. However, suppose that the convolution is to be computed by means of fields defined over the transformed domain . We are given the transformed and want to compute so that . We can get this by changing the integration variable of the convolution (assumed to be scalar to keep the example simple):

If we define the connectivity field [ A_uv = [h^-1(u)-h^-1(v)]h'[h^-1(v)] ,] then the convolution integral becomes [ _u = _' A_uv _v v, ] which is the integral operator, . This is a linear operator, but not a convolution, which means that it is still implemented by a simple pattern of connectivity, but that it is not a single pattern duplicated throughout the region. (If, as is often the case, the transformation h is a homeomorphism, then it will preserve the topology of , which means that a local convolution in will translate into local connections A in .)

We remark without proof that if the domains are of more than one dimension, then the connectivity kernel is defined [ A_uv = [h^-1(u)-h^-1(v)]; J[h^-1(v)] , ] where is the Jacobian of evaluated at v.

Now, conversely, suppose we do a convolution in the transformed coordinates; what is its effect in the original coordinates? By a similar derivation we find that where the kernel is defined [ C_xy = [h(x)-h(y)];J[h(y)] . ] In effect, the convolution kernel is projected backward through the transformation h. For example, if, like the logmap transformation, h expands the space in the center of the visual field and compresses it at the periphery, then the back-transformation of will result in a C that defines small receptive fields near the center of the visual field, and large ones near its periphery.


next up previous
Next: Field Dynamics Up: Multilinear Operations Previous: Convolution Connections

Bruce MacLennan
Wed Oct 2 16:55:07 EDT 1996