specreader/meshscanplot.m at master · ennogra/specreader · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
function meshscanplot
hFigSpecr = findall(0,'Tag','specr_Fig');
settings = getappdata(hFigSpecr,'settings');
file = settings.file;
scan = settings.scan;
if isempty(file) | isempty(scan) | ~isfield(scan,'selectionIndex') | isempty(scan.selectionIndex)
    return;
end

% --- get z data coloum
hPopupmenuY = findall(hFigSpecr,'Tag','specr_PopupmenuY');
popupmenuYString = get(hPopupmenuY,'String');
popupmenuYValue = get(hPopupmenuY,'Value');

% --- load data and plot
[~,filename,~] = fileparts(file);
for iscan = 1:length(scan.selectionIndex)
    scanParts = regexp(scan.head{scan.selectionIndex(iscan)}, '\ *', 'split');
    if ~strcmpi(scanParts{3},'mesh')
        continue;
    end
    xdata = scan.selection{iscan}.colData(:,1);
    ydata = scan.selection{iscan}.colData(:,2);
    zdata = scan.selection{iscan}.colData(:,popupmenuYValue);
    ptsX = str2double(scanParts{7})+1;
    ptsY = floor(length(xdata)/ptsX);
    nOfPts = ptsX*ptsY;
    if nOfPts == 0
        continue;
    end
    meshdata.X=reshape(xdata(1:nOfPts), ptsX, ptsY)';
    meshdata.Y=reshape(ydata(1:nOfPts), ptsX, ptsY)';
    meshdata.Z=reshape(zdata(1:nOfPts), ptsX, ptsY)';
    if settings.meshScanPlotScale == 1
        meshdata.Z = log10(meshdata.Z);
    end
    title_str = [filename,' - Scan #',num2str(scan.selectionNumber(iscan))];
    figure('Name',title_str);
    switch settings.meshScanPlotType
        case 1
            contour(meshdata.X,meshdata.Y,meshdata.Z);
            colormap;
        case 2
            contourf(meshdata.X,meshdata.Y,meshdata.Z);
            colormap(flipud(colormap('gray')));
        case 3
            mesh(meshdata.X,meshdata.Y,meshdata.Z);
            colormap;
        case 4
            surf(meshdata.X,meshdata.Y,meshdata.Z);
            colormap;
        case 5
            imagesc(meshdata.X(1,:),meshdata.Y(:,1),meshdata.Z);
            colormap;
            set(gca,'ydir','normal');
    end
    xlabel(titlestr(scan.selection{iscan}.colHead{1}));
    ylabel(titlestr(scan.selection{iscan}.colHead{2}));
    %set(gca,'String',popupmenuYString(popupmenuYValue));
    title(titlestr(title_str));
    hC = colorbar;
    ylabel(hC,titlestr(popupmenuYString{popupmenuYValue}));
    if settings.meshScanExport2WS == 1
        varname = ['specrMeshScan',num2str(scan.selectionNumber(iscan))];
        assignin('base',varname,meshdata);
    end
end
return;