function [] = bead_landscape(cycle_startnum, experiment_path, display_flag, save_flag)

fret_num = cycle_startnum;
cy5_num = cycle_startnum+1;
cy3_num = cycle_startnum+2;
txred_num = cycle_startnum+3;

fret_name = sprintf('%03d', fret_num);
cy5_name = sprintf('%03d', cy5_num);
cy3_name = sprintf('%03d', cy3_num);
txred_name = sprintf('%03d', txred_num);

fret_bgfilename = [experiment_path '/PROC/BEADMEAN/' fret_name '.beadmean'];
cy5_bgfilename = [experiment_path '/PROC/BEADMEAN/' cy5_name '.beadmean'];
cy3_bgfilename = [experiment_path '/PROC/BEADMEAN/' cy3_name '.beadmean'];
txred_bgfilename = [experiment_path '/PROC/BEADMEAN/' txred_name '.beadmean'];

fid = fopen(fret_bgfilename);
data_col_vector = uint16(fread(fid, inf, 'uint16'));
bkgrnd(:,1) = data_col_vector';
fclose(fid);

fid = fopen(cy5_bgfilename);
data_col_vector = uint16(fread(fid, inf, 'uint16'));
bkgrnd(:,2) = data_col_vector';
fclose(fid);

fid = fopen(cy3_bgfilename);
data_col_vector = uint16(fread(fid, inf, 'uint16'));
bkgrnd(:,3) = data_col_vector';
fclose(fid);

fid = fopen(txred_bgfilename);
data_col_vector = uint16(fread(fid, inf, 'uint16'));
bkgrnd(:,4) = data_col_vector';
fclose(fid);

pos = load([experiment_path '/PROC/DATASETS/position_list.dat']);
pos(:,2) = pos(:,2).*-1;

if(display_flag==1)
    signal_vs_position(pos(:,2:3), bkgrnd(:,1),1,0);
    colorbar;
xlabel('Bead array X (left to right); -1 * stage X coord');
ylabel('Bead array Y (top to bottom); stage Y coord');
    title(['Mean bead fluorescence for cycle ' fret_name ' (FRET)']);
view([0 90]);
    if(save_flag==1)
        savefilename = [experiment_path '/QC/qc_beadland_' fret_name];
        saveas(gcf, savefilename, 'bmp');
    end

    signal_vs_position(pos(:,2:3), bkgrnd(:,2),1,0);
    colorbar;
xlabel('Bead array X (left to right); -1 * stage X coord');
ylabel('Bead array Y (top to bottom); stage Y coord');
    title(['Mean bead fluorescence for cycle ' cy5_name ' (Cy5)']);
view([0 90]);
    if(save_flag==1)
        savefilename = [experiment_path '/QC/qc_beadland_' cy5_name];
        saveas(gcf, savefilename, 'bmp');
    end

    signal_vs_position(pos(:,2:3), bkgrnd(:,3),1,0);
    colorbar;
xlabel('Bead array X (left to right); -1 * stage X coord');
ylabel('Bead array Y (top to bottom); stage Y coord');
    title(['Mean bead fluorescence for cycle ' cy3_name ' (Cy3)']);
view([0 90]);
    if(save_flag==1)
        savefilename = [experiment_path '/QC/qc_beadland_' cy3_name];
        saveas(gcf, savefilename, 'bmp');
    end

    signal_vs_position(pos(:,2:3), bkgrnd(:,4),1,0);
    colorbar;
xlabel('Bead array X (left to right); -1 * stage X coord');
ylabel('Bead array Y (top to bottom); stage Y coord');
title(['Mean bead fluorescence for cycle ' txred_name ' (TxRed)']);
view([0 90]);
    if(save_flag==1)
        savefilename = [experiment_path '/QC/qc_beadland_' txred_name];
        saveas(gcf, savefilename, 'bmp');
    end
end

    

