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

79 xlabel('Time [min]','fontsize',14);ylabel('Counts','fontsize',14)
80 axis([0 Time(end) 0 max(DTc99m)]);
81 set(get(h2,'CurrentAxes'),'fontsize',14)
82 legend('Uncorrected Data','Final data');
83
84 save eks cor Time STc99m
Listing A.2: load plot count.m
1 function C = load plot count(Data,TID)
2 % This function extracts the 4 images in one dicom−file and flips the
3 % posterior images. Data is corrected using anterior and posterior image.
4 % The output of this function is a vector containing counts from the two
5 % isotopes and the time of examination.
6 % Copyrigth Sabrina Wendt, DTU, 31/3−2009
7
8 Data 1 = Data.Images(1:end,1:end,1); % ANT Tc99m.
9 Data 2 = Data.Images(1:end,end:−1:1,2); % POST Tc99m, is mirrored
10 % compaired to ANT Tc99m.
11 Data 3 = Data.Images(1:end,1:end,3); % ANT In111.
12 Data 4 = Data.Images(1:end,end:−1:1,4); % POST In111, is mirrored
13 % compaired to ANT In111.
14
15 Data 1d = double(Data 1)+1; % Data is converted to double.
16 Data 2d = double(Data 2)+1; % Data is converted to double.
17 Data 3d = double(Data 3)+1; % Data is converted to double.
18 Data 4d = double(Data 4)+1; % Data is converted to double.
19
20 Data Tc99md = sqrt(Data 1d.*Data 2d); % Geometric correction.
21 Data In111d = sqrt(Data 3d.*Data 4d); % Geometric correction.
22
23 Data Tc99m = uint16(Data Tc99md−1); % Data is converted to uint16.
24 Data In111 = uint16(Data In111d−1); % Data is converted to uint16.
25
26 LP = [1 2 1;2 4 2;1 2 1]/16; % Defining a lowpas filter.
27 Data Tc99m LP = conv2(Data Tc99m,LP,'same'); % Filtering Tc99m data.
28
29 Borderline = find border(Data Tc99m LP); % Finding the area of the stomach.
30
31 S Tc99md = Borderline .* Data Tc99md; % Using counts from the stomach.
32 S In111d = Borderline .* Data In111d; % Using counts from the stomach.
33 S Tc99m = uint16(S Tc99md−1); % Data is converted to uint16.
34 S In111 = uint16(S In111d−1); % Data is converted to uint16.
35
36 c Tc = sum(sum(S Tc99m)); % Calculating total number of counts
37 % within the stomach.
38 c In = sum(sum(S In111)); % Calculating total number of counts
39 % within the stomach.
40
41 % Finding time for count instance.
42 t = Data.DicomHeader.ContentTime;
43 time sec = (t(1)*10 + t(2)*1)*60*60 + (t(3)*10 + t(4)*1)*60 + ...
44 (t(5)*10 + t(6)*1);
45 time min = time sec/60;
85