How to plot an ortographic projection of the celestial sphere with equatorial coordinates in python, for a given latitude?

Question

I am trying to obtain an ortographic projection of the celestial sphere, with equatorial coordinates, as seen from a certain latitude, as in the following picture:

Grid o