next previous
Up: On the automatic folding


9 The Mk III auto-folder

The Mk III auto-folder is based on the working hypothesis that the folding process should only use the exterior parts of the rotation curve, so that hole-cutting becomes an active part of this process. However, there is one immediate practical difficulty: the hole-cutting algorithm requires rotation curves to be already folded. The resolution of this problem is embodied in the algorithmic definition of the Mk III auto-folder, described in Sect. 9.1, and the solution arising from it applied over the PS data is shown in Fig. 5 right.
  \begin{figure}\includegraphics{fig5.eps}\end{figure} Figure 5: Comparison of Mk II and Mk III auto-folder solutions. The vertical lines indicate the positions of the A, B, C, D and E peaks in the PS solution

9.1 Algorithmic definition of the Mk III auto-folder

The Mk III auto-folder is defined by the following algorithm:

9.2 Partition of Fourier region

The application of the hole-cutting algorithm of Roscoe 1999A cuts a hole out of the rotation curve which is centred on the estimated dynamical centre of the galaxy; consequently, instead of representing a continuous curve, the reduced rotation curve data represents two separated continuous sections with a gap in between. This means that the Fourier cosine coefficients cannot be computed using (2), but must be computed by an integral of the form
$\displaystyle A(m) = {1 \over 2\,X} \int_{-X}^{-X'} f(x) \cos (m \pi x/X) ~{\rm d}x ~+~$      
$\displaystyle {1 \over 2\, X} \int_{X'}^X f(x) \cos (m \pi x/X) ~{\rm d}x,$     (3)

where X' is some positive number representing the boundaries of the cut-out hole. In practice, of course, the hole-cutting strategy will generally remove a section which is non-symmetric about the current estimate of $O_{{\rm dyn}}$, and so the data on one of the two sections of the reduced rotation curve must be interpolated/extrapolated to ensure that the hole which has been effectively cut out of the rotation curve is centred exactly on the current estimate of $O_{{\rm dyn}}$. Once the removed section is centred in this way, and the Fourier modes computed using (3), rather than (2), the process is as before.
next previous
Up: On the automatic folding

Copyright The European Southern Observatory (ESO)