[529] | 1 | Duchamp version 1.1.7 (2009/01/16) |
---|
[520] | 2 | ---------------------------------- |
---|
| 3 | |
---|
| 4 | Improvements: |
---|
| 5 | |
---|
[522] | 6 | * The mask file that can be created with flagMask can now specify |
---|
| 7 | which object a given detected pixel belongs to, by taking the object |
---|
| 8 | ID as its value. This applies if you set flagMaskWithObjectNum=true |
---|
| 9 | (it defaults to false, which means each detected pixel has the value |
---|
| 10 | 1). Non-detected pixels still have the value 0. |
---|
[520] | 11 | |
---|
[536] | 12 | * The names of the output FITS files are now able to be specified |
---|
| 13 | directly by input parameters. These will override the default names |
---|
| 14 | that Duchamp provides. |
---|
| 15 | |
---|
[520] | 16 | * The Karma annotation file now specifies a font - previously it did |
---|
| 17 | not and so was at the mercy of whatever was specified in other |
---|
| 18 | annotation files being viewed at the same time. Please let me know |
---|
| 19 | if the font I've chosen (lucidasans-12) is not suitable. |
---|
| 20 | |
---|
| 21 | * The default flag, used in the absence of 'E', 'S' or 'N', is now '-' |
---|
| 22 | rather than blank. This will ensure that each column in the results |
---|
| 23 | file has something in it. |
---|
| 24 | |
---|
[522] | 25 | * The configure script is somewhat more robust to non-standard |
---|
| 26 | systems, particularly when the X11 libraries are in unusual |
---|
| 27 | places. It can also deal with the slightly different way that |
---|
| 28 | wcslib-4.3 installs itself. |
---|
| 29 | |
---|
[520] | 30 | * The installation part of the Makefile was altered so that the |
---|
| 31 | symbolic link Duchamp that points to Duchamp-1.1.7 does not have the |
---|
| 32 | absolute path encoded in it. This will enable the files to be moved |
---|
| 33 | around at a later date. |
---|
| 34 | |
---|
| 35 | Bug Fixes: |
---|
| 36 | |
---|
| 37 | * Fixed a bug that caused problems when the flux units were less than |
---|
| 38 | 5 characters long (e.g. 'K') |
---|
| 39 | |
---|
| 40 | * The new velocity width values were not getting proper headers in the |
---|
| 41 | VOTable output. This has been remedied. |
---|
| 42 | |
---|
| 43 | * Minor improvements to help non-pgplot use |
---|
| 44 | |
---|
| 45 | * Minor fix in WCS functionality. |
---|
| 46 | |
---|
| 47 | * Minor fixes to the FitsHeader class. |
---|
| 48 | |
---|
| 49 | * Improved FITS data reading, that doesn't require a nullarray |
---|
| 50 | array. This will reduce memory usage. |
---|
| 51 | |
---|
| 52 | * Improved subsection parsing, particularly for non-standard WCS |
---|
| 53 | systems. |
---|
| 54 | |
---|
| 55 | * Improving the memory allocation of the Cube class. |
---|
| 56 | |
---|
| 57 | * Fixed a bug that meant the growthThreshold values were sometimes |
---|
| 58 | being lost. |
---|
| 59 | |
---|
| 60 | * Improving the way the growing and merging is done, particularly in |
---|
| 61 | the case of non-adjacent merging. |
---|
| 62 | |
---|
[529] | 63 | * Fixed a scaling bug in the GaussSmooth function. |
---|
| 64 | |
---|
[520] | 65 | Other: |
---|
| 66 | |
---|
| 67 | * Numerous changes to the code that do not affect how Duchamp itself |
---|
| 68 | runs, but were necessary for code development with the Duchamp |
---|
| 69 | libraries (primarily as a part of the ASKAP code development). |
---|
| 70 | |
---|
[522] | 71 | |
---|
[482] | 72 | Duchamp version 1.1.6 (2008/09/08) |
---|
[480] | 73 | ---------------------------------- |
---|
| 74 | |
---|
| 75 | New features: |
---|
| 76 | |
---|
| 77 | * The velocity widths can now be quoted in two additional ways: the |
---|
| 78 | full width at 50% of the peak flux (ie. the FWHM) and the full width |
---|
| 79 | at 20% of the peak flux (w_50 and w_20 respectively). The w_50 |
---|
| 80 | parameter is now quoted in the results printed to the screen |
---|
| 81 | (instead of w_VEL), but all three are written to the results file. |
---|
| 82 | |
---|
| 83 | * The facility now exists to bypass the searching step and go straight |
---|
| 84 | to the output stage using a list of previously detected objects. The |
---|
| 85 | list is extracted from a previously-created log file. The user can |
---|
| 86 | change the way the spectra are plotted (e.g. choose to plot |
---|
| 87 | integrated rather than peak spectra), and select a subset of the |
---|
| 88 | detected objects to be printed to individual postscript files. Use |
---|
| 89 | the input parameter usePrevious (set to true) to use this facility, |
---|
| 90 | and provide a list of objects using the objectList parameter |
---|
| 91 | (eg. 1,2,4,7-10,15) |
---|
| 92 | |
---|
[482] | 93 | * A new parameter "growthThreshold" has been created, that goes with |
---|
| 94 | "threshold". If you manually specify a threshold, and you wish to |
---|
| 95 | grow detections, you need to manually specify a growthThreshold (in |
---|
| 96 | flux units). |
---|
| 97 | |
---|
[480] | 98 | Bug fixes: |
---|
| 99 | |
---|
[482] | 100 | * Cubes that do not have a spectral axis (such as Stokes cubes) were |
---|
| 101 | causing memory leaks and crashing, due to some assumptions of the |
---|
| 102 | existence of the spectral axis. These bugs have been fixed, so that |
---|
| 103 | Duchamp will run on such cubes, although the WCS transforms may not |
---|
| 104 | always work appropriately. Please let me know via bug reports of |
---|
| 105 | difficulties you encounter! |
---|
[480] | 106 | |
---|
[482] | 107 | * The growing function has had a few bugs fixed and is now more |
---|
| 108 | efficient. |
---|
| 109 | |
---|
| 110 | * The FDR method was previously overriding any manually-specified |
---|
| 111 | threshold. This has been fixed so that if the user specifies |
---|
| 112 | threshold=X, then the FDR method is not called. |
---|
| 113 | |
---|
| 114 | * Some undefined things were happening with textual output when no |
---|
| 115 | objects were detected. This has been fixed. |
---|
| 116 | |
---|
| 117 | * Memory leaks caused by FITS WCS irregularities have been fixed |
---|
| 118 | (specifically, if a FITS header had NAXIS=2 but more than two WCS |
---|
| 119 | axes were defined). |
---|
| 120 | |
---|
| 121 | |
---|
| 122 | |
---|
[438] | 123 | Duchamp version 1.1.5 (2008/03/20) |
---|
| 124 | ---------------------------------- |
---|
| 125 | |
---|
| 126 | Bug fixes: |
---|
| 127 | * A major bug was discovered in the calculation of integrated |
---|
| 128 | fluxes. The flux was incorrectly being weighted by the spatial size |
---|
| 129 | of the detection, resulting in a (often considerable) overestimation |
---|
[447] | 130 | of the integrated flux. This has been corrected. |
---|
[438] | 131 | |
---|
[447] | 132 | * The beam size calculation was also incorrect, being out by a factor |
---|
| 133 | of 1/ln(2). |
---|
[438] | 134 | |
---|
[447] | 135 | * The "make install" command previously failed to work properly if the |
---|
| 136 | library had not been created. This has been fixed so that the |
---|
| 137 | installation works even if there is no library file. |
---|
| 138 | |
---|
| 139 | New Features: |
---|
| 140 | * You can now provide a new flux unit string, so that the pixel fluxes |
---|
| 141 | will be converted to these new fluxes before any calculations are |
---|
| 142 | done. This is done with the newFluxUnits input parameter. The new |
---|
| 143 | units should be directly convertible from the old ones (given by the |
---|
| 144 | BUNIT header -- e.g. if BUNIT = Jy/beam, you can give newFluxUnits |
---|
| 145 | as mJy/beam, but not mJy). |
---|
| 146 | |
---|
| 147 | * You can now specify the precision used in the output files for flux, |
---|
| 148 | velocity and S/N ratio columns, via new input parameters precFlux, |
---|
| 149 | precVel and precSNR. These will apply to the VOTable output as well |
---|
| 150 | as the text table and information on the spectral plots. |
---|
| 151 | |
---|
| 152 | * The Karma annotation files previously provided just a circle at the |
---|
| 153 | spatial location of the detected object. You now have the option to |
---|
| 154 | draw a border around the detected pixels (in a similar manner to the |
---|
| 155 | PGPLOT moment maps). Karma does not do this perfectly, so the lines |
---|
| 156 | may not line up exactly with the pixel borders, so the option |
---|
| 157 | remains to keep using the circles. This is controlled by the |
---|
| 158 | parameter "annotationType", with possible options being "borders" |
---|
| 159 | (the default) or "circles". |
---|
| 160 | |
---|
| 161 | * There is a new function getMetadata() that loads the header and WCS |
---|
| 162 | data, but no pixel information. This is only of interest if you are |
---|
| 163 | doing code development with the Duchamp library -- it does not |
---|
| 164 | affect the running of Duchamp itself. |
---|
| 165 | |
---|
[425] | 166 | Duchamp version 1.1.4 (2008/02/15) |
---|
| 167 | ---------------------------------- |
---|
| 168 | |
---|
| 169 | Most of the changes incorporated in this update result from ASKAP code |
---|
| 170 | development (that makes use of some of the Duchamp functionality), |
---|
| 171 | but do not actually affect the way Duchamp runs (other than, in a |
---|
| 172 | couple of cases, making it a bit faster). Changes that are |
---|
| 173 | relevant for regular users of Duchamp are as follows: |
---|
| 174 | |
---|
| 175 | * A new option has been made available (thanks to Mary Putman et al |
---|
| 176 | for the suggestion!) to save the spectra of all detected objects in |
---|
| 177 | a text file. This has the form of one column for the spectral |
---|
| 178 | coordinates and one column of fluxes for each object. This is |
---|
| 179 | enabled by the input parameters flagTextSpectra and spectraTextFile. |
---|
| 180 | |
---|
| 181 | * Slight changes to the output include he formatting of the way the |
---|
| 182 | detected objects are written to the log file. This meant the |
---|
| 183 | verification files were also updated. Also the cube statistics are |
---|
| 184 | not printed to the results file if a user threshold is specified (as |
---|
| 185 | they are not calculated!). |
---|
| 186 | |
---|
| 187 | * The determination of beam size has been made more reliable, as has |
---|
| 188 | the BUNIT determination. |
---|
| 189 | |
---|
| 190 | * Fixed a bug that meant source merging was not done for 2D images |
---|
| 191 | (not usually a problem, as the source-finding takes care of this, |
---|
| 192 | but came up in parallelisation development). |
---|
| 193 | |
---|
[406] | 194 | Duchamp version 1.1.3 (2008/01/08) |
---|
| 195 | ---------------------------------- |
---|
[404] | 196 | |
---|
| 197 | A minor update that just improves the build process. The location of |
---|
| 198 | the wcslib include files was not specified correctly, so that some |
---|
| 199 | machines were unable to find them. The configure script has been fixed |
---|
| 200 | so that the compilation should work in all cases. |
---|
| 201 | |
---|
[405] | 202 | Some minor changes have also been made to the code in the PixelMap |
---|
[412] | 203 | directory, and to param.cc, to fix warnings or errors that may appear |
---|
| 204 | in compilation. |
---|
[404] | 205 | |
---|
[405] | 206 | However, there are no changes to the actual workings of the code in |
---|
| 207 | this minor version, so if you have a working version of Duchamp 1.1.2 |
---|
| 208 | there is probably no reason to update. |
---|
| 209 | |
---|
[389] | 210 | Duchamp version 1.1.2 (2007/12/07) |
---|
| 211 | ---------------------------------- |
---|
[349] | 212 | |
---|
[363] | 213 | Compilation: |
---|
| 214 | |
---|
[349] | 215 | * The compilation had problems working correctly when the cfitsio and |
---|
[373] | 216 | wcs libraries were not in standard locations. The configure script |
---|
[349] | 217 | has been tweaked to make sure the libraries are found. |
---|
| 218 | |
---|
| 219 | * Problems with compiling in the absence of pgplot have been fixed, |
---|
| 220 | and it should now compile fine. |
---|
| 221 | |
---|
[362] | 222 | * Improved the configure script for use with Mac OS X and gfortran. |
---|
[349] | 223 | |
---|
[363] | 224 | |
---|
| 225 | Bugfixing: |
---|
| 226 | |
---|
[362] | 227 | * A major problem with the text-based output was fixed. The widths of |
---|
| 228 | the columns were not being calculated correctly, and were running |
---|
[373] | 229 | together. This should no longer be the case. |
---|
[349] | 230 | |
---|
[374] | 231 | * There was a bug in the way Duchamp interpreted spectral axes, |
---|
| 232 | particularly AIPS types such as FELO-HEL, leading to incorrect |
---|
| 233 | scaling of the velocity axis. This has been fixed, so that |
---|
| 234 | non-standard types are now dealt with appropriately and the velocity |
---|
| 235 | scaling should be accurate. |
---|
[362] | 236 | |
---|
[376] | 237 | * Another bug meant the beam correction was always applied to the |
---|
| 238 | integrated spectra, regardless of whether the flux units involved |
---|
| 239 | the beam (e.g. Jy/beam). This has been fixed. Note that the values |
---|
| 240 | of the integrated flux presented in the results table were always |
---|
| 241 | calculated correctly. |
---|
| 242 | |
---|
[389] | 243 | * Other minor changes to improve the speed and efficiency, especially |
---|
| 244 | of the merging code. |
---|
| 245 | |
---|
[363] | 246 | New user parameters: |
---|
| 247 | |
---|
[381] | 248 | * There is now the possibility to save a FITS file that contains a |
---|
| 249 | mask array, indicating where the detected objects are. This is |
---|
| 250 | governed by the flagOutputMask parameter. |
---|
| 251 | |
---|
[363] | 252 | * The header information that is printed at the start of the results |
---|
| 253 | file -- that is, the list of parameters, the statistics of the cube, |
---|
| 254 | and the number of detections -- can now be written to a separate |
---|
| 255 | file. Set the flagSeparateHeader parameter to true or 1 and it will |
---|
| 256 | go to the file indicated by the headerFile parameter (the default is |
---|
| 257 | duchamp-Results.hdr). |
---|
| 258 | |
---|
[362] | 259 | * A "scaleMax" parameter has been added, so that the wavelet |
---|
| 260 | reconstruction can include only a specified range of scales (for |
---|
[373] | 261 | instance, if one wants to search just for narrow features such as |
---|
[362] | 262 | absorption lines). |
---|
| 263 | |
---|
[389] | 264 | * A "flagRobustStats" parameter has been added. This allows the user |
---|
| 265 | to elect to use robust statistics (median and median absolute |
---|
| 266 | deviation from the median) to estimate the noise parameters. If set |
---|
| 267 | to false, the mean and standard deviation will be used. The default |
---|
| 268 | value is true. |
---|
[362] | 269 | |
---|
[369] | 270 | New Feature: |
---|
[362] | 271 | |
---|
[369] | 272 | * When Duchamp is run on a 2D image (i.e. with no spectral output), it |
---|
| 273 | now produces a postscript file that shows information each |
---|
| 274 | individual source plus the image cutout of the source. This is |
---|
| 275 | similar to the spectral output file, but without the spectra. The |
---|
| 276 | filename is given by the spectraFile parameter. |
---|
| 277 | |
---|
[381] | 278 | * Ability to save a FITS file containing a mask array, with pixel |
---|
| 279 | values of 1 for pixels in a detected object, and 0 elsewhere. |
---|
[369] | 280 | |
---|
[381] | 281 | |
---|
[328] | 282 | Duchamp version 1.1.1 (2007/07/20) |
---|
| 283 | ------------------------------------ |
---|
| 284 | |
---|
[329] | 285 | Minor update to the release, with just a few changes: |
---|
[328] | 286 | |
---|
| 287 | * The verification script was generating diffs due to differences in |
---|
| 288 | the precision of some calculations when done on different machines. |
---|
| 289 | One of these was the integrated flux, and this has been corrected so |
---|
| 290 | that the same value should appear on all machines. |
---|
| 291 | The verification script has also changed so that just the number of |
---|
| 292 | sources and their positions are tested, as well as the |
---|
| 293 | logfiles. This should avoid the reporting of spurious differences. |
---|
| 294 | |
---|
| 295 | * The grey-scale used in the image cutouts (in the spectral plots) has |
---|
| 296 | been inverted. This makes it consistent with the scale used in the |
---|
| 297 | maps, and hopefully easier to see when printing. |
---|
| 298 | |
---|
[338] | 299 | * The configure/make system has been altered to make it easier to |
---|
| 300 | install, and to specify different external libraries. In particular, |
---|
| 301 | it is now not necessary to have PGPLOT installed on your system to |
---|
| 302 | run Duchamp. In the absence of PGPLOT, Duchamp will run as normal, |
---|
| 303 | just without any graphical output. Some source code needed to be |
---|
[344] | 304 | restructured to make this work. Please read the README file, or the |
---|
| 305 | Installation appendix in the User's Guide, for complete details. |
---|
[328] | 306 | |
---|
[329] | 307 | * The time of completion is now written to the log file, so that the |
---|
| 308 | user can have some idea of how long the program took. |
---|
[328] | 309 | |
---|
| 310 | |
---|
[310] | 311 | Duchamp version 1.1 (2007/05/17) |
---|
[259] | 312 | ------------------------------------ |
---|
| 313 | |
---|
| 314 | Changes made since the last numbered release: |
---|
| 315 | |
---|
[295] | 316 | EFFICIENCY OF EXECUTION: |
---|
| 317 | |
---|
[259] | 318 | * Major change to the way detections are stored. They are now recorded |
---|
| 319 | via run-length encoding in the x-y plane, and stored as collections |
---|
| 320 | of channel-maps. This improves the memory usage and the speed of |
---|
[298] | 321 | execution (particularly for routines such as the merging |
---|
| 322 | routines). The memory usage has also been improved by changing the |
---|
| 323 | way the list of detected objects is stored. |
---|
[259] | 324 | |
---|
[295] | 325 | SMOOTHING: |
---|
| 326 | |
---|
[279] | 327 | * Addition of ability to spatially smooth the cube before |
---|
| 328 | searching. This is done with a Gaussian kernel, the size and shape |
---|
[287] | 329 | of which can be specified by the user, using the parameters kernMaj |
---|
| 330 | & kernMin (which are the FWHM for the major & minor axes -- if |
---|
| 331 | kernMin is not given it is assumed to match kernMaj) and kernPA (the |
---|
| 332 | position angle of the major axis). There is a new parameter |
---|
[279] | 333 | smoothType that distinguishes this method from the spectral (Hanning) |
---|
| 334 | smoothing already in use. |
---|
| 335 | |
---|
[295] | 336 | DETECTIONS: |
---|
| 337 | |
---|
[264] | 338 | * Changed the way detection is done -- now only search in 2D images, |
---|
| 339 | keeping all detections (ie. no minimum-pixel requirement). This |
---|
| 340 | should be the same as doing both the 1D and 2D searches with the |
---|
| 341 | minimum pixel requirement in, but simpler. The rejection is done at |
---|
[295] | 342 | the end after all the merging and growing has been done. This leads |
---|
| 343 | to a larger number of intermediate detections (ie. the number |
---|
| 344 | before merging), which may increase the run-time slightly, although |
---|
| 345 | this is offset by the changes to the efficiency (see above). |
---|
[264] | 346 | |
---|
| 347 | * Increased the range over which objects can be grown to match the |
---|
| 348 | threshSpatial & threshVelocity parameters, as well as rejecting new |
---|
| 349 | pixels based on the BLANK and MW tests. |
---|
| 350 | |
---|
[295] | 351 | INPUTS: |
---|
[259] | 352 | |
---|
[287] | 353 | * There has been a subtle change in the way BLANK pixels are dealt |
---|
| 354 | with. To enable trimming of BLANK pixels, users should use the new |
---|
| 355 | parameter flagTrim. The parameters flagBlankPix and blankPixVal are |
---|
[295] | 356 | no longer accessible via the parameter file -- warning messages are |
---|
| 357 | given if they are provided in the parameter file. |
---|
[287] | 358 | |
---|
[295] | 359 | * Added ability to specify a statsec, so that the statistics are only |
---|
| 360 | calculated from a given subsection of the cube, rather than the full |
---|
| 361 | cube. |
---|
| 362 | |
---|
| 363 | * The subsections can also specify the ability to remove a border of n |
---|
| 364 | pixels from the edges of an axis (ie. n from the start and from the |
---|
| 365 | end) -- as requested in Trac ticket #5. |
---|
| 366 | |
---|
| 367 | * Changed the default value of the reconDim parameter from 3 to 1, and |
---|
| 368 | the default value of the flagATrous parameter to false. |
---|
| 369 | |
---|
| 370 | OUTPUTS: |
---|
| 371 | |
---|
[279] | 372 | * The output now includes three different estimates of the "centre" of |
---|
| 373 | the detection: peak location, centroid (flux-weighted), or average. |
---|
| 374 | |
---|
[295] | 375 | * Some slight improvements to the user feedback during execution. |
---|
| 376 | |
---|
[282] | 377 | * Spectral output now includes the baseline, if one was fitted. The |
---|
| 378 | format of the plot headers is slightly changed, to stop information |
---|
| 379 | falling off the edge of the page. |
---|
[279] | 380 | |
---|
[295] | 381 | * Added command-line option -x to directly disable the PGPLOT X-window |
---|
| 382 | map plotting. This overrides the setting of flagXOutput in the |
---|
| 383 | parameter file. |
---|
[259] | 384 | |
---|
[295] | 385 | * Fallback spectral units (in case the FITS header fails to record |
---|
| 386 | them) changed to SPC, but code has been fixed so that any problem |
---|
| 387 | should be caught before this becomes an issue. |
---|
| 388 | |
---|
| 389 | BUGFIXES: |
---|
| 390 | |
---|
[259] | 391 | * Fixed bug (Trac ticket #4) that meant the velocity units were not |
---|
| 392 | correctly dealt with. Duchamp should be more robust for different |
---|
| 393 | velocity axis setups now. |
---|
| 394 | |
---|
| 395 | * Fixed bug (ticket #9) with the VOTable output, where the tags were |
---|
| 396 | not closing. |
---|
| 397 | |
---|
[305] | 398 | * Fixed a bug that meant the integrated fluxes were not calculated |
---|
| 399 | properly -- no allowance was made for the number of spatial |
---|
| 400 | pixels. This should also work properly for 2D images now as well. |
---|
| 401 | |
---|
[264] | 402 | * Fixed minor bug in calculating the flux threshold from the p-value |
---|
| 403 | threshold when the FDR method is used. This did not affect the |
---|
| 404 | actual thresholding, merely the quoted flux threshold in the results |
---|
| 405 | file. Also, the MW channels are kept out of the p-value threshold |
---|
| 406 | determination if they are being flagged. |
---|
| 407 | |
---|
[259] | 408 | |
---|
| 409 | |
---|
[211] | 410 | Duchamp version 1.0.7 (2006/11/22) |
---|
| 411 | ---------------------------------- |
---|
| 412 | |
---|
| 413 | One new feature, and some improvements and bug-fixes. |
---|
| 414 | |
---|
[264] | 415 | * The user can now elect to Hanning-smooth the spectral direction of |
---|
| 416 | the cube before the searching is done. This can be done instead of the |
---|
| 417 | wavelet reconstruction (although if both are selected, the |
---|
| 418 | reconstruction takes precendence). This functionality is switched on |
---|
| 419 | and off by the flagSmooth parameter, and the width of the Hanning |
---|
| 420 | kernel is governed by the hanningWidth parameter. |
---|
[211] | 421 | |
---|
[264] | 422 | * The smoothed array can also be written to a FITS file and read in at |
---|
| 423 | the start in the same manner as the reconstructed array. |
---|
[211] | 424 | |
---|
[264] | 425 | * Several improvements have been made to the memory |
---|
| 426 | management. Several memory leaks were fixed, and the memory usage of |
---|
| 427 | Duchamp is considerably less, which will be a help for running it on |
---|
| 428 | large cubes. The speed of some routines has also been improved, by |
---|
| 429 | removing unnecessarily repetitive loops. |
---|
[211] | 430 | |
---|
[287] | 431 | * The sorting routines have been changed, and should now be slightly quicker |
---|
[215] | 432 | (and a bit more robust). |
---|
| 433 | |
---|
[287] | 434 | * Enabled the scale bar (plotted on the 0th moment map in the spectral |
---|
| 435 | output file) to be able to deal with much finer spatial scales (as |
---|
| 436 | may be expected from VLBI data). |
---|
[211] | 437 | |
---|
[215] | 438 | * Other minor changes to code and presentation. |
---|
[211] | 439 | |
---|
| 440 | |
---|
[194] | 441 | Duchamp version 1.0.6 (2006/11/01) |
---|
| 442 | ---------------------------------- |
---|
| 443 | |
---|
| 444 | Some incremental improvements and changes to the code: |
---|
| 445 | |
---|
| 446 | * The statistics used for the detection are now calculated once for |
---|
| 447 | the entire cube, rather than being calculated separately for each |
---|
| 448 | channel map and spectrum. This results in a uniform threshold value |
---|
| 449 | across the entire cube. The way the statistics calculations are |
---|
| 450 | implemented in the code was changed, using a new statistics class. |
---|
| 451 | |
---|
| 452 | * This means we can report a peak signal-to-noise value for each |
---|
| 453 | detection. This is given in the text-based output, as well as in the |
---|
| 454 | summary information above the spectral output for each source. |
---|
| 455 | |
---|
[195] | 456 | * The user can now elect not to have the moment map displayed in a |
---|
| 457 | PGPlot window, through the use of the new flagXOutput |
---|
| 458 | parameter. This should help in cases of running Duchamp on a dumb |
---|
| 459 | terminal (with no X-display capabilities), or as part of a batch |
---|
| 460 | script where the continual popping-up of the PGPlot window might |
---|
| 461 | become annoying. |
---|
| 462 | |
---|
[265] | 463 | * There are different Flags reported for the case of an object lying on |
---|
[194] | 464 | the spatial edge and spectral edge of the cube (E and S respectively). |
---|
| 465 | |
---|
| 466 | * A few minor bug fixes: |
---|
| 467 | - When providing a subsection, any axes with the range given by * |
---|
| 468 | had the pixel values offset incorrectly. This has been fixed. |
---|
| 469 | - Negative searches now correctly invert the reconstructed array |
---|
| 470 | when it is read in from a file. |
---|
| 471 | |
---|
| 472 | * Other minor changes to the code, to improve memory handling and |
---|
| 473 | remove redundant arrays and variables. |
---|
| 474 | |
---|
| 475 | |
---|
[146] | 476 | Duchamp version 1.0.5 (2006/09/06) |
---|
[142] | 477 | ---------------------------------- |
---|
| 478 | |
---|
[146] | 479 | Since version 1.0, there have been a range of mostly minor |
---|
| 480 | changes. These are generally bug-fixes, thanks to feedback from a |
---|
[157] | 481 | number of users, as well as a few extra features. |
---|
[146] | 482 | |
---|
[157] | 483 | Here is a summary of the changes since version 1.0: |
---|
[142] | 484 | |
---|
[157] | 485 | USER INTERFACE: |
---|
[146] | 486 | |
---|
[142] | 487 | * The user is now able to enter the value of flag parameters in the |
---|
| 488 | parameter file as strings, as well as integers (so one can enter |
---|
| 489 | true or 1 to give the same result). |
---|
| 490 | |
---|
[157] | 491 | GRAPHICAL OUTPUT: |
---|
| 492 | |
---|
[142] | 493 | * Added a new function that draws the edge of the blank region on the |
---|
[146] | 494 | maps and image cutouts. There is a new flag parameter -- |
---|
| 495 | drawBlankEdges -- that controls the use of this function (the |
---|
| 496 | default is true). |
---|
| 497 | Also, the edge of the field is plotted on the image cutouts as a |
---|
| 498 | yellow line. This distinguishes it from the purple blank pixel |
---|
| 499 | dividing line. |
---|
[142] | 500 | |
---|
[146] | 501 | * The size of the tick mark in the image cutout (indicating the |
---|
| 502 | spatial scale of the image) is now adaptable, so that it shows the |
---|
| 503 | scale that best fits the image (ranging from 1 arcsec to 15 |
---|
| 504 | degrees). (This is a change from the previous standard length of 15 |
---|
| 505 | arcmin.) |
---|
[142] | 506 | |
---|
[113] | 507 | * The zoomed spectrum had its flux scale determined by all points in |
---|
[146] | 508 | the zoom box,including those flagged by the Milky Way (MW) |
---|
| 509 | range. This has been fixed so that any MW channels will not affect |
---|
| 510 | the flux scaling, potentially improving the appearance of objects |
---|
| 511 | detected close in velocity to strong MW emission or absorption. |
---|
[112] | 512 | |
---|
[157] | 513 | * Improved Karma Annotation file output, so that files without good |
---|
| 514 | WCS are dealt with appropriately. |
---|
| 515 | |
---|
| 516 | TEXT-BASED OUTPUT: |
---|
| 517 | |
---|
[113] | 518 | * The position resolution in the IAU-format names has been increased, |
---|
| 519 | so that distinct objects close to one another are more readily |
---|
| 520 | distinguished. Overlapping sources at different velocities may |
---|
| 521 | still have the same name, but this is a more general problem. |
---|
[112] | 522 | |
---|
[146] | 523 | * The presentation of the results in tabular format has been improved, |
---|
| 524 | so that the precision of faint sources is sufficient to show their |
---|
| 525 | flux, for instance, is not 0. This also applies to the information |
---|
| 526 | presented above the spectral plots. |
---|
[112] | 527 | |
---|
[146] | 528 | * Improved error and warning reporting, so that the formatting and |
---|
| 529 | presentation of information is clearer. |
---|
| 530 | |
---|
[157] | 531 | FITS I/O RELATED: |
---|
| 532 | |
---|
[172] | 533 | * The previous version was did not deal well with 4 (or greater) |
---|
| 534 | dimensions in the FITS file -- the memory allocation was badly done, |
---|
| 535 | and the axes needed to be in a particular order. This has been |
---|
| 536 | fixed, so that the FITS I/O routines can now handle any sized FITS |
---|
| 537 | file, with axes in any order (as long as there are two spatial and one |
---|
| 538 | spectral axes). |
---|
[157] | 539 | |
---|
| 540 | * When the FITS file does not have a rest frequency defined, there is |
---|
| 541 | no way to transform the frequency axis into velocity. In this case, |
---|
| 542 | the spectral axis is left as frequency, and put into units of MHz. |
---|
| 543 | |
---|
[172] | 544 | * If the beam size is not indicated by the FITS header (through the |
---|
| 545 | BMAJ and BMIN keywords), the user can manually specify the size |
---|
| 546 | through the parameter file (using parameter beamSize). |
---|
| 547 | |
---|
| 548 | * The FITS I/O routines make use of the command fits_file_exists. This |
---|
| 549 | was introduced into the CFITSIO package from version 2.5, so Duchamp |
---|
| 550 | will only work with this version or higher. Version 3+ is preferred. |
---|
| 551 | |
---|
[157] | 552 | OTHER BUG-FIXING: |
---|
| 553 | |
---|
| 554 | * Several bugs related to the blank-pixels were dealt with. |
---|
| 555 | The correct values were not propagating through correctly, meaning |
---|
| 556 | that the default values were being used. Also, there was the |
---|
| 557 | potential for blank pixels to be detected when the FDR method was |
---|
| 558 | used. Finally, the image cutout displays occasionally did not |
---|
| 559 | obscure the blank pixels. All these problems should be fixed now. |
---|
| 560 | |
---|
| 561 | * The FDR setup routine was occasionally failing with a segmentation |
---|
| 562 | fault, due to it accessing unallocated memory. This has been fixed. |
---|
| 563 | |
---|
| 564 | * Fixed bugs that affected the way 2D images (rather than 3D cubes) |
---|
| 565 | were dealt with. The minChannels test failed when minChannels=0. |
---|
| 566 | |
---|
[165] | 567 | * If there is no beam information in the FITS header, the user can now |
---|
| 568 | enter a beam size as a parameter. |
---|
| 569 | |
---|
[146] | 570 | * Other minor fixes to the code, to improve its structure and |
---|
| 571 | robustness. |
---|
| 572 | |
---|
| 573 | |
---|
[112] | 574 | Duchamp version 1.0 (2006/06/30) |
---|
| 575 | -------------------------------- |
---|
| 576 | |
---|
[113] | 577 | The first full release version. |
---|
| 578 | |
---|
| 579 | Several major changes compared to previous version: |
---|
| 580 | |
---|
| 581 | * Incorporated a FitsHeader class, to keep track of header information |
---|
| 582 | (including WCS) from the FITS file. All units-related calculations |
---|
| 583 | are done with this class, such as conversion to velocity, and |
---|
| 584 | pixel-to-WCS conversion. The main advantage of this is improved |
---|
| 585 | speed of execution, more adaptable to different units, and improved |
---|
| 586 | code (easier for me!). |
---|
| 587 | |
---|
| 588 | * Able to select the dimensionality of the reconstruction, rather than |
---|
| 589 | just doing it in 3-D. |
---|
| 590 | |
---|
| 591 | * Robust text output, so that the columns are sized correctly given |
---|
| 592 | the data being written, and the relevant units (derived from the |
---|
| 593 | FITS header) are included in the output table. |
---|
| 594 | |
---|
| 595 | * Improved the way the reconstructed image is saved and read back in, |
---|
| 596 | with a better naming scheme and a better set of new FITS header |
---|
| 597 | keywords. The MW headers are now not written. |
---|
| 598 | |
---|
| 599 | Other, more minor changes include: |
---|
| 600 | |
---|
| 601 | * MW channels now not set to 0, but just ignored for searching |
---|
| 602 | purposes (and for the scaling of the spectral output -- they are |
---|
| 603 | plotted though and highlighted with a hashed box). |
---|
| 604 | |
---|
| 605 | * The scale bar in the moment map in the spectral output file has its |
---|
| 606 | length indicated. |
---|
| 607 | |
---|
| 608 | * Improved VOTable output, with more information about the way Duchamp |
---|
| 609 | was run. |
---|
| 610 | |
---|
| 611 | * Made sure all statistics calculations in the reconstruction are |
---|
| 612 | based on robust (ie. median) methods. |
---|
| 613 | |
---|
| 614 | Also moved to an autoconf based configure/make compilation, for ease of use. |
---|
| 615 | |
---|
| 616 | |
---|
| 617 | Duchamp version 0.9.2 (2006/06/27) |
---|
| 618 | -------------------------------- |
---|
| 619 | |
---|
| 620 | * Added the use of warning flags for detections: for either edge |
---|
| 621 | location or negative enclosed flux. |
---|
| 622 | |
---|
| 623 | * A new command line option added so that the user can specify a FITS |
---|
| 624 | file to be searched with the default parameters, rather than giving |
---|
| 625 | a full parameter file. |
---|
| 626 | |
---|
| 627 | * The scale bar on the moment cutout now works in general case, using |
---|
| 628 | an iterative process to get to the correct length. |
---|
| 629 | |
---|
| 630 | * RA axes are now done correctly on the full-scale maps. |
---|
| 631 | |
---|
| 632 | * Various coding fixes, improving the readability of the code, |
---|
| 633 | particularly regarding the declaration of temporary variables. |
---|
| 634 | |
---|
| 635 | * Changed the way the time of execution is called : now using time.h |
---|
| 636 | functions, rather than a system call. |
---|
| 637 | |
---|
| 638 | * Changed output in Detection/outputDetection.cc for the pixel |
---|
| 639 | information of detected objects (the info printed above the spectral |
---|
| 640 | plots) -- now print out number of voxels in object, and have |
---|
| 641 | improved the formatting a little. |
---|
| 642 | |
---|
| 643 | * Improved the functionality of saving the FITS file containing the |
---|
| 644 | reconstructed array. Added ability to save and read back in the MW |
---|
| 645 | parameters when saving a reconstructed cube (if a reconstructed cube |
---|
| 646 | has different MW range to that being considered, it would cause |
---|
| 647 | problems). ReadRecon now just returns a FAILURE and will reconstruct |
---|
| 648 | the cube anew. |
---|
| 649 | |
---|
| 650 | * Changed default value of flagGrowth to 0. |
---|
| 651 | |
---|
| 652 | |
---|
| 653 | |
---|
| 654 | Duchamp version 0.9 (2006/05/05) |
---|
| 655 | -------------------------------- |
---|
| 656 | |
---|
| 657 | First public release version. |
---|