Как сделать алгоритм восстановления фазы Нокса-Томпсона - PullRequest
Новая
       84

Как сделать алгоритм восстановления фазы Нокса-Томпсона

0 голосов
/ 03 июля 2019

Метод Нокса-Томпсона (КТ) 1974 года - одно из первых исследований, которое дало хорошие результаты в восстановлении фаз объектов, наблюдаемых через атмосферную турбулентность.Метод КТ использует общий момент второго порядка или кросс-спектр (двойную корреляцию) вместо спектра мощности.

O (u) O * (u + Δu) = среднее (I (u) I* (u + Δu))

где I - данные изображения (комплексные) в пространстве-частоте Фурье, а * - комплексное сопряжение.

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

Если I - матрица MxN, и есть смещение значений K, будет ли операция такой?? 

Let O[M,N,K,K]
O = 0  

for x, 0 to M
{  
 for y, 0 to N
 {  
  for dx, 1 to K
  {  
   for dy, 1 to K
   {  
    O[x, y, dx, dy] += I[x, y] * cmplx_conj(I[x + dx, x + dy])  
   }  
  }  
 }  
}

O / = (K * K)

или я полностью ошибаюсь?

Тогда как мне это сделать?получить фазу от I, таким образом?:

фаза (u, Δu) = фаза (O (u)) - фаза (O (u + Δu))

Как это происходит вseudocode?

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

...