- unsigned char* pPix;
- pPix = (unsigned char*)image->GetScalarPointer();
-
- for(int i=0;i<1041;i=i+3){
- int p1 = pPix[i];
- int p2 = pPix[i+1];
- int p3 = pPix[i+2];
-
- double* vecProj;
- vecProj = vectorProjection (Px, Py, Pz, Nx, Ny, Nz, p1, p2, p3);
- double a, b;
-
- //double* cutPoint;
- //cutPoint = intersectionPlaneLine(pPlane[0], pPlane[1], pPlane[2], pPlane[3], pPlane[4], pPlane[5], pPlane[6], pPlane[7], pPlane[8], px1, py1, pz1, px2, py2, pz2);
-
- a = alpha (p1, p2, p3,
- p1, p2, p3, //?? es el vector?
- vecProj[0], vecProj[1], vecProj[2]);
-
- }
-
-
-}
+ alphaImage = image;
+ alphaImage->SetExtent( image->GetExtent() );
+ alphaImage->SetOrigin( image->GetOrigin() );
+ alphaImage->SetSpacing( image->GetSpacing() );
+ alphaImage->SetScalarTypeToUnsignedChar();
+ alphaImage->SetNumberOfScalarComponents( image->GetNumberOfScalarComponents() );
+ alphaImage->AllocateScalars();
+
+ betaImage = image;
+ betaImage->SetExtent( image->GetExtent() );
+ betaImage->SetOrigin( image->GetOrigin() );
+ betaImage->SetSpacing( image->GetSpacing() );
+ betaImage->SetScalarTypeToUnsignedChar();
+ betaImage->SetNumberOfScalarComponents( image->GetNumberOfScalarComponents() );
+ betaImage->AllocateScalars();