{"id":2191,"date":"2024-04-29T16:17:43","date_gmt":"2024-04-29T08:17:43","guid":{"rendered":"https:\/\/sanlangcode.com\/?p=2191"},"modified":"2024-04-29T16:17:43","modified_gmt":"2024-04-29T08:17:43","slug":"vtk%e5%ae%9e%e7%8e%b0%e8%a7%a3%e6%9e%90dcm%e5%ae%9e%e7%8e%b0%e5%90%8c%e7%aa%97%e5%8f%a3%e5%9b%9b%e8%a7%86%e5%9b%be%e5%92%8c%e4%b8%89%e7%bb%b4%e9%87%8d%e5%bb%ba","status":"publish","type":"post","link":"https:\/\/sanlangcode.com\/index.php\/2024\/04\/29\/vtk%e5%ae%9e%e7%8e%b0%e8%a7%a3%e6%9e%90dcm%e5%ae%9e%e7%8e%b0%e5%90%8c%e7%aa%97%e5%8f%a3%e5%9b%9b%e8%a7%86%e5%9b%be%e5%92%8c%e4%b8%89%e7%bb%b4%e9%87%8d%e5%bb%ba\/","title":{"rendered":"VTK\u5b9e\u73b0\u89e3\u6790dcm\u5b9e\u73b0\u540c\u7a97\u53e3\u56db\u89c6\u56fe\u548c\u4e09\u7ef4\u91cd\u5efa"},"content":{"rendered":"\n

VTK\uff08Visualization Toolkit\uff09\u662f\u4e00\u4e2a\u7528\u4e8e 3D \u8ba1\u7b97\u673a\u56fe\u5f62\u5b66\u548c\u56fe\u50cf\u5904\u7406\u7684\u5f00\u6e90\u8f6f\u4ef6\u7cfb\u7edf\u3002\u5b83\u63d0\u4f9b\u4e86\u5e7f\u6cdb\u7684\u529f\u80fd\uff0c\u5305\u62ec\u6570\u636e\u53ef\u89c6\u5316\u3001\u56fe\u50cf\u5904\u7406\u3001\u4f53\u7ed8\u5236\u3001\u8868\u9762\u91cd\u5efa\u3001\u4f53\u7ed8\u5236\u3001\u4f53\u6570\u636e\u7684\u4ea4\u4e92\u5f0f\u5207\u7247\u548c\u4f53\u79ef\u6e32\u67d3\u7b49\u3002VTK \u901a\u5e38\u7528\u4e8e\u533b\u5b66\u56fe\u50cf\u5904\u7406\u3001\u79d1\u5b66\u53ef\u89c6\u5316\u3001\u5de5\u7a0b\u548c\u5730\u7403\u79d1\u5b66\u7b49\u9886\u57df\u3002<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\u5b9e\u73b0\u89e3\u6790 DICOM\uff08Digital Imaging and Communications in Medicine\uff09\u5e76\u5728\u540c\u4e00\u4e2a\u7a97\u53e3\u4e2d\u663e\u793a\u56db\u4e2a\u89c6\u56fe\uff08\u51a0\u72b6\u9762\u3001\u77e2\u72b6\u9762\u3001\u6a2a\u65ad\u9762\u548c\u4e09\u7ef4\u91cd\u5efa\uff09\u7684\u57fa\u672c\u539f\u7406\u5982\u4e0b\uff1a<\/p>\n\n\n\n

  1. DICOM \u6570\u636e\u8bfb\u53d6\uff1a<\/strong> \u9996\u5148\uff0c\u9700\u8981\u4f7f\u7528 VTK \u7684 DICOMImageReader \u7c7b\u6216\u5176\u4ed6\u7c7b\u4f3c\u7684\u5de5\u5177\u8bfb\u53d6 DICOM \u6570\u636e\u96c6\u3002DICOM \u662f\u533b\u5b66\u56fe\u50cf\u7684\u6807\u51c6\u683c\u5f0f\uff0c\u5176\u4e2d\u5305\u542b\u4e86\u6765\u81ea\u533b\u5b66\u5f71\u50cf\u8bbe\u5907\uff08\u5982 CT\u3001MRI\u3001X\u5c04\u7ebf\u7b49\uff09\u7684\u56fe\u50cf\u6570\u636e\u4ee5\u53ca\u76f8\u5173\u7684\u5143\u6570\u636e\u4fe1\u606f\u3002<\/li>
  2. \u5e73\u9762\u90e8\u4ef6\u8bbe\u7f6e\uff1a<\/strong> \u63a5\u4e0b\u6765\uff0c\u9700\u8981\u521b\u5efa\u56db\u4e2a\u5e73\u9762\u90e8\u4ef6\uff0c\u5206\u522b\u7528\u4e8e\u663e\u793a\u51a0\u72b6\u9762\u3001\u77e2\u72b6\u9762\u548c\u6a2a\u65ad\u9762\u7684\u56fe\u50cf\u3002\u6bcf\u4e2a\u5e73\u9762\u90e8\u4ef6\u90fd\u4e0e\u76f8\u5e94\u7684 DICOM \u5207\u7247\u65b9\u5411\u5bf9\u9f50\uff0c\u5e76\u5c06\u5176\u5173\u8054\u5230\u6e32\u67d3\u5668\u4e2d\u3002\u8fd9\u6837\uff0c\u5c31\u53ef\u4ee5\u5728\u540c\u4e00\u4e2a\u7a97\u53e3\u4e2d\u663e\u793a\u56db\u4e2a\u89c6\u56fe\u3002<\/li>
  3. \u4ea4\u4e92\u5f0f\u64cd\u4f5c\uff1a<\/strong> \u8bbe\u7f6e\u4ea4\u4e92\u5f0f\u64cd\u4f5c\uff0c\u4f7f\u7528\u6237\u80fd\u591f\u5728\u89c6\u56fe\u4e2d\u8fdb\u884c\u5e73\u79fb\u3001\u7f29\u653e\u548c\u65cb\u8f6c\u7b49\u64cd\u4f5c\u3002\u8fd9\u901a\u5e38\u901a\u8fc7\u8bbe\u7f6e vtkRenderWindowInteractor \u548c\u76f8\u5e94\u7684\u4e8b\u4ef6\u76d1\u542c\u5668\u6765\u5b9e\u73b0\u3002\u7528\u6237\u53ef\u4ee5\u4f7f\u7528\u9f20\u6807\u6216\u952e\u76d8\u4e0e\u56fe\u50cf\u4ea4\u4e92\uff0c\u4ee5\u4fbf\u67e5\u770b\u611f\u5174\u8da3\u7684\u533a\u57df\u3002<\/li>
  4. \u4e09\u7ef4\u91cd\u5efa\uff1a<\/strong> \u4f7f\u7528\u4f53\u7ed8\u5236\u6216\u8868\u9762\u91cd\u5efa\u7b49\u6280\u672f\u521b\u5efa\u4e09\u7ef4\u6a21\u578b\u3002\u5728\u8fd9\u79cd\u60c5\u51b5\u4e0b\uff0c\u53ef\u4ee5\u4f7f\u7528 vtkMarchingCubes \u6216\u5176\u4ed6\u7c7b\u4f3c\u7684\u7b97\u6cd5\u4ece\u4f53\u6570\u636e\u4e2d\u63d0\u53d6\u7b49\u503c\u9762\u3002\u7136\u540e\uff0c\u5c06\u4e09\u7ef4\u6a21\u578b\u6dfb\u52a0\u5230\u53e6\u4e00\u4e2a\u6e32\u67d3\u5668\u4e2d\uff0c\u4ee5\u4fbf\u5728\u540c\u4e00\u4e2a\u7a97\u53e3\u4e2d\u663e\u793a\u4e09\u7ef4\u91cd\u5efa\u7684\u7ed3\u679c\u3002<\/li>
  5. \u8bbe\u7f6e\u89c6\u53e3\u548c\u6e32\u67d3\u5668\uff1a<\/strong> \u6700\u540e\uff0c\u5c06\u56db\u4e2a\u6e32\u67d3\u5668\u653e\u7f6e\u5728\u540c\u4e00\u4e2a vtkRenderWindow \u4e2d\uff0c\u5e76\u8bbe\u7f6e\u5b83\u4eec\u7684\u89c6\u53e3\uff0c\u4ee5\u4fbf\u5b83\u4eec\u5728\u7a97\u53e3\u4e2d\u6b63\u786e\u6392\u5217\u3002\u53ef\u4ee5\u8bbe\u7f6e\u6e32\u67d3\u5668\u7684\u80cc\u666f\u989c\u8272\u548c\u76f8\u673a\u53c2\u6570\uff0c\u4ee5\u83b7\u5f97\u66f4\u597d\u7684\u89c6\u89c9\u6548\u679c\u3002<\/li><\/ol>\n\n\n\n

    \u901a\u8fc7\u8fd9\u4e9b\u6b65\u9aa4\uff0c\u5c31\u53ef\u4ee5\u5b9e\u73b0\u5728\u540c\u4e00\u4e2a\u7a97\u53e3\u4e2d\u663e\u793a\u56db\u4e2a\u89c6\u56fe\uff08\u51a0\u72b6\u9762\u3001\u77e2\u72b6\u9762\u3001\u6a2a\u65ad\u9762\u548c\u4e09\u7ef4\u91cd\u5efa\uff09\u7684\u529f\u80fd\u3002\u7528\u6237\u53ef\u4ee5\u901a\u8fc7\u4ea4\u4e92\u5f0f\u64cd\u4f5c\u6d4f\u89c8\u56fe\u50cf\uff0c\u5e76\u4e14\u5728\u540c\u4e00\u4e2a\u7a97\u53e3\u4e2d\u540c\u65f6\u67e5\u770b\u4e0d\u540c\u65b9\u5411\u7684\u5207\u7247\u548c\u4e09\u7ef4\u91cd\u5efa\u7ed3\u679c\u3002<\/p>\n\n\n\n

    #include \"vtkSmartPointer.h\"\n#include \"vtkActor.h\"\n#include \"vtkCamera.h\"\n#include \"vtkCellPicker.h\"\n#include \"vtkCommand.h\"\n#include \"vtkImageActor.h\"\n#include \"vtkImageReslice.h\"\n#include \"vtkInteractorStyleImage.h\"\n#include \"vtkImageMapToColors.h\"\n#include \"vtkImagePlaneWidget.h\"\n#include \"vtkImageReader.h\"\n#include \"vtkInteractorEventRecorder.h\"\n#include \"vtkLookupTable.h\"\n#include \"vtkOutlineFilter.h\"\n#include \"vtkDICOMImageReader.h\"\n#include \"vtkPolyDataMapper.h\"\n#include \"vtkProperty.h\"\n#include \"vtkRenderWindow.h\"\n#include \"vtkRenderWindowInteractor.h\"\n#include \"vtkRenderer.h\"\n#include \"vtkImageData.h\"\n#include \"vtkPointData.h\"\n#include \"vtkPlaneSource.h\"\n#include \"vtkPlane.h\"\n#include \"vtkResliceCursorActor.h\"\n#include \"vtkResliceCursorPolyDataAlgorithm.h\"\n#include \"vtkResliceCursor.h\"\n#include \"vtkResliceCursorWidget.h\"\n#include \"vtkResliceCursorLineRepresentation.h\"\n#include \"vtkBiDimensionalWidget.h\"\n#include\"vtkAutoInit.h\"\n#include\"vtkAxesActor.h\"\n#include\"vtkTransform.h\"\n#include\"vtkTextActor.h\"\n#include\"vtkProperty2D.h\"\n\nVTK_MODULE_INIT(vtkRenderingOpenGL2);\nVTK_MODULE_INIT(vtkInteractionStyle);\nVTK_MODULE_INIT(vtkRenderingFreeType);\nVTK_MODULE_INIT(vtkRenderingVolumeOpenGL2);\n\n#include \"vtkSmartPointer.h\"\n#include \"vtkDICOMImageReader.h\"\n#include \"vtkMarchingCubes.h\"\n#include \"vtkPolyDataMapper.h\"\n#include \"vtkActor.h\"\n#include \"vtkRenderer.h\"\n#include \"vtkRenderWindow.h\"\n#include \"vtkRenderWindowInteractor.h\"\n#include <vtkOrientationMarkerWidget.h>\n#include \"Reconstruct3D.cpp\"\n\nclass vtkResliceCursorCallback : public vtkCommand\n{\npublic:\n    static vtkResliceCursorCallback* New()\n    {\n        return new vtkResliceCursorCallback;\n    }\n    void Execute(vtkObject* caller, unsigned long, void* callData) override\n    {\n        vtkImagePlaneWidget* ipw = dynamic_cast<vtkImagePlaneWidget*>(caller);\n        if (ipw)\n        {\n            double* wl = static_cast<double*>(callData);\n\n            if (ipw == this->IPW[0])\n            {\n                this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);\n                this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);\n            }\n            else if (ipw == this->IPW[1])\n            {\n                this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);\n                this->IPW[2]->SetWindowLevel(wl[0], wl[1], 1);\n            }\n            else if (ipw == this->IPW[2])\n            {\n                this->IPW[0]->SetWindowLevel(wl[0], wl[1], 1);\n                this->IPW[1]->SetWindowLevel(wl[0], wl[1], 1);\n            }\n        }\n        vtkResliceCursorWidget* rcw = dynamic_cast<vtkResliceCursorWidget*>(caller);\n        if (rcw)\n        {\n            vtkResliceCursorLineRepresentation* rep = dynamic_cast<vtkResliceCursorLineRepresentation*>(rcw->GetRepresentation());\n            vtkResliceCursor* rc = rep->GetResliceCursorActor()->GetCursorAlgorithm()->GetResliceCursor();\n            for (int i = 0; i < 3; i++)\n            {\n                vtkPlaneSource* ps = static_cast<vtkPlaneSource*>(this->IPW[i]->GetPolyDataAlgorithm());\n                ps->SetNormal(rc->GetPlane(i)->GetNormal());\n                ps->SetCenter(rc->GetPlane(i)->GetOrigin());\n                this->IPW[i]->UpdatePlacement();\n            }\n        }\n        this->RCW[0]->Render();\n    }\n    vtkResliceCursorCallback() {}\n    vtkImagePlaneWidget* IPW[3];\n    vtkResliceCursorWidget* RCW[3];\n};\n\nint FourViews()\n{\n    vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();\n    reader->SetDirectoryName(\"C:\/Users\/DeepA\/Desktop\/CTA_Cardio\");\n    \/\/reader->SetDirectoryName(\"C:\/Users\/DeepA\/Desktop\/RadioPic\/DICOM\/Heart\");\n\n    reader->Update();\n   \n    vtkSmartPointer<vtkPolyDataMapper> outlineMapper = vtkSmartPointer<vtkPolyDataMapper>::New();\n    outlineMapper->SetInputConnection(reader->GetOutputPort());\n\n    vtkSmartPointer<vtkActor> outlineActor = vtkSmartPointer<vtkActor>::New();\n    outlineActor->SetMapper(outlineMapper);\n\n    vtkSmartPointer<vtkRenderer> ren[4];\n    vtkSmartPointer<vtkRenderWindow> renWin = vtkSmartPointer<vtkRenderWindow>::New();\n    renWin->SetMultiSamples(0);\n\n    for (int i = 0; i < 4; i++)\n    {\n        ren[i] = vtkSmartPointer<vtkRenderer>::New();\n        renWin->AddRenderer(ren[i]);\n    }\n\n    vtkSmartPointer<vtkRenderWindowInteractor> iren = vtkSmartPointer<vtkRenderWindowInteractor>::New();\n    iren->SetRenderWindow(renWin);\n\n    vtkSmartPointer<vtkCellPicker> picker = vtkSmartPointer<vtkCellPicker>::New();\n    picker->SetTolerance(0.005);\n\n    vtkSmartPointer<vtkProperty> ipwProp = vtkSmartPointer<vtkProperty>::New();\n\n    vtkSmartPointer<vtkImagePlaneWidget> planeWidget[3];\n    int imageDims[3];\n    reader->GetOutput()->GetDimensions(imageDims);\n\n    for (int i = 0; i < 3; i++)\n    {\n        planeWidget[i] = vtkSmartPointer<vtkImagePlaneWidget>::New();\n        planeWidget[i]->SetInteractor(iren);\n        planeWidget[i]->SetPicker(picker);\n        planeWidget[i]->RestrictPlaneToVolumeOn();\n        double color[3] = { 0, 0, 0 };\n        color[i] = 1;\n        planeWidget[i]->GetPlaneProperty()->SetColor(color);\n        planeWidget[i]->SetTexturePlaneProperty(ipwProp);\n        planeWidget[i]->TextureInterpolateOff();\n        planeWidget[i]->SetResliceInterpolateToLinear();\n        planeWidget[i]->SetInputConnection(reader->GetOutputPort());\n        planeWidget[i]->SetPlaneOrientation(i);\n        planeWidget[i]->SetSliceIndex(imageDims[i] \/ 2);\n        planeWidget[i]->DisplayTextOn();\n        planeWidget[i]->SetDefaultRenderer(ren[3]);\n        planeWidget[i]->SetWindowLevel(1358, -27);\n        planeWidget[i]->On();\n        planeWidget[i]->InteractionOn();\n    }\n\n    planeWidget[1]->SetLookupTable(planeWidget[0]->GetLookupTable());\n    planeWidget[2]->SetLookupTable(planeWidget[0]->GetLookupTable());\n\n    vtkSmartPointer<vtkResliceCursorCallback> cbk = vtkSmartPointer<vtkResliceCursorCallback>::New();\n    vtkSmartPointer< vtkResliceCursor > resliceCursor = vtkSmartPointer< vtkResliceCursor >::New();\n    resliceCursor->SetCenter(reader->GetOutput()->GetCenter());\n    resliceCursor->SetThickMode(0);\n    resliceCursor->SetThickness(10, 10, 10);\n    resliceCursor->SetImage(reader->GetOutput());\n\n    vtkSmartPointer< vtkResliceCursorWidget > resliceCursorWidget[3];\n    vtkSmartPointer< vtkResliceCursorLineRepresentation > resliceCursorRep[3];\n\n    double viewUp[3][3] = { { 0, 0, -1 }, { 0, 0, 1 }, { 0, 1, 0 } };\n    for (int i = 0; i < 3; i++)\n    {\n        resliceCursorWidget[i] = vtkSmartPointer< vtkResliceCursorWidget >::New();\n        resliceCursorWidget[i]->SetInteractor(iren);\n\n        resliceCursorRep[i] = vtkSmartPointer< vtkResliceCursorLineRepresentation >::New();\n        resliceCursorWidget[i]->SetRepresentation(resliceCursorRep[i]);\n        resliceCursorRep[i]->GetResliceCursorActor()->GetCursorAlgorithm()->SetResliceCursor(resliceCursor);\n        resliceCursorRep[i]->GetResliceCursorActor()->GetCursorAlgorithm()->SetReslicePlaneNormal(i);\n\n        const double minVal = reader->GetOutput()->GetScalarRange()[0];\n        if (vtkImageReslice* reslice = vtkImageReslice::SafeDownCast(resliceCursorRep[i]->GetReslice()))\n        {\n            reslice->SetBackgroundColor(minVal, minVal, minVal, minVal);\n        }\n\n        resliceCursorWidget[i]->SetDefaultRenderer(ren[i]);\n        resliceCursorWidget[i]->SetEnabled(1);\n\n        ren[i]->GetActiveCamera()->SetFocalPoint(0, 0, 0);\n        double camPos[3] = { 0, 0, 0 };\n        camPos[i] = 1;\n        ren[i]->GetActiveCamera()->SetPosition(camPos);\n        ren[i]->GetActiveCamera()->ParallelProjectionOn();\n        ren[i]->GetActiveCamera()->SetViewUp(viewUp[i][0], viewUp[i][1], viewUp[i][2]);\n        ren[i]->ResetCamera();\n        cbk->IPW[i] = planeWidget[i];\n        cbk->RCW[i] = resliceCursorWidget[i];\n        resliceCursorWidget[i]->AddObserver(vtkResliceCursorWidget::ResliceAxesChangedEvent, cbk);\n        double range[2];\n        reader->GetOutput()->GetScalarRange(range);\n        resliceCursorRep[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) \/ 2.0);\n        planeWidget[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) \/ 2.0);\n        resliceCursorRep[i]->SetLookupTable(resliceCursorRep[0]->GetLookupTable());\n        planeWidget[i]->GetColorMap()->SetLookupTable(resliceCursorRep[0]->GetLookupTable());\n    }\n  \n    ren[0]->SetBackground(0, 0, 0);\n    ren[1]->SetBackground(0, 0, 0);\n    ren[2]->SetBackground(0, 0, 0);\n    ren[3]->AddActor(outlineActor);\n    ren[3]->SetBackground(0, 0, 0);\n    renWin->SetSize(900, 900);\n\n    ren[0]->SetViewport(0, 0, 0.5, 0.5);\n    ren[1]->SetViewport(0.5, 0, 1, 0.5);\n    ren[2]->SetViewport(0, 0.5, 0.5, 1);\n    ren[3]->SetViewport(0.5, 0.5, 1, 1);\n    renWin->Render();\n      \n    ren[3]->GetActiveCamera()->Elevation(110);\n    ren[3]->GetActiveCamera()->SetViewUp(0, 0, -1);\n    ren[3]->GetActiveCamera()->Azimuth(45);\n    ren[3]->GetActiveCamera()->Dolly(1.15);\n    ren[3]->ResetCameraClippingRange();\n\n    vtkSmartPointer< vtkInteractorStyleImage > style = vtkSmartPointer< vtkInteractorStyleImage >::New();\n\n    iren->SetInteractorStyle(style);\n    iren->Initialize();\n    iren->Start();\n    return EXIT_SUCCESS;\n}\n\n\nint FourViewsAnd3D()\n{\n    vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();\n    reader->SetDirectoryName(\"C:\/Users\/DeepA\/Desktop\/CTA_Cardio\");\n    reader->Update();\n\n    vtkSmartPointer<vtkRenderer> ren[5];\n    vtkSmartPointer<vtkRenderWindow> renWin = vtkSmartPointer<vtkRenderWindow>::New();\n    renWin->SetMultiSamples(0);\n\n    for (int i = 0; i < 5; i++)\n    {\n        ren[i] = vtkSmartPointer<vtkRenderer>::New();\n        renWin->AddRenderer(ren[i]);\n    }\n\n    vtkSmartPointer<vtkRenderWindowInteractor> iren = vtkSmartPointer<vtkRenderWindowInteractor>::New();\n    iren->SetRenderWindow(renWin);\n\n    vtkSmartPointer<vtkCellPicker> picker = vtkSmartPointer<vtkCellPicker>::New();\n    picker->SetTolerance(0.005);\n\n    vtkSmartPointer<vtkProperty> ipwProp = vtkSmartPointer<vtkProperty>::New();\n\n    vtkSmartPointer<vtkImagePlaneWidget> planeWidget[3];\n    int imageDims[3];\n    reader->GetOutput()->GetDimensions(imageDims);\n\n    vtkSmartPointer<vtkResliceCursorCallback> cbk = vtkSmartPointer<vtkResliceCursorCallback>::New();\n    vtkSmartPointer< vtkResliceCursor > resliceCursor = vtkSmartPointer< vtkResliceCursor >::New();\n    resliceCursor->SetCenter(reader->GetOutput()->GetCenter());\n    resliceCursor->SetThickMode(0);\n    resliceCursor->SetThickness(10, 10, 10);\n    resliceCursor->SetImage(reader->GetOutput());\n\n    vtkSmartPointer< vtkResliceCursorWidget > resliceCursorWidget[3];\n    vtkSmartPointer< vtkResliceCursorLineRepresentation > resliceCursorRep[3];\n\n    double viewUp[3][3] = { { 0, 0, -1 }, { 0, 0, 1 }, { 0, 1, 0 } };\n    for (int i = 0; i < 3; i++)\n    {\n        planeWidget[i] = vtkSmartPointer<vtkImagePlaneWidget>::New();\n        planeWidget[i]->SetInteractor(iren);\n        planeWidget[i]->SetPicker(picker);\n        planeWidget[i]->RestrictPlaneToVolumeOn();\n        double color[3] = { 0, 0, 0 };\n        color[i] = 1;\n        planeWidget[i]->GetPlaneProperty()->SetColor(color);\n        planeWidget[i]->SetTexturePlaneProperty(ipwProp);\n        planeWidget[i]->TextureInterpolateOff();\n        planeWidget[i]->SetResliceInterpolateToLinear();\n        planeWidget[i]->SetInputConnection(reader->GetOutputPort());\n        planeWidget[i]->SetPlaneOrientation(i);\n        planeWidget[i]->SetSliceIndex(imageDims[i] \/ 2);\n        planeWidget[i]->DisplayTextOn();\n        planeWidget[i]->SetDefaultRenderer(ren[3]);\n        planeWidget[i]->SetWindowLevel(1358, -27);\n        planeWidget[i]->On();\n        planeWidget[i]->InteractionOn();\n\n        resliceCursorWidget[i] = vtkSmartPointer<vtkResliceCursorWidget>::New();\n        resliceCursorWidget[i]->SetInteractor(iren);\n\n        resliceCursorRep[i] = vtkSmartPointer<vtkResliceCursorLineRepresentation>::New();\n        resliceCursorWidget[i]->SetRepresentation(resliceCursorRep[i]);\n        resliceCursorRep[i]->GetResliceCursorActor()->GetCursorAlgorithm()->SetResliceCursor(resliceCursor);\n        resliceCursorRep[i]->GetResliceCursorActor()->GetCursorAlgorithm()->SetReslicePlaneNormal(i);\n\n        const double minVal = reader->GetOutput()->GetScalarRange()[0];\n        if (vtkImageReslice* reslice = vtkImageReslice::SafeDownCast(resliceCursorRep[i]->GetReslice()))\n        {\n            reslice->SetBackgroundColor(minVal, minVal, minVal, minVal);\n        }\n\n        resliceCursorWidget[i]->SetDefaultRenderer(ren[i]);\n        resliceCursorWidget[i]->SetEnabled(1);\n\n        ren[i]->GetActiveCamera()->SetFocalPoint(0, 0, 0);\n        double camPos[3] = { 0, 0, 0 };\n        camPos[i] = 1;\n        ren[i]->GetActiveCamera()->SetPosition(camPos);\n        ren[i]->GetActiveCamera()->ParallelProjectionOn();\n        ren[i]->GetActiveCamera()->SetViewUp(viewUp[i][0], viewUp[i][1], viewUp[i][2]);\n        ren[i]->ResetCamera();\n        cbk->IPW[i] = planeWidget[i];\n        cbk->RCW[i] = resliceCursorWidget[i];\n        resliceCursorWidget[i]->AddObserver(vtkResliceCursorWidget::ResliceAxesChangedEvent, cbk);\n        double range[2];\n        reader->GetOutput()->GetScalarRange(range);\n        resliceCursorRep[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) \/ 2.0);\n        planeWidget[i]->SetWindowLevel(range[1] - range[0], (range[0] + range[1]) \/ 2.0);\n        resliceCursorRep[i]->SetLookupTable(resliceCursorRep[0]->GetLookupTable());\n        planeWidget[i]->GetColorMap()->SetLookupTable(resliceCursorRep[0]->GetLookupTable());\n    }\n\n    \/\/ Setup 3D reconstruction\n    vtkSmartPointer<vtkMarchingCubes> marchingCubes = vtkSmartPointer<vtkMarchingCubes>::New();\n    marchingCubes->SetInputConnection(reader->GetOutputPort());\n    marchingCubes->SetValue(0, 450); \/\/ Set iso-surface threshold\n\n    vtkSmartPointer<vtkPolyDataMapper> mapper3D = vtkSmartPointer<vtkPolyDataMapper>::New();\n    mapper3D->SetInputConnection(marchingCubes->GetOutputPort());\n\n    vtkSmartPointer<vtkActor> actor3D = vtkSmartPointer<vtkActor>::New();\n    actor3D->SetMapper(mapper3D);\n    ren[4]->AddActor(actor3D);\n\n    \/\/ Set up viewport for the four DICOM slice views and the 3D reconstruction view\n    ren[0]->SetViewport(0.0, 0.5, 0.25, 1.0); \/\/ \u5de6\u4e0a\n    ren[1]->SetViewport(0.25, 0.5, 0.5, 1.0); \/\/ \u53f3\u4e0a\n    ren[2]->SetViewport(0.0, 0.0, 0.25, 0.5); \/\/ \u5de6\u4e0b1\n    ren[3]->SetViewport(0.25, 0.0, 0.5, 0.5); \/\/ \u5de6\u4e0b2\n    ren[4]->SetViewport(0.5, 0.0, 1.0, 1.0); \/\/ \u53f3\u4e0b\n\n    \/\/ \u8bbe\u7f6e\u80cc\u666f\u989c\u8272\n    ren[0]->SetBackground(0.3, 0.0, 0.0); \/\/ \u7ea2\u8272\n    ren[1]->SetBackground(0.0, 0.8, 0.0); \/\/ \u7eff\u8272\n    ren[2]->SetBackground(0.0, 0.0, 0.3); \/\/ \u84dd\u8272\n    ren[3]->SetBackground(0.3, 0.3, 0.0); \/\/ \u9ec4\u8272\n    ren[4]->SetBackground(0, 0, 0); \/\/ \u9ed1\u8272\n\n    renWin->SetSize(1300, 700);\n    renWin->Render();\n\n    \/\/ Start interaction\n    iren->Start();\n\n    return EXIT_SUCCESS;\n}\n\nvoid Reconstruct3D()\n{\n    vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();\n    reader->SetDirectoryName(\"C:\/Users\/DeepA\/Desktop\/CTA_Cardio\");\n    reader->Update();\n\n    vtkSmartPointer<vtkMarchingCubes> marchingCubes = vtkSmartPointer<vtkMarchingCubes>::New();\n    marchingCubes->SetInputConnection(reader->GetOutputPort());\n    marchingCubes->SetValue(0, 450); \/\/ \u8bbe\u7f6e iso-surface \u7684\u9608\u503c\n\n    vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();\n    mapper->SetInputConnection(marchingCubes->GetOutputPort());\n\n    vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();\n    actor->SetMapper(mapper);\n\n    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();\n    renderer->AddActor(actor);\n\n    renderer->SetBackground(0.2, 0.2, 0.2); \/\/ \u9ed1\u7070\u8272\u80cc\u666f\n\n    \/\/ \u8ba1\u7b97\u6a21\u578b\u4e2d\u5fc3\u4f4d\u7f6e\n    double center[3];\n    actor->GetCenter();\n\n    \/\/ \u521b\u5efa\u4e00\u4e2a\u5750\u6807\u8f74\uff0c\u5c06\u5176\u653e\u7f6e\u5728\u6a21\u578b\u4e2d\u5fc3\n    vtkSmartPointer<vtkAxesActor> axes = vtkSmartPointer<vtkAxesActor>::New();\n    axes->SetTotalLength(100, 100, 100);\n    axes->SetPosition(center);\n\n    renderer->AddActor(axes);\n\n    vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();\n    renderWindow->AddRenderer(renderer);\n    renderWindow->SetSize(800,480);\n    renderWindow->SetWindowName(\"VTK Window\");\n\n    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor = vtkSmartPointer<vtkRenderWindowInteractor>::New();\n    renderWindowInteractor->SetRenderWindow(renderWindow);\n\n    \/\/ \u8bbe\u7f6e\u6e32\u67d3\u7a97\u53e3\u8fb9\u6846\u5c5e\u6027\n    vtkSmartPointer<vtkProperty> borderProp = vtkSmartPointer<vtkProperty>::New();\n    borderProp->SetColor(1, 1, 1); \/\/ \u8bbe\u7f6e\u8fb9\u6846\u989c\u8272\u4e3a\u7ea2\u8272\n    borderProp->SetLineWidth(10);   \/\/ \u8bbe\u7f6e\u8fb9\u6846\u7ebf\u6761\u7c97\u7ec6\u4e3a2\n\n    \/\/ \u6e32\u67d3\u5e76\u663e\u793a\u7a97\u53e3\n    renderWindow->Render();\n    renderWindowInteractor->Start();\n}\n\nint main() {\n    \/\/FourViews();\n    FourViewsAnd3D();\n    \/\/Reconstruct3D();\n    \/\/InitializeReconstruct3D();\n    return EXIT_SUCCESS;\n}<\/code><\/pre>\n","protected":false},"excerpt":{"rendered":"
    VTK\uff08Visualizatio...<\/p>\n","protected":false},"author":1,"featured_media":2192,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[45,42,43],"tags":[],"class_list":["post-2191","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-c-2","category-vtk","category-43"],"_links":{"self":[{"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/posts\/2191","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/comments?post=2191"}],"version-history":[{"count":0,"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/posts\/2191\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/"}],"wp:attachment":[{"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/media?parent=2191"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/categories?post=2191"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sanlangcode.com\/index.php\/wp-json\/wp\/v2\/tags?post=2191"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}