2013-05-15

Great results from Skyloop

This is just a quick post to show the clear benefits of using skyloop when reducing your own data (see our earlier blog or read: http://www.starlink.ac.uk/docs/sc21.htx/node51.html). Below is an image composed of several observations reduced individually using the iterative map maker. "Blooms" (bright false emission) are clearly evident at the edges of maps where high noise has made its way into the ast model.

 
Reducing the same data using skyloop helps to identify what is noise and what is true emission-particularly in the overlap regions. The benefit in this situation is clearly seen. The cost to the user is in computing power and time but the result is clearly worth it.


2013-05-03

Makemap now gives slightly different numbers

From today, people using a cutting edge starlink rsynced from /stardev at JAC may notice small changes in the maps created by MAKEMAP compared to previous versions.This is caused by a change to the way in which the (RA,Dec) of each bolometer sample is calculated. The change was needed to support FTS-2.

The new scheme is mathematically equivalent to the old scheme, but groups the required transformations a little differently, resulting in different (not worse) numerical rounding errors. This has no effect on the majority of samples, but some samples that are very close to the edge of a pixel will be moved just that little tiny amount needed to push them over into the adjacent pixel. This results in noticeably different (not worse) noise in the final map.

The black curve in the following plot is a histogram of the differences between the two maps produced by the two versions of makemap, for a small Uranus observation (S8A only). The red curve is a histogram of Error component (i.e. sqrt(Variance) ) for the old map. Note the logarithmic Y axis scale. 


Numerically, the RMS difference between the data values in the two maps is 8.57E-5, the square root of the mean Variance value in the old image 53.0E-5, and the  square root of the mean Variance value in the new image is 52.9E-5.

So these changes are much smaller than the noise level, and so should not be a concern.