betaImage->SetNumberOfScalarComponents( image->GetNumberOfScalarComponents() );
betaImage->AllocateScalars();
- unsigned char* pPix;
- pPix = (unsigned char*)image->GetScalarPointer();
- double a;
- double b;
+ unsigned char* pPix;
+ pPix = (unsigned char*)image->GetScalarPointer();
+ double a;
+ double b;
- int numberOfPoints = dim[0]*dim[1]*dim[2];
+ int numberOfPoints = dim[0]*dim[1]*dim[2];
- vtkDoubleArray *array;
+ vtkDoubleArray *array;
array = (vtkDoubleArray*) image->GetPointData()->GetVectors();
//for(int i = 0 ; i < 8600 ; i++) {
// std::cout<<" arreglo #"<< i << " "<< array->GetValue(i)<<" tuple3 "<< array->GetTuple3(i)[0]<<" "<< array->GetTuple3(i)[1]<<" "<<array->GetTuple3(i)[2]<< std::endl;
//}
- int i, j, k;
- int numTuple=0;
+ int i, j, k;
+ int numTuple=0;
- for( k = 0 ; k < dim[2] ; k++ )
+ for( k = 0 ; k < dim[2] ; k++ )
{
for( j = 0 ; j < dim[1] ; j++ )
+ {
+ for( i = 0 ; i < dim[0] ; i++ )
{
- for( i = 0 ; i < dim[0] ; i++ )
- {
- //unsigned char* pPix = (unsigned char *)image->GetScalarPointer( i , j , k );
-
- double p1 = array->GetTuple3(numTuple)[0];
- double p2 = array->GetTuple3(numTuple)[1];
- double p3 = array->GetTuple3(numTuple)[2];
- numTuple++;
+ //unsigned char* pPix = (unsigned char *)image->GetScalarPointer( i , j , k );
+
+ double p1 = array->GetTuple3(numTuple)[0];
+ double p2 = array->GetTuple3(numTuple)[1];
+ double p3 = array->GetTuple3(numTuple)[2];
+ numTuple++;
- //std::cout << " p1: " << p1 << " p2: " << p2 <<" p3: " << p3 << std::endl;
- //std::cout << " Px: " << Px << " Py: " << Py <<" Pz: " << Pz << std::endl;
- //std::cout << " p2x: " << p2x << " p2y: " << p2y <<" p2z: " << p2z << std::endl;
-
- double* vecProj;
- vecProj = vectorProjection (Nx, Ny, Nz,
- p1, p2, p3);
- //std::cout << " vecProj[0]: " << vecProj[0] << " vecProj[1]: " << vecProj[1] <<" vecProj[2]: " << vecProj[2] << std::endl;
+ //std::cout << " p1: " << p1 << " p2: " << p2 <<" p3: " << p3 << std::endl;
+ //std::cout << " Px: " << Px << " Py: " << Py <<" Pz: " << Pz << std::endl;
+ //std::cout << " p2x: " << p2x << " p2y: " << p2y <<" p2z: " << p2z << std::endl;
+
+ double* vecProj;
+ vecProj = vectorProjection (Nx, Ny, Nz, p1, p2, p3);
+ //std::cout << " vecProj[0]: " << vecProj[0] << " vecProj[1]: " << vecProj[1] <<" vecProj[2]: " << vecProj[2] << std::endl;
a = alpha (Px, Py, Pz,
p1, p2, p3,
vecProj[0], vecProj[1], vecProj[2]);
//if(a>=90&&a<=360) { std::cout << " numTuple: " << numTuple << " a: " << a << " b: " << b << std::endl; }
if(p1==0&&p2==0&&p3==0){
unsigned char *zPtr1 = (unsigned char *) alphaImage->GetScalarPointer( i , j , k );
- *zPtr1 = (unsigned char)255;
+ *zPtr1 = (unsigned char)0;
unsigned char *zPtr2 = (unsigned char *) betaImage->GetScalarPointer( i , j , k );
- *zPtr2 = (unsigned char)255;
+ *zPtr2 = (unsigned char)0;
}else{
unsigned char *zPtr1 = (unsigned char *) alphaImage->GetScalarPointer( i , j , k );
*zPtr1 = (unsigned char)a;
} // for j
} // for k
-
std::cout << "CFT creaVtkHeartAngles::calculateImages End"<<std::endl;
}