[转载]VTK类视频切片移动直线获取灰度值

首页 > 代码库 > [转载]VTK类视频切片移动直线获取灰度值

[转载]VTK类视频切片移动直线获取灰度值

2024-08-10 08:07:28 218人阅读

#include "stdafx.h"
#include "vtkActor.h"
#include "vtkCamera.h"
#include "vtkCommand.h"
#include "vtkImageData.h"
#include "vtkImagePlaneWidget.h"
#include "vtkInteractorEventRecorder.h"
#include "vtkKochanekSpline.h"
#include "vtkOutlineFilter.h"
#include "vtkParametricSpline.h"
#include "vtkPlaneSource.h"
#include "vtkPointData.h"
#include "vtkPolyData.h"
#include "vtkPolyDataMapper.h"
#include "vtkProbeFilter.h"
#include "vtkProperty.h"
#include "vtkProperty2D.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkSplineWidget.h"
#include "vtkTextProperty.h"
#include "vtkVolume16Reader.h"
#include "vtkXYPlotActor.h"

#include "vtkTestUtilities.h"
#include "vtkRegressionTestImage.h"
#include "vtkDebugLeaks.h"

char TSWeventLog[] =
"# StreamVersion 1\n"
"CharEvent 133 125 0 0 98 1 i\n"
"KeyReleaseEvent 133 125 0 0 98 1 i\n"
"MouseMoveEvent 133 125 0 0 0 0 i\n"
"RightButtonPressEvent 133 125 0 0 0 0 i\n"
"MouseMoveEvent 133 123 0 0 0 0 i\n"
"MouseMoveEvent 133 119 0 0 0 0 i\n"
"MouseMoveEvent 132 115 0 0 0 0 i\n"
"MouseMoveEvent 132 111 0 0 0 0 i\n"
"MouseMoveEvent 132 107 0 0 0 0 i\n"
"RightButtonReleaseEvent 132 107 0 0 0 0 i\n"
"MouseMoveEvent 132 129 0 0 0 0 i\n"
"LeftButtonPressEvent 132 129 0 0 0 0 i\n"
"MouseMoveEvent 132 130 0 0 0 0 i\n"
"MouseMoveEvent 132 135 0 0 0 0 i\n"
"MouseMoveEvent 132 143 0 0 0 0 i\n"
"MouseMoveEvent 131 152 0 0 0 0 i\n"
"MouseMoveEvent 130 159 0 0 0 0 i\n"
"MouseMoveEvent 129 165 0 0 0 0 i\n"
"MouseMoveEvent 127 170 0 0 0 0 i\n"
"MouseMoveEvent 125 176 0 0 0 0 i\n"
"MouseMoveEvent 124 181 0 0 0 0 i\n"
"MouseMoveEvent 122 183 0 0 0 0 i\n"
"LeftButtonReleaseEvent 122 183 0 0 0 0 i\n"
"MouseMoveEvent 133 163 0 0 0 0 i\n"
"MiddleButtonPressEvent 133 163 0 0 0 0 i\n"
"MouseMoveEvent 132 161 0 0 0 0 i\n"
"MouseMoveEvent 128 158 0 0 0 0 i\n"
"MouseMoveEvent 124 155 0 0 0 0 i\n"
"MouseMoveEvent 120 151 0 0 0 0 i\n"
"MouseMoveEvent 116 147 0 0 0 0 i\n"
"MouseMoveEvent 118 146 0 0 0 0 i\n"
"MouseMoveEvent 121 148 0 0 0 0 i\n"
"MouseMoveEvent 123 150 0 0 0 0 i\n"
"MouseMoveEvent 125 154 0 0 0 0 i\n"
"MouseMoveEvent 129 158 0 0 0 0 i\n"
"MouseMoveEvent 132 161 0 0 0 0 i\n"
"MouseMoveEvent 134 165 0 0 0 0 i\n"
"MouseMoveEvent 136 168 0 0 0 0 i\n"
"MiddleButtonReleaseEvent 136 168 0 0 0 0 i\n"
"MouseMoveEvent 178 186 0 0 0 0 i\n"
"KeyPressEvent 178 186 -128 0 0 1 Control_L\n"
"MiddleButtonPressEvent 178 186 8 0 0 0 Control_L\n"
"MouseMoveEvent 178 185 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 183 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 181 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 179 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 177 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 175 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 173 8 0 0 0 Control_L\n"
"MouseMoveEvent 179 171 8 0 0 0 Control_L\n"
"MouseMoveEvent 177 169 8 0 0 0 Control_L\n"
"MouseMoveEvent 176 167 8 0 0 0 Control_L\n"
"MouseMoveEvent 174 165 8 0 0 0 Control_L\n"
"MouseMoveEvent 172 164 8 0 0 0 Control_L\n"
"MouseMoveEvent 171 163 8 0 0 0 Control_L\n"
"MiddleButtonReleaseEvent 171 163 8 0 0 0 Control_L\n"
"KeyReleaseEvent 171 163 0 0 0 1 Control_L\n"
"MouseMoveEvent 170 167 0 0 0 0 Control_L\n"
"MiddleButtonPressEvent 170 167 0 0 0 0 Control_L\n"
"MouseMoveEvent 172 167 0 0 0 0 Control_L\n"
"MouseMoveEvent 176 167 0 0 0 0 Control_L\n"
"MouseMoveEvent 181 167 0 0 0 0 Control_L\n"
"MouseMoveEvent 188 167 0 0 0 0 Control_L\n"
"MouseMoveEvent 198 165 0 0 0 0 Control_L\n"
"MouseMoveEvent 205 163 0 0 0 0 Control_L\n"
"MouseMoveEvent 211 161 0 0 0 0 Control_L\n"
"MouseMoveEvent 216 160 0 0 0 0 Control_L\n"
"MouseMoveEvent 222 158 0 0 0 0 Control_L\n"
"MiddleButtonReleaseEvent 222 158 0 0 0 0 Control_L\n"
"MouseMoveEvent 230 158 0 0 0 0 Control_L\n"
"MiddleButtonPressEvent 230 158 0 0 0 0 Control_L\n"
"MouseMoveEvent 229 156 0 0 0 0 Control_L\n"
"MouseMoveEvent 228 153 0 0 0 0 Control_L\n"
"MouseMoveEvent 226 150 0 0 0 0 Control_L\n"
"MouseMoveEvent 224 148 0 0 0 0 Control_L\n"
"MouseMoveEvent 222 145 0 0 0 0 Control_L\n"
"MouseMoveEvent 220 141 0 0 0 0 Control_L\n"
"MouseMoveEvent 216 135 0 0 0 0 Control_L\n"
"MouseMoveEvent 214 129 0 0 0 0 Control_L\n"
"MouseMoveEvent 212 123 0 0 0 0 Control_L\n"
"MouseMoveEvent 209 118 0 0 0 0 Control_L\n"
"MouseMoveEvent 207 113 0 0 0 0 Control_L\n"
"MouseMoveEvent 204 109 0 0 0 0 Control_L\n"
"MouseMoveEvent 202 105 0 0 0 0 Control_L\n"
"MouseMoveEvent 200 103 0 0 0 0 Control_L\n"
"MouseMoveEvent 198 99 0 0 0 0 Control_L\n"
"MouseMoveEvent 196 97 0 0 0 0 Control_L\n"
"MouseMoveEvent 194 93 0 0 0 0 Control_L\n"
"MiddleButtonReleaseEvent 194 93 0 0 0 0 Control_L\n"
"MouseMoveEvent 254 98 0 0 0 0 Control_L\n"
"MiddleButtonPressEvent 254 98 0 0 0 0 Control_L\n"
"MouseMoveEvent 254 100 0 0 0 0 Control_L\n"
"MouseMoveEvent 254 104 0 0 0 0 Control_L\n"
"MouseMoveEvent 255 108 0 0 0 0 Control_L\n"
"MouseMoveEvent 255 112 0 0 0 0 Control_L\n"
"MouseMoveEvent 255 116 0 0 0 0 Control_L\n"
"MouseMoveEvent 255 120 0 0 0 0 Control_L\n"
"MouseMoveEvent 256 124 0 0 0 0 Control_L\n"
"MouseMoveEvent 257 128 0 0 0 0 Control_L\n"
"MouseMoveEvent 257 132 0 0 0 0 Control_L\n"
"MouseMoveEvent 257 136 0 0 0 0 Control_L\n"
"MouseMoveEvent 258 141 0 0 0 0 Control_L\n"
"MouseMoveEvent 258 146 0 0 0 0 Control_L\n"
"MouseMoveEvent 258 151 0 0 0 0 Control_L\n"
"MouseMoveEvent 258 157 0 0 0 0 Control_L\n"
"MouseMoveEvent 258 159 0 0 0 0 Control_L\n"
"MiddleButtonReleaseEvent 80 206 0 0 0 0 Control_L\n"
;

// Callback for the image plane widget interaction
class vtkIPWCallback : public vtkCommand
{
public:
static vtkIPWCallback *New()
    { return new vtkIPWCallback; }
virtual void Execute(vtkObject *caller, unsigned long, void*)
    {
      vtkImagePlaneWidget *planeWidget = reinterpret_cast<vtkImagePlaneWidget*>(caller);
      if(planeWidget->GetPlaneOrientation() == 3)
        {
        Spline->SetProjectionPosition(0);
        }
      else
        {
        Spline->SetProjectionPosition(planeWidget->GetSlicePosition());
        }
      Spline->GetPolyData(Poly);
    }
vtkIPWCallback():Spline(0),Poly(0){};
vtkSplineWidget* Spline;
vtkPolyData* Poly;
};

// Callback for the spline widget interaction
class vtkSWCallback : public vtkCommand
{
public:
static vtkSWCallback *New()
    { return new vtkSWCallback; }
virtual void Execute(vtkObject *caller, unsigned long, void*)
    {
      vtkSplineWidget *spline = reinterpret_cast<vtkSplineWidget*>(caller);
      spline->GetPolyData(Poly);
    }
vtkSWCallback():Poly(0){};
vtkPolyData* Poly;
};

int main( int argc, char *argv[] )
{
// char* fname = vtkTestUtilities::ExpandDataFileName(argc, argv, "Data/headsq/quarter");

vtkVolume16Reader* v16 = vtkVolume16Reader::New();
    v16->SetDataDimensions( 64, 64);
    v16->SetDataByteOrderToLittleEndian();
    v16->SetImageRange( 1, 93);
    v16->SetDataSpacing( 3.2, 3.2, 1.5);
v16->SetFilePrefix("e://VTK 5.6/VTKDATA/Data/headsq/quarter");
    v16->SetDataMask( 0x7fff);
    v16->Update();

// delete[] fname;

vtkRenderer* ren1 = vtkRenderer::New();
vtkRenderer* ren2 = vtkRenderer::New();
vtkRenderWindow* renWin = vtkRenderWindow::New();
    renWin->SetMultiSamples(0);
    renWin->AddRenderer( ren1);
    renWin->AddRenderer( ren2);
vtkRenderWindowInteractor* iren = vtkRenderWindowInteractor::New();
    iren->SetRenderWindow( renWin);

vtkOutlineFilter* outline = vtkOutlineFilter::New();
outline->SetInputConnection(v16->GetOutputPort());

vtkPolyDataMapper* outlineMapper = vtkPolyDataMapper::New();
outlineMapper->SetInputConnection(outline->GetOutputPort());

vtkActor* outlineActor = vtkActor::New();
outlineActor->SetMapper(outlineMapper);

vtkImagePlaneWidget* ipw = vtkImagePlaneWidget::New();
    ipw->DisplayTextOn();
    ipw->TextureInterpolateOff();
    ipw->UserControlledLookupTableOff();
    ipw->SetInput(v16->GetOutput());
    ipw->KeyPressActivationOn();
    ipw->SetKeyPressActivationValue(‘x‘);
    ipw->SetResliceInterpolateToNearestNeighbour();
    ipw->SetInteractor(iren);
    ipw->SetPlaneOrientationToXAxes();
    ipw->SetSliceIndex(32);
    ipw->GetPlaneProperty()->SetColor(1,0,0);

vtkSplineWidget* spline = vtkSplineWidget::New();
    spline->SetInteractor( iren);
    spline->SetInput(v16->GetOutput());
    spline->SetPriority(1.0);
    spline->KeyPressActivationOff();
    spline->PlaceWidget();
    spline->ProjectToPlaneOn();
    spline->SetProjectionNormal(0);
    spline->SetProjectionPosition(102.4); //initial plane oriented position
    spline->SetProjectionNormal(3); //allow arbitrary oblique orientations
    spline->SetPlaneSource(
      static_cast<vtkPlaneSource*>(ipw->GetPolyDataAlgorithm()));

// Specify the type of spline (change from default vtkCardinalSpline)
vtkKochanekSpline* xspline = vtkKochanekSpline::New();
vtkKochanekSpline* yspline = vtkKochanekSpline::New();
vtkKochanekSpline* zspline = vtkKochanekSpline::New();

vtkParametricSpline* para = spline->GetParametricSpline();

para->SetXSpline(xspline);
para->SetYSpline(yspline);
para->SetZSpline(zspline);

xspline->Delete();
yspline->Delete();
zspline->Delete();

vtkPolyData* poly = vtkPolyData::New();
spline->GetPolyData(poly);

vtkProbeFilter* probe = vtkProbeFilter::New();
probe->SetInput(poly);
probe->SetSource(v16->GetOutput());

vtkIPWCallback* ipwcb = vtkIPWCallback::New();
ipwcb->Spline = spline;
ipwcb->Poly = poly;

ipw->AddObserver(vtkCommand::InteractionEvent,ipwcb);

vtkSWCallback* swcb = vtkSWCallback::New();
swcb->Poly = poly;

spline->AddObserver(vtkCommand::InteractionEvent,swcb);

vtkImageData* data = http://www.mamicode.com/v16->GetOutput();
double* range = data->GetPointData()->GetScalars()->GetRange();

vtkXYPlotActor* profile = vtkXYPlotActor::New();
    profile->AddInput(probe->GetOutput());
    profile->GetPositionCoordinate()->SetValue( 0.05, 0.05, 0);
    profile->GetPosition2Coordinate()->SetValue( 0.95, 0.95, 0);
    profile->SetXValuesToNormalizedArcLength();
    profile->SetNumberOfXLabels( 6);
    profile->SetTitle( "Profile Data ");
    profile->SetXTitle( "s");
    profile->SetYTitle( "I(s)");
    profile->SetXRange( 0, 1);
    profile->SetYRange( range[0], range[1]);
    profile->GetProperty()->SetColor( 0, 0, 0);
    profile->GetProperty()->SetLineWidth( 2);
    profile->SetLabelFormat("%g");
    vtkTextProperty* tprop = profile->GetTitleTextProperty();
    tprop->SetColor(0.02,0.06,0.62);
    tprop->SetFontFamilyToArial();
    profile->SetAxisTitleTextProperty(tprop);
    profile->SetAxisLabelTextProperty(tprop);
    profile->SetTitleTextProperty(tprop);

ren1->SetBackground( 0.1, 0.2, 0.4);
ren1->SetViewport( 0, 0, 0.5, 1);
ren1->AddActor(outlineActor);

ren2->SetBackground( 1, 1, 1);
ren2->SetViewport( 0.5, 0, 1, 1);
ren2->AddActor2D( profile);

renWin->SetSize( 600, 300);

ipw->On();
ipw->SetInteraction(0);
ipw->SetInteraction(1);
spline->On();
spline->SetNumberOfHandles(4);
spline->SetNumberOfHandles(5);
spline->SetResolution(399);

// Set up an interesting viewpoint
vtkCamera* camera = ren1->GetActiveCamera();
camera->Elevation(110);
camera->SetViewUp(0, 0, -1);
camera->Azimuth(45);
camera->SetFocalPoint(100.8,100.8,69);
camera->SetPosition(560.949, 560.949, -167.853);
ren1->ResetCameraClippingRange();

// Position the actors
// renWin->Render();
// iren->SetEventPosition(200,200);
// iren->SetKeyCode(‘r‘);
// iren->InvokeEvent(vtkCommand::CharEvent,NULL);
// ren1->ResetCameraClippingRange();
// renWin->Render();
// iren->SetKeyCode(‘t‘);
// iren->InvokeEvent(vtkCommand::CharEvent,NULL);

// Playback recorded events
vtkInteractorEventRecorder *recorder = vtkInteractorEventRecorder::New();
recorder->SetInteractor(iren);
recorder->ReadFromInputStringOn();
recorder->SetInputString(TSWeventLog);

// Test On Off mechanism
ipw->SetEnabled(0);
spline->EnabledOff();
ipw->SetEnabled(1);
spline->EnabledOn();

// Test Set Get handle positions
double pos[3];
int i;
for(i=0;i<spline->GetNumberOfHandles();i++)
    {
    spline->GetHandlePosition(i,pos);
    spline->SetHandlePosition(i,pos);
    }

// Test Closed On Off
spline->ClosedOn();
spline->ClosedOff();

// Render the image
iren->Initialize();
renWin->Render();
//ren1->ResetCamera();
recorder->Play();

// Remove the observers so we can go interactive. Without this the "-I"
// testing option fails.
recorder->Off();

int retVal = vtkRegressionTestImage( renWin );
if ( retVal == vtkRegressionTester::DO_INTERACTOR)
    {
    iren->Start();
    }

// Clean up
recorder->Off();
recorder->Delete();

outlineActor->Delete();
outlineMapper->Delete();
outline->Delete();
renWin->Delete();
ren1->Delete();
ren2->Delete();
iren->Delete();

ipw->RemoveObserver(ipwcb);
ipw->Delete();
ipwcb->Delete();
spline->RemoveObserver(swcb);
spline->Delete();
swcb->Delete();
poly->Delete();
probe->Delete();
profile->Delete();
v16->Delete();

return !retVal;
}

[转载]VTK类视频切片移动直线获取灰度值

声明：以上内容来自用户投稿及互联网公开渠道收集整理发布，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任，若内容有误或涉及侵权可进行投诉：投诉/举报工作人员会在5个工作日内联系你，一经查实，本站将立刻删除涉嫌侵权内容。

联系
我们

首页 > 代码库 > [转载]VTK类视频 切片移动直线获取灰度值

[转载]VTK类视频 切片移动直线获取灰度值

看完仍有疑问？有类似问题直接问程序猿

首页 > 代码库 > [转载]VTK类视频切片移动直线获取灰度值

[转载]VTK类视频切片移动直线获取灰度值