如何读取.mat文件中特定经纬度范围内的sst数据

Black_kiki

我有一个.mat文件，里面有sst, lon, lat这三个变量数据，我想读取出0-50N, 100-150N区域范围内的sst数据，并且将读取出的数据重新保存为还是这三个参数的.mat数据。
新人求教，非常感谢。

二爷名声在外

load
save

Lancelot

load
save     

Black_kiki

二爷名声在外 发表于 2017-11-27 13:42
load
save

我就是用的load,save
可是我只会提取经纬度数值，没提出来来经纬度对应的sst

Black_kiki

二爷名声在外 发表于 2017-11-27 13:42
load
save

我再研究研究，谢谢。

sunzy025

我有点笨，我感觉实在不行的话，就新建个变量，然后if 经纬度 then SSt赋给新变量。
你看呢？但是确实有点笨啊

Black_kiki

哎，我这个菜鸟
sst和qual_qual 是个8192*4096的矩阵，lon和lat 分别是8192*1 和4096*1的向量。
琢磨了一下午，写了个下面的代码。
生成的lon_new和lat_new 分别是:1137*1和1138*1 的向量。生成的sst_new是：129306*1的向量
现在还在琢磨怎么把sst_new 生成为 1137*1138的矩阵？？


    load (sstfilename, 'sst', 'qual_qual', 'lon', 'lat');
    lon_new = find(lon>100 & lon<150);
    lat_new = find(lat >0 & lat <50);
    [X, Y]=meshgrid(lon,lat);
    sst_new = sst( X>100 & X<150 & Y >0 & Y <50);
    qual_qual_new = qual_qual( X>100 & X<150 & Y >0 & Y <50);

Black_kiki

sunzy025 发表于 2017-11-27 17:38
我有点笨，我感觉实在不行的话，就新建个变量，然后if 经纬度 then SSt赋给新变量。
你看呢？但是确实有点 ...

哎，我这个菜鸟
sst和qual_qual 是个8192*4096的矩阵，lon和lat 分别是8192*1 和4096*1的向量。
琢磨了一下午，写了个下面的代码。
生成的lon_new和lat_new 分别是:1137*1和1138*1 的向量。生成的sst_new是：129306*1的向量
现在还在琢磨怎么把sst_new 生成为 1137*1138的矩阵？？


    load (sstfilename, 'sst', 'qual_qual', 'lon', 'lat');
    lon_new = find(lon>100 & lon<150);
    lat_new = find(lat >0 & lat <50);
    [X, Y]=meshgrid(lon,lat);
    sst_new = sst( X>100 & X<150 & Y >0 & Y <50);
    qual_qual_new = qual_qual( X>100 & X<150 & Y >0 & Y <50);

Black_kiki

Black_kiki 发表于 2017-11-27 18:21
哎，我这个菜鸟
sst和qual_qual 是个8192*4096的矩阵，lon和lat 分别是8192*1 和4096*1的向量。
琢磨了 ...

不知用reshape对还是用for循环对

sunzy025

Black_kiki 发表于 2017-11-27 18:21
哎，我这个菜鸟
sst和qual_qual 是个8192*4096的矩阵，lon和lat 分别是8192*1 和4096*1的向量。
琢磨了 ...

和你讲个超笨的办法，特别慢
[m,n]=size(sst);
for i=1:m
   for j=1:n
      if lon（j）>100 & lon（j）<150 & lat(i) >0 & lat(i)<50
           x=x+1;y=y+1
          new_sst(x,y)=sst(i,j);
      end
   end
end



什么i,j,m,n啥的，你看看是不是反了，我感觉差不多能出来