如何计算每个网格单元的面积？

Question

我有空气温度网格数据，空间分辨率为1.25 x 1.25度（lon-lat）。数据涵盖北半球，第一纬度为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

Answer 1

如果您不介意非matlab解决方案，最简单的方法可能是从命令行中简单地使用cdo：

cdo gridarea yourdata.nc gridarea.nc 

然后netcdf文件gridarea.nc将包含m ** 2中每个单元格的gridarea。然后，您可以将其读入您的matlab程序。

在ubuntu中可以安装cdo

sudo apt-get install cdo 

文档在这里：https://code.mpimet.mpg.de/projects/cdo/