• 12
• 12
• 9
• 10
• 13

# Decomposition of a 2D spectrum

This topic is 3082 days old which is more than the 365 day threshold we allow for new replies. Please post a new topic.

## Recommended Posts

I have a problem for those with experience in 2D signal processing. Suppose I have a large 2D hermitian spectrum representing a 2D real signal (S). The straightforward way to calculate the signal S is by applying an IFFT to the entire spectrum. What I am trying to do instead is to subdivide the spectrum into non-overlapping bands, calculate the corresponding signals (Si) again with an IFFT, and then build S by combining Si. A schematic representation of the spectrum, subdivided into numbered bands, is (excuse my bad ASCII art):
^
|
|-----------|
| 4   |  3  |
|_____|_____|
|     |     |
|  1  |  2  |
|--------------------->
|     |     |
|  5  |  6  |
|-----------|
|  7  |  8  |
|-----------|


The left side of the spectrum is the conjugate of the right side and is ignored. Let's start with the bands (1,5). We can calculate the signal (S1+S5) by zeroing all other bands and calculating an IFFT of the spectrum (1,5), with half the size as the full one. This is the simple part. The problem comes in applying the same approach to the remaining couples of high-frequency bands - for example { (4,7), (2,6), (3,8) } . First, they have to be shifted to the low-frequencies (previously occupied by (1,5)) in order to use a half-sized IFFT. Second, the corresponding signal has to be brought back to the correct frequencies before being combined with the others to reconstruct the original signal S. But how?? I can't figure out the correct way to do this, considering the implicit periodic nature of the spectrum and signals, the hermitian requirement, etc. I am 100% positive this problem is shown and solved, at least in the 1D case, in my books on digital processing but unfortunately they're in other country. Can anybody help me in this or maybe point me to another forum? Thanks! Stefano Lanza

##### Share on other sites
I am not really an expert and I don't have the time to work out the details in your case, but it looks like the Wikipedia page on the discrete Fourier transform has a theorem that might help.