A new formulation for light propagation in geometric optics by means of the Bi-local Geodesic Operators is considered. We develop BiGONLight Mathematica package, uniquely designed to apply this framework compute optical observables Numerical Relativity. Our package can be used on a wide range scales and redshifts accepts numerical as well analytical input spacetime metric. In paper we focus two...