Małgorzata Curyło, Eric Thrane, Paul D. Lasky, Dawson S. Gaynor

We present a practical implementation of a phase-coherent mapping technique for pulsar timing arrays that resolves the full complex polarisation state of the gravitational-wave sky as a function of direction and frequency. Unlike standard cross-correlation methods, this approach preserves the amplitude, phase, and polarisation of the signal in every sky pixel. The resulting maps constitute a compact, minimally processed summary of the data from which all subsequent analyses -- characterisation of a stochastic background, searches for anisotropy, and identification of individual sources -- can be derived within a single unified framework. Our implementation is fully compatible with established pulsar timing data analysis methods. We validate the framework through a series of realistic simulations with varying array configurations, noise properties, and signal types. We demonstrate robust recovery of source amplitudes and sky locations across different scenarios, and discuss the impact of polarisation leakage, noise, and direction-dependent array sensitivity on the recovery of astrophysical signals.