https://github.com/Chaogan-Yan/REST
Tip revision: 5f9de23b90b10fd5eca9e9bc1016bbde75d43ab6 authored by Chaogan-Yan on 16 June 2013, 16:06:36 UTC
Fixed a bug in temporal correlation of two groups of images in Image Calculator;The midline of VMHC results were set to zero.
Fixed a bug in temporal correlation of two groups of images in Image Calculator;The midline of VMHC results were set to zero.
Tip revision: 5f9de23
rest_Corr2FisherZ.m
function rest_Corr2FisherZ(ACorrMapFilename, AZScoreMapFilename, AMaskFilename)
% Fisher Z score transformation
%Usage:
% rest_Corr2FisherZ(ACorrMap, AZScoreMap, AMaskFilename)
% ACorrMapFilename, Input, Original Pearson Product Moment correlation map
% AZScoreMapFilename, Output, Fisher Z score map
% AMaskFilename, Input, Mask file
%
%-----------------------------------------------------------
% Copyright(c) 2007~2010
% State Key Laboratory of Cognitive Neuroscience and Learning in Beijing Normal University
% Written by Xiao-Wei Song
% http://resting-fmri.sourceforge.net
% Mail to Authors: <a href="Dawnwei.Song@gmail.com">Xiaowei Song</a>; <a href="ycg.yan@gmail.com">Chaogan Yan</a>
% Version=1.3;
% Release=20090321;
% Revised by YAN Chao-Gan, 080610. NIFTI compatible
% Last Revised by YAN Chao-Gan, 090321. Result data will be saved in the format 'single'.
%-----------------------------------------------------------
if ~(nargin==3) error(' Error using ==> rest_Corr2FisherZ. 3 arguments wanted.'); end
%Load the original ReHo map file
[BrainMap,VoxelSize, Header]=rest_readfile(ACorrMapFilename);
nDim1 = size(BrainMap,1); nDim2 = size(BrainMap,2); nDim3 = size(BrainMap,3);
BrainSize = [nDim1 nDim2 nDim3]; VoxelSize =VoxelSize';
mask=rest_loadmask(nDim1, nDim2, nDim3, AMaskFilename);
pos=find(mask);
BrainMap(pos) =0.5 * log((1 +BrainMap(pos))./(1- BrainMap(pos)));
rest_writefile(single(BrainMap),AZScoreMapFilename,BrainSize,VoxelSize,Header, 'single'); %Revised by YAN Chao-Gan, 090321. Result data will be stored in 'single' format. %'double');