A system, device and a method to perform bilateral filtering using linear convolution by way of an FFT or a recursive sequence method. Proper selection of functions for the photometric and spatial components of a bilateral filter may reduce the computational cost of the bilateral filter while preserving the bilateral filter de-noising and edge detection capabilities. Such functions may reduce the computational cost of a bilateral filter to substantially O(1).