A2 A program to clean ithe artifacts

Current computing technology makes it possible to esti-x A mate and remove most of the scanning and reconstruction H5 artifacts (Byrne, 2002). A two-stage process is used to clean pen the images in this book. First, the systematic artifacts such Ap as white lines between framelets and brightness variations across framelets are measured and removed. Second, a filter suppresses remaining variable striping patterns. The cleaned image is shown in Figure A.2.

The software programs for each of the stages were written in Visual Basic® and run in a Windows 98 environment. They use utility routines from "Visual Basic Graphics® Programming, a book by Rod Stevens (Stevens, 1999).

Extrapolate or interpolate those framelet edges that were not identified by search processes Scan each framelet edge and, for each horizontal pixel coordinate, estimate the excess brightness by comparison with the immediately adjacent pixels. Then remove the excess brightness from the pixels of the affected scan lines

Determine the average linearized brightness ratio of the pixels, as a function of their relative position between fra-melet edges, for the entire mosaic image. A model of the nonlinear contrast function is used in the linearization process

8. Correct the brightness of each pixel for the average normalized brightness ratio corresponding to its position between framelet edges. The contrast model is used in the correction.

Specific Techniques of Stage 1

The program addresses several problems that are specific to the mosaics.

Finding the edges of framelets. Usually, framelet edges are represented by bright lines caused by light shining between the framelet filmstrips as the mosaics were laid up. Although they are narrower than the distance between scans of the digital images, the lines appear in either two or three adjacent lines of pixels because the scan spot used in creating the digital images was larger than the distance between scans. The scan lines are usually tilted with respect to the framelet edges by two or more pixels across a frame and have curvature of the order of one pixel across a frame. In many frames, the white lines are obscured by saturation or overcome by image signal in large parts of the frame. Development artifacts and some valid signals can produce false segments of bright lines that are not at framelet edges.

The approach taken to find these lines is to search in a band of scan lines 17 pixels high, looking in 17 x 21 blocks of pixels for a bright white line segment of either two or three pixels in width. A quadratic best-fit line is calculated using the centers of the line segments. Successive fits are calculated, eliminating those centers that are most off the line until all remaining centers are within half a pixel of the line. If the remaining centers are at least 25% of the possible centers, the

Steps of Stage 1

The steps performed by the program that measures and removes systematic artifacts (Stage 1) are the following:

1. Search from the top of the mosaic image for the first identifiable bright line between framelets (a framelet edge)

2. Jump half a typical 40-pixel framelet width and search for the next identifiable framelet edge

3. Repeat the jump-and-search process until the bottom of the image is reached

4. Calculate the precise framelet width of the image (which varies between 38 and 43 pixels on the currently sampled images) by averaging the framelet width between successive identified framelet edges. Successive edges are those separated by about one typical framelet width

Figure A.1. LO4-187H2, Mare Orientale and Montes Rook.

40.00

40.00

Period (Pixels)

Period (Pixels)

Figure A.3. Spectral density of L04-140H3 for periods near the framelet width (about 40 pixels). The solid line shows the spectrum of the raw photo and the dashed line shows the spectrum after removal of the systematic artifacts.

line is accepted as a framelet edge. If not, the band is moved down the picture to find an edge.

Removing each framelet edge line. For each horizontal pixel index, the pixels above and below the calculated vertical value of the best-fit line are examined to determine the probable value of the excess brightness of the line and a weighted value is subtracted from the brightness of the two or three pixels near the line. This process preserves more detail in the image than simply averaging the brightness around the line.

Contrast model. The streaking artifact was applied in the spacecraft on a low contrast image, essentially linear in its relation between photographed brightness and density.

Figure A.2. LO4-187H2 after removal of systematic artifacts and filtering.

However, the atlas images were printed at high contrast and are therefore nonlinear. Thus, in order to determine and compensate for streaking, it is necessary to linearize the brightness by reversing the contrast function. An empirical ^J^ study was conducted, varying the assumed contrast function LJ5J and minimizing the residual fundamental and second and A third harmonics of the artifacts at the framelet frequency. £ The resulting contrast function has a linear contrast gain nd of 3 between output brightness values of 0.1-0.9 of the full a range and exponential curves at each end of the brightness < range.

Special handling. A few images require special techniques. For example, a feature has been added to compensate for an alternating pattern of darker and lighter framelets in some of the images.

Quantitative Results of Stage 1

A spectral analysis program has been written to compare the spectra of the input and output images to provide a quantitative measure of improvement. Figure A.1 shows the fine-scale normalized vertical spectral density of L04-140H3 before and after processing.

The removal of systematic artifacts greatly improved the appearance of the pictures. The venetian blind effect disappears, providing a subjective effect of more direct visual contact with the lunar surface. However, traces of the framelet artifacts remain in parts of some of the images.

Steps of Stage 2

Stage 2 is a two-dimensional filter process that performs the following steps:

1. Generate the two-dimensional Fourier transform of the image o

2. Multiply the transform by a two-dimensional filter that suppresses streaks, that is, patterns that extend horizontally more than vertically. This filter is formed from the combination of a low-pass filter and a high-pass filter, as suggested by Lisa Gaddis of USGS (Gaddis, 2001).

3. Return to the spatial domain by an inverse two-dimensional Fourier transform.

The remaining artifacts are substantially removed by application of the two dimensional filter (see Figure A.2). Careful comparison of filtered images with the originals verifies that shadows, topographic variations, and albedo variations such as associated with rays are preserved. This is true even when the digital images are examined under magnification that resolves the pixels.

Specific Techniques of Stage 2

Fourier Transform and inverse transform: A fast-Fourier-transform algorithm runs in the vertical direction and stores arrays of the real and imaginary parts of the transform. The same algorithm runs in the horizontal direction on these arrays to create the real and imaginary two-dimensional arrays that are now in the frequency domain. The same process is used for the inverse transform. The fast-Fourier-transform algorithm used is limited to image arrays that have dimensions M and N that must be integer powers of two; therefore, the dimensions of the actual images are built up to the next power of two by adding virtual blank brightness data.

Filter characteristic: The filter has a low-pass characteristic in the vertical (north-south) direction to suppress the frame-let frequency and its harmonics. In the horizontal direction (east-west), it has a high-pass characteristic to reject patches of horizontal stripes while permitting topographic detail to remain.

Both filters are of the second-order Darlington type. A second-order low-pass Darlington filter has the form

A high-pass Darlington filter has the form 1

Examination of the two-dimensional Fourier transform of an input image (whose systematic artifacts have been removed) shows strong noise at horizontal frequencies whose period is less than 150 pixels and vertical frequencies whose period is more than 25 pixels.

A compound (two-dimensional) "Normal" filter with the following form is used to suppress the noise:

where fx is the frequency in the x direction (east-west) and fy is the frequency in they direction (north-south).

The "Normal" filter is used as the default, and is effective for most of the pictures. If streaks remain after application of the "Normal" filter, a "Strong" filter is used instead. This happens most often when large areas of dark sky or shadow beyond the terminator appear in an image. The form of the "Strong" filter is

A.3. Summary where f is the frequency and T is the period of the frequency that reduces the filter function G(f) to 0.5.

The Lunar Orbiter photographs have been used as the major reference for topographic information for the last 40 years. They commonly appear in papers about the Moon even today in their original state, artifacts included.

Throughout this time, researchers have learned to distinguish the artifacts from the topographic information. It is my hope that these cleaned photographs will make the nature of the Moon more available to both nonprofessionals and to new students of the Moon.

Was this article helpful?

0 0

Post a comment