A bilateral filtering method includes decomposing an image patch by pixel intensity to form a stack of patches; computing spatial filtering response of each intensity; multiplying the spatial filtering response of each intensity by corresponding intensity, thereby resulting in multiplied spatial filtering response of each intensity; and summing up the multiplied spatial filtering responses of different intensities weighted with corresponding range weights.