Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#62 2024-05-13 16:42:56 +02:00
Optimize execution speed

Thought about using benchmark_raw_images_loading.py:

#!/usr/bin/env python

from tqdm import tqdm

IMAGES_CAMERAS_FOLDER = {
    'RAISE': 'flat-field/nef',
    'Rafael…
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#72 2024-05-13 16:41:17 +02:00
Correctly implement iterative mean for camera attribution
class iterativeMean:
    mean = None
    numberOfElementsInMean = 0

    def add(self, element):
        if self.mean is None:
            self.mean = element
        else:
    
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#62 2024-05-13 16:37:48 +02:00
Optimize execution speed

It does not seem that can only load part of a raw image with rawpy, if even possible at all.

In getRawColorChannel:

rawImageVisible = raw.raw_image_visible.copy()[:100,
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#72 2024-05-13 16:21:34 +02:00
Correctly implement iterative mean for camera attribution

50 images for both cameras take about 22 GB of memory.

Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#66 2024-05-13 15:58:23 +02:00
Implement correlation

Should implement an enum to choose between RMS and correlation.

Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#72 2024-05-13 14:50:27 +02:00
Correctly implement iterative mean for camera attribution

Implementation

So now let us correct the implementation to implement:

  • $prnu_{camera}^l = \text{mean}(image_{camera}^{training_j} − denoiser(image_{camera}^{training_j}, l, camera) \text{…
Correctly implement iterative mean for camera attribution
5e6702e007 #63: use scipy.stats.pearsonr
1a910d2e2c Add a try with numpy.corrcoef
6ed57bc477 Precise purpose of a part of code
Compare 2 commits »
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#63 2024-05-13 12:05:49 +02:00
Attribute source camera

I doubt that current iterative mean leverages the best its current knowledge of the training set.

More precisely we add to the estimated PRNU mean the image substracted the average of already…

3f7b6bf21a Replace _ with for end-user color channel name
b2ffc20473 Add comments to attribute_source_camera.py
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#66 2024-05-13 02:34:42 +02:00
Implement correlation

Unclear if correlation should in absolute be close to 1 or just close to 1.

Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#63 2024-05-13 02:32:56 +02:00
Attribute source camera

Is it the most appropriate, while staying in the use-case context, to compute {min,max}Color only on training images. Do not forget all denoisers.

226936020e Remove comments of {min,max}Color
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#68 2024-05-13 01:43:33 +02:00
Use same color map for different images to more accurately compare them

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…

Benjamin_Loison opened issue Benjamin_Loison/linux#23 2024-05-13 00:45:47 +02:00
How to write ISO to USB key like Disk Image Writer?
Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#68 2024-05-12 03:15:38 +02:00
Use same color map for different images to more accurately compare them

Could locate such maximums to potentially see visually other outliers.

Benjamin_Loison commented on issue Benjamin_Loison/Robust_image_source_identifica...#68 2024-05-12 03:10:39 +02:00
Use same color map for different images to more accurately compare them
>>> index = -100; np.partition(npArray.flatten(), index)[index]
27.28
>>> index = 100; np.partition(npArray.flatten(), index)[index]
-26.35