Path: news.mathworks.com!newsfeed-00.mathworks.com!newsfeed2.dallas1.level3.net!news.level3.com!postnews.google.com!w35g2000prg.googlegroups.com!not-for-mail
From: aheartthatlovesisalwaysyoung@gmail.com
Newsgroups: comp.soft-sys.matlab
Subject: How to deal with this transposition in inverse problem ie. image 
	restoration
Date: Thu, 2 Apr 2009 23:04:12 -0700 (PDT)
Organization: http://groups.google.com
Lines: 27
Message-ID: <471e1973-38cd-4ba9-b081-d9a74c7b8d9c@w35g2000prg.googlegroups.com>
NNTP-Posting-Host: 60.191.99.5
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-7
Content-Transfer-Encoding: quoted-printable
X-Trace: posting.google.com 1238738652 26484 127.0.0.1 (3 Apr 2009 06:04:12 GMT)
X-Complaints-To: groups-abuse@google.com
NNTP-Posting-Date: Fri, 3 Apr 2009 06:04:12 +0000 (UTC)
Complaints-To: groups-abuse@google.com
Injection-Info: w35g2000prg.googlegroups.com; posting-host=60.191.99.5; 
	posting-account=C4OUQwoAAAASp-XNuyTH5t7Fw4jpI80V
User-Agent: G2/1.0
X-HTTP-UserAgent: Mozilla/4.0 (compatible; MSIE 6.0; Windows NT 5.1; SV1; 
	Mozilla/4.0 (compatible; MSIE 6.0; Windows NT 5.1; SV1) ; 
	GreenBrowser),gzip(gfe),gzip(gfe)
Xref: news.mathworks.com comp.soft-sys.matlab:529942


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

For example, when minimizing || g - Hf ||^2 + =EB||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 =3D H^Tg/(H^TH+ =EBQ^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 AHeartThatLovesIsAlwaysYoung
20090403