Initial code to run script.

run.m

@@ -0,0 +1,59 @@
+% each scene
+scene_files = dir('scenes/*.jpg');
+scene_files = {scene_files.name};
+
+% color spaces
+color_spaces = {};
+color_spaces{end + 1} = @(img) img;
+color_spaces{end + 1} = @im_rgb2lab;
+color_spaces{end + 1} = @im_rgb2upvpl;
+color_spaces{end + 1} = @im_rgb2uvl;
+color_spaces{end + 1} = @im_rgb2xyy;
+color_spaces{end + 1} = @im_rgb2xyz;
+color_spaces{end + 1} = @im_rgb2lab;
+color_spaces{end + 1} = @(img) select_channel(im_rgb2lab(img), 2:3);
+color_spaces{end + 1} = @(img) select_channel(im_rgb2upvpl(img), 1:2);
+color_spaces{end + 1} = @(img) select_channel(im_rgb2uvl(img), 1:2);
+color_spaces{end + 1} = @(img) select_channel(im_rgb2xyz(img), [1 3]);
+color_spaces{end + 1} = @(img) select_channel(im_rgb2lab(img), [1 3]);
+
+% tuned for my laptop (3 workers)
+parpool('local', 3);
+
+% build
+for i = 1:length(scene_files)
+    % file name
+    scene_file = scene_files{i};
+
+    % build scene and target
+    [scene, target] = build_scene(scene_file);
+
+    % save
+    save(['output/' scene_file '.mat'], 'scene', 'target');
+
+    fprintf('Scene %s...\n', scene_file);
+
+    % for each color space
+    for j = 1:length(color_spaces)
+        % get callback
+        color_cb = color_spaces{j};
+
+        % transform image
+        img = color_cb(scene);
+
+        fprintf('Color %d...\n', j);
+
+        % run four algorithms
+        img_dwest = dwest(img);
+        save(sprintf('output/%s-%d-dwest.mat', scene_file, j), 'img_dwest');
+
+        img_nswtd = nswtd(img);
+        save(sprintf('output/%s-%d-nswtd.mat', scene_file, j), 'img_nswtd');
+
+        img_mwnswtd = mwnswtd(img);
+        save(sprintf('output/%s-%d-mwnswtd.mat', scene_file, j), 'img_mwnswtd');
+
+        img_pcag = pcag(img);
+        save(sprintf('output/%s-%d-pcag.mat', scene_file, j), 'img_pcag');
+    end 
+end