# 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(