{"id":1987,"date":"2014-02-05T17:32:41","date_gmt":"2014-02-05T17:32:41","guid":{"rendered":"https:\/\/www.anagram.at\/en\/diplomarbeit\/phase-difference\/"},"modified":"2014-02-05T17:32:41","modified_gmt":"2014-02-05T17:32:41","slug":"phase-difference","status":"publish","type":"page","link":"https:\/\/www.anagram.at\/en\/diplomarbeit\/phase-difference\/","title":{"rendered":"Phase Difference"},"content":{"rendered":"<p><body><br \/>\n<!--Navigation Panel--><br \/>\n<b> Next:<\/b> <a name=\"tex2html415\" href=\"https:\/\/www.anagram.at\/diplomarbeit\/feature-based-correspondence-analysis\/\">Feature-based Correspondence Analysis<\/a><br \/>\n<b> Up:<\/b> <a name=\"tex2html411\" href=\"https:\/\/www.anagram.at\/diplomarbeit\/intensity-based-correspondence-analysis\/\">Intensity-based correspondence analysis<\/a><br \/>\n<b> Previous:<\/b> <a name=\"tex2html407\" href=\"https:\/\/www.anagram.at\/diplomarbeit\/correlation\/\">Correlation<\/a><br \/>\n<!--End of Navigation Panel--><\/p>\n<h3><a name=\"SECTION00341200000000000000\"><br \/>\nPhase Difference<\/a><br \/>\n<\/h3>\n<p>\nAnother property to match is local frequency components [<a href=\"node47.html#Sang88\">San88<\/a>,<a href=\"node47.html#Fleet91\">FJJ91<\/a>]. If a function <img loading=\"lazy\" width=\"41\" height=\"37\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img145.png\" alt=\"$ f(x)$\"\/> with fourier transform <img loading=\"lazy\" width=\"44\" height=\"37\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img146.png\" alt=\"$ F(u)$\"\/> is shifted by an amount of <img loading=\"lazy\" width=\"14\" height=\"20\" align=\"BOTTOM\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img33.png\" alt=\"$ d$\"\/> then the resulting fourier transform of the shifted function <img loading=\"lazy\" width=\"74\" height=\"37\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img147.png\" alt=\"$ f(x+d)$\"\/> is <!-- MATH\n $e^{jdu}F(u)$\n --><br \/>\n<img loading=\"lazy\" width=\"75\" height=\"40\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img148.png\" alt=\"$ e^{jdu}F(u)$\"\/>. The shift in the spatial domain is equivalent to a phase shift in the frequency domain. It is possible to determine the disparity if the phase differences are found. Since the shift in the spatial domain is not equal for different regions of the image, for example the disparity differs for different objects that are mapped onto the image plane, a local frequency filter is needed to determine the phase differences. The <i>Gabor filter<\/i> [<a href=\"node47.html#Gabor46\">Gab46<\/a>], which is a bandpass filter with limited bandwidth, can be used for this. Equation <a href=\"#Gaborfilter\">2.27<\/a> shows the filter\n<\/p>\n<p\/>\n<div align=\"CENTER\"><a name=\"Gaborfilter\"\/><!-- MATH\n begin{equation}\ng_{w_0}(x) = gauss(x)e^{jw_0x}\nend{equation}\n --><\/p>\n<table cellpadding=\"0\" width=\"100%\" align=\"CENTER\">\n<tr valign=\"MIDDLE\">\n<td nowrap=\"nowrap\" align=\"CENTER\"><img loading=\"lazy\" width=\"193\" height=\"41\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img149.png\" alt=\"$displaystyle g_{w_0}(x) = gauss(x)e^{jw_0x}$\"\/><\/td>\n<td nowrap=\"nowrap\" width=\"10\" align=\"RIGHT\">\n(2.27)<\/td>\n<\/tr>\n<\/table>\n<\/div>\n<p><br clear=\"ALL\"\/><\/p>\n<p\/>\n<p\/>\n<div align=\"CENTER\"><!-- MATH\n begin{equation}\ngauss(x) = frac{1}{sigmasqrt{2pi}}e^{frac{-x^2}{2sigma^2}}\nend{equation}\n --><\/p>\n<table cellpadding=\"0\" width=\"100%\" align=\"CENTER\">\n<tr valign=\"MIDDLE\">\n<td nowrap=\"nowrap\" align=\"CENTER\"><img loading=\"lazy\" width=\"190\" height=\"59\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img150.png\" alt=\"$displaystyle gauss(x) = frac{1}{sigmasqrt{2pi}}e^{frac{-x^2}{2sigma^2}}$\"\/><\/td>\n<td nowrap=\"nowrap\" width=\"10\" align=\"RIGHT\">\n(2.28)<\/td>\n<\/tr>\n<\/table>\n<\/div>\n<p><br clear=\"ALL\"\/><\/p>\n<p\/>\nand Equation <a href=\"#gabor_fourier\">2.29<\/a> <\/p>\n<p\/>\n<div align=\"CENTER\"><a name=\"gabor_fourier\"\/><!-- MATH\n begin{equation}\nG_{w_0}(w) = e^{-frac{(w-w_0)^2}{2tau^2}}\nend{equation}\n --><\/p>\n<table cellpadding=\"0\" width=\"100%\" align=\"CENTER\">\n<tr valign=\"MIDDLE\">\n<td nowrap=\"nowrap\" align=\"CENTER\"><img loading=\"lazy\" width=\"163\" height=\"60\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img151.png\" alt=\"$displaystyle G_{w_0}(w) = e^{-frac{(w-w_0)^2}{2tau^2}}$\"\/><\/td>\n<td nowrap=\"nowrap\" width=\"10\" align=\"RIGHT\">\n(2.29)<\/td>\n<\/tr>\n<\/table>\n<\/div>\n<p><br clear=\"ALL\"\/><\/p>\n<p\/>\nshows its fourier transform. The first part of the <i>Gabor filter<\/i> is the <i>Gaussian function<\/i><\/p>\n<p>\n<img loading=\"lazy\" width=\"16\" height=\"19\" align=\"BOTTOM\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img152.png\" alt=\"$ sigma$\"\/> is the filter width and <img loading=\"lazy\" width=\"26\" height=\"33\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img153.png\" alt=\"$ w_0$\"\/> is the filter frequency for solving the correspondence problem. The product <!-- MATH\n $tausigma$\n --><br \/>\n<img loading=\"lazy\" width=\"26\" height=\"19\" align=\"BOTTOM\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img154.png\" alt=\"$ tausigma$\"\/> is one, which is the theoretical minimum of any linear complex filter [<a href=\"node47.html#Gabor46\">Gab46<\/a>]. Convolving <img loading=\"lazy\" width=\"31\" height=\"33\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img155.png\" alt=\"$ g_{w_0}$\"\/> with the image intensities <img loading=\"lazy\" width=\"42\" height=\"35\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img156.png\" alt=\"$ I_l,I_r$\"\/> yields a joint spatial and frequency representation of an image [<a href=\"node47.html#Daugman85\">Dau85<\/a>]:\n<\/p>\n<p\/>\n<div align=\"CENTER\"><!-- MATH\n begin{equation}\nbegin{aligned}\nc_l(x_0) = int I_l(x_0)g_{w_0}(x_0-x)dx = rho_l(x_0)e^{jphi_l(x_0)} \nc_r(x_0) = int I_r(x_0)g_{w_0}(x_0-x)dx = rho_r(x_0)e^{jphi_r(x_0)}\nend{aligned}\nend{equation}\n --><\/p>\n<table cellpadding=\"0\" width=\"100%\" align=\"CENTER\">\n<tr valign=\"MIDDLE\">\n<td nowrap=\"nowrap\" align=\"CENTER\"><img loading=\"lazy\" width=\"385\" height=\"108\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img157.png\" alt=\"begin{equation*}begin{aligned}c_l(x_0) = int I_l(x_0)g_{w_0}(x_0-x)dx = rho_...&#10;...(x_0)g_{w_0}(x_0-x)dx = rho_r(x_0)e^{jphi_r(x_0)} end{aligned}end{equation*}\"\/><\/td>\n<\/tr>\n<\/table>\n<\/div>\n<p><br clear=\"ALL\"\/><\/p>\n<p\/>\n<p>\nAs said before, a shift in the spatial domain is represented as a phase shift in the frequency domain, this gives an already estimation of the disparity <img loading=\"lazy\" width=\"14\" height=\"20\" align=\"BOTTOM\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img33.png\" alt=\"$ d$\"\/>.<\/p>\n<\/p>\n<p\/>\n<div align=\"CENTER\"><!-- MATH\n begin{equation}\nf(x+d) rightarrow F(u)e^{jdu}\nend{equation}\n --><\/p>\n<table cellpadding=\"0\" width=\"100%\" align=\"CENTER\">\n<tr valign=\"MIDDLE\">\n<td nowrap=\"nowrap\" align=\"CENTER\"><img loading=\"lazy\" width=\"174\" height=\"42\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img158.png\" alt=\"$displaystyle f(x+d) rightarrow F(u)e^{jdu}$\"\/><\/td>\n<td nowrap=\"nowrap\" width=\"10\" align=\"RIGHT\">\n(2.31)<\/td>\n<\/tr>\n<\/table>\n<\/div>\n<p><br clear=\"ALL\"\/><\/p>\n<p\/>\nThis theorem states that a spatial shift <img loading=\"lazy\" width=\"14\" height=\"20\" align=\"BOTTOM\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img33.png\" alt=\"$ d$\"\/> corresponds to a frequency shift <img loading=\"lazy\" width=\"25\" height=\"20\" align=\"BOTTOM\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img159.png\" alt=\"$ du$\"\/>. A suitable approximation of the local image shift is the normalized phase difference.<\/p>\n<p\/>\n<div align=\"CENTER\"><!-- MATH\n begin{equation}\nd(x) = frac{phi_r(x)-phi_l(x)}{w_0}\nend{equation}\n --><\/p>\n<table cellpadding=\"0\" width=\"100%\" align=\"CENTER\">\n<tr valign=\"MIDDLE\">\n<td nowrap=\"nowrap\" align=\"CENTER\"><img loading=\"lazy\" width=\"176\" height=\"63\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img160.png\" alt=\"$displaystyle d(x) = frac{phi_r(x)-phi_l(x)}{w_0}$\"\/><\/td>\n<td nowrap=\"nowrap\" width=\"10\" align=\"RIGHT\">\n(2.32)<\/td>\n<\/tr>\n<\/table>\n<\/div>\n<p><br clear=\"ALL\"\/><\/p>\n<p\/>\nThe estimation can be further improved using local frequencies instead of the filter mid-frequency <img loading=\"lazy\" width=\"26\" height=\"33\" align=\"MIDDLE\" border=\"0\" src=\"https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/https:\/\/www.anagram.at\/app\/uploads\/2014\/02\/img153.png\" alt=\"$ w_0$\"\/> [<a href=\"node47.html#Fleet91\">FJJ91<\/a>].<br \/>\nAdvantages of the phase difference method are that it is not that sensitive to noise and that the correspondence analysis works with subpixel accuracy. <\/p>\n<hr\/>\n<p><!--Navigation Panel--><b> Next:<\/b> <a name=\"tex2html415\" href=\"https:\/\/www.anagram.at\/diplomarbeit\/feature-based-correspondence-analysis\/\">Feature-based Correspondence Analysis<\/a><br \/>\n<b> Up:<\/b> <a name=\"tex2html411\" href=\"https:\/\/www.anagram.at\/diplomarbeit\/intensity-based-correspondence-analysis\/\">Intensity-based correspondence analysis<\/a><br \/>\n<!--End of Navigation Panel--><\/p>\n<p><\/body><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Phase Difference<\/p>\n","protected":false},"author":1,"featured_media":0,"parent":1946,"menu_order":0,"comment_status":"open","ping_status":"open","template":"","meta":{"_genesis_hide_title":false,"_genesis_hide_breadcrumbs":false,"_genesis_hide_singular_image":false,"_genesis_hide_footer_widgets":false,"_genesis_custom_body_class":"","_genesis_custom_post_class":"","_genesis_layout":""},"categories":[],"featured_image_src":null,"featured_image_src_square":null,"_links":{"self":[{"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/pages\/1987"}],"collection":[{"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/comments?post=1987"}],"version-history":[{"count":0,"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/pages\/1987\/revisions"}],"up":[{"embeddable":true,"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/pages\/1946"}],"wp:attachment":[{"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/media?parent=1987"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.anagram.at\/en\/wp-json\/wp\/v2\/categories?post=1987"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}