We present a flexible linear optimization model for correcting multi-angle curtaining effects in plasma focused ion beam scanning electron microscopy (PFIB-SEM) images produced by rocking-polishing schemes. When PFIB-SEM is employed in a serial sectioning tomography workow, it is capable of imaging large three-dimensional volumes quickly, providing rich information in the critical 10–100 nm feature length scale. During tomogram acquisition, a “rocking polish” is often used to reduce straight-line “curtaining” gradations in the milled sample surface. While this mitigation scheme is effective for deep curtains, it leaves shallower line artifacts at two discretized angles. Segmentation and other automated processing of the image set requires that these artifacts be corrected for accurate microstructural quantification. Our work details a new Fourier-based linear optimization model for correcting curtaining artifacts by targeting curtains at two discrete angles. We demonstrate its capabilities by processing images from a tomogram from a multiphase, heterogeneous concrete sample. We present methods for selecting the parameters which meet the user’s goals most appropriately. Compared to previous works, we show that our model provides effective multi-angle curtain correction without introducing artifacts into the image, modifying non-curtain structures or causing changes to the contrast of voids. Our algorithm can be easily parallelized to take advantage of multi-core hardware.