因此,我正在尝试执行一个不适用于热带纬度的海洋学计算,因此,我需要通过将30N和30S之间的任何纬度变为零来删除它们。我试着这样做如下:
lat = th_lim(1):1/5.999:th_lim(2);
lat_no_tropics(30>lat>-30) = 0;
lat_no_tropics(isinf(lat_no_tropics)|isnan(lat_no_tropics)) = 0;
th_lim是我使用的数据中的纬度限制(-90.0833和90.00833(。选择1/5.999以给出与将在计算中使用的其他掩模相同的尺寸。然而,它的输出是lat_no_tropics(尽管它有正确的维度(只是一系列零,仅此而已。
如果有人能告诉我哪里出了问题,我将不胜感激
您可以通过以下方式将范围内的值设置为零:
% create data
lat = (rand(100,1)*100)-50;
% thresholds:
th_lower = -30;
th_upper = 30;
% set data within thresholds to 0
lat((th_upper > lat) & (lat > th_lower)) = 0;