SPLIT участков устойчивого направления 'вниз-вверх'

[nickleb]

Здравствуйте. С Наступающим. Всего хорошего и желаемого.
Каким быстрым алгоритмом определить участки устойчивого направления "вниз"-"вверх"
у такого вида наблюдений как зондирование толщи вод?
Зонд даёт показания глубины (горизонта наблюдений) своего нахождения и др. показатели воды с определенной дискретностью по времени.
Пять этапов зондирования:
1)Зонд вначале на подповерхностном горизонте некоторое время - выдерживаются датчики (глубина то чуть выше, то чуть ниже подповерхностного горизонта).
2)Затем его опускаем до придонного горизонта.
3)На придонном горизонте тоже некоторое время зонд выдерживается (глубина то чуть выше, то чуть ниже придонного горизонта).
4)Зонд поднимаем.
5)На подповерхностом горизонте зонд также некоторое время даёт показания, пока его датчики не выключат и не вынут из воды (глубина то чуть выше, то чуть ниже подповерхностного горизонта).
Надо выделить 2)-ой и 4)-ый этапы.

[gamm]

пока ничего не понятно. У вас датчик работает как высотомер, дает расстояние до дна? с каким дискретом пишутся данные, сколько отметок глубины приходится на каждый режим (дно/поверхность)? насколько велика вариабельность глубин в обоих режимах (в первую очередь по сравнению с размахом между режимами)?

самый простой вариант, если имеем не менее (например) 20 отсчетов в каждом режиме, то гоним локальное окно размером 40, вычитая из медианы левых 20 отсчетов медиану правых 20 отсчетов. Получите отчетливые пики вверх и вниз при переключении режимов. Сглаживаем полученную картину Савицким-Галаем, и берем первую и вторую производную, на которых все пики как на ладони (смотрим что-нибудь про поиск пиков в спектрографии и хроматографии).

Если копать глубже, то нужно хоть что-то знать о распределении всех величин, в том числе и о длительности каждого режима. Тогда получается типичная state-space model, копать в сторону HMM (hidden Markov model)

[nickleb]

пока ничего не понятно. У вас датчик работает как высотомер, дает расстояние до дна? с каким дискретом пишутся данные, сколько отметок глубины приходится на каждый режим (дно/поверхность)? насколько велика вариабельность глубин в обоих режимах (в первую очередь по сравнению с размахом между режимами)?

самый простой вариант, если имеем не менее (например) 20 отсчетов в каждом режиме, то гоним локальное окно размером 40, вычитая из медианы левых 20 отсчетов медиану правых 20 отсчетов. Получите отчетливые пики вверх и вниз при переключении режимов. Сглаживаем полученную картину Савицким-Галаем, и берем первую и вторую производную, на которых все пики как на ладони (смотрим что-нибудь про поиск пиков в спектрографии и хроматографии).

Если копать глубже, то нужно хоть что-то знать о распределении всех величин, в том числе и о длительности каждого режима. Тогда получается типичная state-space model, копать в сторону HMM (hidden Markov model)

gamm, здравствуйте. спасибо за подсказку направления куда двигаться.

-датчик даёт расстояние от поверхности до своего положения
-опрос всех датчиков идёт с частотой 16 Hz
-скорость "майны-виры" неодинакова и неравномерна (т.е. зонд тоже слегка может "тюлюпаться" по глубине)
-изменчивость глубины 1)-го и 5)-го этапов низка, а 3)-го чуть выше

[gamm]

-опрос всех датчиков идёт с частотой 16 Hz

это "ни о чем", нужны временнЫе характеристики (сколько времени идет режим, сколько подъем, сколько спуск - типичные времена и разброс, не "на глаз" же вы прибор опускаете/поднимаете, небось лебедка с командоаппаратом есть)

[nickleb]

у меня в распоряжении только набор файлов с набором отсортированных по времени записей показаний прибора... я лично это не наблюдал... связаться с наблюдателем ныне возможности нет..

[nickleb]

http://www.oceanscience.com/products/un ... y-ctd.aspx

[gamm]

у меня в распоряжении только набор файлов с набором отсортированных по времени записей показаний прибора... я лично это не наблюдал... связаться с наблюдателем ныне возможности нет..

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

кстати, раз это стандартные данные, то и программа обработки там же должна быть, на сайте, и все характеристики. Ищите там ссылки

[nickleb]

ладно, gamm... я малёхо ещё сам покумекаю... Вас обременять пока не хочется... если, что: вышлю, напишу... спасибо за уделённое мне время... R-Cran библиотеки, м.б, порекомендуете какие по Вашему первому сообщению?

[gamm]

R-Cran библиотеки, м.б, порекомендуете какие по Вашему первому сообщению?

медианы и сравнение - руками (скриптик написать), производные - savitzkyGolay {prospectr}, HMM лучше поискать по ключевым словам, в том числе на r-bloggers (http://www.r-bloggers.com/)

P.S. если готовых программ для обработки не найдете, и конечная задача содержательная (предполагается статья в приличном издании), то можно дальше пообсуждать