Thread Subject: HowTo deal with this transposition in inverse problem ie. image restoration

Subject: HowTo deal with this transposition in inverse problem ie. image restoration

From: timedcy Asayaku

Date: 3 Apr, 2009 07:22:02

Message: 1 of 1

How to deal with this transposition in inverse problem ie. image restoration

For example, when minimizing || g – Hf ||^2 + λ||Qf||^2 , where image matrix f & g are the true and observed value in column-lexically heaped vectors, H is block-toeplitz matrix of the convolution kernel.

Then the solution is f = H^Tg/(H^TH+ λQ^TQ) , here H^T means the conjugate transpose matrix of H.

In programming, we usually compute it in the Fourier domain.
My question is how to deal with H^T in convolution ie H & x. conv2(H,x) (here H, x is a image matrix)?

And what is the relation between (H^T x) and (H x) in frequency domain?

ifftn( conj(fftn(H)) .* fftn(x) )?
or conv2(rot90(H,2), x)

this problem confused me very long long time!

Can you help me ?

Sincerely timedcy
20090403

Tags for this Thread

Everyone's Tags:

Add a New Tag:

Separated by commas
Ex.: root locus, bode

What are tags?

A tag is like a keyword or category label associated with each thread. Tags make it easier for you to find threads of interest.

Anyone can tag a thread. Tags are public and visible to everyone.

Tag Activity for This Thread
Tag Applied By Date/Time
convolution timedcy timedcy 3 Apr, 2009 03:25:07
fft timedcy timedcy 3 Apr, 2009 03:25:07
transposition timedcy timedcy 3 Apr, 2009 03:25:07
inverse problem timedcy timedcy 3 Apr, 2009 03:25:07
image restoration timedcy timedcy 3 Apr, 2009 03:25:07
deblur timedcy timedcy 3 Apr, 2009 03:25:07
regularization timedcy timedcy 3 Apr, 2009 03:25:07
fourier domain timedcy timedcy 3 Apr, 2009 03:25:07
rssFeed for this Thread

Contact us at files@mathworks.com