Almost all of the ingredient needed to add opaque rings already exist: in the sky plane, you can trace out all the relevant regions by traveling along ellipses, so we can reuse all of the machinery in polynomial_limb_darkened_transit. We'd need to solve for intersections between the rings and the planet, then correctly choose which curves to integrate along. We'd also would probably want some new visualization tools and tests.
Almost all of the ingredient needed to add opaque rings already exist: in the sky plane, you can trace out all the relevant regions by traveling along ellipses, so we can reuse all of the machinery in
polynomial_limb_darkened_transit. We'd need to solve for intersections between the rings and the planet, then correctly choose which curves to integrate along. We'd also would probably want some new visualization tools and tests.