Pengolahan Citra Digital

       r = sqrt((x-xc)^2+(y-yc)^2);

       z = sqrt(rmaks^2-r^2);

       bx = (1 - 1/rho) * asin((x-xc)/...

sqrt((x-xc)^2+z^2));

       by = (1 - 1/rho) * asin((y-yc)/...

sqrt((y-yc)^2+z^2));

       if r <= rmaks

           x2 = x - z * tan(bx);

           y2 = y - z * tan(by);

       else

           x2 = x;