Return to Deconvolution

I ran the v0.7 reductions with deconvolution on for 50 iterations. I had cut out deconvolution originally because of the funky noise maps, but that was partly an error on my part. There is also an issue with bright sources being largely left over in the noise maps.

The deconvolver does MUCH better at filtering out the fuzzy atmospheric emission, so I want to use it. It leaves some flux from bright point sources behind, though, so I decided to try to make the deconvolution kernel smaller to see if that recovers more of the pointlike flux.

planet fluxes

; ephemerides from the JCMT
; MARS:
; June 30 2005: 730.14 Jy UT:53551
; July 15 2005: 872.83 Jy UT:53566
; Sept 10 2005: 1941.72 Jy UT:53623
; June 5 2006: 553.13 Jy UT:53891
; June 23 2006: 674.14 Jy UT:53544
; Sept 10 2006: 135.79 Jy UT:53896
; July 20 2007: 381.18 Jy UT:54301
; Sept 10 2007: 597.87 Jy UT:54353
;
; URANUS:
; June 30 2005: 43.43 Jy UT:53551
; July 15 2005: 44.35 Jy UT:53566
; Sept 10 2005: 45.78 Jy UT:53623
; June 5 2006: 41.71 Jy UT:53891
; June 23 2006: 42.96 Jy UT:53544
; Sept 10 2006: 41.62 Jy UT:53896
; July 20 2007: 43.90 Jy UT:54301
; Sept 10 2007: 45.57 Jy UT:54353
;
; NEPTUNE:
; June 30 2005: 17.42 Jy UT:53551
; July 15 2005: 17.58 Jy UT:53566
; Sept 10 2005: 17.50 Jy UT:53623
; June 5 2006: 17.04 Jy UT:53891
; June 23 2006: 17.33 Jy UT:53544
; Sept 10 2006: 17.09 Jy UT:53896
; July 20 2007: 17.59 Jy UT:54301
; Sept 10 2007: 17.56 Jy UT:54353

Rant: calc_beam_locations

It took me a few days to figure this out, but "calc_beam_locations" is about 800 lines of wasted space. It does nothing substantive until line 335. Everything to that point is parameter parsing. But there doesn't need to be any of that crap, really, and it should have been outsourced to functions to begin with.

NCDF files are read to get the rotation angle - JUST as an error check! There is no a priori reason to include it.

All that the code does is read in a centroid file (a list of x,y offsets), rotate them, and output them as r,theta,error. Sure, there's a bunch of automated outlier rejection etc, but... seriously?! We don't have enough observations to hold up the statistics necessary for that to begin with! NO ONE would if each observation takes an hour. It's absurd. Odd as it is coming from me, manual rejection makes a lot more sense in this case.

Now, I still have to understand WHY the beam locations are rotated by the fiducial array angle.

distortion mapping done?

Created 'beam_locations_0707.txt' from uranus_070702_o42 with a few contributions from g34.3_070630_o34.

The rest were created by averaging over all of the beam location files

0707 distortion maps

They're consistent but not very close to each other.

sign errors

pinned down the problem. Was a minor sign error in the offsets. Why is it that simple sign errors are ALWAYS the hardest things to track down?

Now, open question: should x,y scaling be free parameters or not? What I mean is, when I measure the positions of bolometers on the array using the planet map, should I allow the X and Y stretch (the bolometer spacing) to change? Should it be a uniform stretch in X and Y or should it be allowed to 'distort' too? My opinion is, none of the above: I'm measuring their actual positions (in terms of a fixed spacing) and therefore stretching or distorting to match the nominal positions is not necessary.

BL/noBL Mars

Mars using uranus_070702_o42 distortion map:

BL070706_o38_raw_ds1.nc_indiv0    134.8757324     -0.2325330      0.0721872   3282.5439453      2.0985703      2.1959550   3147.6682129      2.3311033      2.1237679
BL070706_o39_raw_ds1.nc_indiv0    188.0588379     -0.0947621     -0.0643225   3782.3984375      2.1416662      1.9909470   3594.3395996      2.2364283      2.0552695
BL070706_o40_raw_ds1.nc_indiv0    451.6955566     -0.0427303     -0.3303463   3461.9191895      2.1925442      2.0726480   3010.2236328      2.2352746      2.4029944
BL070713_o36_raw_ds1.nc_indiv0    294.5898438     -0.2068889      0.0493107   5975.9482422      1.9935303      2.1527038   5681.3583984      2.2004192      2.1033931
BL070713_o37_raw_ds1.nc_indiv0    813.5678711     -0.4597032      0.0361557   5515.4501953      2.0334179      2.1772604   4701.8823242      2.4931211      2.1411047
BL070713_o38_raw_ds1.nc_indiv0    286.2934570     -0.1904230      0.0520797   5967.7011719      1.9853363      2.1624830   5681.4077148      2.1757593      2.1104033
BL070713_o39_raw_ds1.nc_indiv0    773.2104492     -0.0816238     -0.3387871   5409.8344727      2.0317609      2.1736624   4636.6240234      2.1133847      2.5124495
BL070714_o76_raw_ds1.nc_indiv0    244.2031250     -0.2079357      0.0572422   5312.0419922      1.9982606      2.1597764   5067.8388672      2.2061963      2.1025343
BL070714_o77_raw_ds1.nc_indiv0    688.7365723     -0.1176772     -0.3138330   4706.1503906      2.0379941      2.1989338   4017.4138184      2.1556714      2.5127668
BL070715_o45_raw_ds1.nc_indiv0    108.5583496     -0.0652246     -0.0818899   2397.9819336      2.1608708      1.9907858   2289.4235840      2.2260954      2.0726757
BL070715_o46_raw_ds1.nc_indiv0    268.4321289     -0.2927990     -0.0692098   2236.4533691      2.1387808      2.0159581   1968.0212402      2.4315798      2.0851679
BL070717_o47_raw_ds1.nc_indiv0    372.6528320     -0.2188239      0.0690587   7070.0649414      1.9820960      2.1670835   6697.4121094      2.2009199      2.0980248
BL070717_o48_raw_ds1.nc_indiv0    498.4628906     -0.2993519      0.0410304   6636.1396484      2.0418782      2.1088474   6137.6767578      2.3412302      2.0678170
BL070718_o41_raw_ds1.nc_indiv0    287.1591797     -0.2016468      0.0641356   6684.1699219      1.9891322      2.1470959   6397.0107422      2.1907790      2.0829604
BL070718_o42_raw_ds1.nc_indiv0    841.6083984     -0.0809386     -0.3511798   5934.2832031      2.0332539      2.1427581   5092.6748047      2.1141925      2.4939380
BL070725_o34_raw_ds1.nc_indiv0    245.6381836     -0.0810674     -0.0582435   5918.6728516      1.9807132      2.1547492   5673.0346680      2.0617807      2.2129927
BL070725_o35_raw_ds1.nc_indiv0    697.7910156      0.1512802     -0.6174591   4873.1469727      2.3222449      2.0499101   4175.3559570      2.1709647      2.6673691
BL070726_o33_raw_ds1.nc_indiv0    216.2124023     -0.0697458     -0.0627950   5056.3364258      2.1531394      1.9873902   4840.1240234      2.2228851      2.0501852
BL070726_o34_raw_ds1.nc_indiv0    433.0415039     -0.0154374     -0.3017542   4586.9204102      2.1483676      2.0385454   4153.8789062      2.1638050      2.3402996
BL070730_o23_raw_ds1.nc_indiv0    105.4489746      0.1025155     -0.2241490   2679.8452148      2.1567283      1.9786894   2574.3962402      2.0542128      2.2028384
BL070730_o24_raw_ds1.nc_indiv0    440.7482910     -0.3968222     -0.0447321   2595.1044922      2.0618360      2.1782801   2154.3562012      2.4586582      2.2230122