Abstract-This paper considers estimation of pixel-wise monotonic increasing (or decreasing) data from a time series of noisy blurred images. The motivation comes from estimation of mechanical structure damage that accumulates irreversibly over time. We formulate a Maximum A posteriory Probablity (MAP) estimation problem and find a solution by direct numerical optimization of a log-likelihood index. Spatial continuity of the damage is modeled using a Markov Random Field (MRF). The MRF prior includes the temporal monotonicity constraints. We tune the MRF prior, using a spatial frequency domain loopshaping technique to achieve a tradeoff between noise rejection and signal restoration properties of the estimate.The MAP optimization is a large-scale Quadratic Programming (QP) problem that could have more than a million of decision variables and constraints. We describe and implement an efficient interior-point method for solving such optimization problem. The method uses a preconditioned conjugate gradient method to compute the search step. The developed QP solver relies on the special structure of the problem and can solve the problems of this size in a few tens of minutes, on a PC.The application example in the paper describes structural damage images obtained using a Structural Health Monitoring (SHM) system. The damage signal is distorted by environmental temperature that varies for each acquired image in the series. The solution for the experimental data is demonstrated to provide an excellent estimate of the damage accumulation trend while rejecting the spatial and temporal noise.Index Terms-Optimal estimation, spatio-temporal filtering, regularization, interior-point method, isotonic regression, Markov Random Field.