Currently based on the local coordinate system and the surface normal. This has the following issues:
- There are by construction singular points in the map if surface with non-zero Euler number. For the sphere, these have been removed, but if actual defects are in this region, they get removed too.
- For this to work on other surfaces, the normal and the coordinate system need to match, which is not quite true for all of our surfaces.