У меня есть данные температуры воздуха в сетке с пространственным разрешением 1,25 х 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