Jeff Mather

Help!
I'm trying to use spectrumLabel for a calculated absorption spectrum.
It works fine from the command window, but I'm not able to put the color axis in a plot which is rendered to a GUI.
Here is the piece of code that I put in the .m file. The plot of the extinction, scattering and absorption spectra are sent to the axes style GUI handles.axes1:
----------------------------
axes(handles.axes1);
plot(wave,qext,wave,qsca,wave,qabs);
title('Spectrum');
xlabel('Wavelength (nm)');
ylabel('Arbitrary units');
spectrumLabel(handles.axes1);
-------------------------------
Apparently, the spectrumLabel is put at the right position but with zero or very small vertical dimension.
Thank you!

It works well with the modification of Jeff in 29 Apr 2011. Thank you!

Thanks to this excellent work, I've been able to image a spectrum plot with both the characteristic peaks marked in colors and intensity plot shown. The code is as following.

%%%%%%%%%%%%%%%%%%%%%%
function SpecRGB_Intensity(wavelength,intensity, peak_threshold)
% wavelength and intensity are vectors
% the characteristic peaks are chosen if the value of the peak
% is bigger than the peak_threshold.
% by springface

sRGB = spectrumRGB(wavelength);

%% Transform spectrum plot into an image, marked by colors
len_wave = length(wavelength);
len_inten = 1000;
intensity_plot = max(intensity):-max(intensity)/(len_inten-1):0;
wavelength_plot = wavelength;

intensity_num = ones(size(intensity));

for cc = 1: len_wave
    if intensity(cc)>=0
        [C,I] = min(abs(intensity(cc)- intensity_plot));
        intensity_num(cc) = len_inten-I;
    end
end

im_plot = ones(len_inten,len_wave,3);

for cc = 1: len_wave
    im_plot(1:intensity_num(cc),cc,:) = repmat(sRGB(1,cc,:),intensity_num(cc),1);
end

%% mark the non-peak area as black
spec_idx = intensity < peak_threshold;
sRGB(1,spec_idx,:) = zeros(1,sum(spec_idx),3);

%% Plot figures
figure;

whitegap_len = 2;
RGBspec_len = 100;

RGBspec_range = max(intensity)+whitegap_len+RGBspec_len+1:max(intensity)+1;

imagesc(wavelength_plot,[RGBspec_range,intensity_plot],...
    cat(1,repmat(sRGB,RGBspec_len,1),ones(whitegap_len,len_wave,3),im_plot));
set(gca,'YTick',[])
%%%%%%%%%%%%%%%%%%%%%%