#!/usr/bin/env python
# coding: utf-8

# In[1]:


import numpy as np
import matplotlib.pyplot as plt
get_ipython().run_line_magic('matplotlib', 'inline')

import dicom


# In[2]:


dcm = dicom.read_file("star.dcm", force=True)
dcm


# In[3]:


block_data = dcm.BeamSequence[0].BlockSequence[0].BlockData
block_data


# In[9]:


energy = float(dcm.BeamSequence[0].ControlPointSequence[0].NominalBeamEnergy)
energy


# In[14]:


ssd = float(dcm.BeamSequence[0].ControlPointSequence[0].SourceToSurfaceDistance) / 10
ssd


# In[17]:


applicator = dcm.BeamSequence[0].ApplicatorSequence[0].ApplicatorID
applicator


# In[5]:


x = np.array(block_data[0::2]).astype(float)
y = np.array(block_data[1::2]).astype(float)

plt.plot(x, y)
plt.axis("equal")


# In[ ]: