Source
88
88
imagename.model already exists, deconvolve will exit with an error.
89
89
90
90
A common usage is to set this parameter equal to a single dish image
91
91
92
92
Negative components in the model image will be included as is.
93
93
94
94
</description>
95
95
<value type="string"/>
96
96
</param>
97
97
98
-
<!-- TODO in CAS-13570: fix mtmfs description and allow mtmfs value once test_multirun_mtmfs3x passes -->
99
98
<!-- TODO in CAS-13570: fix asp description and allow asp value once asp is working -->
100
99
<param type="string" name="deconvolver">
101
100
<shortdescription>Minor cycle algorithm (hogbom,clark,multiscale,mem,clarkstokes)</shortdescription>
102
101
<!-- <shortdescription>Minor cycle algorithm (hogbom,clark,multiscale,mtmfs,mem,clarkstokes,asp)</shortdescription> -->
103
102
<description>Name of minor cycle algorithm (hogbom,clark,multiscale,mem,clarkstokes)
104
103
105
104
Each of the following algorithms operate on residual images and psfs
106
105
from the gridder and produce output model and restored images.
107
106
Minor cycles stop when threshold or niter are reached. For all
108
107
methods, components are picked from the entire extent of the image or
145
144
146
145
mem : Maximum Entropy Method [Cornwell and Evans, 1985]
147
146
- Iteratively solve for values at all individual pixels via the
148
147
MEM method. It minimizes an objective function of
149
148
chi-square plus entropy (here, a measure of difference
150
149
between the current model and a flat prior model).
151
150
152
151
(Note : This MEM implementation is not very robust.
153
152
Improvements will be made in the future.)
154
153
155
-
mtmfs : This deconvolver currently has incorrect end-of-minor-cycle residual
156
-
calculations and is therefore disabled. Please choose a different
157
-
deconvolver.
158
-
154
+
mtmfs : Multi-term (Multi Scale) Multi-Frequency Synthesis [Rau and Cornwell, 2011]
155
+
- Smooth each Taylor residual image to multiple scale sizes
156
+
- Solve a NTxNT system of equations per scale size to compute
157
+
Taylor coefficients for components at all locations
158
+
- Compute gradient chi-square and pick the Taylor coefficients
159
+
and scale size at the location with maximum reduction in
160
+
chi-square
161
+
- Add multi-scale components to each Taylor-coefficient
162
+
model image
163
+
- Subtract scaled,smoothed,shifted PSF (within a small patch size
164
+
per scale) from all smoothed Taylor residual images
165
+
- Repeat from step 2
166
+
159
167
asp : This deconvolver is not currently considered reliable when used with
160
168
task deconvolve and is therefore disabled. Please choose a different
161
169
deconvolver.
162
170
163
171
164
172
165
173
</description>
166
174
<value type="string">hogbom</value>
167
175
<allowed kind="enum">
168
176
<value>hogbom</value>
169
177
<value>clark</value>
170
178
<value>clarkstokes</value>
171
179
<value>clark_exp</value>
172
180
<value>clarkstokes_exp</value>
173
181
<value>multiscale</value>
174
-
<!-- <value>mtmfs</value> -->
182
+
<value>mtmfs</value>
175
183
<value>mem</value>
176
184
<!-- <value>asp</value> -->
177
185
178
186
</allowed>
179
187
</param>
180
188
181
189
<param type="any" name="scales" subparam="true">
182
190
<shortdescription>List of scale sizes (in pixels) for multi-scale algorithms</shortdescription>
183
191
<description>List of scale sizes (in pixels) for multi-scale and mtmfs algorithms.
184
192
This set of scale sizes should represent the sizes
194
202
smaller than the image (or mask) size and smaller than or
195
203
comparable to the scale corresponding to the lowest measured
196
204
spatial frequency (as a scale size much larger than what the
197
205
instrument is sensitive to is unconstrained by the data making
198
206
it harder to recovery from errors during the minor cycle).
199
207
</description>
200
208
<type>intVec</type><type>doubleVec</type>
201
209
<value type="intVec"/>
202
210
</param>
203
211
204
-
<!-- TODO in CAS-13570: uncomment once test_multirun_mtmfs3x passes -->
205
-
<!-- <param subparam="true" name="nterms" type="int">
212
+
<param subparam="true" name="nterms" type="int">
206
213
<shortdescription>Number of Taylor coefficients in the spectral model</shortdescription>
207
214
<description>Number of Taylor coefficients in the spectral model
208
215
209
216
- nterms=1 : Assume flat spectrum source
210
217
- nterms=2 : Spectrum is a straight line with a slope
211
218
- nterms=N : A polynomial of order N-1
212
219
213
220
From a Taylor expansion of the expression of a power law, the
214
221
spectral index is derived as alpha = taylorcoeff_1 / taylorcoeff_0
215
222
242
249
from the propagation of error during the division of
243
250
two noisy numbers (alpha = xx.tt1/xx.tt0) where the
244
251
'error' on tt1 and tt0 are simply the values picked from
245
252
the corresponding residual images. The absolute value
246
253
of the error is not always accurate and it is best to interpret
247
254
the errors across the image only in a relative sense.)
248
255
249
256
250
257
</description>
251
258
<value type="int">2</value>
252
-
</param> -->
259
+
</param>
253
260
254
261
<param type="double" name="smallscalebias" subparam="true">
255
262
<shortdescription>Biases the scale selection when using multi-scale or mtmfs deconvolvers </shortdescription>
256
263
<description>A numerical control to bias the scales when using multi-scale or mtmfs algorithms.
257
264
The peak from each scale's smoothed residual is
258
265
multiplied by ( 1 - smallscalebias * scale/maxscale )
259
266
to increase or decrease the amplitude relative to other scales,
260
267
before the scale with the largest peak is chosen.
261
268
Smallscalebias can be varied between -1.0 and 1.0.
262
269
A score of 0.0 gives all scales equal weight (default).
745
752
746
753
<when param="deconvolver">
747
754
<equals type="string" value="hogbom">
748
755
</equals>
749
756
<equals type="string" value="clark">
750
757
</equals>
751
758
<equals type="string" value="multiscale">
752
759
<default param="scales"><value type="intVec"/></default>
753
760
<default param="smallscalebias"><value>0.0</value></default>
754
761
</equals>
755
-
<!-- TODO in CAS-13570: uncomment once test_multirun_mtmfs3x passes -->
756
-
<!-- <equals type="string" value="mtmfs">
762
+
<equals type="string" value="mtmfs">
757
763
<default param="scales"><value type="intVec"/></default>
758
764
<default param="nterms"><value type="int">2</value></default>
759
765
<default param="smallscalebias"><value>0.0</value></default>
760
-
</equals> -->
766
+
</equals>
761
767
<!-- TODO in CAS-13570: uncomment once asp is working -->
762
768
<!-- <equals type="string" value="asp">
763
769
<default param="fusedthreshold"><value>0.0</value></default>
764
770
<default param="largestscale"><value>-1</value></default>
765
771
</equals> -->
766
772
</when>
767
773
768
774
<when param="restoration">
769
775
<equals type="bool" value="True">
770
776
<default param="restoringbeam"><value type="stringVec"/></default>