|
Last change
on this file since 752 was
752,
checked in by phi196, 18 years ago
|
|
Added velocity_align
|
-
Property svn:eol-style set to
native
-
Property svn:executable set to
*
-
Property svn:keywords set to
Author Date Id Revision
|
|
File size:
993 bytes
|
| Line | |
|---|
| 1 | #!/usr/bin/env python |
|---|
| 2 | from asap import * |
|---|
| 3 | |
|---|
| 4 | rcParams['verbose'] = 0 |
|---|
| 5 | |
|---|
| 6 | # Don't plot to the screen... |
|---|
| 7 | del plotter |
|---|
| 8 | plotter = asapplotter(False) |
|---|
| 9 | |
|---|
| 10 | print "Test of Tidbinbilla" |
|---|
| 11 | |
|---|
| 12 | # Create the quotient spectra |
|---|
| 13 | data = scantable('data/tid-t002.rpf') |
|---|
| 14 | s = data.get_scan('*[^we]') |
|---|
| 15 | r = data.get_scan('*[we]') |
|---|
| 16 | q = quotient(s,r) |
|---|
| 17 | |
|---|
| 18 | # Set the restfreq for each IF |
|---|
| 19 | q.set_restfreqs(freqs= [23694.4700e6,23722.6336e6]) |
|---|
| 20 | q.set_unit('km/s') |
|---|
| 21 | q.set_freqframe('LSRK') |
|---|
| 22 | |
|---|
| 23 | # Align frequencies - Tid doppler tracks, so this isn't really necessary |
|---|
| 24 | |
|---|
| 25 | q.freq_align(perif=True) |
|---|
| 26 | q.set_unit('km/s') |
|---|
| 27 | |
|---|
| 28 | # Average in time |
|---|
| 29 | av = average_time(q) |
|---|
| 30 | |
|---|
| 31 | # Do some random processing, just to test these functions |
|---|
| 32 | av.smooth('gauss',5) |
|---|
| 33 | av.scale(1.05) |
|---|
| 34 | av.add(0.05) |
|---|
| 35 | |
|---|
| 36 | # Baseline |
|---|
| 37 | msk=av.create_mask([-70,-50],[40,60]) |
|---|
| 38 | av.poly_baseline(msk,1) |
|---|
| 39 | |
|---|
| 40 | plotter.plot(av) |
|---|
| 41 | plotter.set_mode('i') |
|---|
| 42 | plotter.save('output/tid.png') |
|---|
| 43 | |
|---|
| 44 | # These are currently broken as Tid does not set the elevation correctly |
|---|
| 45 | av.gain_el() |
|---|
| 46 | av.opacity(0.075) |
|---|
| 47 | |
|---|
| 48 | print "Tid test finished successfully" |
|---|
Note: See
TracBrowser
for help on using the repository browser.
