[Insight-users] 64 Bit windows problem loading 2D images...

[Neil.Burdett at csiro.au]

Hi,
    We have a product (written in C++) that runs on Linux and Windows (XP) 32 bit. I am currently trying to make it work on a 64 bit windows machine. The application can load a number of different types of medical image formats (i.e. Analyze, DiCOM etc) and displays the results on a WxWidget pane that has 4 quadrants (axial, sagittal, coronal and 3D view). The 3D view allows us to rotate the image, whereas the others are essentially 2D and allow us to move through the image slice by slice. A works okay on 32 bit windows and Linux but on 64 bit windows (XP and Windows 7) the 3D view of the image is the only one that can be loaded, the other screens remain greyed out. So I'm just very confused by this...

We are using :

Visual Studio 9 2008
ITK 3.20.0
VTK 5.6.0
WxWidgets 2.8.10
SigC++ 2.2.3
Glew 1.5.2

The images are loaded using a number of different methods and classes, I'll submit as much as I think is relevant:
void milxImageLoaderTask::Run(void)
{
  milxImageLoader loader;
  m_image = loader.Load(m_filename);
}

milxImage *milxImageLoader::Load( const wxString &filename ) const
{
  //milxImage::ComponentType compType;

  milxImageSource source(filename);
  milxImage3D *image = new milxImage3D();

  try {
    milxImageSource::ImageType::Pointer itkImage = source.GetOutput();
    if( itkImage.IsNull() ) {
      delete image;
      //throw milx::Exception( FILE_READ_ERROR, "Error reading file '%s'", std::string(filename.mb_str()).c_str());
      return 0;
      }
    image->SetImage3D( itkImage );
    image->SetName(filename);
    image->SetFileName(filename);
    itkImage->Update();
    } catch( std::exception & err ) {
    delete image;
    //milxErrorMacro(<< "Exception caught while reading !");
    //milxErrorMacro(<< err.what());
    // Output the error to the terminal
    std::cerr << "milxImage *milxImageLoader: Exception caught while reading: " << err.what() << std::endl;
    //throw milx::Exception( FILE_READ_ERROR, "Error reading file '%s': %s", std::string(filename.mb_str()).c_str() );
    return 0;
    }

  if( image ) {
    image->SetName(filename);
    //image->SetComponentType(compType);
    }
  return image;
}

void milxImageLoaderTask::Finalize(bool complete)
{
  if( complete ) {
    if( m_image == NULL ) {
      // Read failed. Display a semi-useless error message for now until
      // we can get some meaningful context.
      wxMessageBox(_("Unable to load image file: '") + m_filename + _("'"), _("Error loading image"), wxOK|wxICON_ERROR);
      } else {
      milxImageDataMap * imageDataMap = milxViewController::GetInstance()->AddImage(m_image);
      if( m_openPage ) {
        milxViewController::GetInstance()->AddPage(imageDataMap);
        }
      }
    }
}

Where milxImageDataMap is initialised too:
void milxImageDataMap::Init( milxImage *image )
{
  PutImage(IMAGE_KEY_IMAGE, image);

  // Get number of dimensions
  int dimensions[3];
  double spacing[3];
  double origin[3];
  image->GetVTKImageData()->GetSpacing(spacing);
  image->GetVTKImageData()->GetOrigin(origin);
  image->GetVTKImageData()->GetDimensions(dimensions);

  PutFloat(IMAGE_KEY_VOXEL_SIZE_X, spacing[0]);
  PutFloat(IMAGE_KEY_VOXEL_SIZE_Y, spacing[1]);
  PutFloat(IMAGE_KEY_VOXEL_SIZE_Z, spacing[2]);

  // Store the original image dimensions
  PutInt(IMAGE_KEY_SIZE_X, dimensions[0]);
  PutInt(IMAGE_KEY_SIZE_Y, dimensions[1]);
  PutInt(IMAGE_KEY_SIZE_Z, dimensions[2]);

  // Store the original vtk image data bounds (equiv to: dimensions * voxel size)
  vtkFloatingPointType bounds[6];
  image->GetVTKImageData()->GetBounds(bounds);
  PutFloat(IMAGE_KEY_BOUNDS_MIN_X, bounds[0]);
  PutFloat(IMAGE_KEY_BOUNDS_MAX_X, bounds[1]);
  PutFloat(IMAGE_KEY_BOUNDS_MIN_Y, bounds[2]);
  PutFloat(IMAGE_KEY_BOUNDS_MAX_Y, bounds[3]);
  PutFloat(IMAGE_KEY_BOUNDS_MIN_Z, bounds[4]);
  PutFloat(IMAGE_KEY_BOUNDS_MAX_Z, bounds[5]);

  int midpos[3] = { dimensions[0]/2, dimensions[1]/2, dimensions[2]/2 };

  PutInt(IMAGE_KEY_SLICE_X, midpos[0]);
  PutInt(IMAGE_KEY_SLICE_Y, midpos[1]);
  PutInt(IMAGE_KEY_SLICE_Z, midpos[2]);

  // Extract default colour ranges etc
  // TODO optimise this?
  double intensityRange[2];
  image->GetVTKImageData()->GetScalarRange(intensityRange);

  PutFloat(IMAGE_KEY_INTENSITY_LOWER, intensityRange[0]);
  PutFloat(IMAGE_KEY_INTENSITY_UPPER, intensityRange[1]);
  PutInt(IMAGE_KEY_COLOUR_RAMP, VTK_RAMP_SCURVE);
  PutVector4( IMAGE_KEY_COLOUR_LOWER, Vector4(0.0,0.0,0.0,1.0) );
  PutVector4( IMAGE_KEY_COLOUR_UPPER, Vector4(1.0,1.0,1.0,1.0) );
  PutInt(IMAGE_KEY_VOXEL_ACCURACY, 100);
  PutInt(IMAGE_KEY_INTERPOLATION_MODE, -1);

  // set some default mesh-related values
  PutBool(MESH_KEY_COLOUR_BAR, false);
  PutInt(MESH_KEY_SELECT, -1);
  PutInt(MESH_KEY_SHADING, 1);

  // Setup the derived data
}
vtkImageData *milxImageDataMap::GetVtkOutputData() const
{
  if( reslice == NULL ) {
    return GetImage(IMAGE_KEY_IMAGE)->GetVTKImageData();
    } else {
    return reslice->GetOutput();
    }
}

milxImage3D::Image3DType::Pointer milxImageDataMap::GetITKOutputData() const
{
  milxImage *image = GetImage(IMAGE_KEY_IMAGE);
  vtkMatrix4x4 *inverse = GetVtkMatrix4x4(IMAGE_KEY_TRANSFORM_MATRIX);
  milxImage3D *image3D = dynamic_cast<milxImage3D *>(image);
  if( image3D != NULL ) {
    vtkMatrix4x4 *matrix = vtkMatrix4x4::New();
    vtkMatrix4x4::Invert(inverse, matrix);
    typedef itk::ChangeInformationImageFilter<milxImage3D::Image3DType> InfoType;

    InfoType::Pointer info = InfoType::New();
    info->SetInput(image3D->GetImage3D());
    InfoType::OutputImageDirectionType dir = info->GetOutputDirection();
    for( int i=0; i<3; i++ ) {
      for( int j=0; j<3; j++ ) {
        dir(j,i) = matrix->GetElement(j,i);
        }
      }
    info->SetOutputDirection(dir);
    info->ChangeDirectionOn();
    matrix->Delete();

    const Vector3 &trans = GetVector3(IMAGE_KEY_TRANSLATE_VECTOR, Vector3(0,0,0));
    milxImage3D::Image3DType::PointType sourceOrigin = image3D->GetImage3D()->GetOrigin();
    double origin[3] = {sourceOrigin[0] + trans[0],
    sourceOrigin[1] + trans[1],
    sourceOrigin[2] + trans[2]};
    info->SetOutputOrigin(origin);
    info->ChangeOriginOn();
    info->Update();
    return info->GetOutput();
    }
  return NULL;
}

I was just wondering if anyone else had problems with 64 bit windows version (64 bit Linux works fine)


Also, when I tried load a DICOM series I get an Output window shown, but the window is too small to display all the warning information. And when I try to pause the visual studio debugger the window always remains in the background and I'm unable to maximise it. Do these warnings get written somewhere? As all I can read at the moment is:

WARNING in ..\..\..\InsightToolkit-3.20.0\Code\IO\itkGDCMImageIO.cxx line
GDCMImageIO [000000000BC3AF40] the Dicom File (then I suspect the location of the dicom file)

There are hundreds of these warnings scrolling down the list.

Any help much appreciated as I'm not sure how to proceed

Cheers
Neil
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20100915/f08b3183/attachment.htm>