OpenCV has a function that does this. Looking at the implementation (available on GitHub) we have given an image with dimensions w x h and a camera matrix fov_y = 2arctan(h/2f_y), fov_x = 2arctan(w/2f