-
Notifications
You must be signed in to change notification settings - Fork 1
/
code_script_example.m
44 lines (34 loc) · 1.29 KB
/
code_script_example.m
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
load 13C-15N %load msi data (.mat)
%%
% define peak from a peak list
pks=table2struct(readtable('list00.xlsx'));
pkid=10; %pick a peak id (10th peak in the list)
mypk=Mzpk(pks(pkid));
mypk.ppm=5;
mypk.offset=0;
%%
padding=[20,20,20,20];
msi=msi_get_metadata(msi,padding); %initialize, get metadata,ref and alphadata
msi=msi_get_idata(msi,mypk); % get intensity data from mypk
msi=msi_update_imgdata(msi); % update imgdata
%%
figure
imobj=imshow(msi.imgdata,msi.ref);
%%
msi.currentID=1;
msi=msi_get_ms(msi);
msi.ref=imref2d(size(msi.imgdata),msi.res,msi.res); %create axis reference
pad=[1,1,1,1]; %padding
XLim=[min(msi.xx)*msi.res-pad(1),max(msi.xx)*msi.res+pad(3)]; %crop,padding, change XY limits
YLim=[min(msi.yy)*msi.res-pad(2),max(msi.yy)*msi.res+pad(4)];
colorscale=[0,1]; %adjust color brightness
CLim=(max(msi.idata)+1e-9)*colorscale;
% axes1 ----------initial drawing of image
figure
imobj=msi_ini_draw(gca,msi.imgdata,msi.ref,msi.alphadata,parula,CLim,'k',XLim,YLim); %%%%%%%initialize
%imobj.ButtonDownFcn = @axes1_ButtonDownFcn;% ----for ButtonDown action on the image
%%
mypk.M=[0,1,0]
msi=msi_get_idata(msi,mypk);
msi=msi_get_imgdata(msi); %get imgdata from idata
imobj.CData=msi.imgdata; %update image object CData;