B-algorithm was applied to process GPR data. The algorithm detects heterogeneities correlated with region boundaries and “regime” shifts in GPR data paths, and plots their distribution. Implementation of the algorithm in the MATLAB system simplifies and accelerates the calculations, bringing the estimates closer to the optimum solutions. The algorithm, tested in application to radargrams from a model medium and a winter lake, proved to be sensitive to “regime” shifts and capable of processing data automatically.