CN105467393B - A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system - Google Patents

A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system Download PDF

Info

Publication number
CN105467393B
CN105467393B CN201610009058.1A CN201610009058A CN105467393B CN 105467393 B CN105467393 B CN 105467393B CN 201610009058 A CN201610009058 A CN 201610009058A CN 105467393 B CN105467393 B CN 105467393B
Authority
CN
China
Prior art keywords
radial velocity
data
dimensional wind
velocity data
millimeter wave
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610009058.1A
Other languages
Chinese (zh)
Other versions
CN105467393A (en
Inventor
冯凯
周亭亭
魏艳强
李岩松
王志锐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace new weather Technology Co., Ltd
Original Assignee
Beijing Institute of Radio Measurement
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Radio Measurement filed Critical Beijing Institute of Radio Measurement
Priority to CN201610009058.1A priority Critical patent/CN105467393B/en
Publication of CN105467393A publication Critical patent/CN105467393A/en
Application granted granted Critical
Publication of CN105467393B publication Critical patent/CN105467393B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention relates to a kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system, wherein, this method includes：Control millimeter wave cloud radar to press the PPI scan patterns of the default angle of pitch and carry out meteorological target acquisition, signal transacting is carried out to target acquisition data and exports original radial velocity data；The original radial velocity data of signal processing module output are gathered, and are wrapped into being transferred to main control module after processing；Original radial velocity data after packing is handled are preserved；Original radial velocity data of the reading and saving in main control module are simultaneously carried out with pretreatment to it, obtain being available for the radial velocity data of dimensional wind inverting；The radial velocity data for being available for dimensional wind inverting are read, and inversion equation is built to carry out dimensional wind inverting to corresponding target acquisition data under default solution criterion, obtain dimensional wind speed data.By the present invention, solve the problems, such as that accurately and efficiently dimensional wind product can not be obtained using millimeter wave cloud radar.

Description

A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system
Technical field
The present invention relates to Radar Technology field, more particularly to a kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and it is System.
Background technology
Wind field information is one of main meteorological factor of Atmospheric Survey.The meteorological support of large-scale activity, aerodrome guarantee, All there is tight demand in the fields such as weather forecast to the fine observation of wind field.Except instruments such as anemoclinograph, wind profile radars Directly survey outside wind, it is also a kind of to measure wind field information indirectly using the tracer grain (such as water dust, raindrop, smog) of float Important method.Just because of this, the weather radar output using Wind-field Retrieval information as a routine of current operation operation Product.
But be present many weak points in weather radar institute inversion of three dimensional wind field, be mainly manifested in：First, weather radar is anti- The tracer for drilling wind Place Attachment is raindrop and stronger water dust particle, and the probability of these particles middle appearance on high is not ten Divide it is numerous, which results in the most of the time effectively survey wind data blank；Second, because the object of observation is different, day The common angle of pitch setting of gas radar is relatively low, and this make it that horizontal radius reaches several hundred kilometers during circular scanning, it is difficult to meet linear Uniform Wind-field Retrieval is it is assumed that cause larger Wind-field Retrieval error；3rd, limited by missing object, weather radar institute energy The altitude range of enough WIND FIELDSs is smaller.
Compared to weather radar, the shorter millimeter wave cloud radar of booming wavelength can detect aerial weaker in recent years Particle, therefore possess more preferable potentiality in terms of dimensional wind inverting.Millimeter wave cloud radar is mainly used in cloud, mist, weak precipitation Etc. the detection of target, using scattering process of the small particles to electromagnetic wave, the echo-signal of continuous measurement website overhead meteorology target, The meteorological target information of high-spatial and temporal resolution is obtained, there is comprehensive, round-the-clock observing capacity, be widely used in atmospheric science Research, weather modification, cloud automatic observation, aerodrome guarantee, military meteorological support etc. and field.
Due to itself system (wave beam is narrow, correlation is weak) and the characteristic of measured target, (intensity is weak, environment and clutter are disturbed It is sensitive), the detection data scale of millimetre-wave radar is big, individual data randomness is strong, continuity and stability are poor.These features Also the method for solving to dimensional wind inverting is set and inverting accuracy causes certain difficulty.
The content of the invention
The technical problems to be solved by the invention are in view of the shortcomings of the prior art, there is provided a kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system, solve accurately and efficiently can not obtain dimensional wind product using millimeter wave cloud radar Problem.The technical scheme that the present invention solves above-mentioned technical problem is as follows：
A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method, comprises the following steps：
Step 1：Control millimeter wave cloud radar to press the PPI scan patterns of the default angle of pitch and carry out meteorological target acquisition, and it is raw Into target acquisition data；
Step 2：Signal transacting is carried out to the target acquisition data and exports original radial velocity data；
Step 3：The original radial velocity data of the signal processing module output are gathered, and are carried out packing processing；
Step 4：Original radial velocity data after packing is handled are preserved；
Step 5：The original radial velocity data of reading and saving simultaneously carry out eliminating millimetre-wave radar echo-signal to it Instable pretreatment, obtain being available for the radial velocity data of dimensional wind inverting；
Step 6：The radial velocity data of dimensional wind inverting are available for described in reading, and are built instead under default solution criterion Equation is drilled to carry out dimensional wind inverting to the corresponding target acquisition data, obtains dimensional wind speed data.
The beneficial effects of the invention are as follows：The inverting of dimensional wind has been carried out using millimeter wave cloud radar, has taken full advantage of milli The advantages of VHF band radar is strong to weak signal target particle detection ability, the reliability and validity for surveying wind data are enhanced, is enriched The data product type of millimeter wave cloud detection radar.In addition, inversion equation is established from data prediction and in terms of solving criterion two, So as to ensure that validity and accuracy using original radial data inversion of three dimensional wind field.
On the basis of above-mentioned technical proposal, the present invention can also do following improvement.
On the basis of above-mentioned technical proposal, the present invention can also do following improvement.
Further, the horizontal sweep speed of millimeter wave cloud radar is 1~2 degrees second, and the maximum of the default angle of pitch is no more than 82 degree.
Beneficial effect using above-mentioned further scheme is that rational sweep speed both ensure that wave beam had to same target Enough residence times, the overlong time needed for refutation process will not be made again.The reasonable setting of the angle of pitch had then both improved simultaneously To the ability of higher vacancy Wind-field Retrieval, turn avoid because the angle of pitch is excessive and the measurement error of Horizontal Winds is larger.
Further, the pretreatment of the step 5 is implemented as：
Step 5.1：Ascending order arrangement is carried out by azimuthal size to the original radial velocity data；
Step 5.2：The predetermined threshold value of one radial velocity data is set to remove in the original radial velocity data Singular point；
Step 5.3：Using 5 × 1 filter window by azimuth sequence and the original radial velocity to eliminating singular point Data carry out smothing filtering, obtain associated parameter data.
Beneficial effect using above-mentioned further scheme is by pre-processing the useful letter in original radial velocity data Breath extracts, and strengthens its availability, ensure that follow-up solution equation can solve and result is accurate.
Further, the step 5.3 is implemented as：
Step 5.3.1:If it is M that scan control module, which carries out total the points of measurement that meteorological target acquisition obtains, total iteration is secondary Number is N；
Step 5.3.2:Since 1, for i-th point, the average m of the radial velocity data of 5 observation stations around it is calculated；
Step 5.3.3:If the difference of one of them described observation station and the average m absolute value>=0.5dBZ, then remove , otherwise retain at described i-th point at described i-th point；
Step 5.3.4:Next point of the observation station is calculated, M is arrived until calculating；
Step 5.3.5:Next iteration is carried out, until iterating to times N.
Using it is above-mentioned further scheme beneficial effect be by using single radial velocity data in its orientation it is adjacent The information of data, the data larger to the error caused by the reasons such as spectrum discrimination error, the disturbance of environment clutter are removed, ensured The precision of follow-up Wind-field Retrieval.
Further, the step 6 is implemented as：
Step 6.1：Assuming that scan control module carries out the azimuth difference for the different observation stations that meteorological target acquisition obtains For β1,β2,......,βM, M is the total number of observation station, establishes the inversion matrix of M rows：
Wherein, VRFor precipitation particles radial velocity, VAFor the horizontal component of precipitation particles radial velocity, VEFor precipitation particles The vertical component of radial velocity, horizontal component further can be analyzed to two parallel to direct north and perpendicular to direct north Component VANAnd VAE, β0For horizontal wind direction and the angle of direct north, α is antenna elevation angle, according to relational expression XV=F, you can obtain Dimensional wind speed V.
Step 6.2：Establish inversion equation：
Wherein, VR(βj) it is that azimuth is βjActual measurement radial velocity, VR(βj) ' obtained after least square fitting In βjThe radial velocity at place；
After solving equation (1), you can solve VAN、VAEAnd VEWind speed component on this three directions.
Beneficial effect using above-mentioned further scheme has fully been excavated in all radial velocities obtained after pretreatment The wind field information contained, and consider that composition error of all single radial velocity data in the case where expected wind field is assumed is minimum so that The wind field precision that institute's inverting obtains is higher.
Another technical scheme that the present invention solves above-mentioned technical problem is as follows：
A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval system, including scan control module, signal processing module, information are adopted Collect module, main control module and data inversion module, wherein, the data inversion module includes pretreatment module and data solve mould Block；
The scan control module, for controlling millimeter wave cloud radar to enter promoting the circulation of qi by the PPI scan patterns of the default angle of pitch Signal processing module is transferred to as target acquisition, and by caused target acquisition data；
The signal processing module, for carrying out signal transacting to the target acquisition data and exporting original radial velocity Data；
Described information acquisition module, for gathering the original radial velocity data of the signal processing module output, and will Main control module is transferred to after its processing of packing；
The main control module, preserved for the original radial velocity data after packing is handled；
The pretreatment module, for original radial velocity data of the reading and saving in the main control module and it is entered Row obtains being available for the radial velocity of dimensional wind inverting to eliminate the instable pretreatment of millimetre-wave radar echo-signal Data；
The data solve module, for reading the radial velocity data for being available for dimensional wind inverting, and default Solve and inversion equation is built under criterion to carry out dimensional wind inverting to the corresponding target acquisition data, obtain dimensional wind Speed data.
The beneficial effects of the invention are as follows：The inverting of dimensional wind has been carried out using millimeter wave cloud radar, has taken full advantage of milli The advantages of VHF band radar is strong to weak signal target particle detection ability, the reliability and validity for surveying wind data are enhanced, is enriched The data product type of millimeter wave cloud detection radar.In addition, inversion equation is established from data prediction and in terms of solving criterion two, So as to ensure that validity and accuracy using original radial data inversion of three dimensional wind field.
Further, the pretreatment module includes：
Sequencing unit, for carrying out ascending order arrangement by azimuthal size to the original radial velocity data；
Unit is deleted, for setting the predetermined threshold value of a radial velocity data to remove the original radial velocity data In singular point；
Filter unit, for the filter window using 5 × 1 by azimuth sequence and the original radial direction to eliminating singular point Speed data carries out smothing filtering, obtains associated parameter data.
Beneficial effect using above-mentioned further scheme is by pre-processing the useful letter in original radial velocity data Breath extracts, and strengthens its availability, ensure that follow-up solution equation can solve and result is accurate.
Further, the filter unit includes：
Assuming that subelement, total the points of measurement that meteorological target acquisition obtains is carried out as M for setting scan control module, is always changed Generation number is N；
First computation subunit, for since 1, for i-th point, calculating the radial velocity number of 5 observation stations around it According to average m；
Judgment sub-unit, for judging whether the difference of absolute value of one of them described observation station and the average m is more than Equal to 0.5dBZ, if it is, removing at described i-th point, otherwise retain at described i-th point；
Second computation subunit, for calculating next point of the observation station, M is arrived until calculating；
Iteration subelement, for carrying out next iteration, until iterating to times N.
Using it is above-mentioned further scheme beneficial effect be by using single radial velocity data in its orientation it is adjacent The information of data, the data larger to the error caused by the reasons such as spectrum discrimination error, the disturbance of environment clutter are removed, ensured The precision of follow-up Wind-field Retrieval.
Further, the data solve module and included：
Inversion matrix unit, for carrying out the side for the different observation stations that meteorological target acquisition obtains according to scan control module The total number of parallactic angle and observation station, establishes inversion matrix；
Inversion equation unit, for establishing inversion equation, and after equation is solved, you can solve VAN、VAEAnd VEThis three Wind speed component on direction.
Beneficial effect using above-mentioned further scheme has fully been excavated in all radial velocities obtained after pretreatment The wind field information contained, and consider that composition error of all single radial velocity data in the case where expected wind field is assumed is minimum so that The wind field precision that institute's inverting obtains is higher.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of millimeter wave cloud radar three-dimensional Wind-field Retrieval method of the present invention；
Fig. 2 is the structural representation of millimeter wave cloud radar three-dimensional Wind-field Retrieval system of the present invention；
In accompanying drawing, the list of parts representated by each label is as follows：
1st, scan control module, 2, signal processing module, 3, information acquisition module, 4, main control module, 5, data inversion mould Block, 6, pci bus, 7, shared drive.
Embodiment
The principle and feature of the present invention are described below in conjunction with accompanying drawing, the given examples are served only to explain the present invention, and It is non-to be used to limit the scope of the present invention.
Fig. 1 is the schematic flow sheet of millimeter wave cloud radar three-dimensional Wind-field Retrieval method of the present invention；As shown in figure 1, this method Comprise the following steps：
Step 1：Millimeter wave cloud radar is controlled to press PPI (Plan Position Indicator, the abbreviation of the default angle of pitch For " circumferential plane position display ") scan pattern carries out meteorological target acquisition, and generates target acquisition data；
Step 2：Signal transacting is carried out to target acquisition data and exports original radial velocity data；
Step 3：The original radial velocity data of signal processing module output are gathered, and are carried out packing processing；
Step 4：Original radial velocity data after packing is handled are preserved；
Step 5：The original radial velocity data of reading and saving simultaneously carry out eliminating millimetre-wave radar echo-signal to it Instable pretreatment, obtain being available for the radial velocity data of dimensional wind inverting；
Step 6：The radial velocity data for being available for dimensional wind inverting are read, and in default structure inverting side under solving criterion Journey obtains dimensional wind speed data to carry out dimensional wind inverting to corresponding target acquisition data.
The present invention has carried out the inverting of dimensional wind by using millimeter wave cloud radar, takes full advantage of millimere-wave band radar The advantages of strong to weak signal target particle detection ability, the reliability and validity for surveying wind data are enhanced, enrich millimeter wave and survey cloud The data product type of radar.In addition, inversion equation is established from data prediction and in terms of solving criterion two, so as to ensure that Utilize the validity and accuracy of original radial data inversion of three dimensional wind field.
Preferably, suitable pulse accumulation number is set to ensure that the residence time of wave beam is less than beam angle divided by radar water Simple scan speed, the maximum for presetting the angle of pitch are no more than 82 degree, and set default angle of pitch minimum value should ensure that at a distance Velocity ambiguity does not occur for the radial velocity of place's radar detection.
Set rational sweep speed both to ensure that wave beam there are enough residence times to same target, inverting will not be made again Overlong time needed for process.The reasonable setting of the angle of pitch had then both improved the ability to higher vacancy Wind-field Retrieval simultaneously, again Avoid because the angle of pitch is excessive and the measurement error of Horizontal Winds is larger.
Preferably, the pretreatment of step 5 is implemented as：
Step 5.1：Ascending order arrangement is carried out by azimuthal size to original radial velocity data；
Step 5.2：Set the predetermined threshold value of a radial velocity data unusual in original radial velocity data to remove Point；Under normal circumstances, lower threshold takes -100m/s, and upper limit threshold takes+100m/s；
Step 5.3：Using 5 × 1 filter window by azimuth sequence and the original radial velocity to eliminating singular point Data carry out smothing filtering, obtain associated parameter data.
The useful information in original radial velocity data can be extracted by pretreatment, strengthen its availability, ensured Follow-up solution equation can solve and result is accurate.
Preferably, step 5.3 is implemented as：
Step 5.3.1:If it is M that scan control module, which carries out total the points of measurement that meteorological target acquisition obtains, total iteration is secondary Number is N；
Step 5.3.2:Since 1, for i-th point, the average m of the radial velocity data of 5 observation stations around it is calculated；
Step 5.3.3:If the difference of one of observation station and average m absolute value>=0.5dBZ, then remove at i-th point, no Then retain at i-th point；
Step 5.3.4:Next point of calculating observation point, M is arrived until calculating；
Step 5.3.5:Next iteration is carried out, until iterating to times N.
By using the information of single radial velocity data adjacent data in its orientation, to because of spectrum discrimination error, environment The larger data of error caused by reason such as clutter disturbance are removed, and ensure that the precision of follow-up Wind-field Retrieval.
Preferably, step 6 is implemented as：
Step 6.1：Assuming that scan control module carries out the azimuth difference for the different observation stations that meteorological target acquisition obtains For β1,β2,......,βM, M is the total number of observation station, establishes the inversion matrix of M rows：
Wherein, VRFor precipitation particles radial velocity, VAFor the horizontal component of precipitation particles radial velocity, VEFor precipitation particles The vertical component of radial velocity, horizontal component preferably can be analyzed to two parallel to direct north and perpendicular to direct north Component VANAnd VAE, β0For horizontal wind direction and the angle of direct north, α is antenna elevation angle, according to relational expression XV=F, you can obtain Dimensional wind speed V.
Step 6.2：Establish inversion equation：
Wherein, VR(βj) it is that azimuth is βjActual measurement radial velocity, VR(βj) ' obtained after least square fitting In βjThe radial velocity at place；
After solving equation (1), you can solve VAN、VAEAnd VEWind speed component on this three directions.
So, the wind field information contained in all radial velocities obtained after pretreatment has fully been excavated, and has been considered all Composition error of the single radial velocity data in the case where expected wind field is assumed is minimum so that the wind field precision that institute's inverting obtains is higher.
Fig. 2 is the structural representation of millimeter wave cloud radar three-dimensional Wind-field Retrieval system of the present invention；As shown in Fig. 2 the system Including scan control module 1, signal processing module 2, information acquisition module 3, main control module 4 and data inversion module 5, wherein, Data inversion module 5 includes pretreatment module and data solve module；
Scan control module 1, for controlling millimeter wave cloud radar to carry out meteorological mesh by the PPI scan patterns of the default angle of pitch Mark detection, and caused target acquisition data are transferred to signal processing module 2；
Signal processing module 2, for carrying out signal transacting to target acquisition data and exporting original radial velocity data；
Information acquisition module 3, for gathering the original radial velocity data of the output of signal processing module 2, and it is wrapped into Main control module 4 is transferred to after processing；
Main control module 4, preserved for the original radial velocity data after packing is handled；
Pretreatment module, for original radial velocity data of the reading and saving in main control module 4 and it is carried out to The instable pretreatment of millimetre-wave radar echo-signal is eliminated, obtains being available for the radial velocity data of dimensional wind inverting；
Data solve module, the radial velocity data of dimensional wind inverting are available for for reading, and solve criterion default Lower structure inversion equation obtains dimensional wind speed data to carry out dimensional wind inverting to corresponding target acquisition data.
Wherein, scan control module 1, signal processing module 2, information acquisition module 3 and data inversion module 5 respectively with master Control module 4 is connected by pci bus 6, and main control module 4 includes pci bus 6 and shared drive 7, and scan control module 1 passes through PCI Target acquisition data are transferred to signal processing module 2 by bus 6；Signal processing module 2 is carried out at signal to target acquisition data Manage and export；Information acquisition module 3 gathers the target acquisition data that signal processing module 2 exports, and leads to after being wrapped into processing Pci bus 6 is crossed to be saved in shared drive 7；Data inversion module 5 is read in shared drive 7 by packing by pci bus 6 Target acquisition data after processing simultaneously carry out dimensional wind inverting, obtain dimensional wind speed data.
When radar is in detection process, information acquisition module 3 gathers the target acquisition data that signal processing module 2 exports, And extract radial velocity data therein, it is saved in after packing processing on the main control module 4 of industrial control computer.One After the completion of group PPI scanning observations, terminate the data storage operations, and observation data are saved as into a data file.
The present invention has carried out the inverting of dimensional wind by using millimeter wave cloud radar, takes full advantage of millimere-wave band radar The advantages of strong to weak signal target particle detection ability, the reliability and validity for surveying wind data are enhanced, enrich millimeter wave and survey cloud The data product type of radar.In addition, inversion equation is established from data prediction and in terms of solving criterion two, so as to ensure that Utilize the validity and accuracy of original radial data inversion of three dimensional wind field.
Preferably, pretreatment module includes：
Sequencing unit, for carrying out ascending order arrangement by azimuthal size to original radial velocity data；
Unit is deleted, for setting the predetermined threshold value of a radial velocity data to remove in original radial velocity data Singular point；
Filter unit, for the filter window using 5 × 1 by azimuth sequence and the original radial direction to eliminating singular point Speed data carries out smothing filtering, obtains associated parameter data.
The useful information in original radial velocity data is extracted by pre-processing, strengthens its availability, ensure that Follow-up solution equation can solve and result is accurate.
Preferably, filter unit includes：
Assuming that subelement, total the points of measurement that meteorological target acquisition obtains is carried out as M for setting scan control module 1, always Iterations is N；
First computation subunit, for since 1, for i-th point, calculating the radial velocity number of 5 observation stations around it According to average m；
Judgment sub-unit, for judging whether the one of observation station and average m difference of absolute value is more than or equal to 0.5dBZ, if it is, removing at i-th point, otherwise retain at i-th point；
Second computation subunit, for next point of calculating observation point, M is arrived until calculating；
Iteration subelement, for carrying out next iteration, until iterating to times N.
By using the information of single radial velocity data adjacent data in its orientation, to because of spectrum discrimination error, environment The larger data of error caused by reason such as clutter disturbance are removed, and ensure that the precision of follow-up Wind-field Retrieval.
Preferably, data solve module and included：
Inversion matrix unit, for carrying out different observation stations that meteorological target acquisition obtains according to scan control module 1 Azimuth and the total number of observation station, establish inversion matrix；
Inversion equation unit, for establishing inversion equation, and after equation is solved, you can solve VAN、VAEAnd VEThis three Wind speed component on direction.
So, the wind field information contained in all radial velocities obtained after pretreatment has fully been excavated, and has been considered all Composition error of the single radial velocity data in the case where expected wind field is assumed is minimum so that the wind field precision that institute's inverting obtains is higher.
The foregoing is only presently preferred embodiments of the present invention, be not intended to limit the invention, it is all the present invention spirit and Within principle, any modification, equivalent substitution and improvements made etc., it should be included in the scope of the protection.

Claims (8)

  1. A kind of 1. millimeter wave cloud radar three-dimensional Wind-field Retrieval method, it is characterised in that comprise the following steps：
    Step 1：Control millimeter wave cloud radar to press the PPI scan patterns of the default angle of pitch and carry out meteorological target acquisition, and generate mesh Mark detection data；
    Step 2：Signal transacting is carried out to the target acquisition data and exports original radial velocity data；
    Step 3：The original radial velocity data of the signal processing module output are gathered, and are carried out packing processing；
    Step 4：Original radial velocity data after packing is handled are preserved；
    Step 5：Original radial velocity data of the reading and saving in main control module simultaneously carry out eliminating millimetre-wave radar to it The instable pretreatment of echo-signal, obtain being available for the radial velocity data of dimensional wind inverting；
    Step 6：The radial velocity data of dimensional wind inverting are available for described in reading, and in default structure inverting side under solving criterion Journey obtains dimensional wind speed data to carry out dimensional wind inverting to the corresponding target acquisition data；
    The step 6 is implemented as：
    Step 6.1：Assuming that the azimuth that scan control module carries out the different observation stations that meteorological target acquisition obtains is respectively β1, β2,......,βM, M is the total number of observation station, establishes the inversion matrix of M rows：
    Wherein, VRFor precipitation particles radial velocity, VAFor the horizontal component of precipitation particles radial velocity, VEFor precipitation particles radially The vertical component of speed, horizontal component further can be analyzed to two components parallel to direct north and perpendicular to direct north VANAnd VAE, β0For horizontal wind direction and the angle of direct north, α is antenna elevation angle, according to relational expression XV=F, you can obtain three-dimensional Wind field speed V;
    Step 6.2：Establish inversion equation：
    Wherein, VR(βj) it is that azimuth is βjActual measurement radial velocity, VR(βj) ' it is to be obtained after least square fitting in βj The radial velocity at place；
    After solving equation (1), you can solve VAN、VAEAnd VEWind speed component on this three directions.
  2. 2. millimeter wave cloud radar three-dimensional Wind-field Retrieval method according to claim 1, it is characterised in that millimeter wave cloud radar Horizontal sweep speed is 1~2 degrees second, and the maximum for presetting the angle of pitch is no more than 82 degree, and minimum value is not it is ensured that radial velocity is sent out It is raw fuzzy.
  3. 3. millimeter wave cloud radar three-dimensional Wind-field Retrieval method according to claim 1, it is characterised in that the step 5 it is pre- Processing is implemented as：
    Step 5.1：Ascending order arrangement is carried out by azimuthal size to the original radial velocity data；
    Step 5.2：Set the predetermined threshold value of a radial velocity data unusual in the original radial velocity data to remove Point；
    Step 5.3：Using 5 × 1 filter window by azimuth sequence and the original radial velocity data to eliminating singular point Smothing filtering is carried out, obtains associated parameter data.
  4. 4. millimeter wave cloud radar three-dimensional Wind-field Retrieval method according to claim 3, it is characterised in that the step 5.3 has Body is embodied as：
    Step 5.3.1:If it is M that scan control module, which carries out total the points of measurement that meteorological target acquisition obtains, total iterations is N；
    Step 5.3.2:Since 1, for i-th point, the average m of the radial velocity data of 5 observation stations around it is calculated；
    Step 5.3.3:If the difference of one of them described observation station and the average m absolute value>=0.5dBZ, then described in removal , otherwise retain at described i-th point at i-th point；
    Step 5.3.4:Next point of the observation station is calculated, M is arrived until calculating；
    Step 5.3.5:Next iteration is carried out, until iterating to times N.
  5. 5. a kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval system, it is characterised in that including scan control module, signal transacting mould Block, information acquisition module, main control module and data inversion module, wherein, the data inversion module include pretreatment module and Data solve module；
    The scan control module, for controlling millimeter wave cloud radar to carry out meteorological mesh by the PPI scan patterns of the default angle of pitch Mark detection, and caused target acquisition data are transferred to signal processing module；
    The signal processing module, for carrying out signal transacting to the target acquisition data and exporting original radial velocity number According to；
    Described information acquisition module, for gathering the original radial velocity data of signal processing module output, and by its dozen Main control module is transferred to after bag processing；
    The main control module, preserved for the original radial velocity data after packing is handled；
    The pretreatment module, for original radial velocity data of the reading and saving in the main control module and it is used To eliminate the instable pretreatment of millimetre-wave radar echo-signal, obtain being available for the radial velocity number of dimensional wind inverting According to；
    The data solve module, are solved for reading the radial velocity data for being available for dimensional wind inverting, and default Inversion equation is built under criterion to carry out dimensional wind inverting to the corresponding target acquisition data, obtains dimensional wind speed Data；
    The data, which solve module, to be included：
    Inversion matrix unit, for carrying out the azimuth for the different observation stations that meteorological target acquisition obtains according to scan control module With the total number of observation station, inversion matrix is established；
    Wherein, VRFor precipitation particles radial velocity, VAFor the horizontal component of precipitation particles radial velocity, VEFor precipitation particles radially The vertical component of speed, horizontal component further can be analyzed to two components parallel to direct north and perpendicular to direct north VANAnd VAE, β0For horizontal wind direction and the angle of direct north, α is antenna elevation angle, according to relational expression XV=F, you can obtain three-dimensional Wind field speed V；
    Wherein：The azimuth that scan control module carries out the different observation stations that meteorological target acquisition obtains is respectively β1, β2,......,βM, M is the total number of observation station；
    Inversion equation unit, for establishing inversion equation；
    Wherein, VR(βj) it is that azimuth is βjActual measurement radial velocity, VR(βj) ' it is to be obtained after least square fitting in βj The radial velocity at place；
    And after equation is solved, you can solve VAN、VAEAnd VEWind speed component on this three directions.
  6. 6. millimeter wave cloud radar three-dimensional Wind-field Retrieval system according to claim 5, it is characterised in that the pretreatment module Including：
    Sequencing unit, for carrying out ascending order arrangement by azimuthal size to the original radial velocity data；
    Unit is deleted, for setting the predetermined threshold value of a radial velocity data to remove in the original radial velocity data Singular point；
    Filter unit, for the filter window using 5 × 1 by azimuth sequence and the original radial velocity to eliminating singular point Data carry out smothing filtering, obtain associated parameter data.
  7. 7. millimeter wave cloud radar three-dimensional Wind-field Retrieval system according to claim 6, it is characterised in that the filter unit bag Include：
    Assuming that subelement, total the points of measurement that meteorological target acquisition obtains is carried out as M for setting scan control module, total iteration is secondary Number is N；
    First computation subunit, for since 1, for i-th point, calculating the radial velocity data of 5 observation stations around it Average m；
    Judgment sub-unit, for judging whether the difference of absolute value of one of them described observation station and the average m is more than or equal to 0.5dBZ, if it is, removing at described i-th point, otherwise retain at described i-th point；
    Second computation subunit, for calculating next point of the observation station, M is arrived until calculating；
    Iteration subelement, for carrying out next iteration, until iterating to times N.
  8. 8. millimeter wave cloud radar three-dimensional Wind-field Retrieval system according to claim 7, it is characterised in that the data solve mould Block includes：
    Inversion matrix unit, for carrying out the azimuth for the different observation stations that meteorological target acquisition obtains according to scan control module With the total number of observation station, inversion matrix is established；
    Inversion equation unit, for establishing inversion equation, and after equation is solved, you can solve VAN、VAEAnd VEThis three directions On wind speed component.
CN201610009058.1A 2016-01-06 2016-01-06 A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system Active CN105467393B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610009058.1A CN105467393B (en) 2016-01-06 2016-01-06 A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610009058.1A CN105467393B (en) 2016-01-06 2016-01-06 A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system

Publications (2)

Publication Number Publication Date
CN105467393A CN105467393A (en) 2016-04-06
CN105467393B true CN105467393B (en) 2017-12-15

Family

ID=55605293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610009058.1A Active CN105467393B (en) 2016-01-06 2016-01-06 A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system

Country Status (1)

Country Link
CN (1) CN105467393B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106199605B (en) * 2016-07-06 2018-08-24 西南技术物理研究所 Wind field error correcting method
CN106526709B (en) * 2016-09-26 2020-02-21 北京无线电测量研究所 Airport cloud and mist comprehensive early warning system and method based on multi-source equipment networking
CN107741587A (en) * 2017-09-15 2018-02-27 北京无线电测量研究所 A kind of meteorological detection method and system of dimensional wind
CN107843895B (en) * 2017-10-20 2019-10-18 厦门市气象灾害防御技术中心(海峡气象开放实验室、厦门市避雷检测技术中心) A kind of Dual-Doppler weather radar dimensional wind inversion method
CN107831481A (en) * 2017-10-30 2018-03-23 厦门市气象灾害防御技术中心(海峡气象开放实验室厦门市避雷监测技术中心) The lower three Doppler radar three-dimensional wind field inversion methods of dynamic earth coordinate system
CN108269303B (en) * 2017-12-22 2021-08-24 广东纳睿雷达科技股份有限公司 Three-dimensional meteorological radar display method
CN109283535B (en) * 2018-11-16 2022-12-20 北京无线电测量研究所 Detection method of non-uniform wind field
CN109849930B (en) * 2018-12-29 2020-08-04 百度在线网络技术（北京）有限公司 Method and device for calculating speed of adjacent vehicle of automatic driving automobile
CN110288856A (en) * 2019-06-21 2019-09-27 中国民用航空总局第二研究所 The Scheduled Flight monitoring system and method for fine forecast based on wind
CN111624623B (en) * 2020-06-02 2023-03-10 青岛镭测创芯科技有限公司 Wind field inversion method based on laser radar non-uniform scanning

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006133203A (en) * 2004-11-09 2006-05-25 Nec Engineering Ltd Wind profiler system
CN101988963A (en) * 2010-04-19 2011-03-23 南京恩瑞特实业有限公司 Method for acquiring three-dimensional wind field by using wind profiler radar
CN102508222A (en) * 2011-11-18 2012-06-20 中国科学技术大学 Middle and upper atmospheric wind field retrieval method
CN102890272A (en) * 2012-11-05 2013-01-23 中国航天科工集团第二研究院二十三所 Method for processing millimeter wave cloud radar signal
CN103197299A (en) * 2013-03-25 2013-07-10 南京信息工程大学 Extraction and quantitative analysis system of weather radar radial wind information
CN104345312A (en) * 2014-10-27 2015-02-11 北京无线电测量研究所 Data fusion method and system of millimeter-wave cloud radar

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006133203A (en) * 2004-11-09 2006-05-25 Nec Engineering Ltd Wind profiler system
CN101988963A (en) * 2010-04-19 2011-03-23 南京恩瑞特实业有限公司 Method for acquiring three-dimensional wind field by using wind profiler radar
CN102508222A (en) * 2011-11-18 2012-06-20 中国科学技术大学 Middle and upper atmospheric wind field retrieval method
CN102890272A (en) * 2012-11-05 2013-01-23 中国航天科工集团第二研究院二十三所 Method for processing millimeter wave cloud radar signal
CN103197299A (en) * 2013-03-25 2013-07-10 南京信息工程大学 Extraction and quantitative analysis system of weather radar radial wind information
CN104345312A (en) * 2014-10-27 2015-02-11 北京无线电测量研究所 Data fusion method and system of millimeter-wave cloud radar

Also Published As

Publication number Publication date
CN105467393A (en) 2016-04-06

Similar Documents

Publication Publication Date Title
CN105467393B (en) A kind of millimeter wave cloud radar three-dimensional Wind-field Retrieval method and system
Sicard et al. Mixed-layer depth determination in the Barcelona coastal area from regular lidar measurements: methods, results and limitations
CN102628944B (en) Stratus cloud and convective cloud automatic recognition method based on Doppler radar data
CN104345312B (en) Data fusion method and system of millimeter-wave cloud radar
CN113868971B (en) Airport area three-dimensional refined wind field reconstruction method based on numerical simulation model and historical wind field characteristics
Wang et al. Lake water footprint identification from time-series ICESat/GLAS data
CN105487076A (en) Millimeter-wave cloud radar dense fog visibility inversion method and system
CN103985127B (en) The detection method of small target of a kind of intensive star background and device
CN113325422B (en) Space-based rain radar target positioning and rainfall information three-dimensional processing method and system
CN106772387A (en) A kind of wind shear recognition methods
CN104200082A (en) Typhoon landing prediction method
CN109407177A (en) Dense fog identifying system and methods for using them based on machine learning and conventional meteorological observation
CN105468899B (en) Based on MICAPS wind field information low-level jet stream automatic identification and method for drafting
CN113866770A (en) Hail cloud early identification method and storage medium
CN110764083B (en) Anti-intrusion data fusion method and system for millimeter wave radar
CN104008403A (en) Multi-target identification and judgment method based on SVM mode
CN111736156B (en) Upwind area identification method and device based on weather radar
CN202433521U (en) Wind profile radar turbulence target detection processing plate
CN113219463A (en) Radar ground object echo identification method and system for power system
CN114047508B (en) Sea radar detection range statistical evaluation method and system
CN106203518B (en) A kind of reservoir watershed Rainstorm Feature factor discrimination method based on space clustering
CN113868970B (en) Airport area multi-source wind field fusion method based on numerical simulation model and mesoscale meteorological model
CN114779278A (en) CNN-RNN airport runway wind shear early warning method based on laser radar
CN111340908B (en) Rail electronic map generation method
CN113281716A (en) Photon counting laser radar data denoising method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191219

Address after: 214000 28 Weiming Road, Binhu District, Wuxi City, Jiangsu Province

Patentee after: Aerospace new weather Technology Co., Ltd

Address before: 100854, 32, Yongding Road, 50, Beijing, Haidian District

Patentee before: Beijing Institute of Radio Measurement