For initial simulations, we assumed a pyroelectric distribution g(x,z) with a "soft bi-morph" shape in the z direction and a sudden transition from g ≠ 0 for X < 0 to g = 0 for X > 0 (Fig. 2a). This distribution was substituted into (2) and the pyroelectric current I(X,ω) was calculated. Preliminary results show good agreement of the reconstructed profile (Fig. 2b) with the original distribution. Most importantly, the new approach reduces the smoothing in the lateral direction observed for n = 0 (Fig. 2c). The present approach deals with a two-dimensional pyroelectric distribution only. However, extension to three dimensions is straightforward, although some numerical optimizations will need to be made to keep the computational speed within practical limits for the increased number of neighbors.