9 require './ctsim.conf';
12 $ENV{'PATH'} = "$::bindir:$::lamrundir:/bin";
14 my $fromhost = $ENV{'REMOTE_HOST'};
16 my $gmt_start = time();
24 # Phantom_Name, Phantom_Nx, Phantom_Ny, Phantom_NSample
25 # PJ_NDet, PJ_NRot, PJ_NRay, PJ_RotAngle,
26 # IR_Nx, IR_Ny, IR_Filter, IR_Filter_Param
30 if ($::single_process_only > 0) {
31 my @processes = '/bin/ps -A';
33 foreach $p (@processes) {
34 if ($p ~= /ctsimtext/) {
39 $error .= "Another online simulation is currently running.<p>Click your browser's <i>Refresh</i> button to resubmit your simulation<br>";
44 my $Phantom_Name = FilterMetaChars($in{'Phantom_Name'});
45 my $Phantom_Nx = FilterToNumber($in{'Phantom_Nx'});
46 my $Phantom_Ny = FilterToNumber($in{'Phantom_Ny'});
47 my $Phantom_NSample = FilterToNumber($in{'Phantom_NSample'});
48 $error .= "Phantom name must not be blank<br>" if ($Phantom_Name eq "");
49 $error .= "Phantom NX and NY must be between 5 and 512<br>" if ($Phantom_Nx < 5 || $Phantom_Nx > 512 || $Phantom_Ny < 5 || $Phantom_Ny > 512);
50 $error .= "Phantom NSample must be between 1 and 5<br>" if ($Phantom_NSample < 1 || $Phantom_NSample > 5);
52 my $PJ_NDet = FilterToNumber($in{'PJ_NDet'});
53 my $PJ_NRot = FilterToNumber($in{'PJ_NRot'});
54 my $PJ_NRay = FilterToNumber($in{'PJ_NRay'});
55 my $PJ_RotAngle = FilterToNumber($in{'PJ_RotAngle'});
56 $error .= "Projection NDet must be between 5 and 1000<br>" if ($PJ_NDet < 5 || $PJ_NDet > 1000);
57 $error .= "Projection NRot must be between 5 and 1000<br>" if ($PJ_NRot < 5 || $PJ_NRot > 1000);
58 $error .= "Projection NRay must be between 1 and 5<br>" if ($PJ_NRay < 1 || $PJ_NRay > 5);
59 $error .= "Projection RotAngle must be between 0.1 and 2<br>" if ($PJ_RotAngle < 0.1 || $PJ_RotAngle > 2);
61 #my $IR_Nx = FilterToNumber($in{'IR_Nx'});
62 #my $IR_Ny = FilterToNumber($in{'IR_Ny'});
63 my $IR_Nx = $Phantom_Nx;
64 my $IR_Ny = $Phantom_Ny;
65 my $IR_FilterMethod = FilterMetaChars($in{'IR_FilterMethod'});
66 my $IR_Zeropad = FilterToNumber($in{'IR_Zeropad'});
67 my $IR_Filter = FilterMetaChars($in{'IR_Filter'});
68 my $IR_Filter_Param = FilterToNumber($in{'IR_Filter_Param'});
69 my $IR_Interp = FilterMetaChars($in{'IR_Interp'});
70 my $IR_Backproj = FilterMetaChars($in{'IR_Backproj'});
72 my $Disp_Min = "auto";
73 my $Disp_Max = "auto";
74 $Disp_Min = FilterToNumber($in{'Disp_Min'}) if ($in{'Disp_Min'} ne "auto" && $in{'Disp_Min'} ne "");
75 $Disp_Max = FilterToNumber($in{'Disp_Max'}) if ($in{'Disp_Max'} ne "auto" && $in{'Disp_Max'} ne "");
76 if ($Disp_Min ne 'auto' && ! ($Disp_Min =~ /^[\d\.\-]+$/)) {
77 $error .= "Display Minimum must be 'auto' or numeric (received $Disp_Min) <br>";
79 if ($Disp_Max ne 'auto' && ! ($Disp_Max =~ /^[\d\.\-]+$/)) {
80 $error .= "Display Maximum must be 'auto' or numeric (received $Disp_Max) <br>";
83 my $MPI_Str = FilterMetaChars($in{'MPI'});
85 $MPI = 1 if ($MPI_Str eq "yes" && $::mpi_enable ne "");
87 $error .= "IR Nx and Ny must be between 5 and 1024<br>" if ($IR_Nx < 5 || $IR_Nx > 1024 || $IR_Ny < 5 || $IR_Ny > 1024);
88 $error .= "IR Filter Parameter must be between 0 and 1<br>" if ($IR_Filter_Param < 0 || $IR_Filter_Param > 1);
92 my $auto_window_img = "std0.1";
93 my $auto_window_diff = "std1";
94 my $auto_window_pj = "full";
95 my $logfile = "$::jobdir/ctsim.log";
97 my $result_fname = "$::datadir/result-$tmpid.html";
98 my $phantom_fname = "$::datadir/phantom-$tmpid.if";
99 my $pj_fname = "$::datadir/pj-$tmpid.pj";
100 my $ir_fname = "$::datadir/ir-$tmpid.if";
101 my $pj_if_fname = "$::datadir/pj-$tmpid.if";
102 my $sub_fname = "$::datadir/sub-$tmpid.if";
103 my $phantom_png = "$::datadir/phantom-$tmpid.png";
104 my $ir_png = "$::datadir/ir-$tmpid.png";
105 my $pj_png = "$::datadir/pj-$tmpid.png";
106 my $sub_png = "$::datadir/sub-$tmpid.png";
108 my $result_url = "$::url_datadir/result-$tmpid.html";
109 my $phantom_png_url = "$::url_datadir/phantom-$tmpid.png";
110 my $ir_png_url = "$::url_datadir/ir-$tmpid.png";
111 my $pj_png_url = "$::url_datadir/pj-$tmpid.png";
112 my $sub_png_url = "$::url_datadir/sub-$tmpid.png";
114 my $pjrec_ver = "$::bindir/pjrec";
115 my $phm2pj_ver = "$::bindir/phm2pj";
116 my $phm2if_ver = "$::bindir/phm2if";
117 my $diff_ver = "$::bindir/if2";
119 $pjrec_ver = "/opt/lam/bin/mpirun N N $::lamrundir/pjrec-lam" if $MPI;
120 $phm2pj_ver = "/opt/lam/bin/mpirun N N $::lamrundir/phm2pj-lam" if $MPI;
121 $phm2if_ver = "/opt/lam/bin/mpirun N N $::lamrundir/phm2if-lam" if $MPI;
123 my $gp_cmd = "$phm2if_ver $phantom_fname $Phantom_Nx $Phantom_Ny --phantom $Phantom_Name --nsample $Phantom_NSample";
124 my $pj_cmd = "$phm2pj_ver $pj_fname $PJ_NDet $PJ_NRot --phantom $Phantom_Name --nray $PJ_NRay --rotangle $PJ_RotAngle";
125 my $pj_if_cmd = "$::bindir/pj2if $pj_fname $pj_if_fname";
126 my $pjrec_cmd = "$pjrec_ver $pj_fname $ir_fname $IR_Nx $IR_Ny --filter $IR_Filter --filter-param $IR_Filter_Param --interp $IR_Interp --backproj $IR_Backproj --filter-method $IR_FilterMethod --zeropad $IR_Zeropad";
127 my $sub_cmd = "$diff_ver $phantom_fname $ir_fname $sub_fname --sub";
128 my $diff_cmd = "$diff_ver $phantom_fname $ir_fname --comp";
130 my $window_options = "--auto $auto_window_img";
131 if ($Disp_Min ne 'auto') {
132 $window_options .= " --min $Disp_Min";
134 if ($Disp_Max ne 'auto') {
135 $window_options .= " --max $Disp_Max";
138 my $png1_cmd = "$::bindir/ifexport $phantom_fname $phantom_png $window_options --stats --format png";
139 my $png2_cmd = "$::bindir/ifexport $ir_fname $ir_png $window_options --stats --format png";
140 my $png3_cmd = "$::bindir/ifexport $pj_if_fname $pj_png --auto $auto_window_pj --stats --format png";
141 my $png4_cmd = "$::bindir/ifexport $sub_fname $sub_png --auto $auto_window_diff --stats --format png";
143 my $title = "CT Simulation Results";
146 $out .= "<HTML> <HEAD><TITLE> $title </TITLE></HEAD>\n";
147 $out .= "<BODY BGCOLOR=\"#FFFFFF\" TEXT=\"#000000\" VLINK=\"#800020\" LINK=\"#0000FF\">\n";
148 $out .= "<H1>$title</H1><HR>\n";
151 $out .= "<H2>Commands</H2>\n";
152 $out .= "$gp_cmd<br>\n$pj_cmd<br>\n$pj_if_cmd<br>\n$pjrec_cmd<br>\n$diff_cmd<br>\n$png1_cmd<br>\n$png2_cmd<br>\n" .
153 "$png3_cmd<br>\n$png4_cmd<br>\n";
158 $out .= "<FONT COLOR=\"#FF0000\">The following errors were present in your request.<br>\n";
159 $out .= "Please correct them and try submitting your request again.</FONT><br>\n";
160 $out .= "<i>$error</i>";
173 if (-s $phantom_fname) {
174 $pj_out .= `$pj_cmd`;
175 $png_gp_out .= `$png1_cmd`;
177 $pj_if_out .= `$pj_if_cmd`;
178 $png_pj_out .= `$png3_cmd`;
179 $pjrec_out .= `$pjrec_cmd`;
181 $png_pjrec_out .= `$png2_cmd`;
182 $sub_out .= `$sub_cmd`;
183 $diff_out .= `$diff_cmd`;
184 $png_sub_out .= `$png4_cmd`;
189 $cmdout = "$gp_cmd\n $gp_out $pj_cmd\n $pj_out $pj_if_cmd\n $pj_if_out $pjrec_cmd\n $pjrec_out $diff_cmd\n $diff_out $png1_cmd\n $png_gp_out $png2_cmd\n $png_pjrec_out $png3_cmd\n $png_pj_out $png4_cmd\n $png_sub_out";
190 if (open(LOGFILE,">> $logfile")) {
191 flock(LOGFILE,LOCK_EX);
193 print LOGFILE "Job $tmpid\n";
194 print LOGFILE $cmdout;
195 print LOGFILE "----------------------------------------------------\n";
196 flock(LOGFILE,LOCK_UN);
200 $execution_time = $gmt_end - $gmt_start;
202 $out .= "<H3>Command Output</H3>$cmdout<HR>\n";
204 my $png_gp_out_html = $png_gp_out;
205 my $png_pjrec_out_html = $png_pjrec_out;
206 my $png_pj_out_html = $png_pj_out;
207 my $png_sub_out_html = $png_sub_out;
208 $png_gp_out_html =~ s/\n/<br>/gms;
209 $png_pjrec_out_html =~ s/\n/<br>/gms;
210 $png_pj_out_html =~ s/\n/<br>/gms;
211 $png_sub_out_html =~ s/\n/<br>/gms;
212 $out .= "<TABLE><TR><TD>Phantom Image</TD><TD>Reconstructed Image</TD></TR>\n";
213 $out .= "<TR><TD><IMG SRC=\"$phantom_png_url\"><br><FONT SIZE=1>$png_gp_out</FONT></TD>\n";
214 $out .= "<TD><IMG SRC=\"$ir_png_url\"><br><FONT SIZE=1>$png_pjrec_out</FONT></TD></TR>\n";
215 $out .= "<TR><TD>Projection Sinusoid</TD><TD>Phantom/Reconst Error</TD></TR>\n";
216 $out .= "<TR><TD><IMG SRC=\"$pj_png_url\"><br><FONT SIZE=1>$png_pj_out</FONT></TD>\n";
217 $out .= "<TD><IMG SRC=\"$sub_png_url\"><br><FONT SIZE=2>$sub_out</FONT><br><FONT SIZE=1>$png_sub_out</FONT></TD></TR>\n";
219 $out .= "Execution time: $execution_time seconds\n";
220 $out .= "<p>\nStatistics<br>";
225 $out .= "Return to <A HREF=\"/ctsim/index.html\">CT Simulation Designer</a>.<br>\n";
226 $out .= "Return to <A HREF=\"www.med-info.com\">Med-Info</a>.\n";
227 $out .= "</BODY> </HTML>";
230 if (open(OUTFILE,"> $result_fname"))
232 flock(OUTFILE,LOCK_EX);
234 flock(OUTFILE,LOCK_UN);
236 print "Location: $result_url\n\n";
240 print "Content-type: text/plain\n\n";
241 print "The simulator was unable to create an result file.\n";
243 if (open(JOBFILES,"> $::jobdir/$tmpid"))
245 flock(JOBFILES,LOCK_EX);
246 print JOBFILES "gmt_start=$gmt_start\n";
247 print JOBFILES "gmt_end=$gmt_end\n";
248 print JOBFILES "execution_time=$execution_time\n";
249 print JOBFILES "Phantom_Name=$Phantom_Name\n";
250 print JOBFILES "Phantom_Nx=$Phantom_Nx\n";
251 print JOBFILES "Phantom_Ny=$Phantom_Nx\n";
252 print JOBFILES "Phantom_NSample=$Phantom_NSample\n";
253 print JOBFILES "PJ_NDet=$PJ_NDet\n";
254 print JOBFILES "PJ_NRot=$PJ_NRot\n";
255 print JOBFILES "PJ_NRay=$PJ_NRay\n";
256 print JOBFILES "PJ_RotAngle=$PJ_RotAngle\n";
257 print JOBFILES "IR_Nx=$IR_Nx\n";
258 print JOBFILES "IR_Ny=$IR_Ny\n";
259 print JOBFILES "IR_Interp=$IR_Interp\n";
260 print JOBFILES "IR_Filter=$IR_Filter\n";
261 print JOBFILES "IR_Filter_Param=$IR_Filter_Param\n";
262 print JOBFILES "IR_Backproj=$IR_Backproj\n";
263 print JOBFILES "Disp_Min=$Disp_Min\n";
264 print JOBFILES "Disp_Max=$Disp_Max\n";
265 print JOBFILES "MPI=$MPI\n";
266 print JOBFILES "Files=$result_fname,$phantom_fname,$pj_fname,$ir_fname,$phantom_png,$ir_png,$pj_if_fname,$pj_png\n" if ($error eq "");
267 printf JOBFILES "cmdout=$cmdout\n";
268 flock(JOBFILES,LOCK_UN);
279 $out .= "<H1>Internal error</H1>
280 Please notify <A HREF=mailto:webmaster\@med-info.com>webmaster\@med-info.com</A>
288 # "Content: text/html\n\n";
302 $var =~ /^(\-*\d*\.*\d*)$/;