Extracting 2d array from numpy.pixel_array dicom image

  • Python
  • Thread starter ProPatto16
  • Start date
  • #1
326
0
Hi guys,

I have a dicom image from a QC spect acquisition.
reading the file in with dicom.read_file and pulling pixel data using numpy.pixel_array gives array of shape [2, 1024, 1024]

The 2 represents two images, one from each detector head on the spect scanner and the 1024x1024 are the arrays i want to separate and pull out individually, i.e. so i can manipulate the pixel data, remove buffering 0's, plot histograms of the pixel values etc etc. i have artifacts in one detector head so i need to be able to separate the two slices

looking at numpy package there are ways to "slice" my array but the orientations are confusing and it doesnt seperate them into their own arrays.

This should be quite elementary i would of thought

Thanks
 

Answers and Replies

  • #2
326
0
so it looks like its this easy..

'array' is my array of shape [2,1024,1024]

to get the two 1024x1024 arrays out is simply

1stArray = array[0] and 2ndArray = array[1]
..
ill know for sure when i map them
 
  • #3
87
138
so it looks like its this easy..

'array' is my array of shape [2,1024,1024]

to get the two 1024x1024 arrays out is simply

1stArray = array[0] and 2ndArray = array[1]
..
Yes, it is correct though it doesn't look intuitive. :smile:
In C++, I would probably use i.e a tuple then tie its values to declared variables. Other than that I might use a struct of e.g 3 elements to store the returned values. You can do the same to make more sense of your current approach.
 

Related Threads on Extracting 2d array from numpy.pixel_array dicom image

  • Last Post
Replies
4
Views
774
Replies
7
Views
7K
  • Last Post
Replies
4
Views
14K
  • Last Post
Replies
8
Views
42K
Replies
12
Views
8K
  • Last Post
Replies
4
Views
5K
Replies
11
Views
1K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
1
Views
596
Replies
16
Views
4K
Top