Matlab Code

From Psych 221 Image Systems Engineering
Jump to navigation Jump to search

ieInit;

%% Capture the Image

% Import the image I = im2double(imread('Fig1.JPG'));

% Choose which section of the image you want to look at xStart = 1300; xEnd = 2100; yStart = 1202; yEnd = 2225;

% Create the cropped image Ir = I(:,:,1)'; Ib = I(:,:,2)'; Ig = I(:,:,3)';

Icrp(:,:,1) = Ir(yStart:yEnd,xStart:xEnd); Icrp(:,:,2) = Ib(yStart:yEnd,xStart:xEnd); Icrp(:,:,3) = Ig(yStart:yEnd,xStart:xEnd);

Icrp = fliplr(Icrp);


%% Create image display display1 = displayCreate('LCD-Apple');

%% Create an ISET scene scene = sceneFromFile(Icrp, 'rgb', [], display1);

%% Capture the scene with the camera % Create the cameras camera1 = cameraCreate; camera2 = cameraCreate; camera3 = cameraCreate; camera4 = cameraCreate; camera5 = cameraCreate;

camera6 = cameraCreate; camera7 = cameraCreate;

% Adjust camera settings camera1.sensor.pixel.readNoise = 1e-3; camera2.sensor.pixel.readNoise = 1e-2; camera3.sensor.pixel.readNoise = 1e-1; camera4.sensor.pixel.readNoise = 5e-1; camera5.sensor.pixel.readNoise = 1e0;

camera6.sensor.pixel.readNoise = 1.5e-2; camera7.sensor.pixel.readNoise = 1.5e-2;

camera6.sensor.rows = 720; camera6.sensor.cols = 1280;

camera7.sensor.rows = 720; camera7.sensor.cols = 1280; camera7.sensor.pixel.height = 1.87e-6; camera7.sensor.pixel.width = 1.87e-6;

% Capture the images with the cameras cameraCapture1 = cameraCompute(camera1,scene); cameraCapture2 = cameraCompute(camera2,scene); cameraCapture3 = cameraCompute(camera3,scene); cameraCapture4 = cameraCompute(camera4,scene); cameraCapture5 = cameraCompute(camera5,scene);

cameraCapture6 = cameraCompute(camera6,scene); cameraCapture7 = cameraCompute(camera7,scene);


%% Take a look at the image on your display sRGB1 = cameraGet(cameraCapture1,'ip data srgb'); sRGB2 = cameraGet(cameraCapture2,'ip data srgb'); sRGB3 = cameraGet(cameraCapture3,'ip data srgb'); sRGB4 = cameraGet(cameraCapture4,'ip data srgb'); sRGB5 = cameraGet(cameraCapture5,'ip data srgb'); sRGB6 = cameraGet(cameraCapture6,'ip data srgb'); sRGB7 = cameraGet(cameraCapture7,'ip data srgb');

vcNewGraphWin; imagesc(sRGB1);

vcNewGraphWin; imagesc(sRGB2);

vcNewGraphWin; imagesc(sRGB3);

vcNewGraphWin; imagesc(sRGB4);

vcNewGraphWin; imagesc(sRGB5);

vcNewGraphWin; imagesc(sRGB6);

vcNewGraphWin; imagesc(sRGB7);