compMultispFlash/s_simNatural100.m at master · hblasins/compMultispFlash · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
% Simulate the appearance of a Natural100 chart under ambient illuminant
% and a number of narrowband LED lights. Use ambient illuminant with different
% CIE xy chromaticity coordinates.
%
% Copyright, Henryk Blasinski 2017

close all;
clear all;
clc;

ieInit;

%%

wave = 400:5:700;
wave = wave(:);
nWaves = length(wave);

% Define multispectral flash LED spectra
fName = fullfile(cmfRootPath,'Parameters','ximeaLights');
flashSpd = ieReadSpectra(fName,wave);
flashSpd = Energy2Quanta(wave,flashSpd);
flashNorm = flashSpd/max(flashSpd(:));

nChannels = size(flashSpd,2);

% Define the chromaticity coordinates of the illuminant
xVec = linspace(0,1,10);
yVec = linspace(0,1,10);
nX = length(xVec);
nY = length(yVec);

% Define camera models

cameras = {'AptinaMT9M031','AptinaMT9M131',...
           'Canon1DMarkIII','Canon5DMarkII','Canon20D','Canon40D','Canon50D','Canon60D','Canon300D','Canon500D','Canon600D',...
           'HasselbladH2',...
           'NikonD1','NikonD3','NikonD3X','NikonD40','NikonD50','NikonD70','NikonD80','NikonD90','NikonD100','NikonD200',...
           'NikonD200IR','NikonD300s','NikonD700','NikonD5100',...
           'NokiaN900',...
           'OlympusE-PL2',...
           'PentaxK-5','PentaxQ',...
           'PhaseOne',...
           'PointGreyGrasshopper50S5C','PointGreyGrasshopper214S5C',...
           'SONYNEX-5N'};

nCameras = length(cameras);

%% Camera simulation

measurement = cell(nX,nY,nCameras);
for xx=1:nX
    for yy=1:nY

        [ spd ] = xy2Spectrum( xVec(xx), yVec(yy), wave );
        if sum(isnan(spd)) > 0
            continue,
        end

        for cam=1:nCameras

            sensor = sensorCreate('bayer (rggb)');
            sensor = sensorSet(sensor,'wave',wave);
            sensor = sensorSet(sensor,'noise flag',0);

            fName = fullfile(isetRootPath,'data','sensor','colorfilters',cameras{cam});
            cameraResp = ieReadSpectra(fName,wave);
            cameraResp(isnan(cameraResp)) = 0;
            sensor = sensorSet(sensor,'filter transmissivities',cameraResp);

            measurement{xx,yy,cam} = renderData(wave, sensor, spd, flashNorm, 'compact',true);

            % Clear ISET data
            vcDeleteSomeObjects('scene',1:length(vcGetObjects('scene')));
            vcDeleteSomeObjects('oi',1:length(vcGetObjects('oi')));
            vcDeleteSomeObjects('sensor',1:length(vcGetObjects('sensor')));
            vcDeleteSomeObjects('ip',1:length(vcGetObjects('ip')));

        end

    end
end

%%
fName = fullfile(cmfRootPath,'Results','simNatural100V2.mat');
save(fName);