1 program roi0;
2 uses mtx09;
3 var
4 a, r: mtx;
5 begin
6 equ(r,'7.1,10.2;9.3,12.5;82.9,85.1'); // define ROIs by increasing brackets,
7 mwr(r); // overlap is ok
8 zer(a,101,2); //
9 lsp(col(a,1),0.0,100.0); // a(1)=0... a(101)=100 (i.e. 100 segments)
10 elm(col(a,2),col(a,1)); // a(:,2) = a(:,1)
11 dvd(col(a,2),100.0); // a(:,2) = a(:,2) / 100.0
12 roi(a,r); // select ROIs...
13 mwr(a); // a is trimmed of all data not in the ROIs
14 end.
15 {
16 [jroux@ami ex]$ ./roi0
17 8 0.0800000
18 9 0.0900000
19 10 0.100000
20 11 0.110000
21 12 0.120000
22 83 0.830000
23 84 0.840000
24 85 0.850000
25
26 [jroux@ami ex]$
27 }
28 {
29 7.10000 10.2000
30 9.30000 12.5000
31 82.9000 85.1000
32
33 8 0.0800000
34 9 0.0900000
35 10 0.100000
36 11 0.110000
37 12 0.120000
38 83 0.830000
39 84 0.840000
40 85 0.850000
41
42 }