In the attached images, the "Data.jpg" image is a spectrogram for a type-II radio burst with frequency on the y-axis and time on the x-axis and the intensity of the signal is labeled with colors. There is also a band-splitting in the signal.
I need to trace the maximum intensities of the signal as shown in "DATA_fitln.jpg" image for both bands (red line & black line) and store the values of (f, t) for each fitting line in a matrix. How can I do that in Matlab or IDL?
I have no experience with Matlab. Probably you are searching for
something like center of mass. The maximum is not unique enough.
https://en.wikipedia.org/wiki/Center_of_mass
Maybe in combination with clustering before.
https://en.wikipedia.org/wiki/K-means_clustering
Regards,
Joachim
Couldn't you slice the figure by column in Matlab, fit each column by two 1d guassians to find peak locations, and overlay the result on the original picture?
Could you please explain in details or send some resources describing this method?
Thank you
There is an iteration fitting method used for type III: Article A type III radio burst automatic analysis system and statist...
But you can modify the parameters (the slope range) to let it fit the type II.
This algorithm has an open sourced implementation:
https://github.com/Pjer-zhang/TypeIIIRadioBurstRecognition
- Badges
- Science topic