MC.zip_VTK python_python vtk_python 医学图像_vtk_医学图像

import pydicom import vtk from vtk import vtkInteractionImage from vtk import vtkInteractionStyle import time # Create the RenderWindow, Renderer and both Actors # coneActor = vtk.vtkActor() # coneActor.SetMapper( coneMapper ) ren1 = vtk.vtkRenderer() renWin1 = vtk.vtkRenderWindow() renWin1.AddRenderer(ren1) iren = vtk.vtkRenderWindowInteractor() iren.SetRenderWindow(renWin1) style = vtk.vtkInteractorStyleTrackballCamera() iren.SetInteractorStyle(style) # set the UI style lgt = vtk.vtkLight() # input the .dicom v16 = vtk.vtkDICOMImageReader() v16.SetDataByteOrderToLittleEndian() v16.SetDirectoryName('D:\work\SRS00003') v16.SetDataSpacing(3.2, 3.2, 0.5) locator = vtk.vtkMergePoints() locator.SetDivisions(32, 32, 46) locator.SetNumberOfPointsPerBucket(100) locator.AutomaticOff() iso = vtk.vtkMarchingCubes() iso.SetInputConnection(v16.GetOutputPort()) iso.SetValue(0, 150) iso.ComputeGradientsOn() iso.ComputeScalarsOff() iso.SetLocator(locator) # smooth smoother = vtk.vtkWindowedSincPolyDataFilter() smoother.SetInputConnection(iso.GetOutputPort()) smoother.SetNumberOfIterations(90) normals = vtk.vtkPolyDataNormals() normals.SetInputConnection(smoother.GetOutputPort()) # normals.SetInputConnection(extrator.GetOutputPort()) normals.SetFeatureAngle(60.0) # Rang = gradient.GetOutputPort().GetScalarRange() isoMapper = vtk.vtkDataSetMapper() isoMapper.SetInputConnection(normals.GetOutputPort()) isoMapper.ScalarVisibilityOff() isoMapper.SetScalarRange(0, 1200) isoMapper.SetColorModeToDefault() isoMapper.ImmediateModeRenderingOn() isoActor = vtk.vtkActor() isoActor.SetMapper(isoMapper) isoActor.GetProperty().SetDiffuseColor(.1, .94, .52) isoActor.GetProperty().SetSpecular(.3) isoActor.GetProperty().SetSpecularPower(20) isoActor.GetProperty().SetColor(1, 0.1, .2)# R G B outline = vtk.vtkOutlineFilter() outline.SetInputConnection(v16.GetOutputPort()) outlineMapper = vtk.vtkPolyDataMapper() outlineMapper.SetInputConnection(outline.GetOutputPort()) outlineActor = vtk.vtkActor() outlineActor.SetMapper(outlineMapper) # outlineProp = set(outlineActor.GetProperty()) lgt.SetPosition(1, 1, 1) ren1.AddActor(outlineActor) ren1.AddActor(isoActor) ren1.SetBackground(1, 1, 1) ren1.AddLight(lgt) renWin1.SetSize(250, 250) ren1.SetBackground(1, 1, 1) aCamera = vtk.vtkCamera() aCamera.SetViewUp(0, 0, 1) aCamera.SetPosition(0, 1, 0) aCamera.SetFocalPoint(0, 0, 0) # cam1 = set(ren1.GetActiveCamera) ren1.ResetCamera() aCamera.Dolly(1.5) ren1.ResetCameraClippingRange() renWin1.Render() # show the window iren.Initialize() iren.Start()