function [C,p]=pattern_correlation(map1,map2)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  PATTERN CORRELATION
%  [C,p]=pattern_correlation(map1,map2)
% 
%    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[nx,ny]=size(map2); % x or y are interchangeable
X1=reshape(map1,[nx*ny,1]);
X2=reshape(map2,[nx*ny,1]);
[Cm,Pm]=corrcoef(X1,X2);
C=Cm(1,2); % either off-diagonal element is fine
p=Pm(1,2);  % P-value for that correlation