本发明提供一种基于多变量核密度欧拉解概率密度成像方法，首先通过测量或位场变换获得重力矢量及重力梯度张量数据，构建三维重力/重力梯度张量欧拉反褶积方程，逐一移动滑动窗口，获得一系列含异常源空间位置[x，y，z，N]欧拉解集，针对常规欧拉反褶积仅通过欧拉反褶积解的构造指数的“色谱”难于评价、分离异常，而不综合考虑欧拉解的空间分布特性及构造指数的相关性。采用三维核密度估计对欧拉反褶积解集综合分析，从而实现了欧拉反褶积对相邻异常有效分离、快速定位。