One source of strange shading is bad VertexNormals. They have to be estimated from the data, and if they "wobble" from true normal, the surface's shading changes slightly.
On the test example, we "know" what the surface normal should be. We fix it to show that it fixes the shading. For real data, it seems unusual to know how calculate exactly the normals.
GraphicsRow[{
foo = ListContourPlot3D[Table[x^2 + y^2 + z^2,
{x, -1, 1, 0.02}, {y, -1, 1, 0.02}, {z, -1, 1, 0.02}],
Contours -> {1}, Mesh -> None],
foo /.
GraphicsComplex[p_, a___, HoldPattern[VertexNormals -> _], b___] :>
GraphicsComplex[p, a,
VertexNormals -> Transpose[Transpose@p - 50], b]
},
ImageSize -> 800]
These are easier ways to fix it, but (update) they just remove the vertex normals, which I didn't notice because of the fineness of the mesh:
DeleteCases[foo, HoldPattern[VertexNormals -> _], Infinity]
RegionPlot3D[DiscretizeGraphics[foo]]
(* each produces graphics like this: *)
A way to get the normals from this answer is not so simple, but closed surfaces even when their equations are unknown.
mesh = BoundaryDiscretizeGraphics[foo];
Graphics3D[
GraphicsComplex[
MeshCoordinates[mesh],
{First[
"DefaultPlotStyle" /. (Method /.
Charting`ResolvePlotTheme[Automatic, RegionPlot3D])],
EdgeForm[], Thread[MeshCells[mesh, 2], Polygon]},
VertexNormals -> Region`Mesh`MeshCellNormals[mesh, 0]
]
]
