PAM 2500 荧光数据导出数据整合脚本

Posted 2021-01-31 liulele

PAM 2500导出的数据是以分号为分隔符的纯文本，且一次测量单独保存在文件里，默认以日期加括号里的序号为文件名，我们只需要其中的Fo、Fm和最后一次的数据。

原始数据如下：

 1 "t";"Date";"Time";"No.";"ML";"Temp.";"PAR";"F";"Fm‘";"~Fo‘";"Y(II)";"Y(NPQ)";"Y(NO)";"NPQ";"qN";"qP";"qL";"ETR";"";
 2        7;28.12.18;17:46:40;Type: IC 
 3        9;28.12.18;17:46:41;default_90.par*
 4        7;28.12.18;17:46:40;1;2; 0.0;0;Fo: 0.548;Fm: 3.018
 5        7;28.12.18;17:46:40;1;2; 0.0;0;0.548;3.018;0.548;0.818;0.000;0.182;0.000;0.000;1.000;1.000;  0.0;;
 6       27;28.12.18;17:47:00;2;2; 0.0;319;2.926;2.979;0.547;0.018;0.013;0.970;0.013;0.071;0.023;0.005;  2.4;;
 7       47;28.12.18;17:47:20;3;2; 0.0;319;2.733;2.891;0.544;0.055;0.040;0.906;0.044;0.182;0.078;0.025;  7.3;;
 8       67;28.12.18;17:47:40;4;2; 0.0;319;2.724;2.797;0.540;0.026;0.071;0.903;0.079;0.279;0.041;0.015;  3.5;;
 9       87;28.12.18;17:48:00;5;2; 0.0;319;2.563;2.661;0.535;0.037;0.114;0.849;0.134;0.274;0.055;0.019;  4.9;;
10      107;28.12.18;17:48:20;6;2; 0.0;319;2.370;2.529;0.529;0.063;0.152;0.785;0.193;0.273;0.089;0.027;  8.4;;
11      127;28.12.18;17:48:40;7;2; 0.0;319;2.168;2.396;0.523;0.095;0.186;0.718;0.260;0.299;0.132;0.040; 12.7;;
12      147;28.12.18;17:49:00;8;2; 0.0;319;1.992;2.293;0.518;0.131;0.209;0.660;0.316;0.325;0.181;0.057; 17.6;;
13      167;28.12.18;17:49:20;9;2; 0.0;319;1.848;2.213;0.514;0.165;0.223;0.612;0.364;0.357;0.230;0.078; 22.1;;
14      187;28.12.18;17:49:40;10;2; 0.0;319;1.729;2.157;0.511;0.198;0.229;0.573;0.399;0.377;0.278;0.099; 26.6;;
15      207;28.12.18;17:50:00;11;2; 0.0;319;1.642;2.124;0.509;0.227;0.229;0.544;0.421;0.394;0.322;0.123; 30.4;;
16      221;28.12.18;17:50:14;File: I_181228_174633.PWS

我们需要以提取第3行和第15行的数据整合成一个新的文件，我们把所有需要整合的文件放在同一文件夹下，运行脚本如下：

use strict;
use warnings;


my @dir;
my @myfiles;
open OUT, ">result.csv";
opendir (DIR, "./") or die "cannot open this directory";
@dir = readdir DIR;
foreach my $file (@dir){
    if ($file =~ /.*.CSV/){
        push @myfiles, $file;
    }
}
close(DIR);

foreach my $file (@myfiles){
    open IN, "<".$file or die "cannot open $file";
    my @lines = <IN>;
    chomp @lines;
    my @F_start = split(/;/, $lines[4]);
    my @F_end = split(/;/, $lines[15]);
    print OUT $file.",".substr($F_start[7],4,).",".substr($F_start[8],4,).",".join(",",@F_end)."
";
    close(IN);
}

close(OUT);

手动加入表头行即可，新结果如下：

1 FileName,No,Fo,Fv,Fv/Fm,t,Date,Time,No.,ML,Temp.,PAR,F,Fm‘,~Fo‘,Y(II),Y(NPQ),Y(NO),NPQ,qN,qP,qL,ETR
2 20181227_2(7).CSV,NX08,0.493,2.65,0.813962264,207,27.12.18,20:07:20,11,2,0,444,1.466,2.036,0.467,0.28,0.167,0.553,0.302,0.325,0.391,0.155,52.2
3 20181227_2(8).CSV,NX08,0.512,2.751,0.81388586,207,27.12.18,20:11:29,11,2,0,444,1.652,2.169,0.488,0.238,0.161,0.601,0.268,0.297,0.329,0.119,44.4
4 20181227_2(9).CSV,NX08,0.502,2.596,0.806625578,207,27.12.18,20:15:46,11,2,0,444,1.502,1.982,0.474,0.242,0.179,0.579,0.31,0.327,0.341,0.13,45.2
5 20181228(13).CSV,NX09,0.499,2.817,0.8228612,207,28.12.18,17:26:43,11,2,0,319,1.296,2.091,0.47,0.38,0.16,0.46,0.347,0.34,0.519,0.224,50.9
6 20181228(14).CSV,NX09,0.512,2.713,0.811279027,207,28.12.18,17:30:48,11,2,0,319,1.523,2.111,0.486,0.279,0.16,0.561,0.285,0.314,0.389,0.154,37.3
7 20181228(15).CSV,NX09,0.505,2.823,0.821112292,207,28.12.18,17:35:16,11,2,0,319,1.371,2.21,0.481,0.38,0.135,0.486,0.277,0.299,0.516,0.22,50.9

随后可以导入到R中进行统计学分析，比如求各组的平均值，如下

mean <- tapply(data$Fv.Fm,data$No,mean)