@likanzhan, you may check other options in this other thread, in particular GMT.jl and DIVAnd.jl.
NB:
GMT’s splines in tension method also produces “surfaces that always pass through the data points”
1 Like