#include <pic2vtk.h>
Static Public Member Functions | |
| static vtkImageData * | loadVtkImage (char *fileName) |
| static vtkImageData * | loadImage (char *fileName, bool littleEndian=true) |
| static vtkImageData * | convert (mitkIpPicDescriptor *pic) |
| static vtkImageData * | convert (mitkIpPicDescriptor *pic, unsigned long mask, int move) |
| static mitkIpPicDescriptor * | convert (vtkImageData *vtkImage) |
| static mitkIpPicDescriptor * | convertSigned2Unsigned (mitkIpPicDescriptor *, unsigned long *offset=NULL) |
| static vtkImageData * | convertVectorImage (mitkIpPicDescriptor *pic) |
| static mitkIpPicDescriptor * | convertVectorImage (vtkImageData *vtkImage) |
Definition at line 33 of file pic2vtk.h.
| vtkImageData * Pic2vtk::convert | ( | mitkIpPicDescriptor * | pic ) | [static] |
konvertiert ein pic-Image in ein vtkImageReader
Definition at line 530 of file pic2vtk.cpp.
Referenced by mitk::ImageMapperGL2D::GenerateData(), mitk::ExtractDirectedPlaneImageFilter::GenerateData(), and loadImage().
{
if ( pic == NULL )
return NULL;
vtkImageData *inData = vtkImageData::New();
if ( inData == NULL )
return NULL;
unsigned long size = 0;
if ( pic->dim == 1 ) {
inData->SetDimensions( pic->n[0] -1, 1, 1);
size = pic->n[0];
inData->SetOrigin( ((float) pic->n[0]) / 2.0f, 0, 0 );
} else if ( pic->dim == 2 ) {
inData->SetDimensions( pic->n[0] , pic->n[1] , 1 );
size = pic->n[0] * pic->n[1];
inData->SetOrigin( ((float) pic->n[0]) / 2.0f, ((float) pic->n[1]) / 2.0f, 0 );
} else if ( pic->dim >= 3 ) {
inData->SetDimensions( pic->n[0], pic->n[1], pic->n[2] );
size = pic->n[0] * pic->n[1] * pic->n[2];
// Test
//inData->SetOrigin( (float) pic->n[0] / 2.0f, (float) pic->n[1] / 2.0f, (float) pic->n[2] / 2.0f );
inData->SetOrigin( 0, 0, 0 );
} else {
inData->Delete () ;
return NULL;
}
inData->SetNumberOfScalarComponents(1);
if ( ( pic->type == mitkIpPicInt || pic->type == mitkIpPicUInt ) && pic->bpe == 1 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
size = size>>3;
register unsigned char* s = (unsigned char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 8 ) {
inData->SetScalarType( VTK_CHAR );
inData->AllocateScalars();
register char* s = (char*) pic->data;
register char* a = (char*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 8 ) {
inData->SetScalarType( VTK_UNSIGNED_CHAR );
inData->AllocateScalars();
register unsigned char* s = (unsigned char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 24 ) {
inData->SetNumberOfScalarComponents(3);
inData->SetScalarType( VTK_UNSIGNED_CHAR );
inData->AllocateScalars();
register unsigned char* s = (unsigned char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size*3 ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 16 ) {
inData->SetScalarType( VTK_SHORT );
inData->AllocateScalars();
register short* s = (short*) pic->data;
register short* a = (short*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 16 ) {
inData->SetScalarType( VTK_UNSIGNED_SHORT );
inData->AllocateScalars();
register unsigned short* s = (unsigned short*) pic->data;
register unsigned short* a = (unsigned short*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 32 ) {
inData->SetScalarType( VTK_INT );
inData->AllocateScalars();
register int* s = (int*) pic->data;
register int* a = (int*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 32 ) {
inData->SetScalarType( VTK_UNSIGNED_INT );
inData->AllocateScalars();
register unsigned int* s = (unsigned int*) pic->data;
register unsigned int* a = (unsigned int*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 64 ) {
inData->SetScalarType( VTK_LONG );
inData->AllocateScalars();
register long* s = (long*) pic->data;
register long* a = (long*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 64 ) {
inData->SetScalarType( VTK_UNSIGNED_LONG );
inData->AllocateScalars();
register unsigned long* s = (unsigned long*) pic->data;
register unsigned long* a = (unsigned long*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 32 ) {
inData->SetScalarType( VTK_FLOAT );
inData->AllocateScalars();
register float* s = (float*) pic->data;
register float* a = (float*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 64 ) {
inData->SetScalarType( VTK_DOUBLE );
inData->AllocateScalars();
register double* s = (double*) pic->data;
register double* a = (double*) inData->GetScalarPointer();
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = *s++;
} else {
inData->Delete();
return NULL;
}
return inData;
}
| vtkImageData * Pic2vtk::convert | ( | mitkIpPicDescriptor * | pic, |
| unsigned long | mask, | ||
| int | move | ||
| ) | [static] |
konvertiert ein pic-Image in ein vtkImageReader mask erwartet ein Bittmuster, das im Zielbild eingeblendet wird. move gibt die Verschiebung an, mit der das Bitmuster im Zielbilt plaziert wird.
Definition at line 995 of file pic2vtk.cpp.
{
if ( pic == NULL )
return NULL;
vtkImageData *inData = vtkImageData::New();
if ( inData == NULL )
return NULL;
unsigned long size = 0;
if ( pic->dim == 1 ) {
inData->SetDimensions( pic->n[0] -1, 1, 1);
size = pic->n[0];
inData->SetOrigin( ((float) pic->n[0]) / 2.0f, 0, 0 );
} else if ( pic->dim == 2 ) {
inData->SetDimensions( pic->n[0] , pic->n[1] , 1 );
size = pic->n[0] * pic->n[1];
inData->SetOrigin( ((float) pic->n[0]) / 2.0f, ((float) pic->n[1]) / 2.0f, 0 );
} else if ( pic->dim >= 3 ) {
inData->SetDimensions( pic->n[0], pic->n[1], pic->n[2] );
size = pic->n[0] * pic->n[1] * pic->n[2];
inData->SetOrigin( 0, 0, 0 );
} else {
inData->Delete () ;
return NULL;
}
inData->SetNumberOfScalarComponents(1);
if ( ( pic->type == mitkIpPicInt || pic->type == mitkIpPicUInt ) && pic->bpe == 1 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
size = size>>3;
register unsigned char* s = (unsigned char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicInt && pic->bpe == 8 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register char* s = (char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 8 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register unsigned char* s = (unsigned char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicInt && pic->bpe == 16 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register short* s = (short*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 16 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register unsigned short* s = (unsigned short*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicInt && pic->bpe == 32 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register int* s = (int*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 32 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register unsigned int* s = (unsigned int*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicInt && pic->bpe == 64 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register long* s = (long*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 64 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
register unsigned long* s = (unsigned long*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
register unsigned char v;
for ( register unsigned long i = 0; i < size ; i++ ) {
v = 0;
for ( int j = 0; j < 8; j++ )
v |= (((mask & (*s++))>>move)<<j);
*a++ = v;
}
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 32 ) {
inData->Delete();
return NULL;
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 64 ) {
inData->Delete();
return NULL;
} else {
inData->Delete();
return NULL;
}
return inData;
}
| mitkIpPicDescriptor * Pic2vtk::convert | ( | vtkImageData * | vtkImage ) | [static] |
konvertiert ein vtkImageReader in ein pic-Image
Definition at line 11 of file pic2vtk.cpp.
References mitkIpPicDescriptor.
{
if ( vtkImage == NULL )
return NULL;
mitkIpPicDescriptor* pic = mitkIpPicNew();
if ( pic == NULL )
return NULL;
int dim[3];
vtkImage->GetDimensions(dim);
if ( dim[2] > 1 ) {
pic->dim = 3;
} else if( dim[1] > 1 ) {
pic->dim = 2;
} else if ( dim[0] > 1 ) {
pic->dim = 1;
} else {
return NULL;
}
pic->n[0] = dim[0];
pic->n[1] = dim[1];
pic->n[2] = dim[2];
register unsigned long size = dim[0] * dim[1] * dim[2] * vtkImage->GetNumberOfScalarComponents();
switch ( vtkImage->GetScalarType () ) {
case VTK_BIT: {
pic->type = mitkIpPicInt;
pic->bpe = 8;
unsigned int vtkSize = size>>3;
if (vtkSize%8)
vtkSize++;
register char *s = (char*) vtkImage->GetScalarPointer();
register char *a = (char*) malloc(sizeof(char)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( register unsigned int i = 0; i < vtkSize; i++ ) {
for ( int k=0; k < 8; k++ )
*a++ = 0x80 & (*s<<k);
s++;
}
break;
}
case VTK_CHAR: {
pic->type = mitkIpPicInt;
pic->bpe = 8*vtkImage->GetNumberOfScalarComponents();
register char *s = (char*) vtkImage->GetScalarPointer();
register char *a = (char*) malloc(sizeof(char)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_UNSIGNED_CHAR: {
pic->type = mitkIpPicUInt;
pic->bpe = 8*vtkImage->GetNumberOfScalarComponents();
register unsigned char *s = (unsigned char*) vtkImage->GetScalarPointer();
register unsigned char *a = (unsigned char*) malloc(sizeof(unsigned char)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_SHORT: {
pic->type = mitkIpPicInt;
pic->bpe = 16;
register short *s = (short*) vtkImage->GetScalarPointer();
register short *a = (short*) malloc(sizeof(short)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_UNSIGNED_SHORT: {
pic->type = mitkIpPicUInt;
pic->bpe = 16;
register unsigned short *s = (unsigned short*) vtkImage->GetScalarPointer();
register unsigned short *a = (unsigned short*) malloc(sizeof(unsigned short)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_INT: {
pic->type = mitkIpPicInt;
pic->bpe = 32;
register int *s = (int*) vtkImage->GetScalarPointer();
register int *a = (int*) malloc(sizeof(int)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_UNSIGNED_INT: {
pic->type = mitkIpPicUInt;
pic->bpe = 32;
register unsigned int *s = (unsigned int*) vtkImage->GetScalarPointer();
register unsigned int *a = (unsigned int*) malloc(sizeof(unsigned int)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_LONG: {
pic->type = mitkIpPicInt;
pic->bpe = 64;
register long *s = (long*) vtkImage->GetScalarPointer();
register long *a = (long*) malloc(sizeof(long)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_UNSIGNED_LONG: {
pic->type = mitkIpPicUInt;
pic->bpe = 64;
register unsigned long *s = (unsigned long*) vtkImage->GetScalarPointer();
register unsigned long *a = (unsigned long*) malloc(sizeof(unsigned long)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_FLOAT: {
pic->type = mitkIpPicFloat;
pic->bpe = 32;
register float *s = (float*) vtkImage->GetScalarPointer();
register float *a = (float*) malloc(sizeof(float)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
case VTK_DOUBLE: {
pic->type = mitkIpPicFloat;
pic->bpe = 64;
register double *s = (double*) vtkImage->GetScalarPointer();
register double *a = (double*) malloc(sizeof(double)* size);
if(a==NULL)
{
mitkIpPicFree( pic );
return NULL;
}
pic->data = a;
for ( unsigned long i = 0; i < size; i++ )
*a++ = *s++;
break;
}
default:
mitkIpPicFree( pic );
return NULL;
}
return pic;
}
| mitkIpPicDescriptor * Pic2vtk::convertSigned2Unsigned | ( | mitkIpPicDescriptor * | pic, |
| unsigned long * | offset = NULL |
||
| ) | [static] |
Definition at line 1217 of file pic2vtk.cpp.
References mitkIpPicDescriptor.
{
if ( pic == NULL )
return NULL;
if ( pic->type == mitkIpPicUInt )
return pic;
mitkIpPicDescriptor* newPic = mitkIpPicCopyHeader( pic, NULL );
const unsigned int size = _mitkIpPicElements( pic );
if ( pic->type == mitkIpPicInt && pic->bpe == 8 ) {
newPic->type = mitkIpPicUInt;
newPic->data = new unsigned char[ size ];
register char* s = (char*) pic->data;
register unsigned char* a = (unsigned char*) newPic->data;
register unsigned short localOffset = 2^7;
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = (unsigned char) ((*s++) + localOffset);
if ( offset )
*offset = localOffset;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 16 ) {
newPic->type = mitkIpPicUInt;
newPic->data = new unsigned short[ size ];
register short* s = (short*) pic->data;
register unsigned short* a = (unsigned short*) newPic->data;
register unsigned int localOffset = 1<<15;
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = (unsigned short) ((*s++) + localOffset);
if ( offset )
*offset = localOffset;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 32 ) {
newPic->type = mitkIpPicUInt;
newPic->data = new unsigned int[ size ];
register int* s = (int*) pic->data;
register unsigned int* a = (unsigned int*) newPic->data;
register unsigned long localOffset = 1<<31;
for ( register unsigned long i = 0; i < size ; i++ )
*a++ = (unsigned int)((*s++) + localOffset);
if ( offset )
*offset = localOffset;
} else if ( pic->type == mitkIpPicInt && pic->bpe == 64 ) {
assert(false);
}
return newPic;
}
| vtkImageData * Pic2vtk::convertVectorImage | ( | mitkIpPicDescriptor * | pic ) | [static] |
konvertiert ein vector pic-Image in ein vector vtkImageData
Definition at line 700 of file pic2vtk.cpp.
References COPYVECTORDATAPIC2VTK.
{
std::cout << "Pic2vtk::convertVectorImage()... converting vector pic image to vector vtkImage " << std::endl;
if ( pic == NULL )
return NULL;
vtkImageData *inData = vtkImageData::New();
if ( inData == NULL )
return NULL;
unsigned long size = 0;
unsigned long volumeSize = 0;
int dim = pic->dim - 1;
unsigned int vectorDim = pic->n[dim];
std::cout << " pic dimension = " << dim << std::endl;
std::cout << " pic vector dim = " << vectorDim << std::endl;
if ( dim == 1 ) {
inData->SetDimensions( pic->n[0] -1, 1, 1);
volumeSize = pic->n[0];
size = volumeSize * vectorDim;
inData->SetOrigin( ((float) pic->n[0]) / 2.0f, 0, 0 );
} else if ( dim == 2 ) {
inData->SetDimensions( pic->n[0] , pic->n[1] , 1 );
volumeSize = pic->n[0]* pic->n[1];
size = volumeSize * vectorDim;
//inData->SetOrigin( ((float) pic->n[0]) / 2.0f, ((float) pic->n[1]) / 2.0f, 0 );
inData->SetOrigin( 0,0, 0 );
} else if ( dim >= 3 ) {
inData->SetDimensions( pic->n[0], pic->n[1], pic->n[2] );
volumeSize = pic->n[0]* pic->n[1]* pic->n[2];
size = volumeSize * vectorDim;
inData->SetOrigin( 0, 0, 0 );
} else {
inData->Delete () ;
return NULL;
}
// std::cout << " size = " << size << std::endl;
// std::cout << " volumeSize = " << volumeSize << std::endl;
//inData->SetNumberOfScalarComponents(vectorDim);
inData->SetNumberOfScalarComponents(3);
if ( ( pic->type == mitkIpPicInt || pic->type == mitkIpPicUInt ) && pic->bpe == 1 ) {
inData->SetScalarType( VTK_BIT );
inData->AllocateScalars();
size = size>>3;
std::cout << "WARNING: 1 bit vector not supported yet ..." << std::endl;
// register unsigned char* s = (unsigned char*) pic->data;
// register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
//
// for ( register unsigned long i = 0; i < size ; i++ )
// for ( register unsigned int j = 0; j<vectorDim; j++)
// {
// *a++ = *(s + i*size);
// }
} else if ( pic->type == mitkIpPicInt && pic->bpe == 8 ) {
inData->SetScalarType( VTK_CHAR );
inData->AllocateScalars();
register char* s = (char*) pic->data;
register char* a = (char*) inData->GetScalarPointer();
#define COPYVECTORDATAPIC2VTK \
for ( register unsigned long i = 0; i < volumeSize ; i++ ) \
{ \
for ( register unsigned int j = 0; j<3; j++) \
{ \
if (j<vectorDim) *a++ = *(s + j*volumeSize); else *a++=0; \
} \
s++; \
}
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 8 ) {
inData->SetScalarType( VTK_UNSIGNED_CHAR );
inData->AllocateScalars();
register unsigned char* s = (unsigned char*) pic->data;
register unsigned char* a = (unsigned char*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicInt && pic->bpe == 16 ) {
inData->SetScalarType( VTK_SHORT );
inData->AllocateScalars();
register short* s = (short*) pic->data;
register short* a = (short*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 16 ) {
inData->SetScalarType( VTK_UNSIGNED_SHORT );
inData->AllocateScalars();
std::cout << " memory = " << inData->GetActualMemorySize() << std::endl;
fflush(stdout);
register unsigned short* s = (unsigned short*) pic->data;
register unsigned short* a = (unsigned short*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicInt && pic->bpe == 32 ) {
inData->SetScalarType( VTK_INT );
inData->AllocateScalars();
register int* s = (int*) pic->data;
register int* a = (int*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 32 ) {
inData->SetScalarType( VTK_UNSIGNED_INT );
inData->AllocateScalars();
register unsigned int* s = (unsigned int*) pic->data;
register unsigned int* a = (unsigned int*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicInt && pic->bpe == 64 ) {
inData->SetScalarType( VTK_LONG );
inData->AllocateScalars();
register long* s = (long*) pic->data;
register long* a = (long*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 64 ) {
inData->SetScalarType( VTK_UNSIGNED_LONG );
inData->AllocateScalars();
register unsigned long* s = (unsigned long*) pic->data;
register unsigned long* a = (unsigned long*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 32 ) {
inData->SetScalarType( VTK_FLOAT );
inData->AllocateScalars();
register float* s = (float*) pic->data;
register float* a = (float*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 64 ) {
inData->SetScalarType( VTK_DOUBLE );
inData->AllocateScalars();
register double* s = (double*) pic->data;
register double* a = (double*) inData->GetScalarPointer();
COPYVECTORDATAPIC2VTK
} else {
inData->Delete();
return NULL;
}
return inData;
}
| mitkIpPicDescriptor * Pic2vtk::convertVectorImage | ( | vtkImageData * | vtkImage ) | [static] |
konvertiert ein vektorwertiges vtkImageData in ein vector pic-Image
Definition at line 276 of file pic2vtk.cpp.
References COPYVECTORDATAVTK2PIC, and mitkIpPicDescriptor.
{
std::cout << "Pic2vtk::convertVectorImage()... converting vector vtkImage to vector pic image" << std::endl;
if ( vtkImage == NULL )
return NULL;
std::cout << " vtk dimension: " << vtkImage->GetDataDimension() << std::endl;
std::cout << " vtk vector dim = " << vtkImage->GetNumberOfScalarComponents() << std::endl;
mitkIpPicDescriptor* pic = mitkIpPicNew();
if ( pic == NULL )
return NULL;
int dim[3];
vtkImage->GetDimensions(dim);
unsigned int vectorDim = vtkImage->GetNumberOfScalarComponents();
unsigned long size = 0;
unsigned long volumeSize = 0;
if ( dim[2] > 1 ) {
pic->dim = 3 + 1;
pic->n[0] = dim[0];
pic->n[1] = dim[1];
pic->n[2] = dim[2];
pic->n[3] = vectorDim;
volumeSize = dim[0] * dim[1] * dim[2];
size = volumeSize * vectorDim;
} else if( dim[1] > 1 ) {
pic->dim = 2 + 1;
pic->n[0] = dim[0];
pic->n[1] = dim[1];
pic->n[2] = vectorDim;
volumeSize = dim[0] * dim[1];
size = volumeSize * vectorDim;
} else if ( dim[1] > 1 ) {
pic->dim = 1 + 1;
pic->n[0] = dim[0];
pic->n[1] = vectorDim;
volumeSize = dim[0];
size = volumeSize * vectorDim;
}
switch ( vtkImage->GetScalarType () ) {
case VTK_BIT: {
std::cout << "WARNING: VTK_BIT vector images not supported yet ... " << std::endl;
// pic->type = mitkIpPicInt;
// pic->bpe = 8;
// unsigned int vtkSize = size>>3;
//
// if (vtkSize%8)
// vtkSize++;
//
// register unsigned char *s = (unsigned char*) vtkImage->GetScalarPointer();
// register unsigned char *a = new unsigned char[size];
// pic->data = a;
//
// for ( register unsigned int i = 0; i < vtkSize; i++ ) {
//
// for ( int k=0; k < 8; k++ )
// *a++ = 0x80 & (*s<<k);
//
// s++;
// }
break;
}
case VTK_CHAR: {
pic->type = mitkIpPicInt;
pic->bpe = 8;
register char *s = (char*) vtkImage->GetScalarPointer();
register char *a = (char*) malloc(sizeof(char)* size);;
pic->data = a;
#define COPYVECTORDATAVTK2PIC \
for ( unsigned long i = 0; i < volumeSize; i++ ) \
{ \
for ( unsigned int j = 0; j<vectorDim; j++) \
{ \
*(a+ j*volumeSize) = *s++; \
} \
a++; \
}
COPYVECTORDATAVTK2PIC
break;
}
case VTK_UNSIGNED_CHAR: {
pic->type = mitkIpPicUInt;
pic->bpe = 8;
register unsigned char *s = (unsigned char*) vtkImage->GetScalarPointer();
register unsigned char *a = (unsigned char*) malloc(sizeof(unsigned char)* size);;
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_SHORT: {
pic->type = mitkIpPicInt;
pic->bpe = 16;
register short *s = (short*) vtkImage->GetScalarPointer();
register short *a = (short*) malloc(sizeof(short)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_UNSIGNED_SHORT: {
pic->type = mitkIpPicUInt;
pic->bpe = 16;
register unsigned short *s = (unsigned short*) vtkImage->GetScalarPointer();
register unsigned short *a = (unsigned short*) malloc(sizeof(unsigned short)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_INT: {
pic->type = mitkIpPicInt;
pic->bpe = 32;
register int *s = (int*) vtkImage->GetScalarPointer();
register int *a = (int*) malloc(sizeof(int)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_UNSIGNED_INT: {
pic->type = mitkIpPicUInt;
pic->bpe = 32;
register unsigned int *s = (unsigned int*) vtkImage->GetScalarPointer();
register unsigned int *a = (unsigned int*) malloc(sizeof(unsigned int)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_LONG: {
pic->type = mitkIpPicInt;
pic->bpe = 64;
register long *s = (long*) vtkImage->GetScalarPointer();
register long *a = (long*) malloc(sizeof(long)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_UNSIGNED_LONG: {
pic->type = mitkIpPicUInt;
pic->bpe = 64;
register unsigned long *s = (unsigned long*) vtkImage->GetScalarPointer();
register unsigned long *a = (unsigned long*) malloc(sizeof(unsigned long)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_FLOAT: {
pic->type = mitkIpPicFloat;
pic->bpe = 32;
register float *s = (float*) vtkImage->GetScalarPointer();
register float *a = (float*) malloc(sizeof(float)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
case VTK_DOUBLE: {
pic->type = mitkIpPicFloat;
pic->bpe = 64;
register double *s = (double*) vtkImage->GetScalarPointer();
register double *a = (double*) malloc(sizeof(double)* size);
pic->data = a;
COPYVECTORDATAVTK2PIC
break;
}
default:
mitkIpPicFree( pic );
return NULL;
}
// add tag tp result pic
mitkIpPicTSV_t *tsv;
tsv = (mitkIpPicTSV_t *) malloc( sizeof(mitkIpPicTSV_t) );
strcpy( tsv->tag, "VectorValued" );
tsv->type = mitkIpPicASCII;
tsv->bpe = 8;
tsv->dim = 1;
tsv->value = strdup( "." );
tsv->n[0] = strlen((char *) tsv->value);
mitkIpPicAddTag( pic, tsv );
return pic;
}
| vtkImageData * Pic2vtk::loadImage | ( | char * | fileName, |
| bool | littleEndian = true |
||
| ) | [static] |
läd ein unkomprimiertes pic file
Definition at line 881 of file pic2vtk.cpp.
References convert(), and mitkIpPicDescriptor.
{
if ( fileName == NULL )
return NULL;
mitkIpPicDescriptor* pic = mitkIpPicGetHeader( fileName, NULL );
if ( pic == NULL )
return NULL;
if ( pic->info->pixel_start_in_file == 0 ) {
// komprimierte Date
mitkIpPicFree( pic );
pic = NULL;
pic = mitkIpPicGet( fileName, NULL );
if ( pic == NULL )
return NULL;
vtkImageData* vtkImage = convert( pic );
mitkIpPicFree( pic );
return vtkImage;
}
vtkImageReader *reader = vtkImageReader::New();
if ( reader == NULL )
return NULL;
reader->ReleaseDataFlagOff();
if ( littleEndian == true )
reader->SetDataByteOrderToLittleEndian();
else
reader->SetDataByteOrderToBigEndian();
if ( pic->dim == 1 ) {
reader->SetFileDimensionality(1);
reader->SetDataExtent ( 0, pic->n[0] -1, 0, 0, 0, 0);
reader->SetDataOrigin( ((float) pic->n[0]) / 2.0f, 0, 0 );
} else if ( pic->dim == 2 ) {
reader->SetFileDimensionality(2);
reader->SetDataExtent ( 0, pic->n[0] -1, 0, pic->n[1] -1, 0, 0);
reader->SetDataOrigin( ((float) pic->n[0]) / 2.0f, ((float) pic->n[1]) / 2.0f, 0 );
} else if ( pic->dim >= 3 ) {
reader->SetFileDimensionality(3);
reader->SetDataOrigin( (float) pic->n[0] / 2.0f, (float) pic->n[1] / 2.0f, (float) pic->n[2] / 2.0f );
reader->SetDataExtent( 0, pic->n[0] - 1, 0, pic->n[1] - 1, 0, pic->n[2] - 1 );
} else {
reader->Delete () ;
mitkIpPicFree( pic );
return NULL;
}
// Datatype
if ( ( pic->type == mitkIpPicInt || pic->type == mitkIpPicUInt ) && pic->bpe == 1 ) {
reader->SetDataScalarType( VTK_BIT );
} else if ( pic->type == mitkIpPicInt && pic->bpe == 8 ) {
reader->SetDataScalarType( VTK_CHAR );
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 8 ) {
reader->SetDataScalarType( VTK_UNSIGNED_CHAR );
} else if ( pic->type == mitkIpPicInt && pic->bpe == 16 ) {
reader->SetDataScalarType( VTK_SHORT );
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 16 ) {
reader->SetDataScalarType( VTK_UNSIGNED_SHORT );
} else if ( pic->type == mitkIpPicInt && pic->bpe == 32 ) {
reader->SetDataScalarType( VTK_INT );
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 32 ) {
reader->SetDataScalarType( VTK_UNSIGNED_INT );
} else if ( pic->type == mitkIpPicInt && pic->bpe == 64 ) {
reader->SetDataScalarType( VTK_LONG );
} else if ( pic->type == mitkIpPicUInt && pic->bpe == 64 ) {
reader->SetDataScalarType( VTK_UNSIGNED_LONG );
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 32 ) {
reader->SetDataScalarType( VTK_FLOAT );
} else if ( pic->type == mitkIpPicFloat && pic->bpe == 64 ) {
reader->SetDataScalarType( VTK_DOUBLE );
} else {
reader->Delete();
mitkIpPicFree( pic );
return NULL;
}
reader->SetDataMask(0xffff);
reader->SetHeaderSize ( pic->info->pixel_start_in_file );
reader->SetFileName( fileName );
reader->UpdateWholeExtent();
mitkIpPicFree( pic );
vtkImageData* vtkImage = reader->GetOutput();
reader->Delete();
return vtkImage;
}
| vtkImageData * Pic2vtk::loadVtkImage | ( | char * | fileName ) | [static] |
Definition at line 1285 of file pic2vtk.cpp.
{
vtkStructuredPointsReader* reader = vtkStructuredPointsReader::New();
// vtkDataReader* reader = vtkDataReader::New();
reader->SetFileName( fileName );
reader->Update();
if ( reader->OpenVTKFile() == 0 )
{
std::cout << "something happend ..." << std::endl;
return NULL;
}
// int test = reader->IsFileStructuredPoints(); /* \todo remove line! */
//int a = reader->GetFileType(); /* \todo remove line! */
//char * kk = reader->GetHeader(); /* \todo remove line! */
vtkStructuredPoints* points = reader->GetOutput();
std::cout << "dim: " << ((vtkImageData *)points)->GetDataDimension() << std::endl;
std::cout << "size: " << ((vtkImageData *)points)->GetScalarSize() << std::endl;
std::cout << "min: " << ((vtkImageData *)points)->GetScalarTypeMin () << std::endl;
std::cout << "min: " << ((vtkImageData *)points)->GetScalarTypeMax () << std::endl;
std::cout << "scalars: " << ((vtkImageData *)points)->GetNumberOfScalarComponents () << std::endl;
((vtkImageData *)points)->SetNumberOfScalarComponents (3);
std::cout << "scalars: " << ((vtkImageData *)points)->GetNumberOfScalarComponents () << std::endl;
int dims[3];
((vtkImageData *)points)->GetDimensions (dims);
std::cout << "dim1 = " << dims[0] << " dims2 = " << dims[1] << " dims3 = " << dims[2] << std::endl;
dims[2] = 2;
((vtkImageData *)points)->SetDimensions (dims);
((vtkImageData *)points)->GetDimensions (dims);
std::cout << "dim1 = " << dims[0] << " dims2 = " << dims[1] << " dims3 = " << dims[2] << std::endl;
std::cout << "dim: " << ((vtkImageData *)points)->GetDataDimension() << std::endl;
int extent[6];
((vtkImageData *)points)->GetExtent (extent);
std::cout << "e1 = " << extent[0] << " e2 = " << extent[1] << " e3 = " << extent[2] << std::endl;
std::cout << "e4 = " << extent[3] << " e5 = " << extent[4] << " e6 = " << extent[5] << std::endl;
((vtkImageData *)points)->GetUpdateExtent (extent);
std::cout << "e1 = " << extent[0] << " e2 = " << extent[1] << " e3 = " << extent[2] << std::endl;
std::cout << "e4 = " << extent[3] << " e5 = " << extent[4] << " e6 = " << extent[5] << std::endl;
((vtkImageData *)points)->GetExtent (extent);
((vtkImageData *)points)->SetUpdateExtent (extent);
std::cout << "e1 = " << extent[0] << " e2 = " << extent[1] << " e3 = " << extent[2] << std::endl;
std::cout << "e4 = " << extent[3] << " e5 = " << extent[4] << " e6 = " << extent[5] << std::endl;
((vtkImageData *)points)->GetExtent (extent);
((vtkImageData *)points)->SetWholeExtent (extent);
std::cout << "e1 = " << extent[0] << " e2 = " << extent[1] << " e3 = " << extent[2] << std::endl;
std::cout << "e4 = " << extent[3] << " e5 = " << extent[4] << " e6 = " << extent[5] << std::endl;
vtkFloatingPointType s[3];
((vtkImageData *)points)->GetSpacing (s);
std::cout << "s1 = " << s[0] << " s2 = " << s[1] << " s3 = " << s[2] << std::endl;
((vtkDataObject*)points)->Update();
return (vtkImageData *)points;
//return NULL;
}
1.7.2