C++ (Cpp) RegionType::GetUpperIndex Examples

Example #1

File: Mesh.cpp Project: H-Snoussi/spinalcordtoolbox

void Mesh::save(string filename, ImageType::Pointer image_ref)
{
	vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
	vtkSmartPointer<vtkDoubleArray> pointNormalsArray = vtkSmartPointer<vtkDoubleArray>::New();
	pointNormalsArray->SetNumberOfComponents(3); //3d normals (ie x,y,z)
	pointNormalsArray->SetNumberOfTuples(points_.size());
	CVector3 p, n;
	for (unsigned int i=0; i<points_.size(); i++) {
		p = points_[i]->getPosition();
		n = points_[i]->getNormal();
		points->InsertNextPoint(p[0],p[1],p[2]);
		pointNormalsArray->SetTuple3(i, n[0], n[1], n[2]) ;
	}
	vtkSmartPointer<vtkPolyData> source = vtkSmartPointer<vtkPolyData>::New();
	source->SetPoints(points);
	source->GetPointData()->SetNormals(pointNormalsArray);
	source->Allocate(triangles_.size());
	for (unsigned int i=0; i<triangles_.size(); i+=3) {
		vtkIdType pts[3] = {triangles_[i],triangles_[i+1],triangles_[i+2]};
		source->InsertNextCell(VTK_TRIANGLE,3,pts);
	}
    
    if (image_ref.IsNotNull())
    {
        ImageType::RegionType largestRegion = image_ref->GetLargestPossibleRegion();
        ImageType::IndexType downSliceIndex, downSliceMIndex, upperSliceIndex, upperSliceMIndex;
        downSliceIndex.Fill(0);
        downSliceMIndex.Fill(0);
        downSliceMIndex[1] = 1;
        upperSliceIndex = largestRegion.GetUpperIndex();
        upperSliceMIndex = upperSliceIndex; upperSliceMIndex[1] = upperSliceMIndex[1]-1;
        PointType downSlicePoint, downSliceMPoint, upperSlicePoint, upperSliceMPoint;
        image_ref->TransformIndexToPhysicalPoint(downSliceIndex, downSlicePoint);
        image_ref->TransformIndexToPhysicalPoint(downSliceMIndex, downSliceMPoint);
        image_ref->TransformIndexToPhysicalPoint(upperSliceIndex, upperSlicePoint);
        image_ref->TransformIndexToPhysicalPoint(upperSliceMIndex, upperSliceMPoint);
        
        CVector3 downSliceNormal = CVector3(downSlicePoint[0]-downSliceMPoint[0],downSlicePoint[1]-downSliceMPoint[1],downSlicePoint[2]-downSliceMPoint[2]).Normalize(), upperSliceNormal = CVector3(upperSlicePoint[0]-upperSliceMPoint[0],upperSlicePoint[1]-upperSliceMPoint[1],upperSlicePoint[2]-upperSliceMPoint[2]).Normalize();
        
        vtkSmartPointer<vtkPlane> downPlane = vtkSmartPointer<vtkPlane>::New(), upperPlane = vtkSmartPointer<vtkPlane>::New();
        downPlane->SetOrigin(downSlicePoint[0],downSlicePoint[1],downSlicePoint[2]);
        downPlane->SetNormal(downSliceNormal[0],downSliceNormal[1],downSliceNormal[2]);
        upperPlane->SetOrigin(upperSlicePoint[0],upperSlicePoint[1],upperSlicePoint[2]);
        upperPlane->SetNormal(upperSliceNormal[0],upperSliceNormal[1],upperSliceNormal[2]);
    
        vtkSmartPointer<vtkClipPolyData> downClipper = vtkSmartPointer<vtkClipPolyData>::New();
        downClipper->SetInputData(source);
        downClipper->SetClipFunction(downPlane);
        downClipper->InsideOutOn();
        downClipper->Update();
    
        vtkSmartPointer<vtkClipPolyData> upperClipper = vtkSmartPointer<vtkClipPolyData>::New();
        upperClipper->SetInputData(downClipper->GetOutput());
        upperClipper->SetClipFunction(upperPlane);
        upperClipper->InsideOutOn();
        upperClipper->Update();
        source = upperClipper->GetOutput();
    }
    
	vtkSmartPointer<vtkPolyDataWriter> writer = vtkSmartPointer<vtkPolyDataWriter>::New();
	writer->SetFileName(filename.c_str());
	writer->SetInputData(source);
	writer->Write();
}