# How split image data with plane form DICOM images

**URL:** https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955
**Category:** Support
**Created:** [August 12, 2020, 6:46am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955 "2020-08-12T06:46:26Z")
**Posts on this page:** 6
**Page:** 1

<div class="post-metadata">

### Author: ![11112](https://discourse.vtk.org/letter_avatar_proxy/v4/letter/1/d78d45/32.png) [@11112](https://discourse.vtk.org/u/11112)
#### Post date: [August 12, 2020, 6:46am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955/1 "2020-08-12T06:46:26Z")

</div>

I already have a series of DICOM images as CT of verts, I read it with vtkDICOMReader.  
And use volume ray cast to transform as 3D image.  
Now I need to split each vert by plane, and get slice of 3 dimensions.  
How can I do that?

---

<div class="post-metadata">

### Author: ![flaviu2](https://discourse.vtk.org/letter_avatar_proxy/v4/letter/f/74df32/32.png) [@flaviu2](https://discourse.vtk.org/u/flaviu2)
#### Post date: [August 12, 2020, 8:24am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955/2 "2020-08-12T08:24:11Z")

</div>

Are you try to say to do multi planar reconstruction ? If you do, you can find several documentation regarding this, here is a good one: [https://github.com/srinivasrvaidya/Multi-Planar-Reconstruction-using-VTK/blob/master/MPR/src/slicer.py](https://github.com/srinivasrvaidya/Multi-Planar-Reconstruction-using-VTK/blob/master/MPR/src/slicer.py)  
Of course, it’s in Python, but you can easily translate it in C++ if you wish.

P.S\> Or here: [https://gitlab.kitware.com/vtk/vtk/blob/77895235b0c99b5cf2bb98af8b1e923209e30525/Examples/ImageProcessing/Cxx/ImageSlicing.cxx](https://gitlab.kitware.com/vtk/vtk/blob/77895235b0c99b5cf2bb98af8b1e923209e30525/Examples/ImageProcessing/Cxx/ImageSlicing.cxx)

---

<div class="post-metadata">

### Author: ![11112](https://discourse.vtk.org/letter_avatar_proxy/v4/letter/1/d78d45/32.png) [@11112](https://discourse.vtk.org/u/11112)
#### Post date: [August 12, 2020, 8:30am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955/3 "2020-08-12T08:30:32Z")

</div>

I’m sorry for my bad expression.  
I know how to slice actually.  
What I want to ask is “How to split imagedata into more object, I have plane data from normals and origins of normals”

---

<div class="post-metadata">

### Author: ![11112](https://discourse.vtk.org/letter_avatar_proxy/v4/letter/1/d78d45/32.png) [@11112](https://discourse.vtk.org/u/11112)
#### Post date: [August 13, 2020, 9:40am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955/4 "2020-08-13T09:40:33Z")

</div>

Hey Guys!!  
I am not clip my volume with two vtkPlane,  
And new question is coming, The origin of this volume still not change.  
That is weird to rotate my result.  
**Hot set origin to the center of the part I clipped??**

---

<div class="post-metadata">

### Author: ![flaviu2](https://discourse.vtk.org/letter_avatar_proxy/v4/letter/f/74df32/32.png) [@flaviu2](https://discourse.vtk.org/u/flaviu2)
#### Post date: [August 13, 2020, 9:50am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955/5 "2020-08-13T09:50:43Z")

</div>

With _vtkObject-\>SetOutputOrigin(…)_

P.S. See also this point: [vtkImageReslice Output Appears to be "Shifted Up and to the Left" by one Pixel](https://discourse.vtk.org/t/vtkimagereslice-output-appears-to-be-shifted-up-and-to-the-left-by-one-pixel/3884/4)

---

<div class="post-metadata">

### Author: ![11112](https://discourse.vtk.org/letter_avatar_proxy/v4/letter/1/d78d45/32.png) [@11112](https://discourse.vtk.org/u/11112)
#### Post date: [August 14, 2020, 10:01am UTC](https://discourse.vtk.org/t/how-split-image-data-with-plane-form-dicom-images/3955/6 "2020-08-14T10:01:54Z")

</div>

Thanks!  
But a new question comes, the slice I need is the area after clipping.  
I’m trying to get use position of imagedata between the area.  
then get the slice.  
Is there any better way??
