denoising.py/code/Noise.py

100 lines
2.7 KiB
Python
Raw Permalink Normal View History

2015-09-09 13:52:20 +00:00
# ~*~ encoding: utf-8 ~*~ #
2015-09-14 14:59:35 +00:00
import random
import time
2015-09-14 14:59:35 +00:00
# import des libs internes
from utility import SaltAndPepper_Noise, Additive_Noise, Multiplicative_Noise, Color_Noise, Gaussian_Noise
from utility import Filter, Shape
2015-09-28 14:54:23 +00:00
from BMPFile import RGBPixel
class Noise:
# instanciation des classes qui sont dans utility/
def __init__(self):
self.SaltAndPepper = SaltAndPepper_Noise.SaltAndPepper_Noise();
self.Additive = Additive_Noise.Additive_Noise();
self.Multiplicative = Multiplicative_Noise.Multiplicative_Noise();
self.Gaussian = Gaussian_Noise.Gaussian_Noise();
self.Color = Color_Noise.Color_Noise();
self.Filter = Filter.Filter();
2015-09-28 14:54:23 +00:00
self.Shape = Shape.Shape();
# retourne le SNR d'une image bruitée par rapport à sa référence
# @param uRef matrice de pixels de l'image d'origine
# @param uNoisy matrice de pixels de l'image bruitée
#
# @return SNR retourne le SNR associé
#
def SNR(self, uRef, uNoisy, grayscale=True):
width = len( uRef[0] )
height = len( uRef )
# si les images n'ont pas les mêmes tailles
if len(uNoisy) != height or len(uNoisy[0]) != width:
print "Tailles différentes";
exit();
# on calcule la puissance du signal
P_uRef = self.Power( uRef ) if not grayscale else self.Power_grayscale( uRef )
# on calcule la matrice de = l'image bruitée - l'image d'origine = le bruit uniquement
uNoise = []
# pour chaque pixel
for lineNoisy, lineRef in zip(uNoisy, uRef):
uNoise.append( [] );
for pixNoisy, pixRef in zip(lineNoisy, lineRef):
2015-10-09 12:55:49 +00:00
pix = RGBPixel(0,0,0,pixRef.x, pixRef.y, pixRef.bpp );
pix.r = pixNoisy.r-pixRef.r;
pix.g = pixNoisy.r-pixRef.g;
pix.b = pixNoisy.r-pixRef.b;
uNoise[pixRef.y].append( pix);
2015-09-28 14:54:23 +00:00
# on calcule la puissance du bruit
P_uNoise = self.Power( uNoise ) if not grayscale else self.Power_grayscale( uNoise )
if P_uNoise == 0:
return 0;
else:
return 1.0 * P_uRef / P_uNoise;
# retourne la puissance d'une image (somme du carré de ses valeurs) teinte de gris
# @param pixelMap matrice de pixels de l'image de laquelle on veut la puissance
#
# @return power puissance de l'image
#
def Power_grayscale(self, pixelMap):
power = 0;
for line in pixelMap:
for pix in line:
power += pix.grayscale() ** 2
return power;
# retourne la puissance d'une image (somme du carré de ses valeurs) teinte de gris
# @param pixelMap matrice de pixels de l'image de laquelle on veut la puissance
#
# @return power puissance de l'image
#
def Power(self, pixelMap):
power = 0;
for line in pixelMap:
for pix in line:
power += pix.r ** 2 + pix.g ** 2 + pix.b ** 2
return power;