bin lat lon data and average

Question

akk 2019-8-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/475462-bin-lat-lon-data-and-average

回答： Kavya Vuriti 2019-8-14

Hi, I am trying to bin data by lat and long, and then index the bins. I can get this far:

minlat=min(lat);
minlong=min(long);
gridspacing_lat=6;
gridspacing_lon=5.2
latidx = 1 + floor((lat - minlat) ./ gridspacing_lat);
longidx = 1 + floor((long - minlong) ./ gridspacing_lon);
num=accumarray( [latidx(:), longidx(:)], 1) %finds the # measurements in each grid

This gives me a total of 6 bins with a different number of measurements in each bin. I also have temperature, month and year within this dataset, and I would like to find the monthly mean temperature within each year, and within each bin. I know I can find the means across the whole data set using

G=findgroups(year,month);
temp_mean=splitapply(@nanmean,temp,G)

However, I don't know how to index each bin number. For example, how do I find the monthly mean of each year within bin 1?

Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Kavya Vuriti 2019-8-14

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/475462-bin-lat-lon-data-and-average#answer_387556

From question, I understood that there is a dataset containing columns of latitude, longitude, year, month and temperature.

The accumarray function returns a matrix whose dimensions depend on the maximum value of “latidx” and “longidx”. For example, the element num(1,1) gives the number of occurrences of the row [1 1] in the matrix passed as input to accumarray function. As mentioned in the question, there are 6 bins which means there are 6 non-zero elements in “num” array. The indices of these non-zero elements gives the value of latitude and longitude in each bin. Now that we have the values of latitude and longitude of a bin, extract the rows corresponding to that latitude and longitude. Loop through these rows to obtain temperatures of particular month in a year.

For more understanding on accumarray function, refer to the following link:https://www.mathworks.com/matlabcentral/answers/299710-do-not-understand-accumarray-command-in-maltab