Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://jet.sao.ru/hq/sun/manual/WorkScan/WRKSCANR.html
Дата изменения: Unknown
Дата индексирования: Tue Oct 2 06:00:08 2012
Кодировка: Windows-1251
Поисковые слова: п п п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р р п
|
Ссылка на эту страницу
| Где скачать программу
История изменений программы
11 юиля 2007Добавлена
функция коррекции скана с помощью деления фурье образа скана на образ диаграммы.Scan->Correction->FFT Beam Convolution (На основе алгоритма
Коржавина А.Н. "shaber" (доступен
в виде программ на Фортране) )Основные параметры:
- размер участка (Scan width) в солнечных радиусах
(т.е. какой участок скана (от центра) будет
использоваться в FFT.
для примера если он равен 2, значит будет использоваться + - 2 радиуса
Солнца. Если этот параметр отрицательный, то
коррекция положения источника (за счет выноса из фокуса) не будет производиться.
Если равен нулю, то для FFT будет использоваться размер
участка, отображенный в данный момент на экране. Нежелательно использовать
меньше 128 точек, т.к. уменьшение количества точек ведет к ухудшению качества
восстановления изображения.
202 FFTbeamCorrection Arg1 -R pol.dfx; Arg2 -L pol.dfx; Arg3
- sun size(>0 shifted;=0 screen;<0 not shifted); Arg4 - noise level
- Величина выноса из фокуса правой и левой поляризаций в мм
(в используемых в 2005-2007 гг спиральных плоских облучателях расстояние между
осями составляет 35 мм. т.е. вынос каждой поляризации при наблюдениях Солнца
теоретически составляет 17.5 мм. для точечных источников:
http://www.spbf.sao.ru/data/utils/oblnew-schema.jpg). При коррекции следует
учитывать, что реальный вынос из фокуса может отличаться от расчетного. так, для
точечного источника 3с273 от 25 марта 2007 21:31 на волне 6.18 см оптимальное
значение смещения для источника слева составит около -115, а не -154, как
следует по расчетам.
- Уровень FFT шума. Т.е.
уровень FFT диаграммы ниже которого не будет деления скана на нее
(примерно 0.001)
По умолчанию (при первом запуске)
предлагаются параметры такие же, как в программе Коржавина.Если отмечена опция "FFT
Test option" , то при выполнении функции в окне GraphView будут выводиться промежуточные результаты вычислений (форма
диаграммы в линейном и логарифмическом масштабах и фурье образы скана и
диаграммы). Также будет выводится окно сообщения, которое
надо закрыть, чтобы двигаться дальше.В пакетном режиме данная функция тоже работает, но
сначала в диалоговой моде надо настроить параметры. Основные
параметры можно задать.
Добавлена возможность в пакетном режиме в имена
файлов добавлять дату (%y) и время (%t) из текущего скана.
Обновления можно скачать по ссылке
WS-new
26 марта 2007
- добавлена функция меню Correction-> FITS Param.
- При запуске запрашивает файл с расширением *.cor.
- В файле можно задать различные параметры фитс- шапки, в зависимости от выбора:
- 1-полуширина д.н.а,
- 2- температура спокойного солнца,
- 3- коэффициент перехода от Т к потокам в загруженных файлах.
- Формат файла должен быть следующим:
- первая колонка-частота, вторая колонка - меняемый параметр (в промежутке делается линейная
интерполяция).
вводимые из файла параметры могут быть применены к одному скану ии ко всем сканам.
при записи в фитс-файл эти параметры заносятся во вторую таблицу.
в IDL парметры из второй таблицы могут быть прочитаны с помощью процедуры mrdfits:
rat=mrdfits(nfile,0,hdr) & ext1=mrdfits(nfile,1,hdr1); freq=ext1[i].freq
. Запись обратно в фитс файл после исправлений:
mwrfits, rat, nfile, hdr, /create & mwrfits, ext1,
nfile, hdr1
Также был обновлен файл commpaсk.dat, в него включена функция калибровки по ГШ, что
дает возможность использовать эту функцию для пакетной обработки.
Пример использования:
- скрипт
- [Commands]
- ViewWorkI=1;Отображать воркшит I
- ViewWorkV=1;Отображать воркшит V
- SetRangePoint=100 130;Вывести участок от 100-й до 130-ой точки
- SubNullLevel=0; Задать уровень нул равным выделенному участку
- ViewWorkV=0; Не отображать воркшит V
- SetRangePoint=73 120; Вывести участок от 73-й до 120-ой точки (запись ГШ в данном случае)
- CalibrNoiseGenerator=0;Откалиброватьс по ГШ
- ViewWorkV=1
- ViewWorkI=0
- SetRangePoint=60 100
- CalibrNoiseGenerator=0
- ViewWorkI=1
- обновления см в каталоге WS-new
Изменения в функции Scan->Calibration->Qiuet
point Ta
crab0koef.dat
2.80 |
6115.00 |
0.040000 |
3.00 |
6111.34 |
0.042087 |
3.18 |
6105.86 |
0.046247 |
3.37 |
6100.08 |
0.050639 |
3.57 |
6094.00 |
0.055263 |
3.75 |
6088.52 |
0.059423 |
4.00 |
6080.92 |
0.065202 |
4.20 |
6074.83 |
0.069825 |
4.40 |
6068.75 |
0.074449 |
4.50 |
6065.71 |
0.076760 |
4.70 |
6059.62 |
0.081383 |
4.90 |
6053.54 |
0.086007 |
5.10 |
6047.45 |
0.090630 |
5.30 |
6041.37 |
0.095253 |
5.50 |
6035.28 |
0.099876 |
5.70 |
6029.20 |
0.10450 |
6.00 |
6020.07 |
0.11143 |
6.15 |
6015.51 |
0.11490 |
6.30 |
6010.95 |
0.11837 |
6.65 |
6000.50 |
0.12631 |
6.80 |
5996.06 |
0.12969 |
6.95 |
5991.50 |
0.13316 |
7.10 |
5987.00 |
0.13658 |
7.30 |
5981.06 |
0.14109 |
7.60 |
5972.09 |
0.14791 |
7.75 |
5967.55 |
0.15136 |
7.90 |
5962.92 |
0.15489 |
8.10 |
5956.95 |
0.15944 |
8.25 |
5952.51 |
0.16282 |
8.40 |
5947.96 |
0.16629 |
8.70 |
5939.14 |
0.17299 |
8.90 |
5933.04 |
0.17763 |
9.05 |
5928.49 |
0.18107 |
9.20 |
5923.99 |
0.18448 |
9.35 |
5919.57 |
0.18781 |
9.50 |
5914.89 |
0.19136 |
9.70 |
5909.02 |
0.19582 |
10.0 |
5900.00 |
0.20275 |
10.15 |
5895.59 |
0.20617 |
10.30 |
5891.16 |
0.20959 |
10.50 |
5885.06 |
0.21416 |
10.60 |
5881.84 |
0.21648 |
10.80 |
5875.30 |
0.22094 |
10.90 |
5871.82 |
0.22317 |
11.30 |
5855.32 |
0.23236 |
11.40 |
5850.45 |
0.23473 |
11.60 |
5839.95 |
0.23926 |
11.70 |
5834.02 |
0.24153 |
11.90 |
5821.06 |
0.24610 |
12.20 |
5797.55 |
0.25290 |
12.40 |
5778.05 |
0.25756 |
12.50 |
5767.35 |
0.25984 |
12.70 |
5743.29 |
0.26441 |
12.90 |
5716.09 |
0.26883 |
13.05 |
5691.28 |
0.27228 |
13.35 |
5635.77 |
0.27919 |
13.50 |
5606.54 |
0.28249 |
13.70 |
5564.56 |
0.28706 |
14.00 |
5499.85 |
0.29404 |
14.15 |
5469.35 |
0.29740 |
14.30 |
5439.43 |
0.30076 |
14.50 |
5399.20 |
0.30535 |
14.60 |
5379.67 |
0.30760 |
15.00 |
5300.43 |
0.31678 |
15.30 |
5240.01 |
0.32372 |
15.45 |
5209.90 |
0.32715 |
15.60 |
5180.05 |
0.33055 |
15.75 |
5150.73 |
0.33388 |
15.90 |
5120.21 |
0.33735 |
16.10 |
5079.93 |
0.34193 |
16.25 |
5050.18 |
0.34531 |
16.40 |
5020.56 |
0.34869 |
16.50 |
5000.68 |
0.35095 |
16.90 |
4919.69 |
0.36019 |
17.05 |
4889.32 |
0.36366 |
17.20 |
4860.11 |
0.36700 |
17.35 |
4830.65 |
0.37036 |
17.50 |
4800.93 |
0.37376 |
17.70 |
4760.63 |
0.37836 |
17.80 |
4740.09 |
0.38070 |
18.20 |
4657.92 |
0.39008 |
18.30 |
4630.00 |
0.39200 |
|
|
В новой версии изменена функция
калибровки по точке на спокойном Солнце. Теперь Температура
Qsun и коэффициенты из файлов tab_sun.dat
и sh-fits.cfg не используются
в этой функции. Сканы масштабируются на температуру, зашитую в
фитс-файлах во второй шапке, в поле tqsun.
Эти температуры заносятся туда при преобразовании dat
файлов с помощью pas2fits, и
берутся из файла
crab0koef.dat:
Данные
температуры получены с помощью калибровки через ГШ по Луне. Т.е.
полагая, что температура Луны на всех волнах 220 град, калибруем ГШ.
Затем по ГШ калибруем запись Солнца, и снимаем спектр точки на
спокойном Солнце (вторая колонка в файле).
В файле tab_sun.dat третья
колонка, показывающая какая часть Солнца попадает в вертикальную
диаграмму, приравнена к 1,
по рекомендации Коржавина А.Н. и Боровик В.Н., из
соображений, что в вертикальную диаграмму1 попадает
больше Солнца, чем в гауссиану соответствующего размера.
(Согласно рекомендацим семинара от 17.03.2011, калибровать по полному потоку можно только волны длинее 8 см)
В файле
sh-fits.cfg
соответственно изменилась колонка
с tqsun. Но этот файл используется только для
калибровки по полному потоку (с Нобеяма)
Frq -
частота в [GHz] |
Coeff2 -
коэффициент перевода антенных температур в поток
радиоизлучения |
Tqsun -
антенная температура спокойного Солнца на данной частоте в
[К] |
Theta -
полуширина диаграммы направленности |
Tann -
шумовая температура антенны в [K]. |
35.0 |
1.105 |
6000 |
6.885 |
90.316 |
30.0 |
0.804 |
6000 |
7.087 |
71.857 |
25.0 |
0.620 |
6000 |
8.489 |
56.114 |
20.0 |
0.411 |
6000 |
11.894 |
44.899 |
15.0 |
0.267 |
6000 |
18.101 |
40.023 |
14.0 |
0.246 |
6000 |
19.742 |
39.933 |
13.0 |
0.226 |
6000 |
21.491 |
39.978 |
12.0 |
0.203 |
6000 |
23.275 |
40.006 |
11.0 |
0.181 |
6000 |
24.907 |
40.008 |
10.0 |
0.164 |
6000 |
26.193 |
39.991 |
9.00 |
0.147 |
6000 |
27.024 |
39.915 |
8.00 |
0.130 |
6000 |
27.679 |
39.904 |
7.00 |
0.113 |
6000 |
29.036 |
40.315 |
6.00 |
0.096 |
6000 |
33.618 |
41.370 |
5.00 |
0.079 |
6000 |
45.550 |
43.137 |
4.50 |
0.071 |
6000 |
55.584 |
44.307 |
4.00 |
0.064 |
6000 |
69.029 |
45.678 |
3.50 |
0.057 |
6000 |
86.408 |
47.278 |
3.00 |
0.052 |
6000 |
108.14 |
50.451 |
2.50 |
0.048 |
6000 |
134.47 |
58.681 |
2.00 |
0.045 |
6000 |
165.32 |
76.394 |
1.50 |
0.044 |
6000 |
199.68 |
110.054 |
1.00 |
0.043 |
6000 |
| | |