function IntensityAnalysis
workingfolder = cd;
cd Time' lapse'/mVenus/

filelist = dir('*.TIF');

for i=1:length(filelist)
        filename = [filelist(i).name];
        image = imread(filename);
        x(i) = sum(sum(image))/(1040*1388);
end

t = 0:10:10*(length(filelist)-1);
subplot(2,1,1)
x(145:150) = NaN;
x(103) = NaN;
x(1) = NaN;
plot(t,x,'r+')
xlim([0 1600])
xlabel('time [min]')
ylabel('mVenus concentration')

cd(workingfolder)
cd Time' lapse'/mCherry/

filelist = dir('*.TIF');

for i=1:length(filelist)
        filename = [filelist(i).name];
        image = imread(filename);
        x(i) = sum(sum(image))/(1040*1388);
end

t = 0:10:10*(length(filelist)-1);
x(145:150) = NaN;
x(103) = NaN;
x(1) = NaN;
subplot(2,1,2)
plot(t,x,'b+')

xlabel('time [min]')
ylabel('mCherry concentration')
xlim([0 1600])

cd(workingfolder)
