function IntensityAnalysis
workingfolder = cd;
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)
plot(t,x,'r+')
xlim([0 1600])
ylim([10 18])
xlabel('time [min]')
ylabel('mVenus intensity')

cd(workingfolder)
cd ..
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);
subplot(2,1,2)
plot(t,x,'b+')
ylim([7 9])

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

cd(workingfolder)

