1. A method of reducing noise in an image comprising: decomposing the image to generate wavelet coefficients at different scales; modifying the wavelet coefficients based on the energy of the wavelet coefficients at said different scales; and reconstructing the image based on the modified wavelet coefficients; wherein during the decomposing, said image is decomposed into wavelet detail spaces at different scales, each wavelet detail space comprising a plurality of wavelet detail subbands and each wavelet detail subband comprising a matrix of wavelet coefficients; and wherein during the modifying, the feature energy of the wavelet coefficients is determined according to: FeatureEnergy j ⁡ ( p ) = { NISE j ⁡ ( p ) + NeighborEnergy j ⁡ ( p ) j < J 0 SSE j ⁡ ( p ) + NeighborEnergy j ⁡ ( p ) j = J 0 where NISE j (p) is the normalized inter-scale energy, SSE j (p) is the single scale energy, NeighborEnergy j (p) is the weighted average single scale energy, j is the scale and p is location (x,y) in each wavelet detail subband.