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 << " 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;