Center Unif ROI

center_comparison.py

import pylinac

ct = pylinac.CatPhan504.from_demo_images()
ct.analyze()
data = ct.results_data()
center_unif_roi = data.ctp486.rois['Center']
for name, roi in data.ctp486.rois.items():
    if name == 'Center':
        continue
    if abs(roi.value - center_unif_roi.value) > 30:
        print(f"ROI {name} was over the threshold of 30 HU")