clear all % Data extraction: directory = 'C:\Users\Usuario\Documents\Carlos AV\UPV\iGEM\Modelling\Experimentos\Spectrum\sfGFP Spectrum\'; sheet = 'Plate 1 - Sheet1'; fname = 'sfGFP_Spectrum_data'; FOD_Em = xlsread([directory fname '.xlsx'],sheet,'L353:N748'); wavelengths_em = xlsread([directory fname '.xlsx'],sheet,'B353:B748'); FOD_Ex = xlsread([directory fname '.xlsx'],sheet,'F353:H748'); wavelengths_ex = xlsread([directory fname '.xlsx'],sheet,'B513:B673'); save([directory fname '.mat'],'wavelengths_em','FOD_Em','wavelengths_ex','FOD_Ex'); %% GRAPHICS clear all directory = 'C:\Users\Usuario\Documents\Carlos AV\UPV\iGEM\Modelling\Experimentos\Spectrum\sfGFP Spectrum\'; fname = 'sfGFP_Spectrum_data'; load([directory fname '.mat']) % FOD Vectors FOD_Em_mean = mean(FOD_Em,2); FOD_Ex_mean = mean(FOD_Ex,2); % Maximum peaks indexmax_em = find(max(FOD_Em_mean) == FOD_Em_mean); indexmax_ex = find(max(FOD_Ex_mean) == FOD_Ex_mean); wlmax_em = wavelengths_em(indexmax_em); wlmax_ex = wavelengths_ex(indexmax_ex); FODmax_em = FOD_Em_mean(indexmax_em); FODmax_ex = FOD_Ex_mean(indexmax_ex); % Normalised data FOD_Em_norm = (FOD_Em_mean/FODmax_em)*100; FOD_Ex_norm = (FOD_Ex_mean/FODmax_ex)*100; FODmax_em_norm = FOD_Em_norm(indexmax_em); FODmax_ex_norm = FOD_Ex_norm(indexmax_ex); % RGB Colours yellow = [1 0.843137264251709 0]; green = [0.400000005960464 1 0]; color_489 = [0 251 255]/255; color_498 = [0 255 169]/255; color_504 = [0 255 97]/255; color_513 = [21 255 0]/255; color_514 = [26 255 0]/255; color_515 = [31 255 0]/255; color_516 = [36 255 0]/255; color_529 = [90 255 0]/255; color_586 = [255 236 0]/255; color_611 = [255 152 0]/255; title_spec = 'sfGFP Spectrum'; % Spectrum % Create figure fig_spectrum = figure('Color','w'); % Create axis axes = axes('Parent',fig_spectrum); hold(axes,'on'); % Areas/plots area_em = area(wavelengths_em,FOD_Em_norm,'FaceColor',color_611,'EdgeColor',color_611,'FaceAlpha',0.6,'LineWidth',1); hold on area_ex = area(wavelengths_ex,FOD_Ex_norm,'FaceColor',color_586,'EdgeColor',color_586,'FaceAlpha',0.5,'LineWidth',1); line([wlmax_em wlmax_em],get(axes,'YLim'),'Color','k','LineStyle','--') hold on line([wlmax_ex wlmax_ex],get(axes,'YLim'),'Color','k','LineStyle','--') hold on % Plot peaks plot(wlmax_em,FODmax_em_norm,'k.','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k') hold on plot(wlmax_ex,FODmax_ex_norm,'k.','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k') % Text descr = {['Excitation maximum: ' int2str(wlmax_ex) ' nm']; ['Emission maximum: ' int2str(wlmax_em) ' nm']}; text(0.9,0.9,descr,'Units','normalized','HorizontalAlignment','right','FontName','Lato','FontSize',11) % Other properties title(title_spec,'FontName','Lato','FontSize',14) xlabel('Wavelength (nm)','FontName','Lato') ylabel('Relative intensity (%)','FontName','Lato') legend({'Emission Spectrum','Excitation Spectrum'},'FontName','Lato','Location','east'); wlMIN = min([min(wavelengths_em) min(wavelengths_ex)]); wlMAX = max([max(wavelengths_em) max(wavelengths_ex)]); xlim([wlMIN wlMAX]) box(axes,'on'); set(axes,'FontName','Lato','XGrid','on','YGrid','on');