" if ($Phantom_Name eq ""); -$error .= "Phantom NX and NY must be between 5 and 1024
" if ($Phantom_Nx < 5 || $Phantom_Nx > 1024 || $Phantom_Ny < 5 || $Phantom_Ny > 1024); -$error .= "Phantom NSample must be between 1 and 10
" if ($Phantom_NSample < 1 || $Phantom_NSample > 10); - -my $RS_NDet = FilterToNumber($in{'RS_NDet'}); -my $RS_NRot = FilterToNumber($in{'RS_NRot'}); -my $RS_NRay = FilterToNumber($in{'RS_NRay'}); -my $RS_RotAngle = FilterToNumber($in{'RS_RotAngle'}); -$error .= "Raysum NDet must be between 5 and 1800
" if ($RS_NDet < 5 || $RS_NDet > 1800); -$error .= "Raysum NRot must be between 5 and 2048
" if ($RS_NRot < 5 || $RS_NRot > 2048); -$error .= "Raysum RotAngle must be between 0.1 and 2
" if ($RS_RotAngle < 0.1 || $RS_RotAngle > 2); - -#my $IR_Nx = FilterToNumber($in{'IR_Nx'}); -#my $IR_Ny = FilterToNumber($in{'IR_Ny'}); -my $IR_Nx = $Phantom_Nx; -my $IR_Ny = $Phantom_Ny; -my $IR_Filter = FilterMetaChars($in{'IR_Filter'}); -my $IR_Filter_Param = FilterToNumber($in{'IR_Filter_Param'}); -my $IR_Interp = FilterMetaChars($in{'IR_Interp'}); -my $IR_Backproj = FilterMetaChars($in{'IR_Backproj'}); - -my $Disp_Min = "auto"; -my $Disp_Max = "auto"; -$Disp_Min = FilterToNumber($in{'Disp_Min'}) if ($in{'Disp_Min'} ne "auto" && $in{'Disp_Min'} ne ""); -$Disp_Max = FilterToNumber($in{'Disp_Max'}) if ($in{'Disp_Max'} ne "auto" && $in{'Disp_Max'} ne ""); -if ($Disp_Min ne 'auto' && ! ($Disp_Min =~ /^[\d\.]+$/)) { - $error .= "Display Minimum must be 'auto' or numeric (received '$Disp_Min')
"; -} -if ($Disp_Max ne 'auto' && ! ($Disp_Max =~ /^[\d\.]+$/)) { - $error .= "Display Maximum must be 'auto' or numeric (received '$Disp_Max')
"; -} - -my $MPI_Str = FilterMetaChars($in{'MPI'}); -my $MPI = 0; -$MPI = 1 if ($MPI_Str eq "yes" && $::mpi_enable ne ""); - -$error .= "IR Nx and Ny must be between 5 and 1024
" if ($IR_Nx < 5 || $IR_Nx > 1024 || $IR_Ny < 5 || $IR_Ny > 1024); -$error .= "IR Filter Parameter must be between 0 and 1
" if ($IR_Filter_Param < 0 || $IR_Filter_Param > 1); - - -my $tmpid = $$; -my $auto_window_img = "std0.1"; -my $auto_window_diff = "std1"; -my $auto_window_rs = "full"; -my $logfile = "$::jobdir/ctsim.log"; - -my $result_fname = "$::datadir/result-$tmpid.html"; -my $phantom_fname = "$::datadir/phantom-$tmpid.sdf"; -my $rs_fname = "$::datadir/rs-$tmpid.rs"; -my $ir_fname = "$::datadir/ir-$tmpid.sdf"; -my $rs_sdf_fname = "$::datadir/rs-$tmpid.sdf"; -my $diff_fname = "$::datadir/diff-$tmpid.sdf"; -my $phantom_png = "$::datadir/phantom-$tmpid.png"; -my $ir_png = "$::datadir/ir-$tmpid.png"; -my $rs_png = "$::datadir/rs-$tmpid.png"; -my $diff_png = "$::datadir/diff-$tmpid.png"; - -my $result_url = "$::url_datadir/result-$tmpid.html"; -my $phantom_png_url = "$::url_datadir/phantom-$tmpid.png"; -my $ir_png_url = "$::url_datadir/ir-$tmpid.png"; -my $rs_png_url = "$::url_datadir/rs-$tmpid.png"; -my $diff_png_url = "$::url_datadir/diff-$tmpid.png"; - -my $ctrec_ver = "$::bindir/ctrec"; -my $phm2rs_ver = "$::bindir/phm2rs"; -my $phm2sdf_ver = "$::bindir/phm2sdf"; -my $diff_ver = "$::bindir/sdf-2"; -$ctrec_ver = "/opt/lam/bin/mpirun N N $::bindir/ctrec-lam" if $MPI; -$phm2rs_ver = "/opt/lam/bin/mpirun N N $::bindir/phm2rs-lam" if $MPI; -$phm2sdf_ver = "/opt/lam/bin/mpirun N N $::bindir/phm2sdf-lam" if $MPI; - -my $gp_cmd = "$phm2sdf_ver $phantom_fname $Phantom_Nx $Phantom_Ny --phantom $Phantom_Name --nsample $Phantom_NSample"; -my $rs_cmd = "$phm2rs_ver $rs_fname $RS_NDet $RS_NRot --phantom $Phantom_Name --nray $RS_NRay --rotangle $RS_RotAngle"; -my $rs_sdf_cmd = "$::bindir/rs2sdf $rs_fname $rs_sdf_fname"; -my $ir_cmd = "$ctrec_ver $rs_fname $ir_fname $IR_Nx $IR_Ny --filter $IR_Filter --filter-param $IR_Filter_Param --interp $IR_Interp --backproj $IR_Backproj"; -my $diff_cmd = "$diff_ver $phantom_fname $ir_fname $diff_fname --comp"; - -my $window_options = "--auto $auto_window_img"; -if ($Disp_Min ne 'auto') { - $window_options .= " --min $Disp_Min"; -} -if ($Disp_Max ne 'auto') { - $window_options .= " --max $Disp_Max"; -} - -my $png1_cmd = "$::bindir/sdf2img $phantom_fname $phantom_png $window_options --stats --format png"; -my $png2_cmd = "$::bindir/sdf2img $ir_fname $ir_png $window_options --stats --format png"; -my $png3_cmd = "$::bindir/sdf2img $rs_sdf_fname $rs_png --auto $auto_window_rs --stats --format png"; -my $png4_cmd = "$::bindir/sdf2img $diff_fname $diff_png --auto $auto_window_diff --stats --format png"; - -my $title = "CT Simulation Results"; - -my $out = head(); -$out .= "
$title
\n"; - -if ($opt_d) { - $out .= "
Commands
\n"; - $out .= "$gp_cmd\n$rs_cmd
\n$rs_sdf_cmd
\n$ir_cmd
\n$diff_cmd
\n$png1_cmd
\n$png2_cmd
\n" . - "$png3_cmd
\n$png4_cmd
\n"; -} - -my $cmdout = ""; -if ($error ne "") { - $out .= "The following errors were present in your request.
\n"; - $out .= "Please correct them and try submitting your request again.
\n"; - $out .= "$error"; -} else { - my $gp_out; - my $rs_out; - my $rs_sdf_out; - my $ir_out; - my $diff_out; - my $png_gp_out; - my $png_ir_out; - my $png_rs_out; - my $png_diff_out; - if (1) { - $gp_out = `$gp_cmd`; - $rs_out .= `$rs_cmd`; - $rs_sdf_out .= `$rs_sdf_cmd`; - $ir_out .= `$ir_cmd`; - $diff_out .= `$diff_cmd`; - $png_gp_out .= `$png1_cmd`; - $png_ir_out .= `$png2_cmd`; - $png_rs_out .= `$png3_cmd`; - $png_diff_out .= `$png4_cmd`; - } else { - $gp_out = `$gp_cmd 2>&1`; - $rs_out .= `$rs_cmd 2>&1`; - $rs_sdf_out .= `$rs_sdf_cmd 2>&1`; - $ir_out .= `$ir_cmd 2>&1`; - $diff_out .= `$diff_cmd 2>&1`; - $png_gp_out .= `$png1_cmd 2>&1`; - $png_ir_out .= `$png2_cmd 2>&1`; - $png_rs_out .= `$png3_cmd 2>&1`; - $png_diff_out .= `$png4_cmd 2>&1`; - } - $cmdout = "$gp_cmd\n $gp_out $rs_cmd\n $rs_out $rs_sdf_cmd\n $rs_sdf_out $ir_cmd\n $ir_out $diff_cmd\n $diff_out $png1_cmd\n $png_gp_out $png2_cmd\n $png_ir_out $png3_cmd\n $png_rs_out $png4_cmd\n $png_diff_out"; - if (open(LOGFILE,">> $logfile")) { - flock(LOGFILE,LOCK_EX); - seek(LOGFILE, 0, 2); - print LOGFILE "Job $tmpid\n"; - print LOGFILE $cmdout; - print LOGFILE "----------------------------------------------------\n"; - flock(LOGFILE,LOCK_UN); - close(LOGFILE); - } - $gmt_end = time(); - $execution_time = $gmt_end - $gmt_start; - if ($opt_d) - { - $out .= "
Command Output
$cmdout\n"; - } - my $png_gp_out_html = $png_gp_out; - my $png_ir_out_html = $png_ir_out; - my $png_rs_out_html = $png_rs_out; - my $png_diff_out_html = $png_diff_out; - $png_gp_out_html =~ s/\n/
/gms; - $png_ir_out_html =~ s/\n/
/gms; - $png_rs_out_html =~ s/\n/
/gms; - $png_diff_out_html =~ s/\n/
/gms; - $out .= "
| Phantom Image | Reconstructed Image |
$png_gp_out | \n";
- $out .= "$png_ir_out |
| Raysum Sinusoid | Phantom/Reconst Error |
$png_rs_out | \n";
- $out .= "$diff_out $png_diff_out |
\n"; -$out .= "Return to CT Simulation Designer.
\n"; -$out .= "Return to Med-Info.\n"; -$out .= " "; -$out .= "\n"; - -if (open(OUTFILE,"> $result_fname")) -{ - flock(OUTFILE,LOCK_EX); - print OUTFILE $out; - flock(OUTFILE,LOCK_UN); - close OUTFILE; - print "Location: $result_url\n\n"; -} -else -{ - print "Content-type: text/plain\n\n"; - print "The simulator was unable to create an result file.\n"; -} -if (open(JOBFILES,"> $::jobdir/$tmpid")) -{ - flock(JOBFILES,LOCK_EX); - print JOBFILES "gmt_start=$gmt_start\n"; - print JOBFILES "gmt_end=$gmt_end\n"; - print JOBFILES "execution_time=$execution_time\n"; - print JOBFILES "Phantom_Name=$Phantom_Name\n"; - print JOBFILES "Phantom_Nx=$Phantom_Nx\n"; - print JOBFILES "Phantom_Ny=$Phantom_Nx\n"; - print JOBFILES "Phantom_NSample=$Phantom_NSample\n"; - print JOBFILES "RS_NDet=$RS_NDet\n"; - print JOBFILES "RS_NRot=$RS_NRot\n"; - print JOBFILES "RS_NRay=$RS_NRay\n"; - print JOBFILES "RS_RotAngle=$RS_RotAngle\n"; - print JOBFILES "IR_Nx=$IR_Nx\n"; - print JOBFILES "IR_Ny=$IR_Ny\n"; - print JOBFILES "IR_Interp=$IR_Interp\n"; - print JOBFILES "IR_Filter=$IR_Filter\n"; - print JOBFILES "IR_Filter_Param=$IR_Filter_Param\n"; - print JOBFILES "IR_Backproj=$IR_Backproj\n"; - print JOBFILES "Disp_Min=$Disp_Min\n"; - print JOBFILES "Disp_Max=$Disp_Max\n"; - print JOBFILES "MPI=$MPI\n"; - print JOBFILES "Files=$result_fname,$phantom_fname,$rs_fname,$ir_fname,$phantom_png,$ir_png,$rs_sdf_fname,$rs_png\n" if ($error eq ""); - printf JOBFILES "cmdout=$cmdout\n"; - flock(JOBFILES,LOCK_UN); - close JOBFILES; -} - - -exit; - - -sub internal_error -{ - my $out = head(); - $out .= "