A simple and robust imaging technique for two-dimensional magnetotelluric interpretations has been developed following the well known Niblett-Bostick transformation for one-dimensional profiles. The algorithm processes series and parallel magnetotelluric impedances and their analytical influence functions using a regularized Hopfield artificial neural network. The adaptive, weighted average approximation preserves part of the nonlinearity of the original problem, yet no initial model in the usual sense is required for the recovery of the model; rather, the built-in relationship between model and data automatically and concurrently considers many half spaces whose electrical conductivities vary according to the data. The use of series and parallel impedances, a self-contained pair of invariants of the impedance tensor, avoids the need to decide on best angles of rotation for identifying TE and TM modes. Field data from a given profile can thus be fed directly into the algorithm without much processing. The solutions offered by the regularized Hopfield neural network correspond to spatial averages computed through rectangular windows that can be chosen at will. Applications of the algorithm to simple synthetic models and to the standard COPROD2 data set illustrate the performance of the approximation.