From ee16bcd4245f521c938990af29be1e4915bdbe5a Mon Sep 17 00:00:00 2001
From: Challen <273239067@qq.com>
Date: Thu, 15 Jun 2017 13:55:49 +0800
Subject: [PATCH] New ROI trace back based on M2 segmention

---
 M2lookuptable.mat                           |  Bin 0 -> 11193 bytes
 ROI_volmue_283_MNI.m                        |    9 +-
 level5_to_level3.m => level5_to_level3_M1.m |    0
 level5_to_level3_M2.m                       |  105 ++
 level5_to_level4.m                          | 1252 ------------------
 level5_to_level4_M1.m                       | 1287 ++++++++++++++++++
 level5tolevel4_manual_select.m              | 1288 +++++++++++++++++++
 look_up_table_M2_level3.mat                 |  Bin 0 -> 3019 bytes
 8 files changed, 2685 insertions(+), 1256 deletions(-)
 create mode 100644 M2lookuptable.mat
 rename level5_to_level3.m => level5_to_level3_M1.m (100%)
 create mode 100644 level5_to_level3_M2.m
 delete mode 100644 level5_to_level4.m
 create mode 100644 level5_to_level4_M1.m
 create mode 100644 level5tolevel4_manual_select.m
 create mode 100644 look_up_table_M2_level3.mat

diff --git a/M2lookuptable.mat b/M2lookuptable.mat
new file mode 100644
index 0000000000000000000000000000000000000000..36082a5dc422ae9c75ed9940a1a07c9422d5d7d8
GIT binary patch
literal 11193
zcma)?cT`hbx3@tA6psX?Nec)nA{~?(K#HI!ML_AGf+EsM=#bD6QMxpNfGEWRf=GwZ
zB=mqZkrH}BFG&cbUXJH|-+S+O|M~V<W9^KLz2};1{^s6mt|u|pvAnCJD|%g4Ues8}
zQs$|f#}jE$Gml4*r`~>^>Y`?PR<}*0MfLoiJc2xN5%u;`7qvP*8~J;QDqIs)P*Yb>
zR##LORg_my5&i!@82)`Q-Zlmte^!uUU<k^~>tC~a#mLP%ANX{DyOTrLK=NIWVuq%I
zTUC1X#hcG0?!4r+3%Ex?I$Gq9?zP8b3UC_48FYi`)l*S@9@4jBG_>zt%!VI!<g+Kb
zR3x9AT`-C5jj^GaQg9+RX$Y>t+u(Z;1xWw!kC98mD@FZ=bwaajklPTckuSPgx^)7x
z0uW~G#5=oqFU>|>h~TiNe%)yY&c;FPPi$?Wlz~;0msF!U<O#O6XgaAvXqFwKF>=4~
zlYZTqSue<a>_OplqL)Y`6|0D4FVuKZA>eSfQ4h;rtP!<+rO_FiRICwM!8>~!;*T{c
z(s(VKSOJ{nhn&R<j9f2rHwfgNRfAwg9Bu@1&R&9KjX3B8@;IDnbjR9`d@pKB+Gc5d
zi}fmUHw+YXU~SaEs*m^=H6?BX8@(XJ5zQM!4hKM^BP1M48aZ1isIL;U&Dn^=#*Jv|
z5V>cSAgGbog`tsN+_UE!v#^duf`%%QUgu`H8z-@Ig`tVt!i_|1$w>GOBIm3G1U?e3
zL*$uNfuKhQ3(1k&5{)(xFcw-wriu_*85rpF&m)s289L`#=0$!vO!06^|MKe41Khs5
zd-7bhV}xNr4)%N5i^5;+ab?`0R$r|)&u%`xdGFpmd_1Llu=`oq!!$Xi^c>7!R!$Jd
z$)iUrE@Q(zXRQ)D8{|LdJi7EQGyw~u@ztH2H9rn$X3-~uI-%nz+Crp;zvx@Vd+{`m
zV<!sj61{{Q30URAPll5bUgC|$SmomBzqK)dOb=f_*2asqF&~>d^1aag=46}~U*qh5
zYxB1<+Z>IZ*rY;@m)qwXm#`*<?uL_*+ZP)xuqMSC(c99EVC+zFQ^dAd;~lJQaZ^+U
z*X(`B<gkNbAk(Zdq;7cDAdq#|0YV&}HJpsC5S(R#07f)*iRWjzAf_XldPJdPE%>mG
z??Pj@uQ#e-=iY7HB(l$j9;>5ONPe-+)d>A}kQo}iv8}~|`a~vs=?pzVUp1`X$VyRS
zQDX6(LcwTR5r=aQY_rzKIoVw#rLO`!PQ}L%u3=++mD95>kW1K-!nZHI1RB?{PxZYX
z#3fSOe<<1CtAX1oyE+ETE5eh)*0jSE`z6GkY?MRm;~#SS{iLU+_x0Vrsg}qqn@#^z
z#EE$i!}O3JxKL51$Isxgwl*VPZ{@i^m}mKqPk4o|yKed>Uz*cR@rk-yqiGOR_*ZXa
zPZxknp1N`FaS?M^S<LT5$WYbJ_vG-kF5f&>!i)w16@R(ao!06;2lITaj+?kkg_@(@
zhFqq%YYQAAS?Ef6ss3h&_~2i$uFz=A+xB8-r#<d?%Zq79F)7ndf!E3354{F)Q}dx(
zdkMw?uY#{P`UY879<`hY$(=>rbMb9$GWqo^>`Mr(zFkG(5F#{?SpejUaXkoPTOZA&
z5d{y)Qk0g?jOt5CQ%Ne)r>>#d0LvQOg{|01^h*3a>t4N}xOvGkQ#q2}<ZJhaIeV9u
z(_-x*cgp6Z<k3=_l+XCd_^47sXiPjmNEui3Q0q!IK;2vmd+TVoJ?v0w0MMaUQn~kn
zR{4YN=VXvR>~htSJ7^hsXILaNTZ@?>S+)~RqWz?Seo90idoX{ZA1C*02vH|kFJ}OS
zNWa{_^0+nNv>JwgMQvre1TDF8IYCXXZz4O7^JtiMdvEKfeSx70J7D9-aRf&8YB{XR
z1tR0~gFaJQzbAR@g;N*X+ctio?u;LflTWqTp(SRiWcCKr2Zsg3-eeX>ByQ8^O0K_K
zvg^c>MAeKKu5!06Ts_0kgDTm)6}%7|K>Om;LTTD4XdnsVg5>Vt$&)$p1Wd4E(A8w9
zHusjo>yA#+5&hG#C@R_|Vp3h#9vmo3={IuvW={GGq*hu{8&?}r65d|)L)->mak2tM
zz`jD7ZZXiy_NbKNepoSR`g9vYu-RYp)<>gUKUpGBbA@yDwDb-<6ROLgT{KlYb%79?
z?b=41l-0h>IYD;=jltH$B7(~3CyqlWa!#5q03g8TDbxu{{I<BjuPyCZrdF8K5TkpT
zDVpRm?|*mn%QET2Y_nnTZ)&D7)xM*U<PNH(D^O330d%MwW5hP-kl--Nf1$LH#yysN
zgZeB2ih<Io@p{n0y*7I1Oq+xrp3e9u)pP8;f6EK~h&HNooJguLC382mZ2At^iP}nw
z#kY&Ixjneac9LK5I~p2rj?J21Rm(v4qZ?}H7qU*<4O#A}xd**87mgNJT2-xE#8eVo
z`4(rt(6^ht0@u<0+~J}Gtuj^$KNfXB(JmMc1mGPor9R4I3i-4GDe4<f#VhfjuFuie
z_}_~}(eym&lP?`>XG;MK)_8_drQ)Li_`s?i_7!FtY-l!To;(mJ|1<?L=1rYm<&Aq9
zd?w#_?m~9psa0W;?5Z5TBshIQ+Pg>!$8Lh1J2t*0>g0W{V3iVDlccDE>CH1-s~%`8
zM?OmLKd48W_@J-YPnC>3PJL?A1((WTPD=JM_+%`Q9L(+g*|*U5C_Seiw{p1;PX?F!
z$Vv@?U_n~etfi`xY?#JHeZd0W#N79t4fj>+cwD)WYgfM*mUK4UQ&r(X+R-GfR{OPB
z!4n0Q&T4g2d`R%q1HY#Mescx!mxWNcr?B)7=WX*Yv72E`N;oYg?r55qnwCfn!(tEP
zw==|O9~{A&_L8JeU$lIB4{k)N-Gz_VI2jFHMWyj#e)<l1p>vk+v>z48td_b7eV`$|
zqI`yAlNU}8X~c?87c4J-ib2NKzaJ&g8ml_otDqrlmBfxEaIi_=oXw(q{@$!q9$CgP
zFKuLhPgFh>L?z7PTTp?9D<_cn0c%qo*zXR+o*I&28vVf$NAX^VrEoo?yqN>hBcf^c
zXv7|kt&J=<j6&QyzVzi0v7IpKT}6}SZt+Mwl5zpuXqf=Vc2l|^(wsM;Unu=XMHHdM
zuX;)MS{PP!+CgO>{%#<7fcL8{)4Ddw9Zf&2*KJgV4vHl9*d72+ya78gq|Jyyv+l6M
z3Sq$J`}pcX93jpyCjx)EeB)$Z=c^60zd~}cbpK>%c$(a?|2u|(Y&#W;nwR@^*YY(T
z-(v0=qf&u>P>y^ac60i<ue6Izq)`{4iwnJ=-s;Ci2vZZku$JR@udvUM(7DD{45bp?
zN<UpM@t&qsMt|%WZ+>oIu_0L4b|c5FW&gMwqtiESrk^j#a}4VK-kLMgexzNLOs+4t
z2|cT<Y=r%EPiT}a?DWf)J9MaK!20Em&rgK)1O=k6Sg2wKUcFCMY?Rtx4(<LRO&h{C
zRUnf?>?f8!0&zBEhs`YjD{QT4YPy2QU{8Jj21uNOVx+h-o7T!CHrqdNTOMpw6$O@{
zaQ7E3OdIW+@7Lg-odA~f<BS(BgdSx8KRjatvUJ_YeNOW@gTX|w4x09SxOKV56j}Th
zF2@F#I@{YG)g7FzY>ICEs!uaB{)$j9;3R*Ns}N245fDXb-o;`t=^c6pBaaIr;$P1<
zJP7Kk4f?+D<lGzkYbH4)=R^+~W)YsrLhZkZtV5_`>8{KB+v$oEgDBDDr}p9g11-rk
zvvNM)ToCd)B+2nA$c0`ly`u(n290()T?wkrsvR9SwE!{g?bx7M2~@7hK%H>nXY?AH
zeI8?YcqHxfn~Y#g%8vpnjX8&EriTdeW=4<USSoE$p?l}XvPh&Bg8a!1w^I{fkul;R
zakRsyMJ9!C`3qa#4qNUHPxcOvD!TIX-%CuwMpI#{hA!QxJ$-8&A*b&#hy1<t{Q&L0
zL0myqe5G&!f5RB-=_h=s(Il<1ohXvmL+G5@ib@y%Y$>P{Hwea)y27aiCG-q)i4z-w
zExTCa>mgbm^$4zQuygsSqJ(kI!M1rX+&5KIQ<V6Fr83xwQa0a9lc}A47r@Y|<V-KG
zOz1+-N0R-l-wX%q1dZZ4P7kgJzE*-_>@@c?T*ukN<FP>o&3ZFeK-abplS%_Fq<>!)
z^TsJKCO?1<(O$uccRdaYYUuI72A|<|tAG{_3e%s&%W`n<kmrt*C?f~VdGJjq|B?+X
z4=7`fJs003J@x9dnMHy+<*2{v`)J|&kazwm9Sy3B<h(jJzclWl#nP)w263Gy$u5sn
z1&>ZbD|^a%uk_ca+`CJ2GTPUDC*^^GRQpwMvo9Av+?9OR=KPY8U7<@Y0&aD-^s|pB
zQA~?b?uup=!L$f!T#T^9rvy$t<c))<36BU>cl3sW_A8S4k$PbOVOI-yQP?^nSj52t
zVWieW%PEDs$1)wM=6>ey0N?5x8APOYL^iLeo{>1Ose={lA)R%h0;{3kAO5y$flmpg
zHUiLT38LqICh)UQHHnOCR79`2p3y9Ui*1BT5&W0x-SSz(p1pT^FXTd3e-95|Z5r+P
z!Z++Z@oDtA_(q-5q_!14MFoRuA-d#mS0+zl5)GW}NBKTX>}-~cm#We~i{PN_!E<iz
zssJ`tKU1$;OWv2CVefH%l|1Y0Z4k^)4!^&ka}664H+bvoa;E2tbhj3&p0)wneNsy~
z3t)h8gCxK7?%4oonfqh8FFU~Z%1^cG);6~8FKS8J_q43=`<ziV-7fL<IohpXy0eej
zR&Ko*xNU;=jHhlqSkT>Wte}63d9rMItV#ZI_=|nK?UT54U%1+VrfFUBoRm-x{a~cz
zMNp%R;0k9aXahz1M*FONeA4(3S}Kg#Y<Bunj&A1I-RYF)ju0Op3u&qCS{-7G6;n@h
z-BDh3PEdL->VD_hO?u01Z+E>S8MTM35**#ZJ;*FQ6O+DdeD*oe(|y8JY0TYI7r>gy
zc+dz~zuGq&$r|Fgf_Q4bT7`V`qP70<{kap_@ld8+yFE|ttEvb#0CgBZ&7+p^dfP6w
zcqogV`@ziAUyO~1hHk5HwY!Sh3`9~sg7A6>I<d)p!C1N=!PS0GMfOScEa8Hy$()F{
zt!lA=t3DL}tJL2t*l*RWE)3u);+dzE=)7pi!?LHC+8xWmz`ixaZ<puWSSdPj(^i1d
z=b6<|><@*#eX9xclkG+9<ex$~hJYc{*ri-JK3YIj`%<YpbVw#C;aeqnHZlM+Jd{_`
zc&q8`jqo>lK*D%UCEov33S77}6l!vOk$XIUgkGqRz|Q^5=8d{Y5v2=><fGbHu!rwm
z5Z^uS)g2yJ!1{V>3%EIxFgX#zg=VdSs>MDpdV&vcRG$3!>z1P?s?Z)u+p9T1oOxhn
zV$qT?-wE5T2A=I$tEw!)<tqbL9VM7cP4f!@q`w<VK0KFXvL~12PQR)aI)jMI6v6XX
zphW(WFx&H3LJp+FfP%;xa&HAv_?Nnn&hAQ2V_9U4<RFP^?D0{Fdvts8=W{m0h7RgB
z8XwXS0Q#6@;b!3`@pFh^Bb<5wzm<zS@y*}f0g<Vkt~}zqR3Wtdf(uUx)l}gf8k@hO
zU9uK;dE@-~WKxPqO%iZ@sbl^MUE%gKlR$XXybw43q^1e}5Q@;Fr6AlXVqBLx<1v$a
z(Jgf6?MmbDW%e@zK%WD)O+zOkYuH~Q(%=EXr>Zi{LBNk{!mC9+^ZJDFSzMF{+ud&;
zlZ1V2$H8d)Uqw*!#&IbnzfMM|_t6;G%_l=!Vn3CuyRSw}ZN5v-0q(p}wFTpfE~fCH
zv_dXi0mf3od$hlG!TI<4oYFjd#~8{h_0*mvFNJD`Vd6h8d4m9grR4zhQw?)7N08h!
z@zEL98o9$U@Dra5#Ny&xQ9hxuw9#vUgTbj;kCXV9bhq-8R9R8!Iw+KXeUD?_7!Jsa
zhcavGnA(|)+S^i1%g0sc88BOgrHByH0~E{#rSY#P>W>@|R%K%E6!FQ2TqTcfilQKe
zB6rlqPF8dOgbcvEb)A0cUU(x^%x_f>u*!{>n6zhdt9e%WAxT5H{IDqP{yw?5XOPM&
zEG0at?gi<Z#aTH3{g=cNablLJ-U!w7@NJo-UhPjL{ZmotrpwaDDQRWSRbZKxnI?bF
z<whPwH~a>MusPo)kr(7UoQ10Y?mSbx;<^6{`^4aokH@H5rqq+~`^MQ|rTLYf=*2HW
z%~d--xg|eZ<BgZ+wRj$c53~Mht_524u-N+|_|yHa9>TzURL~dKGB6OY12+p_o!G%f
zYiViK@gkS`r{8!?2l6psL@>&KzH~+(sa$ott>zs7B@c1B{t6tTI?c~46|OJ$^3P3w
z80h@LzY(whK-W`Xd)kT#*;E;m76UL~$ScxvXT|g`W{2VWiLYM7jQfkPk#+oIqn1|M
zn&qjq{}o@Jk^#P%VeKWuDpel!Hy#7Ym!vR%y?47kt#JIpvU1>KGtP>7Yd?R*u(uOD
zSuBz^HxiOxhW(`az}G^HAdq>qoG8c7Bc55(5EFl&H)WCf=A{Nt=?cy9>`ur08s)VS
z%hW|oB$^XBhnbnj&Nx3#^@hx-bro1xO}d}2ejiQqP(7fkDXFT%bu%n~FAZFfQT6wq
z9nIJ`yvqBAHlyrAVll*MC9X>U6&9Fm3=;;un4^z<rT_j_35Guvm4rw3Ts)&@iK*vu
zObl+l7Nia<VD1;m<T4D4aH$Ftsg&A`(DP8!i!SbDIt?DbBA3bkd2_WVO7CCF{G-qR
z(__Adv9P$BEi^A*aE<i#XP+9lz#3f8;7!$m9}t~JP5az3N-K!Z5P7^rc7L$-4C>Y5
zJJrhpaM$Q>I66iy{J;~suIhYCCnlu~C31!p?l}b8?^Rz3`=@A_&^Z$&&&|5aXFXQ5
zUySE`r&rn>ROT!lTqu8Sm7TCP{VZM}HY$3I1_17#q*H1mrTMC_qF8M?M<(zC_!)>k
zDz`ECamZFd?PC=QdvMO)bwF-<Qcp2c@HLP)vCA=OYk6p&hx>|z@FERAI9i9Z4N(iw
zT8UlNe>T{+$^PBe%ECn83M(*ZKJwokbnOTNKquJtgnu}Er$4fwwX;&whyo0AK^)%3
z<zZzc-n}iCMPF~=E0`WKl>L5Hv3Fb3F9qkV7EVs?WCDP8+x~R8S=>LFm5E}Dy1@~z
zB0o~N&lwVvBph8a(f%&!v0}<sYY(Zz!!K(^yzL*}(1yJqrEspk^~QW{l64-iSoj)p
z_H9_<<C(BeJlBpEGH%;Sn5VrEe`9JXeZQVJ@2doR*gTWSVH>N~ZX_V=x0=Okixu~@
zj8cf&B}$%@xpv;}P)?vvo<vfo(RXNhMWxylSnN|UY+8#~SjcqA@~C#+&G)rlKj{>8
zM1)FjqE{EUZVB0?78c(+EV&i@GziEdAQ3eGc}UmXW(PTk+;mjgRmer%I=-{OPqI2=
z^fQfal2qwSF=SW5P&z`r44bn;k8*INCXPQ2{l_;HCfw_vK>Vs=>ySTxbHT~XSAmc!
zB(!BdbJOTXk);kFdM<`#Y;m4rMc3p)hID-8?b}5|0~3!1I;~?AzAxl{IJ)@i>g>{Y
zBV7JY7@E|-H=IwVPS|Ifq?bpt0)U{P-TC(D(u>K0{)=JYkfr6dYdPw%p<{FKl`<~O
zcZ8P4!7TXActe8!*o5Ji;#wz&r3QTHh(rtdN~<fb)p{=cooY^QU26Z2vEt366?)?_
zHT`=-mD6(UHn)Md8#f;e_`~pubsi*PTkM?KmGPqJWN#xTW!W{3OPO5toidst5;Bo~
zc5`>b6_Pg^8Yd-WsyR|SCdMWPGSj?_?cuirO;;-_?VS$50jYkEkM!pz)l-U_g)W&k
zwf3$6v-(%R3NN2#oy5FksN<%G)lf$)=NJ>xi}%`DMZ$LfPY1s-L=i9DB9;VWdcHOs
z4(E-nE><;Efp%cihF5XEs;c!0R9B-d@(nsAy!p>mnr-fRzm?~BYWw@cXy$jiOLC|*
zjzy>Av|_VijJ4!B#^soJ+TV!KCsjW1*{uiQTj;>nZ%5mk71I4ANo29m9W~wuI8e=j
zgql_zY#b*ob54E2hz1pIr(zcbG#2?$T#l_jY*Ke9IYHKgO1#K-$ZubDT`^3)cNd_q
z=QAI`=K3;!wOx4abh034dFytk_VvT+p8A0>+{0k&c-QAZNqPEMBvP&lw$6N`8@BR@
zE2!W)@H)uVD#QNbecnX;%(62Ju$c1sTEYE(@q@KXLt`ETivlZ;2Ac?{P>J2RpS$+@
z0vM!wl=@LYbY&CN8NC$UUFw0JngDmq^PezHrMV&LW9n_T)Pf}E`6Cu)ogbfev%^?V
ztT}Nkg~7@+7{59)zXsmzzE`OIeAN`qr|hFW6-QnhS+9(*UJ0?Z%4b^Q_t6x82o=|-
zW@&Bo&9)lu%ICZK3A1*aa_)`hES4&D?JzD8#lF$NX0Qoa6U0SeF&w7C*@x)RBaD_m
zJjztfhI&P!iM*@HC#c860I(<~l6V{8xZxbKwG)(nmhL=&!KCcCF^Mm}<r{B@vt(jO
zEjLIj;P0l3(o#eDYBoyy#%kJUXPIubilaVFZ>}ypG612k0iTGljocak3$)qC^wT|R
z(<&ht(LA=>cd-8r<kr5UVcUGKMQ|OpIZEanMUMwPjibL?c++CiU;T`DNeHmxecN7G
zes40U^N^Kj>fgSloe)0}$`|D5$D~B1eH^r!8nMg!QMw**FLL$KDoBqgzd@|V4!Q9H
zSr9uU;X5{9wsxYJ3HUHPtDC8H=eK1-{>5=O6xb3P0y@?5kn;65N65(hApY?!?@o{=
zTa05}P%+{B$IzJmODi@H^0)l~z2Hxs>12aKC`Sx*7(Wc<9}K-_YN#zrx}!T2Y`ma=
z`C%q=;#_j0(@j$8EjnZSiqk?R&B<^#nkYmud-5^T9qk-!Mwi=<#Mp_<X|w!jk&W1k
zBkaR4BGqt@4m0w-1>cYeoQN2e_4w|eh(r*o;dk)%KgI5Vjl(fh<Tk$FZR`JGOKu8C
zR9zo@D9M`|J<rLUdq4GV)L%Kj5OCvjyr0W_w_lHm-F)~SoV2rtYzc@fz|3M@=Xv0!
zPskrB;rD5D>st6}n<G-_d-x8C|HSTzNI)v-5i7MQR%Pp9RYnQZmF9>gq?Y)3Di^Rw
zM3lhvzV*g>VlJ(Ymok4sdl%~VMAYO4Mk;51wW-~CLNMoUD#^l?isGHCS=)10oz!ew
zV1CxT#;H6fzzXy>aJ_;!3sjnGw=cw>ZT-2kZg;+eGNJkvec|U0$66rPBCoYuk9zss
zYF&3eg*FluT0?R#z&sPVom(FC<=sUYeCWLII*Kn5Src3{d)eGzeW%x;k`;(pTwC9`
z2a{L5g`oa;`B7R^e;<aJ4)T&In2vJIZ^>xa*q^Zh1u~mL%x<^J<k$xp@S@7*b3B!M
zWhKSI3k`Ka+h>58w50oK<v8`5q)iZNZr=vK`VEV69Z5!wLs=Fcz4}7yDrff;cLv-q
zHSAH>)*YWpTwyOqL{BtaU|xN+o*D1m+)vA`?PjU@T*(<!3L37dxCYJYtskWkcum8@
z0|=>1R17x|(M<!LWl*HOEbXU_WuW;To-|o|`sJIaYc`vF4a%dq{l@&2Z2Q(>2+z{$
z<hYx$0sh^AbJ+v&fI^ZJe@^xl=Ay%Q$B_W@EK*P7Z^Makp>mWv$#Y|ubON@Ou;2V!
zSYjSi@mHYf!`(7)vF^Im32haYwlE%UvW}||>$<OsW~k$<NSxcIDqLWL^`@UC_0jGy
z$(bD=nAAhJ#bq?~(>XyMt8014$UU$=rMdJl7_g`GkmA8^hq}f(hg)f*I<w@Vm9l}y
z#tMN<>k5TH#v~EOB)E3unp@r|$wPheqX*24^g^?uh2F7mdBkfjw~ZrQPE(woBrOwq
zCRZrC5BD^=!zz7q;h~d&WCk?7*=4=sBQR))c&jkv%n|T<*&w`DHFrqo;)6Z)OhvS8
z!$!+Jsbsnqu|j04U~*id6Za20%(43$;Y29;3>ouRQ&Tm!M~s^z6>rCfN>^p_-;(Qs
z)WZ<n&&NO8$w}z(OP*AV4WowSuD-8I;&|TmB`$5rPtd7rD}XGOe4xDDIFsM>>gL>Z
zS9fnieEm3VVJu1YdQB2wsxp04j!0Ts7RruJ^!4ISdPN59#_}kvMWwK?%MQ6o=1r`&
zhP9u*^ZJ_4Zj#c#!i=O4za4zDmroqVqtNjN<hHP|*c=)g2;VJfUvDbZ0-6K|ql(IR
zb5B>+@lLq-%rh^wg_XouSP5ZHcS~-k6Xqc0bf8-h-P|OBD7Vq-C=?$ju>TFijY56z
zNZtFi5zD&92=+N`vf6y-9W!Tz%&8|$M`VR~@CP3k(ECT{yBK7E=hq`Yg-QqC;5!If
zWn;BOv8yp}*i1BQ=T3?Et(UTcIX^x?hdPAXSBX@pYpD>@+mJ^|7Ekj8MV6eyUM26T
znZiF49II&;=lz}gB&)ox9VuwhXAPlpCxy1t{f;jfSxzVBPnxMJ0xi5v@|T-fhDG3=
z1f5_|`{AK4EM{8*q{><=_cZJkAm-Z?%*D@;ZS@=GyvR~9vec;MHjSa&v9gwmdgoX`
zBIW1NZ2TwA<(=aU9uCR(i`ZyD;3b<-pzH_J@FenmCvBP2;*K>ZD&YZn=TOC8d|W~=
zOLdSK=_mdaJK3|(J4MgmP{Q`nro8I1eyzwhJ;`7;KgHhI|IqGC)cE7HDkS{pHTO%`
z2Ct2Kv0OYl0W{UR`TAMMg%Y<gu9NLIGFq>=zQae2``dO)yqn1+7KUZ6Fh5QCi{DU+
zUA{}L3>AAE^UxuMfyA2vm2>Xk2%u}E>yxK{`6m6s7%lXu>c8vfys&XsaJ~zf3)drq
zuKSSK47*Pk=Rih20>a3$HQLY0U=?Fod#483bmv1<pM%?2kO5jAb(QHzSyFL=g$n9M
zZfRvN@>XE$N7Y<hsWJRctoYKU{@NJ4+Cy)X*;Xl+FdLT(7h|6IZC&d+B0ztVH-Yil
zL$6<9UVc+3Knth_9uLSl$XS{a<CDX4(sNMna!{Lm<sI!WP;8+#J<9YI9~WqxSgtyU
zaXPA3E@(+^B7avmIiY5tcCk^5r!+R%9ph<;b;dUV1yx&QOymyFWtp_RFPdvvlv#3Z
z&~AT)G@H4=FP;kqQKt@nEo)t|=a*R2e#-aja-A%IGhaj*(WX`B0FCtb`UG`X5OI>S
z=Qq5U0j+ox>Su39wAUyt(|`;38K`7F@-%r%SKe6)Z%!uq_F@G8DOc~gM;~(=4WHs^
zBmRwPtsg+*uaPfy?TQJ-q1;kar}SmYleVYj${%=M#@4rGILxPTZDY`g*N#6MUKDp>
z3N(AVt-CF|b6yqwC8Z6^#YiX);fe*Ug0)a;T^txWy!|twHyW)EYs5|sT)-wv!puzR
zgyx?<TxaSFqM}30P}<|*a$if@xsx89eWw8!j<M*FSJ(#k`M<6w^?CAE*vwD0g+xsv
zH%tZBk+v;;u9m~sS3qacnd_w<G0v_SFG1A_PNcx*GS(tr65b|s0@vhD9O)|GkUQAW
zR$(SeGGq9{`))10Iw0~~s-BhP=ATbt!ZbvJzxY~iL-)v6Z1kry$hx@RQ^_?L<7;yl
zm(=zSN^gKc)@ae-<fbkMEk<%&8q;X{7pvs(DBs^3<!@D#c*EbbtU9$TRGG-noA$cx
zn9kBNZ0VMpjNn4Dy_Z-v{QV<Bunw6oC8-Y=cc5dUM}7UE((g#IIO<2o!O93mitzgh
zW9fa^E5Ar|){j0i=<9F4b&Bgq6LH=J?IWHC++i&S9tXk<u)bDOP(h_f(;%iOm4EgZ
z*`9v{6ZZJ2xrm(V_v$ccS3DivB3(XMd*klKi5+dFrQggs6a4~bIHU4>5aiH;tWo(>
zQ8q0@(W~{t2R1Dg_+!>qYauN;^|;n3cO#IHJ1Hg9ual|DHG(^)=eOz$M3N?#2-ZKr
z=*}E|p(hg8oJKv~Fi=Oh7o+#Fg<2J7ld7GnU-3dFcjtwc2|=0<@byqcc$f&?jLs}4
z6^BGl(j-FZ$Qh?5QjfT-LTERI=~r^VF?7R+gsf(p4){2}nGY^o5t0YoWRZBjvvK=Y
zb!PT{Ns^2Y%Q$WfpQlS*h7>Dt{aRo^)Y-2$D5&-&-KO@OcD_HCBJGoAF^H-n;J08!
z?e_w%9tO&tJQ5dyvF%|rX<2cI>MR!vh3hxgH_X_7FVzplZyct8Xh<A(2|`$4ZM~Gn
zc2p}g7ayux7~cM&z%IIY(=r2ba%#RXGOwaRinqF~@4V7wqwj&Xi?0eKS_%*wKJDl)
z|1h^;O-GcO$)x>nimBVU^)YkvtaNP*@fgGL{cGE|b0<drd<0KLJ=71Mg|A9};B)c3
zu+)0-T^spq+^^sfI6Cv`9LLkm6~)k4Vy9Fqs=)Y%XPgd4-ukI%lJ&J<X{VgU4<QYv
z32Lwk@t?kB&K1rR_ORjQqPOV`zL!_PGEy9HcGs_{svinRFRu<=z4}~`ft`QI##)F=
z1)}$(+EUsgqyHM?VoCxj;$Y)vd(M{6)AlB<u=Q;rf4<DRUapVu(5nSGpSe-Y8(MWe
zf{e+cfR}t~L9V>0WnZCEZ>O)3c~#-1B+~b*DjK>RLW(J17vH<BmLE35lv~tJzKGU3
zrfSCxfl;dU=S>2T?#^Nocs9(B$TdR8Fb=j3e(ZafmVyZzL^uV(*mXqkjC&1!KC8!Q
z?zEg@!6dl=NH{r-JP<bwN4Nfxv5<Dw^*Do|-xfA$GOgPQr-V=HyiS)WZSl5Pi%%R?
zK^)rCHN9Fek}?=pTm}u=Yh)T|E6y+k7lyYINdMkzPaAdK;E3a~TxXQgdY$!HN2-g-
z{h3^+vcq(wv4C;{YR5~!*t+ZTw%A!p<b}HrwOVrqW@E&9=ZyH7S*J1P$S+A;!s5b)
z1mcpq_L#3NVo<39_jUCvFb|yQyl+);i~Dt|M*P_c%)l`Y)Lr}NL6vu9-xxY7ugp>{
z8R!fXFd-gMYU9b(s6Y^AKeF`GM3`5AkHPVj884lAS==^ToKqRiszewSCQLqP>1A>}
z^<(Ak``CBZFAspB3xr>oGl{PINE`g#55L9tJGvtV$Ay=0$h<t?kmbb;SNEcrx2B2z
z4Q>9zNC5w!9M=dftK*@#jIFQ8!uze$9_)8M%1(|yYAdBvrR>4I>573Z#o^$P8R2eT
zq`p7A0c|;W@AuZwXk=#`7KbL|vNWu9`#%;JO~U@Qn*U8)IO4=D?~VO5^o3-wEf(Hk
z5xz6V#Pr2TpdG1F&Z+Q%$JT5ec8tM&M4T+JJ{O!l9ydBSx~o5|d)_posOe*Bm)dV(
zUQ={U)Ae*6=8^x<;t+q9Pa_bydUspds7H@$wST*6XK-3*dci)usXMt98OmOn@b?@`
zZ)@o#JL5v+uFJ$uyk{#`R~7nSWXLVTblKMfDaN<9SyuMz#Kpv1)#f6t@2!(%F<c<u
z69EeMpHvJqVS=Zd-$g(0dz+@g%9;mTqWp$YCc_-JCNUYklp2Yf$<&aWCLhAJ{{zF0
zT2VZI7Q)Q_WA4TsXZUp(yC{x7C1c0tX={BZ+V*8{1e0?NA9hc}F(vmD^Ubfc$Zx7l
z7!CA5n~!w&ue`@7@zfbTAQ<%_Lq^i$rDUv5(KkQikM^I|Yo4P?%c`2b|5G=m$8cfJ
z3j03k_16|~_tCZyw5odOuz-bw6LBon0x?~g=RDn}z->|hSLAm|z(76+9^<wVI)XfX
zmLT;(3zZ0yK3GvW^S)5dx96yE3MmX3%9hLnY2)k!5^wDLH)vl(TB@!}#ec_09~*Jf
zP2G=%mw;q}p8-kJX-aWN=l)krBT<1vTz_UiUd47?21Gn!s<lr2X?V^-UG?b1hgc?N
z@plK*mWP+yjW2#c_;fN;k9yuEOkj%A*FZmNM~8c^RQBbg{r&PknY`m76Xl>PG>vcQ
za=`jAVEYegV@07tyxt~FC+^K2j7m?y>0pUj$c*FqS=jrH)&q*kGooHRyRbT>gM#=U
zOpPno_0RuzHgiT;j(;1BGrgXlldIJSEXtVrQ8B@KUfiQ*;OXXG;A4diuiuHUuC@B;
zzeDfbIrWBHta_}k#t7Z?;An5?Ad5`!@II&y#JPX5bze|hLdcb1gifq<t#XAqyiFBS
z>QUGqJ1~_M@oBiW6_6$&G^1jb?4zVXJe-dPX#Jlqy#}=|fqX4!tU6T@ir+a(QLCCw
zj9)dM+G+SAo=N^2qPfiFJSJ-&F4u>L|MB11`QgKPhR3`tryj=c^G9a$8%0Mjq+o<L
wyW3b@1Q`?rg$ZgL!^r4&fDYoV!9i)+o?9|TuZfHa7VUYWDXI*PPgX1b3s1EU2><{9

literal 0
HcmV?d00001

diff --git a/ROI_volmue_283_MNI.m b/ROI_volmue_283_MNI.m
index e2dfe62..b9bdfe9 100644
--- a/ROI_volmue_283_MNI.m
+++ b/ROI_volmue_283_MNI.m
@@ -1,16 +1,17 @@
 %function:统计一个目录下所有分割结果中283个脑区的ROI大小，放到一个txt文档中
 clc;clear all;close all;
-pathname_G1 = 'I:\NC30';
+pathname_G1 = 'G:\m2question';
 cd(pathname_G1);
 dirs_G1=dir([pathname_G1, '\*.*']);
 dircell_G1=struct2cell(dirs_G1)' ;
-filenames_G1=dircell_G1(3:length(dircell_G1),1);
+filenames_G1=dircell_G1(3:length(dircell_G1)-1,1);
 %file_path_G1=fullfile(pathname_G1,filenames_G1); %保存各网络参数的全路径
 dirs_G2=dir([pathname_G1, '\*.*']);
 dircell_G2=struct2cell(dirs_G2)' ;
-filenames_G2=dircell_G2(3:length(dircell_G2),1);
+filenames_G2=dircell_G2(3:length(dircell_G2)-1,1);
 savemat = zeros(length(filenames_G1),283);
  for i=1:length(filenames_G1)
+    cd(char(filenames_G1(i)));
     cd(char(filenames_G1(i)));
     cd output;
     clear dirs_G2;clear dircell_G2;clear filenames_G2;
@@ -21,7 +22,7 @@
     for j = 1:len
        
         filenamesIS = char(filenames_G2(j));
-        jamar = strfind(filenamesIS,'283Labels_M2_MNI.img');
+        jamar = strfind(filenamesIS,'283Labels_MNI.img');
         if length(jamar) ~= 0 
            
             [hdr,filetype,fileprefix,machine] = load_nii_hdr(filenamesIS);
diff --git a/level5_to_level3.m b/level5_to_level3_M1.m
similarity index 100%
rename from level5_to_level3.m
rename to level5_to_level3_M1.m
diff --git a/level5_to_level3_M2.m b/level5_to_level3_M2.m
new file mode 100644
index 0000000..b65046b
--- /dev/null
+++ b/level5_to_level3_M2.m
@@ -0,0 +1,105 @@
+%This programme trace back ROI information based on M2 segmention via a
+%lookuptable.this programme can also generate a new lookuptable.
+%Author: Challen zhangchanglehit@163.com
+%Date: Jun-14th-2017
+clc;
+clear all;
+close all;
+
+load('M2lookuptable.mat'); %para named 'a'
+
+b = a(:,9);
+hashset = unique(b);
+
+res = a;
+row_num = zeros(53,1);
+cnt = 1;
+for z =1:283
+%     flag = 0;
+    dim = size(hashset);
+    if(dim(2)~=0)
+        for j = 1:dim(1)
+            if isequal(res(z,9),hashset(j))
+%                 flag = 1;
+                row_num(cnt)=z;
+                hashset(j)=[];
+                cnt=cnt+1;
+                break;
+            end
+        end
+    end
+%     if(flag==0)
+%         res(i,:)=[];
+%     end
+end
+%after that we got important row num in array row_num
+
+cnt = 53;
+for z=283:-1:1
+    if(z~=row_num(cnt))
+        res(z,:)=[];
+    else
+        cnt=cnt-1;
+    end
+end
+%开始进行回溯
+path = 'G:\test';
+cd(path)
+dirs_G1=dir([path, '\*.*']);
+dircell_G1=struct2cell(dirs_G1)' ;
+filenames_G1=dircell_G1(3:length(dircell_G1),1);
+for m=1:length(filenames_G1)
+    floder_name = char(filenames_G1(m));
+    cd(floder_name);
+    disp(floder_name);
+    cd('output');
+    dirs_G2=dir;
+    dircell_G2=struct2cell(dirs_G2)' ;
+    filenames_G2=dircell_G2(3:length(dircell_G2),1);
+    len = length(filenames_G2);
+    for j = 1:len
+        
+        filenamesIS = char(filenames_G2(j));
+        jamar = strfind(filenamesIS,'283Labels_M2_MNI.img');
+        if length(jamar) ~= 0
+            
+            [hdr,filetype,fileprefix,machine] = load_nii_hdr(filenamesIS);
+            [img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine);
+            %开始回溯步骤
+            for indexer=1:53
+                for img_index = 1:283
+                    if isequal(res(indexer,9),a(img_index,9))
+                         if indexer~= 37
+                            index=(img==img_index); %index是图像中对应label number的体素序号
+                            for t=img_index
+                                index1=(img==t);
+                                index=index|index1;
+                            end
+                            img(index)=indexer;
+                         else
+                             index=(img==img_index); %index是图像中对应label number的体素序号
+                            for t=img_index
+                                index1=(img==t);
+                                index=index|index1;
+                            end
+                            img(index)=0;
+                         end
+                    end
+                end
+            end
+            %从1到137排列，因为总共有136个脑区
+            Lables=unique(img);
+            for z=1:53
+                temp=Lables(z); %label number
+                index=(img==temp); %index是图像中对应label number的体素序号
+                img(index)=z-1;
+            end
+            
+            CC_286_L=make_ana(img);     %savemat
+            savename = strcat(floder_name,'_M2_MNI_level3');
+            save_untouch_nii(CC_286_L,savename);
+            cd(path);    
+        end
+    end
+end
+            
diff --git a/level5_to_level4.m b/level5_to_level4.m
deleted file mode 100644
index 71d05c7..0000000
--- a/level5_to_level4.m
+++ /dev/null
@@ -1,1252 +0,0 @@
-% Author: Chenfei Ye upgraded by ZCL
-% 2015-11-21    (2016 Mar)
-%funtion: 将level5下的分割结果回溯到level4
-clc
-clear all;
-close all;  %initialization
-%cd 'F:\ADNI_AD\ADT1SEG\AD16\output'
-% name='Zhuxiaolan';
-[filename,pathname]=uigetfile({'*.*'},'choose any dcm file');
-cd(pathname);
-[hdr,filetype,fileprefix,machine] = load_nii_hdr(filename);
-[img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine);
-% img1=double(img);
-% img2=imrotate(img1,180);
-
-%AG_L:
-temp1=[31]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%AG_R:
-temp2=temp1+1; %label number
-index=(img==temp2(1)); %index是图像中对应label number的体素序号
-for t=temp2
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp2(1);
-
-%ALIC_L:
-temp1=[131]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%ALIC_R:
-temp2=temp1+1; %label number
-index=(img==temp2(1)); %index是图像中对应label number的体素序号
-for t=temp2
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp2(1);
-
-%Amyg_L:
-temp1=[73]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Amyg_R:
-temp2=temp1+1; %label number
-index=(img==temp2(1)); %index是图像中对应label number的体素序号
-for t=temp2
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp2(1);
-
-%ant_DPWM_L:
-temp1=[117,119,153]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%ant_DPWM_R:
-temp2=temp1+1; %label number
-index=(img==temp2(1)); %index是图像中对应label number的体素序号
-for t=temp2
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp2(1);
-
-%AnteriorLateralVentricle_L:
-temp1=[171,172,285,176]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%AnteriorLateralVentricle_R:
-temp1=[177,286]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%BasalForebrain_L:
-temp1=[85:2:89,157:2:165,168]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%BasalForebrain_R:
-temp1=[86:2:90,158:2:166,167]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%BCC_L:
-temp1=[125]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%BCC_R:
-temp1=[126]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Caud_L:
-temp1=[77]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Caud_R:
-temp1=[78]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Caudate_tail_L:
-temp1=[281]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Caudate_tail_R:
-temp1=[283]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CentralSul_L:
-temp1=[265]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CentralSul_R:
-temp1=[266]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Cerebellum_L:
-temp1=[96]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Cerebellum_R:
-temp1=[95]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CerebellumWM_L:
-temp1=[253,258]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CerebellumWM_R:
-temp1=[252,259]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CGC_L:
-temp1=[139]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CGC_R:
-temp1=[140]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CGH_L:
-temp1=[141]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CGH_R:
-temp1=[142]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Cingulate_L:
-temp1=[61:2:69]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Cingulate_R:
-temp1=[62:2:70]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CinguSul_L:
-temp1=[275]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%CinguSul_R:
-temp1=[276]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Cu_L:
-temp1=[57]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Cu_R:
-temp1=[58]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%FrontSul_L:
-temp1=[263]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%FrontSul_R:
-temp1=[264]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%FuG_L:
-temp1=[49]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%FuG_R:
-temp1=[50]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Fx/ST_L:
-temp1=[143]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Fx/ST_R:
-temp1=[144]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Fx_L:
-temp1=[145]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Fx_R:
-temp1=[146]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%GCC_L:
-temp1=[123]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%GCC_R:
-temp1=[124]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%GP_L:
-temp1=[81]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%GP_R:
-temp1=[82]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Hippo_L:
-temp1=[75,282]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Hippo_R:
-temp1=[76,284]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%IFG_L:
-temp1=[11:2:15]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%IFG_R:
-temp1=[12:2:16]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%IIIandIV_Ventricle:
-temp1=[181]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%inf_DPWM_L:
-temp1=[135:2:137,147:2:151]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%inf_DPWM_R:
-temp1=[136:2:138,148:2:152]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%InferiorLateralVentricle_L:
-temp1=[175]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%InferiorLateralVentricle_R:
-temp1=[180]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Insula_L:
-temp1=[71]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Insula_R:
-temp1=[72]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%IOG_L:
-temp1=[55]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%IOG_R:
-temp1=[56]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%ITG_L:
-temp1=[43]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%ITG_R:
-temp1=[44]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%LG_L:
-temp1=[59]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%LG_R:
-temp1=[60]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Limbic_L:
-temp1=[45:2:47]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Limbic_R:
-temp1=[46:2:48]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Medulla_L:
-temp1=[109,115]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Medulla_R:
-temp1=[110,116]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%MFG_L:
-temp1=[7,9]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%MFG_R:
-temp1=[8,10]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%midbrain_L:
-temp1=[91:2:93,97:2:99]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%midbrain_R:
-temp1=[92:2:94,98:2:100]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%MOG_L:
-temp1=[53]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%MOG_R:
-temp1=[54]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%MTG_L:
-temp1=[39:2;41]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%MTG_R:
-temp1=[40:2:42]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%OcciptSul_L:
-temp1=[277]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%OcciptSul_R:
-temp1=[278]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%OG_L:
-temp1=[17,19]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%OG_R:
-temp1=[18,20]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%ParietSul_L:
-temp1=[273]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%ParietSul_R:
-temp1=[274]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralCingulateWM_L:
-temp1=[242:2:250]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralCingulateWM_R:
-temp1=[243:2:251]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralFrontalWM_L:
-temp1=[186:2:206,210]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralFrontalWM_R:
-temp1=[187:2:207,211]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralOccipitalWM_L:
-temp1=[232:2:240]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralOccipitalWM_R:
-temp1=[233:2:241]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralParietalWM_L:
-temp1=[155,208,212:2:218]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralParietalWM_R:
-temp1=[156,209,213:2:219]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralTemporalWM_L:
-temp1=[220:2:230]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PeripheralTemporalWM_R:
-temp1=[221:2:231]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PLIC_L:
-temp1=[133]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PLIC_R:
-temp1=[134]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PoCG_L:
-temp1=[23]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PoCG_R:
-temp1=[24]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Pons_L:
-temp1=[101:2:107,111,113,261]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Pons_R:
-temp2=temp1+1; %label number
-index=(img==temp2(1)); %index是图像中对应label number的体素序号
-for t=temp2
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp2(1);
-
-%post_DPWM_L:
-temp1=[121]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%post_DPWM_R:
-temp1=[122]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PosteriorLateralVentricle_L:
-temp1=[173,174]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PosteriorLateralVentricle_R:
-temp1=[178,179]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PrCG_L:
-temp1=[25]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PrCG_R
-temp1=[26]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PrCu_L:
-temp1=[33]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PrCu_R:
-temp1=[34]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Put_L:
-temp1=[79]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Put_R:
-temp1=[80]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PVA_anterior_L:
-temp1=[182]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PVA_anterior_R:
-temp1=[183]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PVA_posterior_L:
-temp1=[129,184]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%PVA_posterior_R:
-temp1=[130,185]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%RG_L:
-temp1=[21]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%RG_R:
-temp1=[22]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SCC_L:
-temp1=[127]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SCC_R:
-temp1=[128]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SFG_L:
-temp1=[1:2:5]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SFG_R:
-temp1=[2:2:6]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SMG_L:
-temp1=[29]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SMG_R:
-temp1=[30]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SOG_L:
-temp1=[51]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SOG_R:
-temp1=[52]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SPG_L:
-temp1=[27]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SPG_R:
-temp1=[28]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%STG_L:
-temp1=[35,37]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%STG_R:
-temp1=[36,38]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SylFrontSul_L:
-temp1=[267]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SylFrontSul_R:
-temp1=[268]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SylParieSul_L:
-temp1=[271]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SylParieSul_L:
-temp1=[272]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SylTempSul_L:
-temp1=[269]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%SylTempSul_R:
-temp1=[270]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%TempSul_L:
-temp1=[279]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%TempSul_R:
-temp1=[280]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);  
-
-%Thalamus_L:
-temp1=[83]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-%Thalamus_R:
-temp1=[84]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=temp1(1);
-
-
-%xxxxx:
-temp1=[169,170,254,255,256,257,260]; %label number
-index=(img==temp1(1)); %index是图像中对应label number的体素序号
-for t=temp1
-    index1=(img==t);
-    index=index|index1;
-end
-img(index)=0;
-
-%从1到137排列
-Lables=unique(img);
-for i=1:137
-    temp=Lables(i); %label number
-    index=(img==temp); %index是图像中对应label number的体素序号
-    img(index)=i-1;
-end
-
-CC_286_L=make_ana(img);     %savemat
-
-%cd(['H:\DTI_file\',name,'\fiber']) %自行更改路径
-save_untouch_nii(CC_286_L,'AD01_M1_MNI_level4');%save_untouch_nii(CC_286_R,'CC_286_R');
\ No newline at end of file
diff --git a/level5_to_level4_M1.m b/level5_to_level4_M1.m
new file mode 100644
index 0000000..ae71d5f
--- /dev/null
+++ b/level5_to_level4_M1.m
@@ -0,0 +1,1287 @@
+% Author: Chenfei Ye upgraded by ZCL, corresponding author GuoXin
+% 2015-11-21    (2016 Mar, May 2017)
+%funtion: 将level5下的分割结果回溯到level4,并实现了批处理功能。
+clc
+clear all;
+close all;  %initialization
+%cd 'F:\ADNI_AD\ADT1SEG\AD16\output'
+% name='Zhuxiaolan';
+
+% 手动选择方法
+% [filename,pathname]=uigetfile({'*.*'},'choose any dcm file');
+% cd(pathname);
+% [hdr,filetype,fileprefix,machine] = load_nii_hdr(filename);
+% [img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine);
+%end of 手动选择方法
+
+% img1=double(img);
+% img2=imrotate(img1,180);
+
+%批处理使用
+
+path = 'J:\宣武医院硬盘备份\PD\seg';
+cd(path)
+dirs_G1=dir([path, '\*.*']);
+dircell_G1=struct2cell(dirs_G1)' ;
+filenames_G1=dircell_G1(3:length(dircell_G1),1);
+for i=1:length(filenames_G1)
+    floder_name = char(filenames_G1(i));
+    cd(floder_name);
+    disp(floder_name);
+    cd('output');
+    dirs_G2=dir;
+    dircell_G2=struct2cell(dirs_G2)' ;
+    filenames_G2=dircell_G2(3:length(dircell_G2),1);
+    len = length(filenames_G2);
+    for j = 1:len
+        
+        filenamesIS = char(filenames_G2(j));
+        jamar = strfind(filenamesIS,'286Labels_MNI.img');
+        if length(jamar) ~= 0
+            
+            [hdr,filetype,fileprefix,machine] = load_nii_hdr(filenamesIS);
+            [img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine);
+            %开始回溯步骤
+            %AG_L:
+            temp1=[31]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %AG_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %ALIC_L:
+            temp1=[131]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ALIC_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %Amyg_L:
+            temp1=[73]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Amyg_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %ant_DPWM_L:
+            temp1=[117,119,153]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ant_DPWM_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %AnteriorLateralVentricle_L:
+            temp1=[171,172,285,176]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %AnteriorLateralVentricle_R:
+            temp1=[177,286]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BasalForebrain_L:
+            temp1=[85:2:89,157:2:165,168]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BasalForebrain_R:
+            temp1=[86:2:90,158:2:166,167]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BCC_L:
+            temp1=[125]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BCC_R:
+            temp1=[126]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caud_L:
+            temp1=[77]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caud_R:
+            temp1=[78]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caudate_tail_L:
+            temp1=[281]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caudate_tail_R:
+            temp1=[283]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CentralSul_L:
+            temp1=[265]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CentralSul_R:
+            temp1=[266]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cerebellum_L:
+            temp1=[96]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cerebellum_R:
+            temp1=[95]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CerebellumWM_L:
+            temp1=[253,258]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CerebellumWM_R:
+            temp1=[252,259]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGC_L:
+            temp1=[139]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGC_R:
+            temp1=[140]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGH_L:
+            temp1=[141]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGH_R:
+            temp1=[142]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cingulate_L:
+            temp1=[61:2:69]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cingulate_R:
+            temp1=[62:2:70]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CinguSul_L:
+            temp1=[275]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CinguSul_R:
+            temp1=[276]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cu_L:
+            temp1=[57]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cu_R:
+            temp1=[58]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FrontSul_L:
+            temp1=[263]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FrontSul_R:
+            temp1=[264]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FuG_L:
+            temp1=[49]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FuG_R:
+            temp1=[50]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx/ST_L:
+            temp1=[143]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx/ST_R:
+            temp1=[144]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx_L:
+            temp1=[145]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx_R:
+            temp1=[146]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GCC_L:
+            temp1=[123]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GCC_R:
+            temp1=[124]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GP_L:
+            temp1=[81]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GP_R:
+            temp1=[82]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Hippo_L:
+            temp1=[75,282]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Hippo_R:
+            temp1=[76,284]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IFG_L:
+            temp1=[11:2:15]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IFG_R:
+            temp1=[12:2:16]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IIIandIV_Ventricle:
+            temp1=[181]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %inf_DPWM_L:
+            temp1=[135:2:137,147:2:151]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %inf_DPWM_R:
+            temp1=[136:2:138,148:2:152]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %InferiorLateralVentricle_L:
+            temp1=[175]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %InferiorLateralVentricle_R:
+            temp1=[180]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Insula_L:
+            temp1=[71]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Insula_R:
+            temp1=[72]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IOG_L:
+            temp1=[55]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IOG_R:
+            temp1=[56]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ITG_L:
+            temp1=[43]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ITG_R:
+            temp1=[44]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %LG_L:
+            temp1=[59]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %LG_R:
+            temp1=[60]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Limbic_L:
+            temp1=[45:2:47]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Limbic_R:
+            temp1=[46:2:48]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Medulla_L:
+            temp1=[109,115]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Medulla_R:
+            temp1=[110,116]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MFG_L:
+            temp1=[7,9]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MFG_R:
+            temp1=[8,10]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %midbrain_L:
+            temp1=[91:2:93,97:2:99]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %midbrain_R:
+            temp1=[92:2:94,98:2:100]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MOG_L:
+            temp1=[53]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MOG_R:
+            temp1=[54]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MTG_L:
+            temp1=[39:2:41]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MTG_R:
+            temp1=[40:2:42]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OcciptSul_L:
+            temp1=[277]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OcciptSul_R:
+            temp1=[278]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OG_L:
+            temp1=[17,19]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OG_R:
+            temp1=[18,20]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ParietSul_L:
+            temp1=[273]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ParietSul_R:
+            temp1=[274]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralCingulateWM_L:
+            temp1=[242:2:250]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralCingulateWM_R:
+            temp1=[243:2:251]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralFrontalWM_L:
+            temp1=[186:2:206,210]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralFrontalWM_R:
+            temp1=[187:2:207,211]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralOccipitalWM_L:
+            temp1=[232:2:240]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralOccipitalWM_R:
+            temp1=[233:2:241]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralParietalWM_L:
+            temp1=[155,208,212:2:218]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralParietalWM_R:
+            temp1=[156,209,213:2:219]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralTemporalWM_L:
+            temp1=[220:2:230]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralTemporalWM_R:
+            temp1=[221:2:231]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PLIC_L:
+            temp1=[133]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PLIC_R:
+            temp1=[134]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PoCG_L:
+            temp1=[23]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PoCG_R:
+            temp1=[24]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Pons_L:
+            temp1=[101:2:107,111,113,261]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Pons_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %post_DPWM_L:
+            temp1=[121]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %post_DPWM_R:
+            temp1=[122]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PosteriorLateralVentricle_L:
+            temp1=[173,174]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PosteriorLateralVentricle_R:
+            temp1=[178,179]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCG_L:
+            temp1=[25]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCG_R
+            temp1=[26]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCu_L:
+            temp1=[33]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCu_R:
+            temp1=[34]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Put_L:
+            temp1=[79]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Put_R:
+            temp1=[80]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_anterior_L:
+            temp1=[182]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_anterior_R:
+            temp1=[183]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_posterior_L:
+            temp1=[129,184]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_posterior_R:
+            temp1=[130,185]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %RG_L:
+            temp1=[21]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %RG_R:
+            temp1=[22]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SCC_L:
+            temp1=[127]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SCC_R:
+            temp1=[128]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SFG_L:
+            temp1=[1:2:5]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SFG_R:
+            temp1=[2:2:6]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SMG_L:
+            temp1=[29]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SMG_R:
+            temp1=[30]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SOG_L:
+            temp1=[51]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SOG_R:
+            temp1=[52]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SPG_L:
+            temp1=[27]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SPG_R:
+            temp1=[28]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %STG_L:
+            temp1=[35,37]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %STG_R:
+            temp1=[36,38]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylFrontSul_L:
+            temp1=[267]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylFrontSul_R:
+            temp1=[268]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylParieSul_L:
+            temp1=[271]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylParieSul_L:
+            temp1=[272]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylTempSul_L:
+            temp1=[269]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylTempSul_R:
+            temp1=[270]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %TempSul_L:
+            temp1=[279]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %TempSul_R:
+            temp1=[280]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Thalamus_L:
+            temp1=[83]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Thalamus_R:
+            temp1=[84]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            
+            %xxxxx:
+            temp1=[169,170,254,255,256,257,260]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=0;
+            
+            %从1到137排列，因为总共有136个脑区
+            Lables=unique(img);
+            for i=1:136
+                temp=Lables(i); %label number
+                index=(img==temp); %index是图像中对应label number的体素序号
+                img(index)=i-1;
+            end
+            
+            CC_286_L=make_ana(img);     %savemat
+            savename = strcat(floder_name,'_M1_MNI_level4');
+            %cd(['H:\DTI_file\',name,'\fiber']) %自行更改路径
+            save_untouch_nii(CC_286_L,savename);
+            cd(path);     
+        end
+        
+    end
+    
+end
diff --git a/level5tolevel4_manual_select.m b/level5tolevel4_manual_select.m
new file mode 100644
index 0000000..d3c6ce6
--- /dev/null
+++ b/level5tolevel4_manual_select.m
@@ -0,0 +1,1288 @@
+% Author: Chenfei Ye upgraded by ZCL, corresponding author GuoXin
+% 2015-11-21    (2016 Mar, May 2017)
+%funtion: 将level5下的分割结果回溯到level4,仅实现手动点击功能。
+clc
+clear all;
+close all;  %initialization
+%cd 'F:\ADNI_AD\ADT1SEG\AD16\output'
+% name='Zhuxiaolan';
+
+% 手动选择方法
+ [filename,pathname]=uigetfile({'*.*'},'choose any dcm file');
+ cd(pathname);
+ [hdr,filetype,fileprefix,machine] = load_nii_hdr(filename);
+ [img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine);
+%end of 手动选择方法
+
+% img1=double(img);
+% img2=imrotate(img1,180);
+
+%批处理使用
+
+% path = 'G:\玄武医院数据\韩璎课题组数据\NC\NC_分割结果';
+% cd(path)
+% dirs_G1=dir([path, '\*.*']);
+% dircell_G1=struct2cell(dirs_G1)' ;
+% filenames_G1=dircell_G1(3:length(dircell_G1),1);
+% for i=1:length(filenames_G1)
+%     floder_name = char(filenames_G1(i));
+%     cd(floder_name);
+%     disp(floder_name);
+%     cd('output');
+%     dirs_G2=dir;
+%     dircell_G2=struct2cell(dirs_G2)' ;
+%     filenames_G2=dircell_G2(3:length(dircell_G2),1);
+%     len = length(filenames_G2);
+%     for j = 1:len
+%         
+%         filenamesIS = char(filenames_G2(j));
+%         jamar = strfind(filenamesIS,'286Labels_MNI.img');
+%         if length(jamar) ~= 0
+%             
+%             [hdr,filetype,fileprefix,machine] = load_nii_hdr(filenamesIS);
+%             [img,hdr] = load_nii_img(hdr,filetype,fileprefix,machine);
+            %开始回溯步骤
+            %AG_L:
+            temp1=[31]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %AG_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %ALIC_L:
+            temp1=[131]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ALIC_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %Amyg_L:
+            temp1=[73]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Amyg_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %ant_DPWM_L:
+            temp1=[117,119,153]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ant_DPWM_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %AnteriorLateralVentricle_L:
+            temp1=[171,172,285,176]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %AnteriorLateralVentricle_R:
+            temp1=[177,286]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BasalForebrain_L:
+            temp1=[85:2:89,157:2:165,168]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BasalForebrain_R:
+            temp1=[86:2:90,158:2:166,167]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BCC_L:
+            temp1=[125]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %BCC_R:
+            temp1=[126]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caud_L:
+            temp1=[77]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caud_R:
+            temp1=[78]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caudate_tail_L:
+            temp1=[281]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Caudate_tail_R:
+            temp1=[283]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CentralSul_L:
+            temp1=[265]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CentralSul_R:
+            temp1=[266]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cerebellum_L:
+            temp1=[96]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cerebellum_R:
+            temp1=[95]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CerebellumWM_L:
+            temp1=[253,258]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CerebellumWM_R:
+            temp1=[252,259]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGC_L:
+            temp1=[139]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGC_R:
+            temp1=[140]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGH_L:
+            temp1=[141]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CGH_R:
+            temp1=[142]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cingulate_L:
+            temp1=[61:2:69]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cingulate_R:
+            temp1=[62:2:70]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CinguSul_L:
+            temp1=[275]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %CinguSul_R:
+            temp1=[276]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cu_L:
+            temp1=[57]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Cu_R:
+            temp1=[58]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FrontSul_L:
+            temp1=[263]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FrontSul_R:
+            temp1=[264]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FuG_L:
+            temp1=[49]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %FuG_R:
+            temp1=[50]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx/ST_L:
+            temp1=[143]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx/ST_R:
+            temp1=[144]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx_L:
+            temp1=[145]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Fx_R:
+            temp1=[146]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GCC_L:
+            temp1=[123]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GCC_R:
+            temp1=[124]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GP_L:
+            temp1=[81]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %GP_R:
+            temp1=[82]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Hippo_L:
+            temp1=[75,282]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Hippo_R:
+            temp1=[76,284]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IFG_L:
+            temp1=[11:2:15]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IFG_R:
+            temp1=[12:2:16]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IIIandIV_Ventricle:
+            temp1=[181]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %inf_DPWM_L:
+            temp1=[135:2:137,147:2:151]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %inf_DPWM_R:
+            temp1=[136:2:138,148:2:152]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %InferiorLateralVentricle_L:
+            temp1=[175]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %InferiorLateralVentricle_R:
+            temp1=[180]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Insula_L:
+            temp1=[71]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Insula_R:
+            temp1=[72]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IOG_L:
+            temp1=[55]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %IOG_R:
+            temp1=[56]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ITG_L:
+            temp1=[43]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ITG_R:
+            temp1=[44]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %LG_L:
+            temp1=[59]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %LG_R:
+            temp1=[60]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Limbic_L:
+            temp1=[45:2:47]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Limbic_R:
+            temp1=[46:2:48]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Medulla_L:
+            temp1=[109,115]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Medulla_R:
+            temp1=[110,116]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MFG_L:
+            temp1=[7,9]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MFG_R:
+            temp1=[8,10]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %midbrain_L:
+            temp1=[91:2:93,97:2:99]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %midbrain_R:
+            temp1=[92:2:94,98:2:100]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MOG_L:
+            temp1=[53]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MOG_R:
+            temp1=[54]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MTG_L:
+            temp1=[39:2:41]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %MTG_R:
+            temp1=[40:2:42]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OcciptSul_L:
+            temp1=[277]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OcciptSul_R:
+            temp1=[278]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OG_L:
+            temp1=[17,19]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %OG_R:
+            temp1=[18,20]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ParietSul_L:
+            temp1=[273]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %ParietSul_R:
+            temp1=[274]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralCingulateWM_L:
+            temp1=[242:2:250]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralCingulateWM_R:
+            temp1=[243:2:251]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralFrontalWM_L:
+            temp1=[186:2:206,210]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralFrontalWM_R:
+            temp1=[187:2:207,211]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralOccipitalWM_L:
+            temp1=[232:2:240]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralOccipitalWM_R:
+            temp1=[233:2:241]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralParietalWM_L:
+            temp1=[155,208,212:2:218]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralParietalWM_R:
+            temp1=[156,209,213:2:219]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralTemporalWM_L:
+            temp1=[220:2:230]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PeripheralTemporalWM_R:
+            temp1=[221:2:231]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PLIC_L:
+            temp1=[133]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PLIC_R:
+            temp1=[134]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PoCG_L:
+            temp1=[23]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PoCG_R:
+            temp1=[24]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Pons_L:
+            temp1=[101:2:107,111,113,261]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Pons_R:
+            temp2=temp1+1; %label number
+            index=(img==temp2(1)); %index是图像中对应label number的体素序号
+            for t=temp2
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp2(1);
+            
+            %post_DPWM_L:
+            temp1=[121]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %post_DPWM_R:
+            temp1=[122]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PosteriorLateralVentricle_L:
+            temp1=[173,174]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PosteriorLateralVentricle_R:
+            temp1=[178,179]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCG_L:
+            temp1=[25]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCG_R
+            temp1=[26]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCu_L:
+            temp1=[33]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PrCu_R:
+            temp1=[34]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Put_L:
+            temp1=[79]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Put_R:
+            temp1=[80]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_anterior_L:
+            temp1=[182]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_anterior_R:
+            temp1=[183]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_posterior_L:
+            temp1=[129,184]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %PVA_posterior_R:
+            temp1=[130,185]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %RG_L:
+            temp1=[21]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %RG_R:
+            temp1=[22]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SCC_L:
+            temp1=[127]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SCC_R:
+            temp1=[128]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SFG_L:
+            temp1=[1:2:5]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SFG_R:
+            temp1=[2:2:6]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SMG_L:
+            temp1=[29]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SMG_R:
+            temp1=[30]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SOG_L:
+            temp1=[51]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SOG_R:
+            temp1=[52]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SPG_L:
+            temp1=[27]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SPG_R:
+            temp1=[28]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %STG_L:
+            temp1=[35,37]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %STG_R:
+            temp1=[36,38]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylFrontSul_L:
+            temp1=[267]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylFrontSul_R:
+            temp1=[268]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylParieSul_L:
+            temp1=[271]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylParieSul_L:
+            temp1=[272]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylTempSul_L:
+            temp1=[269]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %SylTempSul_R:
+            temp1=[270]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %TempSul_L:
+            temp1=[279]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %TempSul_R:
+            temp1=[280]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Thalamus_L:
+            temp1=[83]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            %Thalamus_R:
+            temp1=[84]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=temp1(1);
+            
+            
+            %xxxxx:
+            temp1=[169,170,254,255,256,257,260]; %label number
+            index=(img==temp1(1)); %index是图像中对应label number的体素序号
+            for t=temp1
+                index1=(img==t);
+                index=index|index1;
+            end
+            img(index)=0;
+            
+            %从1到137排列，因为总共有136个脑区
+            Lables=unique(img);
+            for i=1:136
+                temp=Lables(i); %label number
+                index=(img==temp); %index是图像中对应label number的体素序号
+                img(index)=i-1;
+            end
+            
+            CC_286_L=make_ana(img);     %savemat
+           % savename = strcat(floder_name,'_M1_MNI_level4');
+           savename = 'AD28_M1_MNI_level4';
+            %cd(['H:\DTI_file\',name,'\fiber']) %自行更改路径
+            save_untouch_nii(CC_286_L,savename);
+            cd(path);     
+        %end
+        
+    %end
+    
+%end
diff --git a/look_up_table_M2_level3.mat b/look_up_table_M2_level3.mat
new file mode 100644
index 0000000000000000000000000000000000000000..4b8ac6f63869655c8250bb584673e48193c0630c
GIT binary patch
literal 3019
zcma)+dpr{g8^>MB<&CivDPxgz5jt{g7(&M_NeW@)xZh_(=CXzEWs^e3W^|Frr9073
zVW>0Yy4Z$c)37wQwwb-o`+na4-{-IAzu)uyJm1gH>-gD|$DP1dW)@(t<7baV#U@0W
zf=?%erA8$q6K%k!o&DWUnSz~>kzuKk5#Zz`8}QkfG_Xfn68NYU_^74LF{q6d1pJ%D
zQETx3J#zmWd3P_3&E0vMoLmNRaV#kLo5e1r3qrZ7s>sD``!qJHY3<Kay6JJF!90bB
zr3bFW6htU^l(|byB=0^e=SKSNn-tzk-+E@ac?O+GNxyyA>G~_C(+*a9)vrL{koVq$
zJc>X1r_2i2{NQP)Fl&Y<y!B1!e?}+Ak7e89I(djQ>evKmQQ+LHoREX7rV=VA6yUn4
z6O|M4@U5x-j{aNe=LX=4^t1uG9NlYxu0T&5pv%+Y19U}t&Hzk~?mPffpob5@<mtWx
zOL@z>98&)BR!&IX@=i`j{<1oUk+%%uF!GlbIre$WS{#9+P&YNyQMfZzazj#r<mxsQ
zB1axiHjHGcKWS3W@!_9MDwusNZ<=iRjTt6r4`;|snKrW{-XvdaYwsd<=>0~~d`*HV
z>@hJUe<3k-!_8zjVxHM)!slJcWq!3QwvEOiulhiAx7=~YC6asxkx#QTj_OxU-1Edu
zizf$=PZL`6S(r&pQqwfX%^KkiuX-3J?MWr-%<rN$)vfKH0F7oE_oG7tUUt=1*i~=v
zzD8|^ueC*uGXES6EiF3dfF8xrVLo-zGQ05hlC%(6Z@(D(ydycD)K`$ysC|1K<?6*8
zht7^af!YQFijRE$E=5ZD9o)lX7AEns`^TPu4ldgt>wnj2e{U^&-&(H^Uwf`Fe4Y3T
z`Aw&m1)wYn9wP%-`s~-l4pg*ypAjV$9#PP`x?bS};75~CH$=Xp$fuKwLUKx+6Cahp
zBG#oewa83|=FsYiGOS%=WV!Wj#*bkxi(TBgr`*4U{bmk*q}+D|$q`Rq#?=<gd^$0e
zv5^{#?0GSWQ70e&x}CxQi~Sty`h^OAFLfxkWXfI$D-3=zDljoNSZgcqIiri32Ln4k
zuD^?#TbGxUlgSqEGm{AxniT8U`j%z?*W@+x1(r=qgj<+%Qn9^ZX-mCvcwug>lozd4
z+5Ph7Cuq*N<I9J)X7O$8)=w8*W}if<nQ;CVh)l;Wz5~qJ6mj7a+41=HVs5kaT>ry2
zM;nXyzYckfk{~_h{0~zHT6;d{f+5Ol#33Lx_k)k%kjif3LTS$+-n<(VNg=#)5jk!7
z`^wpwYx(;nT6GuL4+Iz4K>VV^QNdx!?PPdBxabz`>JN{Bdq5bLjJ(*x_Ez$zh3(p}
z0s3+${8?jF4?4}qGNxQW1g(o>5cfa_DlE7=(Nka*_RI}<QJoxj5JJsNNGT=m+U8+%
z+#b@v0|reh_VMC0L_?@cJds1J$F?I6ruAmOyy${Ws`U-u;4?O_t}&l9u)MEGW)fxx
zx8y=?JL%m%UB%m)Nu);$x~+Je8t&!TvDPQCFO`;EE7z;~<?~RJ5)i>*uX#Mp@02}0
zT{g@AmFxM%?><w$pY<1OBjou!4HgA?=80%aCxX65hat4fLq7)>#^)YncK0e$4u<^R
z{#)kIq1<xsetDtsSL2TX<sZbV1w%ru@HNCgWbF67paIM=ss4>`P++FP39WgriFfHa
z>DE)|sY~M5&4e3)|M2CM)tlu9??}_V(Eanw9LWUdn-e54#Bhy6c2l4sZ-t$Lwkz01
z4)mL3eoTe^Z1W|asHZv%J8YYH9je`wtmH~Vp317cZytMoG;ck-?)65Oh6Z|m$(tFV
zH!oUEoA>3C9m`7ay^;rp_)C(=o|Kl*DA+3!)*ZiOn?(`9t8Xi`jcGt*NTH$gCvpBR
z9v^+rxCxPf>tls1ZFAe|Z;Hw@WA>0F3~fdq27I#ZXJ<;EK@c`4&HM27X^-@!N|6t{
znTmf!lKJ<CZdU6*#`<j5swTy^IGx^-UMx5v9;ck5U2=<X0X)R}XxCX45h|xnvL^b+
zwz!#I+*NS`ZShqV-YH|-vR;+j^2_?mpYI9t?uS-~^+D^m8?Fb1;Xj&PfvfI$-QFG}
zwG1<u;HfLGLyu$>N8f(i?|p~={2IzQr2`*dGV&}UZ&=H^ZI6sh!A?GEI)@M42l!>#
z7jtKD4%OFZ_*cM;rVjoVrHfBpne@>6(at6#77^EtsdrfK?m(ulrM^)_S9eU~SY&DY
zHEvRvUMVVWZFweTW!;=4ff4!o>`6ZB!4R~32;^?fSx^eGOJJC*rv~!g=)mAC*D8%M
zi~6=NL1M4=?^DguS4C7!5ze(a`w9cXtO;xL<+q(UT^CnUCOi^I9{p8Xq@V3juS3Xc
zdlACN$lSl?J7_8YDOM;a%tOnZSE038><qGL>(TSmnq-j_oYG;Ex?&m61I=fWCH>Lv
z8g1T(#kYR*Yd>ZgpKD1&ocum<%}ce0vdVLAd8)f3L8AA0F9_#u&&<v%(-_dYZ103K
z*h?`GHS>Y1Ev#pP_Lz-M3kmyK2nu_;fekdn4sb{FYa2}X5dp<py=e^%cl4HW{urUC
z<bg#)3ta?~N#lo3t4}1)AgVmOC{wlrT9vTp&Q@wViIN2gH3s}zDLXND(|QaaXYC`%
zaLH<GFlMBElB|?crk`r*;2s&^<3PXAE3I#zlim^)NADhLj4he-WXPBZJ*mj2{__RK
z;UE_3O_T0ki0r=_l_thdY-(gX?jdxeF_!KwMm8lRsRJyjF`-IU=Yn!g_t5kwnkvlR
zM!H@P6x=AvRWl5BbU1afdPr(?{(9Mu#PqKigW;Cb+_UwTZ$i>#?cW!5;>Zh_X<Odn
zJ-0vP?`6zHK&a+rrB^y7kDx}_wvz~|8em-Yb22T702W4@u`}Wql_l&XipTH!RQh{Q
zuUtPdf_Q1CPsdbejQ?z@F*ZuO8_sn`Jzo>ML3#=Qf;YmVuR5@mlt-=40QGkkF+4uH
z5KJ_PU7|i3oULkRUuJ1!&)bYRT#Ivo@+^6|cUOru5=0!L<J7hnjiBL;FH3rj1=0(R
z3-C|VP8H3RwSxcw%<v%t<DEP@o;gU=UHg+ui1YzLF0&)3dK2W}oi7t=cTT^5S}Cwg
z(Pk#^Unb9bq6($WIms&<G1Q)l)I2_I>XP)_@Y$D6`(32lF<w_PGt!ATDZj5A*e+Dk
zS%l!4+Ou{UVt`+&2}&?ZETDz=%o5lD``Izbql|1}-~Aw;#V@U{HGa)1)Ki$lvl-C|
zQjZnPdW|*LGlIDXvb-_={X$o#GwT8An6k)hk-Z?cNQU6DvJ=vDLk)r_x_BAG3m%Jh
znp<;ruPzu;WskfQ+2h<No`(8`E4v`gZMv(m1Jy7S&JbEWtUZQVUfy|C8HXYz@@>>|
z!mni`esiV1Sn%Zr@7sin?@XbNp}0c_=|2WHx$c>*ecYD0N?YUXpo|)X>4UB)|NYlc
zbRyr|wvn|q(xEEmkFQo2piZFtit%xoALzG$%2uvhv3U*@PQSx{G+Qe5iT?!Fq=1(G
z7)ydXco+x+0uSW$2rf_r0SSj(Ea+}ugcZ?|;oGBP0KK17T_Mvlexv&;(X!j`_vM6>
z3Abrgbn$w3_g$=0(64Nx!|1xpm_%efK$q|ckzqnTHs!eCr*fq4$y>i*?%&gN=&ehN
zFSOqbOn6&E)9fSM^Z-P*O7wOF*nFfHXtDYNZMP+}hn`guXgwDh351;Yr<zVVwq5NV
z{v%UO;hIs(R|#4TQi$Er*P|15yN+n21~74$HOZ}zTv8OZ+7UeF0>vA9kqxv`VF7h|
zm|N2>4t}r_n_ai!4}^`%#{K*kBWk`ke)Q9!b8o6t`<FXodNf3P9S*L)UV7#rRAE*W
v4W*mQ$igbNT4K)CZryLflbD`<ET#lV<!D5A{w_JP;O1(xg>L<^E7JPkSi3am

literal 0
HcmV?d00001