小说阅读网,武道至尊帝临小说

新聞中心

這里有您想知道的互聯(lián)網(wǎng)營銷解決方案

VTK繪制等值線

??????VTK 繪制等值線

??使用vtk來繪制等值線，需要使用到vtkContourFilter，它所處理的數(shù)據(jù)需要是規(guī)則網(wǎng)格數(shù)據(jù)，如規(guī)則的正方形網(wǎng)格等。
?通過使用vtkContourFilter::SetValue(i,value)方法，i表示第幾條等值線，從0開始計(jì)算，value表示等值線的值。多邊形的每個(gè)點(diǎn)都有一個(gè)屬性值，兩個(gè)相鄰的屬性值之間要先通過插值，再把這些屬性值作為該函數(shù)的輸入。如：
???????????SetValue(0,5);
畫一條等值為5的等值線。兩個(gè)相鄰的點(diǎn)的屬性值為3和6，則先插值為（3,4,5,6），所以該先應(yīng)該經(jīng)過這條邊的2/3處。
?也可以使用vtkContourFilter：：GenerateValues(count, start, end);來生成一些列的等值線，其中count參數(shù)指定等值線的條數(shù)，而start指定第一條等值線的值，end指定最后一條等值線的值，中間等值線的值則通過計(jì)算得到。

成都創(chuàng)新互聯(lián)公司IDC提供業(yè)務(wù):多線服務(wù)器托管,成都服務(wù)器租用,多線服務(wù)器托管,重慶服務(wù)器租用等四川省內(nèi)主機(jī)托管與主機(jī)租用業(yè)務(wù);數(shù)據(jù)中心含:雙線機(jī)房,BGP機(jī)房,電信機(jī)房,移動(dòng)機(jī)房,聯(lián)通機(jī)房。

????vtkContourFilter::SetValue(…)、GenerateValues(…)、GetNumberOfContours(…)

解釋：
?SetValue和GenerateValues均為設(shè)置等值線值的函數(shù)；這兩個(gè)函數(shù)的不同之處在于SetValue設(shè)置一條等值線值，而GenerateValues對(duì)應(yīng)n條等值線的值，同時(shí)數(shù)量n包括兩個(gè)邊界值，其余n-2條等值線的值自動(dòng)按照線性關(guān)系產(chǎn)生。
?SetValue和GenerateValues函數(shù)根據(jù)先后調(diào)用關(guān)系產(chǎn)生覆蓋，即后調(diào)用的函數(shù)產(chǎn)生的等值線可能覆蓋掉先調(diào)用函數(shù)產(chǎn)生的等值線，這與等值線的條數(shù)(GetNumberOfContours)相關(guān)。
?例如：首先調(diào)用GenerateValues(3, 100, 300)產(chǎn)生3條等值線，分別為100、200、300；然后調(diào)用SetValue(0, 125)，則最終效果為共有3條等值線，分別為125、200、300；如果改為調(diào)用SetValue(3, 400)，則最終效果為共有4條等值線，分別為100、200、300、400。
?SetValue一般用于覆蓋某一條已經(jīng)存在的等值線，或者增加一條等值線；GenerateValues一般重新設(shè)置等值線的條數(shù)。

示例代碼：

#include "vtkAutoInit.h"
VTK_MODULE_INIT( vtkRenderingOpenGL2 );
VTK_MODULE_INIT( vtkInteractionStyle );

#include "vtkSmartPointer.h"
#include "vtkContourFilter.h"
#include "vtkInteractorStyleUser.h"
#include "vtkProperty.h"
#include "vtkOutlineFilter.h"
#include "vtkCommand.h"
#include "vtkSliderWidget.h"
#include "vtkSliderRepresentation.h"
#include "vtkSliderRepresentation3D.h"
#include "vtkImageData.h"
#include "vtkCellArray.h"
#include "vtkFloatArray.h"
#include "vtkPointData.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkXMLImageDataWriter.h"
#include "vtkNamedColors.h"

namespace
{

void CreateData( vtkImageData* data );

class vtkSliderCallback : public vtkCommand
{
public:
    static vtkSliderCallback* New()
    {
        return new vtkSliderCallback;
    }

    virtual void Execute( vtkObject *caller, unsigned long eventId, void *callData )
    {
        vtkSliderWidget* sliderWidget = reinterpret_cast(caller);

        double value = static_cast(sliderWidget->GetRepresentation())->GetValue();
        this->ContourFilter->GenerateValues( 5, 10.0, value );
        //5: 5contour line 
        //10.0 * value: start value
        //value: end value
    }

    vtkSliderCallback() : ContourFilter( NULL ){}
    vtkContourFilter* ContourFilter;
};

}

int main()
{
vtkSmartPointer colors = vtkSmartPointer::New();

vtkSmartPointer imgData = vtkSmartPointer::New();
CreateData( imgData );

//Create an iso surface
vtkSmartPointer contourFilter = vtkSmartPointer::New();
contourFilter->SetInputData( imgData );
contourFilter->GenerateValues( 5, 10, 10 );

//Map the contours to graphical primitives
vtkSmartPointer contourMapper = vtkSmartPointer::New();
contourMapper->SetInputConnection( contourFilter->GetOutputPort() );

//Create an actor for the contours
vtkSmartPointer contourActor = vtkSmartPointer::New();
contourActor->SetMapper( contourMapper );
contourActor->GetProperty()->SetLineWidth( 5 );

//Create the outline
vtkSmartPointer outlineFilter = vtkSmartPointer::New();
outlineFilter->SetInputData( imgData );

vtkSmartPointer outlineMapper = vtkSmartPointer::New();
outlineMapper->SetInputConnection( outlineFilter->GetOutputPort() );
vtkSmartPointer outlineActor = vtkSmartPointer::New();
outlineActor->SetMapper( outlineMapper );
outlineActor->GetProperty()->SetColor( colors->GetColor3d( "Gray" ).GetData() );
outlineActor->GetProperty()->SetLineWidth( 3 );

//Visualize
vtkSmartPointer renderer = vtkSmartPointer::New();
vtkSmartPointer renderWindow = vtkSmartPointer::New();
renderWindow->AddRenderer( renderer );
vtkSmartPointer interactor = vtkSmartPointer::New();
interactor->SetRenderWindow( renderWindow );

renderer->AddActor( contourActor );
renderer->AddActor( outlineActor );
renderer->SetBackground( colors->GetColor3d( "Burlywood" ).GetData() );

vtkSmartPointer sliderRep = vtkSmartPointer::New();
sliderRep->SetMinimumValue( 0.0 );
sliderRep->SetMaximumValue( 30.0 );
sliderRep->SetValue( 10.0 );
sliderRep->SetTitleText( "Contour value" );
sliderRep->SetPoint1InWorldCoordinates( -20, -40, 0 );
sliderRep->SetPoint2InWorldCoordinates( 0, -40, 0 );
sliderRep->SetSliderWidth( 0.2 );
sliderRep->SetLabelHeight( 0.1 );

vtkSmartPointer sliderWidget = vtkSmartPointer::New();
sliderWidget->SetInteractor( interactor );
sliderWidget->SetRepresentation( sliderRep );
sliderWidget->SetAnimationModeToAnimate();
sliderWidget->EnabledOn();

vtkSmartPointer callback = vtkSmartPointer::New();
callback->ContourFilter = contourFilter;

sliderWidget->AddObserver( vtkCommand::InteractionEvent, callback );

vtkSmartPointer interactStyle = vtkSmartPointer::New();
interactor->SetInteractorStyle( interactStyle );

renderWindow->SetSize( 500, 500 );
renderWindow->Render();
interactor->Start();

return EXIT_SUCCESS;

}

namespace
{
    void CreateData( vtkImageData* pData )
    {
        pData->SetExtent( -25, 25, -25, 25, 0, 0 );
        pData->AllocateScalars( VTK_DOUBLE, 1 );
        int* extent = pData->GetExtent();

    for( int y = extent[ 2 ]; y <= extent[ 3 ]; y++ )
    {
        for( int x = extent[ 0 ]; x <= extent[ 1 ]; x++ )
        {
            double* pixel = static_cast(pData->GetScalarPointer( x, y, 0 ));
            pixel[ 0 ] = sqrt( pow( x, 2.0 ) + pow( y, 2.0 ) );
        }
    }

    vtkSmartPointer writer =
        vtkSmartPointer::New();
    writer->SetFileName( "data.vti" );
    writer->SetInputData( pData );
    writer->Write();
    writer->Update();
}

}
結(jié)果如下圖所示：

VTK 繪制等值線

新聞名稱：VTK繪制等值線
文章出自：http://fisionsoft.com.cn/article/jsjspc.html

新聞中心

??????VTK 繪制等值線

其他資訊