Extract sub-matrix compute correlation coefficient 2D case iteratively

2 次查看（过去 30 天）

Joshua 2020-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/506221-extract-sub-matrix-compute-correlation-coefficient-2d-case-iteratively

评论： darova 2020-2-18

Hello friends, I will like to extract 2x2 sub-matrices from a large matrix and calculate correlation coefficient for each of the 2x2. Loop the procedure over the main matrix and store the correlation coefficient in output file. the for loop gave some error, I will appreciate guide to fix the code. Thank you.

% infile
file1='d.dat';   %
file2='b.dat';   %
% Outfile
correlout='correl.dat';  %
numcols=20;           % 
numrows=20;
%open file
fid1=fopen(file1);
f1=fscanf(fid1,'%f',[numcols numrows]);    
fclose(fid1);
fid2=fopen(file2);
f2=fscanf(fid2,'%f',[numcols numrows]);  
fclose(fid2);
%reshape
Refid1=reshape(fid1,[numcols,numrows]);
Refid2=reshape(fid2,[numcols,numrows]);
%pre-allocate
ccoef=zeros(numrows,numcols);
cRes=zeros(numrows,numcols);
subf1(i)=zeros();
subf2(i)=zeros();
%loop over 
for i=1:numcols
    for j=1:numrows
      subf1(i) = Refid1(i:i+1, j:j+1);
        subf2(i) = Refid2(i:i+1, j:j+1);
        ccoef(i)=corrcoef(subf1(i),subf2(i));
      cRes(i)=ccoef(1,2);
    end
end
%outfile
correlout=cRes(i);
%write to file
fid=fopen(correlout, 'w');
fprintf(fid,'%6.2f\n',correlout);
fclose(fid);