binind2latlon文档
binind2latlon将被索引值的正弦geocoordinate网格。
关于正弦网格的更多信息在这里。
内容
语法
(纬度、经度)= binind2latlon (binind)(纬度、经度)= binind2latlon(…,“行”,NumberOfRows)
描述
(纬度、经度)= binind2latlon (binind)给本指标的地理坐标binind。
(纬度、经度)= binind2latlon(…,“行”,NumberOfRows)网格中的指定的行数。默认情况下,binind2latlon将尝试自动算出的行数,但谁也不能保证它会工作。如果你知道的行数,指定他们只是可以肯定的。根据决议表金宝搏官方网站在这里的行,常见的数字如下:
- 180行1学位网格
- 2160行5分钟网格
- 4320行网格1.5分钟
示例1:一个小网格
让我们重新创建简单的18行看到的例子在这里。首先绘制国界与我边界功能:
边界(“国家”,“k -”)轴紧包含“经”ylabel“纬度”
在美国宇航局的网站的例子包含412箱,纬度跨度18行:
%的本数字列表:垃圾箱= (1:412)';% geocoordinate对应每一本什么号码吗?(纬度、经度)= binind2latlon(垃圾箱,“行”,18);%标签每一本以其geocoordinate:文本(经度、纬度、num2str(箱),“颜色”,“r”,“水平的”,“中心”,“字形大小”6)
示例2:实际数据
这个示例使用一些月度排名数据在这里。开始读数据。虽然这些文件在.nc结束,ncread函数不能识别变量的数据类型我们想阴谋。所以我们必须使用h5read而不是:
一个= h5read (“A20021822002212.L3b_MO_CHL.nc”,“/ level-3_binned_data / chlor_a”);B = h5read (“A20021822002212.L3b_MO_CHL.nc”,“/ level-3_binned_data / BinList”);%为方便起见,把z和垃圾箱作为变量:z =双(A.sum_squared);垃圾箱= B.bin_num;
每一个的geocoordinatez值可以很容易地找到:
(纬度、经度)= binind2latlon(箱);
注意这里有超过1100万的测量,因此我们可能想看测量在一个小区域,说印度洋,在该地区(40 s-20n 35 e - 110 e)。这些坐标与对应的指标geomask:
印第安纳州= geomask(纬度,经度,20 [-40],[35 110]);
地理区域有1625033点,不幸的是他们不做出一个完美的网格,所以如果你想要你必须使用网格数据griddata或scatteredInterpolant。或者你可以做一个数据的散点图。我使用的旧版本的Matlab窒息当你给散射大量的数据点,所以下面我用Aslak Grinstedfastscatter函数。我也在使用cmoceancolormap (Thyng et al ., 2016)。
图fastscatter(经度(印第安纳州),lat(印第安纳州)、z(印第安纳州),“markersize”1)轴紧caxis ([0 4]) cb = colorbar;ylabel (cb、的平方之和)cmocean藻类边界(“国家”,“facecolor”,0.6 * (1 1 1))
作者信息
这个函数的一部分气候数据为Matlab工具箱。这个函数和支持文档写的乍得a德克萨斯金宝app大学奥斯汀分校的格林。