HTML automatically generated with rman
Table of Contents
testbsc - testing Beam Smearing Correction algorithms
testbsc
[parameter=value]
testbsc tests various Beam Smearing Correction
schemes to a velocity field. Currently only Begeman’s correction is available,
but expect to try out Leung’s (2018) method in an upcoming version.
Begeman’s
(1989, Appendix B, eq.(8)
) recipe reads as follows:
v = V + b^2/n * ( Vx*Nx + Vy*Ny + N * (Vxx + Vyy) / 2)
where the correction terms (beeg=t below) have single and double derivatives.
(N,V) are the model density and velocity, (n,v) the observed/smoothed values,
in this case with a circular gaussian beam with sigma (not FWHM) of b (see
beam= below) points per pixel.
The following parameters are recognized
in any order if the keyword is also given:
- inD=
- Input Model Density image
file []
- inV=
- Input Model Velocity image file []
- ind=
- Input Observed
(smoothed) Density image file []
- inv=
- Input Observed (smoothed) Velocity
image file []
- out=
- Output file, depends on the values of diff= and beeg=.
- beam=
- Override sigma_beam/pixel ratio. Recall FWHM = 2.355 * SIGMA. By default
the beam size will be in the header.
- diff=f|t
- Compute a Difference map? (as
opposed to Computed Observed). Normally the Observed velocity field is computed
from the Model and BSC terms, but turning this option on, the input Observed
is subtracted from the computed Observed. [f]
- beeg=t|f
- Use the Begeman correction
terms. This is of course the default, but by turning this off, a simple
different of V-v (Model - Observed) can be computed, as can be seen in the
EXAMPLES below. [t]
ccdmath(1NEMO)
, rotcur(1NEMO)
, image(5NEMO)
In this example we show off the interplay of the diff= and beeg=
keywords:
testbsc N.ccd V.ccd n.ccd v.ccd map1 diff=f beeg=f
testbsc N.ccd V.ccd n.ccd v.ccd map2 diff=t beeg=f
testbsc N.ccd V.ccd n.ccd v.ccd map3 diff=f beeg=t
testbsc N.ccd V.ccd n.ccd v.ccd map4 diff=t beeg=t
Now the four maps will have the following meaning:
map1: V
map2: V - v
map3: V + bsc (should be close to ’v’)
map4: V + bsc - v (should be close to 0)
NEMO/src/image/rotcur
Peter Teuben
3-jul-2019 V0.4 Created/released PJT
Table of Contents