# How to extract multiple components image from dicom series?
**URL:** https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097
**Category:** Beginner Questions
**Created:** [July 20, 2024, 4:34am UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097 "2024-07-20T04:34:52Z")
**Posts on this page:** 6
**Page:** 1
### Author:  [@jinquan\_guan](https://discourse.itk.org/u/jinquan_guan)
#### Post date: [July 20, 2024, 4:34am UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097/1 "2024-07-20T04:34:53Z")
The dicom files can be correctly recognized as multiple components image by `ITK-SNAP`:
The proof is : I can change the component/view of the dicom file:
But when I extract nii.gz file from the dicom series, the views/components are concatenated together, just like:
Are there any problems of my python code? Could you help me to extract the nii.gz file from the dicom series by a correct way?
here is my python code:
```auto
def read_dicom(src_dir, out_dir):
if not os.path.exists(out_dir):
os.mkdir(out_dir)
PixelType = itk.ctype("signed short")
Dimension = 3
ImageType = itk.Image[PixelType, Dimension]
namesGenerator = itk.GDCMSeriesFileNames.New()
namesGenerator.SetUseSeriesDetails(True)
namesGenerator.SetGlobalWarningDisplay(False)
namesGenerator.SetDirectory(src_dir)
seriesUID = namesGenerator.GetSeriesUIDs()
if len(seriesUID) < 1:
print("No DICOMs in: " + src_dir)
sys.exit(1)
for uid in seriesUID:
seriesIdentifier = uid
fileNames = namesGenerator.GetFileNames(seriesIdentifier)
reader = itk.ImageSeriesReader[ImageType].New()
dicomIO = itk.GDCMImageIO.New()
dicomIO.LoadPrivateTagsOn()
reader.SetImageIO(dicomIO)
reader.SetFileNames(fileNames)
reader.ForceOrthogonalDirectionOff()
image = reader.GetOutput()
print("IMage Size: ", image.GetLargestPossibleRegion().GetSize())
writer = itk.ImageFileWriter[ImageType].New()
writer.SetFileName("./results.nii.gz")
writer.UseCompressionOn()
writer.SetInput(image)
writer.Update()
```
Thanks for your help!!! Please!!! 😊
---
### Author:  [@jinquan\_guan](https://discourse.itk.org/u/jinquan_guan)
#### Post date: [July 20, 2024, 7:54am UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097/2 "2024-07-20T07:54:15Z")
I found that `itk-snap` can recognize the number of volume/components.
---
### Author:  [@dchen](https://discourse.itk.org/u/dchen)
#### Post date: [July 20, 2024, 5:31pm UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097/3 "2024-07-20T17:31:20Z")
Does your PixelType need to be a Vector?
I’m not an ITK-Python guy, but here’s an example for creating an image of vectors:
[https://examples.itk.org/src/core/common/createanimageofvectors/documentation](https://examples.itk.org/src/core/common/createanimageofvectors/documentation)
---
### Author:  [@jinquan\_guan](https://discourse.itk.org/u/jinquan_guan)
#### Post date: [July 21, 2024, 5:02pm UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097/4 "2024-07-21T17:02:17Z")
Yes, you are right. Thanks for reminding me about the `PixelType`.
The `PixelType` need to be a Vector. But after I set as follow:
```auto
ComponentType = itk.ctype("float")
PixelType = itk.Vector[ComponentType, 2]
ImageType = itk.Image[PixelType, Dimension]
```
It still didn’t work. It does read pixel as vector, but the value of the vector is the same one.
---
### Author:  [@jinquan\_guan](https://discourse.itk.org/u/jinquan_guan)
#### Post date: [July 21, 2024, 5:10pm UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097/5 "2024-07-21T17:10:49Z")
Now, I think that my problem has become `how to read dicom as image of vector`. That’s right?
---
### Author:  [@matt.mccormick](https://discourse.itk.org/u/matt.mccormick)
#### Post date: [July 23, 2024, 2:53pm UTC](https://discourse.itk.org/t/how-to-extract-multiple-components-image-from-dicom-series/7097/6 "2024-07-23T14:53:31Z")
Hi @jinquan_guan ,
It may be worth trying
```auto
itk.imread
```
A directory with DICOM files can be passed in directly.
The pixel type is auto-detected by default. If that is not desired or possible, a desired type can passed. If series restrictions are needed, the `series_uid` kwarg can be set. More information can be found with `help(itk.imread)`.