To me it seems completely impossible. How exactly do you expect to replace the missing information? Statistically from the neighboring pixels or frames? That works if your goal is to make an image prettier, but if you want to make measurements on the data "inventing" the background is not good enough.
Yes, the sats will occlude things anyway, and their brightness is a big problem, but they are fast moving, they can be mitigated by moving to may short expo images combined. Yes, it's an undue burden, SpaceX should do more, but they can't solve it completely.
In LEO the apparent movement of satellites is around 3.6 arcmin per second, which means for an average sized detector that covers maybe 20 arcmin it takes about 6s to produce a streak that covers the full length of the detector (bit longer depending on the angle). If you want it to cover less than 50% of the length of the detector you have to limit your exposure time to under 3s, which means your duty cycle will be between 60% (if you have a fast camera) and 10%(!).
And that's without even considering the loss of sensitivity you incur, because each image will have lower signal to noise since even with those long readout time the readout noise is still significant.
What am I missing?
Second of all: I highly doubt that this is "trivial", given that observatories are built in the most remote areas of the planet because light pollution was already an enormous problem before these satellites appeared, despite the extra costs that adds to any such project.
A 30 arcsec stripe on an (apparently average) 20 arcmin detector is 2.5% of the total area, so it's not a very large area.
Overall, 20000 satellites at any given time cover 20kπ15acsec² of sky area, which is 0.0005% of total sky area, and over 1 second (assuming a satellite speed of 3.6 arcmin per second) it's 0.02424% of the sky, so roughly 1.45% of the surface area of a 1 minute exposure would have to be thrown out. This assumes that satellite streaks don't overlap, with overlap it would be less.