I attached a video that shows a density evolving with distinct peaks that broaden with increasing temperature. These are density plots of large sets of unit normal vectors: to get here, I have used the calcDensity command and used the standard plot function on the resulting s2FunHarmonic object. I have two questions:
1. I would like to automatically count the number of unique peaks (3 for this simple example). I am trying [~,localMax] = max(density,'numlocal',100) and am enforcing some cutoff below which I don't consider a peak to be unique. In practice, the algorithm usually gives more than 3 peaks with the noisy data. Do you have any suggestions for peak finding here?
2. I would also like to quantify the broadening of the peaks. I see a variety of functions that can be used with a S2FunHarmonic object but do not know what would be appropriate to use here. Any guidance is much appreciated.