projects / ctsim.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (parent: b96676b)
r655: *** empty log message ***
authorKevin M. Rosenberg <kevin@rosenberg.net>
Sat, 31 Mar 2001 04:50:29 +0000 (04:50 +0000)
committerKevin M. Rosenberg <kevin@rosenberg.net>
Sat, 31 Mar 2001 04:50:29 +0000 (04:50 +0000)
cgi-bin/ctsim.cgi.in patch | blob | history
html/simulate.html.in patch | blob | history

diff --git a/cgi-bin/ctsim.cgi.in b/cgi-bin/ctsim.cgi.in
index 4ead487c68ac565c3806fbb42266ccddb62497f5..e978881e0ae0d3b3b16e3f49f888b17ffe4f57da 100755 (executable)
--- a/cgi-bin/ctsim.cgi.in
+++ b/cgi-bin/ctsim.cgi.in
@@ -41,10 +41,12 @@ if ($::single_process_only > 0) {
 my $Phantom_Name = FilterMetaChars($in{'Phantom_Name'});
 my $Phantom_Nx = FilterToNumber($in{'Phantom_Nx'});
 my $Phantom_Ny = FilterToNumber($in{'Phantom_Ny'});
+my $Phantom_ViewRatio = FilterToNumber($in{'Phantom_ViewRatio'});
 my $Phantom_NSample = FilterToNumber($in{'Phantom_NSample'});
 $error .= "Phantom name must not be blank<br>" if ($Phantom_Name eq "");
-$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);
+$error .= "Phantom Size must be between 5 and 512<br>" if ($Phantom_Nx < 5 || $Phantom_Nx > 512 || $Phantom_Ny < 5 || $Phantom_Ny > 512);
 $error .= "Phantom NSample must be between 1 and 5<br>" if ($Phantom_NSample < 1 || $Phantom_NSample > 5);
+$error .= "View Ratio must be between 1 and 100<br>" if ($Phantom_ViewRatio < 1 || $Phantom_ViewRatio > 100);
 
 my $PJ_Geometry = FilterMetaChars($in{'PJ_Geometry'});
 my $PJ_NDet = FilterToNumber($in{'PJ_NDet'});
@@ -52,11 +54,13 @@ my $PJ_NRot = FilterToNumber($in{'PJ_NRot'});
 my $PJ_FocalLength = FilterToNumber($in{'PJ_FocalLength'});
 my $PJ_NRay = FilterToNumber($in{'PJ_NRay'});
 my $PJ_RotAngle = FilterToNumber($in{'PJ_RotAngle'});
+my $PJ_ScanRatio = FilterToNumber($in{'PJ_ScanRatio'});
 $error .= "Projection NDet must be between 5 and 1000<br>" if ($PJ_NDet < 5 || $PJ_NDet > 1000);
 $error .= "Projection NRot must be between 5 and 1000<br>" if ($PJ_NRot < 5 || $PJ_NRot > 1000);
 $error .= "Projection NRay must be between 1 and 5<br>" if ($PJ_NRay < 1 || $PJ_NRay > 5);
 $error .= "Projection RotAngle must be between 0.1 and 2<br>" if ($PJ_RotAngle < 0.1 || $PJ_RotAngle > 2);
 $error .= "Focal length must be between 1.01 and 10<b>" if ($PJ_FocalLength <= 1 || $PJ_FocalLength > 10);
+$error .= "Scan Ratio must be between 0 and 10<br>" if ($PJ_ScanRatio <= 0 || $PJ_ScanRatio > 10);
 
 #my $IR_Nx = FilterToNumber($in{'IR_Nx'});
 #my $IR_Ny = FilterToNumber($in{'IR_Ny'});
@@ -120,8 +124,8 @@ $pjrec_ver = "/opt/lam/bin/mpirun N N $::lamrundir/pjrec-lam" if $MPI;
 $phm2pj_ver = "/opt/lam/bin/mpirun N N $::lamrundir/phm2pj-lam" if $MPI;
 $phm2if_ver = "/opt/lam/bin/mpirun N N $::lamrundir/phm2if-lam" if $MPI;
 
-my $gp_cmd = "$phm2if_ver $phantom_fname $Phantom_Nx $Phantom_Ny --phantom $Phantom_Name --nsample $Phantom_NSample";
-my $pj_cmd = "$phm2pj_ver $pj_fname $PJ_NDet $PJ_NRot --phantom $Phantom_Name --nray $PJ_NRay --rotangle $PJ_RotAngle --geometry $PJ_Geometry --focal-length $PJ_FocalLength";
+my $gp_cmd = "$phm2if_ver $phantom_fname $Phantom_Nx $Phantom_Ny --phantom $Phantom_Name --nsample $Phantom_NSample --view-ratio $Phantom_ViewRatio";
+my $pj_cmd = "$phm2pj_ver $pj_fname $PJ_NDet $PJ_NRot --phantom $Phantom_Name --nray $PJ_NRay --rotangle $PJ_RotAngle --geometry $PJ_Geometry --focal-length $PJ_FocalLength --view-ratio $Phantom_ViewRatio --scan-ratio $PJ_ScanRatio";
 my $pj_if_cmd = "$::bindir/pj2if $pj_fname $pj_if_fname";
 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";
 my $sub_cmd = "$diff_ver $phantom_fname $ir_fname $sub_fname --sub";
@@ -140,10 +144,11 @@ my $png2_cmd = "$::bindir/ifexport $ir_fname $ir_png $window_options --stats --f
 my $png3_cmd = "$::bindir/ifexport $pj_if_fname $pj_png --auto $auto_window_pj --stats --format png";
 my $png4_cmd = "$::bindir/ifexport $sub_fname $sub_png --auto $auto_window_diff --stats --format png";
 
-my $title = "CT Simulation Results";
+my $title = "CTSim Results";
 
 my $out = head();
-$out .= "<HTML> <HEAD><TITLE> $title </TITLE></HEAD>\n";
+$out .= "<HTML> <HEAD><TITLE> $title </TITLE><LINK rel=\"stylesheet\" href=\"http://www.ctsim.org/main.css\" type=\"text/css\">
+</HEAD>\n";
 $out .= "<BODY  BGCOLOR=\"#FFFFFF\" TEXT=\"#000000\"  VLINK=\"#800020\" LINK=\"#0000FF\">\n";
 $out .= "<H1>$title</H1><HR>\n";
 
@@ -218,17 +223,18 @@ if ($error ne "") {
   $out .= "<TABLE><TR><TD>Phantom Image</TD><TD>Reconstructed Image</TD></TR>\n";
   $out .= "<TR><TD><IMG SRC=\"$phantom_png_url\"><br><FONT SIZE=1>$png_gp_out</FONT></TD>\n";
   $out .= "<TD><IMG SRC=\"$ir_png_url\"><br><FONT SIZE=1>$png_pjrec_out</FONT></TD></TR>\n";
-  $out .= "<TR><TD>Projection Sinusoid</TD><TD>Phantom/Reconst Error</TD></TR>\n";
+  $out .= "<TR><TD>Projection Sinusoid</TD><TD>Reconstruction Error</TD></TR>\n";
   $out .= "<TR><TD><IMG SRC=\"$pj_png_url\"><br><FONT SIZE=1>$png_pj_out</FONT></TD>\n";
   $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";
   $out .= "</TABLE>";
   $out .= "Execution time: $execution_time seconds\n";
-  $out .= "<p>\nStatistics<br>";
+  $out .= "<p>\nError Measurements<br>";
   $out .= "$diff_out";
 }
 
 $out .= "<HR>\n";
-$out .= "Return to <A HREF=\"http://www.ctsim.org\">CTSim</a>.<br>\n";
+$out .= "Specify another <a href=\"http://www.ctsim.org/online-sim/simulate.html\">simulation</a>.<br>";
+$out .= "Return to CTSim's <A HREF=\"http://www.ctsim.org\">home</a>.<br>\n";
 $out .= "</BODY> </HTML>";
 $out .= "\n";
 
diff --git a/html/simulate.html.in b/html/simulate.html.in
index 88c7abd27d3a795b20a29ca21adadf553837bb9c..a3790788755cd4103a4b841211044d50a1594803 100644 (file)
--- a/html/simulate.html.in
+++ b/html/simulate.html.in
@@ -1,9 +1,17 @@
+<!DOCTYPE HTML PUBLIC "-//IETF//DTD HTML//EN">
 <HTML>
 <HEAD>
-<TITLE>Perform CT Simulation</TITLE>
+<meta name="copyright" content="Copyright (c) 1983-2001, Kevin M. Rosenberg, M.D.">
+<meta name="description" content="CTsim - The Open-Source Computed Tomography Simulator">
+<meta name="keywords" content="ctsim, computed tomography, ct, open source, scientific software, simulation, c++, online simulation, medical imaging">
+<LINK rel="stylesheet" href="http://www.ctsim.org/main.css" type="text/css">
+<TITLE>CTSim Online Simulation</TITLE>
 </HEAD>
 <BODY BGCOLOR="#FFFFFF">
-<H1>Perform CT Simulation</H1>
+<H1><a href="http://www.ctsim.org"">CTSim</a> Online Simulation</H1>
+Reasonable defaults have been filled into this form. For detailed explanations of the
+parameters on this form, please visit CTSim's online <a href="http://www.ctsim.org/manual/ctsim_contents.html">manual</a>.<br>
+
 <FORM action="@cgibinurl@/ctsim.cgi" method="POST">
 <table border=1 cellpadding="5" bgcolor="#FFFFFF">
 <tr><td valign="top">
@@ -11,11 +19,13 @@
 Phantom Name:<br>
 <INPUT type="radio" name="Phantom_Name" value="herman" checked>G.T. Herman<br>
 <INPUT type="radio" name="Phantom_Name" value="rowland">Rowland<br>
-<INPUT type="radio" name="Phantom_Name" value="browland">Bordered Rowland<br>
 <INPUT type="radio" name="Phantom_Name" value="unitpulse">Unit Pulse<br>
 <p>
-Nx: <input type="text" name="Phantom_Nx" size="4" value="256"><br>
-Ny: <input type="text" name="Phantom_Ny" size="4" value="256"><br>
+Image Size:<br>
+Width: <input type="text" name="Phantom_Nx" size="4" value="128"><br>
+Height: <input type="text" name="Phantom_Ny" size="4" value="128"><br>
+<p>
+View Ratio: <input type="text" name="Phantom_ViewRatio" size="4" value="1.1">
 <p>
 Number of Samples<br>per pixel: <input type="text" name="Phantom_NSample" size="2" value="2"><br>
 <p>
@@ -23,57 +33,57 @@ Display Window:<br>
 Minimum: <input type="text" name="Disp_Min" size=5 value="auto"><br>
 Maximum: <input type="text" name="Disp_Max" size=5 value="auto"><br>
 <p> 
-MPI Supercomputing:<br>
+<!--MPI Supercomputing:<br>
 <INPUT type="radio" name="MPI" value="yes">Yes (Multiple CPUs)<br>
 <INPUT type="radio" name="MPI" value="no" checked>No (Single CPU)<br>
-</td><td valign="top">
+</td>
+-->
+<td valign="top">
 <h3>Simulate X-Ray acquistion</h3>
 Geometry:<br>
 <INPUT type="radio" name="PJ_Geometry" value="parallel" checked>Parallel<br>
 <INPUT type="radio" name="PJ_Geometry" value="equiangular">Equiangular<br>
 <INPUT type="radio" name="PJ_Geometry" value="equilinear">Equilinear<br>
 <p>
-Number of Detectors: <input type="text" name="PJ_NDet" size="4" value="367"><p>
-Number of Rotations: <input type="text" name="PJ_NRot" size="4" value="320"><p>
+Number of Detectors: <input type="text" name="PJ_NDet" size="4" value="200"><p>
+Number of Rotations: <input type="text" name="PJ_NRot" size="4" value="200"><p>
 Focal Length: <input type="text" name="PJ_FocalLength" size="4" value="2.0"><p>
+Scan Ratio: <input type="text" name="PJ_ScanRatio" size="4" value="1.1"><p>
 Number of Rays<br>(samples) per detector: <input type="text" name="PJ_NRay" size="2" value="3"><p>
-Rotation Angle<br>as a multiple of PI: <input type="text" name="PJ_RotAngle" size="3" value="1.0"><br>
+Rotation Angle<br>as a fraction of a cricle: <input type="text" name="PJ_RotAngle" size="3" value="1.0"><br>
 </td><td valign="top">
 <H3>Image Reconstruction</H3>
 Filter Method:<br>
 <input type="radio" name="IR_FilterMethod" value="convolution" checked>Convolution<br>
 <input type="radio" name="IR_FilterMethod" value="fourier">Fourier<br>
 <input type="radio" name="IR_FilterMethod" value="fourier_table">Fourier Table<br>
-<input type="radio" name="IR_FilterMethod" value="fft">FFT (Fast Fourier)<br>
-<input type="radio" name="IR_FilterMethod" value="fftw">FFTW (Fastest Fourier Library)<br>
-<input type="radio" name="IR_FilterMethod" value="rfftw">RFFTW (Real-mode Fastest Fourier Library)<br>
+<input type="radio" name="IR_FilterMethod" value="rfftw">FFT<br>
 <p>
-Zeropad: (frequency-base filtering only)<br>
-<input type="text" name="IR_Zeropad" size="1" value="2"><br>
+Zeropad (applies to<br> Fourier & FFT):
+<input type="text" name="IR_Zeropad" size="1" value="1"><br>
 <p>
 Interpolation:<br>
-<input type="radio" name="IR_Interp" value="linear" checked>Linear<br>
 <input type="radio" name="IR_Interp" value="nearest">Nearest Neighbor<br>
+<input type="radio" name="IR_Interp" value="linear" checked>Linear<br>
+<input type="radio" name="IR_Interp" value="cubic">Cubic<br>
+<p>
+Filter:<br>
+<input type="radio" name="IR_Filter" value="abs_bandlimit" checked>Bandlimit<br>
+<input type="radio" name="IR_Filter" value="abs_cos">Cosine<br>
+<input type="radio" name="IR_Filter" value="abs_hamming">Hamming<br>
+<input type="radio" name="IR_Filter" value="abs_hanning">Hanning<br>
+<input type="radio" name="IR_Filter" value="abs_sinc">Sinc<br>
+<input type="radio" name="IR_Filter" value="shepp">Shepp-Logan
 <p>
-Filter for Backprojection:<br>
-<input type="radio" name="IR_Filter" value="abs_bandlimit" checked>Abs*Bandlimit<br>
-<input type="radio" name="IR_Filter" value="abs_cos">Abs*Cosine<br>
-<input type="radio" name="IR_Filter" value="abs_hamming">Abs*Hamming<br>
-<input type="radio" name="IR_Filter" value="shepp">Shepp-Logan<br>
-<input type="radio" name="IR_Filter" value="bandlimit">Bandlimit<br>
-<input type="radio" name="IR_Filter" value="hamming">Hamming<br>
-<input type="radio" name="IR_Filter" value="cos">Cosine<br>
-<input type="radio" name="IR_Filter" value="sinc">Sinc<br>
-<input type="radio" name="IR_Filter" value="triangle">Triangle<p>
 Hamming alpha factor (0.0-1.0): <input type="text" name="IR_Filter_Param" size="3" value="1.0"><p>
 Backprojection Method:<br>
 <input type="radio" name="IR_Backproj" value="trig">Trig<br>
 <input type="radio" name="IR_Backproj" value="table">Trig Table<br>
 <input type="radio" name="IR_Backproj" value="diff">Difference<br>
-<input type="radio" name="IR_Backproj" value="idiff" checked>Difference Highly-optimized (integer math)<br>
+<input type="radio" name="IR_Backproj" value="idiff" checked>Difference (integer math)<br>
 </td></tr></table>
 <p>
-<input type="submit" value="Simulate"><input type="reset" value="clear fields">
+<input type="submit" value="Simulate"><input type="reset" value="Clear fields">
 </form>
 </BODY>
 </HTML>
Computed Tomography simulator and reconstructor