用Matlab一次性计算不同栅格数据之间相关系数（逐月）

这里计算相关系数的代码是参考了这个文章：不同栅格数据之间的相关系数计算（输出为tif影像）_栅格数据相关系数-优快云博客

因为我需要计算的是a变量跟每个月温度和降水量之间的相关性，如果你每跑完一个月手动改一次代码就很麻烦，所以我就稍微扩展了一下，供大家参考，本人比较菜，用的方法也很机械，但保证好用！

clc;
clear;
tic;我习惯加上这个计时器，可以直观感受到代码效率

[a,R]=geotiffread("F:\paper1_data\paper1\snow_pheno\SCD\SCD_2002_modified.tif"); % 导入投影信息，放自己要分析的数据路径即可
info=geotiffinfo("F:\paper1_data\paper1\snow_pheno\SCD\SCD_2002_modified.tif"); % 同上
[m,n]=size(a);
for months=1:12 %一月到十二月
    month_code=sprintf('%02d',months);
    SCDsum=single(zeros(m*n,21)); % 我这里使用了single类型，因为后面要并行运算，可以缓解内存使用量
    for year=2002:2022 %年份开始和结束时间
        filename=strcat('F:\paper1_data\paper1\snow_pheno\SCD\','SCD_',int2str(year),'_modified.tif'); % 此处需要按你文件名修