“…x + x, y + y,r(∆x, 0, z ; x, y))• exp(j 4π λ cr (∆x, 0, z ; x, y))dydx+ ∆x+xm+1 xm+1yn+1 yn e r (x + x, y + y,r(∆x, 0, z ; x, y))• exp(j 4π λ cr (∆x, 0, z ; x, y))dydx(19) wherer(x +∆x, y , z ; x +x, y +y) =r(∆x, 0, z ; x, y), and e r (x, y,r(x + ∆x, y , z ; x, y)) ≈ e r (x, y,r(x , y , z ; x, y)) is valid due to the short interval (x m , x m+1 ] and (y n , y n+1 ].Compared with the original PCD imaging, an extra independent spatial correlation represented by the integral over the interval (y n , y n+1 ] in y-direction as shown in (19) is necessary before the recursive imaging in x-direction. To further reduce the computational cost, the integral over the area enclosed byx interval (x m , ∆x + x m ] and y interval (y n , y n+1 ] in (19) defined as S m,n (x + ∆x, y , z ) = ∆x+xm xm yn+1 yn e r (x + x, y + y,r(∆x, 0, z ; x, y)) • exp(j 4π λ cr (∆x, 0, z ; x, y))dydx (20) can be also calculated recursively along y-direction.…”