Fixed Noise generation functions

Signed-off-by: Nicola Vigano <nicola.vigano@esrf.fr>

zUtil_NoiseAnalysis/gtGenerateAdditiveWhiteNoise.m

@@ -6,7 +6,7 @@ function n = gtGenerateAdditiveWhiteNoise(dims, total_energy, data_type)
         n = cell(dims);
         summed_energy = 0;
         for ii = 1:numel(n)
-            n{ii} = rand(dims, data_type);
+            n{ii} = rand(dims{ii}, data_type);
             summed_energy = summed_energy + gtGetSignalEnergy(n{ii});
         end
         for ii = 1:numel(n)

zUtil_NoiseAnalysis/gtGenerateExtintionNoise.m

@@ -4,7 +4,7 @@ function n = gtGenerateExtintionNoise(blobs, attenuation_range, num, data_type)
     end
     blob_indxs = round((numel(blobs)-1) * rand(1, num) +1);
 
-    attenuation_vals = (rand(1, num) * (max(attenuation_range) - min(attenuation_range))) - min(attenuation_range);
+    attenuation_vals = (rand(1, num) * (max(attenuation_range) - min(attenuation_range))) + min(attenuation_range);
 
     max_radiuses = cellfun(@(x)(norm(size(x)) / 3), blobs);
     max_radiuses = reshape(max_radiuses, [1 numel(max_radiuses)]) .* rand(1, numel(max_radiuses)) ./ 2;
@@ -17,7 +17,7 @@ function n = gtGenerateExtintionNoise(blobs, attenuation_range, num, data_type)
         n{indx} = blobs{indx};
 
         att_indxes = getIndexes(size(n{indx}), max_radiuses(indx));
-        n{indx}(att_indxes) = n{indx}(att_indxes) .* attenuation_vals(ii);
+        n{indx}(att_indxes) = n{indx}(att_indxes) * (1 - attenuation_vals(ii));
     end
 
     % Now let's convert to noise