p c #B2BBB8",
+",p c #9DA1A1",
+"'p c #A3A9A8",
+")p c #98A29E",
+"!p c #A5A8A8",
+"~p c #A8ACAC",
+"{p c #CED3D2",
+"]p c #D1DADB",
+"^p c #808A8C",
+"/p c #AFB5B2",
+"(p c #A1A9AD",
+"_p c #B2B7B5",
+":p c #B0B3B0",
+"q c #A4AEAB",
+",q c #969C9A",
+"'q c #454B49",
+")q c #8B9691",
+"!q c #7F8B8B",
+"~q c #9FA4A7",
+"{q c #AFB2B1",
+"]q c #ABB0AF",
+"^q c #C8CED0",
+"/q c #939999",
+"(q c #44484B",
+"_q c #C2C5C6",
+":q c #59595B",
+"r c #525656",
+",r c #485052",
+"'r c #858587",
+")r c #9EA39D",
+"!r c #313536",
+"~r c #515257",
+"{r c #6C7271",
+"]r c #70767A",
+"^r c #1B1A1C",
+"/r c #A2A5A6",
+"(r c #666A6C",
+"_r c #919696",
+":r c #485051",
+"s c #B2B9B8",
+",s c #56595C",
+"'s c #57585D",
+")s c #D8D8D9",
+"!s c #4F4F54",
+"~s c #999FA0",
+"{s c #9FA5A1",
+"]s c #737475",
+"^s c #727575",
+"/s c #8E9394",
+"(s c #8B888B",
+"_s c #292E2F",
+":s c #A7B2B2",
+"t c #8B9591",
+",t c #9EA3A7",
+"'t c #46464A",
+")t c #858D91",
+"!t c #C7CBCF",
+"~t c #ABB1AE",
+"{t c #ADB3B2",
+"]t c #ACB1B1",
+"^t c #A5A9AB",
+"/t c #D7DFDC",
+"(t c #757776",
+"_t c #000003",
+":t c #030304",
+"u c #0A0808",
+",u c #303231",
+"'u c #D8DEDF",
+")u c #A5ABB1",
+"!u c #B6BCB7",
+"~u c #B3C0BE",
+"{u c #B4BDBD",
+"]u c #B0BABA",
+"^u c #B4BABD",
+"/u c #ACB5B7",
+"(u c #ABB2AF",
+"_u c #AAAEAF",
+":u c #AEAEB1",
+"v c #A5AFAD",
+",v c #A2ABAB",
+"'v c #A2AEAB",
+")v c #A5ADA8",
+"!v c #A2A9AB",
+"~v c #A0A9AA",
+"{v c #A4ABA9",
+"]v c #A6AAA8",
+"^v c #A0ACAB",
+"/v c #A3A9AB",
+"(v c #A4A8A9",
+"_v c #A0A6A4",
+":v c #A2ACA6",
+"w c #AEB2B3",
+",w c #A9B1B2",
+"'w c #ADAFB0",
+")w c #CDD3D3",
+"!w c #010401",
+"~w c #0C0C0A",
+"{w c #464B48",
+"]w c #414544",
+"^w c #4D4F50",
+"/w c #818789",
+"(w c #A1A4A4",
+"_w c #E3E7E6",
+":w c #E8EDE9",
+"x c #ACAFAF",
+",x c #AAB2AD",
+"'x c #ADB3B1",
+")x c #ABB1AF",
+"!x c #DEE3E2",
+"~x c #313432",
+"{x c #010300",
+"]x c #020401",
+"^x c #191B19",
+"/x c #696A6A",
+"(x c #656260",
+"_x c #5D5D5C",
+":x c #545355",
+"y c #BBC3BF",
+",y c #BAC0C2",
+"'y c #B1B9BA",
+")y c #A7AEAD",
+"!y c #A1ACAC",
+"~y c #A5A7A9",
+"{y c #A3A7A8",
+"]y c #A2AAA7",
+"^y c #A3A7A5",
+"/y c #A1A5A3",
+"(y c #A1A7A4",
+"_y c #A2A7A5",
+":y c #A4A9AB",
+"z c #ABAFAE",
+",z c #ABB3AF",
+"'z c #AAB2AF",
+")z c #E6EAE8",
+"!z c #52514E",
+"~z c #000900",
+"{z c #0F110A",
+"]z c #090C04",
+"^z c #D7D7D7",
+"/z c #B3B3B3",
+"(z c #AAAAAA",
+"_z c #A3A3A3",
+":z c #7C7C7C",
+"A c #9CA2A0",
+",A c #9BA2A0",
+"'A c #9FA4A3",
+")A c #9FA4A2",
+"!A c #9BA0A0",
+"~A c #9AA3A2",
+"{A c #9CA6A4",
+"]A c #9AA4A2",
+"^A c #9EA7A4",
+"/A c #9EA8A2",
+"(A c #A2AAA5",
+"_A c #A2A7A9",
+":A c #A0A9A6",
+"B c #030204",
+",B c #070601",
+"'B c #050600",
+")B c #050602",
+"!B c #030700",
+"~B c #DADDDB",
+"{B c #D5D5D5",
+"]B c #BEBDBE",
+"^B c #B8B8B8",
+"/B c #A9A9A9",
+"(B c #919191",
+"_B c #818182",
+":B c #7A7B7B",
+"C c #3B393B",
+",C c #3B3C3A",
+"'C c #464644",
+")C c #565757",
+"!C c #656566",
+"~C c #666866",
+"{C c #6A6C6A",
+"]C c #737577",
+"^C c #808483",
+"/C c #959694",
+"(C c #97A09D",
+"_C c #9C9E9C",
+":C c #9EA6A2",
+"D c #090F08",
+",D c #B5B5B5",
+"'D c #AEAEAE",
+")D c #A8A9A8",
+"!D c #A0A0A0",
+"~D c #8E8E8E",
+"{D c #848584",
+"]D c #7D7D7D",
+"^D c #6C6C6D",
+"/D c #585959",
+"(D c #535253",
+"_D c #4B4D4C",
+":D c #484848",
+"E c #AEB5B2",
+",E c #7D7D7B",
+"'E c #555659",
+")E c #3D3F40",
+"!E c #313232",
+"~E c #2A2B2C",
+"{E c #2A2A2B",
+"]E c #BCBDBD",
+"^E c #B7B6B6",
+"/E c #B2B1B2",
+"(E c #ACACAC",
+"_E c #A4A4A4",
+":E c #8F8E8E",
+"F c #757675",
+",F c #6E6D6E",
+"'F c #646565",
+")F c #606061",
+"!F c #5C5B5B",
+"~F c #595A59",
+"{F c #545555",
+"]F c #4F4F4F",
+"^F c #454544",
+"/F c #414140",
+"(F c #3E4040",
+"_F c #3B3939",
+":F c #3C3838",
+" , ' ) ! ~ { ~ ~ ] ^ / ( _ ) ) ) ) : < ^ [ } ~ { | 1 | ( ) = * 2 * { @ @ . . . . . . . . . . . . . @ . @ . . . . . . @ * . . . . . . . . . . . . . . . . . . . ",
+". @ @ @ . @ . @ + . . . . . . . . . . . . . . . . . . . % + = * = ~ 3 4 4 5 6 7 8 9 9 9 0 a b c d c e f g h i j j j k l m n f o p q r s t u v w x y l z f k A B n C D E F G H I 9 J ) { { { * = @ . . . . . . . . . @ @ . @ @ . . . . . . . . . . . . . . . . . . . ",
+". . . K L + @ @ + @ @ + & M N O P Q P R S T U V W X w v A Y Z ` . ...+. . .+.+.@.m m m m v #.l h n $.l $.W %.j &.@.n #.h *.=.m -.-.;.>.W y v ,.'.-.).!.~.{.k ].z l z z ^.k &.&./.(._.(.:.(./.<.[.}.|.1.2.3.4.5.= @ + @ @ @ . . . @ . @ . . . . . . . . . . . . . . ",
+"6.7.8.9.0.d $.a.'.m .+.p h b.c.d.e.f. .h .v g.h.i.j.k.*.l.w n m.$.n.x -.o.0.'.p.q.r.l.s.b ^.F }.c F t.u.v.w.F x.y.X z.d w d A.|.B.|.}.C.D.V E.X F.B.G.F.H.I.v.J.w K.E ^.A.}.}.p.n.y.v.L.|.B.v.M.N.v.v./.&.O.P.[.N./.}.Q.R.S.T.+ . . . . . . . . . @ . . . . . . . ",
+"` U...V.W.h X.Y.Z.W y $.v #.`.F c +.+e ++@+#+$+%+B.v.V #+d &+X p.p.*+d =+C.-+;+>+a ,+K.X '+V )+K.C !+|.~+K.!+{+V )+|.{+]+^+d /+K.(+;+|.V _+:+)+V %+<+^+F.|.[+}+a |+)+1+2+[+-+3+4+E.$+&+5+K.[+6+(+;+G.(+{+1+#+7+}.G.8+v.9+[.0+9+Q.a+9 b+* . . . . @ @ . . @ . . . . ",
+"c+d d+@+V &+d e+f+0.'+X /+g+c+h+,+i+^+A.j+$+%+k+l+m+,+{+{+[+n+o+o+p+9.q+r+s+t+p+u+v+w+$+Q.l+x+y+z+A+B+/+^+/+{+C+D+E+F+G+P H+I+J+K+L+M+N+O+P+Q+R+S+T+U+I+V+P W+X+L+L+Y+Z+`+Z+ @.@R +@@@#@$@%@I+&@ @Z+*@=@-@1.q+9.@@3+;@1.3+9.5+%+>@N.,@M.'@)@!@: . . ~@. + @ . . @ . ",
+";+V g+V {@]@K.^+/+^@^+^@R /@(@_@$@^+:@<@[@}@H+5+E+_@|@R ^@1@(@-@A+v+2@K+3@H+4@v+5@M+6@7@8@I+9@6@0@a@b@c@d@e@b@b@f@g@8@0@0@h@0@i@j@k@8.l@=@H+i+m@T n@o@(+p@q@'@r@s@t@5+S ^@E+8.u@A+v@r@w@X+;@X+x@y@z@A@R h@B@C@8@D@E@P C@q@8.F@@@2.k@1.t+Q.Q.}.v.|.G@. $ @ . . . . . ",
+"R H@S+/+p+/+I@J@(@(@(@(@S #@Z+E+K@y+W+L@M@N@O@P@Q@R@J@S@4.T@T@U@U@V@W@X@Y@Z@`@ #=@y@.#+#@#2.B@B@T+##$+$+1+Q.t+$#%#%+&#+*#-+}.=#X+-#9.E.B+r@.#t@;#1.E+A@7@;#8@P@>#O@,#'#)#!#~#{#7@O P ]#~#J@a+(@k@^#/#(#a+_#(#`@g@]#:#<#[#}#7@|#5+O+E+1#Q@<#3.B.:.-+2#. @ . . . . ",
+"x+3#4#^@u@5#'#6#7#8#7.R@8#9#j@{#0#C+}#j@(#a#b#-@R@g@P 1.R /+E.$###>+c#d#^+##e#f#Q.g#$+K.)+h#i#j#'@j#C@<#o+k#k@C+]#'#l#R C+{#0#^#m#k@n#o#S@a+p#A@|#q#r#s#0@]#t#]#E+8#U+u#v#t#m#w#a+v#x#y#z#A#R@R@P B#C#D#Q 7.S@E#h@!#O+3.F#^#8@G ;#k@^#[#G#8.(#|#H#@@I#I#J . . . . . ",
+"G J#K#{#L#M#R@C@3.z@9.-@l+l+e+@+*#E.d#}+##k+N#&+j+g+i#{+{+/+(@ #|#0@0@d@M@K#O#P#Q#R#S#9 T#U#V#6#Q G {#W#X#Y#8#P#K#Q K#T@J@Z#`# $.$O@+$@$'##$Z#$$C@%$&$R@8#O S@G d@*$4.a+Z#7.g@e@9 =$-$;$|#>$S@,$Y@4.0@'${#)$s#S@J#4.S@>#x#R.R.R.R.!$~$M#)#7.s#Q G {$}./.]$C@. . . . ",
+"Z+o+l+^$/$p.p.&+&+m.($_$($3+l+:$<#0@v#Y#<$H#z#y#[$}$|$S#1$2$3$4$5$.$d@c@I >$Y#6$>#7.x#P 7.x#7$@$S# $)#P#d@S#b@U@e@K#8$R#[$x#>#-$%$Z#9$@$u#0$%$K#X@a$`#8$4$R.b$)$c$s#9#4.9 P#*$d$d@e$%$a$f$!@g$@$e@b@8$Y@@$)$a@b@0 c@h$c@7.S#S#i$S#j$k$P#7.[$S#l$I 8$X@5+/.m$3.K @ . ",
+"C n$@@=@O+r#s#z#a$k$N o$p$q$r$s$j$4.t$k$>$5$j$8#9#0@d@9@-$C#c$d@<$d@O#e@e@u$>$%$%$e@@$Y@>$z#e@z#[$c@v$y#w$*$t$)@x$7.Y@y$[$z$L#y#z#0 J#k$*$A$B$N )$0 4$b@C$D$U@%$E$F$G$b$e$a@g$a@8$j$0$d@t$b$H$4$8$c$I$J$b$y#Y@I K$L$I 0 5.9#9 %$9 %$M$N$f@o$O$J$J$S#P$O+t+4+]$Q$. . ",
+"R$S$g$0 %$S#4$S#P#[$[$P#N S#a$K#8$T$r#U$M#I$8$[$9$9$a@y#~$d@e$V$B#S#[$>$*$5$b$8$k$y#8$4$[$7$P#;$t$f@8$k$y#W$Y@X$0 e$)$)$Y@9$I O$e@Y$Y$J$J$t$g$b$L#!@v$[$0 a@Z$5.%$0 J$)@5.`$O$N u$!@0 0 %.%4$+%@%#%P$o$N 0 $%%%&%*%=%0 -%J$;%>%;%,%'%o$)%O$O$!%*%e$1$o#'+~%{%]%R . ",
+"k$M#X@^%M#S#[$[$[$S#8#k$S#N 8$*$8 5$I !@X$P#8$8$!@0 @$z#/%(%0 M#_%g$:%9 I 9$5.g$B$<%N [%0 e$0 P$I I : y#N g$g$N$W$!@e$+%X$X$N J$J$p$p$q$!@}%2$|%9$)@t$l$)@R$I .%l$S#g$4$N$I 8$S#N B$O#%$*$a@G d@%$f@H$M#1%H$1%J#{$b@b@c@2%e$g$*$8#k$*$3%d@X$1%4%c@y#B$5%n+[+6%7%8%e$",
+"5.O$O$S#[$[$5.9%4./%0%b$J$J$8$9#a%)@X$b%c%q$o$1%N 9%P$0 >$}$0 d%e%f%1%)$g%5.O$h%i%O$6.: y#v$j%o$8 H k%d%)%l%0 m%n%5$2%l$i%o%S$l$p%i$O$Y$H$[%e@~$4$@$b$X$q%}$L$8$3$m%r%s%t%u%v%w%x%y%& z%+ A%B%C%D%A%A%A%& + + + + E%> F%& + + ~@+ ~@A%D%G%H%I%I%J%B%. D$2+K%#+L%M%N%",
+"5.N N N N [$%$5.I I O$O$P$!@O%g$D#O$K$P%Q%O$N N R%S%S%T%S%U%V%h%4$h%h%i%W%|$X%d%y$0%O$F$N o%[$g%N o$p%8$O$N y$Y%u$.%H q$[$s$Z%T%} 6.`% &.&+&@&. . . . K K 2##&= $&%&0%O$!@V$d@e$C@/#1.g+&&/$c#B.*&=&-&;&>&,&'&)&)&!&~&{&]&^&/&(&(&'&>&,&>@I.8+G.}.(&@+] 7 '@##_&:&<&",
+"u$6.6.y$o$J$)%[&6.)%6.i%[$[&[&d%H }&h%o$.%Z%|&.%1&6.o$2&2&)%3&4&7 F$5&X%6&7 6&}$T%T%7&8&9&] &0&M } a&b&c&d&e&f&. . . . . . g&T./%!@Z@h&i&3+##;&:&i j&k&l&m&n&o&p&q&r&s&t&u&v&w&x&y&z&A&B&C&D&/&E&F&G&H&A&I&J&K&L&M&N&O&P&Q&R&S&T&T&Y.P&U&V&W&X&V&Y&n&Z&_#`&#+ *d+.*",
+"N y$6.6.T%o$o$f%6.7&`%+*@*#*i$}&4&6.$*y$X%f%8 %*[%N )%o$Z%&***=*-*;*>*,*'*)** . # . . . . . K K !*~*{*)$r#=@]*^*/*(*_*:*<*m&[*}*|*1*2*3*4*5*6*7*8*9*0*a*b*c*d*e*f*g*h*i*j*k*l*m*U%n*o*p*q*r*s*t*u*v*w*x*y*z*A*B*C*C*C*D*E*F*z*G*H*I*J*K*L*M*N*O*P*Q*R*($S*J.T*U*V*W*",
+"6.p$6.} 6.R%6.X*Y*2&g$L$5$5$Z*`*A$ =.=} +=^ '** #&* @ @ . . . . . . K @=#=$=(#'@V*7%%=&=N%*===-=;=>=,='=)=!=~={=~=]=^=/=(=_=:=<=[=}=|=1=2=}=q&3=4=5=>=6=7=8=9=0=a=b=c=d=e=f=g=. h=i=j=i=i=k=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=i=l=m=n=o=@@=%p=q=",
+"`%+*Z%r=s=4$q$y$#=R%} t=u=v=w=+ ~@@ x=K . g&y=z=A=a@=@($M%B=C=D=E=m&F=G=H=I=J=2=/=K=L=M=N=O=P=>=Q=R=S=T=U=V=W=X=Y=Z=`= -.-+-S=@-#-;=3=3=;=$-%-&->=*-=---;-'=>-,-'-)-!-~-{-]-^-/-i=(-_-:-<-[-}-|-1-2-3-4-5-6-7-8-9-0-a-b-5-c-d-e-f-g-h-i-j-k-l-m-n-o-p-q-r-s-t-u-v-*.",
+"5 = + $ . . @ A%w-x-^ = T.: 1.=&^&(*y-z-A-B-C-D-E-F-G-H-I-%-J-H-K-K-L-M-N-O-P-Q-R--=S-T-U-V-W-W-X-Y-Z-S-P-`- ;}*.;+;@;J-#;$;%;&;*;=;-;;;>;,;';);!;|=!;!;~;=-{;];^;/;(;_;:;<;[;i=};|;1;2;3;4;5;6;7;8;9;0;a;b;c;d;e;f;g;h;i;j;k;l;m;n;o;p;q;i;j;r;j;s;t;u;v;w;x;y; #^*",
+"k%H I d@O c#z;A;B;C;D;E;F;G;H;I;J;K;L;E;M;N;O;R-P;Q;R;S;T;S;U;V;W;X;Y;Z;Q-O-N;`; >O;.>+>@>g=#>g=$>M-%>&>*>K;=>->`-;>>>,>'>)>&->>}=!>~>{>%-|=]>^>d=/>*-(>_>:><><=[>}>|>1>2>3>i=4>5>6>7>8>9>0>a>b>c>d>e>f>g>h>i>j>k>l>m>n>o>p>q>r>s>t>u>v>w>a;x>y>z>8>A>B>C>D>E>F>G>F+",
+"H>@>I>@;H-&;U=J>K>L>M>N>O>P>Q>R>S>T>U>V>W>X>Y>Z>`> ,H>.,+,@,L>J>R--=#,R-$,J>J>S-%,&,*,=,-,K>#-%,I>L-K-;,@;>,,,->@;',=;),R=!,~,{,Q=}=],^,;;>=/,>-(,}>M=~>~>_,:,<,[,},|,x=1,i=2,3,4,5,6,7,9>8,9,0,a,b,c,9;k>d,i>e,f,g,h,i,z>j,k,8>h;l,p>r;m,n,o,r;p,q,r,s,t,u,v,w,x,^&",
+"O>R;y,z,A,B,Y>C,A,C,X>Z;D,Y>E,F,y,S;m&G,H,I,J,K,L,M,U;N,O,P,Q,R,M>-=S,R,T,U,V,S-W,K;X,Y,Z,N-`,,>W= ' -.'`-`-';+'*;@'H=#'~,;;$'~,>>%'+'&'Q=*'6=|=a=='-'N=;'%'>',''')'!'~'i={']'^'/'('_':'0;<'r;r>o,0;k,['}'o>|'1'2'm>r,3'4'5'6'7'8'p,k>9'8>0'a'b'c'p,d'e'f'g'i=h'6%i'",
+"j'k'l'A,C,y,X>Y>X>X>m'n'o'T>p'*=q'`>r's'Y-t'u'v'N,w'L,R,P-O,D;P;R,J>R--,x'y'z'A'F-@>P-W,W=%>B'C'`-F-D'$;;;}*E'F'4=G';;H'->G-I'J'K'3=[=K'L'M'N'<=}>0=O'P'Q'R'S'!>T'U'V'W'X'Y'Z'`' ).)+)@)#)$)s;%)q,&)2'*)=)@)-)8';)b->),)')))!)~){)9;h;])}'^)0'n,d'f,/)&)()_):)<)[)})",
+"|)1)2)3)4)5)6)7)8)9) ,L,M,P;L>v'u'S-T-0)Q;a)b)M,c)d)e)f)g)#,h)W-$>X-i)j)J;K>&,k),>l)+>`-`,S-m)L-n);=*>o)=;+'p)q)'>F;r)s)n)%'$-*-s),;&'*;6=|*2=M'I'6=K=>;t)u)!;v)w)x)y)z)A)%)%)B)$)C)^'n>D)E)E)F)r>G)j,H)I)J)e;K)L)M)N)O)P)Q)R)S)T)i>f,U)7'V)V)W)j,X)z>0;Y)w><-i=Z)`)",
+"p' !.!+!@!#!F=M>$!%!&!*!L>=!L>-!;!>!@,#,,!'!)!!!~!{!i)P-R,X-Z-]!W-g=+>^!/!(!%,_!:!L-&>`;.>&;,>>1!I>H-),*;H-{>*>L'>=>=2!,;O'3!%'%'4!{>N'|*<,5!/>6!7!p%8!i=9!0!a!f;b!7>>):':'j,J)c!}'g>d!a'g;e!e!f!g!h!i!j!k!l!m!n!7'o!p!q!7'r!s!t!u!v!w!^)v!x!y!z!A!B!C!",
+"D!W;E!D;&!F!A-G!H!I!J!D;T;K!P;L!M!g)X;'!N!;!O!P!Q!R!S!T!=,U!K;R,W=V!K>J>J>&>}*U=W=W!+>J-',;=X!H;Y!N-Z!`!>,*;); ~.~+~@~[=}=H-#~$'$~2==;a=%~&~*~(,=~2=-~;~$->~,~'~)~!~~~{~b,]~5;Y)^~^)^'/~o,(~0'_~r,:~:'<~[~}~|~1~2~3~4~()5~6~7~8~9~T)0~Z'a~k>w>n;=)b~c~d~9;e~f~g~i=h~",
+"i~j~k~l~0)m~n~+,P;V-4)o~@,p~I,q~r~N,s~o~t~u~-!v~w~x~y~z'S-U,$>#,V!N-E-`-=> -E;z~W!K-&>@;o)I>A~H-)>B~;>I-=;;>o));$'C~5=>=1=%-G-L'!>D~E~!=F~G~H~I~e=:,J~q&K~L~M~N~O~k=P~Q~J)s;0'R~S~S~T~1'U~%)V~9;r;U~W~X~Y~Z~`~ {.{i=+{+{@{#{${%{&{9-*{h,={g,-{y>;{>{,{:''{#)m,){i=!{",
+"~{>!{{]{U;K,^{/{=!({{{p~U-P;_{$,]!N!:{R,]!V,<{[{P,P-W-B'}{-,g=W,U! -+>F-|{B'#,1{V!`;2{`-',Z&3{*;2{4{5{3=.;6{7{)>o)@;,;'>@;L'8{q&D~;;,;9{r)0{#~a{{>b{c{d{e{f{g{h{i{X).)V~j{k{l{/)g>a>m{n{o{p{q{i>&)r{s{t{u{v{w{x{y{g>z{e>A{P~o,B{C{0'D{k>E{Y)F{G{H{I{J{p{K{8>L{g'M{N{",
+"O{Q;P{Q{-=R{P;/{Y-S{T{$,V-P;R-U!U{J>%,O;V{W{X{Y{V-Z{`{`,A' ].]W=+-+]_!U,Z, -X!G=',I>@]E- -=;#]:!H=I>$]%]&]*] ~H-=][=G-H'5{{>,;3=3=I',;b{Q=-]r);]),Q=>],],;']Z#)]i=!]~]{]r>]]#)&)^]_~a'V)a' )0~9;0;/](]_]:]<]J)])g>c,[]}]q;}]Y)2'L{|]1]2]r;3]s;4]5]E{E{9;6]&)c,7]u,8]",
+"9]p~g)H>0]a]D;b]c]w'p~@,v~-,R-R-d]e]N!f]j)W{P,`{5=g]h]i]#,X=P-N-+-',j]X=N-%],>]!o)I>J-[!4=#;m)k].>l]m]s)8=n]*'I-],o]@;p]),),>;3=>=H=*;L'5=a{);q]c{r]]>=;s]t]u]v]w]X'+)x](~b~]]y]9,V)s>z]A]b;U)['5;B]C]D]E]F]G]H]9;E)I]J]K]i>}'L]s;M]N]O]a!P]E{c'g,h,b;$)q{&)g,Q]R]S]",
+"Q;u'S{T]u'U]V-#,V]J!_{W]X]Y]=,R--=Z]S!`]N,S- ^.^+^`;@^y'%,W=K-U,X!,>.]Y=.]W=-;I>#^$^%^L-n];;.]G'.;}=&^8=3=}=*^n]=^-^~,;^>^,^}={>>;'^m)3=!;)^>;$'4{+~H=!^~^{^]^i=(~^^:'o,a>/^(^j,():'b,_^:^<^1'y]]~6'[^}^y]|^m;a;s>1^b;9;b,p,K{c!m,2^3^4^5^6^=)B{=)7^8^9^0^a^3]b^c^d^",
+"e^f^a]g^h^i^j^_{S{W-:!V,b)k^l^V]m^@-@>V-n^,>U!I>N-o^.]p^',N--, -q^P-M-r^q^s^t^u^v^+-F;.;w^I;n]';x^3{k)>>;;>>&;-^q&B~5=;;H-y^),Y!->!,(,!;+~E~ ~~>%',;%'+'z^O A^i=B^t{h;/~p{C^D)f;^)r,j,*)^'D^g,p{E^F^G^H^I^l,J^])K^a~}]4]L^j,M^S)S~c!N^O^P^Q^R^S^j,R~8^T^U^['V^W^X^Y^",
+"Z]Z^Q-`^!!Y-R,s'R-J>J>P-W=w~ /S,./+/@/P-q^.]#>T]y'1>T,G=k)#/$/%/&/*/=/-/;/R=*>,>',M-H;5=D'#^%;(,*,H;>/';-/J-';}=}=2{H-&;n)V=,/],->'/)/!/-/{>%-->E'&;%-%'~/{&{/]/^///6,n{g,(/L^6,m,_/j,6,:'5'o,:/E)r,1/0'h,2/3/4/5;B{:'E{5/6/7/8/d'9/G)q!>{0/q>a/K^U)b/c^c/",
+"H>R!d/J>e/f/J>#,R-s'K;#,b)g/l)K>h/^!i/J>-,w^j/:!k^#;Z,q^W=T=k/l/+-1=p'm/),L-n/@]%]y^>>K>;>),'/B'G;o/J-p/q/o)M-@'-^v^%]F;S=V=m)#^J-;=6{r/J;1{3=3=;=I-@;N=s/t/i=u/v/S)g;$)w/6,x/6^<^H)9;y/z/~)@)A/y>B/C/D/f,c,z/J)[]c,E/V)0'F/K{G/H/9;I/J/K/L/C)Z'A/H/h,_~M/u!%)('N)N/",
+"U!J!g)T-V-V-Q,i)#,R,N-W=K>J>W-W,S-O/A'`,P/&,Q/R/A'S/<{+>+>T/#]'/$~U{],U/t~9{=^;;V/1>s^[!A'I>D'W/X/@'$^ ~Y/[!Z/%]`/p/H'w^2{>>#;>>],&/Z!4= -],&; (I-;=#^.(+(@(i=#($(M/u;L]:^m>g;%(*)g>E{5'h>&(S)q!w>^^&(*(h>G)J]=(>{-({]5;]~p,/)A>;(>)>(,(9>w>'(9;V^H/E{i>)(c'Z'!(k=~(",
+"V-{(Q,:{g=j/](W-K;N-]!J>]!N;$>-,E;P-^(k/:!Z,W=K-@-K-K>!,%-#;/(((_(:(<([(}(:&B 9{%-;/[! ~U=*;),,^|(1(2{5{@'H-&];=&/3=`-H'.;->D'%>%;2(D'>>n))>S=.;k)L'5=s]3(4(i=5(6(R~5;7(s;h;6^q{L^^)5;8(p!8'E{_~9(0(a(b(c(J)d'G/r!d(_/s;c!e(f(|/g(h(i(i>j(k(z]l(O^m(N]c>Z'n(e;o,i=o(",
+"p(o&$,q(K;p~N-S/T-b)b),>V]W=Z,+'4{r(#]s(t(u(#~v(S!,>u)w(x(y(z(A(V,X+2=B(C(D(E(F(`;#;&;*;->G(=]o)R=4=H(U=~,Z!%,2=';`;G=F;;=J-F;G(+'*,p/*;J- -H;J-',',G=,-I(J(K(L(M(o>c(@)N(p,5'm>1'g,c!%)O(8'(/P(|/c'o;z/Q(1]R(1]S(K{T(U(V(I]W(V)=)X(Y(Z(*{`( _K{t!._+_K]U^@_#_$_%_&_",
+"O;V,*_N-U-e]<{=_Q!<{-_f];_>_`,,_'_)_+~!_~_{_j)]_^_a=/_(___:_<_[_}_|_1_2_3_4_5_6_ ,H=@;[!7_G-`-W!8_H(J-X!9_H-0_a_Z&*;Y!b_c_d_e_f_',p/*;q^I;#~@'V!@;g_q/M'M#h_i_<-k>(/h,])j_k_k(c'()|/B{^)u>v!$)l_|^m_u;n_o_={p_t!p_t!q_r_s_V)t_u_o,v_w_x_i>B{6]y_T(z_u_A_B_C_D_E_F_u@",
+"'!G_H_P-.-I_C'J_P-i)K_j/I_3=L_M_N_O_P_Q_R_S_T_U_O>V_W_X_Y_Z_`_ :.:+:@:#:$:%:&:*:,>>>L-`->>#'O/U==:Z!-:L-;:2(@'X!>:x^&;b_ ~,:X!@'G-':%] ~#-):5{K>!:!:',(>~:{:]:^:/:8;(:p,m,j,T)R^K{w/a;J{a;G/V)m;_:::V)<:[:}:r!}:t_1]|:1:s!2:K{3:4:5:6:7:8:q!g~ _u_9:0:a:b:c:d:i=e:f:",
+"P-g:h:i)K-`;j)Z-v^i:K>9_j:;.k:l:m:n:o:p:q:r:s:t:u:P$v:w:x:y:z:A:B:. C:D:E:F:G:H:I:!,y^F-`-%];;y^J:n]G=G=H;W!{,W!',#> -p/ 'K:%]+>.;!:4=w^X!k)L:`-F- ~M:N:,$O:P:Q:R:1^]~0;d'J{()~)M(D^S:o!&(M(z_T:o_U:V:W:X:Y:Z:`: <.<+<@<#<$<%/M-j<5{`-J-,>k<):=^!:v^J;`,l<@;Y/%>k)&/.>H;F;!:H(k)&/-,;>k/N-,>q^muD1]7^]~X(EZ/L=[5=jW=I>1{]<>,-[K:;[r^ -.>I>+>'/>[:!&>I>I>*,,[K-Y/#,i)H;r^=>'[)[![i=B^&(a:~[R^-)5'_/sL-h[i[f j[k[_ ]%p=l[m[n[o[p[q[r[s[t[u[v[w[x[y[z[A[B[C[D[E[F[1{I> ;G[R=K-V=D'K>H[L;N-I>#;Y/H;&>w^j<+-I[V=.>]<#,J[S/%>':K-I>I>K[g[R,b{L[M[i=N[n>v!L{$)a;p{O[R^7,5/C)P[6^3'Q[|^R[}:S[(/T[U[R~V[W[X[Y[Z[`[ }8~.}'(g,+}a(H/@}1]#}Q( <$}%}&}i=*}. =}-}",
+"-,+-1{~/,>;}W!l)#>#;N-l<>}^{,}'})}!}~}{}]}^}/}(}_}:}<}[}}}|}1}2}3}. 4}5}6}7}8}9}0}=:`;*;a}V!;,b}c}4=d}U=.]+>_!@'d/i:l)m]=>U!e}&>f},>g}-,K;M:]!:!`;u~b)5=h}i}i=j}k}C)S~^)e~c>6,&(l}O[R^E{B_c~a(m}B{n}o}s,p{p}q}r}s}t}u}v}w}(/x}y}z}A}B}C}D}k,u_E}F}G}H}I}J}K}L}M}N}O}",
+"P}w^K-g_`-Y/`;F-q^.]*;Q}R}S}T}U}V}W}X}Y}Z}`} |.|+|@|#|$|;.%|&|*|=|-|;|>|,|'|)|!|'/I>;;5=K> -U=K;F;Z,U,K> -H;.>&>V!U,;=`,~/|{N-[!P-S,]!#,T,b_~|N,S-r^g={|]|>&^|/|(|R^*)g,m,6>(/c~_|5':|<|T)f!G)[|v>}|_~-)||o>j,g>1|2|3|K(4|Q~J{:|V^5|6|u!L]7|8|p_a~G}9|}]i=0|a|b|c|d|",
+"1>X!.>.>@]-,e|K-1{y',>,>R,f|g|h| |i|j|k|l|m|n|o|p|3{q|r|s|t|u|L,v|Q=.>w|x|y|%^V/+-c_e|K-;=G[:!I>;=W,z|`-U=J:+>.>W,=,g=U,+-X=+>P-J>+-Q-]! ]:!P-`;P- ]r^O,|!Y*A|B|j-%)-)x>i>C|D|}'j;0'E|@)x/V(F|A_e,f!U)V^G|L^H|~[I|J|K|L|B_M|N|u!O|P|0^Q|R|S|A_T|U|V|W|i=X|. Y|Z|`| 1",
+"W!G=Z{.1`--^1>1> -*;Z-S!/(+1@1#1$1%1&1*1=1-1;1>1,1'1~_)1W/V-]<;;H(q^u^!1~1{1Z&]1,:;/f]H;g=D';=]!J>U,^1,>N-N,S-/1V!,>',x'A'V!B'P-+-W-N-Y-g=#,W-N;]-R,]!q^(1_1:1i=f!<1[1W(a(}'}1D|O(=)g>|11121314151G/61v>y!')y{715'1:t_81>{91<:r!01W[(/51-{a1P|EO-+>V,*!y1/1 ]g[z1M;i)A1B1C1D1E1N(F1^)9;K{a(()G16,H1I16,J1M/K1a'j,5;L1U^M1N1n}O1q!P1>{N]Q13/R1S1H/B{S^0'8:T1U1i=V1W1X1Y1Z1`1 2.2",
+"-,n];=g[g};=+2F;v1%>I>@2#;H;d/b_>^#24{;,@'$2,>@'%2&2*2=2K>`-J:V=Y/G[-2I>#,M-h[;2q^.>`-/1,>x~I[h/`;.>.>S-`;W-y~>2R,@>d/,2K>]!W-P-N,P-S- ]N,-= ]S- ] ]-,W{~/!:'2)2j=H|!2%)~2c!c'T)11L{||{211g,S)~[_:e(7[5'L{v>:'z/]2z_tW!j<@>F-W-B~J<92n/U,`;g=.>o)j^+-t(S-O/S-c_T/]<=,S-_!Y,g/K-02@-+-Z/`;`;+>a2P-~/-=h) ]-,g=]-$,g=U!W=Y{E-R- >P!R,@,b2 > ]:!J>v~c2!:d2=%e2f2g2:/&(6,h2<^L^i2p{p,[]:/N(Z'{]51K^0^j2k2O(l2n(m2n2o2p2R(R^q2p_k,b>}[O]r2s2i=t2u2v2w2x2y2z2A2B2W#",
+"X=*>y^@^5{'>|!X!V=C2D2b{E2Q=a=F2)=G24!u)S'H2H2I2J2K2_,Q'&''}=-!;L2=-&;M2N2F'>>#;e_s)L-,>`;W-&,%>P-g=I>`;J>@-g=Q,J>X-W{<{W,0]N,b2$>U! ]v'T,v~J!J>L,v'b)~/R*J!O2P2Q2R2S2H^_~5'L^@_8'g, _3]T2h>W[F/U2k((/9;=)0^J{S~r>V2W2X2Y2Z2`2 3.3+3p{E/g,K(@3#3. $3%3&3*3=3-3;3>3,3",
+"W,;=F-%];,1>o)%]N-w^'3j) ~Y{b2 >Q,J!J!P})3r~!3+,R-P;Z;~3w'L,Q-Q-@,{3!3N,N,N,V-J>]3@-T-0]W-&,M,M,J>J>J>b2R-F=z'J!F=-=P-$,Y-^3W-N,/3&,M, ][*W-(3_3/3V-Y-H>M,h)V-'}. :3i=<3[3*)}3]~()J{K{&)11|3}'|/13P~23P^g'330(X2435363738393k=i=i=i=i=i=i=03a3. B-b3c3d3e3f3g3R = . ",
+"h3<=<=i3h3h3D~j3!/*_k3l3m$@.%.1+S O J@7.f%m3y$Y*7&7&n3y$@$P R@(@h#1+#.y n y f|o3p3q3E,V-G=W,`-H=V=`-@;p/J-=;I-@;L'*;H-'/;>I-P-P-G[P-W-P-O--=Y-A-/1Q!M>_3 ]q'D;J-@;2$r3s3t30!0's>b;v/u3v3w3x38/b~y3z3A3u,B3C3i=i=i=i=i=D3E3F3G3H3I3J3K3L3M3N3O3P3Q3R3S3T3U3V3>+~ . . ",
+">!W3X3}.V$H : i%Y3k%~ - * ,*Z3+ @ . @ . . . . $ . . . . . @ . I%@ + + v=+ * = { { { ~ J ) J A=d%}&) ) W%5$u@E.n `3r j' 42).4+4m' ,_3#,,>1>;;F-),I>1{@>;=#,W-G[@4#4b28+$4%4&4*4=4-4i=i=%_i_;4>4c^c^,4'4)4!4~4{4]4^4/4c.(4_4:4S$f&<4[4}4|414243444546474f3P/L%H @ @ . ",
+"@ @ @ . . . . . . . . . # . . . . @ @ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . K K $ . . w=I%& A%@ + @ { }&[$R@<+#.+@E.%.84t ]%m&R;+1Z;O'y,9404a4b4c4d4e4f4g4h4i4j4k4l4m4n4_ W%9%Y3o4p4q4r4s4t4u4v4w4x4y4z4G*A4B4C4D4E4F4G4V3L%5.@ . . . ",
+"@ @ . @ . @ . K . . . @ B%@ . . . . . . . . . a|K K . @ . . . . . . . . @ . . @ . . . . . . . . K . @ . . . K . . . x=. K . . . @ . . . . . . . . . . . . F%H4I4|&X#J4K4L4M4N4O4P4Q4R4S4T4U4V4W4X4Y4Z4`4 5.5+5@5#5$5%5&5*5=5-5;5>5,5'5)5!5~5{5]5^5/5(5L%H @ @ . . . ",
+"@ . @ . @ . . K K . . g&g&% . . . . . . . . . 2#. . @ . @ . @ @ @ . . . @ @ g&. L . . . . . . . K . K L @ . K . # . K . . . @ . . . . @ . . @ K @ . . % . K . % . . . _5:5. . . f&<5C=[5}5|515253545556575859505a5b5&3c5d5e5f5g5h5i5j5k5l5m5n5o5p5q5.@J @ @ . @ . . ",
+"# . # # . # . . . . . . . . . . . . . . . . @ . . K . K K # % . # . # . . . @ . g&K . . . . @ L . . . $ . . . . # . . . . . . . . . . . . . . . . . . . . K . $ @ . # . . :5r5s5_5t5r5u5v5&*w5x5y5z5A5B5C5D5E5F5G5X_H5I5J5K5L5M5N5O5P5Q5R5S5T5K[U5V5W5. . . . . . . ",
+"$ . % B%2#x=. . K x=x=x=x=K . # @ . # . # @ . $ x=. K g&. x=. @ . X5K . . @ . @ # K . @ @ . . K . % L % . . . @ . . @ . . . . . . . . . @ . . @ . @ . . . . . @ + @ . @ @ . B%g&2#L . . . . . . . . . 7 F 0+Y5Z5`5 6.6+6@6#6$6%6U3&6*6=6-6)!;6t.d@+&. . . . . . . . ",
+"% . @ . . . . . . . . . . . . . . . . . . . . . K . K . . . L g&x=K # % % . # . % . . K w=K . . . . # . f&# # . . . . . . @ . . . . @ @ . . . + . . @ @ . . . . @ . . @ @ . @ . . +&F%F%F%F%f&L t=>6. . . . . . $ ~@.&1$,6`+'6X3)6!6~6{6L%T I { & . % . . . @ . @ # ",
+". . # 2#K ]6,*~@,*^62 >*. . ~@$ ~@$ y%/6$ I%& # $ 2#(6# # # . . . . . . . . . . . . . . . . . . % _6% % % . g&@ + . . . . . . @ + . @ @ . . . . . @ . . . . . . . . . . . . . . . # # . L f&F%B%:6<6. # # . @ . . . . . # t5. s5r5[64&4 Y3n4,*@ . K }6K g&. @ @ . . ",
+"3&|61626263646566676869606a6b6c6d6e6f6g6h6i6j6k6l6m6l6n6o6p6q6r6O4s6t6N+M@L+O T@{#S#a$.%4&-*v=. . . . . . . . . . . . . . + @ @ @ @ . @ . @ + . . . @ . . . . @ . . . . . . . . @ . . . . . . . _6H%g&. @ I%H4a&u6v6w6x6H%B%g&I%t5r5u5_5y6z6z6r%4 '*A6,32#$ . . . . ",
+"B6C6D6E6F6G6H6I6H6J6K6L6M6N6O6P6Q6R6P6R6S6P6T6U6V6W6X6Y6Z6`6 7.7+7@7#7$7%7&7*7=7-7;7>7,7'7J6)7!7~7~7{7]7N3^79@5$R%$%M 2 . . . . . . . . @ . . . @ . @ @ @ . . . . . . . . . . . . . . . . . . . $ # # @ L # % x%/7n%(7_7:7<7[7}7' H #&. . u5[6|717273747471757+&B%67",
+"77879707a7b754b7b7c7d7e7f7C%g7h7i7j7k7l7m7n7o7p7q7p7r7p7s7t7u7v7t7w7x7y7z7'7A7B7C7D7E7F7G7H7I7J7J7K7L7M7M7N7O7P7P7Q7R7S7T7U7(4V7W7X7Y7Z7. . . . . . . . . . . @ . @ . . . . . . @ . . . . . . K K . . _6. . . . . . . $ }&`7 8.8+8@8[%#8#8 8. . . # /6t5}7$8%8v5&8*8",
+"J|i=i=i=i=i=i=i=i=i=i=i=i==8+=-8;8o<>8,8'8)8!8~8{8]8'8^8/8(8_8D7:8D7-2<8[8}8|8]|B618283848586878%[88a698S708a8b8c8d8e8f8g8h8i8j8k8l8d[m8n8o8p8q866r8~ 2 % . . . . g&. . . . % . :6. K :6g&s8. :6. @ # % . @ . @ + $ . . . . x=x=@ $ '*4&9&w6u#t8u8u%@ v8. . # t5w8_5",
+"$}x8y8z8A8B8C8D8E8F8G873H8i=I8~ J8K8L8M8N8O8P8Q8R8S8T8U8V8W8X8Y8Z8`8 9.9+9@9#9t7$9%9&9*9=9-9;9e8>9,9'9)9!9~9{9]988^9/9(9_9:9<9]|[9}9|91929:83949596979899909a9b9c9d9e9f9* * ~@y%. . . . . . @ @ + . . . . . K x=K . + . a|2#+ . . . . . . K '*^6-*g9`%C$Y7W#6$W5+ I%",
+"#}h9i9j9k9l9C>m9n9o9p9q9r9s9i=t9_5u9v9w9x9y9z9z9A9B9C9D9E9F9G9H9W8I9J9K949L9M9N9N9O9;7P9Q909R9d8S9T9U9T9T9V9V9V9{9S9W9K6X9Y9Z9`9 0S7.0+0K6V9@0|9#0$0.0d8%0&0*0'8=0-0;0>0,0'0)0!0]@7.X$~02 I%. . . . . . . . & @ @ . . @ . # ~@@ . . . . . . . K . . . . * 6&r%{0]0^0",
+"/0(0_0n9:0<0]/[0}0|010l9_[2030i=40f@506070809000a0b0c0B9d0e0f0g0h0K9i0j0k0l0m0n0o0N9p0%1q0r0s0t0d8U9u0&9>9v0d8b829w0K6x0X9.0y0z0A0B0C0.0e8D0E0&9F0u0Q9'958G0H0I0J0K0L0M0_8N0O0P0Q0X6R0W8S0T0U0V0$@.$s$. . . . . K 2#K g&g&. . . . . . . . . . x=. . . . . . . . $ ^6",
+"W0X0Y0X0Z0`0Z0 a.a+a@a#a$a%a&a*ai==a-a;a>a,a'a)a!aa0~a+7{a]a)8^a/a(aK9W8=7K7_a]8p7:a9ca|ada}8eafaga1aha}8,7iaA0B7jaka/a60lamanaoapaM8T8qarasata. . . . . K . @ . . . . . K . @ $ . . . % . K @ @ . . ",
+"uavawaxah9yazaAaBa4:k9CaDaEaFaGaHaIaJaKaLaMaNaOaPaQaRaSaP8Ta00'8UaF9VaWaH9Xa/aYaZaL9p7`a b;7%139D6.b+bR9@b#bP7$bP9%b38K61aR9'7&b*b&9A0=b-bia;b>b,b-:2a&9%b'b39@/,bp7o0}8)b!b~b{b1!&0jaK9]b&0X8^b/b(b_b:b7}8o0#bg8AbBb#b@bCbo0T9f89aDb6aEbzbu7FbGbGbO7HbEb<8IbE7=bJbKb9aLb#b9aMbh849Nb]bJbObWaObPbQbH9,8X8PbRbSbTbz9UbVbc0WbXbYbZb`b c.c+c@cC!~@. . . . . . . . . . x=:62#@ @ . g&",
+"M]#c$cr9%c&c*c=c-c;c>ccb,c'c:[C{)c!c~cx~{c]cO0^c/c(c_c:c7Kb1bkcNbN9jalcmcncObh0ocPb^bU8pc,8X8qcrcM3^csctcucvcwcxcyczcAcBcCc`bDcEcFcGcD 0&E%Hc. . . g&_6. @ @ . . ",
+"IcJcKcLcMcNcOcPcQcRcScTcUcVc:0WcXcYci=Zcf9`c d.d+d@dsc#d$dZ6%d&d*d=d:c-dS8;d,8>d,d'dG9X8n7]b]bAb)dH7H7!d~d0c{d]dt7^d|a:at7/dAbAbf8!b@b|a39(d_d#9.9:d8Wb3d 74d d5d6d7dP68d9d0d}badadbd`bcdddedfdgdV9hdY#s$. . . . . + . ",
+"r9idjdkdldAamdndodi=pdqdrdsdtdudvdwd!cxd1)S7ydzdAdgdBdCd8ucR05d5d>e,e'ena)e!e~e{e]e!e^e/e(e_e:ef,f'f)f!f~fY}{f]f^f/fFe7c(f_f:f",
+"/2g,g'gi=)g!gO8~g{gvc]gDdxc(b^gNa.7 f;0/g(g_g:gL8h,h'h)h 9!hj0~hEe{hmcK<]h^_sf^h bXa(8G9/hl7(h^a_h:hi,i'i)i!i~i{i]i^i/i(iAc_i:i0!8M8qcobXiYi}hR0@dZi`iYg jhi.jmg+jog@j#j$j%j&j*j=j-j;j>j,j'ji=)j!j~j{jTc]j^jvh/j(j_j",
+":j+a:0k,k'k)k!k~k{k",
+"]kQg^k/koh(k_kid:kdiHapd8df80*iMkuchiigNkOkPkQkRkSkTkUk9d*fVkWkXkYkZk`k l.l+l@li=#l$l2jg!.{%l^k)k&l*l=le~",
+"-l;lh9>l,l4:Lc'l)l!l~l{l]l^l/l(l_l:l8WbFhre|l1lP;2l3l4l5l6l7l8l9l0lalblcldlelflglhliljlGfklDdllmlnlolplqlrlsltlulvlwlxle}ylzl57(@AlBlClDlElQ7Fl%iGl90>8O8Hl+hIlHlXiJlKlLl*iMlqeNlkgOlPlQlRlSlTlUloaVlWlXlYlZl&f`l m.m+m@m#mi_$m%m*l:0m7i,m{j'mk9)mHaAa!mi=~m{mi=]m^m_aLl*i/mvcGf(m_m:miSjEmFmGmP6HmHlImJmZjKmLmMmNmOmQ6PmQmRmSmWkTmUmVmWmXmi=:[8iEa|j;kYm|kZmkdudm9`m n",
+"z_.n+n@n#n$nEaWc}03[%n&nAai=B|Ha*n=n-ni=r*;n>n,n'nFh`ese)n!n~n{nl7]n^n/n(nM2_nM,)-:no,o@k'o)o!o~o{o]o^o/obc(o_o:op,p'p)p!p~p{p]p^p/p(pM8Y6_p:pq,q'q)qd Gd!qOk~q{qfmv9]q^q/q(q_q:qr`b,r'r)r!r~r{r]r^r/r(r_r:rjKrLrMrNrOrPrQrRrSrTrUrVrWrXrYrZr`r s.sSn+s@s#s$s%s&sv9^n*s=s@$-s;sh:>s,s'sil)s!s~s{s]s^s/s(s_s:st,t't)t!t~t{txr/i(i]t^t/t(ti=k=6sfq:l%__tP)k=:tq;i=HaB]&gXc",
+"R:7TtUtVtWtXtYtZt`tCq~t u6h.u&j+up6@ui=Qcypk=gq#uu*li=,uj&'u)uzf!u~uzew9{u+o]uZg2hMoT6^uKmZjVmWl/ujooaa35h(u_u:uwnv,v'v)v!p!v~v{v]vBt^v/v(v_v:vw,w'w(e)w!pi=K(!wO)vo@3%_k=i=;4~wk=%_8bK(~j",
+"{w]w^w/w(w;;^+4.'#g@9 _w:wykxLa9hVl[bbh,x'x>w)xws@v0qa3#kdv!x~xi=JuK({x]xJu_t#ui=M{^xtpi=fqdb",
+"/x(x_x:xycg,yZ6Fh'yLh=y@wvx2z*jLa,z)x)xP6]e'z u1sAOq,A'A)A!Atm'A~Avz{A]AYy)AqxFt1}^Aen/A(A@z_A6u:Agzzsw-y*B=B-B;B>Bds,BM{xa'BcyGz)B8spp!B*l",
+". . . . . . . . . . . + * = { 5.9 S {+L%A ~B{B@!(3>_]B^BMz/B8ply(B22_B:B=8+OB&.d.;6PBV3QBRBiySBTB_zUB 1VBWBXBYBZB`B C.C+C@C#C$C%CC,C'C)C!C~C{C]C^C/C(C_C:C+L%pCd.^zqC(3L;KzrCsCtCuCvCwCxC(ByCXBzCACM5m5[z}zq*BCi5CCDCECwAFCGCHCICJCKC54LCMC3zNCSrOCPCQCRCSCTCpu`u}wVlzz;j>A{AUCVCWCOAXCXyYCZC`CXy DYvnmFrwg5yYkDx.DN6rvkt+Dit>z@D#D$D%D&Di=*DKB=D-D;DIz&nqp>D",
+". . . . . . . . . . . . . . . . . . . . . . . . @ . . . . . . @ * = ~ H 9 Q P {+F k j d.^zqCkAV3N2rC,D'D)D!DS5qA~DyC{D]DRzFAEIl;E_p>iRa/r!^S+,E'E2q)E!5!E~E{E",
+". . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . @ . . . . . . . . . . . . . . . . . . @ + * { ~ ) H a+a+R F ^.&.e.iA^z@!kAP/N2]E^E/E(E/B_EOzPz:EF+A.^._.3F;6^zjA4FP/QB5F70,D6FTB_E%67F 1N_nyXBn5;F}E+gF 1hFiFjF^.J . . @ @ @ @ . . . . . . . . . . . . . . . . . . . . . . . @ = = = J 5.9 9 Q R >+^.&.e.e.oDjA@!ZEP/N2RBLz*6'DjyTBeF_E%6S5pD",
+". . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ~ lykF<'lFmFnFoFpFqFQB@ @ . @ . . . . . . . . . . . . . . . @ . . . @ . @ . @ . . . . . . . . @ @ + { { ~ J H I a+a+>+1+F j 3Fd.oD^zrFjAg3kAP/",
+". . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . J ,DsFtFuFvFwFxF
+#include
+
+// Qmitk
+#include "QmitkIGTTrackingDataEvaluationView.h"
+#include "QmitkStdMultiWidget.h"
+
+// Qt
+#include
+#include
+
+// MITK
+#include "mitkNavigationDataCSVSequentialPlayer.h"
+#include
+#include
+#include
+
+//ITK
+#include
+
+//VNL
+#include
+
+//vtk headers
+#include
+#include
+#include
+
+
+const std::string QmitkIGTTrackingDataEvaluationView::VIEW_ID = "org.mitk.views.igttrackingdataevaluation";
+
+QmitkIGTTrackingDataEvaluationView::QmitkIGTTrackingDataEvaluationView()
+: QmitkFunctionality()
+, m_Controls( 0 )
+, m_MultiWidget( NULL )
+, m_scalingfactor(1)
+{
+m_CSVtoXMLInputFilenameVector = std::vector();
+m_CSVtoXMLOutputFilenameVector = std::vector();
+}
+
+QmitkIGTTrackingDataEvaluationView::~QmitkIGTTrackingDataEvaluationView()
+{
+}
+
+
+void QmitkIGTTrackingDataEvaluationView::CreateQtPartControl( QWidget *parent )
+{
+ // build up qt view, unless already done
+ if ( !m_Controls )
+ {
+ // create GUI widgets from the Qt Designer's .ui file
+ m_Controls = new Ui::QmitkIGTTrackingDataEvaluationViewControls;
+ m_Controls->setupUi( parent );
+
+ connect( m_Controls->m_LoadInputFileList, SIGNAL(clicked()), this, SLOT(OnLoadFileList()) );
+ connect( m_Controls->m_StartEvaluation, SIGNAL(clicked()), this, SLOT(OnEvaluateData()) );
+ connect( m_Controls->m_AddToCurrentList, SIGNAL(clicked()), this, SLOT(OnAddToCurrentList()) );
+ connect( m_Controls->m_GeneratePointSetOfMeanPositions, SIGNAL(clicked()), this, SLOT(OnGeneratePointSet()) );
+ connect( m_Controls->m_GenerateRotationLines, SIGNAL(clicked()), this, SLOT(OnGenerateRotationLines()) );
+ connect( m_Controls->m_GeneratePointSet, SIGNAL(clicked()), this, SLOT(OnGenerateGroundTruthPointSet()) );
+ connect( m_Controls->m_Convert, SIGNAL(clicked()), this, SLOT(OnConvertCSVtoXMLFile()) );
+ connect( m_Controls->m_loadCSVtoXMLInputList, SIGNAL(clicked()), this, SLOT(OnCSVtoXMLLoadInputList()) );
+ connect( m_Controls->m_loadCSVtoXMLOutputList, SIGNAL(clicked()), this, SLOT(OnCSVtoXMLLoadOutputList()) );
+ connect( m_Controls->m_OrientationCalculationGenerateReference, SIGNAL(clicked()), this, SLOT(OnOrientationCalculation_CalcRef()) );
+ connect( m_Controls->m_OrientationCalculationWriteOrientationsToFile, SIGNAL(clicked()), this, SLOT(OnOrientationCalculation_CalcOrientandWriteToFile()) );
+ connect( m_Controls->m_GeneratePointSetsOfSinglePositions, SIGNAL(clicked()), this, SLOT(OnGeneratePointSetsOfSinglePositions()) );
+ }
+}
+
+void QmitkIGTTrackingDataEvaluationView::OnOrientationCalculation_CalcRef()
+{
+if(m_FilenameVector.size() != 3)
+ {
+ MessageBox("Need exactly three points as reference, aborting!");
+ return;
+ }
+
+//start loop and iterate through all files of list
+for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+ //check if the stream is valid and skip file if not
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+ continue;
+ }
+ */
+
+ //create evaluation filter
+ mitk::NavigationDataEvaluationFilter::Pointer myEvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+
+ //connect pipeline
+ for (int j=0; jGetNumberOfOutputs(); j++) {myEvaluationFilter->SetInput(j,myPlayer->GetOutput(j));}
+
+
+ //update pipline until number of samlples is reached
+ for (int j=0; jm_NumberOfSamples->value(); j++)
+ {myEvaluationFilter->Update();}
+
+ //store mean position as reference
+ switch (i)
+ {
+ case 0:
+ m_RefPoint1 = myEvaluationFilter->GetPositionMean(0);
+ break;
+ case 1:
+ m_RefPoint2 = myEvaluationFilter->GetPositionMean(0);
+ break;
+ case 2:
+ m_RefPoint3 = myEvaluationFilter->GetPositionMean(0);
+ break;
+ }
+ }
+ MessageBox("Created Reference!");
+}
+
+void QmitkIGTTrackingDataEvaluationView::OnOrientationCalculation_CalcOrientandWriteToFile()
+{
+//start loop and iterate through all files of list
+for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+ //check if the stream is valid and skip file if not
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+ continue;
+ }
+ */
+
+
+
+ //open file header
+ QString outputname = QString(m_FilenameVector.at(i).c_str()) + "_orientationFile.csv";
+ m_CurrentWriteFile.open(outputname.toStdString().c_str(), std::ios::out);
+ if (m_CurrentWriteFile.bad())
+ {
+ MessageBox("Error: Can't open output file!");
+ return;
+ }
+
+ //write header to file
+ m_CurrentWriteFile << "Nr;Calypso_Time;Valid_Reference;MeasureTool_Measurement-Tool[x];MeasureTool_Measurement-Tool[y];MeasureTool_Measurement-Tool[z];MeasureTool_Measurement-Tool[qx];MeasureTool_Measurement-Tool[qy];MeasureTool_Measurement-Tool[qz];MeasureTool_Measurement-Tool[qr]\n";
+
+ //update pipeline until number of samples is reached
+ int step = 0;
+ mitk::Point3D point1,point2,point3;
+ mitk::Quaternion current_orientation;
+
+ for (int j=0; !myPlayer->IsAtEnd(); j++)
+ {
+ myPlayer->Update();
+ mitk::NavigationData::Pointer currentNavData = myPlayer->GetOutput(0);
+ switch (step)
+ {
+ case 0:
+ step++;
+ point1 = currentNavData->GetPosition();
+ break;
+ case 1:
+ step++;
+ point2 = currentNavData->GetPosition();
+ break;
+ case 2:
+ step=0;
+ point3 = currentNavData->GetPosition();
+
+ //compute transform from reference to current points
+ if (point1[0] == 0 &&
+ point1[1] == 0 &&
+ point1[2] == 0 &&
+ point2[0] == 0 &&
+ point2[1] == 0 &&
+ point2[2] == 0 &&
+ point3[0] == 0 &&
+ point3[1] == 0 &&
+ point3[2] == 0
+ ) current_orientation.fill(0);
+ else
+ {
+ /* Drehen um eine Achse um das "Umschlagen" zu vermeiden
+ itk::Matrix rot180degreeAroundY;
+ rot180degreeAroundY.Fill(0);
+ rot180degreeAroundY[0][0] = -1;
+ rot180degreeAroundY[1][1] = 1;
+ rot180degreeAroundY[2][2] = -1;
+ point1 = rot180degreeAroundY * point1;
+ point2 = rot180degreeAroundY * point2;
+ point3 = rot180degreeAroundY * point3;
+ */
+
+ vtkSmartPointer transform = vtkSmartPointer::New();
+ vtkSmartPointer sourcePoints = vtkSmartPointer::New();
+ double sourcepoint1[3] = {point1[0],point1[1],point1[2]};
+ double sourcepoint2[3] = {point2[0],point2[1],point2[2]};
+ double sourcepoint3[3] = {point3[0],point3[1],point3[2]};
+ sourcePoints->InsertNextPoint(sourcepoint1);
+ sourcePoints->InsertNextPoint(sourcepoint2);
+ sourcePoints->InsertNextPoint(sourcepoint3);
+ vtkSmartPointer targetPoints = vtkSmartPointer::New();
+ double targetpoint1[3] = {m_RefPoint1[0],m_RefPoint1[1],m_RefPoint1[2]};
+ double targetpoint2[3] = {m_RefPoint2[0],m_RefPoint2[1],m_RefPoint2[2]};
+ double targetpoint3[3] = {m_RefPoint3[0],m_RefPoint3[1],m_RefPoint3[2]};
+ targetPoints->InsertNextPoint(targetpoint1);
+ targetPoints->InsertNextPoint(targetpoint2);
+ targetPoints->InsertNextPoint(targetpoint3);
+
+ transform->SetSourceLandmarks(sourcePoints);
+ transform->SetTargetLandmarks(targetPoints);
+ transform->Modified();
+ transform->Update();
+
+ mitk::Transform::Pointer newTransform = mitk::Transform::New();
+ newTransform->SetMatrix(transform->GetMatrix());
+ current_orientation = newTransform->GetOrientation();
+
+ //add pointset with the three positions
+ if((j>15) && (j<18))
+ {
+ mitk::DataNode::Pointer newNode = mitk::DataNode::New();
+ mitk::PointSet::Pointer newPointSet = mitk::PointSet::New();
+ newPointSet->InsertPoint(0,point1);
+ newPointSet->InsertPoint(1,point2);
+ newPointSet->InsertPoint(2,point3);
+ QString name = QString(m_FilenameVector.at(i).c_str());
+ newNode->SetName(name.toStdString().c_str());
+ newNode->SetData(newPointSet);
+ newNode->SetFloatProperty("pointsize",0.1);
+ this->GetDataStorage()->Add(newNode);
+ }
+ }
+
+ break;
+ }
+ m_CurrentWriteFile << i << ";";
+ m_CurrentWriteFile << currentNavData->GetTimeStamp() << ";"; //IMPORTANT: change to GetIGTTimeStamp in new version!
+ m_CurrentWriteFile << "true;";
+ m_CurrentWriteFile << currentNavData->GetPosition()[0] << ";";
+ m_CurrentWriteFile << currentNavData->GetPosition()[1] << ";";
+ m_CurrentWriteFile << currentNavData->GetPosition()[2] << ";";
+ m_CurrentWriteFile << current_orientation.x() << ";";
+ m_CurrentWriteFile << current_orientation.y() << ";";
+ m_CurrentWriteFile << current_orientation.z() << ";";
+ m_CurrentWriteFile << current_orientation.r() << ";";
+ m_CurrentWriteFile << "\n";
+ }
+ //close output file
+ m_CurrentWriteFile.close();
+
+ }
+ MessageBox("Finished!");
+}
+
+void QmitkIGTTrackingDataEvaluationView::StdMultiWidgetAvailable (QmitkStdMultiWidget &stdMultiWidget)
+{
+ m_MultiWidget = &stdMultiWidget;
+}
+
+
+void QmitkIGTTrackingDataEvaluationView::StdMultiWidgetNotAvailable()
+{
+ m_MultiWidget = NULL;
+}
+
+void QmitkIGTTrackingDataEvaluationView::OnAddToCurrentList()
+{
+ //read in filename
+ QString filename = QFileDialog::getOpenFileName(NULL,tr("Open Measurement Filename List"), "/", tr("All Files (*.*)"));
+ if (filename.isNull()) return;
+
+ /*
+ //save old locale
+ char * oldLocale;
+ oldLocale = setlocale( LC_ALL, 0 );
+
+ //define own locale
+ std::locale C("C");
+ setlocale( LC_ALL, "C" );
+ */ //TODO: check if this is needed here, and load old locale if yes
+
+ //read file
+ std::ifstream file;
+ file.open(filename.toStdString().c_str(), std::ios::in);
+ if (file.good())
+ {
+ //read out file
+ file.seekg(0L, std::ios::beg); // move to begin of file
+ while (! file.eof())
+ {
+ std::string buffer;
+ std::getline(file,buffer); // read out file line by line
+ if (buffer.size() > 0)
+ {
+ std::string thisFilename = "";
+ if (m_Controls->m_AddPath->isChecked()) thisFilename = m_Controls->m_ListPath->text().toStdString();
+ thisFilename.append(buffer);
+ m_FilenameVector.push_back(thisFilename);
+ }
+ }
+ }
+
+ //fill list at GUI
+ m_Controls->m_FileList->clear();
+ for(unsigned int i=0; im_FileList);}
+
+}
+
+void QmitkIGTTrackingDataEvaluationView::OnLoadFileList()
+ {
+ m_FilenameVector = std::vector();
+ m_FilenameVector.clear();
+ OnAddToCurrentList();
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnEvaluateData()
+{
+
+//open output file
+m_CurrentWriteFile.open(std::string(m_Controls->m_OutputFilename->text().toUtf8()).c_str(), std::ios::out);
+if (m_CurrentWriteFile.bad())
+ {
+ MessageBox("Error: Can't open output file!");
+ return;
+ }
+
+//write output file header
+WriteHeader();
+
+//start loop and iterate through all files of list
+for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+
+ //check if the stream is valid and skip file if not
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+
+ continue;
+ }
+ */
+
+ //create evaluation filter
+ mitk::NavigationDataEvaluationFilter::Pointer myEvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+
+ //connect pipeline
+ for (int j=0; jGetNumberOfOutputs(); j++) {myEvaluationFilter->SetInput(j,myPlayer->GetOutput(j));}
+
+
+
+ //update pipline until number of samlples is reached
+ for (int j=0; jm_NumberOfSamples->value(); j++)
+ {
+ myEvaluationFilter->Update();
+ //Debug output:
+ //std::cout.precision(5);
+ //std::cout << "Euler " << j << ";" << myPlayer->GetOutput()->GetOrientation().rotation_euler_angles()[0] << ";" << myPlayer->GetOutput()->GetOrientation().rotation_euler_angles()[1] << ";" << myPlayer->GetOutput()->GetOrientation().rotation_euler_angles()[2] << "\n";
+ }
+
+ //write result to output file
+ WriteDataSet(myEvaluationFilter,m_FilenameVector.at(i));
+
+ }
+
+//close output file
+m_CurrentWriteFile.close();
+
+//calculate angles if option is on
+if(m_Controls->m_settingDifferenceAngles->isChecked() || m_Controls->m_DifferencesSLERP->isChecked()) CalculateDifferenceAngles();
+
+MessageBox("Finished!");
+
+}
+
+void QmitkIGTTrackingDataEvaluationView::OnGeneratePointSetsOfSinglePositions()
+{
+ m_scalingfactor = m_Controls->m_ScalingFactor->value();
+
+ //start loop and iterate through all files of list
+ for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+
+ //check if the stream is valid and skip file if not
+/*
+if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+
+ continue;
+ }
+*/
+
+ //update pipline until number of samlples is reached and store every single point
+ for (int j=0; jm_NumberOfSamples->value(); j++)
+ {
+ myPlayer->Update();
+ mitk::Point3D thisPoint = myPlayer->GetOutput()->GetPosition();
+ thisPoint[0] *= m_scalingfactor;
+ thisPoint[1] *= m_scalingfactor;
+ thisPoint[2] *= m_scalingfactor;
+ thisPointSet->InsertPoint(j,thisPoint);
+ }
+
+ //add point set to data storage
+ mitk::DataNode::Pointer newNode = mitk::DataNode::New();
+ QString name = this->m_Controls->m_prefix->text() + QString("PointSet_of_All_Positions_") + QString::number(i);
+ newNode->SetName(name.toStdString());
+ newNode->SetData(thisPointSet);
+ this->GetDataStorage()->Add(newNode);
+ }
+
+
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnGeneratePointSet()
+ {
+ m_scalingfactor = m_Controls->m_ScalingFactor->value();
+
+ mitk::PointSet::Pointer generatedPointSet = mitk::PointSet::New();
+
+ //start loop and iterate through all files of list
+ for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+
+ //check if the stream is valid and skip file if not
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+
+ continue;
+ }
+ */
+
+ //create evaluation filter
+ mitk::NavigationDataEvaluationFilter::Pointer myEvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+
+ //connect pipeline
+ for (int j=0; jGetNumberOfOutputs(); j++) {myEvaluationFilter->SetInput(j,myPlayer->GetOutput(j));}
+
+ //update pipline until number of samlples is reached
+ for (int j=0; jm_NumberOfSamples->value(); j++) {myEvaluationFilter->Update();}
+
+ //add mean position to point set
+ mitk::Point3D meanPos = myEvaluationFilter->GetPositionMean(0);
+ if (m_scalingfactor!=1)
+ {
+ meanPos[0] *= m_scalingfactor;
+ meanPos[1] *= m_scalingfactor;
+ meanPos[2] *= m_scalingfactor;
+ }
+ generatedPointSet->InsertPoint(i,meanPos);
+ }
+
+ //add point set to data storage
+ mitk::DataNode::Pointer newNode = mitk::DataNode::New();
+ QString name = this->m_Controls->m_prefix->text() + "PointSet_of_Mean_Positions";
+ newNode->SetName(name.toStdString());
+ newNode->SetData(generatedPointSet);
+ this->GetDataStorage()->Add(newNode);
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnGenerateRotationLines()
+ {
+ m_scalingfactor = m_Controls->m_ScalingFactor->value();
+
+ //start loop and iterate through all files of list
+ for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+
+ //check if the stream is valid and skip file if not
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+ }
+ else
+ */
+ {
+ //create evaluation filter
+ mitk::NavigationDataEvaluationFilter::Pointer myEvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+
+ //connect pipeline
+ for (int j=0; jGetNumberOfOutputs(); j++) {myEvaluationFilter->SetInput(j,myPlayer->GetOutput(j));}
+
+ //update pipline until number of samlples is reached
+ for (int j=0; jm_NumberOfSamples->value(); j++)
+ {
+ myEvaluationFilter->Update();
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+ continue;
+ }
+ */
+ }
+
+ if (!myPlayer->IsAtEnd()) continue;
+
+ //create line from mean pos to a second point which lies along the sensor (1,0,0 in tool coordinates for aurora)
+ mitk::Point3D meanPos = myEvaluationFilter->GetPositionMean(0);
+ if(m_scalingfactor!=1)
+ {
+ meanPos[0] *= m_scalingfactor;
+ meanPos[1] *= m_scalingfactor;
+ meanPos[2] *= m_scalingfactor;
+ }
+ mitk::Point3D secondPoint;
+ mitk::Point3D thirdPoint;
+ mitk::Point3D fourthPoint;
+
+ mitk::FillVector3D(secondPoint,2,0,0); //X
+ vnl_vector secondPointTransformed = myEvaluationFilter->GetQuaternionMean(0).rotation_matrix_transpose().transpose() * secondPoint.Get_vnl_vector() + meanPos.Get_vnl_vector();
+ mitk::Point3D secondPointTransformedMITK;
+ mitk::FillVector3D(secondPointTransformedMITK,secondPointTransformed[0],secondPointTransformed[1],secondPointTransformed[2]);
+
+ mitk::FillVector3D(thirdPoint,0,4,0); //Y
+ vnl_vector thirdPointTransformed = myEvaluationFilter->GetQuaternionMean(0).rotation_matrix_transpose().transpose() * thirdPoint.Get_vnl_vector() + meanPos.Get_vnl_vector();
+ mitk::Point3D thirdPointTransformedMITK;
+ mitk::FillVector3D(thirdPointTransformedMITK,thirdPointTransformed[0],thirdPointTransformed[1],thirdPointTransformed[2]);
+
+ mitk::FillVector3D(fourthPoint,0,0,6); //Z
+ vnl_vector fourthPointTransformed = myEvaluationFilter->GetQuaternionMean(0).rotation_matrix_transpose().transpose() * fourthPoint.Get_vnl_vector() + meanPos.Get_vnl_vector();
+ mitk::Point3D fourthPointTransformedMITK;
+ mitk::FillVector3D(fourthPointTransformedMITK,fourthPointTransformed[0],fourthPointTransformed[1],fourthPointTransformed[2]);
+
+
+ mitk::PointSet::Pointer rotationLine = mitk::PointSet::New();
+ rotationLine->InsertPoint(0,secondPointTransformedMITK);
+ rotationLine->InsertPoint(1,meanPos);
+ rotationLine->InsertPoint(2,thirdPointTransformedMITK);
+ rotationLine->InsertPoint(3,meanPos);
+ rotationLine->InsertPoint(4,fourthPointTransformedMITK);
+
+ mitk::DataNode::Pointer newNode = mitk::DataNode::New();
+ QString nodeName = this->m_Controls->m_prefix->text() + "RotationLineNumber" + QString::number(i);
+ newNode->SetName(nodeName.toStdString());
+ newNode->SetData(rotationLine);
+ newNode->SetBoolProperty("show contour",true);
+ newNode->SetFloatProperty("pointsize",0.5);
+ this->GetDataStorage()->Add(newNode);
+ }
+ }
+
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnGenerateGroundTruthPointSet()
+ {
+ mitk::PointSet::Pointer generatedPointSet = mitk::PointSet::New();
+ int currentPointID = 0;
+ mitk::Point3D currentPoint;
+ mitk::FillVector3D(currentPoint,0,0,0);
+ for (int i=0; im_PointNumber1->value(); i++)
+ {
+ for (int j=0; jm_PointNumber2->value(); j++)
+ {
+ generatedPointSet->InsertPoint(currentPointID,currentPoint);
+ currentPointID++;
+ currentPoint[1] += m_Controls->m_PointDistance->value();
+ }
+ currentPoint[1] = 0;
+ currentPoint[2] += m_Controls->m_PointDistance->value();
+ }
+ mitk::DataNode::Pointer newNode = mitk::DataNode::New();
+ QString nodeName = "GroundTruthPointSet_" + QString::number(m_Controls->m_PointNumber1->value()) + "x" + QString::number(m_Controls->m_PointNumber2->value()) + "_(" + QString::number(m_Controls->m_PointDistance->value()) + "mm)";
+ newNode->SetName(nodeName.toStdString());
+ newNode->SetData(generatedPointSet);
+ this->GetDataStorage()->Add(newNode);
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnConvertCSVtoXMLFile()
+ {
+ if(m_Controls->m_ConvertSingleFile->isChecked())
+ { //convert one file
+
+ int lines = ConvertOneFile(this->m_Controls->m_InputCSV->text().toStdString(),this->m_Controls->m_OutputXML->text().toStdString());
+
+ QString result = "Converted one file with" + QString::number(lines) + " data sets";
+ MessageBox(result.toStdString());
+ }
+ else //converte file list
+ {
+ if(m_CSVtoXMLInputFilenameVector.empty() || m_CSVtoXMLOutputFilenameVector.empty())
+ {
+ MessageBox("Error: one list is not loaded!");
+ return;
+ }
+ else if(m_CSVtoXMLInputFilenameVector.size() != m_CSVtoXMLOutputFilenameVector.size())
+ {
+ MessageBox("Error: lists do not have the same number of files!");
+ return;
+ }
+ for(int i=0; i < m_CSVtoXMLInputFilenameVector.size(); i++)
+ {int lines = ConvertOneFile(m_CSVtoXMLInputFilenameVector.at(i),m_CSVtoXMLOutputFilenameVector.at(i));}
+ QString result = "Converted " + QString::number(m_CSVtoXMLInputFilenameVector.size()) + " files from file list!";
+ MessageBox(result.toStdString());
+ }
+ }
+
+int QmitkIGTTrackingDataEvaluationView::ConvertOneFile(std::string inputFilename, std::string outputFilename)
+ {
+ std::vector myNavigationDatas = GetNavigationDatasFromFile(inputFilename);
+ mitk::NavigationDataRecorderDeprecated::Pointer myRecorder = mitk::NavigationDataRecorderDeprecated::New();
+ myRecorder->SetFileName(outputFilename.c_str());
+ mitk::NavigationData::Pointer input = mitk::NavigationData::New();
+ if (m_Controls->m_ConvertCSV->isChecked()) myRecorder->SetOutputFormat(mitk::NavigationDataRecorderDeprecated::csv);
+ myRecorder->AddNavigationData(input);
+ myRecorder->StartRecording();
+ for (int i=0; iGraft(myNavigationDatas.at(i));
+ myRecorder->Update();
+ }
+ myRecorder->StopRecording();
+ return myNavigationDatas.size();
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnCSVtoXMLLoadInputList()
+ {
+ //read in filename
+ QString filename = QFileDialog::getOpenFileName(NULL,tr("Open Measurement Filename List"), "/", tr("All Files (*.*)"));
+ if (filename.isNull()) return;
+
+ m_CSVtoXMLInputFilenameVector = this->GetFileContentLineByLine(filename.toStdString());
+
+ m_Controls->m_labelCSVtoXMLInputList->setText("READY");
+ }
+
+void QmitkIGTTrackingDataEvaluationView::OnCSVtoXMLLoadOutputList()
+ {
+ //read in filename
+ QString filename = QFileDialog::getOpenFileName(NULL,tr("Open Measurement Filename List"), "/", tr("All Files (*.*)"));
+ if (filename.isNull()) return;
+
+ m_CSVtoXMLOutputFilenameVector = this->GetFileContentLineByLine(filename.toStdString());
+
+ m_Controls->m_labelCSVtoXMLOutputList->setText("READY");
+ }
+
+void QmitkIGTTrackingDataEvaluationView::MessageBox(std::string s)
+ {
+ QMessageBox msgBox;
+ msgBox.setText(s.c_str());
+ msgBox.exec();
+ }
+
+void QmitkIGTTrackingDataEvaluationView::WriteHeader()
+ {
+
+ m_CurrentWriteFile << "Filename;";
+ m_CurrentWriteFile << "N;";
+ m_CurrentWriteFile << "N_invalid;";
+ m_CurrentWriteFile << "Percentage_invalid;";
+
+ if(m_Controls->m_settingPosMean->isChecked())
+ {
+ m_CurrentWriteFile << "Position_Mean[x];";
+ m_CurrentWriteFile << "Position_Mean[y];";
+ m_CurrentWriteFile << "Position_Mean[z];";
+ }
+
+ if(m_Controls->m_settingPosStabw->isChecked())
+ {
+ m_CurrentWriteFile << "Position_StandDev[x];";
+ m_CurrentWriteFile << "Position_StandDev[y];";
+ m_CurrentWriteFile << "Position_StandDev[z];";
+ }
+
+ if(m_Controls->m_settingPosSampleStabw->isChecked())
+ {
+ m_CurrentWriteFile << "Position_SampleStandDev[x];";
+ m_CurrentWriteFile << "Position_SampleStandDev[y];";
+ m_CurrentWriteFile << "Position_SampleStandDev[z];";
+ }
+
+ if(m_Controls->m_settingQuaternionMean->isChecked())
+ {
+ m_CurrentWriteFile << "Quaternion_Mean[qx];";
+ m_CurrentWriteFile << "Quaternion_Mean[qy];";
+ m_CurrentWriteFile << "Quaternion_Mean[qz];";
+ m_CurrentWriteFile << "Quaternion_Mean[qr];";
+ }
+
+ if(m_Controls->m_settionQuaternionStabw->isChecked())
+ {
+ m_CurrentWriteFile << "Quaternion_StandDev[qx];";
+ m_CurrentWriteFile << "Quaternion_StandDev[qy];";
+ m_CurrentWriteFile << "Quaternion_StandDev[qz];";
+ m_CurrentWriteFile << "Quaternion_StandDev[qr];";
+ }
+
+ if(m_Controls->m_settingPosErrorMean->isChecked()) m_CurrentWriteFile << "PositionError_Mean;";
+
+ if(m_Controls->m_settingPosErrorStabw->isChecked()) m_CurrentWriteFile << "PositionError_StandDev;";
+
+ if(m_Controls->m_settingPosErrorSampleStabw->isChecked()) m_CurrentWriteFile << "PositionError_SampleStandDev;";
+
+ if(m_Controls->m_settingPosErrorRMS->isChecked()) m_CurrentWriteFile << "PositionError_RMS;";
+
+ if(m_Controls->m_settingPosErrorMedian->isChecked()) m_CurrentWriteFile << "PositionError_Median;";
+
+ if(m_Controls->m_settingPosErrorMinMax->isChecked())
+ {
+ m_CurrentWriteFile << "PositionError_Max;";
+ m_CurrentWriteFile << "PositionError_Min;";
+ }
+
+ if(m_Controls->m_settingEulerMean->isChecked())
+ {
+ m_CurrentWriteFile << "Euler_tx;";
+ m_CurrentWriteFile << "Euler_ty;";
+ m_CurrentWriteFile << "Euler_tz;";
+ }
+
+ if(m_Controls->m_settingEulerRMS->isChecked())
+ {
+ m_CurrentWriteFile << "EulerErrorRMS (rad);";
+ m_CurrentWriteFile << "EulerErrorRMS (grad);";
+ }
+
+ m_CurrentWriteFile << "\n";
+
+ }
+
+void QmitkIGTTrackingDataEvaluationView::WriteDataSet(mitk::NavigationDataEvaluationFilter::Pointer myEvaluationFilter, std::string dataSetName)
+ {
+ if (myEvaluationFilter->GetNumberOfOutputs()==0) m_CurrentWriteFile << "Error: no input \n";
+ else
+ {
+ m_CurrentWriteFile << dataSetName << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetNumberOfAnalysedNavigationData(0) << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetNumberOfInvalidSamples(0) << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPercentageOfInvalidSamples(0) << ";";
+
+
+ if(m_Controls->m_settingPosMean->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionMean(0)[0] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionMean(0)[1] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionMean(0)[2] << ";";
+ }
+
+ if(m_Controls->m_settingPosStabw->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionStandardDeviation(0)[0] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionStandardDeviation(0)[1] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionStandardDeviation(0)[2] << ";";
+ }
+
+ if(m_Controls->m_settingPosSampleStabw->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionSampleStandardDeviation(0)[0] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionSampleStandardDeviation(0)[1] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionSampleStandardDeviation(0)[2] << ";";
+ }
+
+ if(m_Controls->m_settingQuaternionMean->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionMean(0).x() << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionMean(0).y() << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionMean(0).z() << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionMean(0).r() << ";";
+ }
+
+ if(m_Controls->m_settionQuaternionStabw->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionStandardDeviation(0).x() << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionStandardDeviation(0).y() << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionStandardDeviation(0).z() << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetQuaternionStandardDeviation(0).r() << ";";
+ }
+
+ if(m_Controls->m_settingPosErrorMean->isChecked()) m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorMean(0) << ";";
+ if(m_Controls->m_settingPosErrorStabw->isChecked()) m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorStandardDeviation(0) << ";";
+ if(m_Controls->m_settingPosErrorSampleStabw->isChecked()) m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorSampleStandardDeviation(0) << ";";
+ if(m_Controls->m_settingPosErrorRMS->isChecked()) m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorRMS(0) << ";";
+ if(m_Controls->m_settingPosErrorMedian->isChecked()) m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorMedian(0) << ";";
+ if(m_Controls->m_settingPosErrorMinMax->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorMax(0) << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetPositionErrorMin(0) << ";";
+ }
+
+ if(m_Controls->m_settingEulerMean->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetEulerAnglesMean(0)[0] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetEulerAnglesMean(0)[1] << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetEulerAnglesMean(0)[2] << ";";
+ }
+
+ if(m_Controls->m_settingEulerRMS->isChecked())
+ {
+ m_CurrentWriteFile << myEvaluationFilter->GetEulerAnglesRMS(0) << ";";
+ m_CurrentWriteFile << myEvaluationFilter->GetEulerAnglesRMSDegree(0) << ";";
+ }
+
+ m_CurrentWriteFile << "\n";
+
+ }
+
+
+ }
+
+void QmitkIGTTrackingDataEvaluationView::CalculateDifferenceAngles()
+ {
+
+ std::vector EvaluationDataCollection;
+
+ //start loop and iterate through all files of list: store the evaluation data
+ for(int i=0; iSetFiletype(mitk::NavigationDataCSVSequentialPlayer::ManualLoggingCSV);
+ myPlayer->SetFileName(m_FilenameVector.at(i));
+
+
+
+ //create evaluation filter
+ mitk::NavigationDataEvaluationFilter::Pointer myEvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+
+ //check if the stream is valid and skip file if not
+ /*
+ if (!myPlayer->GetStreamValid())
+ {
+ MITK_ERROR << "Error in file " << m_FilenameVector.at(i) << ": " << myPlayer->GetErrorMessage() << " ; Skipping file!";
+ }
+ else
+ */
+ {
+ //connect pipeline
+ for (int j=0; jGetNumberOfOutputs(); j++) {myEvaluationFilter->SetInput(j,myPlayer->GetOutput(j));}
+ //update pipline until number of samlples is reached
+ for (int j=0; jm_NumberOfSamples->value(); j++) {myEvaluationFilter->Update();}
+ }
+
+ myEvaluationFilter->SetInput(NULL);
+ myPlayer=NULL;
+ EvaluationDataCollection.push_back(myEvaluationFilter);
+ }
+
+ //calculation and writing of output data
+ //open output file
+ m_CurrentAngleDifferencesWriteFile.open(std::string((m_Controls->m_OutputFilename->text() + ".angledifferences.csv").toUtf8()).c_str(), std::ios::out);
+ if (m_CurrentAngleDifferencesWriteFile.bad())
+ {
+ MessageBox("Error: Can't open output file for angle differences calculation!");
+ return;
+ }
+ //write header
+ WriteDifferenceAnglesHeader();
+ //compute angle differences
+ QString pos1 = "invalid";
+ QString pos2 = "invalid";
+ //now iterate through all evaluation data and calculate the angles
+ for(int i=0; im_DifferencesSLERP->isChecked())
+ {
+ //compute slerp average
+ q1 = GetSLERPAverage(EvaluationDataCollection.at(i));
+ q2 = GetSLERPAverage(EvaluationDataCollection.at(j));
+ }
+ else
+ {
+ //compute arithmetic average
+ q1 = EvaluationDataCollection.at(i)->GetQuaternionMean(0);
+ q2 = EvaluationDataCollection.at(j)->GetQuaternionMean(0);
+ }
+
+
+ double AngleBetweenTwoQuaternions = GetAngleBetweenTwoQuaterions(q1,q2);
+
+ //write data set
+ WriteDifferenceAnglesDataSet(pos1.toStdString(),pos2.toStdString(),i,j,AngleBetweenTwoQuaternions);
+ }
+ }
+
+ //close output file
+ m_CurrentAngleDifferencesWriteFile.close();
+ }
+
+void QmitkIGTTrackingDataEvaluationView::WriteDifferenceAnglesHeader()
+ {
+ m_CurrentAngleDifferencesWriteFile << "Name;Idx1;Idx2;Angle [-PI..+PI]; Angle [Degree]\n";
+ }
+
+void QmitkIGTTrackingDataEvaluationView::WriteDifferenceAnglesDataSet(std::string pos1, std::string pos2, int idx1, int idx2, double angle)
+ {
+ double PI = 3.1415926535897932384626433832795;
+ double angle_degree = (angle / PI) * 180;
+ m_CurrentAngleDifferencesWriteFile << "Angle between " << pos1 << " and " << pos2 << ";" << idx1 << ";" << idx2 << ";" << angle << ";" << angle_degree << "\n";
+ }
+
+double QmitkIGTTrackingDataEvaluationView::GetAngleBetweenTwoQuaterions(mitk::Quaternion a, mitk::Quaternion b)
+ {
+ double returnValue;
+
+ /*
+ //another variant
+ mitk::Quaternion combinedRotation = b * a;
+
+ itk::Vector pt1; //caution 5D-Tools: Vector must lie in the YZ-plane for a correct result.
+ pt1[0] = 0.0;
+ pt1[1] = 0.0;
+ pt1[2] = 100000.0;
+
+ itk::Matrix rotMatrixA;
+ for(int i=0; i<3; i++) for(int j=0; j<3; j++) rotMatrixA[i][j] = combinedRotation.rotation_matrix_transpose().transpose()[i][j];
+ itk::Vector pt2 = rotMatrixA*pt1;
+
+ //compute angle between the two vectors
+ returnValue = (pt1[0]*pt2[0]+pt1[1]*pt2[1]+pt1[2]*pt2[2]) / ( sqrt(pow(pt1[0],2)+pow(pt1[1],2)+pow(pt1[2],2)) * sqrt(pow(pt2[0],2)+pow(pt2[1],2)+pow(pt2[2],2)));
+ returnValue = acos(returnValue);
+
+ */
+
+ //variant with double precision
+
+ itk::Vector point; //caution 5D-Tools: Vector must lie in the YZ-plane for a correct result.
+ //TODO: welchen Vektor hier nehmen?
+ point[0] = 0;
+ point[1] = 100000.0;
+ point[2] = 100000.0;
+
+ //OB DAS HILFT?
+ a.normalize();
+ b.normalize();
+
+ itk::Matrix rotMatrixA;
+ for(int i=0; i<3; i++) for(int j=0; j<3; j++) rotMatrixA[i][j] = a.rotation_matrix_transpose().transpose()[i][j];
+
+ itk::Matrix rotMatrixB;
+ for(int i=0; i<3; i++) for(int j=0; j<3; j++) rotMatrixB[i][j] = b.rotation_matrix_transpose().transpose()[i][j];
+
+ itk::Vector pt1 = rotMatrixA * point;
+ itk::Vector pt2 = rotMatrixB * point;
+
+ returnValue = (pt1[0]*pt2[0]+pt1[1]*pt2[1]+pt1[2]*pt2[2]) / ( sqrt(pow(pt1[0],2.0)+pow(pt1[1],2.0)+pow(pt1[2],2.0)) * sqrt(pow(pt2[0],2.0)+pow(pt2[1],2.0)+pow(pt2[2],2.0)));
+ returnValue = acos(returnValue);
+
+
+ /* same code with float precision:
+ mitk::Point3D point;
+ mitk::FillVector3D(point,0,0,100); //caution 5D-Tools: Vector must lie in the YZ-plane for a correct result.
+ vnl_vector pt1 = a.rotate(point.Get_vnl_vector());
+ vnl_vector pt2 = b.rotate(point.Get_vnl_vector());
+
+ //compute angle between the two vectors
+ returnValue = (pt1[0]*pt2[0]+pt1[1]*pt2[1]+pt1[2]*pt2[2]) / ( sqrt(pow(pt1[0],2)+pow(pt1[1],2)+pow(pt1[2],2)) * sqrt(pow(pt2[0],2)+pow(pt2[1],2)+pow(pt2[2],2)));
+ returnValue = acos(returnValue);
+ //angle(pt1,pt2);
+ */
+
+
+ return returnValue;
+ }
+
+std::vector QmitkIGTTrackingDataEvaluationView::GetNavigationDatasFromFile(std::string filename)
+{
+std::vector returnValue = std::vector();
+std::vector fileContentLineByLine = GetFileContentLineByLine(filename);
+for(int i=1; i QmitkIGTTrackingDataEvaluationView::GetFileContentLineByLine(std::string filename)
+{
+std::vector readData = std::vector();
+
+//save old locale
+char * oldLocale;
+oldLocale = setlocale( LC_ALL, 0 );
+
+//define own locale
+std::locale C("C");
+setlocale( LC_ALL, "C" );
+
+//read file
+std::ifstream file;
+file.open(filename.c_str(), std::ios::in);
+if (file.good())
+ {
+ //read out file
+ file.seekg(0L, std::ios::beg); // move to begin of file
+ while (! file.eof())
+ {
+ std::string buffer;
+ std::getline(file,buffer); // read out file line by line
+ if (buffer.size() > 0) readData.push_back(buffer);
+
+ }
+ }
+
+file.close();
+
+//switch back to old locale
+setlocale( LC_ALL, oldLocale );
+
+return readData;
+}
+
+mitk::NavigationData::Pointer QmitkIGTTrackingDataEvaluationView::GetNavigationDataOutOfOneLine(std::string line)
+{
+ mitk::NavigationData::Pointer returnValue = mitk::NavigationData::New();
+
+ QString myLine = QString(line.c_str());
+
+ QStringList myLineList = myLine.split(';');
+
+ mitk::Point3D position;
+ mitk::Quaternion orientation;
+
+ double time = myLineList.at(1).toDouble();
+
+ bool valid = false;
+ if (myLineList.at(2).toStdString() == "1") valid = true;
+
+ position[0] = myLineList.at(3).toDouble();
+ position[1] = myLineList.at(4).toDouble();
+ position[2] = myLineList.at(5).toDouble();
+
+ orientation[0] = myLineList.at(6).toDouble();
+ orientation[1] = myLineList.at(7).toDouble();
+ orientation[2] = myLineList.at(8).toDouble();
+ orientation[3] = myLineList.at(9).toDouble();
+
+ //returnValue->SetTimeStamp(time);
+ returnValue->SetDataValid(valid);
+ returnValue->SetPosition(position);
+ returnValue->SetOrientation(orientation);
+
+ return returnValue;
+}
+
+mitk::Quaternion QmitkIGTTrackingDataEvaluationView::GetSLERPAverage(mitk::NavigationDataEvaluationFilter::Pointer evaluationFilter)
+{
+ mitk::Quaternion average;
+
+
+
+ //build a vector of quaternions from the evaulation filter (caution always takes the first (0) input of the filter
+ std::vector quaternions = std::vector();
+ for(int i=0; iGetNumberOfAnalysedNavigationData(0); i++)
+ {
+
+ mitk::Quaternion currentq = evaluationFilter->GetLoggedOrientation(i,0);
+
+ quaternions.push_back(currentq);
+ }
+
+ //compute the slerp average using the quaternion averaging class
+ mitk::QuaternionAveraging::Pointer myAverager = mitk::QuaternionAveraging::New();
+ average = myAverager->CalcAverage(quaternions);
+
+
+ return average;
+
+}
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingDataEvaluationView.h b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingDataEvaluationView.h
new file mode 100644
index 0000000000..a928c578b3
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingDataEvaluationView.h
@@ -0,0 +1,125 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date$
+Version: $Revision$
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+#ifndef QmitkIGTTrackingDataEvaluationView_h
+#define QmitkIGTTrackingDataEvaluationView_h
+
+#include
+
+#include
+
+#include "ui_QmitkIGTTrackingDataEvaluationViewControls.h"
+
+#include
+
+
+
+/*!
+ \brief QmitkIGTTrackingDataEvaluationView
+
+ \warning This application module is not yet documented. Use "svn blame/praise/annotate" and ask the author to provide basic documentation.
+
+ \sa QmitkFunctionality
+ \ingroup Functionalities
+*/
+class QmitkIGTTrackingDataEvaluationView : public QmitkFunctionality
+{
+ // this is needed for all Qt objects that should have a Qt meta-object
+ // (everything that derives from QObject and wants to have signal/slots)
+ Q_OBJECT
+
+ public:
+
+ static const std::string VIEW_ID;
+
+ QmitkIGTTrackingDataEvaluationView();
+ virtual ~QmitkIGTTrackingDataEvaluationView();
+
+ virtual void CreateQtPartControl(QWidget *parent);
+
+ virtual void StdMultiWidgetAvailable (QmitkStdMultiWidget &stdMultiWidget);
+ virtual void StdMultiWidgetNotAvailable();
+
+ protected slots:
+
+ void OnLoadFileList();
+ void OnAddToCurrentList();
+ void OnEvaluateData();
+ void OnGeneratePointSet();
+ void OnGeneratePointSetsOfSinglePositions();
+ void OnGenerateRotationLines();
+ void OnGenerateGroundTruthPointSet();
+ void OnConvertCSVtoXMLFile();
+ void OnCSVtoXMLLoadInputList();
+ void OnCSVtoXMLLoadOutputList();
+
+ /** Reads in exactly three position files als reference. */
+ void OnOrientationCalculation_CalcRef();
+ /** Uses always three positions (1,2,3: first orientation; 4,5,6: second orientation; and so on) in every file to calcualte a orientation. */
+ void OnOrientationCalculation_CalcOrientandWriteToFile();
+
+
+ protected:
+
+ Ui::QmitkIGTTrackingDataEvaluationViewControls* m_Controls;
+
+ QmitkStdMultiWidget* m_MultiWidget;
+
+ std::vector m_FilenameVector;
+
+ void MessageBox(std::string s);
+
+ std::fstream m_CurrentWriteFile;
+ void WriteHeader();
+ void WriteDataSet(mitk::NavigationDataEvaluationFilter::Pointer evaluationFilter, std::string dataSetName);
+
+ //members for orientation calculation
+ mitk::Point3D m_RefPoint1;
+ mitk::Point3D m_RefPoint2;
+ mitk::Point3D m_RefPoint3;
+
+ double m_scalingfactor; //scaling factor for visualization, 1 by default
+
+ //angle diffrences: seperated file
+ std::fstream m_CurrentAngleDifferencesWriteFile;
+ void CalculateDifferenceAngles();
+ void WriteDifferenceAnglesHeader();
+ void WriteDifferenceAnglesDataSet(std::string pos1, std::string pos2, int idx1, int idx2, double angle);
+
+ //different help methods to read a csv logging file
+ std::vector GetNavigationDatasFromFile(std::string filename);
+ std::vector GetFileContentLineByLine(std::string filename);
+ mitk::NavigationData::Pointer GetNavigationDataOutOfOneLine(std::string line);
+
+ //CSV to XML members
+ std::vector m_CSVtoXMLInputFilenameVector;
+ std::vector m_CSVtoXMLOutputFilenameVector;
+
+ //returns the number of converted lines
+ int ConvertOneFile(std::string inputFilename, std::string outputFilename);
+
+ /** @brief calculates the angle in the plane perpendicular to the rotation axis of the two quaterions. */
+ double GetAngleBetweenTwoQuaterions(mitk::Quaternion a, mitk::Quaternion b);
+
+ /** @brief calculates the slerp average of a set of quaternions which is stored in the navigation data evaluation filter */
+ mitk::Quaternion GetSLERPAverage(mitk::NavigationDataEvaluationFilter::Pointer);
+};
+
+
+
+#endif // _QMITKIGTTRACKINGDATAEVALUATIONVIEW_H_INCLUDED
+
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingDataEvaluationViewControls.ui b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingDataEvaluationViewControls.ui
new file mode 100644
index 0000000000..9095e6a5af
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingDataEvaluationViewControls.ui
@@ -0,0 +1,941 @@
+
+
+ QmitkIGTTrackingDataEvaluationViewControls
+
+
+
+ 0
+ 0
+ 378
+ 955
+
+
+
+
+ 0
+ 0
+
+
+
+ QmitkTemplate
+
+
+ -
+
+
+ 0
+
+
+
+ Evaluation
+
+
+
-
+
+
+ Input File List (recorded NavigationData / *.csv):
+
+
+
+ -
+
+
+ -
+
+
-
+
+
+ Add Path:
+
+
+
+ -
+
+
+ D:/Experimente/
+
+
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-style:italic;">(use text files with a complete filename in every line)</span></p></body></html>
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 110
+ 0
+
+
+
+ Load New List
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 110
+ 0
+
+
+
+ Add To Current List
+
+
+
+
+
+ -
+
+
-
+
+
+ Number of samples to analyze:
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ 1000000
+
+
+ 150
+
+
+
+
+
+ -
+
+
-
+
+
+ Result CSV Filename:
+
+
+
+ -
+
+
+ D:/Experimente/output.csv
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 200
+ 50
+
+
+
+
+ 200
+ 50
+
+
+
+ START EVALUATION
+
+
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ -
+
+
+ Visualization Tools:
+
+
+
+ -
+
+
+ Generate PointSet of Mean Positions
+
+
+
+ -
+
+
+ Generate PointSets of Single Positions
+
+
+
+ -
+
+
+ Generate Lines for Rotation
+
+
+
+ -
+
+
-
+
+
+ Prefix for Data Nodes:
+
+
+
+ -
+
+
+
+
+ -
+
+
+ Qt::Vertical
+
+
+ QSizePolicy::Expanding
+
+
+
+ 20
+ 220
+
+
+
+
+
+
+
+
+ Settings
+
+
+ -
+
+
+ Output Settings
+
+
+
-
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-style:italic; text-decoration: underline;">Position</span></p></body></html>
+
+
+
+ -
+
+
+ Mean (x,y,z)
+
+
+
+ -
+
+
+ Standard Deviation (x,y,z)
+
+
+
+ -
+
+
+ Sample Standard Deviation (x,y,z)
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-style:italic; text-decoration: underline;">Orientation</span></p></body></html>
+
+
+
+ -
+
+
+ Quaternion Mean (qx,qy,qz,qr)
+
+
+
+ -
+
+
+ Quaternion Mean (SLERP)
+
+
+
+ -
+
+
+ Quaternion Standard Deviation (qx,qy,qz,qr)
+
+
+
+ -
+
+
+ Euler Mean (tx,ty,tz)
+
+
+
+ -
+
+
+ Difference Angles to all other Positions
+
+
+
+ -
+
+
+ Difference Angles to all other Positions (SLERP)
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-style:italic; text-decoration: underline;">Position Error</span></p></body></html>
+
+
+
+ -
+
+
+ Mean
+
+
+
+ -
+
+
+ Standard Deviation
+
+
+
+ -
+
+
+ Sample Standard Deviation
+
+
+
+ -
+
+
+ RMS
+
+
+
+ -
+
+
+ Median
+
+
+
+ -
+
+
+ Min/Max
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-style:italic; text-decoration: underline;">Orientation Error</span></p></body></html>
+
+
+
+ -
+
+
+ Euler RMS
+
+
+
+
+
+
+ -
+
+
-
+
+
+ Scaling Factor for Visualization:
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ 1.000000000000000
+
+
+
+
+
+ -
+
+
+ Qt::Vertical
+
+
+
+ 20
+ 344
+
+
+
+
+
+
+
+
+ Tools
+
+
+ -
+
+
+ Point Set Ground Truth Generator
+
+
+
-
+
+
-
+
+
+ Generate
+
+
+
+ -
+
+
+ 1
+
+
+ 999
+
+
+ 10
+
+
+
+ -
+
+
+ X
+
+
+
+ -
+
+
+ 1
+
+
+ 999
+
+
+ 9
+
+
+
+ -
+
+
+ Point Set
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+
+
+ -
+
+
-
+
+
+ Inter Point Distance (in mm):
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ 1
+
+
+ 99999
+
+
+ 50
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Generate
+
+
+
+
+
+
+
+
+ -
+
+
+ Result CSV File to NavigationData Converter
+
+
+
-
+
+
+ Convert Single File
+
+
+ true
+
+
+
+ -
+
+
+ Input CSV Logging File:
+
+
+
+ -
+
+
+ C:/Tools/test.csv
+
+
+
+ -
+
+
+ Output Navigation Data File:
+
+
+
+ -
+
+
+ C:/Tools/testoutput.xml
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ -
+
+
+ Convert File List
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-style:italic;">(use text files with a complete filename in every line)</span></p></body></html>
+
+
+
+ -
+
+
-
+
+
+ not loaded
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 100
+ 0
+
+
+
+ Load Input List
+
+
+
+
+
+ -
+
+
-
+
+
+ not loaded
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 100
+ 0
+
+
+
+ Load Output List
+
+
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Output Format
+
+
+
-
+
+
+ XML
+
+
+ true
+
+
+
+ -
+
+
+ CSV
+
+
+
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Convert
+
+
+
+
+
+
+
+
+ -
+
+
+ Orientation Calculation (out of three positions)
+
+
+
-
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Generate Reference From Current List
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Write Orientation Quaternions To File
+
+
+
+
+
+
+
+
+ -
+
+
+ Qt::Vertical
+
+
+
+ 20
+ 632
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.cpp b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.cpp
new file mode 100644
index 0000000000..9a0c25196c
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.cpp
@@ -0,0 +1,667 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date$
+Version: $Revision$
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+// Blueberry
+#include
+#include
+
+// Qmitk
+#include "QmitkIGTTrackingSemiAutomaticMeasurementView.h"
+#include "QmitkStdMultiWidget.h"
+
+// Qt
+#include
+#include
+#include
+#include
+
+// MITK
+#include
+#include
+#include
+
+// POCO
+#include
+#include
+
+const std::string QmitkIGTTrackingSemiAutomaticMeasurementView::VIEW_ID = "org.mitk.views.igttrackingsemiautomaticmeasurement";
+
+QmitkIGTTrackingSemiAutomaticMeasurementView::QmitkIGTTrackingSemiAutomaticMeasurementView()
+ : QmitkFunctionality()
+ , m_Controls(0)
+ , m_MultiWidget(NULL)
+{
+ m_NextFile = 0;
+ m_FilenameVector = std::vector();
+ m_Timer = new QTimer(this);
+ m_logging = false;
+ m_referenceValid = true;
+ m_tracking = false;
+ m_EvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+}
+
+QmitkIGTTrackingSemiAutomaticMeasurementView::~QmitkIGTTrackingSemiAutomaticMeasurementView()
+{
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::CreateResults()
+{
+ QString LogFileName = m_Controls->m_OutputPath->text() + "_results.log";
+ mitk::LoggingBackend::Unregister();
+ mitk::LoggingBackend::SetLogFile(LogFileName.toStdString().c_str());
+ mitk::LoggingBackend::Register();
+
+ double RMSmean = 0;
+ for (int i = 0; i < m_RMSValues.size(); i++)
+ {
+ MITK_INFO << "RMS at " << this->m_FilenameVector.at(i) << ": " << m_RMSValues.at(i);
+ RMSmean += m_RMSValues.at(i);
+ }
+ RMSmean /= m_RMSValues.size();
+ MITK_INFO << "RMS mean over " << m_RMSValues.size() << " values: " << RMSmean;
+
+ mitk::DataNode::Pointer newNode = mitk::DataNode::New();
+ newNode->SetName("Tracking Results");
+ newNode->SetData(this->m_MeanPoints);
+ this->GetDataStorage()->Add(newNode);
+
+ if (m_MeanPoints->GetSize() == 12)
+ {
+ MITK_INFO << "Computing distance error for Compact FG...";
+ std::vector distances;
+ mitk::Point3D test;
+ //row 1
+ distances.push_back(m_MeanPoints->GetPoint(0).EuclideanDistanceTo(m_MeanPoints->GetPoint(1))); //0
+ distances.push_back(m_MeanPoints->GetPoint(1).EuclideanDistanceTo(m_MeanPoints->GetPoint(2))); //1
+ distances.push_back(m_MeanPoints->GetPoint(2).EuclideanDistanceTo(m_MeanPoints->GetPoint(3))); //2
+ //row 2
+ distances.push_back(m_MeanPoints->GetPoint(4).EuclideanDistanceTo(m_MeanPoints->GetPoint(5))); //3
+ distances.push_back(m_MeanPoints->GetPoint(5).EuclideanDistanceTo(m_MeanPoints->GetPoint(6))); //4
+ distances.push_back(m_MeanPoints->GetPoint(6).EuclideanDistanceTo(m_MeanPoints->GetPoint(7))); //5
+ //row 3
+ distances.push_back(m_MeanPoints->GetPoint(8).EuclideanDistanceTo(m_MeanPoints->GetPoint(9))); //6
+ distances.push_back(m_MeanPoints->GetPoint(9).EuclideanDistanceTo(m_MeanPoints->GetPoint(10))); //7
+ distances.push_back(m_MeanPoints->GetPoint(10).EuclideanDistanceTo(m_MeanPoints->GetPoint(11))); //8
+ //column 1
+ distances.push_back(m_MeanPoints->GetPoint(0).EuclideanDistanceTo(m_MeanPoints->GetPoint(4))); //9
+ distances.push_back(m_MeanPoints->GetPoint(4).EuclideanDistanceTo(m_MeanPoints->GetPoint(8))); //10
+ //column 2
+ distances.push_back(m_MeanPoints->GetPoint(1).EuclideanDistanceTo(m_MeanPoints->GetPoint(5))); //11
+ distances.push_back(m_MeanPoints->GetPoint(5).EuclideanDistanceTo(m_MeanPoints->GetPoint(9))); //12
+ //column 3
+ distances.push_back(m_MeanPoints->GetPoint(2).EuclideanDistanceTo(m_MeanPoints->GetPoint(6))); //13
+ distances.push_back(m_MeanPoints->GetPoint(6).EuclideanDistanceTo(m_MeanPoints->GetPoint(10))); //14
+ //column 4
+ distances.push_back(m_MeanPoints->GetPoint(3).EuclideanDistanceTo(m_MeanPoints->GetPoint(7))); //15
+ distances.push_back(m_MeanPoints->GetPoint(7).EuclideanDistanceTo(m_MeanPoints->GetPoint(11))); //16
+
+ std::vector errors;
+ double meanError = 0;
+ for (int i = 0; i < distances.size(); i++)
+ {
+ double currentError = distances.at(i) - (double)50.0;
+ errors.push_back(currentError);
+ meanError += currentError;
+ MITK_INFO << "Error" << i << " : " << currentError;
+ }
+ meanError /= distances.size();
+ MITK_INFO << "Mean error : " << meanError;
+ }
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::CreateQtPartControl(QWidget *parent)
+{
+ // build up qt view, unless already done
+ if (!m_Controls)
+ {
+ // create GUI widgets from the Qt Designer's .ui file
+ m_Controls = new Ui::QmitkIGTTrackingSemiAutomaticMeasurementViewControls;
+ m_Controls->setupUi(parent);
+
+ //buttons
+ connect(m_Controls->m_LoadMeasurementToolStorage, SIGNAL(clicked()), this, SLOT(OnLoadMeasurementStorage()));
+ connect(m_Controls->m_LoadReferenceToolStorage, SIGNAL(clicked()), this, SLOT(OnLoadReferenceStorage()));
+ connect(m_Controls->m_StartTracking, SIGNAL(clicked()), this, SLOT(OnStartTracking()));
+ connect(m_Controls->m_LoadList, SIGNAL(clicked()), this, SLOT(OnMeasurementLoadFile()));
+ connect(m_Controls->m_StartNextMeasurement, SIGNAL(clicked()), this, SLOT(StartNextMeasurement()));
+ connect(m_Controls->m_ReapeatLastMeasurement, SIGNAL(clicked()), this, SLOT(RepeatLastMeasurement()));
+ connect(m_Controls->m_SetReference, SIGNAL(clicked()), this, SLOT(OnSetReference()));
+ connect(m_Controls->m_UseReferenceTrackingSystem, SIGNAL(toggled(bool)), this, SLOT(OnUseReferenceToggled(bool)));
+ connect(m_Controls->m_CreateResults, SIGNAL(clicked()), this, SLOT(CreateResults()));
+
+ //event filter
+ qApp->installEventFilter(this);
+
+ //timers
+ connect(m_Timer, SIGNAL(timeout()), this, SLOT(UpdateTimer()));
+ }
+
+ //initialize some view
+ m_Controls->m_StopTracking->setEnabled(false);
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::StdMultiWidgetAvailable(QmitkStdMultiWidget &stdMultiWidget)
+{
+ m_MultiWidget = &stdMultiWidget;
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnUseReferenceToggled(bool state)
+{
+ if (state)
+ {
+ m_Controls->m_ReferenceBox->setEnabled(true);
+ m_Controls->m_SetReference->setEnabled(true);
+ }
+
+ else
+ {
+ m_Controls->m_ReferenceBox->setEnabled(false);
+ m_Controls->m_SetReference->setEnabled(false);
+ }
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::StdMultiWidgetNotAvailable()
+{
+ m_MultiWidget = NULL;
+}
+
+mitk::NavigationToolStorage::Pointer QmitkIGTTrackingSemiAutomaticMeasurementView::ReadStorage(std::string file)
+{
+ mitk::NavigationToolStorage::Pointer returnValue;
+
+ //initialize tool storage
+ returnValue = mitk::NavigationToolStorage::New();
+
+ //read tool storage from disk
+ mitk::NavigationToolStorageDeserializer::Pointer myDeserializer = mitk::NavigationToolStorageDeserializer::New(GetDataStorage());
+ returnValue = myDeserializer->Deserialize(file);
+ if (returnValue.IsNull())
+ {
+ QMessageBox msgBox;
+ msgBox.setText(myDeserializer->GetErrorMessage().c_str());
+ msgBox.exec();
+
+ returnValue = NULL;
+ }
+
+ return returnValue;
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnSetReference()
+{
+ //initialize reference
+ m_ReferenceStartPositions = std::vector();
+ m_ReferenceTrackingDeviceSource->Update();
+ QString Label = "Positions At Start: ";
+ for (int i = 0; i < m_ReferenceTrackingDeviceSource->GetNumberOfOutputs(); i++)
+ {
+ mitk::Point3D position = m_ReferenceTrackingDeviceSource->GetOutput(i)->GetPosition();
+ Label = Label + "Tool" + QString::number(i) + ":[" + QString::number(position[0]) + ":" + QString::number(position[1]) + ":" + QString::number(position[1]) + "] ";
+ m_ReferenceStartPositions.push_back(position);
+ }
+ m_Controls->m_ReferencePosAtStart->setText(Label);
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnLoadMeasurementStorage()
+{
+ //read in filename
+ QString filename = QFileDialog::getOpenFileName(NULL, tr("Open Toolfile"), "/", tr("All Files (*.*)"));
+ if (filename.isNull()) return;
+
+ m_MeasurementStorage = ReadStorage(filename.toStdString());
+
+ //update label
+ Poco::Path myPath = Poco::Path(filename.toStdString()); //use this to seperate filename from path
+ QString toolLabel = QString("Tool Storage: ") + QString::number(m_MeasurementStorage->GetToolCount()) + " Tools from " + myPath.getFileName().c_str();
+ m_Controls->m_MeasurementToolStorageLabel->setText(toolLabel);
+
+ //update status widget
+ m_Controls->m_ToolStatusWidget->RemoveStatusLabels();
+ m_Controls->m_ToolStatusWidget->PreShowTools(m_MeasurementStorage);
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnLoadReferenceStorage()
+{
+ //read in filename
+ static QString oldFile;
+ if (oldFile.isNull()) oldFile = "/";
+ QString filename = QFileDialog::getOpenFileName(NULL, tr("Open Toolfile"), oldFile, tr("All Files (*.*)"));
+ if (filename.isNull()) return;
+ oldFile = filename;
+
+ m_ReferenceStorage = ReadStorage(filename.toStdString());
+
+ //update label
+ Poco::Path myPath = Poco::Path(filename.toStdString()); //use this to seperate filename from path
+ QString toolLabel = QString("Tool Storage: ") + QString::number(m_ReferenceStorage->GetToolCount()) + " Tools from " + myPath.getFileName().c_str();
+ m_Controls->m_ReferenceToolStorageLabel->setText(toolLabel);
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnStartTracking()
+{
+ //check if everything is ready to start tracking
+ if (m_MeasurementStorage.IsNull())
+ {
+ MessageBox("Error: No measurement tools loaded yet!");
+ return;
+ }
+ else if (m_ReferenceStorage.IsNull() && m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ MessageBox("Error: No refernce tools loaded yet!");
+ return;
+ }
+ else if (m_MeasurementStorage->GetToolCount() == 0)
+ {
+ MessageBox("Error: No way to track without tools!");
+ return;
+ }
+ else if (m_Controls->m_UseReferenceTrackingSystem->isChecked() && (m_ReferenceStorage->GetToolCount() == 0))
+ {
+ MessageBox("Error: No way to track without tools!");
+ return;
+ }
+
+ //build the first IGT pipeline (MEASUREMENT)
+ mitk::TrackingDeviceSourceConfigurator::Pointer myTrackingDeviceSourceFactory1 = mitk::TrackingDeviceSourceConfigurator::New(this->m_MeasurementStorage, this->m_Controls->m_MeasurementTrackingDeviceConfigurationWidget->GetTrackingDevice());
+ m_MeasurementTrackingDeviceSource = myTrackingDeviceSourceFactory1->CreateTrackingDeviceSource(this->m_MeasurementToolVisualizationFilter);
+ if (m_MeasurementTrackingDeviceSource.IsNull())
+ {
+ MessageBox(myTrackingDeviceSourceFactory1->GetErrorMessage());
+ return;
+ }
+ //connect the tool visualization widget
+ for (int i = 0; i < m_MeasurementTrackingDeviceSource->GetNumberOfOutputs(); i++)
+ {
+ m_Controls->m_ToolStatusWidget->AddNavigationData(m_MeasurementTrackingDeviceSource->GetOutput(i));
+ m_EvaluationFilter->SetInput(i, m_MeasurementTrackingDeviceSource->GetOutput(i));
+ }
+ m_Controls->m_ToolStatusWidget->ShowStatusLabels();
+ m_Controls->m_ToolStatusWidget->SetShowPositions(true);
+ m_Controls->m_ToolStatusWidget->SetShowQuaternions(true);
+
+ //build the second IGT pipeline (REFERENCE)
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ mitk::TrackingDeviceSourceConfigurator::Pointer myTrackingDeviceSourceFactory2 = mitk::TrackingDeviceSourceConfigurator::New(this->m_ReferenceStorage, this->m_Controls->m_ReferenceDeviceConfigurationWidget->GetTrackingDevice());
+ m_ReferenceTrackingDeviceSource = myTrackingDeviceSourceFactory2->CreateTrackingDeviceSource();
+ if (m_ReferenceTrackingDeviceSource.IsNull())
+ {
+ MessageBox(myTrackingDeviceSourceFactory2->GetErrorMessage());
+ return;
+ }
+ }
+
+ //initialize tracking
+ try
+ {
+ m_MeasurementTrackingDeviceSource->Connect();
+ m_MeasurementTrackingDeviceSource->StartTracking();
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ m_ReferenceTrackingDeviceSource->Connect();
+ m_ReferenceTrackingDeviceSource->StartTracking();
+ }
+ }
+ catch (...)
+ {
+ MessageBox("Error while starting the tracking device!");
+ return;
+ }
+
+ //set reference
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked()) OnSetReference();
+
+ //start timer
+ m_Timer->start(1000 / (m_Controls->m_SamplingRate->value()));
+
+ m_Controls->m_StartTracking->setEnabled(false);
+ m_Controls->m_StartTracking->setEnabled(true);
+
+ m_tracking = true;
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnStopTracking()
+{
+ if (this->m_logging) FinishMeasurement();
+ m_MeasurementTrackingDeviceSource->Disconnect();
+ m_MeasurementTrackingDeviceSource->StopTracking();
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ m_ReferenceTrackingDeviceSource->Disconnect();
+ m_ReferenceTrackingDeviceSource->StopTracking();
+ }
+ m_Timer->stop();
+ m_Controls->m_StartTracking->setEnabled(true);
+ m_Controls->m_StartTracking->setEnabled(false);
+ m_tracking = false;
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::OnMeasurementLoadFile()
+{
+ m_FilenameVector = std::vector();
+ m_FilenameVector.clear();
+ m_NextFile = 0;
+
+ //read in filename
+ QString filename = QFileDialog::getOpenFileName(NULL, tr("Open Measurement Filename List"), "/", tr("All Files (*.*)"));
+ if (filename.isNull()) return;
+
+ //save old locale
+ char * oldLocale;
+ oldLocale = setlocale(LC_ALL, 0);
+
+ //define own locale
+ std::locale C("C");
+ setlocale(LC_ALL, "C");
+
+ //read file
+ std::ifstream file;
+ file.open(filename.toStdString().c_str(), std::ios::in);
+ if (file.good())
+ {
+ //read out file
+ file.seekg(0L, std::ios::beg); // move to begin of file
+ while (!file.eof())
+ {
+ std::string buffer;
+ std::getline(file, buffer); // read out file line by line
+ if (buffer.size() > 0) m_FilenameVector.push_back(buffer);
+ }
+ }
+
+ //fill list at GUI
+ m_Controls->m_MeasurementList->clear();
+ for (unsigned int i = 0; i < m_FilenameVector.size(); i++) { new QListWidgetItem(tr(m_FilenameVector.at(i).c_str()), m_Controls->m_MeasurementList); }
+
+ //update label next measurement
+ std::stringstream label;
+ label << "Next Measurement: " << m_FilenameVector.at(0);
+ m_Controls->m_NextMeasurement->setText(label.str().c_str());
+
+ //reset results files
+ m_MeanPoints = mitk::PointSet::New();
+ m_RMSValues = std::vector();
+ m_EvaluationFilter = mitk::NavigationDataEvaluationFilter::New();
+ if (m_MeasurementToolVisualizationFilter.IsNotNull()) m_EvaluationFilter->SetInput(0, m_MeasurementToolVisualizationFilter->GetOutput(0));
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::UpdateTimer()
+{
+ if (m_EvaluationFilter.IsNotNull() && m_logging) m_EvaluationFilter->Update();
+ else m_MeasurementToolVisualizationFilter->Update();
+
+ m_Controls->m_ToolStatusWidget->Refresh();
+
+ //update reference
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ m_ReferenceTrackingDeviceSource->Update();
+ QString Label = "Current Positions: ";
+ bool distanceThresholdExceeded = false;
+ for (int i = 0; i < m_ReferenceTrackingDeviceSource->GetNumberOfOutputs(); i++)
+ {
+ mitk::Point3D position = m_ReferenceTrackingDeviceSource->GetOutput(i)->GetPosition();
+ Label = Label + "Tool" + QString::number(i) + ":[" + QString::number(position[0]) + ":" + QString::number(position[1]) + ":" + QString::number(position[1]) + "] ";
+ if (position.EuclideanDistanceTo(m_ReferenceStartPositions.at(i)) > m_Controls->m_ReferenceThreshold->value()) distanceThresholdExceeded = true;
+ }
+ m_Controls->m_ReferenceCurrentPos->setText(Label);
+ if (distanceThresholdExceeded)
+ {
+ m_Controls->m_ReferenceOK->setText("NOT OK!");
+ m_referenceValid = false;
+ }
+ else
+ {
+ m_Controls->m_ReferenceOK->setText("OK");
+ m_referenceValid = true;
+ }
+ }
+
+ //update logging
+ if (m_logging)
+ {
+ //check for missing objects
+ if (m_MeasurementLoggingFilterXML.IsNull() ||
+ m_MeasurementLoggingFilterCSV.IsNull()
+ )
+ {
+ return;
+ }
+
+ //log/measure
+ m_MeasurementLoggingFilterXML->Update();
+ m_MeasurementLoggingFilterCSV->Update();
+
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked() &&
+ m_ReferenceLoggingFilterXML.IsNotNull() &&
+ m_ReferenceLoggingFilterCSV.IsNotNull())
+ {
+ m_ReferenceLoggingFilterXML->Update();
+ m_ReferenceLoggingFilterCSV->Update();
+ }
+ m_loggedFrames++;
+ LogAdditionalCSVFile();
+
+ //check if all frames are logged ... if yes finish the measurement
+ if (m_loggedFrames > m_Controls->m_SamplesPerMeasurement->value()) { FinishMeasurement(); }
+
+ //update logging label
+ QString loggingLabel = "Collected Samples: " + QString::number(m_loggedFrames);
+ m_Controls->m_CollectedSamples->setText(loggingLabel);
+ }
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::StartNextMeasurement()
+{
+ if (this->m_NextFile >= m_FilenameVector.size())
+ {
+ MessageBox("Last Measurement reached!");
+ return;
+ }
+
+ m_loggedFrames = 0;
+ m_logging = true;
+
+ //check if directory exists, if not create one
+ Poco::File myPath(std::string(m_Controls->m_OutputPath->text().toUtf8()).c_str());
+ if (!myPath.exists()) myPath.createDirectory();
+
+ QString LogFileName = m_Controls->m_OutputPath->text() + QString(m_FilenameVector.at(m_NextFile).c_str()) + ".log";
+ mitk::LoggingBackend::Unregister();
+ mitk::LoggingBackend::SetLogFile(LogFileName.toStdString().c_str());
+ mitk::LoggingBackend::Register();
+
+ //initialize logging filters
+ m_MeasurementLoggingFilterXML = mitk::NavigationDataRecorderDeprecated::New();
+ m_MeasurementLoggingFilterXML->SetRecordingMode(mitk::NavigationDataRecorderDeprecated::NormalFile);
+ m_MeasurementLoggingFilterCSV = mitk::NavigationDataRecorderDeprecated::New();
+ m_MeasurementLoggingFilterCSV->SetRecordingMode(mitk::NavigationDataRecorderDeprecated::NormalFile);
+ m_MeasurementLoggingFilterXML->SetOutputFormat(mitk::NavigationDataRecorderDeprecated::xml);
+ m_MeasurementLoggingFilterCSV->SetOutputFormat(mitk::NavigationDataRecorderDeprecated::csv);
+ QString MeasurementFilenameXML = m_Controls->m_OutputPath->text() + QString(m_FilenameVector.at(m_NextFile).c_str()) + ".xml";
+ QString MeasurementFilenameCSV = m_Controls->m_OutputPath->text() + QString(m_FilenameVector.at(m_NextFile).c_str()) + ".csv";
+ m_MeasurementLoggingFilterXML->SetFileName(MeasurementFilenameXML.toStdString());
+ m_MeasurementLoggingFilterCSV->SetFileName(MeasurementFilenameCSV.toStdString());
+ m_MeasurementLoggingFilterXML->SetRecordCountLimit(m_Controls->m_SamplesPerMeasurement->value());
+ m_MeasurementLoggingFilterCSV->SetRecordCountLimit(m_Controls->m_SamplesPerMeasurement->value());
+
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ m_ReferenceLoggingFilterXML = mitk::NavigationDataRecorderDeprecated::New();
+ m_ReferenceLoggingFilterXML->SetRecordingMode(mitk::NavigationDataRecorderDeprecated::NormalFile);
+ m_ReferenceLoggingFilterCSV = mitk::NavigationDataRecorderDeprecated::New();
+ m_ReferenceLoggingFilterCSV->SetRecordingMode(mitk::NavigationDataRecorderDeprecated::NormalFile);
+ m_ReferenceLoggingFilterXML->SetOutputFormat(mitk::NavigationDataRecorderDeprecated::xml);
+ m_ReferenceLoggingFilterCSV->SetOutputFormat(mitk::NavigationDataRecorderDeprecated::csv);
+ QString ReferenceFilenameXML = m_Controls->m_OutputPath->text() + "Reference_" + QString(m_FilenameVector.at(m_NextFile).c_str()) + ".xml";
+ QString ReferenceFilenameCSV = m_Controls->m_OutputPath->text() + "Reference_" + QString(m_FilenameVector.at(m_NextFile).c_str()) + ".csv";
+ m_ReferenceLoggingFilterXML->SetFileName(ReferenceFilenameXML.toStdString());
+ m_ReferenceLoggingFilterCSV->SetFileName(ReferenceFilenameCSV.toStdString());
+ m_ReferenceLoggingFilterXML->SetRecordCountLimit(m_Controls->m_SamplesPerMeasurement->value());
+ m_ReferenceLoggingFilterCSV->SetRecordCountLimit(m_Controls->m_SamplesPerMeasurement->value());
+ }
+
+ //start additional csv logging
+ StartLoggingAdditionalCSVFile(m_FilenameVector.at(m_NextFile));
+
+ //connect filter
+ for (int i = 0; i < m_MeasurementToolVisualizationFilter->GetNumberOfOutputs(); i++)
+ {
+ m_MeasurementLoggingFilterXML->AddNavigationData(m_MeasurementToolVisualizationFilter->GetOutput(i));
+ m_MeasurementLoggingFilterCSV->AddNavigationData(m_MeasurementToolVisualizationFilter->GetOutput(i));
+ }
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked()) for (int i = 0; i < m_ReferenceTrackingDeviceSource->GetNumberOfOutputs(); i++)
+ {
+ m_ReferenceLoggingFilterXML->AddNavigationData(m_ReferenceTrackingDeviceSource->GetOutput(i));
+ m_ReferenceLoggingFilterCSV->AddNavigationData(m_ReferenceTrackingDeviceSource->GetOutput(i));
+ }
+
+ //start filter
+ m_MeasurementLoggingFilterXML->StartRecording();
+ m_MeasurementLoggingFilterCSV->StartRecording();
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ m_ReferenceLoggingFilterXML->StartRecording();
+ m_ReferenceLoggingFilterCSV->StartRecording();
+ }
+
+ //disable all buttons
+ DisableAllButtons();
+
+ //update label next measurement
+ std::stringstream label;
+ if ((m_NextFile + 1) >= m_FilenameVector.size()) label << "Next Measurement: ";
+ else label << "Next Measurement: " << m_FilenameVector.at(m_NextFile + 1);
+ m_Controls->m_NextMeasurement->setText(label.str().c_str());
+
+ //update label last measurement
+ std::stringstream label2;
+ label2 << "Last Measurement: " << m_FilenameVector.at(m_NextFile);
+ m_Controls->m_LastMeasurement->setText(label2.str().c_str());
+
+ m_NextFile++;
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::RepeatLastMeasurement()
+{
+ m_NextFile--;
+ StartNextMeasurement();
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::MessageBox(std::string s)
+{
+ QMessageBox msgBox;
+ msgBox.setText(s.c_str());
+ msgBox.exec();
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::DisableAllButtons()
+{
+ m_Controls->m_LoadList->setEnabled(false);
+ m_Controls->m_StartNextMeasurement->setEnabled(false);
+ m_Controls->m_ReapeatLastMeasurement->setEnabled(false);
+ m_Controls->m_SamplingRate->setEnabled(false);
+ m_Controls->m_SamplesPerMeasurement->setEnabled(false);
+ m_Controls->m_ReferenceThreshold->setEnabled(false);
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::EnableAllButtons()
+{
+ m_Controls->m_LoadList->setEnabled(true);
+ m_Controls->m_StartNextMeasurement->setEnabled(true);
+ m_Controls->m_ReapeatLastMeasurement->setEnabled(true);
+ m_Controls->m_SamplingRate->setEnabled(true);
+ m_Controls->m_SamplesPerMeasurement->setEnabled(true);
+ m_Controls->m_ReferenceThreshold->setEnabled(true);
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::FinishMeasurement()
+{
+ m_logging = false;
+
+ m_MeasurementLoggingFilterXML->StopRecording();
+ m_MeasurementLoggingFilterCSV->StopRecording();
+ if (m_Controls->m_UseReferenceTrackingSystem->isChecked())
+ {
+ m_ReferenceLoggingFilterXML->StopRecording();
+ m_ReferenceLoggingFilterCSV->StopRecording();
+ }
+ StopLoggingAdditionalCSVFile();
+
+ int id = m_NextFile - 1;
+ mitk::Point3D positionMean = m_EvaluationFilter->GetPositionMean(0);
+ MITK_INFO << "Evaluated " << m_EvaluationFilter->GetNumberOfAnalysedNavigationData(0) << " samples.";
+ double rms = m_EvaluationFilter->GetPositionErrorRMS(0);
+ MITK_INFO << "RMS: " << rms;
+ MITK_INFO << "Position Mean: " << positionMean;
+ m_MeanPoints->SetPoint(id, positionMean);
+ if (m_RMSValues.size() <= id) m_RMSValues.push_back(rms);
+ else m_RMSValues[id] = rms;
+
+ m_EvaluationFilter->ResetStatistic();
+
+ EnableAllButtons();
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::StartLoggingAdditionalCSVFile(std::string filePostfix)
+{
+ //write logfile header
+ QString header = "Nr;MITK_Time;Valid_Reference;";
+ QString tool = QString("MeasureTool_") + QString(m_MeasurementTrackingDeviceSource->GetOutput(0)->GetName());
+ header = header + tool + "[x];" + tool + "[y];" + tool + "[z];" + tool + "[qx];" + tool + "[qy];" + tool + "[qz];" + tool + "[qr]\n";
+
+ //open logfile and write header
+ m_logFileCSV.open(std::string(m_Controls->m_OutputPath->text().toUtf8()).append("/LogFileCombined").append(filePostfix.c_str()).append(".csv").c_str(), std::ios::out);
+ m_logFileCSV << header.toStdString().c_str();
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::LogAdditionalCSVFile()
+{
+ mitk::Point3D pos = m_MeasurementTrackingDeviceSource->GetOutput(0)->GetPosition();
+ mitk::Quaternion rot = m_MeasurementTrackingDeviceSource->GetOutput(0)->GetOrientation();
+ std::string valid = "";
+ if (m_referenceValid) valid = "true";
+ else valid = "false";
+ std::stringstream timestamp;
+ timestamp << m_MeasurementTrackingDeviceSource->GetOutput(0)->GetTimeStamp();
+ QString dataSet = QString::number(m_loggedFrames) + ";" + QString(timestamp.str().c_str()) + ";" + QString(valid.c_str()) + ";" + QString::number(pos[0]) + ";" + QString::number(pos[1]) + ";" + QString::number(pos[2]) + ";" + QString::number(rot.x()) + ";" + QString::number(rot.y()) + ";" + QString::number(rot.z()) + ";" + QString::number(rot.r()) + "\n";
+ m_logFileCSV << dataSet.toStdString();
+}
+
+void QmitkIGTTrackingSemiAutomaticMeasurementView::StopLoggingAdditionalCSVFile()
+{
+ m_logFileCSV.close();
+}
+
+bool QmitkIGTTrackingSemiAutomaticMeasurementView::eventFilter(QObject *obj, QEvent *ev)
+{
+ if (ev->type() == QEvent::KeyPress)
+ {
+ QKeyEvent *k = (QKeyEvent *)ev;
+ bool up = false;
+ bool down = false;
+ if (k->key() == 16777238) up = true; //page up
+ else if (k->key() == 16777239) down = true; //page down
+
+ if (down && m_tracking && !m_logging)
+ {
+ StartNextMeasurement();
+ }
+ }
+
+ return false;
+}
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.h b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.h
new file mode 100644
index 0000000000..85fe08d498
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.h
@@ -0,0 +1,133 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date$
+Version: $Revision$
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+#ifndef QmitkIGTTrackingSemiAutomaticMeasurementView_h
+#define QmitkIGTTrackingSemiAutomaticMeasurementView_h
+
+#include
+
+#include
+
+//QT
+#include
+
+//MITK
+#include
+#include
+#include
+#include
+#include
+
+#include "ui_QmitkIGTTrackingSemiAutomaticMeasurementViewControls.h"
+
+/*!
+ \brief QmitkIGTTrackingSemiAutomaticMeasurementView
+
+ \warning This application module is not yet documented. Use "svn blame/praise/annotate" and ask the author to provide basic documentation.
+
+ \sa QmitkFunctionality
+ \ingroup Functionalities
+ */
+class QmitkIGTTrackingSemiAutomaticMeasurementView : public QmitkFunctionality
+{
+ // this is needed for all Qt objects that should have a Qt meta-object
+ // (everything that derives from QObject and wants to have signal/slots)
+ Q_OBJECT
+
+public:
+
+ static const std::string VIEW_ID;
+
+ QmitkIGTTrackingSemiAutomaticMeasurementView();
+ virtual ~QmitkIGTTrackingSemiAutomaticMeasurementView();
+
+ virtual void CreateQtPartControl(QWidget *parent);
+
+ virtual void StdMultiWidgetAvailable(QmitkStdMultiWidget &stdMultiWidget);
+ virtual void StdMultiWidgetNotAvailable();
+
+ protected slots:
+
+ void OnLoadMeasurementStorage();
+ void OnLoadReferenceStorage();
+ void OnStartTracking();
+ void OnStopTracking();
+ void OnMeasurementLoadFile();
+ void OnSetReference();
+ void StartNextMeasurement();
+ void RepeatLastMeasurement();
+ void UpdateTimer();
+ void CreateResults();
+ void OnUseReferenceToggled(bool state);
+
+protected:
+
+ Ui::QmitkIGTTrackingSemiAutomaticMeasurementViewControls* m_Controls;
+
+ QmitkStdMultiWidget* m_MultiWidget;
+
+ //the tool storages
+ mitk::NavigationToolStorage::Pointer m_MeasurementStorage;
+ mitk::NavigationToolStorage::Pointer m_ReferenceStorage;
+
+ //members for the filter pipeline
+ mitk::TrackingDeviceSource::Pointer m_MeasurementTrackingDeviceSource;
+ mitk::NavigationDataObjectVisualizationFilter::Pointer m_MeasurementToolVisualizationFilter;
+ mitk::NavigationDataRecorderDeprecated::Pointer m_MeasurementLoggingFilterXML;
+ mitk::NavigationDataRecorderDeprecated::Pointer m_MeasurementLoggingFilterCSV;
+ mitk::TrackingDeviceSource::Pointer m_ReferenceTrackingDeviceSource;
+ mitk::NavigationDataRecorderDeprecated::Pointer m_ReferenceLoggingFilterXML;
+ mitk::NavigationDataRecorderDeprecated::Pointer m_ReferenceLoggingFilterCSV;
+
+ //members for file name list
+ std::vector m_FilenameVector;
+ int m_NextFile;
+
+ //help methods
+ mitk::NavigationToolStorage::Pointer ReadStorage(std::string file);
+ void MessageBox(std::string s);
+ void DisableAllButtons();
+ void EnableAllButtons();
+ void FinishMeasurement();
+ void StartLoggingAdditionalCSVFile(std::string filePostfix);
+ void LogAdditionalCSVFile();
+ void StopLoggingAdditionalCSVFile();
+
+ //timer
+ QTimer* m_Timer;
+
+ //memebers for reference checking
+ std::vector m_ReferenceStartPositions;
+ bool m_referenceValid;
+
+ //logging members
+ int m_loggedFrames;
+ bool m_logging;
+ std::fstream m_logFileCSV;
+
+ //event filter for key presses
+ bool eventFilter(QObject *obj, QEvent *ev);
+
+ //results members
+ mitk::PointSet::Pointer m_MeanPoints;
+ std::vector m_RMSValues;
+ mitk::NavigationDataEvaluationFilter::Pointer m_EvaluationFilter;
+
+ bool m_tracking;
+};
+
+#endif // _QMITKIGTTRACKINGSEMIAUTOMATICMEASUREMENTVIEW_H_INCLUDED
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementViewControls.ui b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementViewControls.ui
new file mode 100644
index 0000000000..50168cb767
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementViewControls.ui
@@ -0,0 +1,555 @@
+
+
+ QmitkIGTTrackingSemiAutomaticMeasurementViewControls
+
+
+
+ 0
+ 0
+ 419
+ 931
+
+
+
+
+ 0
+ 0
+
+
+
+ QmitkTemplate
+
+
+ -
+
+
+ 0
+
+
+
+ Tracking Initialization
+
+
+
-
+
+
+ Measurement Tracking System
+
+
+
-
+
+
+ -
+
+
-
+
+
+ Tool Storage: <none>
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Load Tool Storage
+
+
+
+
+
+
+
+
+ -
+
+
+ Reference Trackingsystem
+
+
+
-
+
+
+ -
+
+
-
+
+
+ Tool Storage: <none>
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Load Tool Storage
+
+
+
+
+
+
+
+
+ -
+
+
+ Start Tracking
+
+
+
+ -
+
+
+ Stop Tracking
+
+
+
+ -
+
+
+ Qt::Vertical
+
+
+
+ 20
+ 40
+
+
+
+
+
+
+
+
+ Measurement
+
+
+ -
+
+
+ Measurement List:
+
+
+
+ -
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ Load List
+
+
+
+
+
+ -
+
+
-
+
+
+ Output Path:
+
+
+
+ -
+
+
+ C:/temp/
+
+
+
+
+
+ -
+
+
+ Last Measurement: <none>
+
+
+
+ -
+
+
+ Next Measurement: <none>
+
+
+
+ -
+
+
+ Collected Samples: <none>
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 200
+ 0
+
+
+
+ Start Next Measurement
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 200
+ 0
+
+
+
+ Repeat Last Measurement
+
+
+
+
+
+ -
+
+
-
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+
+ 200
+ 0
+
+
+
+ Create Results
+
+
+
+
+
+ -
+
+
+ Qt::Vertical
+
+
+
+ 20
+ 281
+
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-weight:600; text-decoration: underline;">Measurement Sensors:</span></p></body></html>
+
+
+
+ -
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:8pt; font-weight:600; text-decoration: underline;">Reference Sensors:</span></p></body></html>
+
+
+
+ -
+
+
+ Set Reference
+
+
+
+ -
+
+
+ Reference Postion: <none>
+
+
+
+ -
+
+
+ Current Positions: <none>
+
+
+
+ -
+
+
+ <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0//EN" "http://www.w3.org/TR/REC-html40/strict.dtd">
+<html><head><meta name="qrichtext" content="1" /><style type="text/css">
+p, li { white-space: pre-wrap; }
+</style></head><body style=" font-family:'MS Shell Dlg 2'; font-size:8.25pt; font-weight:400; font-style:normal;">
+<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px;"><span style=" font-size:28pt; font-weight:600;">OK</span></p></body></html>
+
+
+ Qt::AlignCenter
+
+
+
+
+
+
+
+ Settings
+
+
+ -
+
+
-
+
+
+ Sampling Rate (Times Per Second):
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ 999.000000000000000
+
+
+ 15.000000000000000
+
+
+
+
+
+ -
+
+
-
+
+
+ Samples Per Measurement:
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ 1000
+
+
+ 150
+
+
+
+
+
+ -
+
+
-
+
+
+ Threshold For Reference Tools:
+
+
+
+ -
+
+
+ Qt::Horizontal
+
+
+
+ 40
+ 20
+
+
+
+
+ -
+
+
+ 1.000000000000000
+
+
+
+
+
+ -
+
+
+ Use Reference Tracking System
+
+
+ true
+
+
+
+ -
+
+
+ Qt::Vertical
+
+
+
+ 20
+ 40
+
+
+
+
+
+
+
+
+
+
+
+
+
+ QmitkToolTrackingStatusWidget
+ QWidget
+ QmitkToolTrackingStatusWidget.h
+ 1
+
+
+ QmitkTrackingDeviceConfigurationWidget
+ QWidget
+ QmitkTrackingDeviceConfigurationWidget.h
+ 1
+
+
+
+
+
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkNavigationDataCSVSequentialPlayer.cpp b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkNavigationDataCSVSequentialPlayer.cpp
new file mode 100644
index 0000000000..66f36f2fda
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkNavigationDataCSVSequentialPlayer.cpp
@@ -0,0 +1,214 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date: 2009-02-10 18:08:54 +0100 (Di, 10 Feb 2009) $
+Version: $Revision: 16228 $
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+#include "mitkNavigationDataCSVSequentialPlayer.h"
+#include
+#include
+#include
+#include
+
+
+mitk::NavigationDataCSVSequentialPlayer::NavigationDataCSVSequentialPlayer()
+ : mitk::NavigationDataPlayerBase()
+{
+m_NavigationDatas = std::vector();
+m_CurrentPos = 0;
+m_Filetype = mitk::NavigationDataCSVSequentialPlayer::NavigationDataCSV;
+}
+
+
+mitk::NavigationDataCSVSequentialPlayer::~NavigationDataCSVSequentialPlayer()
+{
+
+}
+
+bool mitk::NavigationDataCSVSequentialPlayer::IsAtEnd()
+{
+ if (m_CurrentPos >= m_NavigationDatas.size()) return true;
+ else return false;
+}
+
+void mitk::NavigationDataCSVSequentialPlayer::
+ SetFileName(const std::string& fileName)
+{
+ this->SetNumberOfOutputs(1);
+ FillOutputEmpty(0);
+
+ MITK_INFO << "Reading file: " << fileName;
+ m_NavigationDatas = GetNavigationDatasFromFile(fileName);
+
+ this->Modified();
+}
+
+void mitk::NavigationDataCSVSequentialPlayer::FillOutputEmpty(int number)
+ {
+ this->SetNthOutput(number,GetEmptyNavigationData());
+ }
+
+mitk::NavigationData::Pointer mitk::NavigationDataCSVSequentialPlayer::GetEmptyNavigationData()
+ {
+ mitk::NavigationData::Pointer emptyNd = mitk::NavigationData::New();
+ mitk::NavigationData::PositionType position;
+ mitk::NavigationData::OrientationType orientation(0.0,0.0,0.0,0.0);
+ position.Fill(0.0);
+
+ emptyNd->SetPosition(position);
+ emptyNd->SetOrientation(orientation);
+ emptyNd->SetDataValid(false);
+ return emptyNd;
+ }
+
+void mitk::NavigationDataCSVSequentialPlayer::GenerateData()
+{
+
+ for (unsigned int index = 0; index < this->GetNumberOfOutputs(); index++)
+ {
+
+ mitk::NavigationData* output = this->GetOutput(index);
+
+ if (m_CurrentPos > m_NavigationDatas.size())
+ {
+ FillOutputEmpty(index);
+ return;
+ }
+
+ output->Graft(this->m_NavigationDatas.at(m_CurrentPos));
+ m_CurrentPos++;
+
+ }
+}
+
+void mitk::NavigationDataCSVSequentialPlayer::UpdateOutputInformation()
+{
+ this->Modified(); // make sure that we need to be updated
+ Superclass::UpdateOutputInformation();
+}
+
+std::vector mitk::NavigationDataCSVSequentialPlayer::GetNavigationDatasFromFile(std::string filename)
+{
+std::vector returnValue = std::vector();
+std::vector fileContentLineByLine = GetFileContentLineByLine(filename);
+for(int i=1; (i mitk::NavigationDataCSVSequentialPlayer::GetFileContentLineByLine(std::string filename)
+{
+std::vector readData = std::vector();
+
+//save old locale
+char * oldLocale;
+oldLocale = setlocale( LC_ALL, 0 );
+
+//define own locale
+std::locale C("C");
+setlocale( LC_ALL, "C" );
+
+//read file
+std::ifstream file;
+file.open(filename.c_str(), std::ios::in);
+if (file.good())
+ {
+ //read out file
+ file.seekg(0L, std::ios::beg); // move to begin of file
+ while (! file.eof())
+ {
+ std::string buffer;
+ std::getline(file,buffer); // read out file line by line
+ if (buffer.size() > 0) readData.push_back(buffer);
+
+ }
+ }
+
+file.close();
+
+//switch back to old locale
+setlocale( LC_ALL, oldLocale );
+
+return readData;
+}
+
+mitk::NavigationData::Pointer mitk::NavigationDataCSVSequentialPlayer::GetNavigationDataOutOfOneLine(std::string line)
+{
+ mitk::NavigationData::Pointer returnValue = mitk::NavigationData::New();
+
+ QString myLine = QString(line.c_str());
+
+ QStringList myLineList = myLine.split(';');
+
+ mitk::Point3D position;
+ mitk::Quaternion orientation;
+ bool valid = false;
+ double time;
+
+ if (m_Filetype = mitk::NavigationDataCSVSequentialPlayer::NavigationDataCSV)
+ {
+
+ if (myLineList.size() < 10)
+ {
+ MITK_ERROR << "Error: cannot read line: only found " << myLineList.size() << " fields. Last field: " << myLineList.at(myLineList.size()-1).toStdString() ;
+ returnValue = GetEmptyNavigationData();
+ return returnValue;
+ }
+
+ time = myLineList.at(1).toDouble();
+
+ if (myLineList.at(2).toStdString() == "1") valid = true;
+
+ position[0] = myLineList.at(3).toDouble();
+ position[1] = myLineList.at(4).toDouble();
+ position[2] = myLineList.at(5).toDouble();
+
+ orientation[0] = myLineList.at(6).toDouble();
+ orientation[1] = myLineList.at(7).toDouble();
+ orientation[2] = myLineList.at(8).toDouble();
+ orientation[3] = myLineList.at(9).toDouble();
+ }
+ else
+ {
+ if (myLineList.size() < 10)
+ {
+ MITK_ERROR << "Error: cannot read line: only found " << myLineList.size() << " fields. Last field: " << myLineList.at(myLineList.size()-1).toStdString() ;
+ returnValue = GetEmptyNavigationData();
+ return returnValue;
+ }
+
+ time = myLineList.at(1).toDouble();
+
+ //if (myLineList.at(2).toStdString() == "true")
+ //valid-flag wurde hier nicht gespeichert
+ valid = true;
+
+ position[0] = myLineList.at(3).toDouble();
+ position[1] = myLineList.at(4).toDouble();
+ position[2] = myLineList.at(5).toDouble();
+
+ orientation[0] = myLineList.at(6).toDouble();
+ orientation[1] = myLineList.at(7).toDouble();
+ orientation[2] = myLineList.at(8).toDouble();
+ orientation[3] = myLineList.at(9).toDouble();
+ }
+
+ //returnValue->SetTimeStamp(time); //DOES NOT WORK ANY MORE... CANNOT SET TIME TO itk::timestamp CLASS
+ returnValue->SetDataValid(valid);
+ returnValue->SetPosition(position);
+ returnValue->SetOrientation(orientation);
+
+ return returnValue;
+}
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkNavigationDataCSVSequentialPlayer.h b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkNavigationDataCSVSequentialPlayer.h
new file mode 100644
index 0000000000..50bb773d18
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkNavigationDataCSVSequentialPlayer.h
@@ -0,0 +1,98 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date: 2009-02-10 18:08:54 +0100 (Di, 10 Feb 2009) $
+Version: $Revision: 16228 $
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+
+#ifndef MITKNavigationDataCSVSequentialPlayer_H_HEADER_INCLUDED_
+#define MITKNavigationDataCSVSequentialPlayer_H_HEADER_INCLUDED_
+
+#include
+#include "tinyxml.h"
+
+
+namespace mitk
+{
+
+ /**Documentation
+ * \brief This class is a NavigationDataPlayer which can play CSV formatted
+ * files in sequential order, which means it doesn't
+ * care about timestamps and just
+ * outputs the navigationdatas in their sequential order
+ *
+ * \ingroup IGT
+ */
+ class NavigationDataCSVSequentialPlayer
+ : public NavigationDataPlayerBase
+ {
+ public:
+
+ enum Filetype
+ {
+ NavigationDataCSV,
+ ManualLoggingCSV
+ };
+
+ mitkClassMacro(NavigationDataCSVSequentialPlayer, NavigationDataPlayerBase);
+ itkNewMacro(Self);
+
+ /**
+ * \brief sets the file name and path (if XMLString is set, this is neglected)
+ */
+ void SetFileName(const std::string& _FileName);
+
+ /**
+ * \brief returns the file name and path
+ */
+ itkGetStringMacro(FileName);
+
+ itkSetMacro(Filetype,Filetype);
+
+ bool IsAtEnd();
+
+ /**
+ * \brief Used for pipeline update just to tell the pipeline
+ * that we always have to update
+ */
+ virtual void UpdateOutputInformation();
+
+ protected:
+ NavigationDataCSVSequentialPlayer();
+ virtual ~NavigationDataCSVSequentialPlayer();
+
+ ///
+ /// do the work here
+ ///
+ virtual void GenerateData();
+
+ std::string m_FileName;
+
+ int m_CurrentPos;
+ Filetype m_Filetype;
+
+ //member for the navigation datas which were read (only one output is supported at the moment)
+ std::vector m_NavigationDatas;
+
+ std::vector GetNavigationDatasFromFile(std::string filename);
+ std::vector GetFileContentLineByLine(std::string filename);
+ mitk::NavigationData::Pointer GetNavigationDataOutOfOneLine(std::string line);
+
+ void FillOutputEmpty(int number);
+ mitk::NavigationData::Pointer GetEmptyNavigationData();
+
+ };
+} // namespace mitk
+
+#endif /* MITKNavigationDataCSVSequentialPlayer_H_HEADER_INCLUDED_ */
\ No newline at end of file
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkPluginActivator.cpp b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkPluginActivator.cpp
new file mode 100644
index 0000000000..b1ca328364
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkPluginActivator.cpp
@@ -0,0 +1,38 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date$
+Version: $Revision$
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+#include "mitkPluginActivator.h"
+
+#include
+
+#include "QmitkIGTTrackingSemiAutomaticMeasurementView.h"
+
+namespace mitk {
+ void PluginActivator::start(ctkPluginContext* context)
+ {
+ BERRY_REGISTER_EXTENSION_CLASS(QmitkIGTTrackingSemiAutomaticMeasurementView, context)
+ }
+
+ void PluginActivator::stop(ctkPluginContext* context)
+ {
+ Q_UNUSED(context)
+ }
+}
+
+#if QT_VERSION < QT_VERSION_CHECK(5, 0, 0)
+Q_EXPORT_PLUGIN2(org_mitk_gui_qt_igttrackingsemiautomaticmeasurement, mitk::PluginActivator)
+#endif
\ No newline at end of file
diff --git a/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkPluginActivator.h b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkPluginActivator.h
new file mode 100644
index 0000000000..ecf559f1d4
--- /dev/null
+++ b/Plugins/org.mitk.gui.qt.igt.app.hummelprotocolmeasurements/src/internal/mitkPluginActivator.h
@@ -0,0 +1,44 @@
+/*=========================================================================
+
+Program: Medical Imaging & Interaction Toolkit
+Language: C++
+Date: $Date$
+Version: $Revision$
+
+Copyright (c) German Cancer Research Center, Division of Medical and
+Biological Informatics. All rights reserved.
+See MITKCopyright.txt or http://www.mitk.org/copyright.html for details.
+
+This software is distributed WITHOUT ANY WARRANTY; without even
+the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+
+#ifndef MITKPLUGINACTIVATOR_H
+#define MITKPLUGINACTIVATOR_H
+
+#include
+
+namespace mitk {
+
+class PluginActivator :
+ public QObject, public ctkPluginActivator
+{
+ Q_OBJECT
+ #if QT_VERSION >= QT_VERSION_CHECK(5, 0, 0)
+ Q_PLUGIN_METADATA(IID "org_mitk_gui_qt_igtapphummelprotocolmeasurements")
+ #endif
+ Q_INTERFACES(ctkPluginActivator)
+
+public:
+
+ void start(ctkPluginContext* context);
+ void stop(ctkPluginContext* context);
+
+}; // PluginActivator
+
+}
+
+#endif // MITKPLUGINACTIVATOR_H