VTK 实现MIP（两种） Activiz

2019独角兽企业重金招聘Python工程师标准>>>

 public partial class Form1 : Form
    {

        #region 私有变量
        private Kitware.VTK.vtkRenderer _render = null;
        private Kitware.VTK.vtkRenderWindow _renWin = null;
        private Kitware.VTK.vtkRenderWindowInteractor _iren = null;
        private Kitware.VTK.vtkDICOMImageReader _reader = null;
        private double opacityWindow = 4096;
        private double opacityLevel = 2048;



        #endregion

        #region 属性
        public string FileDir
        {
            get;
            set;
        }
        #endregion
        public Form1()
        {
            InitializeComponent();
        }

        private void btn_open_Click(object sender, EventArgs e)
        {
            FolderBrowserDialog dlg = new FolderBrowserDialog();
            if (DialogResult.OK == dlg.ShowDialog())
            {
                FileDir = dlg.SelectedPath;
            }
            
        }

        private void renderWindowControl1_Load(object sender, EventArgs e)
        {
            _render = renderWindowControl1.RenderWindow.GetRenderers().GetFirstRenderer();
            _renWin = renderWindowControl1.RenderWindow;
            _iren = renderWindowControl1.RenderWindow.GetInteractor();


        }

        private void btn_render_Click(object sender, EventArgs e)
        {
            if (!String.IsNullOrEmpty(FileDir))
            {
                _reader = new Kitware.VTK.vtkDICOMImageReader();
                _reader.SetDirectoryName(FileDir);
                _reader.SetDataByteOrderToLittleEndian();
                _reader.Update();

                vtkPiecewiseFunction opacityFun = new vtkPiecewiseFunction();
                opacityFun.AddSegment(opacityLevel-0.5*opacityWindow,0.0,
                                      opacityLevel+0.5*opacityWindow,1.0);

                Console.WriteLine("opacityLevel-0.5*opacityWindow:" + (opacityLevel - 0.5 * opacityWindow).ToString());
                Console.WriteLine("opacityLevel+0.5*opacityWindow:" + (opacityLevel + 0.5 * opacityWindow).ToString());


                vtkColorTransferFunction colorFun = new vtkColorTransferFunction();
                colorFun.AddRGBSegment(0.0, 0.0, 0.0,0.0, 255.0, 1.0, 1.0, 1.0);

                vtkVolumeProperty property = new vtkVolumeProperty();
                property.SetInterpolationTypeToLinear();
                property.ShadeOn();
                property.SetColor(colorFun);
                property.SetScalarOpacity(opacityFun);
                


               vtkSmartVolumeMapper mapper = new vtkSmartVolumeMapper();
               // vtkVolumeProMapper mapper = new vtkVolumeProMapper();
               // vtkVolumeRayCastMIPFunction

                mapper.SetBlendModeToMaximumIntensity();
               
                mapper.SetInput(_reader.GetOutput());


                vtkVolume volume = new vtkVolume();
                volume.SetProperty(property);
                volume.SetMapper(mapper);


                _render.AddViewProp(volume);
                _render.ResetCamera();

                _renWin.Render();
                _iren.Initialize();
                _iren.Start();
            }
            else
            {
                MessageBox.Show("请先选定文件夹！");
            }
        }
    }

第二种：

 public partial class Form1 : Form
    {
        #region 私有变量
        private Kitware.VTK.vtkRenderer _render = null;
        private Kitware.VTK.vtkRenderWindow _renWin = null;
        private Kitware.VTK.vtkRenderWindowInteractor _iren = null;
        private Kitware.VTK.vtkDICOMImageReader _reader = null;
        private double opacityWindow = 4096;
        private double opacityLevel = 2048;



        #endregion

        #region 属性
        public string FileDir
        {
            get;
            set;
        }
        #endregion

        public Form1()
        {
            InitializeComponent();
        }

        private void button1_Click(object sender, EventArgs e)
        {
            FolderBrowserDialog dlg = new FolderBrowserDialog();
            if (DialogResult.OK == dlg.ShowDialog())
            {
                FileDir = dlg.SelectedPath;
            }
        }

        private void button2_Click(object sender, EventArgs e)
        {
            if (!String.IsNullOrEmpty(FileDir))
            {
                _reader = new Kitware.VTK.vtkDICOMImageReader();
                _reader.SetDirectoryName(FileDir);
                _reader.SetDataByteOrderToLittleEndian();
                _reader.Update();

                _reader.SetDataScalarTypeToShort();
                _reader.UpdateWholeExtent();
                _reader.GetOutput().UpdateInformation();

                double[] range=_reader.GetOutput().GetScalarRange();//获得颜色映射标量值的范围

                double min = range[0];
                double max = range[1];
                double diff = max - min;
                double slope=255.0 / diff;//斜率
                double inter = -slope * min;//截距
                double shift = inter / slope;

                vtkImageShiftScale vtkImageCast = vtkImageShiftScale.New();
                //With vtkImageShiftScale Pixels are shifted 
                //(a constant value added) and then scaled (multiplied by a scalar. As a convenience,
                //this class allows you to set the output scalar type similar to vtkImageCast.
                //This is because shift scale operations frequently convert data types.

                //如果不进行转换，vtkVolumeRayCastMapper会无法渲染short类型的数据，
                //必须要经过转换将数据转换为vtkVolumeRayCastMapper可以渲染的数据
                vtkImageCast.SetInput(_reader.GetOutput());
                vtkImageCast.SetScale(slope);
                vtkImageCast.SetShift(shift);
                vtkImageCast.SetOutputScalarTypeToUnsignedShort();
                vtkImageCast.Update();

                //获得透明度转换函数的参数
                range=vtkImageCast.GetOutput().GetScalarRange();
                double level = 0.5 * (range[1] + range[0]);
                double window = range[1] - range[0];          

                vtkPiecewiseFunction opacityFun = new vtkPiecewiseFunction();
                //opacityFun.AddPoint(level-window/2,0.0);
             //  opacityFun.AddPoint(level+window/2,1.0);
               // opacityFun.AddSegment(0,0.0,255,1.0);//all

                //opacityFun.AddPoint(73.0, 0.0);
                //opacityFun.AddPoint(89, 1.0);
                //opacityFun.AddPoint(104, 0.0);//bone
                //opacityFun.AddPoint(13, 0.0);
                //opacityFun.AddPoint(32, 0.25);
                //opacityFun.AddPoint(51, 0.0);//lung黑色
                opacityFun.AddPoint(33.0, 0.0);
                opacityFun.AddPoint(49, 0.8);//0.3
                opacityFun.AddPoint(64, 0.0);//soft

            

                //opacityFun.AddPoint(opacityLevel-0.5*opacityWindow,0.0);
                //opacityFun.AddPoint(opacityLevel + 0.5 * opacityWindow, 1.0);


                vtkPiecewiseFunction grayTransferFun = new vtkPiecewiseFunction();
              //  grayTransferFun.AddSegment(level-window/2,0.0,level+window/2,1.0); 
                grayTransferFun.AddSegment(0, 0.0, 255, 1.0); 

                Console.WriteLine("level:{0}\n window :{1}",level,window);
                vtkVolumeProperty property = new vtkVolumeProperty();
                property.SetInterpolationTypeToLinear();
                property.SetScalarOpacity(opacityFun);
                property.SetColor(grayTransferFun);

                vtkVolumeRayCastMIPFunction mipFun = new vtkVolumeRayCastMIPFunction();
                mipFun.SetMaximizeMethodToOpacity();

                vtkVolumeRayCastMapper mapper = new vtkVolumeRayCastMapper();
                mapper.SetVolumeRayCastFunction(mipFun);
                mapper.SetInput(vtkImageCast.GetOutput());
               // mapper.SetInput(_reader.GetOutput());
                //会报错：无法渲染非unsigned char或者unsigned short 类型的数据
                //vtkVolumeRayCastMapper (09C92D38): 
                //Cannot volume render data of type short, only unsigned char or unsigned short.


                vtkVolume volume = new vtkVolume();
                volume.SetProperty(property);
                volume.SetMapper(mapper);


                _render.AddViewProp(volume);
                _render.ResetCamera();

                _renWin.Render();
                _iren.Initialize();
                _iren.Start();
            }
            else
            {
                MessageBox.Show("请先选定文件夹！");
            }
        }

        private void renderWindowControl1_Load(object sender, EventArgs e)
        {
            _render = renderWindowControl1.RenderWindow.GetRenderers().GetFirstRenderer();
            _renWin = renderWindowControl1.RenderWindow;
            _iren = renderWindowControl1.RenderWindow.GetInteractor();

        }
    }

转载于:https://my.oschina.net/u/2009228/blog/536368