Python code with comments:

Check if the name of mm.Item_g is 'OR TP Lobe-R' or 'TemporalLobe_R'

if mm['Item_' + str(g)].name == 'OR TP Lobe-R' or mm['Item_' + str(g)].name == 'TemporalLobe_R': dcmOrigin = ak.ImagePositionPatient dcmSpacing = ak.PixelSpacing K = len(CT_Slice[qq].Organ['Item_' + str(g)].ContourData) / 3 x1 = [0] * (K + 1) y1 = [0] * (K + 1) for ii in range(K): x1[ii] = (mm['Item_' + str(g)].ContourData[1 + (ii - 1) * 3][0] - dcmOrigin[0][0]) / dcmSpacing[0][0] y1[ii] = (mm['Item_' + str(g)].ContourData[2 + (ii - 1) * 3][0] - dcmOrigin[1][0]) / dcmSpacing[1][0] x1[K] = x1[0] y1[K] = y1[0] fill(x1, y1, [2/255, 2/255, 2/255]) f = getframe(gcf) gg = f.cdata gg = rgb2gray(gg) a, b = find(gg) len = len(a) for j in range(len): ggsv[a[j]][b[j]] = 2/255 patch(ax, bx, 'k') hold on

Check if the name of mm.Item_g is 'OR Eye-L' or 'Eye_L'

if mm['Item_' + str(g)].name == 'OR Eye-L' or mm['Item_' + str(g)].name == 'Eye_L': dcmOrigin = ak.ImagePositionPatient dcmSpacing = ak.PixelSpacing K = len(CT_Slice[qq].Organ['Item_' + str(g)].ContourData) / 3 x1 = [0] * (K + 1) y1 = [0] * (K + 1) for ii in range(K): x1[ii] = (mm['Item_' + str(g)].ContourData[1 + (ii - 1) * 3][0] - dcmOrigin[0][0]) / dcmSpacing[0][0] y1[ii] = (mm['Item_' + str(g)].ContourData[2 + (ii - 1) * 3][0] - dcmOrigin[1][0]) / dcmSpacing[1][0] x1[K] = x1[0] y1[K] = y1[0] fill(x1, y1, [3/255, 3/255, 3/255]) f = getframe(gcf) gg = f.cdata gg = rgb2gray(gg) a, b = find(gg) len = len(a) for j in range(len): ggsv[a[j]][b[j]] = 3/255 patch(ax, bx, 'k') hold on