Light Pollutoin Meter (LPM) - устройство для измерения световой засветки

[Palomar]

Может тогда предусмотреть в прошивке отключение дисплея на время измерения.. нажал кнопку, экран погас, как измерения произведены, экран включился и отобразил данные?

Да да, так и реализовано.


Михаил опередил))

[Алексей Юдин]

Вы что, даже Сикорука не читали, что Вас это так удивляет?

Сикорук-то здесь причем?

Сикорук на двух страничках описал все основные реалии наблюдателя слабых объектов.

Удивляет то, что пользователи sqm сообщают о результатах 22.0 и темнее, что как-то выходит за мои представления о разумном. Перед замерами крышку снять забыли, или точность такая?
Или под ясным небом в самом деле можно столько намерить?

Прозрачность не измеряется, поэтому смысла в таком замере нет - он не позволяет сравнить реальную видность дипская, что и является действительным качеством неба.
Краткое название ситуации тут:"Самый лучший SQM в туалете..."

[megavoltt71]

Получается, что прибор нужен в основной массе пару - тройку раз для определения какой-то цифры в любимых местах. Можно обойтись камерой с шириком, но она не даст какое-то числовое значение . При том, что основная масса не совсем понимает что показал прибор, но кто-то говорит что 21 это ого-го! 

[Алексей Юдин]

При том, что основная масса не совсем понимает что показал прибор, но кто-то говорит что 21 это ого-го! 

Это почти как Штрель, только с ним ещё хуже, т.к. основная масса его вообще не понимает, несмотря на элементарнейшее определение и прозрачную связь с качеством изображения.

Но сам по себе SQM как метрика качества неба куда хуже, т.к. он откровенно некорректен.

[Palomar]

но кто-то говорит что 21 это ого-го!

Вот у меня как раз 21.04 - 21.10. После долгого городского заточения восхищает, но спустя 20 минут начинаешь испытывать сожаление и вспоминать про 2000-2005 года, когда МП просматривался до горизонта.

[leviathan]

У меня когда был SQM, я как-то вывел такую систему, что выше 21 уже неплохо вполне, 21.2-21.3 очень хорошо, выше 21.4 - отлично.

21 это же только-только зелёная зона, так-себе. Вот 21.3 уже нормально. Лично для меня отлично это выше 21.6, у нас немало таких мест.

[Palomar]

21 это же только-только зелёная зона, так-себе.

Да, к сожалению это текущая ситуация с моего дачного участка. Теперь за хорошим визуалом потребуется перемещаться, искать новую площадку. Последнее время как раз нахожусь в поиске таких мест. Смотрю карту, еду, а на месте уже прибором оцениваю. Пока в окрестностях удалось поймать лишь более менее приличную зеленую зону: 21.42-21.48. Выезжаю по несколько раз, снимаю показания, строю графики и свою карту, составляю впечатление. В общем жить еще можно.

[Владимир ARS]

У меня когда был SQM, я как-то вывел такую систему, что выше 21 уже неплохо вполне, 21.2-21.3 очень хорошо, выше 21.4 - отлично.

21 это же только-только зелёная зона, так-себе. Вот 21.3 уже нормально. Лично для меня отлично это выше 21.6, у нас немало таких мест.

считаю что выше 21,5 уже погрешность прибора и прозрачности...

[Алексей Юдин]

Погрешность прибора тут не при чём, ошибка именно методическая - подразумевается что одним измерением яркости фона можно охарактеризовать качество неба.

[Дмитрий Смирнов]

Для новичков, которые сами не знают пока как понять , какое лучше поле - это или соседнее, прибор весьма полезен будет.

А еще карта засветки. В первую очередь)

[Serj]

Это абсурд и профанация, если так, причём ещё большая чем само использование SQM на небе лучше 21. Небо 22 SQM вполне себе яркое, делать машинный ноль на этой величине немылимо. Надо минимум до 23-й иметь диапазон с приличной точностью.

Уровень естественной яркости неба – очень спорный вопрос. Даже Чинзано приводит в своих работах разные цифры. В ранних работах упоминается 252 мккд/м2 (21,58 зв.вел./угл.сек.2) с оговоркой что это для такой чистоты атмосферы, при которой горизонтальная видимость составляет 26км. В работе The new world atlas of artificial night sky brightness (2016 Fabio Falchi, Pierantonio Cinzano и др.) уже говорится о 170 мккд/м2 (22 зв.вел./угл.сек.2).

[Грехов Михаил]

с оговоркой что это для такой чистоты атмосферы, при которой горизонтальная видимость составляет 26км.

Это еще не предел. Бывает горизонтальная видимость и в 200...250км запросто.  Но больше зависит еще и от Airglow - естественном свечении неба от солнечного ветра и прочих космических лучей. Вот 2 разных неба - одно Кавказкое (Приэльбрусье), где это дело, видимо разбушевалось и в принципе довольно ярокое. А второе-Чилийское на высоте 5640 метров. Ну и видимо геомагнитная обстановка была абсолютно спокойной. Разница ощутима.

ИМХО для более профессиональной оценки нужен сколько не такой прибор, а некую стандартизированную AllSky камеру (оптика, ПЗС, электроника) которая снимет небо на разных выдержках. И потом уже сравнивать как яркость неба, как проницание для снимка по звездам... так и специфику и характер. Статистику набрать....

[Алексей Юдин]

Стандартизировать сам прибор не выйдет - нет тоталитарных механизмов. Все научные оллскайки разные. Можно только надеяться выдачу причесать до стандартных полос. Впрочем, как мы уже обсуждали выше, визуальщику не особо интересны замеры в полосе V с максимумом на 550 нм. Для рынка любителей нужна специализированная система.

Ну и что сильно мешает широкому распространению таких приборов - в случае с SQM-L пусть и платится почти ни за что, зато совсем немного. Тут будет платиться несколько килобаксов, но за ого-го чего! Психологически на это мало кто пойдёт - пусть и приятно иметь в хозяйстве крутой фишай, чб астрокамерку и колесо фильтров, но посмотрим правде в глаза - тут элементарный 16" добик очень немногие покупают...

[ilch63]

   Купил Ваш прибор, вчера его получил на почте и хочу высказать несколько замечаний.
1. Кнопка выступает и легко нажимается в кармане. Это очевидно снижает срок службы батареи.
2. Индикатор не прикрыт защитным стеклом.
   Защитное стекло толщиной 5-6мм по форме лицевой панели решает эти проблемы. Или можно только кнопку защитить кольцом, его можно или приклеить, или привинтить к панели.
3. Яркость индикатора можно сделать поменьше.
4. Защитный светофильтр болтался и выпадал из гнезда. Оказалось что линза не до конца вставлена в гнездо. Я потряс прибор чтобы светофильтр встал на место и нажал пальцем через лист бумаги на линзу. Линза немного прожалась внутрь и светофильтр перестал болтаться и выпадать из гнезда.
   Я законтрил линзу по краям герметиком ВГО-1, чтобы линза не выпала и не двигалась в оправе.
5. При выключении прибора он мог бы сказать OFF. Иначе меряешь темноту в ванной, и не знаешь что прибор отключился.
6. Тучи в С.Петербурге освещены ярче 15. Прибор показывает Е, как и описано. Насколько критично именно значение 15?
7. Угол зрения прибора на глаз около 20-30 градусов и аналогичен он скорее SQM-L, чем SQM.
8. Направил прибор на темные деревья - показал 17, повторил замер - показывает Е. Не могу понять в чем дело. Обычно показания хорошо повторяются.
9. Прибор легкий и удобно лежит в руке.

 К процессору можно подключить датчики положения и компас, можно подключить GPS и запись на флешку. Кому очень хочется, будет Вам обзор неба по всем направлениям и с геолокацией. А если подключить часы реального времени, то Вы можете даже рисовать суточные  и годовые графики.
   Процессор есть, датчик света есть, детали доступны и недороги. Программатор Вам в руки - и вперед.
   

[Palomar]

   Купил Ваш прибор, вчера его получил на почте и хочу высказать несколько замечаний.
1. Кнопка выступает и легко нажимается в кармане. Это очевидно снижает срок службы батареи.
2. Индикатор не прикрыт защитным стеклом.
   Защитное стекло толщиной 5-6мм по форме лицевой панели решает эти проблемы. Или можно только кнопку защитить кольцом, его можно или приклеить, или привинтить к панели.
3. Яркость индикатора можно сделать поменьше.
4. Защитный светофильтр болтался и выпадал из гнезда. Оказалось что линза не до конца вставлена в гнездо. Я потряс прибор чтобы светофильтр встал на место и нажал пальцем через лист бумаги на линзу. Линза немного прожалась внутрь и светофильтр перестал болтаться и выпадать из гнезда.
   Я законтрил линзу по краям герметиком ВГО-1, чтобы линза не выпала и не двигалась в оправе.
5. При выключении прибора он мог бы сказать OFF. Иначе меряешь темноту в ванной, и не знаешь что прибор отключился.
6. Тучи в С.Петербурге освещены ярче 15. Прибор показывает Е, как и описано. Насколько критично именно значение 15?
7. Угол зрения прибора на глаз около 20-30 градусов и аналогичен он скорее SQM-L, чем SQM.
8. Направил прибор на темные деревья - показал 17, повторил замер - показывает Е. Не могу понять в чем дело. Обычно показания хорошо повторяются.
9. Прибор легкий и удобно лежит в руке.

 К процессору можно подключить датчики положения и компас, можно подключить GPS и запись на флешку. Кому очень хочется, будет Вам обзор неба по всем направлениям и с геолокацией. А если подключить часы реального времени, то Вы можете даже рисовать суточные  и годовые графики.
   Процессор есть, датчик света есть, детали доступны и недороги. Программатор Вам в руки - и вперед.
 

• Да, есть такая неприятность, исправлю. В любом случае потребление там мизерное, если неделями включенным лежать не будет, то за батареи можно не волноваться.
• Хорошо, подумаю и над этим. Светодиодный индикатор защищать, может смысла особого и нет (он по умолчанию пленкой был покрыт, но мне как то не понравилось, снял в итоге), а вот с защитой кнопки от случайного срабатывания придумать что то нужно.
• Если вдруг есть или окажется под рукой программатор, могу выслать прошивку с уменьшенной яркостью.
• Странно, вот за этим я следил особо, чтобы блок линзы и светофильтра стоял надежно. Прототип я и тряс и кидал на пол, все было ОК. Отправляемые изделия тоже проверил - фильтры стояли хорошо.
• Это тоже можно решить перепрошивкой, если у Вас будет возможность. Прошивку могу отредактировать.
• Это нормально. Значение 15 не критично. Просто смысла измерять дальше нет - небо однозначно непригодно.
• Совершенно верно. Прибор аналогичен SQM-L, который также с линзой.
• Деревья мерить, конечно, не совсем корректно. Я на небо направил близ торгового центра. Вроде на глаз все ровно, а показания отличаются, что за чудеса... А причина оказалась в разной яркости рекламных роликов, которые крутились на большой телевизионной панели, установленной на крыше ТЦ. В общем в городе может быть много динамических источников, излучение которых может попасть на деревья, стены домов и прочий сооружений, а также отразиться и на небе, как в моем случае при условии низкой облачности. При измерениях в деревне стабильность результатов была на очень хорошем уровне если направлять прибор на один участок неба.

Про GPS я думал, хотел измерять показания вместе с координатами, сохранять, а потом сливать по USB с помощью десктопной утилитки на компьютер в общую базу, с целью дальнейшей привязки и публикации на картах засветки. Эдакий "глобальный проект". Востребованность такого прибора под сомнением. Идеально он подошел бы ЛА, регулярно, интенсивно исследующим разные локации на предмет засветки. Большинство же, на мой взгляд, имеет несколько любимых площадок и изучает небольшое количество новых. А для этого можно пользоваться и этим самым простым решением, просто запомнив местность. Что касается энтузиастов-хардкорщиков, этим людям полезней будет полноценная установка, описанная Алексеем Юдиным, снимающая во всех нужных областях спектра (по основным линиям засветки) и измеряющая прозрачность. Вот к ней можно и GPS приделать. Я же стремился сделать прибор простым и бюджетным.

ilch63, спасибо за отзыв!

[LILO]

Насчет кнопки согласен.
Думаю подобрать колечко подходящее и приклеить вокруг.
В остальном все ОК.

[AlYev]

В заголовке во втором слове у вас опечатка

[ilch63]

Хорошо, что использованы 3 АА батарейки, вместо канадской кроны. Надеюсь, новый приборчик будет не таким хлипким и разболтанным, как североамериканский

прообраз (который как погремушка, с отваливающимся иногда питанием).

P.S. а температуру он меряет?

А абсолютная темнота - ночью в закрытой ванной без света - интересно чего прибор покажет? 23-24m?

Заданы вопросы, смысл их такой: какую часть измеренного значения яркости неба составляет ток утечки датчика TSL237?
Находим даташит датчика и смотрим стр.2. Видим, что при 50 цельсия нам гарантируют частоту не более 2 Гц. Это параметр fD Dark frequency. Датчик кремниевый, у кремния утечка растет в 2 раза при повышении температуры на 6.5 градуса.

Теперь вопрос: сколько это в m/arcSec^2 magnitudes/arcSecond2 ? Находим в интернете описание аналогичного прибора
https://sourceforge.net/projects/arduinomysqmskyqualitymeter/files/SQM_Daniel-Ola-FreqMeasure/SQM_Daniel-Ola-FreqMeasure.ino/download и находим формулу в тексте программы на языке С++:
Msqm = A - 2.5 * log10(frequency); //Frequency to magnitudes/arcSecond2 formula
где
const float A = 20.7;
   Получаем нужные значения для 50 град., заносим в эксел и строим табличку: в следующей строке уменьшаем температуру на 6,5 и делим частоту на 2.
   В результате получаем два столбца следующей таблицы:
   " Msqm = A - 2.5 * log10(frequency); //Frequency to magnitudes/arcSecond2 formula"
   21.45257499

t____   fD _____  T=1/fD    A __   MsqmDark  %     fD+f21.45  T=1/f  A __   Msqm __    dm
50.0    2.00000   0.5000    20.7    19.9474    400.0   2.5000   0.400  20.7   19.7051    1.747
43.5    1.00000   1.0000    20.7    20.7000    200.0   1.5000   0.666  20.7   20.2597    1.192
37.0    0.50000   2.0000    20.7    21.4525    100.0   1.0000   1.000  20.7   20.7000    0.752
30.5    0.25000   4.0000    20.7    22.2051    50.00   0.7500   1.333  20.7   21.0123    0.440
24.0    0.12500   8.0000    20.7    22.9577    25.00   0.6250   1.600  20.7   21.2102    0.242
17.5    0.06250   16.000    20.7    23.7102    12.50   0.5625   1.777  20.7   21.3246    0.127
11.0    0.03125   32.000    20.7    24.4628    6.250   0.5312   1.882  20.7   21.3867    0.065
4.50    0.01562   64.000    20.7    25.2154    3.125   0.5156   1.939  20.7   21.4191    0.033
-2.0    0.00781   128.00    20.7    25.9680    1.562   0.5078   1.969  20.7   21.4357    0.016
-8.5    0.00390   256.00    20.7    26.7205    0.781   0.5039   1.984  20.7   21.4441    0.008
-15.0   0.00195   512.00    20.7    27.4731    0.390   0.5019   1.992  20.7   21.4483    0.004
-21.5   0.00097   1024.0    20.7    28.2257    0.195   0.5009   1.996  20.7   21.4504    0.002
-28.0   0.00048   2048.0    20.7    28.9783    0.097   0.5004   1.998  20.7   21.4515    0.001
-34.5   0.00024   4096.0    20.7    29.7308    0.048   0.5002   1.999  20.7   21.4520    0.0005

           0.50000    2.0000   20.7     21.4525
          500000  0.000002    20.7     6.45
         1000000  0.000001    20.7     5.7

   Частоты в доли Герца воспринимаются тяжело, следующий столбец - это период импульсов датчика в секундах. Дальше считаем по формуле яркость на квадратную секунду.
Проверяем расчёт тока утечки: при нормальной температуре нам ничего не гарантируют, но обещают частоту около 0,1Гц. У нас в таблице при 24 частота 0,125Гц. Вроде ошибки нет, но и технологический запас по утечке у них небольшой.
Следующий столбец: а сколько это по отношению к темному небу? Пусть темное небо даёт 0,5 Гц с датчика. У нас в Еглино SQM показывает 21,4. Вот и будет 21.45257499 яркостью тёмного неба. При 17.5 цельсия утечка даст 12% яркости неба.

Какова ошибка в звездных величинах? Берем яркость неба 0,5Гц и складываем с утечкой. Полученную частоту переводим в звёздные величины, вычитаем из Msqm 21.45257499 и видим погрешность измерения тёмного неба dm. Теперь внимательно смотрим график 7 на стр. 4 и читаем примечание А внизу. Figure 7, NOTE A: Internal offsets that result in dark frequency can be both positive and negative. The dashed line represents the case of negative offset in which an equivalent amount of light signal is required to obtain a non-zero output frequency.

   То есть утечка может как добавляься к сигналу, так и вычитаться из него. То есть если нам попался датчик с отрицательной утечкой, то вот стали мы его в руках греть - а он нам небо всё лучше и лучше показывает. А если мы его над пуншем горящим немного подержим аккуратно, то он нам вообще абсолютную темноту покажет.
   Теперь вопрос такой: в течение какого времени надо измерять частоту чтобы получить яркость с заданной точностью? каждый импульс датчика TSL237 соответствует накоплению в его интеграторе Nимп электронов: накопились электроны, напряжение на конденсаторе достигло порога, схема переключилась и выдала импульс, накапливается следующая порция электронов. и так по кругу.
   Как по количеству электронов сосчитать статистическую погрешность здесь на форуме многие знают. Погрешность равна корню из числа электронов. Измерили частоту по одному импульсу - получили погрешность sqrt(Nимп ). Измерили частоту по М^2 импульсов датчика - погрешность будет в М раз лучше. Nимп - если интересно - можно оценить так же, как шум матрицы. Нужно измерять дисперсию частоты в зависимости от света, по обычной методике.
   Все почему-то измеряют частоту по 30 импульсам, не меньше.
   Подсветка пульсирует с частотой сети, время измерения  должно быть кратно 1/50 сек. При малых засветках это неактуально - оно и так очень большое, при больших и так ясно что всё плохо, но если задать погрешность измерения большой засветки, то надо про такой момент не забыть. 
   Вернёмся к утечкам. SQM использует термодатчик для компенсации погрешности от тока утечки. По температуре они рассчитывают утечку и корректируют показания прибора. Представьте себе: холодный прибор полежал на чём-то теплом, ну вот в руки его взяли например. TSL237 снаружи торчит, он быстро нагреется. Термодатчик - внутри прибора и нагреется позже. Пока температура TSL237 и термодатчика не сравняются корректировка будет не точной. Поэтому датчик надо ставить рядом с TSL237.
  И вопрос с какой погрешностью мы должны измерять яркость неба? Должны ли мы дополнительную засветку измерять точно? Какой диапазон яркости должен показывать прибор?
Диапазон яркости ограничен:
1. Возможностями датчика.
2. Возможностями программы измерения частоты.
3. Снизу диапазон яркости ограничен допустимым максимальным временем измерения.
4. Практически необходимым диапазоном измерения. Под это требование собственно подбирается датчик и программа обработки его сигнала.

   Датчик мы не выбираем - датчик выбрала фирма Unuhedron, он как бы стал стандартом. О нижней границе его возможностей поговорили ранее, верхняя граница возможностей датчика по даташиту см. Figure 1 на стр.3 и Full-scale frequency‡ 500

1000 kHz на стр. 2 , что по формуле соответствует 6,45-5.7 Msqm.

https://sourceforge.net/projects/arduinomysqmskyqualitymeter/?source=directory
https://sourceforge.net/projects/arduinomysqmskyqualitymeter/files/SQM_Daniel-Ola-FreqMeasure/
https://sourceforge.net/projects/arduinomysqmskyqualitymeter/files/CODE%20ARDUINO/OLED/

https://sourceforge.net/projects/arduinomysqmskyqualitymeter/files/latest/download?source=files
Это файл mySQM[1].pdf
Здесь они пишут следующее:
mySQM readings and Frequency
The following table (use column1 for frequency values) shows approximate mySQM readings. For example, if you measured a frequency around 6,300,000 then this should equate to an SQM reading of around 5. You can use the table to quickly check if your mySQM is working - during the day you might get a reading between 2 and 5 depending upon conditions and location, during the day in a dark room it might be somewhere between 8 and 11. If you have an oscilloscope then you could measure the frequency directly by monitoring the signal output of the TLS237 on data pin 5 of the Nano.
Днем обещают от 2 до 5, а в темной комнате от 8 до 11. Там у них в таблице частоты до 36 МГц. Какое-то сомнительное значение.

[ilch63]

Находим даташит датчика

[ilch63]

Следующий столбец: а сколько это по отношению к темному небу? Пусть темное небо даёт 0,5 Гц с датчика. У нас в Еглино SQM показывает 21,4. Вот и будет 21.45257499 яркостью тёмного неба. При 17.5 цельсия утечка даст 12% яркости неба.

Расчётные данные хотелось бы проверить экспериментально, но вот беда: прибор не показывает значения яркости ниже 22 и измеряет в течение 30с, после чего отключается. Однако, у прибора есть особенность: он "измеряет 30с после последней активности" и если во время измерения яркости нажимать кнопку, то это воспринимается как активность и позволяет увеличить время измерения. Надо всего лишь периодически нажимать кнопку прибора.
Теперь можно производить измерения. Забираемся в тёмную ванную, нажимаем кнопки прибора и секундомера и ждём появления показаний. Не забываем периодически нажимать кнопку LPM. Если при нажатии на кнопку LPM не загораются 00,00 - значит прибор не выключился и измерение продолжается. Через 5мин20сек на табло прибора появляется "Е".
Открываем дверь ванной, нажимаем кнопки, через 4,04сек появляется значение 19,02, подзакрываем дверь - за 11,36сек 20,26.
Вот теперь можно и уточнить значение А для данного прибора, и измерить темновой ток при 25 градусах. Считая, что прибор определяет частоту по 30 периодам, получаем
t               fD            T=30/fD    A               Msqm
25                            4.04                          19.02
25            7.42574   4.04         21.196      19.02010028

25                            11.36                        20.26
25            2.64084   11.36       21.315      20.26064269

25                            320                           Е
25            0.09375   320          21.255      23.82604681

Для измерений с известным временем и яркостью определяем А, усредняем значение А и определяем темновой ток в попугаях.
А=21,26
Яркость тока утечки 23,8 при 25 градусах.

Теперь вопрос: сколько это в m/arcSec^2 magnitudes/arcSecond2 ? Находим в интернете описание аналогичного прибора
https://sourceforge.net/projects/arduinomysqmskyqualitymeter/files/SQM_Daniel-Ola-FreqMeasure/SQM_Daniel-Ola-FreqMeasure.ino/download и находим формулу в тексте программы на языке С++:
Msqm = A - 2.5 * log10(frequency); //Frequency to magnitudes/arcSecond2 formula
где
const float A = 20.7;
   Получаем нужные значения для 50 град., заносим в эксел и строим табличку: в следующей строке уменьшаем температуру на 6,5 и делим частоту на 2.
   В результате получаем два столбца следующей таблицы:
   " Msqm = A - 2.5 * log10(frequency); //Frequency to magnitudes/arcSecond2 formula"
   21.45257499

t____   fD _____  T=1/fD    A __   MsqmDark  %     fD+f21.45  T=1/f  A __   Msqm __    dm
50.0    2.00000   0.5000    20.7    19.9474    400.0   2.5000   0.400  20.7   19.7051    1.747
43.5    1.00000   1.0000    20.7    20.7000    200.0   1.5000   0.666  20.7   20.2597    1.192
37.0    0.50000   2.0000    20.7    21.4525    100.0   1.0000   1.000  20.7   20.7000    0.752
30.5    0.25000   4.0000    20.7    22.2051    50.00   0.7500   1.333  20.7   21.0123    0.440
24.0    0.12500   8.0000    20.7    22.9577    25.00   0.6250   1.600  20.7   21.2102    0.242
17.5    0.06250   16.000    20.7    23.7102    12.50   0.5625   1.777  20.7   21.3246    0.127
11.0    0.03125   32.000    20.7    24.4628    6.250   0.5312   1.882  20.7   21.3867    0.065
4.50    0.01562   64.000    20.7    25.2154    3.125   0.5156   1.939  20.7   21.4191    0.033
-2.0    0.00781   128.00    20.7    25.9680    1.562   0.5078   1.969  20.7   21.4357    0.016
-8.5    0.00390   256.00    20.7    26.7205    0.781   0.5039   1.984  20.7   21.4441    0.008
-15.0   0.00195   512.00    20.7    27.4731    0.390   0.5019   1.992  20.7   21.4483    0.004
-21.5   0.00097   1024.0    20.7    28.2257    0.195   0.5009   1.996  20.7   21.4504    0.002
-28.0   0.00048   2048.0    20.7    28.9783    0.097   0.5004   1.998  20.7   21.4515    0.001
-34.5   0.00024   4096.0    20.7    29.7308    0.048   0.5002   1.999  20.7   21.4520    0.0005