Fit multiple lines un function

Dear community,

I have the following binary image (blue scatter) where I would like to detect the two clear branching points. The objective is to fit three different functions for each of the branches. Any ideas?

Best.