Average multiple cells in the same lat and lon

6 views (last 30 days)
Hello,
I have a 3422x5 (double) numeric matrix A = [Latitude Longitude Biomass lower_boundary upper_boundary]
the lower and upper boundary represent depth
I want to be able to average the biomass across any bins that are in the same exact latitude and longitude. There are several measurements that are in the same lat and lon, so I want to be able to get just one data point for each lat and lon. I also want to be able to calculate the standard deviation of each of those calculations in a different column to keep track of the statistics. Along with this, I would like to keep the lower and upper boundary columns intact with respect to the new, averaged dataset. The duplicated bins all have the same correlated upper and lower boundary
What would be the most efficient way of accomplishing this? I am attaching the matlab file for better reference as to what I am talking about
Any help is appreciated.
Thanks,
Melissa

Accepted Answer

Sean de Wolski
Sean de Wolski on 15 Dec 2014
Edited: Sean de Wolski on 15 Dec 2014
% Remove nans
idxKeep = ~any(isnan(A(:,1:3)),2);
Akeep = A(idxKeep,:);
% Unique combinations of lat/lon and where they are
[ulatlon,~,idx] = unique(Akeep(:,1:2),'rows');
% Mean of each combination
meanBm = accumarray(idx,Akeep(:,3),[],@mean);
% Combine unique lat/lon combos and meanBm into table
T = table(ulatlon,meanBm,'VariableNames',{'LatLon','MeanBioMass'});
% Show Results
disp(T)
  5 Comments
Sean de Wolski
Sean de Wolski on 7 Jan 2015
I guess I'm not clear. Isn't that what idx is above? idx, right now is telling you which row each unique lat/lon is?
So for a given array
x = [1;4;4]
[ux,idxf,idxr] = unique(x)
ux =
1
4
idxf =
1
2
idxr =
1
2
2
So idxr (idx in my original) is telling me where each row came from. idxf, is the first match for the unique values. I think all you need to do is group by idxr, but I'm not completely clear.
Could you provide a small example of what you want with five rows, and just lower boundary?
Melissa
Melissa on 8 Jan 2015
Ah, yes! So what I was looking for in your example is idxf! That is the one that you omitted in your original code with ~. So I switched it to
[ulatlon,idxf,idx] = unique(Akeep(:,1:2),'rows');
Now I can use idxf to do this
lowerKeep = lower_boundary(idxf);
upperKeep = upper_boundary(idxf);
Which now allows me to make
A = [Lat Lon meanBm Lower Upper]
Thank you so much for your quick responses!

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!