Beyond diagonal reconfigurable intelligent surface (BD-RIS)s enhance wave manipulation through inter-element couplings but pose significant channel estimation challenges due to cascaded channels and block-Kronecker structures. This paper proposes a compressive sensing framework exploiting sparse Tucker decomposition of the measurement tensor and the Kronecker rank-one structure of channel components. Two algorithms are developed: Sparse Tensor Orthogonal Recovery Method (STORM), which uses orthogonal matching pursuit (OMP) for greedy support recovery, and Sparse Tensor subspace- Aided Recovery (STAR), which leverages subspace-based projection for enhanced noise robustness. Both perform joint sparse support identification, followed by a Kronecker rank-one factorization via singular value decomposition (SVD) to recover the channel parameters. Simulations show that STAR achieves oracle-assisted least squares (LS) performance at moderate-to-high signal-to-noise ratio (SNR) with significantly fewer measurements than baseline methods, enabling practical BD-RIS deployment in next-generation millimeter wave (mmWave)/sub-terahertz (sub-THz) networks.