Skip to content

Commit

Permalink
pubs: Ensure similar configs are used for 2D cf and fw, and describe …
Browse files Browse the repository at this point in the history
…them
  • Loading branch information
David Berry committed Jul 18, 2014
1 parent abddeb0 commit 856ac36
Show file tree
Hide file tree
Showing 9 changed files with 90 additions and 20 deletions.
27 changes: 25 additions & 2 deletions pubs/fellwalker/comp4
Expand Up @@ -7,6 +7,27 @@ set peak_value = 100
echo "# rms deltat nc gain sum sum_target true_data mask" > cf.asc
echo "# rms maxjump mindip nc gain sum sum_target true_data mask" > fw.asc

# Create the clumpfind config.
echo "ClumpFind.Allowedge=1" > cf.conf
echo "ClumpFind.Fwhmbeam=0" >> cf.conf
echo "ClumpFind.Idlalg=1" >> cf.conf
echo "ClumpFind.Maxbad=1.0" >> cf.conf
echo "ClumpFind.Minpix=16" >> cf.conf
echo "ClumpFind.Naxis=2" >> cf.conf
echo "ClumpFind.Tlow=2*RMS" >> cf.conf
echo "ClumpFind.VeloRes=0" >> cf.conf

# Create the fellWalker config.
echo "FellWalker.AllowEdge=1" > fw.conf
echo "FellWalker.CleanIter=1" >> fw.conf
echo "FellWalker.FlatSlope=1*RMS" >> fw.conf
echo "FellWalker.FwhmBeam=0" >> fw.conf
echo "FellWalker.MaxBad=0.05" >> fw.conf
echo "FellWalker.MinHeight=2*RMS" >> fw.conf
echo "FellWalker.MinPix=16" >> fw.conf
echo "FellWalker.Noise=2*RMS" >> fw.conf
echo "FellWalker.VeloRes=0" >> fw.conf


$CUPID_DIR/makeclumps angle=\[0,0\] beamfwhm=0 deconv=no \
fwhm1=\[$fwhm,0\] fwhm2=\[$fwhm,0\] lbnd=\[1,1\] ubnd=\[100,100\] \
Expand Down Expand Up @@ -44,7 +65,8 @@ foreach rms (2 4 6 8 10 12 14 16)
set cf_mask = "comp4_cf_${rms}_${pv}"
$CUPID_DIR/findclumps in=$truedata deconv=no method=clumpfind \
out=$cf_mask outcat=comp4_cf.fit rms=$rms \
config="'^idl_conf,DeltaT=$deltat*RMS'"
config="'^cf.conf,DeltaT=$deltat*RMS'" \
backoff=yes wcspar=no perspectrum=no

set text = `stilts tpipe in=comp4_cf.fit cmd='keepcols SUM;stats median' | grep -P "\| \d+\.\d+ \|" | grep -o -P "\d+\.\d+"`

Expand All @@ -65,7 +87,8 @@ foreach rms (2 4 6 8 10 12 14 16)
set fw_mask = "comp4_fw_${rms}_${maxjump}_${pv}"
$CUPID_DIR/findclumps in=$truedata deconv=no method=fellwalker \
out=$fw_mask outcat=comp4_fw.fit rms=$rms \
config="'MaxJump=$maxjump,MinDip=$mindip*RMS'"
backoff=yes wcspar=no perspectrum=no \
config="'^fw.conf,MaxJump=$maxjump,MinDip=$mindip*RMS'"

set text = `stilts tpipe in=comp4_fw.fit cmd='keepcols SUM;stats median' | grep -P "\| \d+\.\d+ \|" | grep -o -P "\d+\.\d+"`

Expand Down
Binary file modified pubs/fellwalker/comp4_cf.pdf
Binary file not shown.
Binary file modified pubs/fellwalker/comp4_fw_0.pdf
Binary file not shown.
Binary file modified pubs/fellwalker/comp4_fw_1.pdf
Binary file not shown.
Binary file modified pubs/fellwalker/comp4_fw_2.pdf
Binary file not shown.
Binary file modified pubs/fellwalker/comp4_fw_3.pdf
Binary file not shown.
Binary file modified pubs/fellwalker/comp4_fw_4.pdf
Binary file not shown.
Binary file modified pubs/fellwalker/comp4_fw_5.pdf
Binary file not shown.
83 changes: 65 additions & 18 deletions pubs/fellwalker/fellwalker.tex
Expand Up @@ -530,7 +530,7 @@ \subsection{Characterising Each Clump}


\section{Comparing FellWalker and CLUMPFIND}

\label{sec:compare}
\cite{2010Watson} made a detailed comparison of the performance of several
different 2-dimensional clump finding algorithms, including FellWalker and CLUMPFIND.
This concluded that, under the conditions used in the study, FellWalker
Expand Down Expand Up @@ -626,20 +626,11 @@ \subsection{The Clump-finding Parameters}
package follows the description of the algorithm contained in
\cite{1994Williams}. The IDL version of CLUMPFIND distributed by Williams
includes some enhancements to the published algorithm. These are also
available in the CUPID version, but are disabled by default. They have,
however been enabled for the purposes of the current comparison, using
the following extra CUPID parameter settings:

\smallskip
\begin{tabular}{c c c}
ClumpFind.IDLAlg & = & 1 \\
ClumpFind.FwhmBeam & = & 0 \\
ClumpFind.MaxBad & = & 1.0 \\
ClumpFind.MinPix & = & 20 \\
ClumpFind.VeloRes & = & 0 \\
ClumpFind.AllowEdge & = & 1 \\
\end{tabular}
\smallskip
available in the CUPID version, but are disabled by default. They were,
however, enabled for the purposes of the current comparison.

Full lists of parameters used for these tests are given in
appendix~\ref{app:configs}.

\subsection{Measuring Performance}
In each test, the artificial data consists of a collection of identical
Expand Down Expand Up @@ -862,10 +853,66 @@ \section{Acknowledgements}

%% The Appendices part is started with the command \appendix;
%% appendix sections are then done as normal sections
%% \appendix
\appendix

\section{Configuration Parameter Values}
\label{app:configs}
The following two sections shows the CUPID commands that were
used to perform the 2-dimensional tests described in section
section~\ref{sec:compare}. Also shown are the contents of the text files
that specify the values used for the additional configuration parameters
required by CUPID. See the CUPID documentation \citep{SUN255} for
descriptions of these commands and parameters. In both cases,
\texttt{data.sdf} is the data file containing the artificial clumps,
\texttt{mask.sdf} is the generated clump pixel mask and \texttt{cat.fit}
is the generate clump catalogue. For each specific test, different
numerical values were assigned to the variables \texttt{\$rms},
\texttt{\$maxjump}, \texttt{\$mindip} and \texttt{\$deltat}, as described
in section~\ref{sec:compare}.
\subsection{FellWalker}
{\small
\begin{verbatim}
% findclumps in=data.sdf deconv=no method=fellwalker \
out=mask.sdf outcat=cat.fit rms=$rms \
config="'^fw.conf,MaxJump=$maxjump, \
MinDip=$mindip*RMS'" \
wcspar=no backoff=yes perspectrum=no
% cat fw.conf
FellWalker.AllowEdge=1
FellWalker.CleanIter=1
FellWalker.FlatSlope=1*RMS
FellWalker.FwhmBeam=0
FellWalker.MaxBad=0.05
FellWalker.MinHeight=2*RMS
FellWalker.MinPix=16
FellWalker.Noise=2*RMS
FellWalker.VeloRes=0
\end{verbatim}
}
\subsection{ClumpFind}
{\small
\begin{verbatim}
% findclumps in=data.sdf deconv=no method=clumpfind \
out=mask.sdf outcat=cat.fit rms=$rms \
config="'^cf.conf,DeltaT=$deltat*RMS'" \
wcspar=no backoff=yes perspectrum=no
% cat cf.conf
ClumpFind.Allowedge=1
ClumpFind.Fwhmbeam=0
ClumpFind.Idlalg=1
ClumpFind.Maxbad=1.0
ClumpFind.Minpix=16
ClumpFind.Naxis=2
ClumpFind.Tlow=2*RMS
ClumpFind.VeloRes=0
\end{verbatim}
}
For the 3-dimensional tests, the value used for \texttt{ClumpFind.Naxis}
was changed from 2 to 3.


%% \section{}
%% \label{}

%% References
%%
Expand Down

0 comments on commit 856ac36

Please sign in to comment.