2008ApJ...684..794K

Motivated by the limitations encountered with the commonly used direct reconstruction techniques of producing mass maps, we have developed a multiresolution maximum likelihood reconstruction method for producing two-dimensional mass maps using weak gravitational lensing data. To utilize all the shear information, we employ an iterative inverse method with a properly selected regularization coefficient which fits the deflection potential at the position of each galaxy. By producing mass maps with multiple resolutions in the different parts of the observed field, we can achieve a comparable level of signal to noise by increasing the resolution in regions of higher distortions or regions with an overdensity of background galaxies. In addition, we are able to better study the substructure of the massive clusters at a resolution which is not attainable in the rest of the observed field. We apply our method to the simulated data and to a 4 deg² field obtained by the Deep Lens Survey.