Invariant pattern recognition can be achieved by use of harmonic decomposition, for example circular harmonics are used for rotation invariant recognition. A common problem with such methods is that often only a single term of the harmonic decomposition is used, and it does not contain a sufficient amount of the reference energy. Thus discrimination capability is limited, especially in the presence of noise or other disturbances. By using several terms of the harmonic decomposition together this problem can be solved; this can be achieved by the use of code division filter multiplexing. Several harmonic terms are encoded onto a single filter, and the signal is simultaneously correlated with all of them, hence producing enhanced discrimination capabilities. Here two methods are suggested for such encoding. The first involves multiplexing the filters in the Fourier plane, while the second involves multiplexing in the image plane.