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
69
|
function extraction
% Pour extraire l'intensité de 7 bandes Raman, il faut indiquer
% les nombres d'onde précisément qui ont été lus à partir des spectres sous
% Labspec.
% Il faut sauvegarder Le fichier ligne ou cartographie en format TXT et
% aussi un spectre ponctuel en format TXT
nu1=0;
nu2=0;
lgn=0;
fpictname=[];fpixname=[];filename=[];filesave=[];
fpictname = input('Nom du fichier image ','s');
fpixname = input('Nom du fichier spectre du pixel ','s');
filesave = input ('Nom du fichier à sauvegarder ','s');
nu1=input('Nombre d onde nu1 ','s');
nu2=input('Nombre d onde nu2 ','s');
nu3=input('Nombre d onde nu3 ','s');
nu4=input('Nombre d onde nu4 ','s');
nu5=input('Nombre d onde nu5 ','s');
nu6=input('Nombre d onde nu6 ','s');
nu7=input('Nombre d onde nu7 ','s');
lgn=input('Nombre de lignes ','s');
image = dlmread([fpictname,'.txt'],' ');
spectre = dlmread([fpixname,'.txt'],' ');
nu1=str2num(nu1);
nu2=str2num (nu2);
nu3=str2num (nu3);
nu4=str2num (nu4);
nu5=str2num (nu5);
nu6=str2num (nu6);
nu7=str2num (nu7);
lgn=str2num(lgn);
lgn=lgn-1;
idx=find (spectre ==nu1);
idx=idx-1;
p1=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
idx=0;
idx=find(spectre== nu2);
idx=idx-1;
p2=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
idx=0;
idx=find(spectre== nu3);
idx=idx-1;
p3=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
idx=0;
idx=find(spectre== nu4);
idx=idx-1;
p4=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
idx=0;
idx=find(spectre== nu5);
idx=idx-1;
p5=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
idx=0;
idx=find(spectre== nu6);
idx=idx-1;
p6=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
idx=0;
idx=find(spectre== nu7);
idx=idx-1;
p7=dlmread([fpictname,'.txt'],' ',[0 idx lgn idx]);
NU=[nu1,nu2,nu3,nu4,nu5,nu6,nu7];
P=[p1,p2,p3,p4,p5,p6,p7];
P=[NU;P];
save([filesave,'.txt'],'P','-ASCII','-TABS'); |
