Matematisk Model for Mavesækkens Tømning - Danmarks Tekniske ...

46
86 MATLAB kode til kantfindingsprogram
47 minC = min(min(Data Tc99m)); % Calculating limits for the colormap
48 maxC = max(max(Data Tc99m));
49 minC2 = min(min(Data In111));
50 maxC2 = max(max(Data In111));
51
52 figure
53 subplot(2,2,1)
54 plot ven data(Data Tc99m) % Plotting raw Tc99m data.
55 title('Original Tc99m');
56 xlabel('Horizontal [mm]'); ylabel('Vertical [mm]');
57 subplot(2,2,2)
58 plot ven data(S Tc99m) % Plotting the stomach from Tc99m data.
59 title('Stomach Tc99m');
60 xlabel('Horizontal [mm]'); ylabel('Vertical [mm]');
61 caxis([minC,maxC]);
62 subplot(2,2,3)
63 plot ven data(Data In111) % Plotting raw In111 data.
64 title('Original In111');
65 xlabel('Horizontal [mm]'); ylabel('Vertical [mm]');
66 subplot(2,2,4)
67 plot ven data(S In111) % Plotting the stomach from In111 data.
68 title('Stomach In111');
69 xlabel('Horizontal [mm]'); ylabel('Vertical [mm]');
70 caxis([minC2,maxC2]);
71 print('−depsc', '−tiff', ['subplot tid' num2str(TID)])
72 print('−djpeg', '−loose', ['subplot tid' num2str(TID)])
73
74 C = [time min; c Tc; c In]; % Returning number of counts from Tc99m
75 % and In111 and time instance.
Listing A.3: find border.m
1 function Borderline = find border(Data)
2 % This function finds the position of the stomach.
3 % Data − A matrix containing a set of data from a stomach−examination
4 % with gamma−cameras.
5 % The call to the function is: find border(Data).
6 % Copyrigth Sabrina Wendt, DTU, 31/3−2009
7
8
9 si = 58; % Stomach Index, describes the lower bondary of the stomach.
10 cs = (1/15)*max(max(Data(1:si,:))); % Count Sensivity, depends on maximal
11 % number of counts in that image.
12
13 if cs == 0 % To awoid getting 0 counts when the Count Sensitivity is
14 % very low.
15 cs = 1;
16 end
17
18 LD = length(Data);
19 Borderline = zeros(LD);
20 for i = 1:1:si
21 for j = 1:1:LD