Fractional calculus has increased in popularity in recent years, as the number of its applications in different fields has increased. Compared to the traditional operations in calculus (integration and differentiation) which are uniquely defined, the fractional-order operators have numerous definitions. Furthermore, a consensus on the most suitable definition for a given task is yet to be reached. Fractional operators are defined as continuous operators and their implementation requires a discretization step. In this article, we propose a discrete fractional Laplacian as a matrix operator. The proposed operator is real (non-complex) which makes it computationally efficient. The construction of the proposed fractional Laplacian utilizes the DCT transform avoiding the complexity associated with the discretization step which is typical in the constructions based on signal processing. We demonstrate the utility of the proposed operator on a number of data modeling and image processing tasks. INDEX TERMS Fractional-Laplacian, discrete operator, image-processing, trend-filtering, fractional calculus.