RadxDealias James and Houze 2001 - Revision history

Jcdehart at 21:04, 7 May 2021

2021-05-07T21:04:53Z

Jcdehart at 21:03, 7 May 2021

2021-05-07T21:03:45Z

Jcdehart at 21:02, 7 May 2021

2021-05-07T21:02:58Z

Jcdehart at 21:01, 7 May 2021

2021-05-07T21:01:17Z

Jcdehart at 20:38, 7 May 2021

2021-05-07T20:38:28Z

Jcdehart at 20:27, 7 May 2021

2021-05-07T20:27:17Z

Jcdehart: /* Overview */

2021-05-07T20:27:10Z

Overview

Jcdehart at 20:19, 7 May 2021

2021-05-07T20:19:28Z

Jcdehart at 18:51, 7 May 2021

2021-05-07T18:51:47Z

Jcdehart: /* Dealiasing */

2021-05-07T18:45:20Z

Dealiasing

← Older revision		Revision as of 21:04, 7 May 2021
Line 4:		Line 4:

	=== '''Background''' ===		=== '''Background''' ===
−	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 discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one ~~of several algorithms~~ that ~~automate~~ the unfolding process ~~to be used on~~ large radar datasets.	+	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.

← Older revision		Revision as of 21:03, 7 May 2021
Line 4:		Line 4:

	=== '''Background''' ===		=== '''Background''' ===
−	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 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 discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one of several algorithms that automate the unfolding process to be used on large radar datasets.	+	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 discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one of several algorithms that automate the unfolding process to be used on large radar datasets.

← Older revision		Revision as of 21:02, 7 May 2021
Line 4:		Line 4:

	=== '''Background''' ===		=== '''Background''' ===
−	Doppler radars observe atmospheric motions by analyzing frequency shifts of pulses that backscatter off of 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 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 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 discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one of several algorithms that automate the unfolding process to be used on large radar datasets.	+	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 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 discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one of several algorithms that automate the unfolding process to be used on large radar datasets.

@@ Line 17: / Line 17: @@
 ===== '''Spatial Dealiasing''' =====
-Remaining "bad" gates are next compared with the 8 neighboring gates within the same tilt. For this step, the tolerance in terms of horizontal wind shear is <math>0.4V_n</math> 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 <math>0.4V_n</math>. 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 <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.
+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.

@@ Line 4: / Line 4: @@
 === '''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 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 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 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 discontinuities in velocity to diagnose regions where aliasing has occurred. RadxDealias is one of several algorithms that automate the unfolding process to be used on 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 tolerance in terms of horizontal wind shear is <math>0.4V_n</math> 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.
+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 <math>0.4V_n</math>. 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.
+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 <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.

@@ Line 1: / Line 1: @@
 === '''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''' ===

← Older revision		Revision as of 20:19, 7 May 2021
Line 4:		Line 4:

	=== '''Background''' ===		=== '''Background''' ===
−	Dealiasing (or folding) of Doppler velocity data occurs when moving hydrometeors produce a frequency shift too high for a radar's pulse repetition frequency (PRF) to resolve. The ~~Doppler~~ phase shift measured by the radar may differ from the true phase shift by <math>\pm2n\pi</math>, where n is an integer.	+	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.

@@ Line 8: / Line 8: @@
 === '''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 nearby spatial and temporal gates. 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.
+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.
 === '''Sounding Quality''' ===