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 }