Нет, пока в SA ничего там не разглядел.
Имеешь в виду на водопаде? Судя по тому, что я наблюдал в HDSDR - на водопаде можно ничего и не увидеть - во всяком случае в HDSDR я не разглядел никакого намека на сигнал. Надо видимо свертку "вслепую" делать - т.е. вырезать весь 10-минутный кусок на этой частоте и дальше демодулировать и строить свертку

Надо видимо свертку "вслепую" делать - т.е. вырезать весь 10-минутный кусок на этой частоте и дальше демодулировать и строить свертку
Так и делал, плюс ещё производил поиск ACF в выделенной полосе 30-40 герц вокруг 319,6 kHz - ничего не нашёл...

Хм.. неужто спец-автомат тут лучше? Интервал не пробовал ставить известный? (на скриншоте он есть - 4,850 секунды)
Могу вечером загрузить построенную спектру. Ну или можешь попробовать построить - настройки "стандартные" - fft размером примерно в 4-5 герц, оверлэп 90%

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

Впрочем, на текущий момент считаю, что чувствительности (которая уже точно лучше уха) уже хватает для выпуска в "продакшен" первой версии, равно как и для того, чтобы развивать прогу в сторону функционала автоматизации "скимминга".

Могу вечером загрузить построенную спектру
KarapuZ
Вот: https://yadi.sk/d/5H06plhq0T6EYQ

wazzoo
Вот: https://yadi.sk/d/5H06plhq0T6EYQ
Качаю, и как-раз дошли руки до тестов софта - обрабатываю свой IQ файл:

32flt бита теперь пошли, спасибо, а то у меня все архивы в этом формате.

KarapuZ
Ага - узнаю на скрине FLT32, фирменный Максусовский с высокой полкой и крутыми боками :)
Тебе, видимо, нужна еще и функция склейки? Чтобы из "коротких" 7-минутных отрезков клеить 15-20 минутные?

Тебе, видимо, нужна еще и функция склейки? Чтобы из "коротких" 7-минутных отрезков клеить 15-20 минутные?
Не помешает:), 7 минут в FLT32 при полосе 625kHz "весят" 1,95 Гб - размеры файла ограничиваю в настройках
HDSDR, ибо часто их визуализирую в SDR Data файл анализаторе, а у него лимит на размер входных файлов,
не более 2 Гб, зато читает массивы из файлов, при условии что они записаны подряд без временнЫх пауз.

И было-бы удобно, что-бы в имени файлов спектры была указана ширина спектра и центральная частота в привычных герцах. А так, все работает, и эта рисовалка под сонограммой очень удобна!
Маяки ищутся на "ура" одним кликом по сонограмме:

Не помешает
Ага, сделаю

И было-бы удобно, что-бы в имени файлов спектры была указана ширина спектра и центральная частота в привычных герцах.
Я сегодня добавил прописывание в имя файла диапазона - что-то вроде "300.100-350.200" - на работе осталась версия, завтра могу отправить

Маяки ищутся на "ура" одним кликом по сонограмме:
Это радует :)

Ошибка.

KarapuZ
Новая версия вьювера. Изменения только в части функционала (алгоритмы не менялись)

https://yadi.sk/d/FaCQYj09AFc3-g



Добавлены функции на главное окно:
1. Вывод значения частоты на шкале
2. Вывод меток времени на водопаде
3. "Source lines per line": Сжатие водопада по вертикали (по умолчанию - 8 линий усредняются в одну, можно выбрать другое сжатие). Это позволяет кушать меньше памяти, так же позволяет больше обозревать по времени. На отрисовку свертки не влияет (свертка делается по полному водопаду).
4. "Append": Последовательное добавление файлов в анализ. Т.е. - например, есть отрезки по 7 минут - 4 штуки подряд. Можно загрузить их последовательно. Имена загруженных файлов - в тултипе при наведении на лейбл "x files loaded"
5. Окно статистики по спектру теперь не отдельное - а внизу экрана, рядом с графиком уровня выбранного отсчета.
6. Файлы выбираются проводником, при загрузке рисуется бегунок.



Добавлено на окно анализа NDB:
1. Картинка свертки теперь "интерактивная". Можно двигать интервал, потянув за любую точку на картинке. Т.е. нажимаем ЛКМ и, не отпуская, двигаем влево-вправо - соответственно меняя интервал. Это удобно, если на свертке видны несколько маяков - как правило соседние маяки видны сильно косыми линиями - взявшись внизу за конец такой косой линии можно её "выпрямить".

Дальнейшие планы:

1. В рамках версии "1" (т.е. условно "односигнальной") - добавить логирование принятых маяков. В виде "скриншотов" рабочего окна свертки и записи в табличку.

wazzoo
1. Вывод значения частоты на шкале
2. Вывод меток времени на водопаде

Сразу после загрузки файла спектры во вьювер, шкалы времени и частоты отображаются корректно (правда
цвет шрифта "времени" трудно читаем на фоне спектра):

Но, стоит только воспользоваться нижним ползунком прокрутки спектра, шкала времени "исчезает", равно как и шкала частот:

При возврате ползунка в крайнее левое положение шкала "время" появляется, но шкала частот нет...

По остальным пунктам всё отлично, 3. "Source lines per line": Сжатие водопада по вертикали - этот пункт
тоже актаулен, спасибо!

Но, стоит только воспользоваться нижним ползунком прокрутки спектра, шкала времени "исчезает", равно как и шкала частот
То, что шкала времени исчезает - это нормально - она отрисована прямо поверх водопада, т.е. это все одна картинка - соответственно картинка уехала вбок. А вот то, что исчезла напрочь шкала с частотами - довольно странное явление. У меня она просто "ездит" вслед за ползунком. Буду думать, что может быть не так в твоих условиях :)

Nabludatel
В Москве слышно удовлетворительно БПРМ аэродром Каменск-Уральска, Свердловская обл.
http://www.radioscanner.ru/base/freq36173.html#comm_19164
В REU пишут, что есть маяк 1025 G Саранск. Это гораздо ближе к нам.

В REU пишут, что есть маяк 1025 G Саранск. Это гораздо ближе к нам.
Да, верно, ошибся я с идентификацией буквы. Спасибо за поправку.

Unid KKS 216kHz, 1020 Hz USB/LSB, 8,6 sec. Обкатываю прогу wazzoo. Манипулированные поднесущие KKS на частотах 214,978kHz:

и 217,019kHz:


Уточню параметры по ещё одному сигналу Unid 410 AJA, 10 sec:

Пару лет назад сигнал этот уже принимался. Может у кого есть информация по
этому приводу?

wazzoo
Это потрясающе! В записи от 9 января 2019, в районе частоты 318kHz коррелятор детектировал
последовательность BR, манипуляционные поднесущие на частотах:
317,603 kHz -

и 318,399 kHz -

период последовательности ~8,1 sec, а это всё вместе бъётся, как сигнал из Брунея, 8574 km!

KarapuZ
Ничего себе! Это же супер dx, экзотическая локация. Поздравляю!

Вадим, а можно тебя еще попросить исходный IQ выложить (можно децимированный) - мне для тестов алгоритмов, как пример слабого сигнала?

wazzoo
Спасибо! Но, попробую перепроверить результат по более длительной свёртке. Эта 15-ти минутная запись, а
есть ещё полу-часовая.
Обнаружились ещё интересные неопознанные сигналы:
334 NASB, 10 sec, LSB 400 герц

USB 400 герц


320 RZB, 6,5 sec, LSB 1030 герц


Возможно, туркменский KHM, сменил частоту, отмечен на 333, USB 1044 Hz, 11,8 sec:


339 OB, 5,1 sec, USB 409 Hz:


346 AG, 8,15 sec, LSB 1017 Hz:


348 KRD, 11,8 sec, USB 1040 Hz:

исходный IQ выложить (можно децимированный)
Да, конечно выложу. Полу-часовую попробую загрузить полосу, там где этот брунеец предположительный.

Да, конечно выложу. Полу-часовую попробую загрузить полосу
Главное вот этот кусок, с которого свертка получена, что выше в теме. Я хочу покрутить - вдруг есть возможность более четко показывать такие сигналы.

KarapuZ
wazzoo

Друзья! Нет слов- ваши разработки просто чудеса творят! Поздравления!
Вадим, UNID's крайне интересны (в RWW прием некоторых упоминается только из Египте-RZB,KRD), ну а Бруней -это просто финиш!

Виктор, это не финиш, а старт :)
Жаль, что Рапа-Нуи больше нет в эфире...
----------------------------------------------------
"Живые" NDB из Колумбии (от коллеги Рафаэля):

Увеличить

vinipu
Виктор, ты же вроде тоже HDSDR используешь? Значит прога и тебе подойдет :)

Виктор, это не финиш, а старт :)
Это всё вместе )) ИМХО - это финиш ловли на слух, но старт нового этапа более дальнобойных приёмных систем, равно как и новых DX-ов - не всё же скандинавам с их условиями (низкие атмосферные шумы, хорошая серая зона) дальняки ловить.

Если же еще удастся показать превосходство специально пересчитаной TX3A под ДВ (на бумаге оно - есть, а ММАНА в общем не врет в части диаграмм даже для "приземляных" антенн) - то и преимущество скандинавских Бевереджей будет не очень то сильным - точнее - те, что меньше километра, будут проигрывать. Тут правда нужен будет малошумящий усилок - прямо на антенну (может Хайо попросить помочь его спроектировать)

vinipu
Спасибо! Это всё Иван, он не прошёл мимо идеи и очень изящно реализовал метод в софте, это действительно захватывающе работает, "свежая кровь" в хобби. Полагаю, добью таки эту запись от 9 января на предмет Unid-ов прежде всего. Метод позволяет "расширить" горизонты приёма даже на уже имеющихся приёмных трактах.

wazzoo
Главное вот этот кусок, с которого свертка получена, что выше в теме. Я хочу покрутить - вдруг есть возможность более четко показывать такие сигналы.
Иван, вот ссылка на IQ по Брунею, там децимирована из спектра
полоса шириной 19,5kHz, центральная частота 318 длительность около получаса. На таком длительном интервале, свертка по BR ещё более убедительна и похоже, однозначна! Коррелятор четко, на автомате, определил интервал ~8,1 sec. Сигнал, как оказалось, увеличился в уровне как раз во второй половине записи. При построении файлов спектры я увеличил разрешение FFT до ~2-х герц, оверлап 90, попробуй повторить у себя. Приёмник при записи "маяков" калибрую с точностью до ~1-2 герц, ошибка по частоте маловероятна поэтому и все параметры "бъются" под 318 BR:
317,602 BR (LSB)


318,401 BR (USB)

На таком длительном интервале, свертка по BR ещё более убедительна и похоже, однозначна!
Да - на первом отрезке тоже было абсолютно убедительно. Это сказал бы любой, кто поработал с прогой/свертками по кр. мере. Тем и хороши корреляционные механизмы - вероятность ошибиться даже на предельно слабых сигналах типа этого - на порядки ниже, чем при приеме на слух того что "слышно" (а это будут заметно более мощные сигналы, чем детектируемые сверткой)

Приёмник при записи "маяков" калибрую с точностью до ~1-2 герц, ошибка по частоте маловероятна
Только надо уточнить, что есть специфика отображения частоты для данного метода (анализ по спектрограмме). Дело в том, что частота анализируемого фрейма - это округленная до герца частота начала "фрейма". Т.е., к примеру, имеем центральную частоту IQ равную 500000, ширина записи - 400000. При этом у нас, например, размер окна равен 4 (упростим пример) - тогда частоты фреймов будут равны 300000, 400000, 500000, 600000.
Второй момент - сам сигнал (информация о нём) может присутствовать сразу в нескольких фреймах (это зависит в том числе от оконой функции). Поэтому для точности - нам нужен фрейм, где сигнал максимален. Это будет значить, что сигнал по частоте где-то внутри фрейма. В нашем примере, если сигнал сильнее всего во 2-м фрейме, значит он где-то между 400000 и 500000.
Т.е., переводя на уже практически используемые размеры FFT - когда ты выбираешь фрейм, где сигнал максимален - значит его частота лежит где-то между F и F+0.5+fft_size - где F - отображаемая частота в герцах, а fft_size - размер одного фрейма fft в герцах.

Или еще проще, в случае, если FFT до ~2-х герц, оверлап 90, - то реальная частота если и отличается, то не более чем на 2..3 герца от определенной по фрейму.

wazzoo
Спасибо за уточнения!
Продолжаю собирать "урожай" по DX-ам из январского IQ, ещё один, яркий пример работы проги -
375 VM Исландия, приём на обычный "дворовый" наклонный луч 40m, дистанция 4091 km:

Это какая-то революция в приёме NDB!:)

Это какая-то революция в приёме NDB!:)
И сильно напоминает вот это. А ведь он прав.
Всё это прекрасно, если известны исходные параметры маяка. По сути же, метода - подгон решения задачи под заранее известный результат. А вот что-то новое по такой методе едва ли удастся открыть. 😏

KarapuZ
Очень мощно прошла Исландия - на графике свертки четкий сигнал, равно как на агрегированной сонограмме.

RegRX
Всё это прекрасно, если известны исходные параметры маяка
Нет - никаких "исходных" данных для метода не надо. Единственный принцип, на котором он работает - это то, что сигнал маяка - это самокоррелирующий временной ряд (т.е. то, что в эфир передаются одни и те же манипуляции). Т.е. метод автоматически находит маяки и определяет длительность их периода.

Обратите внимание на последний скриншот Исландии - на верхний график. На нём видны "палки" - это места, где была найдена автокорреляция. Ось Х - длительность периода автокорреляции в миллисекундах, ось Y - степень автокорреляции. Нижняя "полка" на графике - это уровень автокорреляции шумов (шум имеет гауссово распределение, подвыборка на основе равномерного распределения - будет повторять это гауссово распределение - что мы и наблюдаем).

Ну или вот - видео первой пробы пера проги (первая версия авто-алгоритма) - https://www.youtube.com/watch?v=7XKklO9Im4o - там я просто прощелкиваю по фреймам спектрограммы, прога выводит автоматом то, что нашла.

подгон решения задачи под заранее известный результат
О - это важный момент! Дело в том, что метода как раз наоборот даёт более объективный результат, чем прием на слух. Именно прием на слух характеризуется наличием мозгового шумодава, который пытается найти известные ему шаблоны у себя в мозгу, подгоняя услышанное под знакомое. Это известный эффект, приводящий к тому, что люди часто в условиях помех слышат совсем не то, что реально звучит. Т.е. да - мы имеем с методом не только заметно лучшую чувствительность, но и на порядки меньше ошибок идентификации.

Еще важное последствие такой методы - любители приема NDB теперь могут легко обмениваться объективными данными по маякам - в виде такой вот свертки - на которой видны все нюансы манипулятора конкретного маяка - длительность периода, длительности точек и тире, пауз между ними, наличие каких-то "косяков" наподобие сдвига периода. У некоторых маяков даже несущие имеют особености - т.е. автокоррелируют по некоему принципу

- никаких "исходных" данных для метода не надо
О, тогда это действительно "бомба"!

Несколько интересных неопознанных сигналов ещё, приняты 2019-01-09 21:50-22:20Z.
374 JRM, 8,6 sec, LSB 1030 Hz:


410 NAJB, 10 sec, LSB 400 Hz, USB 400 Hz:


414 ELIF, 8,05 sec, LSB 1020 Hz, USB 1020 Hz:


Активирован турецкий 413 GEY.