View source for RadxDealias James and Houze 2001

=== '''Overview''' ===
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.


=== '''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, 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.

===== '''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.

===== '''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.

===== '''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.

===== '''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.


=== '''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. [https://doi.org/10.1175/1520-0426(2001)018%3C1674:ARTFDD%3E2.0.CO;2 Link]