Probing the three-dimensional emission source and neutron skin via <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>π</mml:mi> <mml:mtext>−</mml:mtext> <mml:mi>π</mml:mi> </mml:mrow> </mml:math> correlations in heavy-ion c

The Richardson-Lucy algorithm is applied to reconstruct the three-dimensional source function of identical pions from their two-particle correlation functions. The algorithm's performance is first evaluated through simulations with Gaussian-type initial source functions. Its imaging quality and robustness are further demonstrated with experimental data from Au+Au collisions at 1.23 A GeV, collected by the HADES Collaboration. Additionally, using UrQMD simulations of Pb+Pb collisions at 1.5 A GeV, we show that the deblurred source functions exhibit sensitivity to the initial neutron skin thickness of the colliding nuclei. This highlights the potential of the Richardson-Lucy algorithm as a tool for probing the neutron density distribution in heavy nuclei.