To begin with, I would suggest also plotting the line-of-sight velocity and velocity dispersion maps for each inclination. In particular, you should see a clear anti-correlation between the velocity and the skewness, especially in the edge-on case. This is a signature of rotation-dominated kinematics and a flat or slowly changing rotation curve, due to the line-of-sight vetor intersecting (nearly) circular orbits at different angles outside the tangent radius:

Edited to add: The idea is that your line of sight (dashed green line) is tangent to the point A, and so the projected lines-of-sight velocity = the rotation velocity at that radius. The line of sight intersects stellar orbits at larger radii at increasingly large angles, so the projected velocity along the line of sight gets smaller and smaller (if the stellar velocity curve is flat or decreasing or even mildly increasing). Assuming a locally Gaussian distribution of stellar velocities at each radius, the contributions from radii $> A$ show up at smaller velocities (curves in the right-hand figure). Assuming that stellar density decreases at larger radii, the contributions from larger radii are smaller (lower heights in the right-hand figure). The net effect will be a shape with negative skew. For the other side of the galaxy, where velocities are negative (blueshifted), you will get the mirror situation, with a positively skewed LOSVD.
(One small suggestion would be to plot positive quantities with read and negativity quantities with blue -- the opposite of what you're currently doing -- since then positive velocities correspond directly to redshifts and negative velocities to blueshifts.)
The strong pixel-to-pixel variations in the outer part of the disk in your plots is almost certainly just due to low S/N (since the number density of particles is low at large radii).
Note that the standard approach in analyzing galaxy stellar kinematics is to model the line-of-sight velocity distribution using Gauss-Hermite polynomials, where the first- and second-order terms correspond to the mean velocity ($V$) and the velocity dispersion ($\sigma$), and the third- and fourth-order terms ($h_{3}$ and $h_{4}$) to the skewness and kurtosis, respectively. (See, e.g., van der Marel & Franx (1993) and Gerhard (1993).) I point this out mainly in case you want to compare your results with published analyses of other models and real galaxies.