A comparative study was made of three methods of interpolation – inverse distance weighting (IDW), spline and ordinary kriging – after optimization of their characteristic parameters. These interpolation methods were used to represent the electric field levels for three emission frequencies (774kHz, 900kHz, and 1107kHz) and for the electrical stimulation quotient, QE, characteristic of complex electromagnetic environments. Measurements were made with a spectrum analyser in a village in the vicinity of medium-wave radio broadcasting antennas. The accuracy of the models was quantified by comparing their predictions with levels measured at the control points not used to generate the models. The results showed that optimizing the characteristic parameters of each interpolation method allows any of them to be used. However, the best results in terms of the regression coefficient between each model's predictions and the actual control point field measurements were for the IDW method. •Interpolation methods were used to represent the electric field levels in a village.•Temporal variation values were less than spatial variations detected in the area.•Inverse distance weighting, spline, and ordinary kriging were compared.•The accuracy was quantified comparing the predictions with levels at control points.•The best fits between predicted and measured values corresponded to the IDW model.