分享自己编的[x,y,t]三维数据截取区域的FUNCTION

我是柿子

function [area,arealonx,arealaty]=SUBAREA(x,lonx,laty,westb,eastb,northb,southb);
% Program for the area you need only
%  x is input data to divide, lonx is the zonal form of x,laty is the meridional form of x,westb,eastb,northb,southb are the boundaries of the area
%format :: x [lon lay time]
%area :: is the area you need


%找到东西南北四个边界所在位置

[Cw,Iw]=min(abs(lonx(:)-westb));
[Ce,Ie]=min(abs(lonx(:)-eastb));
[Cn,In]=min(abs(laty(:)-northb));
[Cs,Is]=min(abs(laty(:)-southb));

area(:,:,:)=x(Iw:Ie,Is:In,:);
arealonx(:,1)=lonx(Iw:Ie,1);
arealaty(:,1)=laty(Is:In,1);

只需要输入待截取变量，待截取变量的经纬度分布、所需区域的东、西、南、北四个边界即可

mofangbao

额  贴出来了为啥还要卖个贡献值

lqouc

mofangbao 发表于 2014-6-24 16:38
额  贴出来了为啥还要卖个贡献值

她比较呆，就是想把所有的发帖功能都试一下

新八不是废材男

灰常好，谢分享{:5_213:}

二爷名声在外

额……对楼主的创新意识表示不明觉厉……

topmad

lqouc 发表于 2014-6-24 16:43
她比较呆，就是想把所有的发帖功能都试一下

你们有情况....

二爷名声在外

topmad 发表于 2014-6-25 08:12
你们有情况....

老五一语惊醒梦中人啊……

lqouc

二爷名声在外 发表于 2014-6-25 08:14
老五一语惊醒梦中人啊……

怎么个情况？你俩一块睡的呀？我得和老五说说让他晚上别总把你吵醒。

二爷名声在外

lqouc 发表于 2014-6-25 08:35
怎么个情况？你俩一块睡的呀？我得和老五说说让他晚上别总把你吵醒。

歪楼了……

2006202788

不错不错  好东西