基于ENVI的极化SAR数据处理流程介绍

ENVI Tutorial: Polarimetric SAR Processing andAnalysisTable of ContentsO VERVIEW OF T HIS T UTORIAL (2)Background: SIR-C/SAR (2)P REPARE SIR-C D ATA (3)Optional: Read a SIR-C CEOS Data Tape (3)Optional: Multilook SIR-C Data (3)S YNTHESIZE I MAGES (4)Default Polarization Combinations (4)Other Polarization Combinations (4)Display Images (5)Define ROIs for Polarization Signatures (6)Extract Polarization Signatures (6)Adaptive Filters (8)Slant-to-Ground Range Transformation (9)Preview CEOS Header (9)Resample Image (9)Texture Analysis (10)Create Color-coded Texture Map (10)Image-Map Output (11)Overview of This TutorialThis tutorial demonstrates the use of ENVI’s tools for analyzing polarimetric synthetic aperture radar (SAR) data. You will multilook Spaceborne Imaging Radar-C (SIR-C) from Death Valley, California; synthesize images, define ROIs for (and extract) polarization signatures, use adaptive filters, perform slant-to-ground range transformation, use texture analysis, and create an output image-map.Files Used in This TutorialENVI Tutorial Data DVD: envidata\ndv_sircFile Descriptionndv_l.cdp L-band SIR-C subset in ENVI compressed data product (.cdp) formatpol_sig.roi Region of interest (ROI) fileBackground: SIR-C/SARSIR-C is a polarimetric SAR instrument that uses two microwave wavelengths: L-band (24 cm) and C-band (6 cm). The SIR-C radar system was flown as a science experiment on the Space Shuttle Endeavor in April (SRL-1) and October 1994 (SRL-2), collecting high-quality SAR data over many sites around the world. (A second radar system, XSAR, was also flown on this mission, but these data are neither discussed nor processed here.) Additional information about SIR-C is available on the NASA/JPL Imaging Radar Home Page at /.Prepare SIR-C DataThe data used in this tutorial are a subset of L-band Single Look Complex (SLC) SIR-C data that cover the northern part of Death Valley, including Stovepipe Wells, a site of active sand dunes and extensive alluvial fans at the base of mountains. These data were preprocessed by reading and subsetting from tape and multilooking (averaging) to 13 m square pixels. The data are provided in ENVI compressed data product (.cdp) format. This non-image format is similar to the tape format and cannot be viewed until images are synthesized for specific polarizations.The first two functions described in this example—reading the data tape and multilooking—were already applied to the SIR-C data. The steps are provided here only for completeness if you want to learn more about the processes. Skip to the Synthesizing Images section if you are not interested in data preparation.Optional: Read a SIR-C CEOS Data Tape1.From the ENVI main menu bar, select File→Tape Utilities→Read Known Tape Formats→SIR-C CEOS.The SIRC Format - Load Tape dialog appears.2.Enter the tape device name and accept the default record size of 65,536. Click OK. The tape is scanned todetermine what SIR-C files it contains. A dialog appears to let you select the desired datasets. By default, ENVI reads all of the data files on the tape.3.If you do not want to read all of the data files, click Clear, then select the check box next to each desired file.Click OK.4.You can independently subset and multilook the selected data files as they are being read from tape. However,you should perform multilooking on disk (unless you have insufficient disk space) as this function is extremely slow from tape.5.Select a filename, then click Spatial Subset or Multi-Look to enter parameters for the data file. Enter anoutput filename. Each input file must have an output filename. By convention, the output filenames should take the form filename_c.cdp and filename_l.cdp for the C- and L-bands, respectively. The SIR-C data areread from the tape, and one compressed scattering matrix output file is created for selected each dataset.Optional: Multilook SIR-C DataMultilooking is a method for reducing speckle noise in SAR data and for changing the size of a SAR file. You can multilook SIR-C data to a specified number of looks, number of lines and samples, or azimuth and range resolutions. The SIR-C file used in this tutorial was a single-look dataset with a range resolution of 13 m and an azimuth size of 5 m. Multilooking has already been performed in the azimuth direction to make 13 m square pixel sizes. Instructions on multilooking are included here only for completeness.1.From the ENVI main menu bar, select Radar→Polarimetric Tools→Multilook Compressed Data→SIR-C Multilook. An Input Data Product Files dialog appears.2.Click Open File and select an input file. ENVI detects whether the file contains L- or C- band data and displaysthe filename in the appropriate field of the dialog. Click OK.3.Select the file to multilook by selecting the check box next to the name. You can select multiple files.4.Enter any one of three values—number of looks, number of pixels, or pixel size—and the other two are calculatedautomatically. Integer and floating-point number of looks are supported.5.Enter the desired Samples (range) and Lines (azimuth) values.6.Enter a base filename and click OK.Synthesize ImagesThe SIR-C quad-polarization data provided with this tutorial and available on tape from JPL are in a non-image, compressed format. Accordingly, images of the SIR-C data must be mathematically synthesized from the compressed scattering matrix data. You can synthesize images with any transmit and receive polarization combinations you want.1.From the ENVI main menu bar, select Radar→Polarimetric Tools→Synthesize SIR-C Data. An InputProduct Data Files dialog appears.2.Click Open File. A file selection dialog appears.3.Navigate to envidata\ndv_sirc and select ndv_l.cdp. Click Open. When the filename appears in theSelected Files L: field, click OK. The Synthesize Parameters dialog appears.Default Polarization CombinationsFour standard transmit/receive polarization combinations—HH, VV, HV, and TP—are listed in the Select Bands to Synthesize list of the Synthesize Parameters dialog. By default, all of these bands are selected to be synthesized.1.Enter ndv_l.syn in the Enter Output Filename field.2.Click the Output Data Type drop-down list and select Byte. This scales the output data to byte values. (If youwill be performing quantitative analysis, the output should remain in floating-point format.) Click OK. Afterprocessing is complete, four bands corresponding to the four polarization combinations are added to the Available Bands List.Other Polarization CombinationsThe transmit and receive ellipticity and orientation angles determine the polarization of the radar wave used to synthesize an image. The ellipticity angle falls between -45 and 45 degrees and determines the “fatness” of the ellipse. The orientation angle is measured with respect to horizontal and ranges from 0 to 180 degrees. You can synthesize images of non-default polarization combinations by entering the desired parameters as follows.1.From the ENVI main menu bar, select Radar→Polarimetric Tools→Synthesize SIR-C Data. The filendv_l.cdp should still appear in the Selected Files field. Click OK. The Synthesize Parameters dialog appears.2.Enter -45 in both the Transmit Ellip and Receive Ellip fields and 135 in the Transmit Orien and ReceiveOrien fields.3.Click Add Combination. This will produce a right-hand circular polarization image.4.Enter 0 in both the Transmit Ellip and Receive Ellip fields and 30 in the Transmit Orien and Receive Orienfields.5.Click Add Combination. This will produce a linear polarization with an orientation angle of 30 degrees.6.Click Clear under the list of polarization combinations to turn off synthesis of the standard polarization bands,which have already been generated.7.Select the Yes radio button for Output in dB? This will produce images that are in decibels with values typicallybetween –50 and 0.8.In the Enter Output Filename field, enter ndv_l2.syn and click OK. After processing is complete, two bandscorresponding to the polarization combinations are added to the Available Bands List.Display Images1.In the Available Bands List, select [L-TP] under ndv_l.syn and click Load Band. The SIR-C, L-band, total-power image appears in a new display group.2.From the Display group menu bar, select Enhance→Interactive Stretching. A histogram plot windowappears, which shows the current stretch (between the vertical dotted lines on the input histogram) and thecorresponding DN values in the text fields.3.Drag the dotted vertical lines to change the stretch, or enter the desired DN values into the appropriate fields.4.Enter 5 in the left Stretch field and 95 in the right field.5.From the histogram menu bar, select Stretch Type→Gaussian. Click Apply. A Gaussian stretch is appliedwith a 5% low and high cutoff.6.Generate and compare linear and square-root stretches.7.To display a color composite, select the RGB Color radio button in the Available Bands List. Select[L-HH], [L-VV], and [L-HV] in sequential order under ndv_l.syn.8.Click Display #1 and select New Display. Click Load RGB to display the HH band in red, VV in green, and HVin blue. The color variations in the images are caused by variations in the radar reflectivity of the surfaces. The bright areas in the sand dunes are caused by scattering of the radar waves by vegetation (mesquite bushes). The alluvial fans show variations in surface texture due to age and composition of the rock materials.9.Adjust the stretch as desired (Gaussian and square-root stretches work well on all three bands).10.Close the histogram plot window and Display #2 when you are finished. Keep Display #1 open for later exercises.Define ROIs for Polarization SignaturesYou can extract polarization signatures from a SIR-C compressed scattering matrix for a region of interest (ROI) or a single pixel in a polarimetric radar image. Define ROIs by selecting pixels or by drawing lines or polygons within an image.1.From the Display group menu bar, select Overlay→Region of Interest. An ROI Tool dialog appears.2.Four ROIs were previously defined and saved for use in extracting polarization signatures for this tutorial. Fromthe ROI Tool dialog menu bar, select File→Restore ROIs. A file selection dialog appears.3.Select pol_sig.roi. A dialog box appears, stating that the regions were restored. Click OK.4.Regions named veg, fan, sand, and desert pvt appear in the table in the ROI Tool and are drawn in the displaygroup.5.To draw your own ROI, select ROI_Type→Polygon, Polyline, or Point from the ROI Tool menu bar.6.Click New Region, enter a name, and choose a color.Draw polygons by clicking the left mouse button in the display group to select the endpoints of linesegments, or by holding down the left mouse button and moving the cursor for continuous drawing. Clickthe right mouse button once to close the polygon and a second time to accept the polygon.Draw polylines in the same manner as polygons. Click the left mouse button to define the line endpointsand click the right button to end the polyline and a second time to accept the polyline.Point mode is used to select individual pixels. Click the left mouse button to add the pixel currently underthe cursor to the ROI.You can select multiple polygons, lines, and pixels for each ROI.7.Repeat Step 6 to draw a second ROI. You can save the ROIs to a file and restore them later by selecting File →Save ROI from the ROI Tool dialog menu bar.Extract Polarization SignaturesPolarization signatures are 3D representations of the complete radar scattering characteristics of the surface for a pixel or average of pixels. They show the backscatter response at all combinations of transmit and receive polarizations and are represented as either co-polarized or cross-polarized. Co-polarized signatures have the same transmit and receive polarizations. Cross-polarized signatures have orthogonal transmit and receive polarizations. Polarization signatures are extracted from the compressed scattering matrix data using the ROIs for pixel locations. Polarization signatures are displayed in viewer dialogs, as shown on the next page. To extract your own polarization signatures, perform the following steps.1.From the ENVI main menu bar, select Radar→Polarimetric Tools→Extract Polarization Signatures→SIR-C. The filename ndv_l.cdp should appear in the Input Data Product Files dialog. If not, click Open File and select this file. Click OK. The Polsig Parameters dialog appears.2.Select the four ROIs (veg, fan, sand, and desert pvt) by clicking Select All Items.3.Select the Memory radio button and click OK. Four Polarization Signature Viewer dialogs appear, one for eachROI. The polarization signatures are displayed as 3D wire mesh surface plots and as 2D gray scale images. The X and Y axes represent ellipticity and orientation angles, respectively. You can selectively plot the vertical axis as intensity, normalized intensity, or dB by selecting Polsig_Data from the Polarization Signature Viewer dialog menu bar.4.Polarization signature statistics appear at the bottom of each Polarization Signature Viewer dialog. Notice therange of intensity values for the different surfaces. The smoother surfaces (sand and desert pvt) have low Z values. The rough surfaces (fan and veg) have higher Z values. The minimum intensity indicates the pedestal height of the polarization signature. The rougher surfaces have more multiple scattering and therefore higher pedestal heights than the smoother surfaces. The shape of the signature also indicates the scatteringcharacteristics. Signatures with a peak in the middle show a Bragg-type (resonance) scattering mechanism.5.In any given Polarization Signature Viewer dialog, change the Z-axis by selecting Polsig_Data→Normalizedfrom the Polarization Signature Viewer dialog menu bar. This normalizes the signature by dividing by itsmaximum; the signature is plotted between 0 and 1. This representation shows the difference in pedestal heights and shapes better, but it removes the absolute intensity differences.Alternately, select Polsig_Data → Co-Pol and Cross-Pol to toggle between co-polarized and cross-polarized signatures.e the left mouse button to drag a 2D cursor on the polarization signature image on the right side of the plot.Note the corresponding 3D cursor in the polarization plot.7.Click-and-drag any axis to rotate the polarization signature.8.You can optionally output the signatures to a file or printer by selecting File→Save Plot As or File→Printfrom the Polarization Signature Viewer dialog menu bar.9.Close the Polarization Signature Viewer and ROI Tool dialogs when you are finished.Adaptive FiltersAdaptive filters are used to reduce the speckle noise in a radar image while preserving the texture information. Statistics are calculated for each kernel and used as input into the filter, allowing the filter to adapt to different textures within the image.1.From the ENVI main menu bar, select Radar→Adaptive Filters→Gamma. A Gamma Filter Input File dialogappears with a list of open files. You can apply a filter to an entire file or to an individual band.2.In the Gamma Filter Input File dialog, click the Select by toggle button to choose Band.3.Select [L-HH] under ndv_l.syn and click OK. The Gamma Filter Parameters dialog appears.4.Accept the default values, and select the Memory radio button. Click OK.5.In the Available Bands List, click Display #1 and select New Display. Select the Gray Scale radio button,select the new band name (Gamma), and click Load Band.6.From the Display group menu bar, select Enhance→[Image] Square Root.7.In the Available Bands List, click Display #2 and select Display #1. Select [L-HH] under ndv_l.syn, andclick Load Band.8.From the Display #1 menu bar, select Enhance→[Image] Square Root.9.From any Display group menu bar, select Tools→Link→ Link Displays. The Link Displays dialog appears.Click OK to link the gamma-filtered L-HH image (Display #2) with the original L-HH image (Display #1).10.Click in an Image window to toggle between the two images, using the dynamic overlay feature. The figure belowshows a portion of the original image (left) and the gamma-filtered image (right).11.Close Display #2 when you are finished. Leave Display #1 (ndv_l.syn) open for the next exercise.Slant-to-Ground Range TransformationA radar system looks to the side and records the locations of objects using the distance from the sensor to the object along the line of sight, rather than along the surface. An image collected using this geometry is referred to as a slant range image. Slant range radar data have a systematic geometric distortion in the range direction. The true, or ground range, pixel sizes vary across the range direction because of the changing incident angles. This makes the image appear compressed in the near range, relative to what it would look like if all of the pixels covered the same area on the ground. Slant-to-ground range correction for SIR-C is performed on synthesized images. In other words, the correction is not performed on the entire SIR-C compressed data product file. However, this file does store the required information in the CEOS header about the sensor orientation.Preview CEOS Header1.From the ENVI main menu bar, select Radar→Open/Prepare Radar File→View Generic CEOS Header. Afile selection dialog appears. You must select the original unsynthesized data file from which to extract thenecessary information.2.Select ndv_l.cdp and click Open. A CEOS Header Report dialog appears. Scroll down and note that the linespacing (azimuth direction) is 5.2 m, while the pixel spacing (slant range direction) is 13.32 m. Close the CEOS Header Report dialog when you are finished reviewing it.Next, you will use the Slant-to-Ground-Range function to resample the image to square 13.32 m pixels, thusremoving slant range geometric distortion.Resample Image3.From the ENVI main menu bar, select Radar→Slant to Ground Range→SIR-C. A file selection dialogappears.4.Select ndv_l.cdp and click Open. The Slant Range Correction Input File dialog appears.5.Select ndv_l.syn and click OK. The Slant to Ground Range Correction Dialog appears. ENVI automaticallypopulates the Instrument height (km), Near range distance (km), and Slant range pixel size (m) fields withinformation from the CEOS header.6.Enter 13.32 in the Output pixel size (m) field to generate square ground-range pixels.7.From the Resampling Method drop-down list, select Bilinear.8.In the Enter Output Filename field, enter ndv_gr.img. Click OK. The input image is resampled to square13.32 m pixels. Four new bands appear in the Available Bands List. Band 1 of the resampled image correspondsto the L-HH band of the original, slant-range image (ndv_l.syn), Band 2 corresponds to L-VV, etc.9.In the Available Bands List, click Display #1 and select New Display.10.Select a band from the resampled image and click Load Band. The resampled image appears in Display #2.Make sure Display #1 (ndv_l.syn) shows the corresponding polarization band.pare the two images.12.When you are finished comparing images, close Display #2. Keep Display #1 (ndv_l.syn) open for the nextexercise.Texture AnalysisTexture is a measure of the spatial variation in the gray levels in the image, as a function of scale. ENVI calculates texture based on a processing window size you specify. The texture measures demonstrated in this tutorial are Occurrence Measures, including data range, mean, variance, entropy, and skewness. These terms are explained in ENVI Help. Texture is best calculated for radar data with no resampling or filtering applied.1.From the ENVI main menu bar, select Radar→Texture Filters→Occurrence Measures. A Texture InputFile dialog appears.2.Click the Select By toggle button to choose Band. Select [L-HH] under ndv_l.syn and click OK. An OccurrenceTexture Parameters dialog appears.3.Deselect all of the Textures to Compute options except for Data Range.4.Set the Processing Window: Rows and Cols to 7 and 7.5.In the Enter Output Filename field, enter ndv_hh.tex and click OK.Create Color-coded Texture Map6.In the Available Bands List, click Display #1 and select New Display.7.Select Data Range under ndv_hh.tex and click Load Band.8.From the Display #2 menu bar, select Enhance→[Image] Square Root.9.From any Display group menu bar, select Tools→Link→ Link Displays. The Link Displays dialog appears.Click OK to link the original image (Display #1) with the colored texture image (Display #2).10.Click in an Image window to toggle between the two images.11.Double-click inside an Image window to display the Cursor Location/Value tool. Examine the data values in thetextured image, and compare these to the original image.12.From the Display #2 menu bar, select Tools→Color Mapping→Density Slice. A Density Slice Band Choicedialog appears.13.Select the Data Range band and click OK. A Density Slice dialog appears.14.Accept the default ranges by clicking Apply.15.Try creating your own density-sliced image and view the results.16.Keep Display #2 open for the next exercise.Image-Map OutputIn this exercise, you will create a map of your color-coded textured image and add a border and map key.1.From the Display #2 menu bar, select Overlay→Annotation. An Annotation dialog appears.2.From the Annotation dialog menu bar, select Options→Set Display Borders.3.In the Display Borders dialog, enter 100 in the upper field, and leave the remaining fields 0.4.Click Border Color and select Items 1:20→White. Click OK. This adds a 100-pixel white border at the top ofthe image.5.Move the Image box in the Scroll window to the top of the image containing the border.6.Enter a map title in the empty field in the Annotation dialog. Set the Size value to 16. Click the Color box onceto select black.7.Click in the Image window to show the map title, then move it inside the white border to the far left. Right-clickto lock the map title in place. You can place multiple text items on the image in this manner, and you can change their font size, type, color, and thickness as desired.8.From the Annotation dialog menu bar, select Object→Color Ramp.9.Enter Min and Max values of 0 and 255 respectively, set Inc to 4, and set the font Size to 14 to annotate thecolor ramp.10.Click in the Image window to show the map key, move it inside the white border to the far right, then right-clickto lock it in place. The following figure shows a sample map; your results may be different.11.Save the image to a PostScript file by selecting File→Save Image As→Postscript File from the Display #2menu bar. An Output Display to PostScript File dialog appears.12.Leave the default values, and enter an output filename or accept the default name of ndv_hh.ps. Click OK.Or, output the map directly to your printer by selecting File → Print from the Display #2 menu bar.13.When you are finished, select File→Exit from the ENVI main menu bar.。