VTK例子6-样条曲线KSpline

medimage

样条曲线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_MAX

  vtkRenderer *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 spline

int 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)[0],

                                    points->GetPoint(i)[1],points->GetPoint(i)[2]);

}

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();

}

样条曲线KSpline

13636600535

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

ultrasound

想看看 啊    想学学啊

dentaldev

谢谢楼主！

13636600535

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

oevko

这个有用了

medstd

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

13636600230

回复看下，！！

medimagedev

高手，厉害，向你学习。