# visual checking of images and mask
ImSeg.display_overlay_volume(mri_array,
ser_mask_array,
'DCE-MRI image + SER mask',
aspect=ar)
if ii == 0:
# figure out the patient MRN
tag_mrn = dicom.tag.Tag(0x0010, 0x0020)
pt_mrn = dceSeries._DS[tag_mrn].value
tag_acc_num = dicom.tag.Tag(0x0008, 0x0050)
pt_acc_num = dceSeries._DS[tag_acc_num].value
# get VOI info
voi_ixyz, voi_size = ImSeg.GetVOIinfo(
CENTER, HALFLENGTHS, IG_mri_img)
# Get mask VOI (only need to do it once)
ser_mask_sagittal_roi = ImSeg.GetITKVOI(
ser_mask_sagittal, voi_size, voi_ixyz)
ser_mask_array = ITKImageHelper.itkImage_to_ndarray(
ser_mask_sagittal_roi)
# get the image VOI
mri_itk_sagittal_roi = ImSeg.GetITKVOI(mri_itk_sagittal,
voi_size, voi_ixyz)
mri_array = ITKImageHelper.itkImage_to_ndarray(
mri_itk_sagittal_roi)
# rescale the MRI VOI image
mri_roi_rescale_itk = ImSeg.ITKStandardizeImageIntensity(
