VTK例子5-样条曲线CSpline
样条曲线CSpline#include "stdafx.h"#include "vtkSphereSource.h"#include "vtkPolyDataMapper.h"#include "vtkRenderWindow.h"#include "vtkRenderWindowInteractor.h"#include "vtkCamera.h"#include "vtkActor.h"#include "vtkRenderer.h"#include "vtkPolyData.h"#include "vtkBoxWidget.h"#include "vtkCardinalSpline.h"#include "vtkPoints.h"#include "stdio.h"#include <time.h>#include <vtkGlyph3D.h>#include <vtkProperty.h>#include <vtkCellArray.h>#include <vtkTubeFilter.h>#include <vtkPolyLine.h> void main(){ srand(time(NULL)) ; //rand()/RAND_MAX;// 0 - RAND_MAXvtkRenderer *ren = vtkRenderer::New();vtkRenderWindow *renWindow = vtkRenderWindow::New(); renWindow->AddRenderer(ren);renWindow->SetSize( 600, 600 ); vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New(); iren->SetRenderWindow(renWindow); int numberOfInputPoints = 30; vtkCardinalSpline *aSplineX = vtkCardinalSpline::New(); vtkCardinalSpline *aSplineY = vtkCardinalSpline::New(); vtkCardinalSpline *aSplineZ = vtkCardinalSpline::New(); vtkPoints *inputPoints = vtkPoints::New(); double x,y,z;int i;for(i=0;i<numberOfInputPoints;i++){ x = (float)rand()/(float)RAND_MAX; y = (float)rand()/(float)RAND_MAX; z = (float)rand()/(float)RAND_MAX; printf(" x:%f,y:%f,z:%f\n",x,y,z); aSplineX->AddPoint(i, x); aSplineY->AddPoint(i, y); aSplineZ->AddPoint(i, z); inputPoints->InsertPoint(i, x, y, z);} vtkPolyData *inputData =vtkPolyData::New();inputData->SetPoints(inputPoints); vtkSphereSource *balls =vtkSphereSource::New();balls->SetRadius(.01);balls->SetPhiResolution(10);balls->SetThetaResolution(10); vtkGlyph3D *glyphPoints =vtkGlyph3D::New();glyphPoints->SetInput(inputData);glyphPoints->SetSource(balls->GetOutput()); vtkPolyDataMapper *glyphMapper = vtkPolyDataMapper::New();glyphMapper->SetInputConnection(glyphPoints->GetOutputPort()); vtkActor *glyph = vtkActor::New();glyph->SetMapper(glyphMapper);glyph->GetProperty()->SetDiffuseColor(1,0,0);glyph->GetProperty()->SetSpecular(.3);glyph->GetProperty()->SetSpecularPower(30); ///////////////vtkPoints *points = vtkPoints::New(); // Number of points on the splineint numberOfOutputPoints = 400;float t; for(i=0;i<numberOfOutputPoints;i++){ t = (numberOfInputPoints-1.0)/(numberOfOutputPoints-1.0)*i; points->InsertPoint(i, aSplineX->Evaluate(t), aSplineY->Evaluate(t), aSplineZ->Evaluate(t)); printf(" point:%f, %f, %f; \n",points->GetPoint(i), points->GetPoint(i),points->GetPoint(i));} vtkCellArray *lines = vtkCellArray::New(); int a=lines->InsertNextCell(numberOfOutputPoints); for(i=0;i<numberOfOutputPoints;i++) { lines->InsertCellPoint(i); } vtkPolyData *profileData = vtkPolyData::New();profileData->SetPoints(points);profileData->SetLines(lines); // Add thickness to the resulting line.vtkTubeFilter *profileTubes = vtkTubeFilter::New();profileTubes->SetNumberOfSides(8);profileTubes->SetInput(profileData);profileTubes->SetRadius(0.005); vtkPolyDataMapper *profileMapper = vtkPolyDataMapper::New();profileMapper->SetInputConnection(profileTubes->GetOutputPort()); vtkActor *profile = vtkActor::New();profile->SetMapper(profileMapper);profile->GetProperty()->SetDiffuseColor(1,1,0);profile->GetProperty()->SetSpecular(0.3);profile->GetProperty()->SetSpecularPower(30); ////////////// ren->AddActor(glyph);ren->AddActor(profile); iren->Initialize();renWindow->Render(); iren->Start();} 好好好好哈
页:
[1]