Use same color map for different images to more accurately compare them #68

New Issue

Closed

opened 2024-05-09 23:18:36 +02:00 by Benjamin_Loison · 9 comments

Benjamin_Loison commented 2024-05-09 23:18:36 +02:00

Owner

Copy Link

Copy Source

help(plt.imsave)

...
    vmin, vmax : float, optional
        *vmin* and *vmax* set the color scaling for the image by fixing the
        values that map to the colormap color limits. If either *vmin*
        or *vmax* is None, that limit is determined from the *arr*
        min/max value.
    cmap : str or `~matplotlib.colors.Colormap`, default: :rc:`image.cmap`
        A Colormap instance or registered colormap name. The colormap
        maps scalar data to colors. It is ignored for RGB(A) data.
...

```python help(plt.imsave) ``` ``` ... vmin, vmax : float, optional *vmin* and *vmax* set the color scaling for the image by fixing the values that map to the colormap color limits. If either *vmin* or *vmax* is None, that limit is determined from the *arr* min/max value. cmap : str or `~matplotlib.colors.Colormap`, default: :rc:`image.cmap` A Colormap instance or registered colormap name. The colormap maps scalar data to colors. It is ignored for RGB(A) data. ... ```

Benjamin_Loison added the enhancementlow priorityquick labels 2024-05-09 23:18:36 +02:00

Benjamin_Loison commented 2024-05-09 23:28:47 +02:00

Author

Owner

Copy Link

Copy Source

plt.imsave(PREFIX + 'multiple_colors.png', multipleColorsImage, vmin = multipleColorsImage.min(), vmax = multipleColorsImage.max())

is identical (according to diff) to:

plt.imsave(PREFIX + 'multiple_colors.png', multipleColorsImage)

```python plt.imsave(PREFIX + 'multiple_colors.png', multipleColorsImage, vmin = multipleColorsImage.min(), vmax = multipleColorsImage.max()) ``` is identical (according to `diff`) to: ```python plt.imsave(PREFIX + 'multiple_colors.png', multipleColorsImage) ```

Benjamin_Loison commented 2024-05-09 23:51:45 +02:00

Author

Owner

Copy Link

Copy Source

denoiser = bilateral np.min(multipleColorsImage) = -0.04890543970650529 np.max(multipleColorsImage) = 0.047219058067642594
denoiser = wavelet np.min(multipleColorsImage) = -0.05360332845280251 np.max(multipleColorsImage) = 0.03556257229735231
denoiser = mean np.min(multipleColorsImage) = -69.4864864864864 np.max(multipleColorsImage) = 36.67181467181476

If use the same color scale, then only mean is interesting as it enforces its own color scale and 2 other images are monocolor.

Mean:

Bilateral:

When remove mean denoiser:

Bilateral:

Wavelet:

``` denoiser = bilateral np.min(multipleColorsImage) = -0.04890543970650529 np.max(multipleColorsImage) = 0.047219058067642594 denoiser = wavelet np.min(multipleColorsImage) = -0.05360332845280251 np.max(multipleColorsImage) = 0.03556257229735231 denoiser = mean np.min(multipleColorsImage) = -69.4864864864864 np.max(multipleColorsImage) = 36.67181467181476 ``` If use the same color scale, then only `mean` is interesting as it enforces its own color scale and 2 other images are monocolor. Mean:  Bilateral:  When remove mean denoiser: Bilateral:  Wavelet: 

mean_rafael_230424_mean_multiple_colors.png

50 MiB

mean_rafael_230424_bilateral_multiple_colors.png

102 KiB

mean_rafael_230424_wavelet_multiple_colors.png

36 MiB

mean_rafael_230424_bilateral_multiple_colors.png

36 MiB

Benjamin_Loison referenced this issue from a commit 2024-05-10 00:03:13 +02:00

Fix #68: Use same color map for different images to more accurately compare them

Benjamin_Loison closed this issue 2024-05-10 00:03:13 +02:00

Benjamin_Loison commented 2024-05-10 00:27:53 +02:00

Author

Owner

Copy Link

Copy Source


prnu_estimate_distribution_of_values_for_rafael_230424_images.png

47 KiB

Benjamin_Loison commented 2024-05-10 02:25:41 +02:00

Author

Owner

Copy Link

Copy Source

Should make similar plot for RAISE flat-field to maybe understand why Rafael images are absorbing in #63.

Why have such an outlier for mean blue? Maybe because of borders? I would not say so as they do not behave differently as far as I understand. Maybe because of point at (644, 1200)? It seems to:

import numpy as np
import matplotlib.pyplot as plt

fileName = 'means/mean_flat-field_NEF_mean_blue'
npArray = np.load(f'{fileName}.npy')

plt.title('Distribution of PRNU estimation of RAISE flat-field raw blue channel')
plt.hist(npArray.flatten(), 'auto')
plt.show()

>>> index = -1; np.partition(npArray.flatten(), index)[index]
5042.010000000002

>>> index = -2; np.partition(npArray.flatten(), index)[index]
423.28999999999934

So the first and second maximum are 5042 and 423.29.

Should make similar plot for RAISE flat-field to maybe understand why Rafael images are absorbing in #63.    Why have such an outlier for mean blue? Maybe because of borders? I would not say so as they do not behave differently as far as I understand. Maybe because of point at (644, 1200)? It seems to:   ```python import numpy as np import matplotlib.pyplot as plt fileName = 'means/mean_flat-field_NEF_mean_blue' npArray = np.load(f'{fileName}.npy') plt.title('Distribution of PRNU estimation of RAISE flat-field raw blue channel') plt.hist(npArray.flatten(), 'auto') plt.show() ```   ```pycon >>> index = -1; np.partition(npArray.flatten(), index)[index] 5042.010000000002 >>> index = -2; np.partition(npArray.flatten(), index)[index] 423.28999999999934 ``` So the first and second maximum are 5042 and 423.29.

prnu_estimate_distribution_of_values_for_raise_flat-field_images.png

46 KiB

image.png

176 KiB

image.png

199 KiB

distribution_of_prnu_estimation_of_raise_flat-field_raw_blue_channel.png

24 KiB

distribution_of_prnu_estimation_of_raise_flat-field_raw_blue_channel.png

20 KiB

distribution_of_prnu_estimation_of_raise_flat-field_raw_blue_channel.png

49 KiB

Figure_1.png

46 KiB

Benjamin_Loison commented 2024-05-12 02:54:58 +02:00

Author

Owner

Copy Link

Copy Source

Could be interesting to just add a line with raw images.

Could be interesting to just add a line with raw images.

Benjamin_Loison commented 2024-05-12 03:10:39 +02:00

Author

Owner

Copy Link

Copy Source

>>> index = -100; np.partition(npArray.flatten(), index)[index]
27.28
>>> index = 100; np.partition(npArray.flatten(), index)[index]
-26.35

So width of 53.63 and compared to other denoisers width of about 0.1, the width is 563 times wider...

```pycon >>> index = -100; np.partition(npArray.flatten(), index)[index] 27.28 >>> index = 100; np.partition(npArray.flatten(), index)[index] -26.35 ``` So width of 53.63 and compared to other denoisers width of about 0.1, the width is 563 times wider...

Benjamin_Loison commented 2024-05-12 03:15:38 +02:00

Author

Owner

Copy Link

Copy Source

Could locate such maximums to potentially see visually other outliers.

>>> np.unravel_index(npArray.argmax(), npArray.shape)
(1200, 644)
>>> index = -0; np.where(npArray == np.partition(npArray.flatten(), index)[index])
(array([1202]), array([643]))
>>> index = -1; np.array(np.where(npArray == np.partition(npArray.flatten(), index)[index])[::-1]) * 2
array([[1288],
       [2400]])

They are not clearly visible as the first maximum, this is explainable as the intensity is 133 times less.

Could locate such maximums to potentially see visually other outliers. ```pycon >>> np.unravel_index(npArray.argmax(), npArray.shape) (1200, 644) >>> index = -0; np.where(npArray == np.partition(npArray.flatten(), index)[index]) (array([1202]), array([643])) >>> index = -1; np.array(np.where(npArray == np.partition(npArray.flatten(), index)[index])[::-1]) * 2 array([[1288], [2400]]) ``` They are not clearly visible as the first maximum, this is explainable as the intensity is 133 times less.

Benjamin_Loison commented 2024-05-12 03:45:34 +02:00

Author

Owner

Copy Link

Copy Source

Related to Benjamin_Loison/gimp/issues/29.

Related to [Benjamin_Loison/gimp/issues/29](https://gitlab.gnome.org/Benjamin_Loison/gimp/-/issues/29).

Benjamin_Loison commented 2024-05-13 01:43:33 +02:00

Author

Owner

Copy Link

Copy Source

Could compute proportion of pixels for mean outside largest other denoisers band.

But it may be a significant proportion exceeding not much, so should maybe just plot distribution of pixel values and zoom on the largest other denoisers band.

Could compute proportion of pixels for mean outside largest other denoisers band. But it may be a significant proportion exceeding not much, so should maybe just plot distribution of pixel values and zoom on the largest other denoisers band.

Sign in to join this conversation.

No Branch/Tag Specified

Branches Tags

master

prnu_written_as_such

evaluation

Labels

Clear labels

bug

Something is not working

Context-Adaptive Interpolator

duplicate

This issue or pull request already exists

enhancement

New feature

epic

A task that is going to take more than a day to complete.

help wanted

Need some help

high priority

invalid

Something is wrong

left for future work

low priority

medium

A task that should take less than a day to complete.

medium priority

meta

question

More information is needed

quick

A task that should take less than two hours to complete.

wontfix

This won't be fixed

No labels enhancement low priority quick

Milestone

No items

No Milestone

Projects

Clear projects

No projects

Assignees

Clear assignees

Benjamin_Loison (Benjamin Loison)

No Assignees

1 Participants

Notifications

Subscribe

Due Date

No due date set.

Dependencies

No dependencies set.

Reference: Benjamin_Loison/Robust_image_source_identification_on_modern_smartphones#68