matlab如何利用不规则矢量掩膜裁剪栅格数据？

马俊

各位老师，大佬好，
     我目前想用matlab提取两种文件相同地理位置的格网值，并且计算他们的相关系数。这两种文件大小不同，nc文件大，tif文件小，nc文件空间范围包含tif文件。 ps：tif文件本质上是用一个不规则的矢量掩膜在arcgis软件上裁剪大的tif文件（大小和nc文件一样）得到的，是我的真实研究区。 文件的示意图如下：

   我的思路是，先用ind2sub函数提取tif文件有效值的坐标，然后利用得到的坐标基于sub2ind函数得到tif文件和nc文件的索引值。  问题在于 这样操作得到的nc文件值的位置和tif文件的位置是否是一 一对应的？   如果不对应的话，有大神知道matlab如何根据不规则矢量或者栅格掩膜来裁剪，从而使得两种文件的空间范围一致？

我的代码如下：
filename1='D:\modis\result\MOD11C1\day'
filename5='C:\Users\Administrator\Desktop\2004_5km'
file1=dir(fullfile(filename1,'*.tif'))   %tif范围小
file5=dir(fullfile(filename5,'*.nc'))  %nc范围大
Input_length=length(file1)
lst_modis=[]
lst_clm=[]   %定义初始0矩阵
for i=1:Input_length
    [lst_day,R]=geotiffread(fullfile(filename1,file1(i).name))
    a=double(lst_day*0.02)  %尺度因子转换单位
    idx=find(a~=0)  %去除无效值
    [x,y]=ind2sub(size(a),idx)
    A=a(sub2ind(size(a),x,y))

    filename=fullfile(filename5,file5(i).name)
    skin=(ncread(filename,'TSKIN'))'   %需要转置
    SKIN=skin(sub2ind(size(skin),x,y))  %提取同样位置nc文件的温度信息
    lst_modis=[lst_modis;average]
    lst_clm=[lst_clm;SKIN]
end
number=length(lst_clm)
CC=corrcoef(lst_modis,lst_clm); %R 相关性系数

wjy_ecnu

马俊 发表于 2020-6-21 22:15
这个该如何操作.... 大神可以具体点吗 或者有对应的代码？

我也没有具体的代码  文件要是不大 你传给我 我看看 ？  我觉得比较感兴趣  毕竟我没试过 2935356853@qq.com

wjy_ecnu

用经纬度找到对应的格点不是挺好吗

马俊

wjy_ecnu 发表于 2020-6-21 21:29
用经纬度找到对应的格点不是挺好吗

这个该如何操作.... 大神可以具体点吗 或者有对应的代码？   

学习渣渣

论坛搜掩膜，有类似的叭好像