Author: Sebastian Korsak
Special thanks to Max Pierini and his repository: https://github.com/maxdevblock/j_johnson_M, and my prof. V. Costantoudis who gave me the code in MATLAB.
The method we followed is described on the paper: Numerical Simulation of 3D Rough Surfaces and Analysis of Interfacial Contact Characteristics, by Guoqing Yang, Baotong Li, Yang Wang and Jun Hong. Here we have some of the steps as are writen in this paper.
Just import the library,
import non_gaussian_surfaces as seb
and type,
seb.SAimage_fft_2(500,3,1,3,20,20,0.9,800.0,corr=True)