Detailed structural decomposition of galaxy images

We present a two-dimensional fitting algorithm (GALFIT) designed to extract structural components from galaxy images, with emphasis on closely modeling light pro files of spatially well-resolved, nearby galaxies observed with the Hubble Space Telescope. Our algorithm improves on previous techniques in two areas: by being able to simultaneously fit a galaxy with an arbitrary number of components and with optimization in computation speed, suited for working on large galaxy images. We use two-dimensional models such as the Nuker law, the Sersic (de Vaucouleurs) profile, an exponential disk, and Gaussian or Moffat functions. The azimuthal shapes are generalized ellipses that can fit disky and boxy components. Some potential applications of our program include: standard modeling of global galaxy profiles; extracting bars, stellar disks, double nuclei, and compact nuclear sources; and measuring absolute dust extinction or surface brightness fluctuations after removing the galaxy model. When examined in detail, we find that even simple looking galaxies generally require at least three components to be modeled accurately, rather than the one or two components more often employed. Many galaxies with complex isophotes, ellipticity changes, and position angle twists can be modeled accurately in two dimensions. We illustrate this by way of 11 case studies, which include regular and barred spiral galaxies, highly disky lenticular galaxies, and elliptical galaxies displaying various levels of complexities. A useful extension of this algorithm is to accurately extract nuclear point sources in galaxies. We compare two-dimensional and one-dimensional extraction techniques on simulated images of galaxies having nuclear slopes with different degrees of cuspiness, and we then illustrate the application of the program to several examples of nearby galaxies with weak nuclei.