Implement The MATLAB 'fitdist' In Python

I am working on a image processing tool, and I am having some trouble finding a good substitute for matlab's fitdist to a Python soltuion. The matlab code works something like this

Solution 1:

Generally SciPy is useful in your case:

import scipy.stats as st
# KDE
st.gaussian_kde(data)
# Fit to specified distribution (Normal distribution in this example)
st.norm.fit(data)

Full reference is here: https://docs.scipy.org/doc/scipy/reference/stats.html