diff --git a/datasets/raise/extract_noise.py b/datasets/raise/extract_noise.py
index 2dfb0aa..c583f48 100755
--- a/datasets/raise/extract_noise.py
+++ b/datasets/raise/extract_noise.py
@@ -4,9 +4,8 @@ import numpy as np
 import os
 from tqdm import tqdm
 import csv
-from utils import Color, denoise, iterativeMean, isARawImage, escapeFilePath, getColorChannel, saveNpArray, rescaleRawImageForDenoiser, updateExtremes
+from utils import Color, denoise, iterativeMean, isARawImage, escapeFilePath, getColorChannel, saveNpArray, rescaleRawImageForDenoiser, updateExtremes, getColorMeans
 import matplotlib.pyplot as plt
-from scipy.ndimage import gaussian_filter
 
 imagesFolderPath = 'rafael/arw'
 imagesFolderPathFileName = escapeFilePath(imagesFolderPath)
@@ -91,17 +90,12 @@ if (minColor is None or maxColor is None) and denoiser != 'mean':
     print(f'{maxColor=}')
 
 if denoiser == 'mean':
-    colorMeans = {}
-    for color in colors:
-        colorIterativeMean = iterativeMean()
-        for imageFileName in tqdm(imagesFileNames, f'Computing mean of {color} colored images'):
-            imageNpArray = getImageNpArray(imageFileName, False, color)
-            imageNpArray = gaussian_filter(imageNpArray, sigma = 5)
-            colorIterativeMean.add(imageNpArray)
-        colorMeans[color] = colorIterativeMean.mean
+    colorMeans = getColorMeans(imagesFileNames, colors)
+    for color in Color:
+        colorMeans[color] = colorMeans[color]
         fileName = f'mean_{imagesFolderPathFileName}_{color}'
         # Then use `merge_single_color_channel_images_according_to_bayer_filter.py` to consider all color channels, instead of saving this single color channel as an image.
-        saveNpArray(fileName, colorIterativeMean.mean)
+        saveNpArray(fileName, colorMeans[color])
 
 for color in colors:
     estimatedPrnuIterativeMean = iterativeMean()
diff --git a/datasets/raise/utils.py b/datasets/raise/utils.py
index 704411e..7314baa 100644
--- a/datasets/raise/utils.py
+++ b/datasets/raise/utils.py
@@ -7,6 +7,7 @@ from PIL import Image
 from skimage import img_as_float
 from datetime import datetime
 import builtins as __builtin__
+from scipy.ndimage import gaussian_filter
 
 class Color(Enum):
     RED = auto()
@@ -130,3 +131,14 @@ def updateExtremes(imageNpArray, minColor, maxColor):
 
 def print(*toPrint):
     __builtin__.print(datetime.now(), *toPrint)
+
+def getColorMeans(imagesFileNames, colors = Color):
+    colorMeans = {}
+    for color in colors:
+        colorIterativeMean = iterativeMean()
+        for imageFileName in tqdm(imagesFileNames, f'Computing mean of {color} colored images'):
+            imageNpArray = getImageNpArray(imageFileName, False, color)
+            imageNpArray = gaussian_filter(imageNpArray, sigma = 5)
+            colorIterativeMean.add(imageNpArray)
+        colorMeans[color] = colorIterativeMean.mean
+    return colorMeans