Difference between revisions of "RadxDealias James and Houze 2001"

Latest revision as of 21:04, 7 May 2021

Doppler radars observe atmospheric motions by analyzing frequency shifts of pulses that backscatter off of hydrometeors moving toward or away from the radar. In order to reconstruct a wave, one needs to observe both the peaks and troughs of a wave (i.e., the sampling frequency needs to be double the frequency of the wave). Dealiasing (or folding) of Doppler velocity data occurs when the moving hydrometeors produce a frequency shift that is too high for a radar's pulse repetition frequency (PRF) to resolve. The phase shift measured by the radar may differ from the true phase shift by [math]\displaystyle{ \pm2n\pi }[/math], where n is an integer. The maximum unambiguous velocity a radar can measure, the Nyquist velocity ([math]\displaystyle{ V_n }[/math]), is defined as [math]\displaystyle{ V_n=\pm\frac{PRF\lambda}{4} }[/math]. Although there is no way of knowing whether velocity at a single gate is folded, users can use nearby observations and large radial discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one algorithm that automates the unfolding process for large radar datasets.

Dealiasing

Fundamentally, RadxDealias attempts to unfold velocity data by adjusting the velocity value at each gate by Nyquist intervals such that the velocity value at a gate is similar to values at gates that are nearby, both spatially and temporally. Working from higher to lower elevation angles, RadxDealias steps through several processes: quality control, deliasing using velocity data from a prior radar volume (or an environmental sounding, if no prior volume exists) as a reference, dealising using velocity data in the 8 neighboring gates from the same tilt as a reference, dealiasing using velocity data from nearby windows of 11x11 or 21x21 gates as a reference, and finally, comparing any remaining folded gates to the environmental sounding.

Quality Control

First, radar volumes are cleaned by discarding gates with reflectivity values outside a user-defined range (e.g., 0 to 80 dBZ, generally). Next, a Bergen and Albers filter can be applied to remove isolated gates.

Initial Dealiasing

Gates are first compared to the corresponding gates in the previous radar volume in time and in the previous (i.e., higher) elevation angle. All gates that can be dealiased to within [math]\displaystyle{ 0.25V_n }[/math] of those prior gates are modified and marked as "good" data. If a prior radar volume does not exist, the environmental sounding data is used instead. Note: environmental soundings tend to produce more errors than a previous radar volume.

Spatial Dealiasing

Remaining "bad" gates are next compared with the 8 neighboring gates within the same tilt. For this step, the adjusted velocity must agree with all neighboring "good" gates within a horizontal wind shear tolerance of [math]\displaystyle{ 0.4V_n }[/math]. This spatial dealiasing occurs twice, first scanning the tilt angle in the clockwise direction and then in the counterclockwise direction. On the third and final pass, the tolerance is increased to [math]\displaystyle{ V_n }[/math] and the adjusted velocity must agree with a majority of the neighboring "good" gates. All gates that are successfully changed in this step are marked as "good" data.

Window Dealiasing

Remaining "bad" gates are next compared with gates in neighboring azimuth-range windows. If the default 11x11 window has less than 5 gates of "good" data, a larger window of 21x21 is used. If that larger window still has too few "good" data values, the gate is saved for the final dealiasing step. All values that fall within [math]\displaystyle{ 0.4V_n }[/math] of the population mean are adjusted and marked as "good" gates, as long as the standard deviation is within a user-defined threshold.

Auxiliary Dealiasing

Finally, any remaining gates are attempted to be corrected with the environmental sounding data, but only if such data exists. The adjusted value is accepted only if it falls within [math]\displaystyle{ 0.5V_n }[/math] of the environmental wind value. A maximum of 10 passes that alternate scanning in the clockwise and counterclockwise directions are performed and any remaining gates that cannot be unfolded are deleted.

Sounding Quality

In the course of dealiasing radar data, it may become necessary to visualize to what degree the environmental sounding is influencing the final output. RadxDealias allows for the output sounding or Velocity Azimuth Display (VAD) velocities to be saved to the output files instead of the dealiased velocities by setting output_soundVol = TRUE. This option may be helpful in debugging.

References

James, C. N., & Houze , R. A., Jr. (2001). A Real-Time Four-Dimensional Doppler Dealiasing Scheme, Journal of Atmospheric and Oceanic Technology, 18(10), 1674-1683. Link

@@ Line 1: / Line 1: @@
 === '''Overview''' ===
-RadxDealias implements a four-dimensional scheme to correct velocity data when wind speeds exceed the Nyquist velocity. Based on [http://wiki.lrose.net/index.php/RadxDealias_James_and_Houze_2001#References James and Houze 2001], RadxDealias uses a combination of corrected velocity data from a prior radar file and environmental soundings to unfold velocity data.
+RadxDealias implements a four-dimensional scheme to correct folded velocity data. Based on [http://wiki.lrose.net/index.php/RadxDealias_James_and_Houze_2001#References James and Houze 2001], RadxDealias uses a combination of corrected velocity data from a previous radar volume close in time and environmental soundings to help unfold velocity data.
 === '''Background''' ===
-Doppler radars observe atmospheric motions by analyzing the frequency shift of the pulse after backscattering by hydrometeors moving towards or away from the radar. In order to reconstruct a wave, one needs to observe both the peaks and troughs of a wave (i.e., the frequency of sampling needs to be double the frequency of the wave). Dealiasing (or folding) of Doppler velocity data occurs when the moving hydrometeors produce a frequency shift too high for a radar's pulse repetition frequency (PRF) to resolve. The phase shift measured by the radar may differ from the true phase shift by <math>\pm2n\pi</math>, where n is an integer. The maximum unambiguous velocity a radar can measure, the Nyquist velocity (<math>V_n</math>), is defined as <math>V_n=\pm\frac{PRF\lambda}{4}</math>. Although there is no way of knowing whether velocity at a gate is folded on its own, users can use nearby observations and discontinuities in velocity that jump from <math>+V_n</math> to <math>-V_n</math> to diagnose regions where unfolding has occurred. RadxDealias is one of several algorithms that can automate that process on large amounts of data.
+Doppler radars observe atmospheric motions by analyzing frequency shifts of pulses that backscatter off of hydrometeors moving toward or away from the radar. In order to reconstruct a wave, one needs to observe both the peaks and troughs of a wave (i.e., the sampling frequency needs to be double the frequency of the wave). Dealiasing (or folding) of Doppler velocity data occurs when the moving hydrometeors produce a frequency shift that is too high for a radar's pulse repetition frequency (PRF) to resolve. The phase shift measured by the radar may differ from the true phase shift by <math>\pm2n\pi</math>, where n is an integer. The maximum unambiguous velocity a radar can measure, the Nyquist velocity (<math>V_n</math>), is defined as <math>V_n=\pm\frac{PRF\lambda}{4}</math>. Although there is no way of knowing whether velocity at a single gate is folded, users can use nearby observations and large radial discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one algorithm that automates the unfolding process for large radar datasets.
@@ Line 11: / Line 11: @@
 ===== '''Quality Control''' =====
-First, volumes can be prepped by discarding gates with reflectivity values outside a user-defined range (e.g., 0 to 80 dBZ, generally). Next, a Bergen and Albers filter can be applied to remove isolated gates.
+First, radar volumes are cleaned by discarding gates with reflectivity values outside a user-defined range (e.g., 0 to 80 dBZ, generally). Next, a Bergen and Albers filter can be applied to remove isolated gates.
 ===== '''Initial Dealiasing''' =====
-Gates are first compared to the corresponding gate in the previous radar volume in time and the gate in the previous (i.e., higher) elevation angle. All gates that can be dealiased to within a factor of 0.25 of the Nyquist velocity of those prior gates are modified and marked as "good" data. If a prior radar volume does not exist, the environmental sounding data is used instead. Note: environmental soundings tend to produce more errors than a previous radar volume in time.
+Gates are first compared to the corresponding gates in the previous radar volume in time and in the previous (i.e., higher) elevation angle. All gates that can be dealiased to within <math>0.25V_n</math> of those prior gates are modified and marked as "good" data. If a prior radar volume does not exist, the environmental sounding data is used instead. Note: environmental soundings tend to produce more errors than a previous radar volume.
 ===== '''Spatial Dealiasing''' =====
-Remaining "bad" gates are next compared with the 8 neighboring gates within the same tilt. For this step, the allowance in terms of horizontal wind shear is a factor of 0.4 of the Nyquist velocity and must agree with all 8 neighboring gates. This spatial dealiasing occurs twice, first scanning the tilt angle in a clockwise fashion and then counterclockwise. On the third and final pass, the threshold is increased to the Nyquist velocity and changed for any gates that agree with a majority of the 8 neighboring gates. All gates that are successfully changed are marked as "good" data.
+Remaining "bad" gates are next compared with the 8 neighboring gates within the same tilt. For this step, the adjusted velocity must agree with all neighboring "good" gates within a horizontal wind shear tolerance of <math>0.4V_n</math>. This spatial dealiasing occurs twice, first scanning the tilt angle in the clockwise direction and then in the counterclockwise direction. On the third and final pass, the tolerance is increased to <math>V_n</math> and the adjusted velocity must agree with a majority of the neighboring "good" gates. All gates that are successfully changed in this step are marked as "good" data.
 ===== '''Window Dealiasing''' =====
-Remaining "bad" gates are next compared with gates in neighboring azimuth-range windows. If the default 11x11 window has less than 5 gates of good data, a larger window of 21x21 is used. If that larger window still has too few good data values, the gate is saved for the final dealiasing step. The adjustment threshold for this step is the Nyquist velocity. All values that fall within this threshold of the population mean are changed, as long as the standard deviation is within a user-defined threshold.
+Remaining "bad" gates are next compared with gates in neighboring azimuth-range windows. If the default 11x11 window has less than 5 gates of "good" data, a larger window of 21x21 is used. If that larger window still has too few "good" data values, the gate is saved for the final dealiasing step. All values that fall within <math>0.4V_n</math> of the population mean are adjusted and marked as "good" gates, as long as the standard deviation is within a user-defined threshold.
 ===== '''Auxiliary Dealiasing''' =====
-Finally, any remaining gates are attempted to be corrected with the environmental sounding data, but only if such data exists. The adjusted value is saved only if it falls within a factor of 0.5 of the Nyquist velocity of the environmental wind value. A maximum of 10 passes that alternate scanning in the clockwise and counterclockwise directions are performed and any remaining gates that cannot be unfolded are deleted.
+Finally, any remaining gates are attempted to be corrected with the environmental sounding data, but only if such data exists. The adjusted value is accepted only if it falls within <math>0.5V_n</math> of the environmental wind value. A maximum of 10 passes that alternate scanning in the clockwise and counterclockwise directions are performed and any remaining gates that cannot be unfolded are deleted.