In today’s modern wide-field galaxy surveys, there is the necessity for parametric surface brightness decomposition codes characterised by accuracy, small degree of user intervention, and high degree of parallelisation. We try to address this necessity by introducing MORPHOFIT, a highly parallelisable python package for the estimate of galaxy structural parameters. The package makes use of wide-spread and reliable codes, namely, SEXTRACTOR and GALFIT. It has been optimised and tested in both low-density and crowded environments, where blending and diffuse light makes the structural parameters estimate particularly challenging. MORPHOFIT allows the user to fit multiple surface brightness components to each individual galaxy, among those currently implemented in the code. Using simulated images of single Sérsic and bulge plus disk galaxy light profiles with different bulge-to-total luminosity (B/T) ratios, we show that MORPHOFIT is able to recover the input structural parameters of the simulated galaxies with good accuracy. We also compare its estimates against existing literature studies, finding consistency within the errors. We use the package in a companion paper to measure the structural parameters of cluster galaxies in order to study the wavelength dependence of the Kormendy relation of early-type galaxies. The package is available on github¹ and on the Pypi server².