diff --git a/101915/Functional/Movement_Regressors.txt b/101915/Functional/Movement_Regressors.txt new file mode 100644 index 0000000..83a39e2 --- /dev/null +++ b/101915/Functional/Movement_Regressors.txt @@ -0,0 +1,1200 @@ +0.091614 0.11409 0.1303 -0.001428 0.00069 0 +0.054535 0.1271 0.12682 -0.000981 0 0 +0.13167 0.12992 0.15864 -0.001421 0.00053899 0.000555 +0.15006 0.13967 0.25098 -0.002781 0.000545 0.000396 +0.10721 0.12966 0.19533 -0.002462 0.001004 0 +0.06376 0.13955 0.19969 -0.001963 0 0.00033599 +0.097042 0.15591 0.22882 -0.002383 0 0.00046 +0.13926 0.15073 0.30378 -0.003654 0.00071901 0.00037299 +0.072679 0.15863 0.27099 -0.002782 0.00072799 0.000369 +0.080584 0.13767 0.24908 -0.00267 0 0.000271 +0.12384 0.15657 0.2446 -0.002141 0 0.00048799 +0.15055 0.15081 0.32101 -0.003472 0.001051 0.000453 +0.11027 0.15041 0.3056 -0.003389 0.001011 0.00037401 +0.047929 0.14768 0.29313 -0.002622 0 0.00040399 +0.091828 0.15714 0.24736 -0.002248 0 0.000511 +0.1136 0.15938 0.28629 -0.002687 0.000723 0.000513 +0.16402 0.16442 0.35507 -0.004174 0.001143 0 +0.12141 0.1615 0.33004 -0.003674 0.001139 0.000581 +0.053893 0.16213 0.30291 -0.002818 0 0.000568 +0.082924 0.1859 0.26223 -0.002726 0 0.000525 +0.11719 0.16255 0.32638 -0.003513 0.000775 0.000633 +0.14737 0.16587 0.35784 -0.004327 0.001071 0.00057701 +0.063318 0.16845 0.32408 -0.003287 0.001029 0.000538 +0.096242 0.17654 0.30581 -0.002529 0.00074499 0.000733 +0.12355 0.16332 0.32729 -0.003269 0.00071201 0.000761 +0.21345 0.16586 0.38015 -0.004208 0.001517 0.00086101 +0.17575 0.18495 0.36718 -0.004344 0.001228 0.000717 +0.068166 0.17156 0.3434 -0.003183 0.001014 0.00068799 +0.14544 0.17634 0.32408 -0.00343 0 0.000768 +0.14878 0.18805 0.37512 -0.00408 0.001167 0.00075899 +0.15798 0.1963 0.41773 -0.004911 0.001331 0.00078199 +0.095567 0.18912 0.37616 -0.003987 0.001196 0.000737 +0.077707 0.17235 0.35845 -0.003755 0.000805 0.000785 +0.095204 0.16923 0.36628 -0.0036 0.001324 0.000748 +0.14918 0.18037 0.43744 -0.00457 0.001482 0.00088399 +0.17463 0.1955 0.47248 -0.005373 0.001508 0.00079301 +0.08933 0.19607 0.39971 -0.003819 0.001233 0.000714 +0.082703 0.17539 0.36423 -0.003424 0.001296 0.000774 +0.097687 0.18478 0.34232 -0.003179 0.001296 0.00082001 +0.12997 0.1968 0.42512 -0.005004 0.001525 0.00075499 +0.15028 0.21276 0.47847 -0.00542 0.001365 0.00081701 +0.10498 0.2028 0.4618 -0.005531 0.001541 0.00077301 +0.045392 0.19083 0.42478 -0.005123 0.001509 0.000689 +0.072569 0.18839 0.39499 -0.004193 0.001407 0.000697 +0.12907 0.20387 0.44484 -0.00539 0.001442 0.000849 +0.1617 0.21592 0.48245 -0.005889 0.001587 0.00082999 +0.074051 0.21498 0.42116 -0.00508 0.001471 0 +0.1398 0.19635 0.39072 -0.005024 0.001226 0.00081701 +0.11775 0.20407 0.43678 -0.005591 0.001341 0.000824 +0.2092 0.20134 0.48811 -0.006342 0.001722 0.000953 +0.12736 0.2219 0.44961 -0.005713 0.001345 0.000778 +0.094602 0.20984 0.43172 -0.00578 0.001225 0.00082701 +0.1194 0.21633 0.49387 -0.006014 0.001393 0.00086401 +0.1627 0.20368 0.55172 -0.00778 0.001816 0.000841 +0.1712 0.20212 0.52143 -0.006419 0.001566 0.000944 +0.078277 0.22044 0.46946 -0.005333 0.001866 0.000988 +0.086728 0.22256 0.46353 -0.00574 0.001552 0.00068501 +0.15055 0.2011 0.54158 -0.007656 0.001359 0.000778 +0.16207 0.20679 0.57853 -0.008404 0.001685 0.000753 +0.13825 0.21764 0.49063 -0.006376 0.001667 0.000828 +0.14266 0.20054 0.45365 -0.005306 0.001657 0.00096201 +0.18871 0.21555 0.47025 -0.005956 0.00187 0.000984 +0.16093 0.22814 0.56115 -0.007532 0.001883 0.000673 +0.1254 0.22082 0.57542 -0.008411 0.001683 0.00065801 +0.074265 0.21555 0.52567 -0.006901 0.001449 0.0007 +0.081201 0.22844 0.48901 -0.006131 0.001367 0.00064699 +0.13295 0.26125 0.55414 -0.006942 0.001286 0.00060401 +0.15486 0.23416 0.6118 -0.008783 0.001399 0.00065399 +0.09794 0.23355 0.5912 -0.008483 0.001445 0.00051499 +0.058069 0.21097 0.54698 -0.007888 0.000859 0 +0.095727 0.24283 0.56103 -0.006728 0.001201 0.00061401 +0.20933 0.23677 0.60573 -0.008588 0.001373 0.00062699 +0.18802 0.24044 0.60732 -0.008862 0.001395 0.000554 +0.077733 0.22468 0.57002 -0.007832 0.001267 0.000557 +0.085923 0.23618 0.52841 -0.006182 0.001136 0.00069201 +0.13221 0.21808 0.56429 -0.007752 0.000991 0.00068501 +0.1712 0.2369 0.60643 -0.00861 0.001566 0.00067399 +0.13239 0.22914 0.59938 -0.008748 0.00116 0.000679 +0.079153 0.22755 0.5665 -0.007517 0.001287 0.000683 +0.10104 0.24231 0.5135 -0.006167 0.001195 0.000751 +0.12299 0.24517 0.51246 -0.006022 0.001247 0.00086401 +0.15927 0.24405 0.60844 -0.00799 0.001439 0.000798 +0.20654 0.25894 0.65648 -0.008742 0.00162 0.000632 +0.15406 0.24888 0.6381 -0.0085 0.001392 0.00061401 +0.1381 0.23747 0.63538 -0.008433 0.001292 0.00063101 +0.16032 0.24439 0.65207 -0.008517 0.001437 0.000726 +0.20925 0.25968 0.65108 -0.008947 0.001526 0.000602 +0.16557 0.23922 0.66458 -0.009156 0.001608 0.000652 +0.089068 0.24963 0.60998 -0.008225 0.000805 0.000568 +0.11959 0.22882 0.59726 -0.00819 0.001336 0.000574 +0.15181 0.2423 0.64655 -0.009146 0.001582 0.000575 +0.17257 0.23854 0.67997 -0.009878 0.001579 0.00047499 +0.1122 0.23891 0.65233 -0.009311 0.001439 0.000396 +0.084861 0.23446 0.63303 -0.008327 0.001264 0.000525 +0.13521 0.23674 0.65251 -0.008646 0.001324 0.000564 +0.17654 0.25247 0.69142 -0.009599 0.001791 0.000548 +0.14106 0.23651 0.66435 -0.008574 0.001701 0.00050199 +0.086551 0.22374 0.64273 -0.008351 0.001328 0.000464 +0.096584 0.21272 0.63662 -0.008452 0.001438 0.00045799 +0.1792 0.24191 0.70311 -0.008832 0.001672 0.00049201 +0.17033 0.24704 0.72262 -0.009248 0.001619 0.00046799 +0.1228 0.23703 0.68287 -0.008578 0.001315 0.000325 +0.081464 0.24216 0.64306 -0.007982 0.001401 0 +0.13905 0.24369 0.6802 -0.008335 0.001402 0.000454 +0.18288 0.2509 0.723 -0.00945 0.001933 0 +0.1631 0.243 0.6955 -0.009085 0.001585 0.00040399 +0.12761 0.23865 0.62517 -0.008092 0.00093201 0.000447 +0.13008 0.23348 0.64586 -0.008352 0.001491 0.000508 +0.17468 0.24947 0.67777 -0.009217 0.001857 0.00039699 +0.14728 0.25432 0.67618 -0.00948 0.001838 0.000477 +0.087167 0.23362 0.64368 -0.008255 0.001328 0.00051499 +0.15201 0.23239 0.62774 -0.008196 0.00158 0.00058599 +0.18403 0.25731 0.68685 -0.00913 0.002048 0.000601 +0.16578 0.25753 0.69286 -0.00952 0.001851 0.000565 +0.14422 0.24495 0.64701 -0.008826 0.001485 0.00053 +0.15246 0.23981 0.62804 -0.008109 0.001617 0.000574 +0.19499 0.25522 0.68429 -0.008311 0.001656 0.000564 +0.21122 0.26604 0.7545 -0.009896 0.001956 0.000555 +0.15939 0.2695 0.71082 -0.009055 0.001893 0.00041699 +0.11981 0.27026 0.69123 -0.008209 0.001354 0.000491 +0.14835 0.25842 0.71966 -0.009101 0.001846 0.00051901 +0.19854 0.26508 0.75441 -0.009648 0.001952 0.00041699 +0.15913 0.25726 0.73036 -0.00942 0.001875 0.00039401 +0.095724 0.24614 0.69852 -0.008743 0.001783 0.000449 +0.16204 0.24313 0.67626 -0.00842 0.001848 0.000511 +0.17136 0.26488 0.71196 -0.008561 0.002029 0.00046799 +0.18777 0.25475 0.74526 -0.009566 0.002255 0.000437 +0.17183 0.24534 0.71451 -0.00943 0.001771 0.000399 +0.13627 0.24896 0.67148 -0.008546 0.001733 0.000439 +0.13821 0.24186 0.70256 -0.008532 0.002038 0.00044501 +0.20181 0.25482 0.75065 -0.009555 0.001943 0.000473 +0.14144 0.24765 0.73567 -0.009298 0.002131 0.00036999 +0.098549 0.24272 0.6917 -0.008225 0.001754 0.000399 +0.1218 0.25502 0.67293 -0.007385 0.00172 0.00062101 +0.18472 0.25565 0.73655 -0.009092 0.002295 0.00041 +0.17561 0.26411 0.74617 -0.009629 0.00195 0.000369 +0.10912 0.24792 0.72218 -0.009048 0.001787 0.000419 +0.13136 0.24418 0.67641 -0.00815 0.001829 0.000348 +0.13885 0.25162 0.6836 -0.008278 0.001708 0.000415 +0.22153 0.2571 0.75695 -0.008952 0.002463 0.000432 +0.19218 0.25226 0.76469 -0.009349 0.002462 0 +0.11693 0.25733 0.69961 -0.008118 0.001993 0.000403 +0.10115 0.24114 0.69159 -0.008037 0.001805 0.00048901 +0.13472 0.26139 0.71498 -0.008474 0.001943 0.000595 +0.1818 0.26669 0.75809 -0.009135 0.002594 0.00047799 +0.13897 0.26149 0.7291 -0.008624 0.002499 0.00047901 +0.10917 0.25941 0.70721 -0.007718 0.001663 0.00052899 +0.14377 0.2594 0.69565 -0.007689 0.001809 0.00049201 +0.20034 0.26333 0.76026 -0.008847 0.002527 0.000402 +0.22715 0.26834 0.77945 -0.00928 0.002499 0.0005 +0.14631 0.25706 0.74162 -0.008462 0.001995 0.000483 +0.13506 0.23759 0.7139 -0.008239 0.002001 0.00048499 +0.16477 0.25719 0.75219 -0.008349 0.002406 0.00061299 +0.2039 0.26474 0.78819 -0.009389 0.002684 0.000565 +0.13995 0.26348 0.75891 -0.00898 0.002562 0.000446 +0.098415 0.2561 0.71558 -0.007833 0.002067 0.00046501 +0.1441 0.23449 0.71517 -0.008141 0.002151 0.00056701 +0.19801 0.26686 0.77587 -0.008843 0.00266 0.000521 +0.23006 0.27078 0.80159 -0.009585 0.002781 0.000518 +0.17102 0.26605 0.77795 -0.009 0.002661 0.000513 +0.11662 0.28168 0.69385 -0.007186 0.002133 0.000564 +0.16034 0.26687 0.70711 -0.008104 0.002536 0.00058599 +0.22162 0.26118 0.77656 -0.009001 0.002929 0.00064801 +0.22603 0.26213 0.80538 -0.010253 0.002986 0.000633 +0.1361 0.25566 0.76624 -0.008785 0.002748 0.000524 +0.17893 0.25873 0.73179 -0.008096 0.00286 0.00059001 +0.18082 0.27196 0.73233 -0.007643 0.002588 0.00071799 +0.23601 0.28826 0.78846 -0.009087 0.002826 0.000625 +0.19112 0.28078 0.78599 -0.009048 0.003017 0.000592 +0.15425 0.2733 0.75792 -0.00888 0.002413 0.000591 +0.16556 0.2621 0.7526 -0.00883 0.002834 0.000628 +0.22845 0.2691 0.80369 -0.009788 0.002947 0.000598 +0.22005 0.28323 0.81312 -0.0096 0.002995 0.000463 +0.14434 0.2663 0.80411 -0.009679 0.002727 0.00051199 +0.1213 0.24854 0.77394 -0.008573 0.002599 0.00052899 +0.17963 0.26051 0.79989 -0.008866 0.002863 0.00061699 +0.23388 0.26314 0.84422 -0.010382 0.002805 0.00053601 +0.20283 0.25732 0.84515 -0.010363 0.002707 0.00053199 +0.12679 0.26084 0.81875 -0.009163 0.002483 0.00051601 +0.16319 0.25298 0.79126 -0.008199 0.0027 0.000572 +0.22953 0.26893 0.84648 -0.008991 0.002867 0.00057001 +0.26737 0.27684 0.87524 -0.010318 0.003127 0.00050499 +0.17408 0.27013 0.85236 -0.009529 0.002811 0.000361 +0.14872 0.26044 0.81075 -0.00815 0.002719 0.00045201 +0.22165 0.27455 0.84675 -0.009293 0.002925 0.000606 +0.26845 0.28341 0.89663 -0.009757 0.003122 0.00050499 +0.20583 0.27931 0.87571 -0.010009 0.002941 0.00042 +0.14923 0.26982 0.82161 -0.008254 0.002707 0.000486 +0.16189 0.26766 0.81598 -0.0083 0.002851 0.00049201 +0.26346 0.28128 0.8854 -0.009405 0.003119 0.00052899 +0.23876 0.28643 0.89174 -0.009815 0.003042 0.00050499 +0.16365 0.27381 0.86557 -0.00946 0.002884 0.00045 +0.17971 0.27201 0.8232 -0.008118 0.00294 0.000501 +0.2021 0.27147 0.82547 -0.008105 0.003133 0.00060999 +0.27977 0.2813 0.899 -0.009296 0.003413 0.00053601 +0.25764 0.29428 0.92243 -0.010497 0.003288 0.000433 +0.18447 0.28186 0.87386 -0.009441 0.003029 0.000437 +0.17737 0.28364 0.82832 -0.007767 0.003042 0.00061999 +0.21825 0.28892 0.86358 -0.008595 0.003246 0.00050601 +0.27113 0.29901 0.92267 -0.009733 0.003451 0.000437 +0.23034 0.29916 0.91646 -0.009692 0.003334 0.00040101 +0.16066 0.28741 0.84832 -0.008859 0.003037 0.00041699 +0.18191 0.28133 0.79907 -0.00827 0.00265 0.00045 +0.25271 0.29424 0.87382 -0.009047 0.003415 0.00049201 +0.29179 0.30031 0.91328 -0.009867 0.003494 0.00051 +0.20203 0.29317 0.88227 -0.009179 0.003245 0.00044501 +0.16984 0.28631 0.80396 -0.007986 0.003287 0.0005 +0.20733 0.30129 0.79292 -0.007711 0.003219 0.00058001 +0.26219 0.3037 0.88582 -0.009197 0.003428 0.000538 +0.28781 0.31743 0.9247 -0.010229 0.003562 0.00046099 +0.18303 0.31111 0.88304 -0.008956 0.003258 0.000359 +0.20016 0.29654 0.82244 -0.00805 0.003275 0.000611 +0.22339 0.31548 0.84537 -0.008412 0.003299 0.00050199 +0.2752 0.31445 0.94468 -0.010014 0.003492 0.00030899 +0.2753 0.31657 0.95623 -0.010624 0.00354 0.000305 +0.19227 0.30105 0.92289 -0.009733 0.003274 0.000355 +0.20577 0.29333 0.87885 -0.008956 0.003312 0.00047 +0.24322 0.31104 0.8569 -0.00796 0.003938 0.00058701 +0.28562 0.31324 0.95301 -0.009791 0.003646 0.000474 +0.25732 0.30995 0.92779 -0.009292 0.003745 0.00033701 +0.20019 0.29947 0.87782 -0.008491 0.003879 0 +0.21427 0.29575 0.80254 -0.007759 0.003296 0 +0.2385 0.2883 0.85219 -0.007956 0.003698 0 +0.28706 0.30446 0.9582 -0.009069 0.004325 0.000114 +0.25659 0.29997 0.8169 -0.008633 0.003924 -0.00012999 +0.11048 0.25144 0.67107 -0.007708 0.002634 0 +0.18529 0.27166 0.75587 -0.008646 0.003005 0 +0.27657 0.29431 0.89942 -0.01048 0.003692 0 +0.25496 0.30935 0.9365 -0.010962 0.00347 3.002e-06 +0.17604 0.29619 0.89562 -0.010176 0.002866 0.00011299 +0.12978 0.28427 0.85557 -0.009611 0.002835 0 +0.19519 0.29237 0.88209 -0.00991 0.002997 0 +0.26546 0.30635 0.92276 -0.010401 0.003561 0 +0.22394 0.31069 0.93083 -0.01056 0.003143 0 +0.1499 0.29522 0.8876 -0.00972 0.002735 0 +0.17883 0.29054 0.85096 -0.009242 0.002887 0 +0.26599 0.3004 0.92769 -0.010947 0.003278 0.000158 +0.28812 0.31778 0.95524 -0.010633 0.003436 0.000105 +0.1797 0.3135 0.87263 -0.008872 0.002835 9.8995e-05 +0.15482 0.30505 0.82647 -0.008083 0.002765 0.00017 +0.18487 0.3039 0.83325 -0.008028 0.002838 0.00024499 +0.2695 0.31932 0.9265 -0.00951 0.003135 0.00016801 +0.27609 0.32712 0.98682 -0.010711 0.003298 0.000243 +0.18179 0.32105 0.95288 -0.009862 0.003025 0.00021201 +0.16854 0.31815 0.90453 -0.009066 0.00274 0.00026901 +0.20767 0.30927 0.93815 -0.009843 0.002945 0.00020399 +0.28959 0.32636 1.0066 -0.010971 0.00331 0.000196 +0.25508 0.33241 1.0031 -0.010661 0.003286 0.00011701 +0.18175 0.3225 0.97044 -0.010153 0.002694 0.00020801 +0.18438 0.31545 0.89921 -0.008717 0.002716 0.000352 +0.23377 0.32664 0.87015 -0.008139 0.002724 0.000304 +0.29679 0.33204 0.98549 -0.010552 0.003151 0.000457 +0.29545 0.34602 1.0361 -0.011534 0.003328 0.000246 +0.20777 0.34569 1.0067 -0.01079 0.003138 0.000139 +0.19595 0.33289 0.95764 -0.009493 0.003183 0.00019099 +0.20396 0.34033 0.95136 -0.009344 0.003011 0.00025 +0.25408 0.34379 1.0136 -0.010239 0.003102 0.000268 +0.27686 0.3573 1.0463 -0.010777 0.003191 0.00014699 +0.21826 0.34578 1.001 -0.010139 0.003116 0.000121 +0.20034 0.32718 0.99391 -0.010247 0.002923 0.000206 +0.25243 0.33811 1.0069 -0.010191 0.003046 0.000298 +0.29158 0.34876 1.0413 -0.010769 0.003209 0.00020399 +0.24532 0.33728 1.0139 -0.010476 0.003023 0.000163 +0.18997 0.32158 0.96616 -0.010387 0.002922 0.000162 +0.23778 0.32587 0.93183 -0.009286 0.003052 0.000223 +0.28166 0.34594 0.99523 -0.009771 0.003078 0.00023901 +0.29388 0.35751 1.0385 -0.010896 0.003247 0.000173 +0.21203 0.34518 1.001 -0.010278 0.003005 0.000111 +0.20743 0.33228 0.96649 -0.009711 0.002856 0.000236 +0.24903 0.34364 0.98683 -0.009761 0.002697 0.00022501 +0.31388 0.36189 1.0523 -0.010985 0.003021 0.00019399 +0.25659 0.35739 1.0307 -0.010159 0.002895 0.00011001 +0.18366 0.34374 0.95485 -0.009207 0.002918 0 +0.11416 0.34099 0.75051 -0.007922 0.001298 -0.000156 +0.089223 0.32404 0.61182 -0.007573 0 0.00032599 +0.23422 0.32569 0.97094 -0.01114 0.003196 -4.9916e-06 +0.30722 0.32368 1.0466 -0.010927 0.003987 -0.000155 +0.13653 0.24929 0.71755 -0.00875 0.001641 0 +0.0874 0.22166 0.6253 -0.008244 0.00094799 -0.00059701 +0.11133 0.2341 0.52561 -0.007369 0 -0.000846 +0.12435 0.25829 0.61912 -0.008153 0 -0.001091 +0.22692 0.28276 0.82966 -0.01119 0.001021 -0.000734 +0.22302 0.30919 0.91646 -0.012055 0.00093099 -0.00045799 +0.14135 0.31384 0.91696 -0.011725 0 -0.00032599 +0.074454 0.29149 0.83396 -0.010366 0 -0.00030301 +0.08364 0.2912 0.78801 -0.009336 0 -0.00021099 +0.17471 0.326 0.85283 -0.009234 0.000229 -0.000139 +0.23515 0.34885 0.91518 -0.010915 0.000687 0 +0.16473 0.33449 0.93871 -0.010849 0.000477 0 +0.15412 0.31474 0.88966 -0.010211 0.00023799 0 +0.2072 0.32286 0.89153 -0.010467 0.000497 0 +0.2782 0.34465 0.96108 -0.011068 0.000741 0 +0.26432 0.35262 0.98437 -0.011224 0.000792 0 +0.1596 0.35143 0.96365 -0.010537 0 0 +0.15029 0.33002 0.91688 -0.010132 0.00043 0 +0.19111 0.34314 0.9855 -0.010817 0.000626 0 +0.25032 0.35546 1.0404 -0.011656 0.000859 0 +0.18835 0.35278 1.0059 -0.010954 0 -0.00017399 +0.15191 0.34861 0.94804 -0.009661 0 -7.5992e-05 +0.16568 0.32139 0.83482 -0.008926 0 0 +0.1862 0.33737 0.84286 -0.008653 0 0 +0.32342 0.36588 1.0281 -0.010336 0.001635 0 +0.28865 0.3479 0.98981 -0.010207 0.001697 -0.0005 +0.12121 0.30863 0.74104 -0.008079 0.00070099 -0.000605 +0.15471 0.28341 0.63739 -0.007898 0.000331 -0.000635 +0.18338 0.29258 0.8212 -0.009316 0.001572 -0.00086101 +0.31111 0.29376 0.99635 -0.011246 0.002479 -0.000845 +0.27244 0.28074 0.91015 -0.010137 0.002166 -0.000655 +0.25544 0.29311 0.86708 -0.010015 0.002025 0 +0.29762 0.32026 0.9914 -0.012487 0.002024 0 +0.25593 0.33509 1.0318 -0.012637 0.001426 -0.00057299 +0.17703 0.33177 1.009 -0.012826 0.001308 0 +0.1049 0.30974 0.94641 -0.011093 0.00103 0 +0.11247 0.29739 0.86719 -0.010312 0.000892 0 +0.21371 0.32325 0.95069 -0.011338 0.001281 0 +0.27593 0.33004 1.0112 -0.012498 0.001317 0 +0.24494 0.32516 1.0096 -0.012671 0.001093 0 +0.17353 0.31188 0.94741 -0.011534 0.00096201 0 +0.19143 0.31634 0.88697 -0.010081 0.00095201 0 +0.25289 0.32691 0.94397 -0.011105 0.001209 0 +0.3036 0.3445 1.0256 -0.012324 0.001338 -0.000345 +0.22426 0.34753 1.0182 -0.011935 0.001195 -0.00034001 +0.14625 0.32756 0.96609 -0.011035 0.001127 -0.000294 +0.14727 0.31816 0.91479 -0.010742 0.00085699 -7.4997e-05 +0.23199 0.34623 1.0018 -0.011918 0.001246 -0.000247 +0.2403 0.34186 1.036 -0.012073 0.001175 -0.000325 +0.21161 0.34731 1.0238 -0.011949 0.001227 0 +0.14525 0.33694 0.9702 -0.011202 0.001 -0.000192 +0.17236 0.3386 0.89379 -0.009282 0.000672 -3.7996e-05 +0.26395 0.34429 1.003 -0.011328 0.001324 0 +0.30717 0.35778 1.0569 -0.012171 0.001498 0 +0.20741 0.34242 1.0166 -0.011613 0.001251 -0.00023799 +0.16655 0.34197 0.9551 -0.011191 0.000738 -0.00014699 +0.19972 0.34108 0.96788 -0.011348 0.001124 -6.0999e-05 +0.26775 0.35173 1.0466 -0.011911 0.001321 -0.00022201 +0.25279 0.35938 1.066 -0.01237 0.001314 -0.00035401 +0.1816 0.36449 1.0286 -0.011494 0.00109 -0.000335 +0.17043 0.35634 1.0007 -0.011367 0.000676 -0.00038699 +0.20961 0.37343 0.9861 -0.009674 0.00063699 -0.000376 +0.28539 0.38968 1.0834 -0.012728 0.001033 -0.000534 +0.24985 0.37998 1.0853 -0.012493 0.001143 -0.00048799 +0.17213 0.37133 1.0384 -0.011699 0.000652 -0.000501 +0.1827 0.349 0.96999 -0.010547 0.000663 -0.00039699 +0.25229 0.37061 1.0672 -0.011505 0.001172 -0.000413 +0.2795 0.37499 1.1003 -0.01256 0.001254 -0.000508 +0.19133 0.37068 0.98973 -0.009776 0.001089 0 +0.18534 0.37375 1.0025 -0.010219 0.00073499 -0.00035401 +0.24973 0.38462 1.0698 -0.01149 0.001201 0 +0.27635 0.40118 1.1104 -0.012286 0.001222 -0.00048201 +0.20479 0.38131 1.0699 -0.011708 0.001063 -0.000531 +0.1632 0.36498 1.0014 -0.010743 0.000673 -0.000433 +0.19867 0.37315 1.0346 -0.011125 0.000802 -0.00045099 +0.28956 0.39603 1.0912 -0.011177 0.001361 0 +0.24666 0.39424 1.0299 -0.01041 0.001057 0 +0.18383 0.38167 1.0273 -0.010033 0.000797 -0.000477 +0.21611 0.38453 1.0743 -0.010675 0.00119 -0.00045799 +0.29614 0.4038 1.1428 -0.011934 0.001544 -0.000581 +0.28162 0.40364 1.1439 -0.012353 0.001523 -0.000574 +0.18819 0.38217 1.103 -0.011635 0.00095799 -0.00050601 +0.16648 0.36827 1.021 -0.010201 0.001056 -0.000579 +0.22544 0.37657 1.0959 -0.011561 0.001219 -0.000521 +0.295 0.39491 1.1524 -0.012521 0.001402 -0.00058001 +0.24581 0.39762 1.1279 -0.011882 0.001509 -0.00060701 +0.19326 0.3771 1.061 -0.010701 0.001312 -0.00056899 +0.21144 0.36365 0.99964 -0.010159 0.001353 0 +0.21814 0.36581 0.98616 -0.009655 0.00112 -0.000619 +0.3098 0.39963 1.1381 -0.011288 0.001267 0 +0.28386 0.41003 1.1333 -0.011344 0.001182 -0.00065801 +0.16129 0.35734 0.98427 -0.009402 0.000612 -0.000933 +0.13255 0.3373 0.86953 -0.00845 0.000835 -0.001076 +0.19388 0.33661 0.9661 -0.009374 0.00241 -0.001002 +0.29255 0.3556 1.1254 -0.01101 0.003261 -0.001273 +0.26009 0.34176 1.0838 -0.011142 0.0022 -0.00147 +0.13391 0.30184 0.84955 -0.00927 0 -0.001931 +0.084206 0.26009 0.73414 -0.007952 0 -0.002533 +0.11811 0.28837 0.72244 -0.008263 0 -0.002417 +0.22179 0.32535 0.94714 -0.01078 0.001052 -0.002401 +0.21373 0.34604 1.0459 -0.01239 0.001089 -0.002027 +0.14041 0.35453 1.0242 -0.011133 0.00069 -0.001802 +0.083953 0.3476 1.0127 -0.010621 0.000457 -0.00157 +0.085777 0.33048 0.95808 -0.010109 0.000541 -0.001379 +0.17204 0.36064 1.0504 -0.011299 0.000822 -0.00142 +0.22859 0.38075 1.0914 -0.011893 0.000879 -0.001479 +0.17234 0.37485 1.0499 -0.011506 0.000653 -0.001357 +0.15794 0.35836 0.98568 -0.009627 0.00036701 -0.001411 +0.18087 0.34854 0.93225 -0.009462 0.00046099 -0.001278 +0.22494 0.3835 1.021 -0.010119 0.000829 -0.001396 +0.25899 0.38602 1.1191 -0.01168 0.000791 -0.001673 +0.23389 0.39802 1.1224 -0.011871 0.000713 -0.001663 +0.16196 0.39507 1.0998 -0.010818 0 -0.001536 +0.11145 0.37243 1.0414 -0.010107 0 -0.001327 +0.16101 0.38643 1.0558 -0.010382 0 -0.001408 +0.2273 0.3966 1.1196 -0.011769 0.000656 -0.001523 +0.19045 0.39704 1.1313 -0.01203 0.000581 -0.001434 +0.15087 0.37955 1.1065 -0.011323 0 -0.001378 +0.16281 0.36617 1.0502 -0.01083 0 -0.001362 +0.23656 0.39159 1.0929 -0.010923 0 -0.001398 +0.28873 0.39128 1.152 -0.01201 0.00069201 -0.001578 +0.20254 0.39104 1.124 -0.011167 0.000595 -0.00149 +0.17305 0.37228 1.0773 -0.010623 0 -0.001301 +0.20931 0.38511 1.0883 -0.010626 0 -0.001238 +0.26161 0.40861 1.1696 -0.011835 0.00063999 -0.001447 +0.20246 0.39894 1.1387 -0.011292 0.000521 -0.00141 +0.13837 0.39171 1.108 -0.010845 0 -0.001313 +0.14967 0.37167 1.034 -0.009836 0 -0.001162 +0.21058 0.38885 1.103 -0.010879 0.00052601 -0.001296 +0.25663 0.40713 1.1468 -0.011529 0.000618 -0.001453 +0.21847 0.41151 1.0961 -0.010484 0 -0.001351 +0.17418 0.39165 1.0759 -0.010233 0 -0.001295 +0.20315 0.3904 1.0803 -0.010106 0 -0.001267 +0.25852 0.41606 1.1381 -0.011101 0 -0.001313 +0.26305 0.40965 1.171 -0.011753 0.000575 -0.001377 +0.18057 0.4049 1.1175 -0.010561 0 -0.001295 +0.16957 0.39486 1.0748 -0.010293 0 -0.001393 +0.19342 0.39769 1.0332 -0.009033 0 -0.001332 +0.22634 0.404 1.1141 -0.011036 0 -0.001478 +0.25934 0.41296 1.1667 -0.011678 0.000649 -0.001503 +0.23604 0.41179 1.1662 -0.011784 0.000579 -0.001425 +0.15938 0.40422 1.1386 -0.01087 0 -0.001338 +0.18757 0.39965 1.1191 -0.010359 0 -0.001309 +0.21948 0.41622 1.1226 -0.010207 0.000784 -0.001301 +0.28036 0.42663 1.1866 -0.011304 0.000743 -0.001324 +0.22035 0.42978 1.1813 -0.011419 0.00059701 -0.001361 +0.17767 0.41926 1.1277 -0.010219 0 -0.001354 +0.19654 0.41566 1.0665 -0.009388 0 -0.001261 +0.20423 0.41758 1.0677 -0.009551 0.00089801 -0.001333 +0.25191 0.41595 1.1979 -0.010699 0.002221 -0.001954 +0.20027 0.40135 1.0321 -0.009192 0.001173 -0.001791 +0.050058 0.36459 0.80188 -0.007302 0 -0.001651 +0.085013 0.32969 0.73033 -0.005606 0 -0.001547 +0.092794 0.30763 0.67653 -0.005713 0 -0.00149 +0.19082 0.32145 0.85298 -0.007379 0 -0.001651 +0.30919 0.34608 1.1343 -0.009763 0.002281 -0.002044 +0.32184 0.37877 1.1873 -0.010776 0.002081 -0.0019 +0.23376 0.38324 1.1738 -0.01066 0.001552 -0.001599 +0.13037 0.36278 1.1257 -0.009633 0.001301 -0.001432 +0.11775 0.36262 1.0601 -0.008867 0.001282 -0.001357 +0.18446 0.38565 1.1597 -0.00969 0.00177 -0.001536 +0.19899 0.39645 1.2115 -0.010291 0.001653 -0.001531 +0.15776 0.39378 1.2065 -0.010331 0.001482 -0.001415 +0.12262 0.37602 1.1333 -0.008736 0.001531 -0.001342 +0.15588 0.38069 1.0286 -0.007089 0.001259 -0.001468 +0.24303 0.39426 1.1398 -0.00945 0.001499 -0.00165 +0.28034 0.41213 1.2058 -0.009906 0.001767 -0.001628 +0.22658 0.41032 1.1832 -0.009645 0.001683 -0.001642 +0.16413 0.41334 1.1088 -0.00866 0.001332 -0.001673 +0.17743 0.4162 1.1 -0.008279 0.001408 -0.001598 +0.25137 0.40625 1.1889 -0.009104 0.001576 -0.00161 +0.23554 0.44268 1.2165 -0.009156 0.001528 -0.001752 +0.16633 0.42795 1.2014 -0.008909 0.0015 -0.00167 +0.13181 0.40452 1.1325 -0.008187 0.001433 -0.001636 +0.17307 0.42174 1.152 -0.007581 0.00158 -0.001603 +0.2504 0.44603 1.2574 -0.008829 0.001854 -0.00177 +0.23506 0.4352 1.2681 -0.008928 0.002144 -0.001837 +0.1549 0.4271 1.0823 -0.007608 0.00157 -0.002006 +0.087461 0.3888 0.89855 -0.005882 0.001158 -0.00172 +0.20834 0.40858 1.0964 -0.007397 0.002421 -0.001957 +0.26173 0.42389 1.1686 -0.008088 0.002787 -0.002142 +0.11227 0.37232 0.91343 -0.006503 0.001238 -0.002018 +0.065428 0.33485 0.74526 -0.005679 0.00070201 -0.001628 +0.13896 0.3314 0.89099 -0.00743 0.001958 -0.001787 +0.2222 0.34461 1.0456 -0.008921 0.003132 -0.002127 +0.18179 0.32273 0.95374 -0.007771 0.00311 -0.002364 +0.076884 0.29079 0.84247 -0.006636 0.001677 -0.002811 +0.10286 0.30589 0.84769 -0.006368 0.001586 -0.002914 +0.11714 0.31505 0.8283 -0.006453 0.001399 -0.002614 +0.17289 0.31754 0.91644 -0.008077 0.001941 -0.002466 +0.16636 0.34556 1.0393 -0.009494 0.001801 -0.002268 +0.12452 0.35842 1.0775 -0.009579 0.001541 -0.001987 +0.061172 0.35592 1.0302 -0.008875 0.001451 -0.001755 +0.050946 0.33753 0.94083 -0.007463 0.001408 -0.00167 +0.10546 0.3456 0.93439 -0.007485 0.000764 -0.001549 +0.21463 0.36993 1.0284 -0.008486 0.001255 -0.001934 +0.21658 0.37702 1.089 -0.009459 0.00162 -0.001807 +0.16192 0.37211 1.0661 -0.009051 0.001121 -0.001793 +0.11352 0.36468 0.96662 -0.007048 0.001228 -0.001675 +0.14497 0.36237 0.94079 -0.006657 0.001226 -0.001595 +0.2295 0.38314 1.0477 -0.008416 0.001467 -0.00185 +0.21005 0.39357 1.0733 -0.008476 0.00132 -0.002017 +0.10958 0.39631 1.0276 -0.007353 0.001215 -0.0017 +0.10539 0.38754 1.0054 -0.006428 0.001325 -0.001462 +0.16619 0.40653 1.0907 -0.007926 0.001821 -0.001669 +0.22847 0.41654 1.1459 -0.008573 0.001966 -0.001877 +0.18655 0.42015 1.1381 -0.00866 0.001812 -0.001818 +0.12523 0.40908 1.0344 -0.006052 0.001202 -0.00179 +0.15265 0.40805 1.0602 -0.006312 0.001477 -0.001601 +0.22654 0.40899 1.1499 -0.008054 0.001983 -0.001673 +0.26632 0.42773 1.1934 -0.008883 0.002092 -0.001888 +0.19458 0.42155 1.1609 -0.008232 0.001966 -0.001733 +0.15037 0.38987 1.0941 -0.007299 0.001739 -0.001618 +0.17838 0.43768 1.038 -0.005316 0.001412 -0.001672 +0.22498 0.44814 1.1553 -0.007702 0.001979 -0.001768 +0.25122 0.43727 1.2043 -0.008749 0.002013 -0.001806 +0.17543 0.43428 1.1678 -0.008157 0.001722 -0.001676 +0.14439 0.41772 1.0905 -0.006486 0.001523 -0.001586 +0.19785 0.42324 1.1539 -0.007429 0.002026 -0.001644 +0.27134 0.43939 1.2337 -0.008884 0.002211 -0.001702 +0.22665 0.43418 1.214 -0.008735 0.002183 -0.001659 +0.15686 0.41922 1.1477 -0.007365 0.002043 -0.001587 +0.1734 0.42089 1.0996 -0.006672 0.001883 -0.001691 +0.23291 0.42332 1.1927 -0.007522 0.00244 -0.001594 +0.27039 0.46045 1.2207 -0.007945 0.002234 -0.001917 +0.2104 0.42863 1.1971 -0.007466 0.002168 -0.001621 +0.14273 0.43285 1.0969 -0.00595 0.001726 -0.001683 +0.11731 0.40773 0.9284 -0.004787 0.001524 -0.001812 +0.11753 0.42059 1.0124 -0.005708 0.002213 -0.001878 +0.33155 0.4411 1.2042 -0.008714 0.003893 -0.000551 +0.2735 0.38752 1.0353 -0.008438 0.003346 -0.000558 +0.12354 0.3646 0.78382 -0.003869 0 -0.001781 +0.12412 0.34321 0.70618 -0.003609 0.000514 -0.002007 +0.174 0.34832 0.76286 -0.003627 0.00095 -0.002104 +0.2325 0.37061 0.93324 -0.00394 0.0036 -0.00233 +0.27049 0.36935 1.0613 -0.005519 0.0036 -0.002205 +0.26344 0.38938 1.111 -0.007793 0.00328 -0.002038 +0.18119 0.39119 1.1208 -0.007792 0.003017 -0.001805 +0.094057 0.37518 1.0702 -0.006774 0.002759 -0.00149 +0.092936 0.39383 0.95362 -0.003604 0.002121 -0.001349 +0.16888 0.42759 1.0053 -0.004388 0.002107 -0.001473 +0.22481 0.41919 1.1507 -0.007272 0.002943 -0.001535 +0.23362 0.4182 1.1547 -0.007957 0.002799 -0.001651 +0.15638 0.4139 1.0669 -0.005727 0.002399 -0.001387 +0.16246 0.40523 1.0166 -0.005113 0.002297 -0.001409 +0.24245 0.42402 1.111 -0.006359 0.002675 -0.0015 +0.26956 0.44013 1.1631 -0.006955 0.002949 -0.00165 +0.2086 0.44417 1.1196 -0.005565 0.002764 -0.001607 +0.15802 0.44393 1.0375 -0.004515 0.002591 -0.001538 +0.17562 0.45818 1.0112 -0.00424 0.002667 -0.001459 +0.19782 0.46908 1.1115 -0.004879 0.003193 -0.001473 +0.23439 0.47891 1.201 -0.006486 0.002738 -0.001762 +0.2198 0.46361 1.2203 -0.00717 0.002733 -0.001785 +0.14675 0.47504 1.2017 -0.00554 0.002401 -0.001789 +0.13499 0.45559 1.0721 -0.003925 0.002224 -0.001715 +0.15846 0.46558 1.0267 -0.003262 0.002203 -0.001697 +0.2236 0.48748 1.1479 -0.00482 0.002747 -0.001849 +0.23652 0.49111 1.2226 -0.006213 0.002948 -0.001917 +0.16335 0.46588 1.1861 -0.005943 0.002681 -0.00179 +0.14012 0.46939 1.1027 -0.004315 0.002699 -0.0018 +0.17344 0.47591 1.1351 -0.004141 0.002696 -0.001737 +0.21344 0.49619 1.2137 -0.004799 0.0028 -0.002023 +0.21515 0.48497 1.2011 -0.005123 0.00291 -0.001992 +0.13788 0.45528 1.1142 -0.004316 0.002686 -0.001981 +0.12705 0.42475 1.0263 -0.004075 0.003058 -0.00196 +0.17233 0.40224 0.98847 -0.003842 0.002911 -0.002464 +0.24187 0.40368 1.1317 -0.005194 0.00409 -0.002337 +0.36001 0.42296 1.2548 -0.006726 0.005698 -0.002381 +0.28489 0.40384 1.2053 -0.00632 0.004553 -0.002553 +0.1612 0.35921 0.96176 -0.004519 0.003772 -0.002802 +0.13385 0.3546 0.91526 -0.003883 0.002942 -0.002682 +0.12899 0.329 0.84164 -0.004142 0.002501 -0.002909 +0.15389 0.35616 0.99727 -0.004644 0.003168 -0.003109 +0.228 0.38192 1.1834 -0.006218 0.004399 -0.003477 +0.20533 0.37857 1.155 -0.005954 0.004342 -0.003457 +0.1366 0.35231 1.0538 -0.005041 0.003554 -0.003341 +0.11045 0.33976 0.94221 -0.004448 0.002845 -0.003276 +0.010287 0.29439 0.74133 -0.00421 0.002016 -0.003802 +-0.043041 0.31391 0.77775 -0.005181 0.00097601 -0.004615 +0.021874 0.32569 0.97559 -0.007103 0.001006 -0.005059 +0.13111 0.35029 1.1617 -0.00822 0.003218 -0.004786 +0.14282 0.35219 1.1281 -0.007411 0.00364 -0.004601 +0.15219 0.3918 1.1101 -0.007082 0.003688 -0.004448 +0.11137 0.41613 1.1006 -0.006329 0.002511 -0.004219 +0.083074 0.45649 1.096 -0.005079 0.001908 -0.00368 +-0.005676 0.42784 1.0705 -0.005098 0.002233 -0.003036 +-0.040277 0.37715 0.9058 -0.004016 0.002202 -0.002585 +0.031928 0.38708 0.842 -0.002712 0.002428 -0.002819 +0.13564 0.43061 0.75952 0.003334 0.00241 -0.004088 +0.12564 0.4331 0.83321 0.002665 0.002435 -0.003822 +0.067414 0.40951 0.73249 0.002951 0.002658 -0.003521 +0.061024 0.40471 0.64347 0.002337 0.003043 -0.002925 +0.15602 0.40884 0.69009 0.002557 0.00339 -0.003461 +0.21379 0.43569 0.82249 0.001454 0.003849 -0.003436 +0.14148 0.44098 0.81816 0.001526 0.003757 -0.003868 +0.051739 0.41998 0.74244 0.002173 0.003575 -0.00316 +0.076076 0.42541 0.7142 0.00241 0.003281 -0.003098 +0.16366 0.46775 0.82032 0.001537 0.003852 -0.003578 +0.12569 0.44842 0.87724 0.000687 0.004157 -0.003404 +0.064304 0.44117 0.82401 0.001768 0.003668 -0.003189 +0.05466 0.44401 0.75462 0.002955 0.003394 -0.003244 +0.12669 0.45651 0.80066 0.002351 0.003167 -0.003439 +0.20428 0.46561 0.9098 0.000875 0.003862 -0.003875 +0.15332 0.45971 0.92827 0.000852 0.003891 -0.00349 +0.093259 0.45788 0.86544 0.002165 0.003713 -0.003355 +0.10124 0.47666 0.81648 0.003211 0.003592 -0.003324 +0.16026 0.47363 0.92621 0.001576 0.003883 -0.003542 +0.19695 0.48667 0.99445 0.001105 0.004699 -0.003865 +0.10507 0.47682 0.96353 0.001159 0.004056 -0.003445 +0.067366 0.4773 0.8886 0.00173 0.00422 -0.003197 +0.11421 0.47937 0.88604 0.00228 0.003868 -0.003208 +0.19268 0.48474 0.94618 0.00075499 0.004311 -0.003675 +0.16679 0.48844 0.94449 0.001238 0.004263 -0.003585 +0.084298 0.47303 0.86638 0.001842 0.004355 -0.003195 +0.093624 0.4663 0.84511 0.002517 0.003925 -0.00312 +0.15411 0.47669 0.9109 0.001824 0.004581 -0.003518 +0.21448 0.48722 0.9849 0.000821 0.004815 -0.00348 +0.16195 0.48502 0.9786 0.00078599 0.004625 -0.003525 +0.080036 0.48048 0.91742 0.002569 0.004089 -0.003158 +0.093059 0.47656 0.86332 0.002836 0.003877 -0.003132 +0.16051 0.47558 0.92787 0.001592 0.004352 -0.00334 +0.17313 0.49602 0.96093 0.001182 0.004276 -0.003609 +0.11498 0.48345 0.9472 0.001559 0.004389 -0.003279 +0.076875 0.47833 0.89154 0.002343 0.004219 -0.003099 +0.13095 0.48328 0.90893 0.002068 0.004288 -0.003255 +0.19627 0.49089 0.98319 0.000802 0.005015 -0.003652 +0.12664 0.48981 0.95569 0.001282 0.004456 -0.003441 +0.080994 0.48361 0.89122 0.002491 0.004274 -0.003104 +0.11419 0.49937 0.88198 0.002853 0.004263 -0.003179 +0.197 0.50892 1.0258 0.00038699 0.005013 -0.003596 +0.16532 0.51165 1.0042 0.00073 0.005415 -0.003616 +0.090297 0.49795 0.9745 0.001305 0.004382 -0.003211 +0.094869 0.49474 0.91553 0.002383 0.004429 -0.003369 +0.12328 0.49972 0.92179 0.002246 0.00507 -0.003316 +0.19073 0.50951 1.0283 0.00076601 0.005315 -0.003768 +0.17335 0.51461 1.0346 0.00053601 0.005255 -0.003574 +0.081566 0.5022 0.94942 0.001791 0.005032 -0.003211 +0.10162 0.49843 0.90054 0.002544 0.004438 -0.003379 +0.11869 0.49843 0.87451 0.002767 0.004699 -0.003273 +0.18337 0.50678 1.0071 0.000961 0.005143 -0.003581 +0.20162 0.52346 1.0309 0.00030199 0.005458 -0.003518 +0.11035 0.5156 0.99661 0.000649 0.004492 -0.003401 +0.08544 0.5118 0.94224 0.001956 0.005058 -0.00335 +0.099513 0.50257 0.91799 0.002601 0.004868 -0.003221 +0.13257 0.52043 0.96804 0.002134 0.005223 -0.003526 +0.19207 0.52865 1.0432 0.00090801 0.005455 -0.003552 +0.15437 0.52524 1.0403 0.00019099 0.005341 -0.00352 +0.097913 0.51669 0.97684 0.001582 0.005266 -0.003498 +0.087629 0.51348 0.92812 0.002633 0.004983 -0.0035 +0.13445 0.52364 0.92976 0.002936 0.00499 -0.003357 +0.18423 0.53838 1.007 0.00073601 0.005109 -0.003723 +0.13844 0.53037 1.0224 0.00088399 0.004745 -0.003666 +0.071467 0.52143 0.9585 0.002174 0.004531 -0.003419 +0.09417 0.51777 0.92469 0.002629 0.005319 -0.003394 +0.15668 0.52904 0.99011 0.001796 0.005512 -0.003755 +0.1858 0.53553 1.042 0.00077199 0.005434 -0.003581 +0.095012 0.52103 0.96853 0.001876 0.00447 -0.003503 +0.073715 0.51191 0.90569 0.002867 0.004306 -0.003525 +0.1302 0.51986 0.98143 0.001941 0.004293 -0.003926 +0.13781 0.53142 1.0407 0.000858 0.006095 -0.003829 +0.12917 0.51693 1.0719 0 0.005659 -0.003786 +0.11882 0.49919 1.0382 0 0.005865 -0.003911 +0.12795 0.49167 1.0125 0 0.006105 -0.004055 +0.10296 0.46335 1.0425 0 0.006204 -0.004183 +0.17929 0.43882 1.2425 -0.002894 0.007712 -0.004577 +0.24133 0.47048 1.1931 -0.002025 0.007457 -0.004758 +0.133 0.46633 1.1452 -0.001894 0.007263 -0.004229 +0.14882 0.47794 1.1503 -0.002085 0.007522 -0.004076 +0.15653 0.48625 1.1869 -0.003012 0.007549 -0.00415 +0.097475 0.47729 1.1772 -0.002454 0.007448 -0.003978 +0.023815 0.46435 1.157 -0.000933 0.007595 -0.00327 +0.063799 0.45183 1.0843 -0.000639 0.007352 -0.003623 +0.15299 0.47622 1.1245 -0.001499 0.007698 -0.003949 +0.19555 0.48124 1.1874 -0.002572 0.007795 -0.004387 +0.16097 0.47489 1.1746 -0.001875 0.00753 -0.004249 +0.10831 0.47654 1.1178 -0.000476 0.00758 -0.003885 +0.14816 0.46877 1.0998 -0.00084401 0.007268 -0.004128 +0.2245 0.48643 1.1579 -0.002008 0.007329 -0.004477 +0.18884 0.4951 1.19 -0.001943 0.007559 -0.004392 +0.10019 0.49439 1.1752 -0.001165 0.007178 -0.0041 +0.085649 0.48291 1.1298 -0.000513 0.006957 -0.003782 +0.12739 0.48067 1.1642 -0.001282 0.007564 -0.003972 +0.20689 0.49499 1.2569 -0.002628 0.007221 -0.004247 +0.17588 0.49868 1.243 -0.002734 0.007068 -0.004225 +0.10342 0.48814 1.2123 -0.00213 0.007242 -0.003884 +0.11812 0.47138 1.1454 -0.00116 0.007766 -0.00389 +0.16806 0.4818 1.1864 -0.001378 0.008046 -0.004055 +0.2413 0.50101 1.2607 -0.00258 0.008232 -0.004245 +0.18763 0.51147 1.2406 -0.002169 0.007503 -0.004382 +0.11744 0.49158 1.2086 -0.001358 0.007944 -0.003925 +0.12757 0.47705 1.1594 -0.00061401 0.008098 -0.004006 +0.1965 0.49118 1.2317 -0.00157 0.008398 -0.004206 +0.22615 0.50341 1.2953 -0.002547 0.008419 -0.00432 +0.1462 0.49893 1.2661 -0.002167 0.008265 -0.004166 +0.094569 0.4848 1.2202 -0.001154 0.008191 -0.003783 +0.14294 0.48453 1.2125 -0.001252 0.008439 -0.003807 +0.21488 0.50389 1.2597 -0.002431 0.008395 -0.004224 +0.20989 0.50308 1.2914 -0.002814 0.00853 -0.004253 +0.10876 0.49375 1.2438 -0.001629 0.008328 -0.003948 +0.12247 0.4879 1.2128 -0.001029 0.008258 -0.003851 +0.21095 0.49945 1.2744 -0.002063 0.008522 -0.004171 +0.22065 0.50901 1.2999 -0.002796 0.008561 -0.004314 +0.14976 0.50267 1.2901 -0.002275 0.008601 -0.004155 +0.099643 0.49259 1.2271 -0.00088101 0.00854 -0.003877 +0.15193 0.49856 1.236 -0.0012 0.008584 -0.004042 +0.23034 0.51871 1.3186 -0.002715 0.008662 -0.004241 +0.17869 0.51501 1.3099 -0.002721 0.008573 -0.004086 +0.12322 0.50238 1.2818 -0.001596 0.008534 -0.003825 +0.1337 0.49449 1.2162 -0.001155 0.00857 -0.003815 +0.20797 0.51058 1.294 -0.002339 0.008683 -0.003984 +0.2205 0.51795 1.3255 -0.003754 0.008711 -0.004208 +0.13594 0.49785 1.3057 -0.003399 0.008316 -0.003921 +0.12916 0.49414 1.2454 -0.001597 0.008754 -0.003991 +0.17115 0.50461 1.2684 -0.00161 0.008803 -0.004019 +0.25021 0.5231 1.3572 -0.002934 0.008557 -0.004218 +0.21589 0.5259 1.3614 -0.004304 0.008654 -0.004252 +0.13075 0.51235 1.3092 -0.002193 0.0087 -0.004023 +0.13125 0.50603 1.253 -0.001317 0.008829 -0.003901 +0.18323 0.52252 1.2672 -0.001662 0.008806 -0.003928 +0.25338 0.53357 1.3393 -0.003032 0.009069 -0.00386 +0.22018 0.53099 1.3334 -0.003032 0.008955 -0.003776 +0.14089 0.51845 1.2919 -0.001768 0.008785 -0.003811 +0.1545 0.51304 1.2493 -0.001187 0.008771 -0.003761 +0.21868 0.52326 1.322 -0.002404 0.008807 -0.003823 +0.24193 0.52927 1.3565 -0.003252 0.008627 -0.004051 +0.17285 0.52929 1.2938 -0.002535 0.007911 -0.003632 +0.12721 0.51195 1.2604 -0.001207 0.008255 -0.00393 +0.18006 0.5086 1.2833 -0.00184 0.008443 -0.00384 +0.2485 0.53598 1.3535 -0.00328 0.008744 -0.004014 +0.20439 0.52104 1.3422 -0.003168 0.008708 -0.004062 +0.13264 0.51526 1.2985 -0.002006 0.008616 -0.003895 +0.15082 0.50584 1.2575 -0.001451 0.008677 -0.003772 +0.22752 0.51875 1.3177 -0.002408 0.008835 -0.003891 +0.25508 0.53249 1.3606 -0.003202 0.008645 -0.003776 +0.18426 0.53539 1.3519 -0.003135 0.008591 -0.003774 +0.14491 0.5159 1.3009 -0.001692 0.008496 -0.003609 +0.17253 0.52761 1.2542 -0.000535 0.007875 -0.003385 +0.2482 0.52324 1.3392 -0.002233 0.008753 -0.003548 +0.25632 0.53385 1.3838 -0.003166 0.009014 -0.00383 +0.16374 0.51405 1.3403 -0.002004 0.009028 -0.003615 +0.17456 0.51147 1.2977 -0.001311 0.00886 -0.003693 +0.20938 0.51767 1.3291 -0.001714 0.008938 -0.003509 +0.26901 0.53765 1.373 -0.003059 0.008809 -0.003767 +0.20194 0.53963 1.3279 -0.002255 0.008369 -0.003484 +0.15048 0.52147 1.3109 -0.001549 0.008255 -0.003671 +0.17302 0.5128 1.2952 -0.001129 0.00841 -0.003513 +0.22774 0.5275 1.3462 -0.001982 0.008409 -0.003604 +0.25794 0.53119 1.3675 -0.003063 0.008511 -0.00366 +0.17568 0.52407 1.329 -0.001669 0.008377 -0.003484 +0.16128 0.51501 1.2637 -0.00034001 0.008401 -0.003401 +0.21631 0.52397 1.2958 -0.00087101 0.008493 -0.003453 +0.26112 0.5427 1.3293 -0.001739 0.008522 -0.003565 +0.20969 0.52056 1.3392 -0.002316 0.008475 -0.003623 +0.16557 0.51834 1.3034 -0.00089 0.008382 -0.003502 +0.2549 0.52678 1.3546 -0.001929 0.008434 -0.003322 +0.25559 0.52855 1.3932 -0.002949 0.008179 -0.003349 +0.16815 0.51662 1.3455 -0.00148 0.00754 -0.003445 +0.15315 0.5077 1.2776 -0.000497 0.007522 -0.003416 +0.19814 0.51089 1.272 -0.00055301 0.007364 -0.003179 +0.27765 0.52339 1.3558 -0.002143 0.008481 -0.003306 +0.25411 0.52765 1.3685 -0.002434 0.008337 -0.003221 +0.14475 0.5138 1.2902 -0.00081701 0.007612 -0.003365 +0.15699 0.5042 1.1996 9.9484e-07 0.007279 -0.003116 +0.21572 0.51732 1.2779 -0.000951 0.008105 -0.003242 +0.284 0.53235 1.3405 -0.001838 0.008341 -0.003436 +0.22137 0.52962 1.3296 -0.001479 0.007896 -0.003317 +0.14697 0.52157 1.2708 0.000102 0.007534 -0.003304 +0.18768 0.5169 1.255 0.00010001 0.00795 -0.003076 +0.26192 0.53114 1.335 -0.001254 0.008174 -0.003256 +0.2577 0.53485 1.3551 -0.001814 0.007681 -0.003183 +0.16636 0.54185 1.2767 -0.00010299 0.007473 -0.003259 +0.17575 0.5209 1.2491 0.00034001 0.007725 -0.003043 +0.23139 0.5229 1.2939 -0.00026499 0.007981 -0.003164 +0.31157 0.53901 1.3878 -0.002225 0.008369 -0.003181 +0.25771 0.53896 1.3747 -0.002029 0.008141 -0.003111 +0.16862 0.52707 1.2858 -0.00027601 0.007973 -0.003083 +0.16941 0.51091 1.2207 0.000585 0.007856 -0.003108 +0.22731 0.52667 1.2876 -0.00022 0.008124 -0.003187 +0.28878 0.5391 1.3654 -0.00179 0.008192 -0.003188 +0.22583 0.53728 1.3188 -0.001319 0.007556 -0.003179 +0.17465 0.5252 1.2602 3.7996e-05 0.00807 -0.003099 +0.19423 0.52025 1.2491 0.00049 0.008147 -0.003095 +0.27398 0.53314 1.3046 -0.000622 0.008292 -0.002891 +0.28239 0.54113 1.344 -0.001894 0.008306 -0.003225 +0.20353 0.53245 1.3073 -0.00089101 0.008011 -0.003123 +0.16647 0.51353 1.2241 0.00040101 0.007315 -0.003038 +0.21668 0.52528 1.2349 7.1995e-05 0.007429 -0.00297 +0.29674 0.53999 1.3201 -0.001383 0.007831 -0.0031 +0.23893 0.53787 1.3219 -0.001339 0.008182 -0.003092 +0.16887 0.52778 1.2724 -0.00017101 0.007497 -0.002933 +0.16488 0.52357 1.1851 0.00071099 0.007336 -0.002863 +0.23033 0.52431 1.2376 2.0001e-05 0.007925 -0.002892 +0.30156 0.54432 1.3215 -0.001255 0.007731 -0.003141 +0.23847 0.54511 1.3022 -0.00056701 0.007578 -0.003129 +0.17492 0.52138 1.2723 -0.000131 0.008132 -0.003011 +0.20782 0.51766 1.2534 0.00015999 0.008021 -0.002963 +0.2723 0.53626 1.3342 -0.00109 0.008408 -0.003088 +0.29465 0.54441 1.3776 -0.001848 0.008475 -0.003127 +0.18307 0.53969 1.3344 -0.000943 0.007784 -0.003042 +0.196 0.53072 1.2851 -8.8994e-05 0.008165 -0.002833 +0.21937 0.53289 1.2548 0.000186 0.008127 -0.002917 +0.27562 0.551 1.3519 -0.001251 0.008414 -0.003023 +0.29409 0.55588 1.3758 -0.00214 0.008386 -0.003021 +0.22485 0.55106 1.331 -0.001105 0.008305 -0.002949 +0.21935 0.54782 1.2962 -8.4998e-05 0.008285 -0.002946 +0.25138 0.54697 1.3071 -0.000311 0.008483 -0.002997 +0.30052 0.55396 1.3782 -0.001902 0.008557 -0.003017 +0.30058 0.55433 1.3775 -0.002405 0.008568 -0.002989 +0.2118 0.53888 1.3422 -0.001434 0.008478 -0.002966 +0.21081 0.53263 1.3236 -0.00069201 0.008434 -0.002845 +0.26565 0.54254 1.3633 -0.00149 0.008554 -0.002908 +0.30319 0.55943 1.4124 -0.002828 0.008133 -0.003126 +0.25265 0.54582 1.3938 -0.002691 0.008022 -0.003112 +0.18725 0.53766 1.3601 -0.001251 0.008486 -0.002981 +0.20614 0.51704 1.2785 -0.000524 0.008158 -0.00278 +0.26495 0.53584 1.3299 -0.00099 0.008652 -0.003048 +0.3315 0.55181 1.4087 -0.002372 0.008789 -0.003183 +0.28596 0.54934 1.4182 -0.002691 0.008739 -0.003146 +0.21321 0.5369 1.3489 -0.00085699 0.008646 -0.003063 +0.2298 0.52793 1.2947 -0.00020399 0.008798 -0.002977 +0.28238 0.53891 1.3609 -0.001269 0.008948 -0.003018 +0.33145 0.54732 1.4212 -0.002766 0.008917 -0.003068 +0.25318 0.54141 1.3826 -0.002374 0.008393 -0.003091 +0.18811 0.53112 1.2966 -0.000325 0.009021 -0.002993 +0.20929 0.52731 1.2589 0.000125 0.009024 -0.003039 +0.29279 0.54135 1.3576 -0.001243 0.009297 -0.003122 +0.3252 0.56612 1.4057 -0.002715 0.009191 -0.003153 +0.24431 0.54169 1.3563 -0.001622 0.009174 -0.003156 +0.21472 0.52656 1.2688 -0.000304 0.009159 -0.003048 +0.25154 0.53292 1.2827 -0.000146 0.009106 -0.003178 +0.33248 0.55225 1.3912 -0.001745 0.009228 -0.00326 +0.27921 0.5485 1.3849 -0.001994 0.008873 -0.003199 +0.20865 0.53841 1.3381 -0.000966 0.008837 -0.003218 +0.19932 0.52553 1.2903 -0.00019 0.008893 -0.003188 +0.23914 0.54272 1.2928 -5.7002e-05 0.008293 -0.003198 +0.30904 0.54746 1.3762 -0.001857 0.008983 -0.003328 +0.25565 0.54744 1.3648 -0.001892 0.008755 -0.003282 +0.20171 0.5368 1.3128 -0.000689 0.008867 -0.003173 +0.20489 0.52903 1.2772 7.4997e-05 0.008829 -0.003141 +0.29731 0.54753 1.3614 -0.001262 0.009054 -0.003218 +0.30313 0.55655 1.3963 -0.00204 0.008929 -0.003215 +0.2262 0.54986 1.3707 -0.001382 0.00872 -0.00318 +0.21226 0.53418 1.3049 -7.4002e-05 0.008918 -0.003175 +0.23264 0.54025 1.3156 -0.00031599 0.008926 -0.003153 +0.31577 0.55889 1.3887 -0.001675 0.00892 -0.0032 +0.2758 0.56081 1.3839 -0.001761 0.008921 -0.003281 +0.21775 0.55364 1.3157 -5.6008e-05 0.008872 -0.00312 +0.24527 0.56147 1.3294 -0.000186 0.009003 -0.003156 +0.31357 0.57156 1.4254 -0.001968 0.008935 -0.003271 +0.30058 0.575 1.428 -0.002297 0.008951 -0.00333 +0.22515 0.56275 1.3791 -0.001063 0.009064 -0.003223 +0.22317 0.55026 1.315 -6.4001e-05 0.009127 -0.003112 +0.2829 0.56275 1.3699 -0.001041 0.009113 -0.003207 +0.32503 0.57457 1.4235 -0.002365 0.00909 -0.003286 +0.26489 0.57176 1.4025 -0.001764 0.009028 -0.003256 +0.20813 0.55327 1.3241 -0.00023901 0.009155 -0.003206 +0.23531 0.56156 1.2277 0.001402 0.008585 -0.00307 +0.32562 0.57873 1.344 -0.001189 0.008503 -0.003307 +0.32373 0.58341 1.3927 -0.002018 0.008795 -0.003355 +0.25016 0.57644 1.359 -0.0011 0.008735 -0.003322 +0.22191 0.5609 1.2768 -0.00019 0.009371 -0.003292 +0.23339 0.55583 1.2372 0.00017399 0.009068 -0.003257 +0.2986 0.56161 1.3338 -0.00052 0.009248 -0.003374 +0.3211 0.575 1.376 -0.00089 0.00872 -0.003434 +0.22937 0.56364 1.302 0.00027 0.008446 -0.003574 +0.21264 0.55873 1.2553 0.000562 0.0087 -0.003643 +0.24316 0.54964 1.2993 -0.000331 0.009039 -0.003666 +0.31702 0.56846 1.408 -0.001974 0.009256 -0.003676 +0.29303 0.57818 1.4292 -0.002834 0.009224 -0.003758 +0.22268 0.56082 1.4035 -0.0018 0.008574 -0.003635 +0.1986 0.5334 1.3304 -0.00107 0.008794 -0.003357 +0.2225 0.53419 1.3212 -0.00096201 0.009296 -0.003307 +0.33266 0.55772 1.4098 -0.002328 0.008853 -0.003418 +0.33236 0.55954 1.441 -0.002969 0.009075 -0.003455 +0.25771 0.56437 1.3925 -0.00175 0.008933 -0.00347 +0.2299 0.55824 1.354 -0.00087101 0.009173 -0.003341 +0.25582 0.55325 1.3265 -0.000419 0.009493 -0.003384 +0.33271 0.57225 1.407 -0.001749 0.009316 -0.003313 +0.33816 0.58039 1.428 -0.002343 0.008826 -0.003516 +0.2457 0.57134 1.3819 -0.000886 0.009175 -0.003428 +0.24619 0.56257 1.3224 -0.000216 0.008753 -0.003511 +0.2282 0.55067 1.2246 0.001273 0.008612 -0.003489 +0.23669 0.55776 1.2702 0.000639 0.009073 -0.003605 +0.31805 0.56734 1.4034 -0.001866 0.009303 -0.003592 +0.33706 0.59293 1.4517 -0.002943 0.009385 -0.003577 +0.26543 0.57185 1.438 -0.002861 0.009145 -0.003505 +0.17927 0.5475 1.376 -0.001433 0.008962 -0.003475 +0.19006 0.55231 1.2747 -0.000697 0.008931 -0.003474 +0.26543 0.56341 1.3576 -0.001905 0.009166 -0.003472 +0.32352 0.56756 1.4171 -0.002767 0.009223 -0.003563 +0.26235 0.56804 1.3594 -0.001344 0.008541 -0.003617 +0.19018 0.54092 1.3322 -0.000724 0.009067 -0.003483 +0.2063 0.5317 1.2628 -2.2009e-05 0.008421 -0.003388 +0.29895 0.55002 1.369 -0.001558 0.009184 -0.00349 +0.33651 0.56435 1.3977 -0.002468 0.008635 -0.003547 +0.23697 0.55612 1.3747 -0.00166 0.008982 -0.003478 +0.19674 0.54129 1.3418 -0.000575 0.008857 -0.003438 +0.23577 0.55339 1.3116 -0.000307 0.008365 -0.003398 +0.31611 0.55894 1.4016 -0.002186 0.009084 -0.003502 +0.29821 0.57612 1.4216 -0.00291 0.009187 -0.003397 +0.23061 0.55891 1.4094 -0.001858 0.008935 -0.003437 +0.20166 0.54485 1.3591 -0.00069099 0.009077 -0.003285 +0.26118 0.55224 1.3774 -0.00105 0.009052 -0.003304 +0.3342 0.56988 1.4228 -0.002266 0.009199 -0.003337 +0.2913 0.56643 1.4034 -0.001797 0.009258 -0.003335 +0.21723 0.55975 1.3521 -0.0007 0.009053 -0.003348 +0.24336 0.56291 1.306 0.000247 0.008764 -0.00327 +0.33058 0.57017 1.4052 -0.002055 0.009162 -0.003296 +0.31733 0.58022 1.4314 -0.002636 0.009174 -0.00338 +0.24656 0.57196 1.4141 -0.001984 0.009101 -0.003394 +0.20481 0.5546 1.3629 -0.000706 0.009137 -0.003257 +0.23657 0.55462 1.3474 -0.000602 0.009267 -0.003283 +0.33186 0.57622 1.4069 -0.001605 0.009059 -0.003361 +0.31856 0.57462 1.431 -0.00226 0.009163 -0.003383 +0.23118 0.56785 1.4029 -0.001527 0.009143 -0.003315 +0.23567 0.5618 1.3324 -0.00054301 0.009316 -0.003241 +0.27701 0.57834 1.3177 0.00016699 0.009043 -0.003269 +0.34807 0.59145 1.4308 -0.002124 0.009253 -0.003341 +0.32074 0.60066 1.4589 -0.002525 0.009255 -0.003384 +0.24152 0.58907 1.3906 -0.001487 0.00887 -0.003285 +0.23371 0.57488 1.3109 -0.00010701 0.009259 -0.003224 +0.28111 0.59436 1.3854 -0.001165 0.009439 -0.003298 +0.33545 0.60318 1.4534 -0.002244 0.009024 -0.003305 +0.23893 0.58681 1.3561 -0.001253 0.008932 -0.003309 +0.2182 0.57149 1.271 0.00014699 0.009421 -0.00328 +0.25908 0.57795 1.3134 -0.00032301 0.00954 -0.003363 +0.3318 0.60588 1.4359 -0.001797 0.009544 -0.003479 +0.28624 0.58117 1.336 -0.00125 0.008812 -0.003557 +0.17557 0.5455 1.1314 7.8994e-05 0.008561 -0.003522 +0.12659 0.52279 0.93832 0.00097601 0.007324 -0.003626 +0.20608 0.54048 1.1564 0.00015401 0.008652 -0.00372 +0.28785 0.55228 1.2908 -0.002083 0.009173 -0.003663 +0.29498 0.55798 1.3534 -0.003222 0.009144 -0.003497 +0.20513 0.56021 1.3348 -0.002459 0.009052 -0.003383 +0.18 0.52868 1.286 -0.001378 0.009204 -0.003358 +0.20777 0.53609 1.3088 -0.00084701 0.009328 -0.003301 +0.30394 0.55354 1.3921 -0.002211 0.008949 -0.003459 +0.31147 0.57703 1.3701 -0.001268 0.009212 -0.003417 +0.24602 0.56175 1.3394 -0.000785 0.009045 -0.003421 +0.25494 0.56357 1.3023 -0.000148 0.009007 -0.003279 +0.2744 0.56794 1.3461 -0.001132 0.009389 -0.003408 +0.33828 0.5812 1.4099 -0.001775 0.009071 -0.003387 +0.28381 0.57385 1.3407 -0.001379 0.009764 -0.003457 +0.18459 0.54612 1.1697 0.000165 0.008705 -0.003574 +0.04853 0.51985 0.85558 0.00049799 0.004376 -0.003746 +0.16057 0.51772 0.98937 -0.000777 0.006622 -0.003629 +0.25878 0.54088 1.2029 -0.002333 0.007297 -0.003598 +0.25191 0.5575 1.282 -0.002782 0.007344 -0.003773 +0.15011 0.54095 1.2353 -0.001963 0.007111 -0.003489 +0.13989 0.51831 1.1304 -0.00071099 0.007144 -0.003498 +0.18935 0.52412 1.1095 -0.000551 0.00704 -0.003284 +0.28701 0.55035 1.2377 -0.001228 0.007644 -0.003585 +0.28509 0.54543 1.2701 -0.001821 0.007945 -0.003552 +0.19694 0.53216 1.2101 -0.00065699 0.007568 -0.00341 +0.17726 0.51365 1.0857 9.1996e-05 0.007406 -0.003302 +0.22993 0.52754 1.1062 0.000129 0.007372 -0.003268 +0.27704 0.54965 1.2272 -0.001293 0.007617 -0.003381 +0.22349 0.54059 1.2605 -0.001659 0.007537 -0.003324 +0.14259 0.53376 1.217 -0.000396 0.007208 -0.003276 +0.16663 0.52349 1.1495 5.0004e-05 0.006912 -0.003235 +0.2274 0.54501 1.206 -0.000535 0.007296 -0.00328 +0.28878 0.55632 1.2843 -0.001626 0.007931 -0.003393 +0.25439 0.55556 1.2928 -0.001439 0.007693 -0.003446 +0.18483 0.54061 1.2418 -0.000398 0.006979 -0.003293 +0.1906 0.53234 1.1538 0.00052899 0.007134 -0.003273 +0.25411 0.54831 1.2511 -0.000792 0.007534 -0.003328 +0.31095 0.55715 1.3302 -0.00223 0.007981 -0.003328 +0.23116 0.55819 1.3263 -0.001673 0.007354 -0.003295 +0.15688 0.5428 1.2551 -0.00025201 0.007019 -0.003207 +0.17066 0.53555 1.1901 0.00044099 0.007418 -0.003289 +0.26422 0.55707 1.3044 -0.001067 0.007897 -0.003258 +0.28839 0.56585 1.343 -0.001893 0.008059 -0.003372 +0.22733 0.56038 1.3291 -0.001381 0.007895 -0.003213 +0.17016 0.54824 1.2496 1.2008e-05 0.007164 -0.003221 +0.20069 0.55143 1.2285 0.00024201 0.006958 -0.003239 +0.29137 0.57142 1.3159 -0.001157 0.007884 -0.003347 +0.30322 0.58127 1.3467 -0.001704 0.008109 -0.00337 +0.22881 0.57626 1.3248 -0.000869 0.007363 -0.003234 +0.20747 0.58179 1.2367 0.00059299 0.00775 -0.003279 +0.24947 0.58151 1.291 -0.000185 0.007893 -0.003314 +0.31303 0.59296 1.3618 -0.001644 0.008167 -0.003327 +0.25676 0.59169 1.3671 -0.001646 0.008185 -0.003291 +0.19247 0.57853 1.2906 0.00042101 0.007702 -0.003323 +0.20786 0.58123 1.2741 0.000426 0.007729 -0.003304 +0.26751 0.58262 1.351 -0.00099601 0.008193 -0.003372 +0.3431 0.59412 1.41 -0.002107 0.008397 -0.003354 +0.25399 0.58737 1.3738 -0.00156 0.00803 -0.003265 +0.21791 0.58732 1.2805 0.000601 0.007932 -0.003082 +0.23331 0.58621 1.2724 0.000477 0.008086 -0.003242 +0.26949 0.59676 1.3325 -0.000176 0.008169 -0.003361 +0.31991 0.60836 1.3815 -0.001216 0.007525 -0.003435 +0.27048 0.60189 1.3533 -0.00091 0.007538 -0.003432 +0.24792 0.59904 1.328 -4.0998e-05 0.008069 -0.00345 +0.26862 0.60418 1.3769 -0.00061699 0.008101 -0.003413 +0.3261 0.6186 1.4177 -0.001263 0.007735 -0.003495 +0.2464 0.61217 1.3188 0.00026 0.007395 -0.003376 +0.21707 0.62272 1.2755 0.00071099 0.007901 -0.003296 +0.21875 0.6073 1.3031 0.000358 0.007693 -0.003328 +0.26389 0.61879 1.39 -0.001646 0.008054 -0.00336 +0.29113 0.60608 1.4241 -0.00259 0.00808 -0.003356 +0.2357 0.594 1.3984 -0.002382 0.007993 -0.003354 +0.19933 0.57664 1.3382 -0.000966 0.007957 -0.003294 +0.20526 0.57096 1.2954 0.000102 0.007944 -0.003369 +0.28141 0.59266 1.3474 -0.000919 0.008015 -0.003458 +0.33191 0.60287 1.3848 -0.00222 0.008343 -0.003356 +0.25844 0.59936 1.333 -0.00032899 0.00831 -0.003314 +0.22157 0.59381 1.3079 0.00044501 0.007972 -0.003296 +0.26095 0.59427 1.2881 0.000315 0.008156 -0.003302 +0.32069 0.60904 1.3622 -0.001409 0.007882 -0.003459 +0.32958 0.61484 1.3971 -0.002041 0.007925 -0.00331 +0.24748 0.6109 1.3725 -0.001226 0.007839 -0.003355 +0.20097 0.60042 1.2938 0.000464 0.008029 -0.003403 +0.22038 0.59315 1.2621 0.000956 0.007806 -0.003413 +0.26547 0.60602 1.3589 -0.000115 0.008465 -0.003517 +0.31526 0.616 1.3835 -0.001115 0.007884 -0.003604 +0.21329 0.58844 1.1716 0.000697 0.008047 -0.003611 +0.19751 0.56871 1.0895 0.002711 0.008256 -0.003989 +0.20701 0.57259 1.0609 0.003176 0.008241 -0.004296 +0.14041 0.56089 0.95603 0.003368 0.007257 -0.00469 +0.087227 0.5196 0.8594 0.003043 0.006831 -0.0047 +0.23001 0.55221 1.0345 0.001015 0.007552 -0.004632 +0.29938 0.56281 1.2179 -0.00094 0.008266 -0.004587 +0.24614 0.56766 1.2367 -0.00074499 0.008112 -0.004389 +0.13748 0.55852 1.1875 0.00062699 0.007756 -0.004115 +0.14838 0.54653 1.1525 0.001481 0.007198 -0.004068 +0.22296 0.56099 1.2455 0.000206 0.008061 -0.004402 +0.26159 0.58409 1.2791 -0.000578 0.008248 -0.004516 +0.22949 0.57486 1.3186 -0.00067399 0.007829 -0.004297 +0.17385 0.57729 1.278 0.000758 0.007606 -0.004141 +0.23539 0.57016 1.2473 0.00071901 0.008147 -0.004031 +0.31607 0.58484 1.3202 -0.000605 0.008438 -0.004137 +0.30446 0.59103 1.3508 -0.00122 0.008498 -0.004086 +0.21155 0.57949 1.2932 -3.6006e-05 0.007513 -0.004258 +0.19955 0.56859 1.2251 0.00094201 0.008083 -0.004204 +0.25746 0.577 1.2979 -0.00015101 0.008306 -0.004142 +0.28782 0.5834 1.3589 -0.001126 0.008455 -0.004289 +0.22769 0.57924 1.3397 -0.000951 0.008338 -0.004399 +0.18094 0.56232 1.2806 0.00044 0.008281 -0.004249 +0.2039 0.55913 1.252 0.000595 0.008295 -0.004152 +0.29885 0.57476 1.3342 -0.000653 0.008639 -0.004029 +0.31455 0.58689 1.3809 -0.001402 0.008751 -0.004106 +0.23495 0.58239 1.3501 -0.000656 0.008452 -0.004197 +0.21844 0.571 1.2728 0.000568 0.00848 -0.003986 +0.24808 0.57881 1.2871 0.00028601 0.008594 -0.003878 +0.32462 0.59709 1.3795 -0.001075 0.008777 -0.003985 +0.28628 0.59774 1.376 -0.00138 0.008746 -0.003938 +0.20373 0.59415 1.3281 4.0998e-05 0.008409 -0.004073 +0.21348 0.57788 1.3075 0.000332 0.008639 -0.004036 +0.27271 0.58321 1.3684 -0.00044099 0.008824 -0.003965 +0.32039 0.59528 1.4119 -0.001465 0.008851 -0.004346 +0.2431 0.58837 1.3841 -0.001196 0.008754 -0.004369 +0.18791 0.57172 1.3014 0.00022099 0.008768 -0.004222 +0.2111 0.57263 1.2935 0.00012999 0.009042 -0.004308 +0.30412 0.59439 1.4006 -0.001297 0.009014 -0.004364 +0.30908 0.60445 1.417 -0.002042 0.008959 -0.004404 +0.22085 0.59563 1.3798 -0.001242 0.008658 -0.004288 +0.20463 0.58377 1.335 -0.00015699 0.008976 -0.004161 +0.28506 0.59424 1.3955 -0.000812 0.009071 -0.004195 +0.34063 0.60866 1.4416 -0.001958 0.009107 -0.004224 +0.27709 0.60348 1.4152 -0.001863 0.008971 -0.004226 +0.20886 0.58658 1.3509 -0.000639 0.00883 -0.004124 +0.22471 0.58383 1.3253 0.000152 0.008996 -0.004075 +0.30454 0.5994 1.4083 -0.001038 0.009113 -0.003946 +0.31926 0.61759 1.4125 -0.001542 0.008504 -0.004033 +0.2619 0.61114 1.3916 -0.001456 0.008844 -0.004047 +0.22779 0.59117 1.3627 -0.00078199 0.009058 -0.00394 +0.25598 0.58775 1.3611 -0.00085401 0.009198 -0.004014 +0.33798 0.61083 1.4271 -0.001969 0.009219 -0.004037 +0.29953 0.6117 1.4373 -0.002379 0.009071 -0.004039 +0.21927 0.5996 1.3907 -0.001265 0.00903 -0.004207 +0.22292 0.59437 1.3732 -0.000538 0.008945 -0.004143 +0.28921 0.60115 1.4338 -0.001693 0.009 -0.00389 +0.32985 0.6117 1.459 -0.002768 0.009152 -0.004075 +0.25911 0.60967 1.4368 -0.002184 0.009013 -0.003962 +0.22217 0.59566 1.3911 -0.00092799 0.009146 -0.004144 +0.24816 0.59263 1.3736 -0.000582 0.009224 -0.003899 +0.30737 0.61123 1.4323 -0.001713 0.008998 -0.003958 +0.32205 0.6211 1.445 -0.002294 0.009057 -0.004057 +0.24402 0.61193 1.3855 -0.001267 0.008438 -0.00404 +0.23395 0.59566 1.3287 0.000257 0.009362 -0.003885 +0.27547 0.61099 1.3401 -0.000274 0.009261 -0.003913 +0.32875 0.63039 1.4297 -0.001924 0.008657 -0.004052 +0.29667 0.62343 1.4112 -0.001246 0.009192 -0.004096 +0.22778 0.60917 1.3492 -0.000598 0.00864 -0.004027 +0.2201 0.5987 1.2664 -4.9916e-06 0.008946 -0.003679 +0.25967 0.61528 1.3257 -0.00049201 0.009197 -0.00429 +0.32627 0.62392 1.398 -0.001385 0.009787 -0.004375 +0.33457 0.62277 1.437 -0.00239 0.009789 -0.004462 +0.25958 0.60984 1.4059 -0.001934 0.009641 -0.004263 +0.24822 0.59412 1.3406 -0.00088399 0.009752 -0.003961 +0.27042 0.59978 1.3771 -0.001104 0.009685 -0.004064 +0.33073 0.62655 1.4502 -0.001735 0.009671 -0.004108 +0.2977 0.6311 1.4527 -0.002915 0.009283 -0.004066 +0.21069 0.60454 1.4054 -0.001233 0.009538 -0.004161 +0.22832 0.59366 1.3558 -0.000998 0.008987 -0.003971 +0.31025 0.61229 1.392 -0.00071799 0.008979 -0.004161 +0.34341 0.62177 1.4503 -0.002607 0.009319 -0.00412 +0.28937 0.61515 1.428 -0.002222 0.009338 -0.004125 +0.24773 0.60609 1.3575 -0.000947 0.009199 -0.003972 +0.25706 0.60126 1.3388 -0.000114 0.01004 -0.003943 +0.32883 0.63127 1.4203 -0.000835 0.009411 -0.004067 +0.33864 0.64038 1.4641 -0.001789 0.009537 -0.004058 +0.27584 0.62637 1.4014 -0.001104 0.00924 -0.004229 +0.21624 0.58881 1.2588 4.3005e-05 0.009577 -0.004017 +0.14656 0.55653 1.084 0.00065101 0.008304 -0.004825 +0.1575 0.53361 1.1326 -4.1993e-05 0.008825 -0.004574 +0.32685 0.54302 1.4211 -0.002342 0.010854 -0.005282 +0.27584 0.5023 1.2143 -0.001714 0.011328 -0.005914 +0.22313 0.4975 1.1505 -0.001452 0.010402 -0.005221 +0.32601 0.53519 1.3096 -0.002973 0.010963 -0.005123 +0.32484 0.55678 1.386 -0.004039 0.01044 -0.004952 +0.21453 0.55161 1.368 -0.003587 0.010302 -0.004795 +0.14688 0.53459 1.3054 -0.002246 0.010406 -0.00445 +0.14528 0.5329 1.257 -0.00093501 0.010302 -0.004299 +0.24343 0.5611 1.3608 -0.00278 0.010447 -0.00453 +0.27631 0.56518 1.3836 -0.003122 0.010166 -0.004719 +0.21408 0.57134 1.3412 -0.001872 0.009995 -0.00468 +0.22082 0.56019 1.3088 -0.001344 0.010247 -0.004721 +0.2747 0.56706 1.3527 -0.001499 0.01035 -0.00439 +0.3556 0.58668 1.3987 -0.002911 0.01022 -0.004614 +0.33641 0.59555 1.418 -0.002912 0.01001 -0.004611 +0.23186 0.58678 1.3862 -0.001807 0.009975 -0.004472 +0.20817 0.57128 1.3176 -0.000757 0.010235 -0.004478 +0.26358 0.58896 1.3836 -0.00194 0.010072 -0.004557 +0.3025 0.60457 1.4284 -0.002742 0.010074 -0.004494 +0.23983 0.59418 1.4142 -0.00269 0.01005 -0.004452 +0.20199 0.58472 1.3578 -0.00086401 0.010159 -0.004282 +0.22511 0.58424 1.3029 -0.000392 0.01016 -0.004336 +0.31848 0.60324 1.4036 -0.002186 0.01016 -0.004384 +0.33428 0.61468 1.4296 -0.002659 0.010222 -0.004465 +0.26188 0.60685 1.3908 -0.001405 0.010171 -0.004406 +0.24227 0.59698 1.314 -0.000365 0.009919 -0.003865 +0.27982 0.6038 1.3008 0.00027 0.010268 -0.004263 +0.32735 0.63679 1.4181 -0.002234 0.010406 -0.004404 +0.32312 0.6293 1.4245 -0.002012 0.010365 -0.004454 +0.24253 0.61834 1.3416 -0.000963 0.010581 -0.004483 +0.21505 0.60607 1.2523 -0.0002 0.010195 -0.004573 +0.17672 0.57481 1.1177 2.7995e-05 0.00974 -0.004392 +0.23856 0.58556 1.2682 -0.000667 0.010701 -0.004859 +0.32523 0.57713 1.3996 -0.002468 0.011049 -0.005101 +0.26232 0.53377 1.1634 -0.001415 0.009873 -0.004717 +0.069961 0.54096 0.75301 0.004181 0.007419 -0.004779 +0.11451 0.54288 0.74704 0.004148 0.007633 -0.004502 +0.26097 0.55287 1.013 0.002047 0.008992 -0.00484 +0.25531 0.58204 1.1343 0.001752 0.009243 -0.004798 +0.15344 0.56134 1.1232 0.001653 0.009029 -0.00443 +0.12026 0.5456 1.0752 0.002465 0.009006 -0.00428 +0.17446 0.56313 1.0707 0.002874 0.009064 -0.004336 +0.28385 0.57157 1.2054 0.001068 0.009664 -0.004657 +0.26964 0.57426 1.2103 0.000271 0.009638 -0.004454 +0.1942 0.55934 1.1636 0.001227 0.009176 -0.004274 +0.18195 0.55576 1.0857 0.002459 0.009273 -0.004234 +0.2694 0.57005 1.1619 0.001303 0.009588 -0.004324 +0.32449 0.58314 1.2575 0.00036299 0.009504 -0.004453 +0.25844 0.58627 1.2519 0.000331 0.00962 -0.004661 +0.1911 0.58321 1.2079 0.00135 0.009306 -0.004531 +0.21005 0.56902 1.1849 0.001577 0.009778 -0.004391 +0.27975 0.58399 1.243 0.001134 0.009976 -0.004318 +0.29037 0.59149 1.2856 0.000443 0.010042 -0.00445 +0.20577 0.58501 1.2424 0.001269 0.009672 -0.004585 +0.17896 0.57838 1.1851 0.001803 0.009184 -0.004335 +0.23358 0.57258 1.2166 0.001268 0.009793 -0.004457 +0.32183 0.59445 1.289 0.000429 0.010051 -0.004383 +0.2798 0.59147 1.2994 0.000267 0.009967 -0.004508 +0.2031 0.57996 1.2516 0.001519 0.00978 -0.004469 +0.20141 0.58331 1.1836 0.002036 0.009268 -0.004365 +0.26259 0.58681 1.2441 0.00148 0.008951 -0.004482 +0.31323 0.59662 1.3125 0.000273 0.00971 -0.004629 +0.2541 0.59155 1.3154 0.000402 0.009815 -0.004429 +0.16886 0.59808 1.2391 0.001162 0.009007 -0.004244 +0.17459 0.56677 1.1775 0.002249 0.009127 -0.004161 +0.28951 0.59243 1.2807 0.00124 0.009648 -0.004091 +0.3167 0.5995 1.3036 0.000419 0.009945 -0.004108 +0.23828 0.5932 1.2894 0.001181 0.009656 -0.004109 +0.20894 0.58654 1.2634 0.002076 0.009498 -0.003992 +0.2531 0.59019 1.2618 0.001962 0.009585 -0.003883 +0.33223 0.61046 1.3228 0.00069 0.009957 -0.004066 +0.32659 0.6099 1.3343 0.00021 0.010049 -0.003896 +0.23394 0.60653 1.2962 0.001454 0.009684 -0.004121 +0.22822 0.59044 1.2699 0.001875 0.00985 -0.003943 +0.27827 0.59632 1.2885 0.001531 0.009909 -0.003667 +0.32767 0.61079 1.3521 0.00023501 0.010012 -0.003771 +0.30307 0.60939 1.3695 -0.000111 0.010001 -0.003719 +0.21983 0.61098 1.3192 0.00107 0.009632 -0.004013 +0.24907 0.59824 1.3128 0.001443 0.00981 -0.003779 +0.30171 0.60772 1.3476 0.000605 0.009947 -0.003794 +0.35418 0.61781 1.3849 -0.00013299 0.009988 -0.003727 +0.28209 0.61588 1.3622 0.000172 0.009842 -0.003675 +0.22917 0.61342 1.3054 0.001488 0.009508 -0.003582 +0.27391 0.59513 1.2917 0.001337 0.009904 -0.003386 +0.34547 0.60544 1.3781 0.000362 0.009814 -0.003714 +0.35079 0.61581 1.3834 -0.00012701 0.009664 -0.003666 +0.28022 0.60841 1.354 0.000738 0.00945 -0.003644 +0.26285 0.60156 1.3213 0.001433 0.009728 -0.003582 +0.29669 0.60513 1.3401 0.001437 0.009751 -0.003553 +0.36463 0.62264 1.4055 8.1001e-05 0.009819 -0.003579 +0.32988 0.62286 1.3996 -7.1e-05 0.009626 -0.003555 +0.26844 0.61718 1.3724 0.00097499 0.009543 -0.003599 +0.28525 0.61285 1.3628 0.001664 0.00967 -0.003556 +0.31133 0.61627 1.3665 0.001372 0.009782 -0.003498 +0.37392 0.62447 1.4066 0.00016699 0.009702 -0.0036 +0.33163 0.62244 1.404 -0.00035599 0.009612 -0.003661 +0.29212 0.62232 1.3671 0.001029 0.009432 -0.003515 +0.27666 0.60829 1.3662 0.001332 0.00948 -0.003533 +0.30628 0.60147 1.3646 0.00082701 0.009574 -0.003422 +0.37626 0.62003 1.4184 0.00015401 0.009733 -0.003351 +0.37851 0.62713 1.4418 -0.000297 0.009691 -0.003262 +0.31629 0.63278 1.4176 -0.00011299 0.009062 -0.003202 +0.28836 0.59923 1.3599 0.00079001 0.00958 -0.003181 +0.34187 0.60837 1.385 0.000342 0.009752 -0.003178 +0.41177 0.6284 1.441 -0.000197 0.00963 -0.003231 +0.33514 0.60765 1.4228 -0.000264 0.009123 -0.003148 +0.2844 0.60205 1.3641 0.000856 0.009465 -0.003114 +0.33699 0.60537 1.3741 0.000521 0.009645 -0.003023 +0.42059 0.62104 1.45 -0.00047499 0.009602 -0.003036 +0.36527 0.62621 1.4487 0 0.009227 -0.003157 +0.281 0.62067 1.4314 0 0.008886 -0.003045 +0.30957 0.60423 1.3908 0.001005 0.009251 -0.003 +0.34716 0.60805 1.4239 0 0.009217 -0.002934 +0.42369 0.62676 1.4846 -0.00018 0.008678 -0.003025 +0.38676 0.61662 1.4811 -0.00023901 0.008897 -0.002366 +0.31084 0.60334 1.4232 0 0.008545 -0.002898 +0.32581 0.59529 1.3969 0 0.009327 -0.002853 diff --git a/CsfMask_07_121x145x121.nii b/CsfMask_07_121x145x121.nii new file mode 100644 index 0000000..b4c36b7 Binary files /dev/null and b/CsfMask_07_121x145x121.nii differ diff --git a/Script_Preprocessing_HCP.m b/Script_Preprocessing_HCP.m index 3325519..0b5149c 100644 --- a/Script_Preprocessing_HCP.m +++ b/Script_Preprocessing_HCP.m @@ -1,474 +1,472 @@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Authors: Thomas Bolton and Mert Inan % Date: 03/08/2017 - % Description: % This script will perform classical preprocessing steps on the HCP data in % order to analyse it in different ways: % % - Through regional sparse coupled hidden Markov models % - Through co-activation pattern analysis % - Through variance-related assessment of motion % - Through Total activation on graphs %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 1. Setting up path files % Adds the shared folder to do all we want: contains SPM8 PathSPM = '/media/miplab-nas3/Code-NAS1/shared/static_FC_pipeline/spm8'; PathPrepro = '/media/miplab-nas3/Code-NAS1/projects/ABIDE/preprocessing_Thomas/MyPreprocessing'; PathDPARSFA = '/media/miplab-nas3/Code-NAS1/shared/static_FC_pipeline/DPARSF_V2.3_130615'; % If the 'server link' exists, we add the path. Else, we add the computer % path (with Volumes) try addpath(genpath(PathSPM)); catch addpath(genpath('/Volumes/Code-NAS1/shared/static_FC_pipeline/spm8')); end try addpath(genpath(PathPrepro)); catch addpath(genpath('/Volumes/Code-NAS1/projects/ABIDE/preprocessing_Thomas/MyPreprocessing')); end try addpath(genpath(PathDPARSFA)); catch addpath(genpath('/Volumes/Code-NAS1/shared/static_FC_pipeline/DPARSF_V2.3_130615')); end % ------------- DATA PATH TO READ FROM ------------- % **** Changed from /media/miplab-nas2/data/HCP_Data_TA' **** % Need to check that writing to this path is not done! readPath = '/media/miplab-nas3/HCP-Data/HCP_Giulia'; % ------------- DATA PATH TO WRITE TO ------------- writePath = '/media/miplab-nas3/Data/Regional_SCHMM'; %% 2. Setting up options % We change directory towards where we store the RS data cd(readPath); % Choose the wave of the subjects wave = '1st_wave_90subjs'; %Change loop cycle count when you change this %wave ='2nd_wave_158subjs'; % Indices of the subjects to preprocess pathToSubjectsFunc = fullfile(readPath, 'Preproc_sFC', wave); pathToSubjectsStruct = fullfile(readPath,'Preproc_sFC',wave); d = dir(pathToSubjectsFunc); d = d(4:end,:); % Type of acquisition for the subjects (RL or LR): for this data, we only % have the LR session available for now d2 = 'rfMRI_REST1_LR'; %d2 = 'rfMRI_REST1_RL'; % Number of scans to remove removeFirstNscans = 10; %% 3. Preprocessing itself % I iterate through all the HCP folders and create respective folders in % the write directory for i = 1:158 try disp(['Subject ',d(i).name,'...']); %Creating a folder in the write directory if isequal(d2,'rfMRI_REST1_RL') writeFolder = [d(i).name,'_RL']; %Change according to d2 else writeFolder = d(i).name; end cd(writePath); mkdir(writeFolder); cd(writeFolder); - mkdir('Structural'); + %mkdir('Structural'); mkdir('Functional'); + mkdir('CAPs'); + mkdir('VAR'); + mkdir('SCHMM'); cd(readPath); % READ Structural and functional paths - structPath = fullfile(pathToSubjectsStruct,d(i).name,'T1'); + %structPath = fullfile(pathToSubjectsStruct,d(i).name,'T1'); functPath = fullfile(pathToSubjectsFunc,d(i).name,'fMRI'); % WRITE Structural and functional paths - structWritePath = fullfile(writePath,writeFolder,'Structural'); + %structWritePath = fullfile(writePath,writeFolder,'Structural'); functWritePath = fullfile(writePath,writeFolder,'Functional'); % Reads and modifies the text file to have 594 rows and to have proper % units, then copies it from read location to write location. movText = textread(fullfile(functPath,'rpMovement_Regressors.txt')); %Changes degrees to radians - movText = movText(1:594,1:6); + movText = movText(:,1:6); for l = 4:6 movText(:,l) = movText(:,l)*pi/180; end dlmwrite(fullfile(functWritePath,'Movement_Regressors.txt'),movText,' '); %% Fine-Tuning WM and CSF Masks for Each Volume % Gets one volume only [fMeanFN_noPath, foosubdirs] = spm_select('List',functPath,['^' 'frf' '.*\.nii$']); % Gets the path and the volume itself - fMean_fn=fullfile(functWritePath,fMeanFN_noPath(1,:)); + fMean_fn=fullfile(functPath,fMeanFN_noPath(1,:)); % Transforms the WM and CSF masks from the toolbox into native % resolution Masks={'WhiteMask_09_121x145x121.nii','CsfMask_07_121x145x121.nii'}; % I must convert the masks to the right resolution for j=1:length(Masks) tmp=which(Masks{j}); if j == 1 WMmaskLidx = mapVolumeToVolume_TB(tmp,fMean_fn); else CSFmaskLidx = mapVolumeToVolume_TB(tmp,fMean_fn); end end %% 5. EXTRACT FMRI TIMECOURSES disp(['Subject ',d(i).name,'...']); disp('Extracting time courses...'); % Names of all the functional files fVolsFNlist_VX=getImageFNinAcqOrder_TB(functPath,'s5frf','nii',[]); fVolsFNlist_fullpath_VX=cellfun(@(x) fullfile(functPath,x), fVolsFNlist_VX,'UniformOutput',false); % Same for unsmoothed data (needed because we want to average WM and % CSF signals on unsmoothed data) fVolsFNlist_VX_UNSMOOTHED=getImageFNinAcqOrder_TB(functPath,'frf','nii',[]); fVolsFNlist_fullpath_VX_UNSMOOTHED=cellfun(@(x) fullfile(functPath,x), fVolsFNlist_VX_UNSMOOTHED,'UniformOutput',false); % read all headers and files in temporal order into a 4-D array V0i_VX=spm_vol(fVolsFNlist_fullpath_VX); if removeFirstNscans~=0 V0i_VX=V0i_VX((removeFirstNscans+1):end); end % Same for unsmoothed V0i_VX_UNSMOOTHED=spm_vol(fVolsFNlist_fullpath_VX_UNSMOOTHED); if removeFirstNscans~=0 V0i_VX_UNSMOOTHED=V0i_VX_UNSMOOTHED((removeFirstNscans+1):end); end % Indices of the frames of interest (common to smoothed and unsmoothed % cases V0idx=1:length(V0i_VX_UNSMOOTHED); % Creates the 4D volumes V0_VX=zeros(V0i_VX{1}.dim(1),V0i_VX{1}.dim(2),V0i_VX{1}.dim(3),length(V0i_VX),'single'); for iter=1:length(V0i_VX) V0_VX(:,:,:,V0idx(iter))=spm_read_vols(V0i_VX{iter}); end V0_VX_UNSMOOTHED=zeros(V0i_VX_UNSMOOTHED{1}.dim(1),V0i_VX_UNSMOOTHED{1}.dim(2),V0i_VX_UNSMOOTHED{1}.dim(3),length(V0i_VX_UNSMOOTHED),'single'); for iter=1:length(V0i_VX_UNSMOOTHED) V0_VX_UNSMOOTHED(:,:,:,V0idx(iter))=spm_read_vols(V0i_VX_UNSMOOTHED{iter}); end %% 6. DETRENDING disp(['Subject ',d(i).name,'...']); disp('Detrending...'); AllVolume_VX=reshape(V0_VX,[],size(V0_VX,4))'; % time x space AllVolume_VX_UNSMOOTHED=reshape(V0_VX_UNSMOOTHED,[],size(V0_VX_UNSMOOTHED,4))'; % time x space CUTNUMBER=10; % cut space in parts - SegmentLength = ceil(size(AllVolume_VX_UNSMOOTHED,2) / CUTNUMBER); + SegmentLength = ceil(size(AllVolume_VX,2) / CUTNUMBER); for iCut=1:CUTNUMBER if iCut~=CUTNUMBER Segment = (iCut-1)*SegmentLength+1 : iCut*SegmentLength; else Segment = (iCut-1)*SegmentLength+1 : size(AllVolume_VX_UNSMOOTHED,2); end AllVolume_VX(:,Segment) = detrend(AllVolume_VX(:,Segment)); AllVolume_VX_UNSMOOTHED(:,Segment) = detrend(AllVolume_VX_UNSMOOTHED(:,Segment)); end V0_VX=reshape(AllVolume_VX',size(V0_VX,1),size(V0_VX,2),size(V0_VX,3),size(V0_VX,4)); V0_VX_UNSMOOTHED=reshape(AllVolume_VX_UNSMOOTHED',size(V0_VX_UNSMOOTHED,1),size(V0_VX_UNSMOOTHED,2),size(V0_VX_UNSMOOTHED,3),size(V0_VX_UNSMOOTHED,4)); %% 7. COMPUTE AVG CSF TIMECOURSE disp(' Computing average CSF signal...'); % Returns the locations in X, Y and Z of the elements from the mask [Ci,Cj,Ck]=ind2sub(size(V0_VX_UNSMOOTHED(:,:,:,1)),find(CSFmaskLidx)); % How many voxels nCSFvoxels=numel(Ci); WriteInformation_TB(fid,['Voxels CSF: ',num2str(nCSFvoxels),'\t']); % Creates the time course itself tcCSF=zeros(1,1,1,size(V0_VX_UNSMOOTHED,4)); for idx=1:numel(Ci) tcCSF=tcCSF+V0_VX_UNSMOOTHED(Ci(idx),Cj(idx),Ck(idx),:); end CSFavg=squeeze(tcCSF/nCSFvoxels)'; CSFavg_demean=CSFavg-mean(CSFavg); clear tcCSF; %% 8. COMPUTE AVG WM TIMECOURSE disp(' Computing average WM signal...'); [Ci,Cj,Ck]=ind2sub(size(V0_VX_UNSMOOTHED(:,:,:,1)),find(WMmaskLidx)); nWMvoxels=numel(Ci); % also, sum(CSFmaskLidx) WriteInformation_TB(fid,['WM: ',num2str(nWMvoxels),'\t']); tcWM=zeros(1,1,1,size(V0_VX_UNSMOOTHED,4)); for idx=1:numel(Ci) tcWM=tcWM+V0_VX_UNSMOOTHED(Ci(idx),Cj(idx),Ck(idx),:); end WMavg=squeeze(tcWM/nWMvoxels)'; WMavg_demean=WMavg-mean(WMavg); clear tcWM; %% 9. Construction of the matrix of regressors fprintf('Regressing out nuisance covariates...'); % Constant, linear and quadratic regressors X=[ones(numel(CSFavg),1) [1:numel(CSFavg)]'/numel(CSFavg) ... [1:numel(CSFavg)].^2'/(numel(CSFavg)^2)]; % Adds the WM regressor X=[X,WMavg_demean']; % Adds the CSF regressor X=[X,CSFavg_demean']; %% 10. Regression itself disp('Regressing out covariates...'); % Regression for both the data to atlas and the voxelwise data V0Cov_VX_UNSMOOTHED =zeros(size(V0_VX_UNSMOOTHED)); V0Cov_VX =zeros(size(V0_VX_UNSMOOTHED)); for I=1:size(V0_VX_UNSMOOTHED,1) fprintf('.'); for j=1:size(V0_VX_UNSMOOTHED,2) for k=1:size(V0_VX_UNSMOOTHED,3) % We keep the residual from the regression process (i.e. % what is left after we remove the impact from the % covariates) [~,res_VX] = y_regress_ss(squeeze(V0_VX(I,j,k,V0idx)),X); [~,res_VX_UNSMOOTHED] = y_regress_ss(squeeze(V0_VX_UNSMOOTHED(I,j,k,V0idx)),X); % The volumes following regression of covariates are created V0Cov_VX(I,j,k,V0idx)=res_VX+mean(squeeze(V0_VX(I,j,k,V0idx))); V0Cov_VX_UNSMOOTHED(I,j,k,V0idx)=res_VX_UNSMOOTHED+mean(squeeze(V0_VX_UNSMOOTHED(I,j,k,V0idx))); end end end V0Cov_VX(isnan(V0Cov_VX))=0; V0Cov_VX_UNSMOOTHED(isnan(V0Cov_VX_UNSMOOTHED))=0; %% 12. ATLASING (on unsmoothed data) % 1. Get the proper atlas file(s): read content of atlas NIFTI file % (spm_vol and spm_read_vols) and then convert it into the functional data % resolution (mapVolumeToVolume) if isequal(d2,'rfMRI_REST1_RL') chosenDataFile = 'frfMRI_REST1_RL001.nii'; else chosenDataFile = 'frfMRI_REST1_LR001.nii'; end % Path of the atlas is going to be used by the mapVolumeToVolume % function. - atlasDir = '/media/miplab-nas2/Data/HCP_Data_TA/HCP_for_SCHMM/Atlases'; + atlasDir = '/Users/TiBiUan/Desktop/Github/Regional_SCHMM/Atlases'; - for m = 1:2 + n_atlas = 4; + + atlas_names = {'Craddock_349.nii','Craddock_950.nii','AAL_90.nii','HCP-MMP1_MNI152_2mm_360_fnirt2diff.nii'}; + + AllVolume_V0Cov_VX = reshape(V0Cov_VX,[],size(V0Cov_VX,4)); + + AllVolume_V0Cov_VX_UNSMOOTHED = reshape(V0Cov_VX_UNSMOOTHED,[],size(V0Cov_VX_UNSMOOTHED,4)); + + for m = 1:n_atlas %Selection of the number of regions in the atlas. - if m == 1 - regionCnt = '300'; - elseif m==2 - regionCnt = '950'; + switch m + case 1 + regionCnt = 349; + case 2 + regionCnt = 950; + case 3 + regionCnt = 90; + case 4 + regionCnt = 360; end %Feeding into the mapVolumeToVolume function with a chosen 3D volume - mappedAtlas = mapVolumeToVolume(fullfile(atlasDir,['Craddock_',regionCnt,'.nii']),fullfile(functWritePath,chosenDataFile)); + mappedAtlas = mapVolumeToVolume(fullfile(atlasDir,atlas_names{m}),fullfile(functPath,chosenDataFile)); % 2. Atlasing itself with reshaped data and atlas AllVolume_mappedAtlas = reshape(mappedAtlas,[],1); - AllVolume_V0Cov_VX = reshape(V0Cov_VX,[],size(V0Cov_VX,4)); - [atlasedData, n_nc, n_TC] = Make_Atlasing(AllVolume_V0Cov_VX, AllVolume_mappedAtlas); - - %% 13. SCRUBBING - % Converts the rotational components into [mm] - disp('Scrubbing...'); - Mot = movText; - Mot(:,4:6) = 50*Mot(:,4:6); - - % Discard the first n scans if necessary - if removeFirstNscans~=0 - Mot = Mot((removeFirstNscans+1):end,:); - end + [atlasedData{m}, n_nc, n_TC] = Make_Atlasing(AllVolume_V0Cov_VX_UNSMOOTHED, AllVolume_mappedAtlas); + end + + + %% 13. SCRUBBING (of voxelwise and of atlased time courses) + + % Converts the rotational components into [mm] + disp('Scrubbing...'); + Mot = movText; + Mot(:,4:6) = 50*Mot(:,4:6); - % Computes FD - FD = sum(abs([0 0 0 0 0 0; diff(Mot)]),2); + % Discard the first n scans if necessary + if removeFirstNscans~=0 + Mot = Mot((removeFirstNscans+1):end,:); + end - % Locates the time points to scrub (set to 1) - Scrub_mask1 = logical(FD > 0.5); + % Computes FD + FD = sum(abs([0 0 0 0 0 0; diff(Mot)]),2); - % Scrubbing is done for 6 frames after - disp('Scrubbing Masks...'); - Scrub_mask2 = circshift(Scrub_mask1,[1,0]); - Scrub_mask2(1) = 0; + % Locates the time points to scrub (set to 1) + Scrub_mask1 = logical(FD > 0.5); - Scrub_mask3 = circshift(Scrub_mask1,[2,0]); - Scrub_mask3(1) = 0; + % Scrubbing is done for 6 frames after + disp('Scrubbing Masks...'); + Scrub_mask2 = circshift(Scrub_mask1,[1,0]); + Scrub_mask2(1) = 0; - Scrub_mask4 = circshift(Scrub_mask1,[3,0]); - Scrub_mask4(1) = 0; + Scrub_mask3 = circshift(Scrub_mask1,[2,0]); + Scrub_mask3(1) = 0; - Scrub_mask5 = circshift(Scrub_mask1,[4,0]); - Scrub_mask5(1) = 0; + Scrub_mask4 = circshift(Scrub_mask1,[3,0]); + Scrub_mask4(1) = 0; - Scrub_mask6 = circshift(Scrub_mask1,[5,0]); - Scrub_mask6(1) = 0; + Scrub_mask5 = circshift(Scrub_mask1,[4,0]); + Scrub_mask5(1) = 0; - Scrub_mask7 = circshift(Scrub_mask1,[6,0]); - Scrub_mask7(1) = 0; + Scrub_mask6 = circshift(Scrub_mask1,[5,0]); + Scrub_mask6(1) = 0; - %Scrubbing done for 3 frames before. - Scrub_mask8 = circshift(Scrub_mask1,[-1,0]); - Scrub_mask8(end) = 0; + Scrub_mask7 = circshift(Scrub_mask1,[6,0]); + Scrub_mask7(1) = 0; - Scrub_mask9 = circshift(Scrub_mask1,[-2,0]); - Scrub_mask9(end) = 0; + %Scrubbing done for 3 frames before. + Scrub_mask8 = circshift(Scrub_mask1,[-1,0]); + Scrub_mask8(end) = 0; - Scrub_mask10 = circshift(Scrub_mask1,[-3,0]); - Scrub_mask10(end) = 0; + Scrub_mask9 = circshift(Scrub_mask1,[-2,0]); + Scrub_mask9(end) = 0; - Scrub_mask = Scrub_mask1 | Scrub_mask2 | Scrub_mask3 | Scrub_mask4 | Scrub_mask5 | Scrub_mask6 | Scrub_mask7 | Scrub_mask8 | Scrub_mask9 | Scrub_mask10; + Scrub_mask10 = circshift(Scrub_mask1,[-3,0]); + Scrub_mask10(end) = 0; + + Scrub_mask = Scrub_mask1 | Scrub_mask2 | Scrub_mask3 | Scrub_mask4 | Scrub_mask5 | Scrub_mask6 | Scrub_mask7 | Scrub_mask8 | Scrub_mask9 | Scrub_mask10; + + % SCRUBBING ATLASED TIME COURSES + + for m = 1:n_atlas - scrubPrcntg = sum(Scrub_mask)/length(Scrub_mask)*100; - save('Scrubbed_Percentage','scrubPrcntg'); - % Finds the time points that we know - TCon = atlasedData(:,find(~Scrub_mask)); - TCon(isnan(TCon)) = 0; + TCon{m} = atlasedData{m}(:,find(~Scrub_mask)); + TCon{m}(isnan(TCon{m})) = 0; % tinter is all the time points tinter = 1:length(Scrub_mask); % torig is the time points that we know torig = tinter(find(~Scrub_mask)); % We interpolate all time point values (tinter) using the info that we % know (torig,TCon) % atlasedData has size n_vox x n_TP disp('Before interpolation...'); - TC2 = interp1(torig,TCon',tinter,'spline')'; + TC2{m} = interp1(torig,TCon{m}',tinter,'spline')'; % We want to keep the 'NaN's that we had because they are used to % remove bad regional time courses later on - TC2(isnan(atlasedData)) = NaN; + TC2{m}(isnan(atlasedData{m})) = NaN; disp('Reviving all the NaNs...'); + end + + + % Finds the time points that we know + TCon_VX = AllVolume_V0Cov_VX(:,find(~Scrub_mask)); + TCon_VX(isnan(TCon_VX)) = 0; + % We interpolate all time point values (tinter) using the info that we + % know (torig,TCon) + % atlasedData has size n_vox x n_TP + disp('Before interpolation...'); + TC2_VX = interp1(torig,TCon_VX',tinter,'spline')'; + + % We want to keep the 'NaN's that we had because they are used to + % remove bad regional time courses later on + TC2_VX(isnan(TCon_VX)) = NaN; + disp('Reviving all the NaNs...'); + + + + %% 14. FILTERING + + for m = 1:n_atlas + + + Atlas_data{m} = y_IdealFilter(TC2{m}', 0.72, [0.01, 0]); - %% 14. FILTERING - disp(['Subject ',d(i).name,'...']); - disp('Filtering...'); - % In this step, two different filtering strategies are used on the - % atlased data: - % 1) 0.01 Hz - 0.1 Hz Bandpass Filtering (B) - % 2) 0.01 Hz Highpass filtering (H) - - B_Filt_Atl_Data = y_IdealFilter(TC2', 0.72, [0.01, 0.1]); - H_Filt_Atl_Data = y_IdealFilter(TC2', 0.72, [0.01, 0]); - - %Following code generates the same thing, but using a for loop: -% B2_Filt_Atl_Data = zeros(size(TC2,1),size(TC2,2)); -% for m=1:size(TC2,1) -% B2_Filt_Atl_Data(m,:)=y_IdealFilter(TC2(m,:)',0.72,[0.01, 0.1]); -% end -% -% H2_Filt_Atl_Data = zeros(size(TC2,1),size(TC2,2)); -% for m=1:size(TC2,1) -% H2_Filt_Atl_Data(m,:)=y_IdealFilter(TC2(m,:)',0.72,[0.01, 0]); -% end - - %% 15. Saving All Files - disp(['Subject ',d(i).name,'...']); - disp(['Saving All Files for ', regionCnt, '...']); - cd(fullfile(writePath,writeFolder)); - save([regionCnt,'B_Filt_Atl_Data',int2str(i),'.mat'], 'B_Filt_Atl_Data'); -% save([int2str(regionCnt),'B2_Filt_Atl_Data',int2str(i),'.mat'], 'B2_Filt_Atl_Data'); - save([regionCnt,'H_Filt_Atl_Data',int2str(i),'.mat'], 'H_Filt_Atl_Data'); -% save([int2str(regionCnt),'H2_Filt_Atl_Data',int2str(i),'.mat'], 'H2_Filt_Atl_Data'); - - colormap('jet'); - caxis([-1 1]); - - %Plotting the functional connectivity matrix - %Pearson correlation matrix is found by the following: - B_pCorr = corr(real(B_Filt_Atl_Data)); - imagesc(B_pCorr); - print([regionCnt,'B_Filt_Atl_Data_plot'],'-dpng'); -% -% B2_pCorr = corr(real(B2_Filt_Atl_Data)'); -% imagesc(B2_pCorr); -% print([int2str(regionCnt),'B2_Filt_Atl_Data_plot'],'-dpng'); - - H_pCorr = corr(real(H_Filt_Atl_Data)); - imagesc(H_pCorr); - print([regionCnt,'H_Filt_Atl_Data_plot'],'-dpng'); - -% H2_pCorr = corr(real(H2_Filt_Atl_Data)'); -% imagesc(H2_pCorr); -% print([int2str(regionCnt),'H2_Filt_Atl_Data_plot'],'-dpng'); - - %Some statistical work on the Pearson Correlation - %First vectorize the pCorr - %vpCorr = reshape(pCorr,[],1); - %rng_vpCorr = range(pCorr); - %min_vpCorr = min(pCorr); end + H_Filt_VX_Data = y_IdealFilter(TC2_VX', 0.72, [0.01, 0]); + + + % Z-scoring the CAP data + CAP_Data = zscore(H_Filt_VX_Data); + + VAR_Data = var(H_Filt_VX_Data); + + + %% 15. Saving All Files + + save(fullfile(functWritePath,'CAP','CAP_Data'),'CAP_Data','-v7.3'); + save(fullfile(functWritePath,'VAR','VAR_Data'),'VAR_Data','-v7.3'); + save(fullfile(functWritePath,'SCHMM','Atlas_Data'),'Atlas_Data','-v7.3'); + disp(['Finished for Subject ',d(i).name,'...']); catch warning(['Could not do it for the ',int2str(i), '. iteration.']); continue; end -end - -fclose(fid); - -% ** End of Script_Preprocessing_HCP_Mert.m ** \ No newline at end of file +end \ No newline at end of file diff --git a/WhiteMask_09_121x145x121.nii b/WhiteMask_09_121x145x121.nii new file mode 100644 index 0000000..4541eed Binary files /dev/null and b/WhiteMask_09_121x145x121.nii differ