HAZEREM

Removes haze from satellite imagery


EnvironmentsPYTHON :: EASI :: MODELER
Batch ModeYes
Quick linksDescription :: Parameters :: Parameter descriptions :: Details :: References :: Related

Back to top

Description


HAZEREM removes haze from optical imagery using a variation of the Haze Optimized Transformation method. By default, the function uses the metadata in the file to determine the sensor and band positions. If the metadata is not available, the user must provide this information by specifying the Satellite Sensor (ASENSOR), Blue, Red, Near-Infrared Channels (VISIRCH), and Channel Option (CHANOPT) parameters. HAZEREM is typically run after cloud and water masks are generated using MASKING. The output is always TOA Radiance. The input can be either TOA Radiance or TOA Reflecance.
Back to top

Parameters


Name Type Length Value range
Input: Input multispectral image file* Raster port 1 -    
InputPAN: Input panchromatic image Raster port 0 -    
Source background options String 0 -   Default: File Metadata
Source background values Integer 0 - 1024 Default: 0
Satellite sensor String 0 -  
Blue, red, near-infrared channels Integer 0 - 3  
Channel options String 0 -    
InputMask: Input haze, cloud, and water mask file Bitmap port 0 -    
Haze cover (%) Integer 0 - 2  
Haze filter size (pixels) Integer 0 - 1 1 - 301
Default: 7
Output: Output haze-free multispectral image file* Raster port 1 -    
OutputPAN: Output haze-free panchromatic image file Raster port 0 -    

* Required parameter
Back to top

Parameter descriptions

Input: Input multispectral image file

Specifies the name of the input file that contains multispectral imagery from one of the supported sensors (see ASENSOR).

The pixel values in the input raster file (FILI) must be TOA Radiance or TOA Reflectance.

InputPAN: Input panchromatic image

Optionally specifies the name of the input file that contains the panchromatic image associated with the input multispectral image.

The panchromatic image must have been acquired at the exact same time, under the same haze conditions as the multispectral image. The input panchromatic image must be a single-channel file. If the specified file contains more than one channel, only the first channel is processed.

Source background options

Specifies, potentially with Source Background Values, which pixels in the source image are to be considered background (NoData) pixels. In general, if a pixel is considered NoData, the application handles the pixel in a special manner.

Available options are:

See the Source Background Values parameter for specific examples.

Source background values

Optionally specifies the source background value(s) when the following options are specified for the Source Background Options parameter:
Note: Specify multiple values in a comma-delimited list. The first value is applied to the first channel, the second value to the second channel, and so on. If fewer values are specified than the number of input channels, the last value is repeated for all remaining channels. If more values are specified than the number of input channels, the extra values are ignored.
For example:

Satellite sensor

Optionally specifies the name of the sensor used to capture the input image.

See the Function Details section for a list of supported sensors.

If this parameter is not specified, the function checks for the PlatformName metadata tag at the file level.

If this parameter is not specified and the PlatformName metadata tag is not found, the function will error.

Blue, red, near-infrared channels

Optionally specifies the channels that contain the blue, red, and near-infraRed (NIR) bands. If this parameter is not specified, channel numbers are automatically determined from the specified sensor type (ASENSOR). If the NIR band is not available, the algorithm will still run, but the urban areas and bare soil pixels will likely be overcorrected in the output.

Note: If the channel order has been altered from its standard order as supplied by the data vendor, this parameter must be used to define the new order.

Channel options

Optionally specifies how the channels in the file should be processed and written to the output file. Channel options are defined using a string of characters, one character for each channel.

Supported processing options are:

For example, if the input file has four channels and this parameter is either 'p,p,c,s' or 'p,p,c', the first two channels should have haze removal applied, the third should be copied as is, since the input is TOA Radiance and the fourth skipped. The output file would then contain three channels: the first two with haze removed and the third with an exact copy of the third input channel.

If this parameter is not specified, the function uses a default value for the specified sensor. The default applies haze removal to all visible channels (Coastal/Blue/Green/Red/RedEdge/NIR). and copies the rest (SWIR) and thermal IR channels.

Output is always in TOA Radiance

InputMask: Input haze, cloud, and water mask file

Optionally specifies the name of the file that contains haze, cloud, and water masks. These masks are typically generated by MASKING, which is run prior to running HAZEREM.

Supplying haze, cloud, and/or water masks improves haze removal results over land. Haze will still be removed over water, but clouds are left undisturbed. When a cloud mask is not provided but the scene contains clouds, the cloud pixels are assumed to be haze and will receive the maximum possible correction, resulting in an overcorrected image.

Haze cover (%)

Optionally specifies the percentage of the image that is covered in haze. The default value is 50.

This value should be set as low as possible while still allowing haze to be removed. Low percentages cause the least change to the original image values. It may be necessary to run HAZEREM with different haze cover values to achieve the best results.

Haze filter size (pixels)

Optionally specifies the size of a smoothing filter to apply to the haze computations. Valid values are odd integers, from 1 to 301; the default value is 7.

Although large filter sizes tend to generate more natural color, if the haze thickness varies quickly (that is, it is thin, high cloud or fills small valleys), a large filter size will leave small pockets and fringes of haze. This value should be set as high as possible while still allowing all haze to be removed. It may be necessary to run HAZEREM with different haze smoothing filter values to achieve the best results.

Note: The smoothing filter used in internal haze computations is not applied to the actual imagery.

Output: Output haze-free multispectral image file

Specifies the name of the output file to receive the haze-corrected image. The output will be in TOA Radiance.

The specified file must not already exist; a new file will always be created to store the haze-corrected multispectral image.

OutputPAN: Output haze-free panchromatic image file

Optionally specifies the name of the output file to receive the haze-corrected panchromatic image if an input PAN file was specified. The output will be in TOA Radiance. (FILI_PAN).

The specified file must not already exist; a new file will always be created to store the haze-corrected panchromatic image.

Back to top

Details

HAZEREM is the second step in a typical atmospheric correction workflow. This algorithm prepares the image by correcting the radiometric variability caused by haze, thus reducing some of the atmospheric effect on the image in the raw domain. This produces a better-looking image, ready for analysis or operations such as automated GCP collection, color balancing, or mosaicking, without relying on the accuracy of atmospheric lookup tables.

HAZEREM uses the coarse classification masks created by MASKING to identify and derive haze statistics. The identification of haze is based on a transformation of the visible band space, where spectral response to diverse surface cover classes under clear-sky conditions is highly correlated. Bright land surface cover types such as urban and bare soil may be problematic because they are often confused with haze. The algorithm produces the best results when there are a good number of clear and hazy pixels over vegetation.

The output will be in TOA Radiance.

Supported sensors

The following sensors are supported:
  • ALI
  • ALOS Avnir-2
  • Aster
  • Cartosat PAN
  • CBERS-4
  • Deimos-1
  • Deimos-2
  • DMC
  • FASat-Charlie
  • Formosat-2
  • FORMOSAT-5
  • Geoeye-1
  • GF1
  • GF2
  • GF4
  • GF6
  • GF7
  • Gokturk1
  • IRS-1A
  • IRS-1B
  • IRS-1C
  • IRS-1D
  • IRS-P6
  • Ikonos-2
  • Jilin-1
  • KazEOSat-2
  • KOMPSAT-2
  • KOMPSAT-3
  • KOMPSAT-3A
  • Landsat-4 MSS
  • Landsat-5 MSS
  • Landsat-4 TM
  • Landsat-5 TM
  • Landsat-7 ETM+
  • Landsat-8
  • OrbView-3
  • PeruSAT-1
  • PlanetScope
  • Pleiades
  • QuickBird
  • PlanetScope
  • RapidEye
  • Resourcesat-2
  • SAC-C
  • Sentinel-2
  • SPOT-1
  • SPOT-2
  • SPOT-3
  • SPOT-4
  • SPOT-5
  • SPOT-6
  • SPOT-7
  • SuperView
  • Thaichote (THEOS)
  • TripleSat
  • WorldView-2
  • WorldView-3
  • WorldView-4
  • ZY3
  • ZY3-2
Note: When specifying a sensor, use the exact syntax, as shown in the preceding list.
Back to top

References

Y. Zhang, B. Guindon, and J. Cihlar. "An image transform to characterize and compensate for spatial variations in thin cloud contamination of Landsat images", Remote Sensing of Environment 82 (2002), pp. 173-87.

© PCI Geomatics Enterprises, Inc.®, 2024. All rights reserved.