function band_modis1b_scatter_callback

handles = guidata(gcf);
s = get(gcbo,'Label');
[t, r] = strtok(s, ' ()');
s = r;
[t, r] = strtok(s, ' ()');
band = str2num(t);
s = r;
[t, r] = strtok(s, ' ()');
wn = str2num(t);

visible=[1:19,26];

handles.modis.select_name = 'MODIS Level1b';
switch handles.modis.axis_variable
  case 'x'
    handles.modis.x_band = band;
    handles.modis.x_label = ['Band ' num2str(handles.modis.x_band)];
    if ismember(band, visible)
      handles.modis.x_label = ['Band ' num2str(handles.modis.x_band) ' (%)'];
    else
      handles.modis.x_label = ['Band ' num2str(handles.modis.x_band) ' (k)'];
    end
    handles.modis.x_data = read_mod1b(handles.modis.filename.level1b, handles.modis.x_band);
  case 'y'
    handles.modis.y_band = band;
    handles.modis.y_label = ['Band ' num2str(handles.modis.y_band)];
    if ismember(band, visible)
      handles.modis.y_label = ['Band ' num2str(handles.modis.y_band) ' (%)'];
    else
      handles.modis.y_label = ['Band ' num2str(handles.modis.y_band) ' (k)'];
    end
    handles.modis.y_data = read_mod1b(handles.modis.filename.level1b, handles.modis.y_band);
end

handles = get_modis_scatter(handles);
scatter_plot(handles)

% Update handles structure
guidata(gcbo, handles);
