Hi, I am now answering from my office.
I just found out your last post but, in the meantime, I have done a few things on my own.
Because of our security policy it takes me some time to transfer a file from my professional network on this open one, so no file given here.
I mainly worked on smoothing the original matrix NP for a better rendering.
Here his the code I made (the look-up table F is just an example among many matrices we can use ; I'm going to look at what ImageTools proposes)
Filtering := proc(M::Matrix(square))
local N, MF, F, KP, nx, ny:
KP := KroneckerProduct:
N := numelems(M):
MF := copy(M):
F := Matrix(3$2, [[1, 2, 1], [2, 4, 2], [1, 2, 1]]):
for nx from 2 to N-1 do
for ny from 2 to N-1 do
MF[nx, ny] := KP(M[nx-1..nx+1, ny-1..ny+1]; F):
NPF := Filtering(NP)
Used in plots:-surfdata, this filtered matrix NPF looks rather pretty (prtty doesn't mean correct of course).
I will continue to work on this idea with all the material you sent me.
See you later