A method and apparatus for seismic analysis include obtaining an initial geophysical model and seismic data for a subsurface region; producing a subsurface image of the subsurface region with the seismic data and the geophysical model; generating a map of one or more geologic features of the subsurface region by automatically interpreting the subsurface image; and iteratively updating the geophysical model, subsurface image, and map of geologic features by: building an updated geophysical model based on the geophysical model of a prior iteration constrained by one or more geologic features from the prior iteration; imaging the seismic data with the updated geophysical model to produce an updated subsurface image; and automatically interpreting the updated subsurface image to generate an updated map of geologic features. The method and apparatus may also include post-stack migration, pre-stack time migration, pre-stack depth migration, reverse-time migration, gradient-based tomography, and/or gradient-based inversion methods.