[Insight-users] Image registration - 2 different initializations

Radhika Sivaramakrishna radhika.sivaramakrishna at synarc.com
Thu, 15 Jan 2004 09:48:29 -0800


This message is in MIME format. Since your mail reader does not understand
this format, some or all of this message may not be legible.

------_=_NextPart_001_01C3DB8F.C9218CB2
Content-Type: text/plain

Hi Luis,
I appear to be having a strange problem with image registration. Basically I
am using the same code which you finally sent me although I have
combined the previous code you sent me to read a dicom series so as to
directly read a pair of test-retest DICOM series. Since the data are
unsigned short
there is already some problem I am having, I suspect w.r.t to Endianness
which I have not solved. But I am getting some other strange things which
I will try to describe below:
 
When I directly run the registration code on the 2 dicom series, then I get 
 
Transform center before Initializer 0  0  0
 Transform center after Initializer -2.34973  225.964  219.384
Intial Parameters =
[0, 0, 0, -0.0650073, 1.4008, -2.19123]
 Translationscale   = 1000
 
However, when I save the images as Raw format and then convert that to an
Analyze file using ImageJ and then load in the two images again, using the
straight version
of the code you sent me, here is what I get:
 
Transform center before Initializer 0  0  0
 Transform center after Initializer -2.34973  225.964  219.384
Intial Parameters =
[0, 0, 0, -0.0650073, 1.4008, -2.19123]
 Translationscale   = 1000
 
I cant figure out what is going on. It appears the data is being read fine,
because the raw files for fixed and moving image are as they should be when
I view it on ImageJ.
 
Can you help figure out what is happening. 
 
Thanks
Radhika
 
                                                     
-----------------------------------------------------  
Confidentiality Notice. 
This email message is for the sole use of the intended recipient(s) and may
contain confidential and privileged information. Any unauthorized review,
use, disclosure or distribution is prohibited. If you are not the intended
recipient, please contact the sender by reply email and destroy all copies
of the original message. If you are the intended recipient, please be
advised that the content of this message is subject to access, review and
disclosure by the sender's Email System Administrator.

------_=_NextPart_001_01C3DB8F.C9218CB2
Content-Type: text/html
Content-Transfer-Encoding: quoted-printable

<html xmlns:o=3D"urn:schemas-microsoft-com:office:office" =
xmlns:w=3D"urn:schemas-microsoft-com:office:word" =
xmlns=3D"http://www.w3.org/TR/REC-html40">

<head>
<META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
charset=3Dus-ascii">


<meta name=3DProgId content=3DWord.Document>
<meta name=3DGenerator content=3D"Microsoft Word 10">
<meta name=3DOriginator content=3D"Microsoft Word 10">
<link rel=3DFile-List href=3D"cid:filelist.xml at 01C3DB4C.BAF4AAE0">
<!--[if gte mso 9]><xml>
 <o:OfficeDocumentSettings>
  <o:DoNotRelyOnCSS/>
 </o:OfficeDocumentSettings>
</xml><![endif]--><!--[if gte mso 9]><xml>
 <w:WordDocument>
  <w:SpellingState>Clean</w:SpellingState>
  <w:GrammarState>Clean</w:GrammarState>
  <w:DocumentKind>DocumentEmail</w:DocumentKind>
  <w:EnvelopeVis/>
  <w:Compatibility>
   <w:BreakWrappedTables/>
   <w:SnapToGridInCell/>
   <w:WrapTextWithPunct/>
   <w:UseAsianBreakRules/>
  </w:Compatibility>
  <w:BrowserLevel>MicrosoftInternetExplorer4</w:BrowserLevel>
 </w:WordDocument>
</xml><![endif]-->
<style>
<!--
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
	{mso-style-parent:"";
	margin:0in;
	margin-bottom:.0001pt;
	mso-pagination:widow-orphan;
	font-size:12.0pt;
	font-family:"Times New Roman";
	mso-fareast-font-family:"Times New Roman";}
a:link, span.MsoHyperlink
	{color:blue;
	text-decoration:underline;
	text-underline:single;}
a:visited, span.MsoHyperlinkFollowed
	{color:purple;
	text-decoration:underline;
	text-underline:single;}
span.EmailStyle17
	{mso-style-type:personal-compose;
	mso-style-noshow:yes;
	mso-ansi-font-size:10.0pt;
	mso-bidi-font-size:10.0pt;
	font-family:Arial;
	mso-ascii-font-family:Arial;
	mso-hansi-font-family:Arial;
	mso-bidi-font-family:Arial;
	color:windowtext;}
span.SpellE
	{mso-style-name:"";
	mso-spl-e:yes;}
span.GramE
	{mso-style-name:"";
	mso-gram-e:yes;}
 at page Section1
	{size:8.5in 11.0in;
	margin:1.0in 1.25in 1.0in 1.25in;
	mso-header-margin:.5in;
	mso-footer-margin:.5in;
	mso-paper-source:0;}
div.Section1
	{page:Section1;}
-->
</style>
<!--[if gte mso 10]>
<style>
 /* Style Definitions */=20
 table.MsoNormalTable
	{mso-style-name:"Table Normal";
	mso-tstyle-rowband-size:0;
	mso-tstyle-colband-size:0;
	mso-style-noshow:yes;
	mso-style-parent:"";
	mso-padding-alt:0in 5.4pt 0in 5.4pt;
	mso-para-margin:0in;
	mso-para-margin-bottom:.0001pt;
	mso-pagination:widow-orphan;
	font-size:10.0pt;
	font-family:"Times New Roman";}
</style>
<![endif]-->
</head>

<body lang=3DEN-US link=3Dblue vlink=3Dpurple =
style=3D'tab-interval:.5in'>

<div class=3DSection1>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Hi Luis,<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>I appear to be having a strange problem with image
registration. Basically I am using the same code which you finally sent =
me
although I have<o:p></o:p></span></font></p>

<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>combined</span></font></spa=
n><font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> the
previous code you sent me to read a <span class=3DSpellE>dicom</span> =
series so
as to directly read a pair of test-retest DICOM series. Since the data =
are
unsigned short<o:p></o:p></span></font></p>

<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>there</span></font></span><=
font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> is already
some problem I am having, I suspect <span class=3DSpellE>w.r.t</span> =
to <span
class=3DSpellE>Endianness</span> which I have not solved. But I am =
getting some
other strange things which<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>I will try to describe =
below:<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>When I directly run the registration code on the 2 =
<span
class=3DSpellE>dicom</span> series, then I get =
<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Transform center before <span =
class=3DSpellE>Initializer</span>
<span class=3DGramE>0<span style=3D'mso-spacerun:yes'>&nbsp; =
</span><span
class=3DSpellE>0</span></span><span style=3D'mso-spacerun:yes'>&nbsp; =
</span><span
class=3DSpellE>0</span><o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><span =
style=3D'mso-spacerun:yes'>&nbsp;</span>Transform center
after <span class=3DSpellE>Initializer</span> -<span =
class=3DGramE>2.34973<span
style=3D'mso-spacerun:yes'>&nbsp; </span>225.964</span><span
style=3D'mso-spacerun:yes'>&nbsp; =
</span>219.384<o:p></o:p></span></font></p>

<p class=3DMsoNormal><span class=3DSpellE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>Intial</span></font></span>=
<font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> Parameters
=3D<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>[0, 0, 0, -0.0650073, 1.4008, =
-2.19123]<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><span style=3D'mso-spacerun:yes'>&nbsp;</span><span
class=3DSpellE>Translationscale</span><span =
style=3D'mso-spacerun:yes'>&nbsp;&nbsp;
</span>=3D 1000<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>However, when I save the images as <span =
class=3DGramE>Raw</span>
format and then convert that to an Analyze file using <span =
class=3DSpellE>ImageJ</span>
and then load in the two images again, using the straight =
version<o:p></o:p></span></font></p>

<p class=3DMsoNormal><span class=3DGramE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>of</span></font></span><fon=
t size=3D2
face=3DArial><span style=3D'font-size:10.0pt;font-family:Arial'> the =
code you sent
me, here is what I get:<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Transform center before <span =
class=3DSpellE>Initializer</span>
<span class=3DGramE>0<span style=3D'mso-spacerun:yes'>&nbsp; =
</span><span
class=3DSpellE>0</span></span><span style=3D'mso-spacerun:yes'>&nbsp; =
</span><span
class=3DSpellE>0</span><o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><span =
style=3D'mso-spacerun:yes'>&nbsp;</span>Transform center
after <span class=3DSpellE>Initializer</span> -<span =
class=3DGramE>2.34973<span
style=3D'mso-spacerun:yes'>&nbsp; </span>225.964</span><span
style=3D'mso-spacerun:yes'>&nbsp; =
</span>219.384<o:p></o:p></span></font></p>

<p class=3DMsoNormal><span class=3DSpellE><font size=3D2 =
face=3DArial><span
style=3D'font-size:10.0pt;font-family:Arial'>Intial</span></font></span>=
<font
size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;font-family:Arial'> Parameters
=3D<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>[0, 0, 0, -0.0650073, 1.4008, =
-2.19123]<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><span style=3D'mso-spacerun:yes'>&nbsp;</span><span
class=3DSpellE>Translationscale</span><span =
style=3D'mso-spacerun:yes'>&nbsp;&nbsp;
</span>=3D 1000<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>I <span class=3DGramE>cant</span> figure out what is =
going on.
It appears the data is being read fine, because the raw files for fixed =
and
moving image are as they should be when I view it on <span =
class=3DSpellE>ImageJ</span>.<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Can you help figure out what is happening. =
<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Thanks<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'>Radhika<o:p></o:p></span></font></p>

<p class=3DMsoNormal><font size=3D2 face=3DArial><span =
style=3D'font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

</div>

</body>

</html>

<P><FONT SIZE=3D2 FACE=3D"Arial">                                       =
              </FONT></P>

<P><FONT SIZE=3D2 =
FACE=3D"Arial">-----------------------------------------------------  =
</FONT></P>

<P><FONT SIZE=3D2 FACE=3D"Arial">Confidentiality Notice. </FONT></P>

<P><FONT SIZE=3D2 FACE=3D"Arial">This email message is for the sole use =
of the intended recipient(s) and may contain confidential and =
privileged information. Any unauthorized review, use, disclosure or =
distribution is prohibited. If you are not the intended recipient, =
please contact the sender by reply email and destroy all copies of the =
original message. If you are the intended recipient, please be advised =
that the content of this message is subject to access, review and =
disclosure by the sender's Email System Administrator.</FONT></P>

------_=_NextPart_001_01C3DB8F.C9218CB2--