medimage 发表于 2015-12-4 22:26:47

VTK例子6-样条曲线KSpline

样条曲线KSpline#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 "vtkKochanekSpline.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 = 10; vtkKochanekSpline*aSplineX = vtkKochanekSpline ::New();vtkKochanekSpline*aSplineY = vtkKochanekSpline ::New();vtkKochanekSpline*aSplineZ = vtkKochanekSpline ::New(); aSplineX-> SetDefaultBias(0);aSplineX-> SetDefaultTension(0);aSplineX-> SetDefaultContinuity(0);aSplineY-> SetDefaultBias(0);aSplineY-> SetDefaultTension(0);aSplineY-> SetDefaultContinuity(0);aSplineZ-> SetDefaultBias(0);aSplineZ-> SetDefaultTension(0);aSplineZ-> SetDefaultContinuity(0); 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();}

13636600535 发表于 2016-5-10 17:23:53

一直在摸索，还是楼主厉害！

ultrasound 发表于 2016-5-12 18:26:06

想看看 啊    想学学啊

dentaldev 发表于 2016-5-13 22:02:30

谢谢楼主！

13636600535 发表于 2016-5-14 04:16:01

正需要，学习学习，谢谢楼主！

oevko 发表于 2016-5-14 18:22:15

这个有用了

medstd 发表于 2016-5-16 02:49:17

正需要，学习学习，谢谢楼主！

13636600230 发表于 2016-5-21 13:30:17

回复看下，！！

medimagedev 发表于 2016-5-24 06:50:37

高手，厉害，向你学习。
页: [1]
查看完整版本: VTK例子6-样条曲线KSpline