The extended use of ray maps in the phase-space representation of optical fields in any state of spatial coherence has promoted the search of efficient approaches for their computation. In this paper, those maps are calculated via a matrix-based approach that i) allows accounting for complex degrees of spatial coherence of any shape and value, ii) models transmittances with real dimensions and not just point objects, iii) accounts for Fresnel phases in the modeled fields, and iv) consumes reasonably short computation time (few seconds for relatively big matrices). Examples illustrate the mentioned features.