pixel array example

gistfile1.py

import os.path as osp
from pydicom import read_file
import numpy as np
import matplotlib.pyplot as plt
# use sample file; cd is in pydicom root folder
test_file = osp.join(osp.dirname(osp.abspath(__file__)), 'tests', 'test_files', 'CT_small.dcm')

dcm = read_file(test_file)

# show pixel array
plt.imshow(dcm.pixel_array)
plt.show()

new_array = np.ones(dcm.pixel_array.shape, dtype=dcm.pixel_array.dtype)
# try to replace pixel array
dcm.pixel_array = new_array
# will show no change
plt.imshow(dcm.pixel_array)
plt.show()