arXiv:2610.31026v1 Announce Type: new Abstract: We introduce a data-dependent compact multinode Shepard method for the approximation of scattered multivariate data. The proposed approach extends the compact multinode Shepard operator by incorporating nonlinear, data-dependent weights inspired by the weighted essentially non-oscillatory (WENO) methodology while preserving the partition of unity structure and the locality of the approximation. The construction is based on compactly supported multinode weight functions associated with local polynomial interpolants defined on unisolvent stencils. A smoothness indicator is employed to adapt the contribution of each stencil according to the local regularity of the underlying function, thereby reducing spurious oscillations near discontinuities while retaining the optimal approximation order in smooth regions. We establish the theoretical approximation properties of the proposed operator: it achieves the order of accuracy of its linear counterpart for sufficiently smooth functions and, under an explicit indicator-separation hypothesis, the same order is retained at every covered evaluation point near a jump discontinuity. Numerical experiments for local degrees up to four illustrate the performance of the method on both smooth and non-smooth data, on scattered as well as on gridded nodes.