<br>Hi Michael,<br><br>Thanks for running these additional tests <br>and for reporting back.<br><br><br>0) You mentioned that you found an error in the Software Guide<br> in page 28.<br><br> Could you please let us know what it was, so we can fix it ?<br>
(Thanks).<br><br><br>1) I have a less optimistic interpretation of the results from your test.<br> If I understand correctly from your email,<br> You inserted the two writers for saving the input images <br> before calling the resampling, and you say that this code<br>
crashes in Release mode, while it works on Debug.<br><br> I would say that this clearly indicates that something is <br> going wrong *before* the resampling.<br><br> I common cause for code to work in Debug and not in Release,<br>
is to have uninitialized variables. You may want to run your code<br> through a dynamic checker such ad Valgrind or Purify.<br><div class="gmail_quote"><br><br>2) I'm not aware of any incompatibilities between ITK and wxWidgets,<br>
and would be very surprised to learn that there is any. <br> GUI libraries should be totally independent.<br><br><br>3) Typical suspects at this point are:<br><br> - Unititialized variables<br><br> - ITK Classes that have been created in a certain scope<br>
and that have been destroyed by the time you get them<br> to the code that performs the resampling. (e.g. check the<br> pointer to Transform). You may have SmartPointer that<br> are pointing to classes that have been already deallocated.<br>
<br> One way of debugging this is to connect an Observer to <br> the suspect objects. Upon destruction ITK classes invoke<br> an itk::DeleteEvent.<br><br><br>4) BTW: Why are you using GetPointer() in the calls:<br>
<br> resampler->SetTransform( transform.GetPointer() );<br> resampler->SetInterpolator( interpolator.GetPointer() );<br><br> That shouldn't be necessary.<br> It you did this in order to get around some compilation error,<br>
it is likely that that compilation error was flagging a fundamental<br> problem with the code.<br><br> Could you please remove the GetPointer() calls ?<br> It should just be:<br><br> resampler->SetTransform( transform );<br>
resampler->SetInterpolator( interpolator );<br>
<br><br>5) And... another big suspect....<br><br> Looking back at your first email, you told us that the crash <br> happens in the execution of the ProcessObject update at:<br><br><br> if ( m_Source ) // <---m_Source is valid, <br>
but when i step over<br>
here, it gets an invalid pointer<br><br> This is probably because when you generate the images<br> with the import filter, you are not disconnecting them from<br> the filter, so they get to the Resampling function believing<br>
that they are still attached to a source filter, when in fact,<br> that import filter was destroyed when reached the end <br> of the scope in the function <br><br> dataInterface2itk()<br>
<br> Please, at the end of that function, replace the code<br><br> import->Update();<br>
image = import->GetOutput();<br> image->Update();<br>
<br> // cout << image;<br> return image;<br> }<br><br> with<br><br> import->Update();<br>
image = import->GetOutput();<br>
image->Update();<br>
<br>
image->DisconnectPipeline(); // IMPORTANT !<br><br>
return image;<br>
}<br>
<br><br><br>I'm suspecting that (5) is the real source of the problem.<br><br><br> Please try this and let us know what you find.<br><br><br> Thanks<br><br><br> Luis<br><br><br>----------------------------------------------------------------<br>
On Thu, Jul 2, 2009 at 2:04 PM, Michael Schildt <span dir="ltr"><<a href="mailto:michael.schildt@ifn-magdeburg.de">michael.schildt@ifn-magdeburg.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello Luis,<br>
<br>
I checked the things, you mentioned:<br>
<br>
<br>
Luis Ibanez schrieb:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Hi Michael,<br>
<br><div class="im">
Thanks for adding the Print() statements<br>
and posting the outcome to the mailing list.<br>
<br>
The image data look ok.<br>
and so it does the filter metadata.<br>
<br>
It is good that you also added the code<br>
that you are using for generating the images.<br>
<br>
A common suspect when you using the<br>
import filter is to get wrong the ownership<br>
of the image pixel buffer.<br>
<br>
You are using<br>
<br>
const unsigned int numberOfPixels = diImage->getSize();<br>
dType* data = diImage->getData();<br>
import->SetImportPointer( data, numberOfPixels, false );<br>
<br>
I'll strongly suggest that you verify that the<br>
"numberOfPixels" is correct here, and also<br>
reconsider whether the flag in SetImportPointer()<br>
should be "false" or "true". <br>
</div></blockquote>
I did this intentionally after consulting documentation to save memory. I have the data sets already in memory and they are displayed correctly in the views of my wxWidgets application. So, i just borrow this memory to ITK to do the registration task, whch by the way give strange results even for registration of one replicated dataset. Afterwards datasets are further used in the program without assitance of ITK. For that reason i do not want ITK to own and release the memory. The number of byte i have verified it is 256x256x256=16777216 so it is correct too.<div class="im">
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
When set to true, the import filter will own the memory<br>
buffer, which means that it will deallocate it when the<br>
importer is destroyed.<br>
<br>
<br>
The process of importing the image is a lot more<br>
suspicious than the resampling. <br>
Please do the following:<br>
<br>
1) Before the resampling code, instantiate two<br>
ImageFileWriters, one for the fixed image,<br>
and the other one for the moving image.<br>
and write the images to a fileformat for which<br>
you have a trusty image viewer.<br>
<br>
(you could use Slicer, VV, ParaView, SNAP...)<br>
</blockquote></div>
I have done this and included write just before resample stuff posted earlier in this thread:<br>
<br>
itk::ImageFileWriter<FloatImage3DType>::Pointer writerfixedImage;<br>
writerfixedImage = itk::ImageFileWriter<FloatImage3DType>::New();<br>
writerfixedImage->SetInput (fixedIm);<br>
writerfixedImage->SetFileName( "c:\\temp\\fixedImage.hdr" );<br>
writerfixedImage->Update();<br>
<br>
itk::ImageFileWriter<FloatImage3DType>::Pointer writermovingImage;<br>
writermovingImage = itk::ImageFileWriter<FloatImage3DType>::New();<br>
writermovingImage->SetInput (movingIm);<br>
writermovingImage->SetFileName( "c:\\temp\\movingImage.hdr" );<br>
writermovingImage->Update(); <br>
By the way, i took the exampe from the ItkSoftwareGuide-2.4.0.pdf on page 28 but it contains several error. Maybe it can be noted for a update of the guide. In Debug mode i get out the data as Analyze and it looks exactly like the original data sets. In Release mode the application crashes already at writerfixedImage->Update(); So the good point is, it does seem to be the resampler.<br>
But it gets even more confusing. I wrote a short test programm using my code to load the data, resample and save via the source above. It runs in Debug and Release mode without crash! So, i copied this small piece into my application to be executed on button press. That crashs again. I copied the piece in the MyApp::OnInit() of my wxWidgets application, it crashes too. Then, I change the application type of the test programm from console to windows program and it still runs fine in Release and Debug. So the only difference seems to be that my application is using wxWidgets, but the test programm does not. In the debugger it looks like the ITK is executed on NULL-Pointer but this makes no sense and could be due to compiler optimizations. At that point i have no idea what to test further.<br>
<br>
Are there any issues or side effects known when using ITK together with wxWidgets?<br>
<br>
With best reguards,<br>
Michael Schildt<div class="im"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
2) Run the code and save the images, no need<br>
to run the resampler.<br>
</blockquote></div>
I did so.<div class="im"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
3) Open the images with the viewer and let us<br>
know if they are in good shape.<br>
<br>
4) I would expect that the images are corrupted<br>
at this point, but we will only know for sure<br>
after you run this test.<br>
</blockquote></div>
At least they look exactly the way they should do, so this doesn't seem to be the point.<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
<br>
<br>
Please let us now what you find.<br>
<br>
<br>
Thanks<br>
<br>
<br>
Luis<br>
<br>
<br>
<br>
----------------------------------------------------------<br></div><div><div></div><div class="h5">
On Wed, Jul 1, 2009 at 8:28 AM, Michael Schildt <<a href="mailto:michael.schildt@ifn-magdeburg.de" target="_blank">michael.schildt@ifn-magdeburg.de</a> <mailto:<a href="mailto:michael.schildt@ifn-magdeburg.de" target="_blank">michael.schildt@ifn-magdeburg.de</a>>> wrote:<br>
<br>
Hello Luis,<br>
<br>
thank you for your suggestions.<br>
<br>
I inserted the Prints in my source code and here ist the output:<br>
--------------------------- fixedIm -------------------------------<br>
Image (00339668)<br>
RTTI typeinfo: class itk::Image<float,3><br>
Reference Count: 3<br>
Modified Time: 20<br>
Debug: Off<br>
Observers:<br>
none<br>
Source: (none)<br>
Source output index: 0<br>
Release Data: Off<br>
Data Released: False<br>
Global Release Data: Off<br>
PipelineMTime: 13<br>
UpdateMTime: 19<br>
LargestPossibleRegion:<br>
Dimension: 3<br>
Index: [0, 0, 0]<br>
Size: [256, 256, 256]<br>
BufferedRegion:<br>
Dimension: 3<br>
Index: [0, 0, 0]<br>
Size: [256, 256, 256]<br>
RequestedRegion:<br>
Dimension: 3<br>
Index: [0, 0, 0]<br>
Size: [256, 256, 256]<br>
Spacing: [1, 1, 1]<br>
Origin: [0, 0, 0]<br>
Direction:<br>
1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
IndexToPointMatrix:<br>
1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
PointToIndexMatrix:<br>
1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
PixelContainer:<br>
ImportImageContainer (01390EA8)<br>
RTTI typeinfo: class itk::ImportImageContainer<unsigned<br>
long,float><br>
Reference Count: 1<br>
Modified Time: 17<br>
Debug: Off<br>
Observers:<br>
none<br>
Pointer: 034B0020<br>
Container manages memory: false<br>
Size: 16777216<br>
Capacity: 16777216<br>
<br>
--------------------------- movingIm ------------------------------<br>
Image (00337178)<br>
RTTI typeinfo: class itk::Image<float,3><br>
Reference Count: 4<br>
Modified Time: 40<br>
Debug: Off<br>
Observers:<br>
none<br>
Source: (none)<br>
Source output index: 0<br>
Release Data: Off<br>
Data Released: False<br>
Global Release Data: Off<br>
PipelineMTime: 33<br>
UpdateMTime: 39<br>
LargestPossibleRegion:<br>
Dimension: 3<br>
Index: [0, 0, 0]<br>
Size: [256, 256, 256]<br>
BufferedRegion:<br>
Dimension: 3<br>
Index: [0, 0, 0]<br>
Size: [256, 256, 256]<br>
RequestedRegion:<br>
Dimension: 3<br>
Index: [0, 0, 0]<br>
Size: [256, 256, 256]<br>
Spacing: [1, 1, 1]<br>
Origin: [0, 0, 0]<br>
Direction:<br>
1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
IndexToPointMatrix:<br>
1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
PointToIndexMatrix:<br>
1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
PixelContainer:<br>
ImportImageContainer (0033BE98)<br>
RTTI typeinfo: class itk::ImportImageContainer<unsigned<br>
long,float><br>
Reference Count: 1<br>
Modified Time: 37<br>
Debug: Off<br>
Observers:<br>
none<br>
Pointer: 074C0020<br>
Container manages memory: false<br>
Size: 16777216<br>
Capacity: 16777216<br>
<br>
--------------------------- resampler -----------------------------<br>
ResampleImageFilter (00338A08)<br>
RTTI typeinfo: class itk::ResampleImageFilter<class<br>
itk::Image<float,3>,cla<br>
s itk::Image<float,3>,double><br>
Reference Count: 2<br>
Modified Time: 64<br>
Debug: Off<br>
Observers:<br>
none<br>
Number Of Required Inputs: 1<br>
Number Of Required Outputs: 1<br>
Number Of Threads: 1<br>
ReleaseDataFlag: Off<br>
ReleaseDataBeforeUpdateFlag: Off<br>
Input 0: (00337178)<br>
Output 0: (00338B88)<br>
AbortGenerateData: Off<br>
Progress: 0<br>
Multithreader:<br>
RTTI typeinfo: class itk::MultiThreader<br>
Reference Count: 1<br>
Modified Time: 49<br>
Debug: Off<br>
Observers:<br>
none<br>
Thread Count: 1<br>
Global Maximum Number Of Threads: 128<br>
Global Default Number Of Threads: 1<br>
DefaultPixelValue: 0<br>
Size: [256, 256, 256]<br>
OutputStartIndex: [0, 0, 0]<br>
OutputOrigin: [0, 0, 0]<br>
OutputSpacing: [1, 1, 1]<br>
OutputDirection: 1 0 0<br>
0 1 0<br>
0 0 1<br>
<br>
Transform: 00337350<br>
Interpolator: 00338860<br>
UseReferenceImage: Off<br>
<br>
I already call Update on the Images but after the Set..- commands. I'm<br>
not sure about the order of commands to import the images to itk. Here<br>
is my import method called to create fixedIm and movingIm, maybe this<br>
helps to localize the problem (MRIScalarLayer is a custom class just<br>
containing the image data and some infos):<br>
<br>
template<typename dType> typename itk::Image<dType,3>::Pointer<br>
dataInterface2itk(MRIScalarLayer<dType> *diImage)<br>
{<br>
typedef itk::Image<dType,3> Image3DType;<br>
typename Image3DType::Pointer image = Image3DType::New();<br>
typedef itk::ImportImageFilter<dType,3> ImportFilterImage3DType;<br>
typename ImportFilterImage3DType::Pointer import =<br>
ImportFilterImage3DType::New();<br>
typename ImportFilterImage3DType::SizeType size;<br>
size[0] = diImage->getSize((unsigned char)0); // size along X<br>
size[1] = diImage->getSize((unsigned char)1); // size along Y<br>
size[2] = diImage->getSize((unsigned char)2); // size along Z<br>
typename ImportFilterImage3DType::IndexType start;<br>
start.Fill( 0 );<br>
typename ImportFilterImage3DType::RegionType region;<br>
region.SetIndex( start );<br>
region.SetSize( size );<br>
import->SetRegion( region );<br>
double spacing[ 3 ];<br>
// hier sollte noch das gap mit beachtet werden, welches es nicht in<br>
itk gibt<br>
spacing[0] = diImage->getSpacing((unsigned char)0); // along X<br>
direction<br>
spacing[1] = diImage->getSpacing((unsigned char)1); // along Y<br>
direction<br>
spacing[2] = diImage->getSpacing((unsigned char)2); // along Z<br>
direction<br>
import->SetSpacing( spacing );<br>
double origin[ 3 ];<br>
origin[0] = diImage->getOrigin((unsigned char)0); // X coordinate<br>
origin[1] = diImage->getOrigin((unsigned char)1); // Y coordinate<br>
origin[2] = diImage->getOrigin((unsigned char)2); // Z coordinate<br>
import->SetOrigin( origin );<br>
<br>
const unsigned int numberOfPixels = diImage->getSize();<br>
dType* data = diImage->getData();<br>
import->SetImportPointer( data, numberOfPixels, false );<br>
import->Update();<br>
image = import->GetOutput();<br>
image->Update();<br>
<br>
// cout << image;<br>
return image;<br>
}<br>
<br>
Best requards,<br>
Michael Schildt<br>
<br>
Luis Ibanez schrieb:<br>
<br>
<br>
Hi Michael,<br>
<br>
<br>
This looks like a bug in your code :-(<br>
<br>
<br>
<br>
Typical suspects are:<br>
<br>
<br>
A) The fixed image is not available yet at<br>
the time you call GetLargestPossibleRegion(),<br>
GetOrigin(), GetSpacing().<br>
<br>
You should call fixedIm->Update()<br>
<br>
before you call any of the above methods.<br>
<br>
<br>
We will learn a lot about other possible causes if<br>
you add to your code the following:<br>
<br>
fixedIm->Print( std::cout );<br>
movingIm->Print( std::cout );<br>
resampler->Print( std::cout );<br>
<br>
before the line<br>
<br>
resampler->Update();<br>
<br>
and if you post the resulting output text to the<br>
mailing list.<br>
<br>
<br>
Please let us know what you find.<br>
<br>
<br>
Thanks<br>
<br>
<br>
Luis<br>
<br>
<br>
----------------------------------------------------------------<br>
On Fri, Jun 26, 2009 at 6:42 AM, Michael Schildt<br>
<<a href="mailto:michael.schildt@ifn-magdeburg.de" target="_blank">michael.schildt@ifn-magdeburg.de</a><br>
<mailto:<a href="mailto:michael.schildt@ifn-magdeburg.de" target="_blank">michael.schildt@ifn-magdeburg.de</a>><br>
<mailto:<a href="mailto:michael.schildt@ifn-magdeburg.de" target="_blank">michael.schildt@ifn-magdeburg.de</a><br>
<mailto:<a href="mailto:michael.schildt@ifn-magdeburg.de" target="_blank">michael.schildt@ifn-magdeburg.de</a>>>> wrote:<br>
<br>
Hello!<br>
<br>
I have some trouble with ResampleImageFilter. I'm not sure<br>
if it is a<br>
bug in my program or in ITK. It works as expected in Debug<br>
configuration<br>
but crashes in Release mode. I traced to the source line in ITK<br>
were it<br>
appear in RelWithDebMode configuration. I attach some<br>
information and<br>
source.<br>
<br>
Best reguards,<br>
Michael Schildt<br>
<br>
The method based on the example in the itkSoftwareGuide:<br>
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////<br>
template <typename dataT, typename TrfType><br>
typename itk::Image<dataT,3>::Pointer transformImage(typename<br>
itk::Image<dataT,3>::Pointer fixedIm,<br>
typename<br>
itk::Image<dataT,3>::Pointer movingIm,<br>
typename<br>
TrfType::ConstPointer transform)<br>
{<br>
typedef itk::Image<dataT, 3> ImageType;<br>
typedef TrfType TransformType;<br>
typedef itk::LinearInterpolateImageFunction<ImageType,double><br>
InterpolatorType;<br>
//typedef<br>
itk::WindowedSincInterpolateImageFunction<ImageType,3><br>
InterpolatorType;<br>
<br>
typedef itk::ResampleImageFilter<ImageType,ImageType><br>
ResamplerType;<br>
<br>
typename InterpolatorType::Pointer interpolator =<br>
InterpolatorType::New();<br>
<br>
typename ResamplerType::Pointer resampler =<br>
ResamplerType::New();<br>
resampler->SetInput( movingIm );<br>
<br>
resampler->SetTransform( transform.GetPointer() );<br>
resampler->SetInterpolator( interpolator.GetPointer() );<br>
resampler->SetSize(<br>
fixedIm->GetLargestPossibleRegion().GetSize() );<br>
resampler->SetOutputOrigin( fixedIm->GetOrigin() );<br>
resampler->SetOutputSpacing( fixedIm->GetSpacing() );<br>
resampler->SetDefaultPixelValue( 0 );<br>
<br>
// resample the moving image<br>
resampler->Update();<br>
<br>
return resampler->GetOutput();<br>
}<br>
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////<br>
<br>
It appears in the Update() method of the resampler. While<br>
debugging it<br>
seems like a problem in itkDataObject.cxx at line marked<br>
with **<br>
<br>
//----------------------------------------------------------------------------<br>
void<br>
DataObject<br>
::UpdateOutputData()<br>
{<br>
// If we need to update due to PipelineMTime, or the fact<br>
that our<br>
// data was released, then propagate the UpdateOutputData<br>
to the<br>
source<br>
// if there is one.<br>
if ( m_UpdateMTime < m_PipelineMTime || m_DataReleased ||<br>
this->RequestedRegionIsOutsideOfTheBufferedRegion() )<br>
{<br>
if ( m_Source ) // <---m_Source is valid, but when i<br>
step over<br>
here, it gets an invalid pointer<br>
{<br>
** m_Source->UpdateOutputData(this); // maybe call on<br>
invalid<br>
pointer?<br>
}<br>
}<br>
}<br>
//----------------------------------------------------------------------------<br>
<br>
Stack backtrace:<br>
<br>
BrainGUI.exe!std::_Vector_iterator<unsigned<br>
int,std::allocator<unsigned int> >::operator+() + 0x17<br>
Bytes C++<br>
<br>
BrainGUI.exe!std::vector<unsigned<br>
int,std::allocator<unsigned int> ::resize(unsigned int<br>
_Newsize=1, unsigned int _Val=0) Zeile 721 +<br>
<br>
0x36 Bytes C++<br>
BrainGUI.exe!std::vector<bool,std::allocator<bool><br>
<br>
::_Insert_x(std::_Vb_const_iterator<unsigned<br>
<br>
int,int,std::vector<bool,std::allocator<bool> > ><br>
_Where=..., unsigned<br>
int _Count=1) Zeile 2312 C++<br>
BrainGUI.exe!std::vector<bool,std::allocator<bool><br>
<br>
::_Insert_n(std::_Vb_const_iterator<unsigned<br>
<br>
int,int,std::vector<bool,std::allocator<bool> > ><br>
_Where=..., unsigned<br>
int _Count=1, bool _Val=false) Zeile 2292 C++<br>
BrainGUI.exe!std::vector<bool,std::allocator<bool><br>
<br>
::resize(unsigned int _Newsize=1, bool _Val=false)<br>
Zeile 2044<br>
C++<br>
<br>
BrainGUI.exe!itk::ProcessObject::CacheInputReleaseDataFlags()<br>
Zeile 1055 C++<br>
BrainGUI.exe!itk::ProcessObject::UpdateOutputData(itk::DataObject<br>
*<br>
__formal=0x02079270) Zeile 964 C++<br>
BrainGUI.exe!itk::DataObject::UpdateOutputData() Zeile<br>
420 + 0x8<br>
Bytes C++<br>
<br>
BrainGUI.exe!transformImage<float,itk::VersorRigid3DTransform<double><br>
<br>
(itk::SmartPointer<itk::Image<float,3> > fixedIm={...},<br>
<br>
itk::SmartPointer<itk::Image<float,3> > movingIm={...},<br>
itk::SmartPointer<itk::VersorRigid3DTransform<double> const ><br>
transform={...}) Zeile 123 C++<br>
BrainGUI.exe!VersorRigid3DRegistration::doRegistration(bool<br>
newVolume=true, void (const double &, const<br>
std::basic_string<char,std::char_traits<char>,std::allocator<char><br>
> &)*<br>
progress=0x004029e1) Zeile 254 + 0x81 Bytes C++<br>
BrainGUI.exe!RegThread::Entry() Zeile 273 + 0xf Bytes C++<br>
BrainGUI.exe!wxThreadInternal::DoThreadStart() + 0xa5<br>
Bytes C++<br>
BrainGUI.exe!wxThreadInternal::WinThreadStart() + 0x6a<br>
Bytes C++<br>
BrainGUI.exe!_callthreadstartex() Zeile 348 + 0x6 Bytes<br>
C<br>
BrainGUI.exe!_threadstartex(void * ptd=0x04ef1118)<br>
Zeile 326 + 0x5<br>
Bytes C<br>
<br>
System:<br>
- Windows XP Prof SP3<br>
- Visual Studio 2008 Express Edition<br>
- ITK 3.14<br>
- CMake 2.6.4<br>
- wxWidgets 2.8.9 (wxPack Installer)<br>
- Pentium 4<br>
- 2GB Ram<br>
<br>
<br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a> <<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><br>
<<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><br>
<br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br>
<br>
<br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a> <<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br>
<br>
</div></div></blockquote>
</blockquote></div><br>