This paper presents efficient algorithm for interpretation of resistivity sounding data. The earth is modelled as a set of horizontal layers. Field observations are carried out at each sounding point using general configuration of four symmetrical electrodes : the Schlumberger and Wenner configurations are special cases of these. Sample values of the resistivity transform are computed from the field apparent resistivity data using an exponential approximation of Stefanescu kernel function. Parameters of the horizontally stratified earth model are obtained by an iterative interpretation of the average resistivity transform of an entire earth. The sample values of the average resistivity transform computed from field observation data are compared with those computed from the layer model parameters.