约化密度矩阵（RDM）在量子纠缠和测量中起着关键作用，因为它能够提取几乎所有与约化自由度相关的物理量。然而，受限于环境的自由度，整个系统的规模往往受到限制，更不用说子系统了。为了解决这一挑战，我们提出了一种技术壁垒较低的量子蒙特卡洛方案，能够精确提取RDM。为了展示该方法的能力，我们展示了约化纠缠谱（ES）的精细结构，即RDM的对数特征值。我们清晰地展示了具有长纠缠边界的1维梯子模型以及具有塔状态的2维海森堡模型的ES。此外，我们提出了从采样的RDM数据高效恢复算符形式的纠缠哈密顿量的方法。我们的模拟结果利用前所未有的大系统规模，建立了基于RDM确定纠缠量（如ES）的实用计算框架，特别是在环境具有大量自由度的场景中。

The reduced density matrix (RDM) plays a key role in quantum entanglement and measurement, as it allows the extraction of almost all physical quantities related to the reduced degrees of freedom. However, restricted by the degrees of freedom in the environment, the total system size is often limited, let alone the subsystem. To address this challenge, we propose a quantum Monte Carlo scheme with a low technical barrier, enabling precise extraction of the RDM. To demonstrate the power of the method, we present the fine levels of the entanglement spectrum (ES), which is the logarithmic eigenvalues of the RDM. We clearly show the ES for a $1$D ladder with a long entangled boundary, and that for the $2$D Heisenberg model with a tower of states. Furthermore, we put forward an efficient way to restore the entanglement Hamiltonian in operator-form from the sampled RDM data. Our simulation results, utilizing unprecedentedly large system sizes, establish a practical computational framework for determining entanglement quantities based on the RDM, such as the ES, particularly in scenarios where the environment has a huge number of degrees of freedom.