fock[n_, q_, p_] := (-1)^n /
   Pi Exp[-q ^2 - p^2] LaguerreL[n, 2 q^2 + 2 p^2];
rq = 4; rp = 4; t = 0.003; FontS = 20;
For[n = 0, n <= 4, n++, plot = Show[Plot3D[fock[n, q, p],
    {q, -rq, rq}, {p, -rp, rp}, ImageSize -> 800,
    Mesh -> {Range[-Floor[rq], Floor[rq] ], 
      Range[-Floor[rp], Floor[rp] ], Range[-1, 1, 1/4]/Pi},
    MeshFunctions -> {#1 &, #2 &, #3 &},
    MeshStyle -> {Directive[Black, Thickness[t] ], 
      Directive[Black, Thickness[t] ], 
      Directive[White, Thickness[t] ]},
    PlotRange -> {-1/Pi, 1/Pi}, PlotPoints -> 81, MaxRecursion -> 4, 
    Method -> {Refinement -> {ControlValue -> 0.02} },
    PerformanceGoal -> "Quality",
    PlotStyle -> Opacity[0.85], Lighting -> "Classic",
    ColorFunction -> ({RGBColor[1, 1, 0.75], Glow[GrayLevel[0.06] ], 
        Specularity[0.5, 60]} &),
    Axes -> False, Boxed -> False,
    ViewPoint -> 
     FromSphericalCoordinates[{Sqrt[229/20], Pi/3, -0.64 Pi}] ],
   Graphics3D[{Thickness -> t, Black, 
     Line[{ {-rq, rp, 0}, {-rq, -rp, 0}, {rq, -rp, 0} }]}],
   Graphics3D[{Thickness -> t, Black, 
     Line[{ {-rq, rp, -1/Pi}, {-rq, rp, 1/Pi} }]}],
   (* q ticks *)
   Sequence @@ 
    Table[Graphics3D[{Thickness -> t, Black, 
       Line[{ {x, -rp, 0}, {x, -0.2 - rp, 0} }]}], {x, -Floor[rq], 
      Floor[rq]}],
   (* p ticks *)
   Sequence @@ 
    Table[Graphics3D[{Thickness -> t, Black, 
       Line[{ {-rq, y, 0}, {-rq - 0.2, y, 0} }]}], {y, -Floor[rp], 
      Floor[rp - 1/2]}],
   (*W ticks *)
   Sequence @@ 
    Table[Graphics3D[{Thickness -> t, Black, 
       Line[{ {-rq, rp, z/(2 Pi)}, {-rq - 0.2, rp, 
          z/(2 Pi)} }]}], {z, -2, 2}],
   (* box *)
   (*Graphics3D[{Thickness\[Rule]t/2,Gray,Line[{ {rq,-rp,1/Pi},{rq,rp,
   1/Pi},{-rq,rp,1/Pi} }],Line[{ {rq,rp,0},{rq,rp,1/Pi} }],
   Line[{ {rq,-rp,0},{rq,-rp,1/Pi} }]}],*)
   (* axes labels *)
   Graphics3D[Text[Style["q", FontS, Black], {0, -rp*1.15, -0.07}] ],
   Graphics3D[Text[Style["p", FontS, Black], {-rq*1.15, 0, -0.07}] ],
   Graphics3D[
    Text[Style["W", FontS, Black], {-rq*0.93, rp*0.93, 0.8/Pi}] ],
   Sequence @@ 
    Table[Graphics3D[{Text[
        Style[TextString[x], FontS, 
         Black], {x, -rp - 0.10 Max[rq, rp], 0}, {0, 
         1}]}], {x, -Floor[rq], Floor[rq]}],
   Sequence @@ 
    Table[Graphics3D[{Text[
        Style[TextString[y], FontS, Black], {-rq - 0.10 Max[rq, rp], 
         y, 0}, {0, 1}]}], {y, -Floor[rp], Floor[rp - 1/2]}],
   Sequence @@ 
    Table[Graphics3D[{Text[
        Style[If[z == 0, "0", ToString[z/(2 Pi), TraditionalForm] ], 
         FontS, Black], {-rq - 0.3, rp, z/2/Pi}, {1, 0}]}], {z, -2, 
      2}],
   BoxRatios -> {Automatic, Automatic, 6}, PlotRange -> All
   ];
 trim = { {0., .25}, {.95, .91} };
 imgname = "Wignerfunction_fock_" <> TextString[n] <> ".png";
 Export[imgname, 
  ImageResize[
   ImageTrim[Image[plot, ImageResolution -> 400], trim, 
    DataRange -> { {0, 1}, {0, 1} }], 2000, Resampling -> "Linear"] ];
 ]