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]

[![enter image description here][1]][1]

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: *)

[![enter image description here][2]][2]

A way to get the normals from [this answer](https://mathematica.stackexchange.com/a/130356/4999) 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]
      ]
     ]

[![enter image description here][3]][3]


  [1]: https://i.sstatic.net/ut3sh.png
  [2]: https://i.sstatic.net/0IeJE.png
  [3]: https://i.sstatic.net/S3dNu.png