% Download and unzip the density grid file for VISp SectionDataSet

    % 100 micron volume size
    sizeGrid = [133 81 115];
    % DENSITY = 3-D matrix of projection density grid volume
    fid = fopen('density.raw', 'r', 'l' );
    DENSITY = fread( fid, prod(sizeGrid), 'float' );
    fclose( fid );
    DENSITY = reshape(DENSITY,sizeGrid);

    % Display one coronal and one sagittal section
    figure;imagesc(squeeze(DENSITY (95,:,:)));caxis([0,1]);colormap(hot);
    figure;imagesc(squeeze(DENSITY (:,:,72)));caxis([0,1]);colormap(hot);