<html>
<head>
<style>
P
{
margin:0px;
padding:0px
}
body
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body>Jakub,<BR>
<BR>
There is one slight problem for when r = 0, this causes the crash.<BR>
The following fixes the crash and also gives a full circular area.<BR>
Note that I used (r, theta) and NOT (theta, r) as the output coordinate.<BR>
<BR>
in function <FONT size=2>TransformPoint:</FONT><BR>
<FONT size=2></FONT> <BR><FONT size=2><FONT size=2>
OutputPointType opoint(point);<BR>
opoint[0] = sqrt( point[0] * point[0] + point[1] * point[1] );<BR>
</FONT><FONT color=#0000ff size=2>if</FONT><FONT size=2>(point[0] != 0)<BR>
{<BR>
opoint[1] = acos( point[0] / opoint[0] );<BR>
}<BR>
</FONT><FONT color=#0000ff size=2>else<BR></FONT><FONT size=2>
{<BR>
opoint[1] = 0;<BR>
}<BR>
</FONT><FONT color=#0000ff size=2>if</FONT><FONT size=2>(point[0] < 0 || point[1] < 0 )<BR>
{<BR>
opoint[1] = 2 * vnl_math::pi - opoint[1];<BR>
}<BR>
</FONT><FONT color=#0000ff size=2>return</FONT><FONT face="" size=2> opoint;</FONT><BR>
<FONT face="" size=2></FONT> <BR>
<FONT face="" size=2>-----</FONT><BR>
<FONT face="" size=2></FONT> <BR><FONT face="" size=2><PRE>Hi.
I checked the code and the paper and i found some typos in the paper.
The code was good. You can compute the angle as acos(x/r), same as
atan(y/x) (better atan2(y,x)). These options are mathematically
equivalent. I modulate negative angle values, so that the angle is
always in [0,360]
I had typo in the paper so that there was acos(x/y) instead of
acos(x/r), but in the code it was correct.
By definition, polar transform is located in the center of cartesian
coordinate system. In the ITK transform layout, geometrical transforms
do not know anything about images so it is not possible to determine
the center of the image. The solution is (as you wrote) to set the
origin to -1*Size*Spacing*0.5 in each dimension. In this way, you can
also set the center to any aother point of the image. I am sure that
this will not crash as i am using the same code in my experiments...
Regards,
Jakub</PRE></FONT></FONT><br /><hr />Be one of the first to try <a href='http://ideas.live.com/programpage.aspx?versionId=5d21c51a-b161-4314-9b0e-4911fb2b2e6d' target='_new'>Windows Live Mail.</a></body>
</html>