statmorph calculates non-parametric morphological diagnostics of galaxy images (e.g., Gini-M_{20} and CAS statistics), and fits 2D Sérsic profiles. Given a background-subtracted image and a corresponding segmentation map indicating the source(s) of interest, statmorph calculates the following morphological statistics for each source: - Gini-M20 statistics; - Concentration, Asymmetry and Smoothness (CAS) statistics; - Multimode, Intensity and Deviation (MID) statistics; - outer asymmetry and shape asymmetry; - Sérsic index; and, - several shape and size measurements associated to the above statistics, such as ellipticity, Petrosian radius, and half-light radius, among others.