Andrea and Steve, Many thanks to for your suggestions. Íve made some tests with siena_flow2tal (it took me a while to figure out that the script needs some of the output from SIENA first) and now I get a threshold number and a *_FLOW_DIL and *_FLOW_TO_TAL images that I dońt know how to interpret. I understand that a posterior voxel-wise analysis is a way to go but I was thinking of a way to segment directly into ROÍs the IM1_HALFWAYTO_IM2_RENDER output from SIENA?? On the other hand, the solution of J Chen is really neat: measuring the atrophy in mm of cortical thickness change. That́s a much more practical result and easy to interpret for a clinician like me. The problem is how to get those numbers from SIENA output?? Does the script SIENA_FLOW2TAL include anything to do that? Many thanks in advance! Santiago ><((((º>`·.¸¸.·´¯`·.¸.·´¯`·...¸><((((º>`·.¸¸.·´¯`·.¸.·´¯`·...¸><((((º> `·.¸¸.·´¯`·.¸.·´¯`·...¸><((((º>`·.¸¸.·´¯`·.¸.·´¯`·...¸><((((º>`·.¸¸.·´ Santiago Reig, PhD. Gregorio Maranon Hospital, Universidad Complutense de Madrid Experimental Surgery and Medicine Lab Dr. Esquerdo 46, 28007 Madrid, SPAIN Fax# 34 91 426 5108, Ph# 34 91 426 5067 [log in to unmask] http://www.hggm.es/image