The SWAP Filter: A Simple Azimuthally Varying Radial Filter for Wide-Field EUV Solar Images

We present the SWAP Filter: an azimuthally varying, radial normalizing filter specifically developed for EUV images of the solar corona, named for the Sun Watcher with Active Pixels and Image Processing (SWAP) instrument on the Project for On-Board Autonomy 2 (PROBA2) spacecraft. We discuss the origins of our technique, its implementation and key user-configurable parameters, and highlight its effects on data via a series of examples. We discuss the filter’s strengths in a data environment in which wide field-of-view observations that specifically target the low signal-to-noise middle corona are newly available and expected to grow in the coming years. Supplementary Information The online version contains supplementary material available at 10.1007/s11207-023-02183-w.


Introduction
A problem of understated difficulty and importance for observational coronal physics is simply finding ways to display and analyze observations so they are D.B. Seaton dseaton@boulder.swri.eduactually useful.The need for such techniques stems from a number of sources, but it generally boils down to a few key facts: the corona has a tremendous dynamic range in brightness, many interesting features and events are inherently faint and hard to detect, and various sources of noise can interfere with both imaging itself and can be amplified by image-processing techniques.Noise sources include instrumental effects and photon shot noise embedded in the signal itself.
This problem is not unique to observations in the extreme ultraviolet (EUV), but it has emerged as an especially important one in this spectral band, particularly in light of the development of new EUV imagers -and new ways of using these imagers -that can observe the corona to heights as large as ≈6 R .Such imagers include the Sun Watcher with Active Pixels and Image Processing (SWAP) on the Project for On-Board Autonomy 2 (PROBA2) spacecraft (Seaton et al., 2013b;Halain et al., 2013), the Extreme-Ultraviolet Imager on Solar Orbiter (Rochus et al. 2020;Auchère et al. 2023; see also the observations of a prominence eruption to > 6 R reported by Mierla et al. 2022), and the Solar Ultraviolet Imager (SUVI) on NOAA's GOES-R line of spacecraft (Seaton et al., 2020;Darnel et al., 2022).Both SWAP and SUVI have fields of view (FOV) that can reveal features in the EUV corona to heights greater than 2 R , and both have been used in off-pointed campaigns that revealed EUV corona structures to substantially larger heights, as great as ≈ 5 R (e.g.Goryaev et al., 2014;O'Hara et al., 2019;Tadikonda et al., 2019;Seaton et al., 2021).Observing the region between 1.5 -6 R -the middle corona (West et al., 2022b) -is important for a wide variety of scientific questions including coronal mass ejection (CME) initiation, solar wind acceleration, and characterizing the processes that give the corona its large-scale structure.
Upcoming and proposed instruments, including the EUV CME and Coronal Connectivity Observatory (ECCCO: Golub et al., 2020) and the Sun's Coronal Eruption Tracker CubeSat (SunCET; Mason et al., 2021Mason et al., , 2022)), should reveal the EUV middle corona in dramatically more detail than has been possible, even with SWAP and SUVI.All the new observations that these instruments yield point to the great need for image-processing techniques that can help improve the accessibility of interesting features and events buried in these datasets.
In fact, this need is also shared by other data sets, including both those from the Association of Spacecraft for Polarimetric and Imaging Investigation of the Corona of the Sun (ASPIICS) coronagraph (Shestov and Zhukov, 2018;Shestov, Zhukov, and Seaton, 2019) on the formation-flying PROBA3 mission and ground-based eclipse observations (Habbal et al., 2007(Habbal et al., , 2011) ) This and other new instruments will also observe the corona in largely the same region, and will likely encounter many of the same challenges, albeit in other wavelength regimes.
In this article we briefly describe our approach for processing these images to manage dynamic range and reveal features that would otherwise be invisible, which we refer to as the SWAP Filter, thanks to its origins.These techniques have been demonstrated extensively on EUV coronal observations, including observations of large-scale EUV structure and evolution (Seaton et al., 2013a), eruptive flares (Seaton and Darnel, 2018), and off-pointed observations (Seaton et al., 2021), but they have not yet been described in detail in the literature.
Although today this problem is typically addressed using software, techniques to manage dynamic range in coronal imaging have their origins many decades ago, during the era of photographic observations.One of the earliest examples was developed by Gordon Newkirk, who found that photographic images lacked sufficient dynamic range to capture the detail in the corona visible by eye during a total solar eclipse.In response, Newkirk devised an optical radially graded neutral-density filter that successfully compensated for the steep brightness gradient of the corona (Eddy, 1989).Newkirk first deployed his filter at an eclipse in 1966 in Bolivia, generating a stunning image of the corona that arguably revealed more detail than was visible to eye (Newkirk, 1967).Subsequently, Newkirk-style filters were widely used for eclipse photography before digital imaging made it possible to assemble high-dynamic range (HDR) composites using multiple exposure times instead of radial-normalizing filters.
Although the origins of our technique are firmly rooted in the era of spacebased electronic imaging, the fundamental motivation is the same as Newkirk's: during our work with SWAP images, we realized that there was more to see in our images than we could display using simple techniques.In this article, we describe the problem and the simple observation that led to the development of our approach (Section 2), describe the approach itself and some of its configurable parameters (Section 3), and finally highlight a few examples of its applications (Section 4) before making a few concluding remarks (Section 6).

The Origins of Our Filter Technique
Early in the SWAP mission, we realized that SWAP's large FOV might reveal features and evolution in the solar corona that had not been well observed previously in these wavelengths.We began to produce movies composed of deepexposure images spanning full Carrington rotations, using median stackingthat is, computing the median of each image pixel in the temporal direction -of many individual SWAP observations to suppress noise and enhance the visibility of faint features in the middle corona (see Section 2.2 in West et al., 2022a).We quickly realized, however, that without a technique to equalize the dynamic range in the corona, there was no path to take full advantage of these new observations, since we could not clearly visualize the large-scale features in the data.Figure 1 highlights this dilemma, showing an example of a deep-stacked SWAP image from 28 Oct 2014 at 13:21:13 UT, displayed using three different datascaling approaches.The left panel uses a linear scale, which results in an image with too much contrast to be shown on typical displays.A common approach to address this problem is to use a logarithmic scaling (center), but this tends to suppress fine variations too much to visualize anything other than the very large-scale features.A third approach is to raise the image values to a power or other nonlinear function, (e.g.Lupton et al., 2004), which reduces dynamic range enough to reasonably display the image without washing out fine features.We found a good value for SWAP observations was the1 ⁄4 power, but found it was still not adequate to provide a detailed view of the entire corona within the FOV.Thus, while these various data scaling approaches are used widely for narrow-field EUV image data, we found that none were adequate to capture both large-scale structure, which can vary in brightness by up to five orders of magnitude over the FOV, and fine features, which often represent variations of <1%.
A straightforward approach to overcoming this challenge is to rescale the data using a symmetrical radial gradient normalizing filter.Much like Newkirk's analog implementation, these filters are used to neutralize the steep radial falloff of coronal brightness as a function of height above the disk.Typically one generates the filter by sampling the coronal brightness in successive rings at greater heights, and then dividing the coronal image by the resulting array to eliminate the gradient with height, leaving only local variations to be displayed 1 .
Although filters like this have proven highly successful for visible-light coronal imaging, we found that this simple approach was not satisfactory for large-scale EUV images due to the large variation in brightness between polar coronal holes and equatorial streamers and other bright features.
We sought another approach to "normalize" -that is, divide -the data that could be adapted to adjust the degree of data scaling as a function of azimuth around the Sun.(Note that throughout this article we adopt the term normalize to refer to a process in which data are rescaled by division by a function or array.)Because at the time we were working with observations that spanned complete Carrington rotations, we had a large amount of data that revealed the evolution of the Sun and corona over a full solar rotation -or even longerreadily available to work with.
As an initial experiment, we adopted a concept similar to the monthly minimum images (e.g.Section 10 in Morrill et al., 2006), which are used to subtract background (stray light and F-corona) from coronagraph observations, to see whether a similar type of stacking might yield an image that we could use to normalize our data.This approach turned out to be very effective, because longterm averages smooth out most (but not all) local variations and leave a largely smooth -but azimuthally variable -fall-off in intensity that can serve as the basis to normalize the height gradient of an image.
Because the large-scale features in the off-disk corona often persist for many months (Seaton et al., 2013a), long-term image stacks preserve some of the variation of coronal brightness as a function of latitude, which leaves a normalizing filter that acts differently at equatorial latitudes, where bright active regions are often present; at mid-latitudes, where extended coronal fans are often observed; and at the poles, where coronal holes are most common.In our case, we found median-stacking the data to be most effective approach to generating these filters, since the median suppresses both temporal noise and transient events that might otherwise skew the behavior of the filter.
Figure 2 shows an example of a long-term median-stack image generated using all SWAP observations during Carrington Rotation 2156, and a sample of an image (the same image shown in Figure 1) processed using this image as the basis for a normalizing filter.In addition to their value in equalizing the brightness in the off-disk region, an important feature of these filters is the bright bands left by active regions as they transit the solar disk.When used to normalize a full-Sun image, these bands suppress a bit of brightness at latitudes where the most prominent active regions appear, and comparatively amplify quiet-Sun features.The net result is a disk image that is normalized so both bright active regions and dark coronal holes and filaments are well scaled and overall contrast is better than what can be achieved with simpler image normalization techniques.It is important to note that the filtered images should not be used to determine absolute brightness or density on their own, but they are excellent for determining where signal is coming from, and characterizing coronal morphology.
Because our experimental strategy leveraged all of the data for a full solar rotation, it turned out that these filters were extremely useful in processing movies of the evolution of the extended EUV corona over one or more solar rotations.(See the animation accompanying Figure 2 for an example.)However, although this experimental approach yielded dramatic results, it suffered from an important drawback: because the method requires a full Carrington rotation of SWAP observations to generate the filter, it was only practical for processing SWAP's long-duration movies, when the data were already available.Thus we sought an alternative approach that could be more easily deployed for shortduration events such as CMEs and solar flares, but would share most of the benefits that we realized from this technique.The Carrington-rotation-filter method yielded filters with several specific, helpful features: First, they capture the bright limb and the gradual falloff in brightness with increasing height, which largely determines the overall dynamic range of the image.Second, they capture gradual changes as a function of azimuthal angle, allowing the filter to adapt to appropriately normalize both bright extended features at lower latitudes and fainter features in the coronal-hole regions near the poles appropriately.Finally, they achieve good normalization across the disk and off-limb regions, specifically improving contrast on the disk because they include band-like features that result from the passage of either active regions or darker features, such as coronal holes and filaments, at different latitudes.In the Section 3 we describe how we translated these naturally derived features into a filtering approach that requires only a single input image and subsequently applied this filter technique to a variety of observations from SWAP and other instruments.

Implementation and Use
Our strategy uses the three key features that we identified in filters based on Carrington-rotation image stacks discussed above: it captures the fall-off in brightness above solar limb, it captures the slow variation of this gradient as a function of azimuth, and captures the latitudinal distribution of brightness on the solar disk.The method to generate filter arrays is implemented in a SolarSoft IDL function, called p2sw image filter.pro in the PROBA2/SWAP software package.
The generation of the filter is mostly implemented in two straightforward steps, both of which leverage transformations of coordinates that facilitate data smoothing and averaging in directions that are not practical in image coordinates.
We generate the off-disk part of the filter by reprojecting the part of the image above the limb into polar coordinates, with azimuthal angle along the horizontal axis and radial height on the vertical.By smoothing the image in the azimuthal direction with a spatial median filter, which does not affect the radial direction, we generate an array that varies slowly in azimuth, but captures the remaining local gradient with good fidelity.We pad the edges of the transformed image with copies of itself to ensure that the processing wraps smoothly around the Sun at the discontinuity in the azimuthal angle (i.e.0 • /360 • ).After the smoothing process, we transform the processed array back from radial coordinates to Cartesian image coordinates.
We generate the on-disk part of the filter by transforming the disk image to a simple cylindrical map projection, which accounts for the observatory's viewing angle in heliographic latitude (that is, the B 0 -angle).We then compute the median value for each latitude (that is, across the longitudinal direction) and replace all the pixels in the map with their appropriate median.Again, we transform this back to image coordinates.
These two processed parts of the filter are merged into a single image and modulated using a set of user-configurable parameters (described below), yielding a basis image to perform radial filtering that shares its broad characteristics with the Carrington-rotation-stack images that gave us the idea for this approach.
Figure 3 shows the resulting filter and processed image, displayed using the same parameters as Figure 2. The method that we describe here is more sensitive to the specific features visible in the image used to manufacture the filter than the long-term stacking approach is, and thus more strongly modulates the characteristics of these specific features.The resulting image has overall less contrast than the processed image in Figure 2, but in each case both the overall global structure of the corona and its local-scale variations, both on-and off-disk, are far more visible than in any of the images in Figure 1.
The filter's overall sensitivity to local features can be adjusted somewhat by altering the width of the median filter used off-disk, which is one of several userconfigurable parameters in the code.Here we used a value of ±15 • , which we found matches the scale of many typical off-disk coronal features, but which can be adjusted to emphasize different scales.Very large values diminish the filter's ability to adapt azimuthally to changing coronal brightness, while very small values make the filter too sensitive to local variations and suppress all but the most fine-scale features.The filter can also be tuned by an overall smoothing parameter, which is useful both for setting the scale of features to emphasize and suppressing small artifacts that sometimes develop during the filter generation.
Here we used a 2D Gaussian blur with a σ value of four pixels on the final filter in Cartesian space, but both the blurring method and width are adjustable depending on the user's needs.
There are two other very important parameters that must be set to obtain a high-quality, low-noise, processed image that retains the overall natural appearance of typical solar EUV images: a bright limb, bright active regions, and a gradual decrease in brightness with height.(Such aesthetic choices are not necessarily requirements of a generic image filter; they help meet one of the Figure 3. Filter image generated by our method displayed at 1 ⁄4 power (left) and the result of using this filter on our test image from 28 Oct 2014 13:21:13 UT (right).Though specific features differ somewhat between this and the processed image in Figure 2, the overall results are comparable without the requirement of large quantities of data.
objectives of our specific filter: namely to generate images and movies for publicoutreach materials.)The filter, as derived, acts uniformly across the entire image.If applied with no adjustments, the brightest features will be re-normalized to have the same brightness as the faintest features, while the faintest features, which are also the noisiest in most images, will be enhanced to match the brightness of all other features.The result is an image that both lacks contrast -or is excessively "flat" in appearance -and excessively noisy in its faintest regions.We therefore must tune the filter slightly to avoid both of these outcomes.This is straightforwardly achieved by modulating the overall filter array [F] to yield a final filter [F ], in the following way, Here, t 0 is a constant offset, which we refer to as the tapering parameter, that depends on the image background brightness (or can be adjusted by the user) and reduces the amount by which the filter enhances the faintest pixelsand specifically the noise in these regions.c 0 is a factor that we refer to as the crush factor and is defined such that 0 < c 0 < 1.This factor reduces the overall normalization effect of the filter, particularly on bright values.Careful tuning of these parameters can significantly improve the overall image result, but any filter that amplifies the faintest part of the image without also applying noise reduction is likely to also enhance noise in faint regions as well, thus there is always a trade-off to make between how aggressively to filter regions where the signal-to-noise ratio is low.
For the SWAP images in both Figures 2 and 3, we set the constants in Equation 1 to be t 0 = 1.5 and c 0 = 0.75.In our experience, the distribution of brightness in typical wide-FOV EUV images is such that the value of t 0 should be roughly the median of all pixels in the filter image, but manual fine-tuning SOLA: manuscript.tex;22 May 2023; 0:33; p. 8 The SWAP Filter for EUV Solar Images  is often required to obtain the user's desired appearance for the image.Note that the value of t 0 differs between the two techniques that we illustrate in this article, so we select t 0 = 1.5, which is larger than the median of the Carringtonrotation-stack and smaller than the median of our method's filter, for both cases to ensure the overall image is treated roughly equally by both methods.
Figure 4 illustrates why it is necessary to rescale the filter using c 0 before processing an image using some cross-sectional cuts through the data in our sample image.If we do not compress (or "crush") the filter dynamic range, dividing by the filter will completely flatten all variation in brightness with height above limb, which results in an image that looks significantly modified from the original observations to many viewers.Crushing this function using c 0 = 0.75 has little effect on how the filter treats faint features, but reduces somewhat how much it flattens bright features, resulting in an image that retains some -but much less -of the gradient from limb to the large heights and is both appropriate for display and more natural in appearance.
Figure 5 illustrates the results of neglecting these two parameters with two processed images.This yields images that either lack contrast (the c 0 = 1.0 case; left) or have excessive enhancement of the background (the c 0 = 0.75, t 0 = 0. case; right).The overall result of appropriate selection of these parameters, as shown in Figure 3, yields a final product that achieves balance between contrast, dynamic range, and resemblance to the original observation without excessively amplifying noise in dark regions.
Two final choices remain before the processed image can be displayed: how to scale the final filtered image for display, if at all, and how to set the appropriate limits to clip the displayed image's dynamic range.We found that our filtered SOLA: manuscript.tex;22 May 2023; 0:33; p. 9 SOLA: manuscript.tex;22 May 2023; 0:33; p. 10 images still have enough variation with height that it is useful to adjust the image scaling, and further, that scaling the image to the 1 ⁄4 power yields a subjectively aesthetically appealing result -though, as we discuss in Section 2, a variety of nonlinear functions can work depending on the specific need.For many of the processed images we have generated, we find that a good display range for the filtered and scaled image is between 0.5 -2.0 renormalized counts.We used this combination for the images in this article, but users are encouraged to experiment and find display approaches that suit their particular needs and data.Some examples in the p2sw image filter.procode itself illustrate how to implement this and display an image in IDL.
In the example above, the fully renormalized image [I ] is then described by where, I 0 is the unprocessed, calibrated image and F is the filter generated in Equation 1.I is re-scaled as a byte-array using the values indicated.

Applications
Although we have refined these filtering techniques somewhat since first developing them, largely to make them more robust to anomalies in image data, they have been used effectively on a variety of EUV image data and applications, including tracking the long-term evolution of large-scale features in SWAP data (Seaton et al., 2013a), observations of eruptive solar flares in AIA (Seaton, Bartz, and Darnel, 2017) and SUVI (Seaton and Darnel, 2018;Veronig et al., 2018), as well as the EUV waves that often accompany these events (de Koning, Pizzo, and Seaton, 2022), and even on simulated observations of eruptions (Mason et al., 2022) from planned instruments (Mason et al., 2021).These techniques have proven especially valuable for ultra-wide-field observations of the EUV corona (Tadikonda et al., 2019;Seaton et al., 2021), where the dynamic range from limb to 3 R can reach almost five orders of magnitude.Figure 6 contrasts the detailed structure visible in wide-field EUV images from SUVI when processed with our filter with the minimal structure visible when these images are displayed using log-scaling.Enhanced observations such as these make it possible to characterize the very faint dynamics associated with the origins of the solar wind's embedded structure (Chitta et al., 2022).
Figure 7 and the accompanying animation show how processing with our filter makes visible multi-scale dynamics in wide-field EUV observations.The movie makes visible both small, ubiquitous, jet-like features and several large eruptions near the east limb.These faint variations would not be detectable compared to the steep radial gradient in the brightness without this processing.(See Seaton et al., 2021, for a complete discussion of the processing required for these SUVI mosaic images and the origins of artifacts in the resulting images and movies.)

Comparison to Other Techniques
Since the development of the SWAP Filter, early in PROBA2's mission, a number of powerful, new, image-processing techniques have also been published (e.g.Druckmüllerová, Morgan, and Habbal, 2011;Druckmüller, 2013;Morgan and Druckmüller, 2014;Auchère et al., 2023).The fundamental idea underpinning our approach was very simple: to adapt well-known radial-normalizing filter approaches -dating all the way back to Newkirk (1967) -to be more responsive to azimuthal variations, which are particularly pronounced in the EUV corona.The filter's primary purpose is thus to reduce dynamic range -as a function of altitude within localized features off-disk, and as a function of latitude on-disk -to improve display of data that span many order of magnitude in radiance.Other methods filter specific spatial frequencies to emphasize and sharpen features with specific characteristics, using, for example, wavelets (e.g.Stenborg andCobelli 2003 andAuchère et al. 2023) or multi-scale Gaussian normalization (Morgan and Druckmüller, 2014).Although these techniques can be tuned to match features on many scales, they are particularly powerful for enhancing the visibility and clarity of low-contrast features on fine scales (see, for example, Figure 3 in Auchère et al. (2023) or Figure 6 in Morgan and Druckmüller (2014)).
In contrast, the SWAP filter's strength lies in its ability to make visible large-scale coherent features without trading off localized contrast -as with logarithmic or power scaling (compare Figure 1 and Figure 3).In this sense, the SWAP filter has much in common with other radial normalizing filters such as the NRGF and FNRGF, which are optimized specifically for off-disk image processing.The SWAP filter marries its off-disk processing technique to a related on-disk filtering technique, but is nonetheless fundamentally similar in approach to these methods.
Figure 8 presents examples of the same SWAP image shown elsewhere in this article processed with the MGN and NRGF filters, implemented in the sunkit-image Python package2 .While the MGN-processed image clearly reveals coherent structures on multiple scales, the overall appearance of the image -particularly off of the solar disk -is somewhat different from the SWAP filter.The spatial filtering of the MGN suppresses the low-frequency background brightness in the image, strongly sharpening features on smaller scales and rendering them with largely uniform intensity across the full FOV.Like the SWAP filter, the processing of the disk image improves the visiblity of loops and large-scale features such as the filament in the northwest part of the Sun, but the MGN makes small-scale variations within these features much more prominent.Furthermore, the MGN can be flexibly tuned to both highlight features on specific scales, by adjusting the specific spatial scales isolated by the filter, and remix the relative importance of each frequency in the image, by adjusting weights for each individual spatial component in the final image.(Here we set the MGN to isolate a range of spatial scales ranging from just over a single pixel to 40 pixels, with the strongest weights for the large-scale features, but alternate scalings are easily achieved via configurable parameters.) The NRGF image, on the other hand, more closely resembles the processed image in Figure 3, because the approach used by the NRGF filter is similar to the off-disk approach used in the SWAP filter.A notable difference is that the NRGF processing -which is described by Morgan, Habbal, and Woo (2006) as being specifically developed to address the steep brightness gradient in coronagraph images -is not optimized for use on disk.This region is often omitted from NRGF-processed images (or, in the case of coronagraphic images, is simply unavailable).Nonetheless, because both the SWAP filter and NRGF were both developed to address the radial gradient in coronal images, it is not surprising that they emphasize many of the same large-scale coronal features.
To more quantitatively highlight differences in the SWAP filter's performance compared to other image scaling and filtering techniques, we compute the azimuthally averaged power spectrum for the images processed with the SWAP filter (see Figure 3) and compare the resulting spectrum to those of the log and 1 ⁄4-power scaling (see Figure 1), as well as the MGN filter (see Figure 8 The MGN enhances structures on relatively small scales, and thus has the most power in this range, but it fully suppresses the largest-scale structures, and has the least power for these scales.The 1 ⁄4-power, log-scaled, and SWAP-filtered images are all are similarly behaved, with the SWAP filter emphasizing small and mid-range scales, and thus having the most power in this range.The log-scaled image suppresses the mid-range frequencies, but equalizes brightness across the whole image, and thus has the most power at large scales.Surprisingly, the logscaled and SWAP-filtered images are very similarly behaved at high frequencies.In fact, this appears to be because they enhance off-disk noise -a very highfrequency feature -compared to the 1 ⁄4-power filter, which attenuates the off-disk brightness much more than these other methods.
Overall, these power spectra show more quantitatively what we have, up until now, mainly illustrated qualitatively: the SWAP filter results in images that more closely resemble scaled, but otherwise unprocessed, solar images, while also improving the image contrast over these simpler methods for features across a wide range of scales.

Conclusions
Here we have described the SWAP Filter: a newly available azimuthally varying radial normalizing filter, specifically designed for wide-field EUV images of the solar corona.Although simple, the filter is highly effective for producing displayable EUV images that preserve both large-scale EUV features and small-scale variations within the global coronal structure.The filter has been demonstrated on a large variety of EUV image data, and is now publicly available via the PROBA2/SWAP SolarSoft IDL software package.Performance tests SOLA: manuscript.tex;22 May 2023; 0:33; p. 15 using an 2.40 GHz Intel Core i9-9980HK demonstrate that filter generation with default parameters typically takes less than 0.3 s for a 1024×1024 SWAP image, so the technique is fast and computationally light compared to other similar techniques, making this a useful approach when computational resources are limited or for large datasets.
Image processing tools such as this one are especially important for studies of the middle corona, as observations that connect this region to the solar disk inherently contain very large brightness gradients that must be reduced in order to effectively display the entire corona.Numerous approaches are available, each of which offers various advantages for specific applications.Our filter is fast to apply and it yields images that do not appear extensively processed, and is well suited for the display of data that emphasizes global-scale features, and for generating outreach images and movies that resemble the appearance of inner coronal images that have been available for decades, and thus feel familiar to viewers.

Figure 2 .
Figure 2. (left) Running median of all SWAP images during Carrington Rotation 2156 displayed at 1 ⁄4 power and (right) the result of using this running median as an image filter to normalize the same test image as in Figure 1 (28 Oct 2014 13:21:13 UT).An accompanying animation in the Electronic Supplementary Material shows the results of applying this type of filter to a full solar rotation's worth of observations.

Figure 4 .
Figure 4. Brightness of a cross-section of our test from South Pole to the edge of the FOV compared to the unprocessed filter in Figure 3 (green) and the same rescaled using c 0 = 0.75 (gold).

Figure 5 .
Figure 5.Comparison of a filtered image with c 0 = 1.0 and t 0 = 1.0 (top left; i.e. no rescaling, but including low-end tapering) resulting in a "flat" looking result and with c 0 = 0.75, but t 0 = 0.0 (top right; rescaling, but no low-end tapering), which leads to an over-enhancement of background brightness and low-level temporal noise, and (bottom) the optimally processed image from Figure 3 (28 Oct 2014 13:21:13 UT).

Figure 6 .
Figure 6.Comparison between log-scaled and filtered data from a SUVI offpoint campaign in 2018.Concentric dotted circles indicate heights of 2 and 3 R .Solar North is indicated by the tick mark near the top of the solar disk.Adapted from Seaton et al. (2021), used with permission.

Figure 7 .
Figure 7. SUVI image of extended EUV corona from a campaign during April 2021, showing an erupting prominence cavity, indicated by arrow near the northeast limb.See also the accompanying animation in the Electronic Supplementary Material.

Figure 8 .
Figure 8. (left) MGN-filtered version and (right) NRGF-filtered version of the image shown in Figure 3, illustrating similarities and differences in the same image when processed with different techniques.

Figure 9 .
Figure9.Azimuthally averaged power spectra of images processed with the SWAP filter compared to log-scaling, 1 ⁄4-Power image scaling, and MGN processing.