A method comprising: determining an image of a subsurface geologic region of interest, the image representing at least in part ghost energy attributable to reflections caused by a reflecting interface; and deghosting the image, characterized in that the deghosting comprises processing data representing the image in a processor-based machine to: determine a first point spread function based at least in part on a first modeling and migration of at least one point scatterer for the region, the first modeling being based at least in part on the reflecting interface; determine a second point spread function based at least in part on a second modeling and migration of at least one point scatterer for the region, the second modeling not being based on the reflecting interface; apply the first and second point spread functions to the data representing the image to generate data indicative of a deghosted image, wherein applying the first point spread function to the data representing the image comprises applying the first point spread function as a blurring operator to generate data representing a reflectivity image and applying the second point spread function to the data representing the image comprises applying the second point spread function to the data representing the reflectivity image to generate the data indicative of a deghosted image; and output the deghosted image. determining an image of a subsurface geologic region of interest, the image representing at least in part ghost energy attributable to reflections caused by a reflecting interface; and deghosting the image, characterized in that the deghosting comprises processing data representing the image in a processor-based machine to: determine a first point spread function based at least in part on a first modeling and migration of at least one point scatterer for the region, the first modeling being based at least in part on the reflecting interface; determine a second point spread function based at least in part on a second modeling and migration of at least one point scatterer for the region, the second modeling not being based on the reflecting interface; apply the first and second point spread functions to the data representing the image to generate data indicative of a deghosted image, wherein applying the first point spread function to the data representing the image comprises applying the first point spread function as a blurring operator to generate data representing a reflectivity image and applying the second point spread function to the data representing the image comprises applying the second point spread function to the data representing the reflectivity image to generate the data indicative of a deghosted image; and output the deghosted image. determine a first point spread function based at least in part on a first modeling and migration of at least one point scatterer for the region, the first modeling being based at least in part on the reflecting interface; determine a second point spread function based at least in part on a second modeling and migration of at least one point scatterer for the region, the second modeling not being based on the reflecting interface; apply the first and second point spread functions to the data representing the image to generate data indicative of a deghosted image, wherein applying the first point spread function to the data representing the image comprises applying the first point spread function as a blurring operator to generate data representing a reflectivity image and applying the second point spread function to the data representing the image comprises applying the second point spread function to the data representing the reflectivity image to generate the data indicative of a deghosted image; and <EOS>