我有空间分辨率为 1.25 x 1.25 度（经纬度）的空气温度网格数据。数据覆盖北半球，第一纬度为90度。
我需要计算每个网格单元的面积，我的方法是从纬度带计算。
另一种方法是直接从 netCDF 文件中读取每个网格单元的区域。这在 MATLAB 中可能吗？
这是我尝试过的代码；我不确定这是否正确。
i=1:72;
j=2:73;
R=6371; % Earth's radius in km^2
deltalon=1.25;
area=(pi/180)*R^2*((sind(latitude(i)) - sind(latitude(j)))*deltalon); % Area of each grid cell in km^2