Research Article | Open Access
Yidong Hou, Biyang Wen, Yonghuai Yang, Jing Yang, Caijun Wang, "Two-Dimensional River Flow Patterns Observed with a Pair of UHF Radar System", International Journal of Antennas and Propagation, vol. 2017, Article ID 4792324, 11 pages, 2017. https://doi.org/10.1155/2017/4792324
Two-Dimensional River Flow Patterns Observed with a Pair of UHF Radar System
A pair of ultrahigh-frequency (UHF) radars system for measuring the two-dimensional river flow patterns is presented. The system consists of two all-digital UHF radars with exactly the same hardware structure, operating separately at 329–339 MHz and 341–351 MHz. The adoption of direct radio frequency (RF) sampling technique and digital pulse compression simplifies the structure of radar system and eliminates the distortion introduced by the analog mixer, which improves the SNR and dynamic range of the radar. The field experiment was conducted at Hanjiang River, Hubei province, China. Over a period of several weeks, the radar-derived surface velocity has been very highly correlated with the measurements of EKZ-I, with a correlation coefficient of 0.958 and a mean square error of 0.084 m/s.
Monitoring the river flow velocity serves as an important part in the field of flood control, waterborne freight, fisheries, and so on. Conventionally, it is measured using the method of river buoy, rotor current meter, ADCP, and electromagnetic flow meter. However, these instruments have to contact with the water directly, bringing much pressure on their placement and maintenance, and require more professional and high technological operators. Besides, these methods are of single-point measurement, which hardly get the flow field of the whole waters. All these challenges have promoted the development of noncontact measurement. In the last decades, radar, an important noncontact measurement method, has been widely used in the measurement of ocean current and river flow. The development of high frequency (HF) and UHF radars has mainly undergone three stages: theoretical basis, monitoring ocean current with HF radar, and river flow measurement using UHF radar.
Crombie  proposed the backscattering mechanism of ocean waves in 1955, and then in 1966, Wait  discovered and quantitatively explained that this backscattering mechanism is produced by the Bragg scattering between electromagnetic waves and half wavelength ocean waves. In 1972, the first- and second-order radar cross-section (RCS) equations were derived by Barrick [3, 4], which laid a theoretical foundation for monitoring ocean and river surface dynamic parameters using radar. After that, Lipa and Barrick [5–7] pioneered efforts to derive ocean information from the HF spectra. Following their lead, HF radar has been widely studied with fruitful results. Several representative HF radar systems including SeaSonde, Wellen Radar (WERA), Ocean Surface Current Radar (OSCR), Ground Wave Radar (GWR), and Ocean State Monitoring and Analysis Radar (OSMAR) have come out one after another. The SeaSonde system using a crossed-loop/monopole as receiving antenna array not only has outstanding performance, but also has the advantages of compact structure, small size, and low cost [8, 9]. Compared with the SeaSonde, the WERA system in Germany uses an array receiving antenna, which has more flexible spatial resolution . The OSCR system in the UK is evolved from the SeaSonde system, which is mainly used to measure the ocean current with high accuracy in the near distance . The GWR system developed in Canada in 1990, using the peak power of 16 KW, is a large array radar which has been not only used to monitor the ocean state, but also mainly used to detect ships, icebergs, and so on . The OSMAR system was developed by Wuhan University in 1987, and so far there are three versions which are array radar (OSMAR2000), portable radar (OSMAR-S), and dual-frequency radar (OSMAR-SD), respectively [13, 14].
The success of HF radar in monitoring the ocean state has promoted the development of the measurement of river surface dynamic parameters with UHF radar. Randell et al.  measured the ocean clutter characteristics using coherent Doppler UHF radar in 2000. Based on the SeaSonde system, CODAR developed the UHF radar system (RiverSonde) to measure the river surface flow and has carried out several field experiments to verify its high performance [16–19]. Wuhan University has been devoted to the research of noncontact measurement for river flow monitoring since 2004 and successfully designed its first UHF radar based on the OSMAR system in 2005 . After years of development, Wuhan University finally designed its first all-digital UHF radar system USDPR in 2013, along with a lot of experiments that demonstrated its excellent performance [21–23].
OSMAR-SU is the latest generation of all-digital UHF radar developed by Wuhan University, whose structure is more simplified and the radar performance gets further improved with the application of direct RF sampling, digital pulse compression, and transmit-receive sharing technology. In this paper, we will describe the structure of OSMAR-SU and the realization of all modules in detail in Section 2. The radar waveform parameters and signal processing methods are presented in Section 3. The field experiment and the comparisons of results are described in Sections 4 and 5, respectively. Conclusion is presented in Section 6.
2. The Radar Structure and Realization of All Modules
Figure 1 shows the structure of OSMAR-SU. There are five sections, antenna system, receiver motherboard, Transmit-Receive switch (T/R switch), small power transmitter, and Industrial Personal Computer (IPC). In this radar system, direct RF sampling is utilized to realize the pulse compression in digital domain, so the structure is greatly simplified.
2.1. Receiver Motherboard
OSMAR-SU employs a more compact structure, as Figure 1 shows. Almost all circuit modules, including analog front end (AFE), analog-to-digital converter (ADC), clock generator, signal generator, serial port, USB module, and FPGA, are integrated on the receiver motherboard. Thereinto, the architecture of AFE and ADC is shown in Figure 2. In the all-digital radar, the AFE only functions as amplifier and filter and no longer contains the mixer, which means that the distortion of nonlinear and intermodulation introduced by the mixer is completely eliminated, and the dynamic range of the receiver will be greatly improved. The linear gain and third-order intermodulation distortion (IMD3) of the AFE are measured as shown in Figure 3, and the measurement results indicate that the gain of the AFE is 40.7 dB, the 1 dB compression point (P1dB) is −21 dBm, the input 3rd-order intercept point (IIP3) is −2 dBm, the sensitivity of AFE is −159 dBm/Hz, and the spurious free dynamic range (SFDR) is up to 108 dB. In addition, the channel isolation and bandwidth of the AFE are 65 dB and 35 MHz, respectively.
Structure of clock and signal generator module is as Figure 4 shows. Temperature compensated crystal oscillator (TCXO) as the clock source of the system output 10 MHz. Two coherent clock signals are generated by digital signal phase-lock loop (DSPLL), and their frequencies are 81.92 MHz and 983.04 MHz. The clock (983.04 MHz) is used as a reference clock for signal generator, while the other clock (81.92 MHz) is driven to 4 channels and is employed as the system clock for FPGA and the sampling clock for ADC separately. The signal generator consisting of direct digital synthesizer (DDS), RF filter, amplifier, and transformer is used to generate radar signal.
2.2. Antenna System and Transmitter
OSMAR-SU antenna system adopts uniform liner array (ULA). The array element spacing equals half wavelength, and the array consists of three five-element Yagi antennas. The beam width of the Yagi antenna is 90 degrees, and the detection range is between −70 and 70 degrees. The antenna, located in the middle of the array, is used as a common antenna for receiving and transmitting, and the radio frequency (RF) switch switches the T/R functions. The transmitter power of the OSMAR-SU is only 250 mW. Table 1 lists more detailed parameters of the antenna and transmitter.
2.3. T/R Switch
T/R switch architecture is shown as Figure 5. Single RF switch difficultly satisfies the requirement of isolation; therefore two RF switches are connected in series to realize a relatively high isolation. Transmitting and receiving status shifting of T/R switch are controlled by a pair of complementary signals of transmitting pulse (TP) signal and bating ground wave (BGW) (defined as receiving pulse) signal. The T/R switch achieves a high isolation of 66 dB, a fast switching time of 30 ns, and a high input power of 40 dBm.
3. Waveform Parameters and Signal Processing Algorithm
3.1. Waveform Parameters Design
OSMAR-SU adopts Frequency Modulated Interrupting Continuous Wave (FMICW) mechanism, and the waveform parameters can be configured according to the requirements of the experimental environment. In order to avoid the mutual interference between the two OSMAR-SUs, the central frequencies of the two radars are 334 MHz and 346 MHz. The signal bandwidth is 10 MHz; that is to say, two radars operate at 329–339 MHz and 341–351 MHz, respectively. All the other waveform parameters are consistent, and detailed parameters are shown in Table 2. Based on Barrick’s theory , it is deduced that 334 MHz and 346 MHz correspond to the water wave velocity of 0.837 m/s and 0.822 m/s; therefore the maximum flow velocities that can be measured by the two radars are, respectively, 4.652 m/s and 4.446 m/s. There is no need to consider the problem of blind area, since the radars are located close to the river bank. In a pulse interrupt cycle, the receiving cycle starts immediately after the end of the transmitting cycle.
3.2. Signal Processing Algorithm
Signal processing flow of OSMAR-SU is shown in Figure 6. The range spectrum is obtained from echo signal through bandpass sampling, deramp (digital mixing), cascaded integrator-comb (CIC) filter, finite impulse response (FIR) filter, and fast Fourier transform (FFT), which is an equivalent pulse compression process. After digital pulse compression and FFT, Doppler spectrum is obtained, which include first-order Bragg peaks and noise. The first-order peaks are separated from noise by absorbance subtraction method; then direction of arrival (DOA) for every Doppler frequency-point in the first-order peaks is estimated through MUSIC algorithm . Next, utilizing the frequency offset of Doppler spectral peaks with Bragg frequency, the actual flow velocity of a certain point in the river is calculated. After interpolating and smoothing, radial flow field is obtained from rough results. Finally, two radial flow fields are applied to synthesize two-dimensional flow pattern. Before DOA estimations, it is necessary to calibrate the array. This paper adopts active calibration method, placing an emission source in front of the array to calibrate its amplitude and phase.
The Doppler spectra are broadened by the radial flow of different directions to varying degrees. The relationship between the maximum Doppler shift and the maximum radial flow can be expressed as where is the maximum Doppler shift, is the maximum radial river flow, is the wavelength of radar signal, and is the Bragg shift. There are three cases in the actual DOAs and radial flow estimation, which are
For case 1, as shown in Figure 7(a), there is no aliasing between the positive and negative Bragg echoes. The first-order Doppler-broadened peaks are selected in and as shown as the dark areas in the Figure 7(a), and then the DOAs of Doppler shifts can be estimated by the classic MUSIC algorithm, and the radial flow can be computed by (1).
When the maximum radial flow is more than 0.82 m/s (Bragg-wave velocity) and less than 1.64 m/s, the positive and negative Bragg echoes overlap, but the aliasing areas do not exceed . We call case 2 partial aliasing, as shown in Figure 7(b). The intervals and , as shown as the dark areas in the Figure 7(b), are selected as the first-order Doppler-broadened peaks, because these two intervals completely contain the flow toward the radar and the flow away from the radar, respectively.
From case 1 and case 2, we can see that the classic MUSIC algorithm can unambiguously estimate the DOAs under the condition that the maximum radial flow is less than 1.64 m/s; nevertheless, the aliasing areas between the positive and negative Bragg echoes in case 3 will exceed , as shown in Figure 7(c). The same intervals as case 2 are still selected as the first-order Doppler-broadened peaks. Obviously, there are aliasing areas in the selected intervals. Since the echoes in the aliasing areas are coherent, the classic MUSIC algorithm cannot unambiguously resolve the direction of the velocity. In view of this situation, the five backup receiving channels reserved in OSMAR-SU will come in handy. The eight-element linear array and the spatial smoothing technique [25, 26] will be used to unambiguously estimate the DOAs of Doppler shifts in the aliasing area.
Here we give the simulation results of case 3 as shown in Figure 8. We assume that the relationship between radial flow and direction of arrival (DOA) is described by where is the radial river flow, is the mean river flow, and is the DOA. Let m/s ( belongs to the case 3) and the simulation of the Doppler spectrum is shown in Figure 8(a), where the red dotted lines represent Bragg shifts, the gray shaded area represents the aliasing region, and the black box areas are selected for compute. We can see that the Doppler spectrum is severely broadened. We give the spatial spectra of three typical Doppler frequencies, as shown in Figures 8(b), 8(c), and 8(d). Thereinto, Figure 8(b) is in the outside of the aliasing region, while Figures 8(c) and 8(d) are in the aliasing region. The simulation results indicate that (1)The results of the classical MUSIC algorithm and those of the space-smoothing MUSIC algorithm are consistent if there is no aliasing.(2)When aliasing occurs, the classical MUSIC algorithm will be invalid, while space-smoothing MUSIC still performs well.(3)In the aliasing area, a relatively large DOA (absolute value) corresponds to a relatively large radial velocity (absolute value).
(a) The simulation of the Doppler spectrum
(b) MUSIC spectra at −4.84 Hz
(c) MUSIC spectra at −1.85 Hz
(d) MUSIC spectra at 1.04 Hz
In 2016, from July to August, a field experiment is implemented in the lower reaches of Hanjiang River. As the largest tributary of the Yangtze River, the width of the reach where the radars are located is about 300 meters. As is shown in Figure 9, the river flows from west to east; the radar field of view of Site A covers the river bend, while Site B monitors the straight channel, and the expected flow direction is shown as the yellow line in Figure 9. During the experiment, the radars successively suffered from concentrative rainfall and continuous high temperature weather but still received high quality data. The preliminary experiments were conducted in flood period, and water level was high in the lower reaches of the Yangtze River, resulting in a relatively small flow velocity; the maximum surface flow velocity was about 0.7 m/s, while in the late experiment the water level dropped, and the maximum surface flow velocity climbed up to 1.2 m/s.
Figure 10 shows where and how the radars are deployed. The radar in Figure 10(a) is fixed in the yard of the hydrologic station close to the river bend with the antenna about 6 m above the water level. The other radar is fixed on the riverbank near the straight channel with the antenna about 5 m above the water level. Both radars are located on the same side of the river, about 155 meters apart. The radars look angle covers −70 to 70 degrees; to prevent the blind area between two sites, the normal of dual antenna array is shifted 20 degrees to the center. As is shown in Figure 9, the red and white sectors mark the field of view of the radars. The flow vector synthesis will be performed in the overlapping area. To compare, the flow measuring equipment EKZ-I records its observations. EKZ-I consists of cableway, lead fish and rotor flow meter. In Figure 9, the black dotted line marks where the cableway is, and the lead fish can cross the river along the cableway, so that it gets flow velocity at different distance from the offshore side.
5. Results and Comparisons
Figure 11 shows the Doppler-Range spectrum and the Doppler spectrum at 200 m, and the SNR is about 40 dB. Based on the first-order scattering equation of rough surface derived by Hickey et al. , the echo signal will be greatly enhanced when water wave and electromagnetic wave satisfy the resonance Bragg scattering. In this paper, the OSMAR-SU operates at 340 MHz, and corresponding Bragg resonance frequency is 1.88 Hz. Due to the existence of water flow, the Doppler spectrum is broadened; the larger the flow velocity is, the more obviously it broadens.
(a) Range-Doppler spectrum
(b) Doppler spectrum
Figure 12 shows the radial flow field of the two sites measured during the same time period. Figure 12(a) shows the result of Site A in the west, located at river bend where the flow direction is roughly from west to east. Due to the inertia of the river flow, the flow direction in this area is not perfectly along the river. Therefore, the radial flow velocity in front of the radar is not equal to 0 m/s. Figure 12(b) shows the result of Site B in the east, where the radar field of view is within the straight channel. The flow direction is from the left to the right of the radar along the river. The radial flow hence faces the radar on the left and the opposite on the right. The radial flow velocity in front of the radar is approximately equal to 0 m/s.
(a) Site A
(b) Site B
The radial flow field of the dual sites in Figure 12 is synthesized into two-dimensional flow pattern. As shown in Figure 13, the arrow direction represents the flow direction, and the arrow length represents the flow velocity. The flow direction obviously changes along the river. The closer to the shore, the smaller the flow velocity, which is in line with the actual situation. Flow velocities measured by EKZ-I and OSMAR-SU are compared as Figure 14 shows. The results of OSMAR-SU and EKZ-I have been highly correlated, with a correlation coefficient of 0.950 and a mean square error of 0.092 m/s. Also the result of dual sites system is better than single site system which has better accuracy in the straight channel (correlation coefficient 0.913, mean square error 0.102 m/s) than in the river bend (correlation coefficient 0.879, mean square error 0.330 m/s). This indicates the measurement accuracy of the dual sites system is better than single site system in the river bend or more complex area.
From 2016-7-15 to 2016-7-30, nine comparison experiments between EKZ-I and OSMAR-SU are conducted and the results are shown in Figure 15, with a correlation coefficient of 0.958, an average error of 0.033 m/s, and a mean square error of 0.084 m/s. Figure 16 is the two-dimensional flow pattern measured with previous generation of UHF radar (USDPR) designed by Wuhan University at 2014-8-13 in the same location. One radar antenna array was split into two units, 57 m apart, using transmission power of 10 W. More details of the experiment can be found in . It can be inferred from Figure 16 that single site measurement results (correlation coefficient 0.952, mean square error 0.091 m/s) are more accurate than the results of the dual sites (correlation coefficient 0.910, mean square error 0.158 m/s), and the detection range of USDPR is much smaller than OSMAR-SU. Table 3 lists some detailed comparisons of OSMAR-SU and USDPR, and the comparison results indicate that OSMAR-SU achieved a more accurate measurement with a lower transmitting power.
(a) Flow velocity comparisons
(b) Two-dimensional flow patterns
This paper presents a pair of all-digital UHF radars system (OSMAR-SU) for measuring the two-dimensional river flow patterns and describes in detail the hardware architecture, waveform parameters, and signal processing method of the radar. The application of direct RF sampling, digital pulse compression, and transmit-receive sharing technology greatly simplifies the structure of the radar and obviously improves dynamic range, sensitivity, and channel isolation of the receiver. In order to verify the performance of the OSMAR-SU, a two-site field experiment was conducted at Hanjiang River, Hubei province, China, in July 2016, and the measurement results of the radar are compared with those of EKZ-I. The correlation is as high as 0.958, and the mean square error is less than 0.084 m/s.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was supported by National Natural Science Foundation of China (Grant no. 61371198) and Natural Science Foundation of Hubei Province (Grant no. 2015CFB355). The authors would like to thank the staff of Hubei Xiantao hydrology station for helping them and providing high quality data.
- D. D. Crombie, “Doppler spectrum of sea echo at 13.56 Mc./s.,” Nature, vol. 175, pp. 681–682, 1955.
- J. R. Wait, “Theories of Prolate Spheroidal Antennas,” Radio Science, vol. 1, no. 4, pp. 475–512, 1966.
- D. E. Barrick, “First-order theory and analysis of MF/HF/VHF scatter from the sea,” IEEE Transactions on Antennas and Propagation, vol. 20, no. 1, pp. 2–10, 1972.
- D. Barrick, “Remote sensing of sea state by radar,” in Proceedings of the IEEE International Conference on Engineering in the Ocean Environment (Ocean '72), pp. 186–192, Newport, RI, USA, September 1972.
- D. E. Barrick, “Extraction of wave parameters from measured hf radar sea-echo doppler spectra,” Radio Science, vol. 12, no. 3, pp. 415–424, 1977.
- B. J. Lipa and D. E. Barrick, “Least-Squares Methods for the Extraction of Surface Currents from CODAR Crossed-Loop Data: Application at ARSLOE,” IEEE Journal of Oceanic Engineering, vol. 8, no. 4, pp. 226–253, 1983.
- B. J. Lipa and D. E. Barrick, “Extraction of sea state from HF radar sea echo: mathematical theory and modeling,” Radio Science, vol. 21, no. 1, pp. 81–100, 1986.
- D. Barrick, P. Lilleboe, B. Lipa, and L. Pederson, “Coastal surface-current mapping to 200 km with Long-Range SeaSondes,” in Proceedings of the MTS/IEEE Conference and Exhibition (OCEANS '00), Conference Proceedings (Cat. No.00CH37158), vol. 1, pp. 555–558, Providence, RI, USA, 2000.
- D. Barrick and B. Lipa, “Using antenna patterns to improve the quality of SeaSonde HF radar surface current maps,” in Proceedings of the IEEE 6th Working Conference on Current Measurement, (Cat. No.99CH36331), pp. 5–8, San Diego, Calif, USA, March 1999.
- T. Helzel, M. Kniephoff, and L. Petersen, “WERA: Remote ocean sensing for current, wave and wind direction,” in Proceedings of the IEEE US/EU Baltic International Symposium (BALTIC '06), pp. 1–8, Klaipėda, Lithuania, May 2006.
- T. M. Hammond, C. B. Pattiaratchi, D. Eccles, M. J. Osborne, L. A. Nash, and M. B. Collins, “Ocean surface current radar (OSCR) vector measurements on the inner continental shelf,” Continental Shelf Research, vol. 7, no. 4, pp. 411–431, 1987.
- K. J. Hickey, E. W. Gill, J. A. Helbig, and J. Walsh, “Measurement of ocean surface currents using a long-range, high-frequency ground wave radar,” IEEE Journal of Oceanic Engineering, vol. 19, no. 4, pp. 549–554, 1994.
- W. Huang, S.-C. Wu, E. Gill, B.-Y. Wen, and J.-C. Hou, “HF radar wave and wind measurement over the Eastern China Sea,” IEEE Transactions on Geoscience and Remote Sensing, vol. 40, no. 9, pp. 1950–1955, 2002.
- Y.-W. Tian, B.-Y. Wen, and H. Zhou, “Measurement of high and low waves using dual-frequency broad-beam HF radar,” IEEE Geoscience and Remote Sensing Letters, vol. 11, no. 9, pp. 1599–1603, 2014.
- C. Randell, D. Power, and R. Khan, “Ocean clutter characteristics for a coherent Doppler UHF radar,” in Proceedings of the MTS/IEEE Conference and Exhibition (OCEANS '00), Conference Proceedings (Cat. No.00CH37158), vol. 3, pp. 1515–1520, Providence, RI, USA, 2000.
- C. C. Teague, D. E. Barrick, P. M. Lilleboe, R. T. Cheng, and C. A. Ruhl, “UHF RiverSonde observations of water surface velocity at Threemile Slough, California,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '05), pp. 4383–4386, Seoul, South Korea, July 2005.
- C. C. Teague, D. E. Barrick, P. M. Lilleboe, R. T. Cheng, and C. A. Ruhl, “Long-term UHF RiverSonde river velocity observations at Castle Rock, Washington and Threemile Slough, California,” in Proceedings of the IEEE/OES 8th Working Conference on Current Measurement Technology, pp. 85–89, June 2005.
- C. C. Teague, D. E. Barrick, P. M. Lilleboe, H. Roarty, D. Holden, and D. Goldinger, “Extended-range RiverSonde operation on the Hudson River,” in Proceedings of the IEEE/OES/CWTM 10th Working Conference on Current, Waves and Turbulence Measurement (CWTM '11), pp. 78–80, Monterey, Calif, USA, March 2011.
- C. C. Teague, D. E. Barrick, P. M. Lilleboe, R. T. Cheng, P. Stumpner, and J. R. Burau, “Dual-RiverSonde measurements of two-dimensional river flow patterns,” in Proceedings of the IEEE/OES/CMTC 9th Working Conference on Current Measurement Technology, pp. 258–263, Charleston, SC, USA, March 2008.
- Z. Ma, B. Wen, H. Zhou, C. Wang, and W. Yan, “UHF surface currents radar hardware system design,” IEEE Microwave and Wireless Components Letters, vol. 15, no. 12, pp. 904–906, 2005.
- S. Wang, B. Wen, C. Wang, Z. Yan, K. Li, and J. Yang, “UHF surface dynamics parameters radar design and experiment,” IEEE Microwave and Wireless Components Letters, vol. 24, no. 1, article 12, pp. 65–67, 2014.
- D. Jiang, B. Wen, C. Wang, K. Li, Z. Yan, and J. Yang, “A new type full digital uhf radar system design,” IEICE Electronics Express, vol. 10, no. 11, pp. 1–6, 2013.
- K. Li, B. Wen, Y. Xu, J. Tan, J. Yang, and Y. Liu, “A novel UHF radar system design for river dynamics monitoring,” IEICE Electronics Express, vol. 12, no. 4, 2015.
- R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986.
- T.-J. Shan, M. Wax, and T. Kailath, “On spatial smoothing for direction-of-arrival estimation of coherent signals,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 33, no. 4, pp. 806–811, 1985.
- R. T. Williams, S. Prasad, A. K. Mahalanabis, and L. H. Sibul, “An improved spatial smoothing technique for bearing estimation in a multipath environment,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 36, no. 4, pp. 425–432, 1988.
Copyright © 2017 Yidong Hou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.