Abstract In this work, we construct spherically symmetric solutions of Bardeen–boson stars within the framework of the Einstein–Klein–Gordon theory coupled to nonlinear electrodynamics (NLED) by employing numerical methods. Considering a thin accretion disk in the equatorial plane as the light source, we systematically investigate the optical appearance of boson stars using the ray-tracing method and the stereographic projection technique. Particular attention is paid to the influence of the initial scalar field <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:msub mml:miϕ</mml:mi> mml:mn0</mml:mn> </mml:msub> </mml:math> , the magnetic charge <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:miG</mml:mi> </mml:math> , and the observation angle <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:msub mml:miθ</mml:mi> mml:mio</mml:mi> </mml:msub> </mml:math> , on the image structure. As compact horizonless objects, boson stars produce optical images dominated by direct emission, while their morphology undergoes significant distortions as <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:msub mml:miθ</mml:mi> mml:mio</mml:mi> </mml:msub> </mml:math> increases. Higher values of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:msub mml:miϕ</mml:mi> mml:mn0</mml:mn> </mml:msub> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> mml:msub mml:miθ</mml:mi> mml:mio</mml:mi> </mml:msub> </mml:math> can give rise to lensed images. For all the parameters, the image center exhibits a brightness depression similar to the inner shadow of black holes, which poses challenges for distinguishing between boson stars and black holes. To address this, we propose two possible approaches: (i) combining the analysis of lensing bands with the effective potential to determine the existence of photon rings; and (ii) examining the polarization effects under synchrotron emission mechanisms. These results provide theoretical support for future high-resolution imaging efforts aimed at discriminating boson stars from black holes.

