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->GetStreamValid()) 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::NavigationDataRecorder::Pointer myRecorder = mitk::NavigationDataRecorder::New();
- myRecorder->SetFileName(outputFilename.c_str());
- mitk::NavigationData::Pointer input = mitk::NavigationData::New();
- if(m_Controls->m_ConvertCSV->isChecked()) myRecorder->SetOutputFormat(mitk::NavigationDataRecorder::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.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingDataEvaluationView.h b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingDataEvaluationView.h
deleted file mode 100644
index a928c578b3..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingDataEvaluationView.h
+++ /dev/null
@@ -1,125 +0,0 @@
-/*=========================================================================
-
-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.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingDataEvaluationViewControls.ui b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingDataEvaluationViewControls.ui
deleted file mode 100644
index 9095e6a5af..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingDataEvaluationViewControls.ui
+++ /dev/null
@@ -1,941 +0,0 @@
-
-
- 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.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.cpp b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.cpp
deleted file mode 100644
index 9a0c25196c..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.cpp
+++ /dev/null
@@ -1,667 +0,0 @@
-/*=========================================================================
-
-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.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.h b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.h
deleted file mode 100644
index 85fe08d498..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementView.h
+++ /dev/null
@@ -1,133 +0,0 @@
-/*=========================================================================
-
-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.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementViewControls.ui b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementViewControls.ui
deleted file mode 100644
index 50168cb767..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/QmitkIGTTrackingSemiAutomaticMeasurementViewControls.ui
+++ /dev/null
@@ -1,555 +0,0 @@
-
-
- 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.igttrackingsemiautomaticmeasurement/src/internal/mitkPluginActivator.cpp b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/mitkPluginActivator.cpp
deleted file mode 100644
index b1ca328364..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/mitkPluginActivator.cpp
+++ /dev/null
@@ -1,38 +0,0 @@
-/*=========================================================================
-
-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.igttrackingsemiautomaticmeasurement/src/internal/mitkPluginActivator.h b/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/mitkPluginActivator.h
deleted file mode 100644
index f1bb717002..0000000000
--- a/Plugins/org.mitk.gui.qt.igttrackingsemiautomaticmeasurement/src/internal/mitkPluginActivator.h
+++ /dev/null
@@ -1,44 +0,0 @@
-/*=========================================================================
-
-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_igttrackingsemiautomaticmeasurement")
- #endif
- Q_INTERFACES(ctkPluginActivator)
-
-public:
-
- void start(ctkPluginContext* context);
- void stop(ctkPluginContext* context);
-
-}; // PluginActivator
-
-}
-
-#endif // MITKPLUGINACTIVATOR_H