

Based on the ETF

Based on gradients



The next step is to control the initial condition and the boundary condition and to use different functions or different diffusion to model the surface.
The original image
3000 iterations, 120 seconds
2000 iterations
After 500 iterations(32 seconds), pixels with height values between 0.1 to 0.000001 are shown in red. It means the narrow areas around the hair have a lot almost-zero heights.
500 iterations