Исследовательский проект "Цифровая оптика" .

[Alex2013]

Решил сделать интегральное продолжение уже начатых мной ранее тем ...




Windows 10 и захват кадра из видеопотока с камеры
Серия проектов "Дополнение к реальности "
Tempus fugit !... (Видео: " ускорение & перестройка" )
Виртуальная и"аргументированная" реальность. Ждем Бума?

Основная цель проекта: создать свой интерфейс на основе технологии "дополненной реальности".
(В чем-то аналогичный интерфейсу ХолаЛенс но в отличии от него полностью программный требующий в идеале одну или в крайнем случае две камеры + есть мысль разобраться с общими возможностями и ограничениями ДР-технологии в целом. )


СМЕЖНЫЕ ТЕМЫ
Ужосы OpenCV...(Работа с библиотекой OpenCV в Лазарусе )
OMarks - библиотека оптических маркеров Как подключить ?

"Полезные темы" (Чтобы не искать по всему форуму )
OpenGl Рисуем в Фоне, возможно ли?
Быстрая отрисовка
Как ускорить прорисовку векторной графики ?
Как оптимизировать определение границ произвольной фигуры ?
(Грандиозный был флейм аж вспомнить приятно.. . но пользы то-же много! )
---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Лабораторный журнал состояние на 12.10.2020
(Журнал обновляется не регулярно (примерно раз в квартал ) так-что если интересно следите за свежими обновлениями в конце темы. )

В процессе исследований мной опробованы и отчасти освоены технологии :
----------------

1 Захвата "живого видео" через ДиректХ

2 "Живого масштабирования" (почти как в шпионских фильмах )

3 Детектора движений (+ получил опыт в написании простой "системы слежения" )\

4 Предпринята попытка добиться Супер Зума
(сборкой кадра с высоким разрешением из нескольких кадров с низким разрешением )

5 Создана экспериментальная библиотека с множеством разных фильтров и эффектов .
(причем в целом достаточно быстро работающая за исключением "сложных" случаев вроде "детектора Канни" )

6 Прикручено "самодельное распознавание меток"...
(расширенная редакция из предыдущей темы Серия проектов "Дополнение к реальности ")

7 "От безысходности" добрался до дико нелюбимого мной OpenCV
(Громоздко, местами явно криво... и вообще на "крестах" то бишь С++ )

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

9 Была идея заменить метки давно "насиженными" 2D-баркодами (QR-code, Aztec, DataMatrx...)...
(В принципе рабочая идея, но найденная мной быстрая библиотека оказалась не очень стабильной, а то что я писал сам тормозило совершенно безбожно )

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

11 Освоены простейшие эффекты в стиле "дополненной реальности"...
(Замена метки на "посторонние" 2D изображение(причем динамическое) с учетом 3D-перспективы )

12 Добавлена простейшая загрузка и анимация 3D-модели(Формат ASC из старого 3DS ) и ее привязка к метке (В "подзорной трубе" пока видно только "каркас" но OpenGL-просмотр в принципе технологически уже доступен )


13 О совсем забыл ! В программе всего кроме прочего есть еще и пока не очень работающее "распознавание ладони" ( то бишь HandsDetect )

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

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

16 А еще есть захват изображений с простых (то есть тех что умеют транслировать по кадрам в обычных jpg ) веб-камер в интернете...

17 В планах добавить скрипты как это было в моей "Лаборатории распознавания образов"
Но разумеется на новом уровне

18 Добавил детектор лиц на каскадах Хаара. (В фас и в профиль )
Не могу сказать что все шикарно работает но интересно то что все действует в комплексе.

19 На конец долгожданный OpenGL ! Нуждается в доработке но уже действует...

20 Добавил экспериментальный механизм планшетного ввода чисел ...

21 Немного вправил мозг "детектору ладони"

22 По случаю попадания в мои "очумелые ручки" шлема (или все-же очков ? ) Lenovo Explorer решил опробовать свою программу и там ... ( Реального стерео-видео-потока с камер шлема пока не получил но есть возможность смотреть на результат работы моей программы в "псевдо-стерео" точнее в
"SBS режиме" )

23. В окно "Цифровые эффекты" на будущее добавил дополнительное "меню выбора библиотеки"
В текущей версии есть Softwrare (27 фильтров) и OpenCV(11 фильтров)

24 Окно "Цифровые эффекты" Software Добавил «Сохранить кадр», «Восстановить кадр», «Кадр в окне» .

25 Окно "Цифровые эффекты" добавил:
"Software" "AntiAlias"
"OpenCV"- Два довольно сложных составных метода cvMorphologiсal и cvMotionDetect
cvMorphologiсal- "морфологический" анализ изображения (по тону насыщенности и яркости )
cvMotionDetect - продвинутый детектор движения отслеживает несколько объектов одновременно.

26 Окно "Цифровые эффекты":Software "Упрощение по цвету"
+ Добавил возможность сохранения и загрузки текущего набора фильтров.
см СтарыйТелевизор.ex1 (входит в комплект)

27 Главное окно: добавил новые "альтернативные режимы" (Загрузка из файла и захват с десктопа и окон) которые становится доступны после включения галки "Сменить источник кадров".(Ранее работал только ввод URL веб-камеры работающей в "JPG-режиме" )

28 Окно "Цифровые эффекты"
"cvStereoTest" попытка работать с "картой глубин".
"cvLkDemo" составной метод "поиск характерных точек" .

29 Окно "Цифровые эффекты"
cvWarpPerspective Искажение и афинные преобразования в 2Д (вписывает кадр в произвольный четырехугольник )

30 Окно "Цифровые эффекты" добавил поддержку библиотеки FastDIB
В текущей версии доступны следующие эффекты и фильтры:
+fdGamma
+fdEmbossEdges
+fdSaturation
+fdFlop
+fdFlip
+fdInvert
+fdSquareWave
+fdRoundWave
+fdMosaic
+fdAutoContrast

31 Окно "Цифровые эффекты" дополнение
(Загрузка и парсинг модели в формате ASC)
Раздел "3D Software" / s3d Загрузка Модели ("Проволочная" модель для превью )
Раздел "3D OpenGL" / ogl Загрузка Модели (OpenGL модель. Пока без текстур )
Оба варианта выводятся на отдельный битмап .

32 Окно "Цифровые эффекты" дополнение
Раздел "OpenCV"
"cvSURF"- Распознавание объектов по методу SURF

33 Окно "слежение за объектами"...
'Отделение фона'- пробная версия.
.
Окно "Цифровые эффекты" дополнение
Раздел "Software"
'Отделение фона' -полная версия .
.
Раздел "OpenCV"
.
'cvPyrMeanShiftFilter' - очень ресурсоемкий метод "сегментации"
(Соседние близкие по цвету пиксели сливаются, а цвет усредняется
- выходит что-то вроде цветной бинаризации )

(Текущая версия (от 12.10.2020) ) Truba_DS_0_0059_14_3_2M0 ТУТ

===========================================================================================

[MysticCoder]

Построение 3д модели пространства в идеале имея на входе одну двигающуюся видеокамеру Можно будет 3д модельки сканировать

[Alex2013]

Убрал из шапки темы ...

И так пока продолжаю разработку внезапно начатую в поезде дальнего следования ....





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

Вообщем "в первом проходе" почти всю функциональность простой системы слежения можно сказать ИССЛЕДОВАЛ !

Далее планирую :
1 Добавить новую редакцию функции "Супер зум "
Что это такое? А вот что !

То есть есть РЯД КАДРОВ низкого разрешения на основе которых я хочу получить один кадр ВЫСОКОГО разрешения .

2 Добавить бар-фильтров
Ну тут думаю все понятно ? Увеличение резкости , размывание, разные цветовые коррекции и т.д. ...
3 Добавить "смесь слежения и распознавания" ... То есть выделяю на кадре участок как образ и слежу социально за ним .

4 Режим "стерео визира " (Изображение с двух камер превращается в стерео картинку) ...

и т.д. постепенно уходя куда-то в строну расширенной и виртуальной реальности ....

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

---------------------------------------------------------------------------------------------------------------------------------------------------------------------

Построение 3д модели пространства в идеале имея на входе одну двигающуюся видеокамеру Можно будет 3д модельки сканировать

Вообще-то возможно ... Microsoft HoloLens вроде умеет ...

https://ru.wikipedia.org/wiki/Microsoft_HoloLens

«голографическая телепортация» (англ. holoportation)[15], позволяющая с помощью специальных камер (как утверждается, минимально необходимое количество камер — 2, использование большего количества камер позволяет улучшить качество) в реальном времени захватить 3D-модель человека и воспроизвести перед носителем HoloLens.

Правда как к этому подступится пока не знаю ... Проще сделать расширенную реальность с маркерами виртуальной клавиатурой и виртуальными экранами.

[Alex2013]

Реальность все расширяется и расширяется...






Начинаю новой этап исследований! Планов как видно по скриншотам громадье ...

Но из новых функций работает только "Супер Зум" в окне которого я победил даже подлое так вредно действующее мне на нервы мерцание !
(Секрет прост нужно рисовать на битмаме а после выводить его на PaintBox в PaintBoxPaint (по OnPaint) )

Truba_DS_0_0058.7z


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

Пока просто беру 2х2 3х3 4х4 и тд кадра и складываю их пиксели по квадратам в один "псевдо пиксель " и строю из них "суммарную картинку" ...
Ясный пень результат полученный таким несложным методом оставляет желать лучшего ... но даже сейчас уже при 64 кадрах видно что в принципе эффект увеличение разрешения получить реально ! Эдакий Anti-Motorn-Blur ( или верно Anti Motion Blur)

Интересно никто про что-то похожее не слышал ? Не может же быть быть что бы я додумался до подобной идеи первым !

[Pavia]

Идея не нова. Но нужно что-бы камера качалась из стороны в сторону.
Мне нравиться идея супер-рес. Там достаточно одного изображения, но никто правда не раскрывает деталей как её сделать. Но я предполагаю сделать так. Нужно восстановить линии(границы) и плавные переходы.
Первое реализуется поиском границ второе функцией размытия. Достаточно обучить НС для подбора коэффициентов и натравить обученную НС на новое изображение.

[Alex2013]

Мне нравиться идея супер-рес...

Вообще-то есть легендарное фрактальное сжатие... (Там можно увеличивать "близко к смыслу" раз в сто и более ) или более банальная трассировка в векторный аналог ( в полоть до построения 3D модели )

Ну пока все читатели темы тихо обтекают от сих "откровений"... я тут уже кое что прикрутил под флажок с многозначительным обозначением "FX "
(На самом деле это пока это довольно примитивный АнтиАлайзинг + Чуть подкручиваю контраст ....)

Кстати, вот код пока немногих переписных для Лазаруса +LCL "Спец эффектов "...

Код: Выделить всё

procedure AntiAliasRect(var Clip: TBitmap; XOrigin,YOrigin,XFinal,YFinal: integer);
var Memo,X,Y: integer; (* Composantes primaires des points environnants *)
    P0,P1,P2: PByteArray;
begin
   if XFinal<XOrigin then
     begin
       Memo:=XOrigin;
       XOrigin:=XFinal;
       XFinal:=Memo;
     end;  (* Inversion des valeurs   *)
   if YFinal<YOrigin then
     begin
       Memo:=YOrigin;
       YOrigin:=YFinal;
       YFinal:=Memo;
     end;  (* si diff‚rence n‚gative*)
   XOrigin:=Max(1,XOrigin);
   YOrigin:=Max(1,YOrigin);
   XFinal:=Min(Clip.Width-2,XFinal);
   YFinal:=Min(Clip.Height-2,YFinal);
   Clip.PixelFormat:=pf24bit;
   for Y:=YOrigin to YFinal do
     begin
     Clip.BeginUpdate(False);
       P0:=Clip.ScanLine[Y-1];
       P1:=Clip.ScanLine[Y];
       P2:=Clip.ScanLine[Y+1];
       for X:=XOrigin to XFinal do
         begin
           P1[X*3]:=(P0[X*3]+P2[X*3]+P1[(X-1)*3]+P1[(x+1)*3]) div 4;
           P1[X*3+1]:=(P0[X*3+1]+P2[X*3+1]+P1[(X-1)*3+1]+P1[(X+1)*3+1]) div 4;
           P1[X*3+2]:=(P0[X*3+2]+P2[X*3+2]+P1[(X-1)*3+2]+P1[(X+1)*3+2]) div 4;
         end;
       Clip.EndUpdate(False);
     end;
end;
function TriangleFilter(Value: single): single;
begin
  if (Value<0.0) then Value:=-Value;
  if (Value<1.0) then Result:=1.0-Value else Result:=0.0;
end;

function Int2Byte(I: integer): byte;
begin
  if I>255 then Result:=255 else
    if I<0 then Result:=0 else Result:=I;
end;

procedure Contrast(var Clip: TBitmap; Amount: integer);
var P0: PByteArray;
    RG,GG,BG,R,G,B,X,Y: integer;
begin
  for Y:=0 to pred(Clip.Height) do
    begin
      Clip.BeginUpdate(False);
      P0:=Clip.ScanLine[Y];
      for X:=0 to pred(Clip.Width) do
        begin
          R:=P0[X*3];
          G:=P0[X*3+1];
          B:=P0[X*3+2];
          RG:=(Abs(127-R)*Amount) div 255;
          GG:=(Abs(127-G)*Amount) div 255;
          BG:=(Abs(127-B)*Amount) div 255;
          if R>127 then R:=R+RG else R:=R-RG;
          if G>127 then G:=G+GG else G:=G-GG;
          if B>127 then B:=B+BG else B:=B-BG;
          P0[X*3]:=Int2Byte (R);
          P0[X*3+1]:=Int2Byte(G);
          P0[X*3+2]:=Int2Byte(B);
        end;
    Clip.EndUpdate(False);
    end;
end;
procedure Lightness(var Clip: TBitmap; Amount: integer);
var P0: PByteArray;
    R,G,B,X,Y: integer;
begin
  for Y:=0 to pred(Clip.Height) do
    begin
      Clip.BeginUpdate(False);
      P0:=Clip.ScanLine[Y];
      for X:=0 to pred(Clip.Width) do
        begin
          R:=P0[X*3];
          G:=P0[X*3+1];
          B:=P0[X*3+2];
          P0[X*3]:=Int2Byte(R+((255-R)*Amount) div 255);
          P0[X*3+1]:=Int2Byte(G+((255-G)*Amount) div 255);
          P0[X*3+2]:=Int2Byte(B+((255-B)*Amount) div 255);
        end;
      Clip.EndUpdate(False);
    end;
end;


А что получилось ? Ну кое что определенно получается ( Видно "просто увеличение" и Супер Зум )


Собственно обновленная версия софтинки ...
Truba_DS_0_0058_1.7z

[Лекс Айрин]

Вообще-то есть легендарное фрактальное сжатие...

Боюсь, что именно легендарное.Как то пробовал фрактальное увеличение...лучше, чем обычное, но до вектора ему еще далеко.

[Alex2013]

Как то пробовал фрактальное увеличение...лучше, чем обычное, но до вектора ему еще далеко.

Я тоже пробовал погонять ... и понял что реализация заточена под максимальное уменьшение размера файла ...
Если сделать иначе (пусть даже размер будет равен несжатому формату BMP или PNG )то эффект будет гораздо лучше .
Там просто ИДЕЯ мощная ... подобрать фрактал описывающий картинку ... а главное свойство фрактала как раз в возможности масштабирования "до потери точности" ...
(Разумеется простеньким "треугольником Серпинского" там не обойтись но способов построить фрактал существует достаточно много )

[Лекс Айрин]

Если сделать иначе (пусть даже размер будет равен несжатому формату BMP или PNG )то эффект будет гораздо лучше .

А смысл?

[Alex2013]

Смысл как раз в возможности получить увеличение .... а если сделать фрактал в "псевдо 3д континууме" где третье измерении не "глубина" а "изменение кадров во времени" там ТАКОЕ можно вытворять ...
Зы
Кстати стало мне интересно как получили найденный мной "букет фрактальных цветов"и наткнулся я на небольшой "шедевр от программирования "
(Между прочим написан на дельфи ... исходники прилагаются )
https://sourceforge.net/projects/apophysis7x/

А это уже я лично "на фракталил" ... " Всем борцам за дело сингулярности посвящается !"
"Чшерский кот снова улыбается ..."

[Лекс Айрин]

Alex2013, смысл фрактальной графики и есть в том чтобы используя минимум формул нарисовать картину. А если размер будет такой, что проще его сохранить в bmp, то смысл исчезает. Идеал это 1% и меньше от bmp
Пример идеала это любой организм, который из одной клетки вырастает в полноценную особь.

[Alex2013]

Есть такой прием "избыточное кодирование" обычно причины его применения в предотвращении потерь ...
Но бывают случаи, когда например с камеры удобнее читать данные в формате RAW то есть вообще не только без сжатия но и с разным техническим мусором ...
Просто потому, что потом данные можно обработать не "экспресс алгоритмами " на слабеньком контролере камеры, а "в полный профиль" медленными но качественными методами .

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

Да, разумеется, сделать фрактально-темпоральное кодирование для меня пока сложновато будет но просто "темпоральное" можно попробовать. Как ? Да просто! Найти движение и отделить его от фона ... Что вроде автоматической генерации спрайтов . Разумеется нужно отслеживать изменения самих движущихся объектов их пресечения и каким-то образом научиться дорисовывать фон заслоняемый пришедшими в движение но ранее неподвижными объектами ...

Но способ кодирования видео фрагментов получится очень интересный и несущий в себе массу ДОПОЛНИТЕЛЬНОЙ информации и удобств для систем "технического зрения" и программ с технологиями дополненной реальности... (то есть фактически при просмотре "TVF-файлов" будет производится "2Д рендаринг" изображений .)
Вообще-то MP4 делает что-то похожее но до "спрайтов" там добраться довольно сложно + в основе лежит "квадратно-гнездовое " JPG - сжатие ...
Зы
Вот еще одни фрактал чисто "для вдохновения" в духе мксимы "Возможно все !"



(Если не строил этот фрактал сам то не поверил бы что это не коллаж в фотожабе ... Но факт ! )

[Лекс Айрин]

Но бывают случаи, когда например с камеры удобнее читать данные в формате RAW то есть вообще не только без сжатия но и с разным техническим мусором ...

Профессионалы всегда считывают в RAW

Есть такой прием "избыточное кодирование" обычно причины его применения в предотвращении потерь ...

вся полезная информация в ДНК занимает не больше 5-10% (точную цифру не помню)... причем, это вместе с зип архивом))). Все остальное это абсолютно левая инфа, например, осколки вирусов.

[Alex2013]

Но бывают случаи, когда например с камеры удобнее читать данные в формате RAW то есть вообще не только без сжатия но и с разным техническим мусором ...

Профессионалы всегда считывают в RAW

Есть такой прием "избыточное кодирование" обычно причины его применения в предотвращении потерь ...

вся полезная информация в ДНК занимает не больше 5-10% (точную цифру не помню)... причем, это вместе с зип архивом))). Все остальное это абсолютно левая инфа, например, осколки вирусов.

+зашита от ошибок копирования она же "код смерти" (теломераза) .... (На самом деле думаю фокус в том что в организме ВСЕГДА есть эталонные "стволовые клетки" а теломераза просто защита от массового "несанкционированного копирования " не из исходного " репозитория" и не там где нужно.... )
А еще всем известно что ДНК "двойная спираль" что тоже как бы намекает на избыточность...

Зы
Но вернусь к "своему барану" то то бишь к "Подзорной трубе" ...
Пробовал интересный способ получения контуров ...

Код: Выделить всё

function Set255(Clr: integer): integer;
asm
  MOV  EAX,Clr  // store value in EAX register (32-bit register)
  CMP  EAX,254  // compare it to 254
  JG   @SETHI   // if greater than 254 then go set to 255 (Max value)
  CMP  EAX,1    // if less than 255, compare to 1
  JL   @SETLO   // if less than 1 go set to 0 (Min value)
  RET           // otherwise it doesn't change, just exit
@SETHI:         // Set value to 255
  MOV  EAX,255  // Move 255 into the EAX register
  RET           // Exit (result value is the EAX register value)
@SETLO:         // Set value to 0
  MOV  EAX,0    // Move 0 into EAX register
end;            // Result is in EAX
type
  TRGBArray = array[0..32767] of TRGBTriple;
  PRGBArray = ^TRGBArray;

procedure ConvolveE(Ray: array of integer; Z: word; var aBmp: TBitmap);

var
  O,T,C,B: pRGBArray;  // Scanlines
  x,y: integer;
  tBufr: TBitmap; // temp bitmap for 'enlarged' image
begin
  tBufr:=TBitmap.Create;
  tBufr.PixelFormat:=pf24bit;
  tBufr.SetSize(aBmp.Width+2,aBmp.Height+2); // Add a box around the outside...

  tBufr.BeginUpdate(False);
   O:=tBufr.ScanLine[0];   // Copy top corner pixels
   T:=aBmp.ScanLine[0];
   O[0]:=T[0];  // Left
   O[tBufr.Width-1]:=T[aBmp.Width-1];  // Right
  tBufr.EndUpdate(False);
  // Copy top lines
  tBufr.Canvas.CopyRect(Rect(1,0,tBufr.Width-1,1),aBmp.Canvas,Rect(0,0,aBmp.Width,1));
  tBufr.BeginUpdate(False);
   O:=tBufr.ScanLine[tBufr.Height-1]; // Copy bottom corner pixels
   T:=aBmp.ScanLine[aBmp.Height-1];
   O[0]:=T[0];
   O[tBufr.Width-1]:=T[aBmp.Width-1];
  tBufr.EndUpdate(False);
  // Copy bottoms
  tBufr.Canvas.CopyRect(RECT(1,tBufr.Height-1,tBufr.Width - 1,tBufr.Height),
         aBmp.Canvas,RECT(0,aBmp.Height-1,aBmp.Width,aBmp.Height));
  // Copy rights
  tBufr.Canvas.CopyRect(RECT(tBufr.Width-1,1,tBufr.Width,tBufr.Height-1),
         aBmp.Canvas,RECT(aBmp.Width-1,0,aBmp.Width,aBmp.Height));
  // Copy lefts
  tBufr.Canvas.CopyRect(RECT(0,1,1,tBufr.Height-1),
         aBmp.Canvas,RECT(0,0,1,aBmp.Height));
  // Now copy main rectangle
  tBufr.Canvas.CopyRect(RECT(1,1,tBufr.Width - 1,tBufr.Height - 1),
    aBmp.Canvas,RECT(0,0,aBmp.Width,aBmp.Height));
  // bmp now enlarged and copied, apply convolve
  for x:=0 to aBmp.Height - 1 do begin  // Walk scanlines
       aBmp.BeginUpdate(False);
    O:=aBmp.ScanLine[x];      // New Target (Original)
    T:=tBufr.ScanLine[x];     //old x-1  (Top)
    C:=tBufr.ScanLine[x+1];   //old x    (Center)
    B:=tBufr.ScanLine[x+2];   //old x+1  (Bottom)
  // Now do the main piece
    for y:=1 to (tBufr.Width - 2) do begin  // Walk pixels
      O[y-1].rgbtRed:=Set255(
          ((T[y-1].rgbtRed*ray[0]) +
          (T[y].rgbtRed*ray[1])+(T[y+1].rgbtRed*ray[2]) +
          (C[y-1].rgbtRed*ray[3]) +
          (C[y].rgbtRed*ray[4])+(C[y+1].rgbtRed*ray[5])+
          (B[y-1].rgbtRed*ray[6]) +
          (B[y].rgbtRed*ray[7])+(B[y+1].rgbtRed*ray[8])) div z
          );
      O[y-1].rgbtBlue:=Set255(
          ((T[y-1].rgbtBlue*ray[0]) +
          (T[y].rgbtBlue*ray[1])+(T[y+1].rgbtBlue*ray[2]) +
          (C[y-1].rgbtBlue*ray[3]) +
          (C[y].rgbtBlue*ray[4])+(C[y+1].rgbtBlue*ray[5])+
          (B[y-1].rgbtBlue*ray[6]) +
          (B[y].rgbtBlue*ray[7])+(B[y+1].rgbtBlue*ray[8])) div z
          );
      O[y-1].rgbtGreen:=Set255(
          ((T[y-1].rgbtGreen*ray[0]) +
          (T[y].rgbtGreen*ray[1])+(T[y+1].rgbtGreen*ray[2]) +
          (C[y-1].rgbtGreen*ray[3]) +
          (C[y].rgbtGreen*ray[4])+(C[y+1].rgbtGreen*ray[5])+
          (B[y-1].rgbtGreen*ray[6]) +
          (B[y].rgbtGreen*ray[7])+(B[y+1].rgbtGreen*ray[8])) div z
          );
    end;
    aBmp.EndUpdate(False);
  end;
  tBufr.Free;

end;
procedure CopyMe(var ToBMP,FromBMP: TBitmap);
begin
  ToBMP.Assign(FromBMP);
end;
{The Ignore (basic) version of a 3 x 3 convolution.

The 3 x 3 convolve uses the eight surrounding pixels as part of the
calculation.  But, for the pixels on the edges, there is nothing to use
for the top row values.  In other words, the leftmost pixel in the 3rd
row, or scanline, has no pixels on its left to use in the calculations.
This version just ignores the outermost edge of the image, and doesn't
alter those pixels at all.  Repeated applications of filters will
eventually cause a pronounced 'border' effect, as those pixels never
change but all others do. However, this version is simpler, and the
logic is easier to follow.  It's the fastest of the three in this
application, and works great if the 'borders' are not an issue. }

procedure ConvolveI(Ray: array of integer; Z: word; var aBmp: TBitmap);
var
  O,T,C,B: pRGBArray;  // Scanlines
  x,y: integer;
  tBufr: TBitmap; // temp bitmap
begin
  tBufr:=TBitmap.Create;
  tBufr.PixelFormat:=pf24bit;
  //tBufr.SetSize(aBmp.Width+2,aBmp.Height+2); // Add a box around the outside...
//  tBufr:=NewBitmap(0,0);
  CopyMe(tBufr,aBmp);
  for x:=1 to aBmp.Height-2 do begin  // Walk scanlines
   aBmp.BeginUpdate(False);
    O:=aBmp.ScanLine[x];      // New Target (Original)
    T:=tBufr.ScanLine[x-1];     //old x-1  (Top)
    C:=tBufr.ScanLine[x];   //old x    (Center)
    B:=tBufr.ScanLine[x+1];   //old x+1  (Bottom)
  // Now do the main piece
    for y:=1 to (tBufr.Width-2) do begin  // Walk pixels
      O[y].rgbtRed:=Set255(
          ((T[y-1].rgbtRed*ray[0])+
          (T[y].rgbtRed*ray[1])+(T[y+1].rgbtRed*ray[2])+
          (C[y-1].rgbtRed*ray[3])+
          (C[y].rgbtRed*ray[4])+(C[y+1].rgbtRed*ray[5])+
          (B[y-1].rgbtRed*ray[6])+
          (B[y].rgbtRed*ray[7])+(B[y+1].rgbtRed*ray[8])) div z);
      O[y].rgbtBlue:=Set255(
          ((T[y-1].rgbtBlue*ray[0])+
          (T[y].rgbtBlue*ray[1])+(T[y+1].rgbtBlue*ray[2])+
          (C[y-1].rgbtBlue*ray[3])+
          (C[y].rgbtBlue*ray[4])+(C[y+1].rgbtBlue*ray[5])+
          (B[y-1].rgbtBlue*ray[6])+
          (B[y].rgbtBlue*ray[7])+(B[y+1].rgbtBlue*ray[8])) div z);
      O[y].rgbtGreen:=Set255(
          ((T[y-1].rgbtGreen*ray[0])+
          (T[y].rgbtGreen*ray[1])+(T[y+1].rgbtGreen*ray[2])+
          (C[y-1].rgbtGreen*ray[3])+
          (C[y].rgbtGreen*ray[4])+(C[y+1].rgbtGreen*ray[5])+
          (B[y-1].rgbtGreen*ray[6])+
          (B[y].rgbtGreen*ray[7])+(B[y+1].rgbtGreen*ray[8])) div z);
    end;
  aBmp.EndUpdate(False);
  end;
  tBufr.Free;
end;
{The mirror version of a 3 x 3 convolution.

The 3 x 3 convolve uses the eight surrounding pixels as part of the
calculation.  But, for the pixels on the edges, there is nothing to use
for the top row values.  In other words, the leftmost pixel in the 3rd
row, or scanline, has no pixels on its left to use in the calculations.
I compensate for this by increasing the size of the bitmap by one pixel
on top, left, bottom, and right.  The mirror version is used in an
application that creates seamless tiles, so I copy the opposite sides to
maintain the seamless integrity.  }
procedure ConvolveM(Ray: array of integer; Z: word; var aBmp: TBitmap);
var
  O,T,C,B: pRGBArray;  // Scanlines
  x,y: integer;
  tBufr: TBitmap; // temp bitmap for 'enlarged' image
begin
//  tBufr:=NewBitmap(aBmp.Width+2,aBmp.Height+2);  // Add a box around the outside...
tBufr:=TBitmap.Create;
tBufr.PixelFormat:=pf24bit;
tBufr.SetSize(aBmp.Width+2,aBmp.Height+2); // Add a box around the outside...
tBufr.BeginUpdate(False);
  O:=tBufr.ScanLine[0];   // Copy top corner pixels
  T:=aBmp.ScanLine[0];
  O[0]:=T[0];  // Left
  O[tBufr.Width-1]:=T[aBmp.Width-1];  // Right
tBufr.EndUpdate(False);
  // Copy bottom line to our top - trying to remain seamless...
  tBufr.Canvas.CopyRect(Rect(1,0,tBufr.Width-1,1),aBmp.Canvas,Rect(0,aBmp.Height-1,aBmp.Width,aBmp.Height-2));
  tBufr.BeginUpdate(False);
  O:=tBufr.ScanLine[tBufr.Height-1]; // Copy bottom corner pixels
  T:=aBmp.ScanLine[aBmp.Height-1];
  O[0]:=T[0];
  O[tBufr.Width-1]:=T[aBmp.Width-1];
tBufr.EndUpdate(False);
  // Copy top line to our bottom
  tBufr.Canvas.CopyRect(Rect(1,tBufr.Height-1,tBufr.Width-1,tBufr.Height),aBmp.Canvas,Rect(0,0,aBmp.Width,1));
  // Copy left to our right
  tBufr.Canvas.CopyRect(RECT(tBufr.Width-1,1,tBufr.Width,tBufr.Height-1),
         aBmp.Canvas,RECT(0,0,1,aBmp.Height));
  // Copy right to our left
  tBufr.Canvas.CopyRect(RECT(0,1,1,tBufr.Height-1),
         aBmp.Canvas,RECT(aBmp.Width - 1,0,aBmp.Width,aBmp.Height));
  // Now copy main rectangle
  tBufr.Canvas.CopyRect(RECT(1,1,tBufr.Width - 1,tBufr.Height - 1),
    aBmp.Canvas,RECT(0,0,aBmp.Width,aBmp.Height));
  // bmp now enlarged and copied, apply convolve
  for x:=0 to aBmp.Height-1 do begin  // Walk scanlines
    aBmp.BeginUpdate(False);
    O:=aBmp.ScanLine[x];      // New Target (Original)
    T:=tBufr.ScanLine[x];     //old x-1  (Top)
    C:=tBufr.ScanLine[x+1];   //old x    (Center)
    B:=tBufr.ScanLine[x+2];   //old x+1  (Bottom)
  // Now do the main piece
    for y:=1 to (tBufr.Width-2) do begin  // Walk pixels
      O[y-1].rgbtRed:=Set255(
          ((T[y-1].rgbtRed*ray[0])+
          (T[y].rgbtRed*ray[1])+(T[y+1].rgbtRed*ray[2])+
          (C[y-1].rgbtRed*ray[3]) +
          (C[y].rgbtRed*ray[4])+(C[y+1].rgbtRed*ray[5])+
          (B[y-1].rgbtRed*ray[6]) +
          (B[y].rgbtRed*ray[7])+(B[y+1].rgbtRed*ray[8])) div z
          );
      O[y-1].rgbtBlue:=Set255(
          ((T[y-1].rgbtBlue*ray[0])+
          (T[y].rgbtBlue*ray[1])+(T[y+1].rgbtBlue*ray[2])+
          (C[y-1].rgbtBlue*ray[3]) +
          (C[y].rgbtBlue*ray[4])+(C[y+1].rgbtBlue*ray[5])+
          (B[y-1].rgbtBlue*ray[6]) +
          (B[y].rgbtBlue*ray[7])+(B[y+1].rgbtBlue*ray[8])) div z
          );
      O[y-1].rgbtGreen:=Set255(
          ((T[y-1].rgbtGreen*ray[0])+
          (T[y].rgbtGreen*ray[1])+(T[y+1].rgbtGreen*ray[2])+
          (C[y-1].rgbtGreen*ray[3]) +
          (C[y].rgbtGreen*ray[4])+(C[y+1].rgbtGreen*ray[5])+
          (B[y-1].rgbtGreen*ray[6]) +
          (B[y].rgbtGreen*ray[7])+(B[y+1].rgbtGreen*ray[8])) div z
          );
    end;
    aBmp.EndUpdate(False);
  end;
  tBufr.Free;
end;

procedure ConvolveFilter(var Src: TBitmap; Filternr,Edgenr: integer);
var
  z: integer;
  ray: array [0..8] of integer;
  OrigBMP: TBitmap;              // Bitmap for temporary use
begin
  z:=1;  // just to avoid compiler warnings!
  case filternr of
    0: begin // Laplace
      ray[0]:=-1; ray[1]:=-1; ray[2]:=-1;
      ray[3]:=-1; ray[4]:=8; ray[5]:=-1;
      ray[6]:=-1; ray[7]:=-1; ray[8]:=-1;
      z:=1;
      end;
    1: begin  // Hipass
      ray[0]:=-1; ray[1]:=-1; ray[2]:=-1;
      ray[3]:=-1; ray[4]:=9; ray[5]:=-1;
      ray[6]:=-1; ray[7]:=-1; ray[8]:=-1;
      z:=1;
      end;
    2: begin  // Find Edges (top down)
      ray[0]:=1; ray[1]:=1; ray[2]:=1;
      ray[3]:=1; ray[4]:=-2; ray[5]:=1;
      ray[6]:=-1; ray[7]:=-1; ray[8]:=-1;
      z:=1;
      end;
    3: begin  // Sharpen
      ray[0]:=-1; ray[1]:=-1; ray[2]:=-1;
      ray[3]:=-1; ray[4]:=16; ray[5]:=-1;
      ray[6]:=-1; ray[7]:=-1; ray[8]:=-1;
      z:=8;
      end;
    4 : begin  // Edge Enhance
      ray[0]:=0; ray[1]:=-1; ray[2]:=0;
      ray[3]:=-1; ray[4]:=5; ray[5]:=-1;
      ray[6]:=0; ray[7]:=-1; ray[8]:=0;
      z:=1;
      end;
    5: begin  // Color Emboss (Sorta)
      ray[0]:=1; ray[1]:=0; ray[2]:=1;
      ray[3]:=0; ray[4]:=0; ray[5]:=0;
      ray[6]:=1; ray[7]:=0; ray[8]:=-2;
      z:=1;
      end;
    6: begin  // Soften
      ray[0]:=2; ray[1]:=2; ray[2]:=2;
      ray[3]:=2; ray[4]:=0; ray[5]:=2;
      ray[6]:=2; ray[7]:=2; ray[8]:=2;
      z:=16;
      end;
    7 : begin  // Blur
      ray[0]:=3; ray[1]:=3; ray[2]:=3;
      ray[3]:=3; ray[4]:=8; ray[5]:=3;
      ray[6]:=3; ray[7]:=3; ray[8]:=3;
      z:=32;
      end;
    8: begin  // Soften less
      ray[0]:=0; ray[1]:=1; ray[2]:=0;
      ray[3]:=1; ray[4]:=2; ray[5]:=1;
      ray[6]:=0; ray[7]:=1; ray[8]:=0;
      z:=6;
      end;
    else exit;
  end;
  OrigBMP:=TBitmap.Create;
  //NewBitmap(0,0);  // Copy image to 24-bit bitmap
  CopyMe(OrigBMP,Src);

  case Edgenr of
    0: ConvolveM(ray,z,OrigBMP);
    1: ConvolveE(ray,z,OrigBMP);
    2: ConvolveI(ray,z,OrigBMP);
  end;
  Src.Assign(OrigBMP);  //  Assign filtered image to Image1
  OrigBMP.Free;
end;



Думал после его применения будет лучше работать датчик движения ...
Увы оказалось это не для моего доморощенного алгоритма . (Или ни какой разницы ли даже хуже... из за более темного изображения и как не странно большей разницы между кадрами. )
Но метод волне рабочий возможно нужно включить сглаживание при его применении но до него я еще не добрался ..

[Лекс Айрин]

Вообще, в организме много защитных механизмов. Начиная от двойного дублирования информации по типу RAID-10 (кстати, избыточной))) ) до уничтожения поврежденных клеток целиком (лейкоцитами). При этом, еще и часть органов продублирована, да и мозг имеет резервы восстановления. В общем, если бы так делали аппаратуру, то она была бы почти вечной. Даже возможность создания дополнительной личности и блокирования вредоносной информации есть.