An extremely large-scale antenna array (ELAA) can potentially provide significantly increased spatial multiplexing and beamforming gains, as well as enhanced localization capability. While presenting new potential, its near-field propagation and spatial non-stationary properties also impose a great challenge on the receiver design. This paper focuses on the receiver design in an uplink orthogonal frequency-division multiplexing system with a planar ELAA deployed at the base station. To solve the challenging problem of 3D user localization and channel estimation with the planar ELAA, a space-frequency user localization and channel reconstruction (SF-ULCR) receiver is proposed. Under the Bayesian framework, an extended probability model is first established, to capture the channel structural information comprehensively, based on which an iterative receiver consisting of three modules is derived: element-wise line spectrum estimation (ELSE), distance parameter estimation (DPE), and near-field localization (NFL). In particular, the ELSE module handles the line spectrum relationships among multiple subcarriers in the frequency domain, the DPE module extracts and integrates the distance information from the line spectrum parameters, and the NFL module utilizes the messages of distances for user localization based on near-field spatial characteristics. Our numerical results demonstrate that the proposed SF-ULCR algorithm outperforms existing baselines in terms of channel estimation and localization performance, and that it approaches the Cramèr–Rao bound.