CN101706571B - Method for processing non-meteorological noise in radar echoes - Google Patents
Method for processing non-meteorological noise in radar echoes Download PDFInfo
- Publication number
- CN101706571B CN101706571B CN200910153740A CN200910153740A CN101706571B CN 101706571 B CN101706571 B CN 101706571B CN 200910153740 A CN200910153740 A CN 200910153740A CN 200910153740 A CN200910153740 A CN 200910153740A CN 101706571 B CN101706571 B CN 101706571B
- Authority
- CN
- China
- Prior art keywords
- clutter
- echo
- zone
- numerical value
- value
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for non-meteorological noise in radar echoes, which comprises the steps of: identifying noise areas existing on an 1.5 DEG elevation plane and an 0.5 DEG elevation plane by simply classifying the radar echoes, then analyzing reflectivity factor values of the echoes and then orderly referring radial velocity values and spectrum width values; performing a best noise rejection strategy on different noises after the noise areas are correctly identified; directly adopting a rejection strategy for noises under the conditions of clear sky and no precipitation echoes at peripheries; and adopting an approximate value filling mechanism for noises which are mixed in the precipitation echoes or are adjacent to the precipitation echoes so as to achieve better rejection effect.
Description
Technical field
The present invention relates to a kind of disposal route of radar return, especially relate to the disposal route of non-meteorological clutter in a kind of radar return.
Background technology
Doppler radar is to utilize the sexual intercourse object that meteorological object is found in electromagnetic back scattering; And measure locus, the strength distribution of meteorological object; Thereby the life of understanding precipitation disappears, develops and shifts to speed etc., but Doppler radar can detect the non-meteorological object equally., will cause the non-meteorological object clutter to pollute when being processed into the radar base data; The non-meteorological object is the non-meteorological clutter; The non-meteorological clutter comprises ground clutter and isolated point and speckle noise clutter; The transition clutter AP that the common ground clutter NP that ground clutter exists when comprising the radar beam normal propagation and the anomalous propagation of radar beam cause, isolated point and speckle noise clutter are mainly caused by the object that moves such as ship, aircraft, flock of birds, group of insects etc.Because the visual product and the corresponding algorithm of all deriving all are to be based upon on the basis of radar base data; The existence of non-meteorological clutter has just directly influenced the visual product of deriving like this; Responsive especially is that precipitation is estimated product, to the visual product influence of deriving bigger mainly be exactly ground clutter.Moszkowicz et al. (1994) and Jossand Lee (1995) have added up the WSR88D precipitation of 2 months and 60 months respectively and have estimated product; Only the precipitation evaluated error that the transition clutter AP that is caused by anomalous propagation causes reaches 15% and 13% respectively, so the removal problem of non-meteorological clutter also is one of emphasis of radar weather scientist research always.In weather radar of new generation; The common ground clutter NP (like mountain, buildings etc.) that the bypass figure that produces with RDASOT exists when suppressing the radar beam normal propagation; This method can be handled common ground clutter NP preferably; To the filtration of the transition clutter AP that mainly caused by the anomalous propagation of radar beam then is to suppress the zone and implement through defining a clutter; This clutter is suppressed regioselective inhibition degree act near the transition clutter AP of all radial velocities zero uniformly; This filter method obviously can reduce this clutter and suppress the reflectivity factor value in the zone, thereby will cause bigger influence to the visual product of deriving, and therefore the SA type radar in the whole nation does not all adopt this filter method that transition clutter AP is filtered basically.Yet the existence of these transition clutters AP usually can influence the confidence level of weather radar data; For this reason; People have studied the recognition methods of multiple ground clutter, and the external research that launches in this respect is a lot, as on WSR-88D radar basis; The application of chaos theory is classified to radar return and is produced serial algorithm etc. and carries out filtering, and the three-dimensional reflectance of Application Radar is carried out clutter filtering and applying neural network and classified automatically and discern or the like.All these methods all are to use the characteristic that the non-meteorological echo is showed on the visual product of deriving, and through algorithm reasonable in design, from software, carry out real-time additional processing filtration again; And the at present domestic research that launches in this respect is also few; On the other hand, these methods are that application chaos theory or neural network are all just rejected processing to clutter, in clutter is mixed in the precipitation zone; When perhaps clutter and precipitation zone is adjacent; Then generally all take not processing policy, so just can not effectively remove the clutter that is mixed in the precipitation zone or with the adjacent clutter in precipitation zone, will cause bigger influence to the visual product of deriving; In addition, for superrefraction clutter and ground clutter, because the interface is with a varied topography under the China, so the existing method of external some is local in the parameter that China all need carry out the great amount of samples statistics.
Summary of the invention
Technical matters to be solved by this invention provides a kind of effectively filtering and mixes adjacent or isolated clutter with true precipitation echo, and the disposal route of non-meteorological clutter in the good radar return of filtration result.
The present invention solves the problems of the technologies described above the technical scheme that is adopted: the disposal route of non-meteorological clutter in a kind of radar return; The radar base data of radar return comprises reflectivity factor numerical value, radial velocity numerical value and spectrum width numerical value, and this disposal route may further comprise the steps:
1. poor based on radar return reflectivity factor numerical value in vertical direction on the size of the reflectivity factor numerical value of radar return on 0.5 ° of elevation angle face and 0.5 ° and 1.5 ° of two elevation angle face; Radar return is divided into clear air echo, precipitation echo and dubiety echo, and precipitation echo comprises weak precipitation echo, medium precipitation echo and precipitation echo;
2. discern the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces; Detailed process is: 2.-1, on 1.5 ° of elevation angle faces in the radar return search area greater than the region area of setting and have the zone of low value spectrum width numerical value; Whether judgement should the zone be low value spectrum width zone, if then continue to carry out; Otherwise, do not handle; 2.-2, judge whether the zone of the reflection radial velocity numerical value that low value spectrum width zone is corresponding is low value radial velocity zone, if, then continue to carry out, otherwise, do not handle; 2.-3, judge whether have on 0.5 ° of elevation angle face of low value spectrum width zone below of 1.5 ° of elevation angle faces whether the zone of zone with strong reflection rate factor numerical value and should the zone corresponding reflection radial velocity numerical value is low value radial velocity zone; If having the zone of zone with strong reflection rate factor numerical value and should the zone corresponding reflection radial velocity numerical value is low value radial velocity zone; Confirm that then this zone with strong reflection rate factor numerical value is strong reflection rate zone; The area in this strong reflection rate zone is greater than the area in low value spectrum width zone, and definite low value spectrum width zone be doubtful superrefraction clutter, continuation execution then; Otherwise, do not handle; 2.-4, judge in the strong reflection rate zone whether have the dubiety echo; If exist; The level of then calculating reflectivity factor numerical value on 0.5 ° of elevation angle face change and 0.5 ° and 1.5 ° of two elevation angle face on the vertical change of reflectivity factor numerical value; The level of remembering reflectivity factor numerical value on 0.5 ° of elevation angle face is changed to Elbz, remembers that the vertical change of reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face is Vtbz, Elbz=Y (r) * | Z
I, j-Z
I, j-1|, Vtbz=Y (r) * (Z
Up-Z
Low), wherein, Y (r) representes the relevant weight of radar radial distance with the echo point position of working as pre-treatment, and its value is radar radial distance * 0.0025 of 1-when the echo point position of pre-treatment, and the unit of radar radial distance is km, Z
I, jRepresent that coordinate is (i, the reflectivity factor numerical value of echo point j), Z on 0.5 ° of elevation angle face
I, j-1Represent that coordinate is (i, the reflectivity factor numerical value of echo point j-1), Z on 0.5 ° of elevation angle face
UpThe reflectivity factor numerical value of representing the echo point on 1.5 ° of elevation angle faces, Z
LowRepresent on 0.5 ° of elevation angle face with 1.5 ° of elevation angle faces on the reflectivity factor numerical value of the corresponding echo point in echo point position; When the value of Elbz is that the value of 25~40db and Vtbz is when being 20~40db; Doubtful superrefraction clutter is labeled as the superrefraction clutter; Otherwise, doubtful superrefraction clutter is not handled;
3.-1, calculate the coefficient of the isolated echo point in the radar return on 1.5 ° of elevation angle faces 3. discern the noise clutter point in the radar return on 1.5 ° of elevation angle faces, detailed process is:, be designated as N
x,
Wherein, Mtotal=5 * 5, M for 5 * 5 region memories that are the center with this isolated echo point the number of non-NULL reflectivity factor numerical value; 3.-2, judgement factor N
xWhether less than the first threshold of setting, if then mark should isolated echo point be a noise clutter point, otherwise, do not handle isolating echo point;
4. superrefraction clutter in the radar return on 1.5 ° of elevation angle faces and noise clutter point are rejected processing; Detailed process is: for the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces; Judge the superrefraction clutter around whether have a precipitation echo; If exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value of this superrefraction clutter, to reject this superrefraction clutter; For the noise clutter point in the radar return on 1.5 ° of elevation angle faces, judge noise clutter point around whether have a precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value of this noise clutter point, to reject this noise clutter point;
5. discern the noise clutter point in the radar return on 0.5 ° of elevation angle face; And reject this noise clutter point, discern with reject 0.5 ° of elevation angle face on the process of noise clutter point in the radar return adopt with discern 1.5 ° of elevation angle faces in the radar return noise clutter point and reject the identical process of this noise clutter point;
6.-1, be outside 25~40db scope and the value of Vtbz is that dubiety sudden change zone is directly confirmed as in zone outside 20~40db scope with the value of Elbz 6. discern on the 0.5 ° of elevation angle face clutter group in the dubiety sudden change zone in the radar return, detailed process is:; 6.-2, adopt the regional filter method of sudden change, judge whether have doubtful clutter group in the dubiety zone, if, then continue execution, otherwise, the dubiety zone is not handled; 6.-3 the mean value Mdsw that, judges all the spectrum width numerical value in this doubtful clutter group whether less than the mean value Mdve of first setting value and all radial velocity numerical value whether less than second setting value; If the mean value Mdsw of all the spectrum width numerical value in this doubtful clutter group less than the mean value Mdve of first setting value and all radial velocity numerical value less than second setting value; Confirm that then it is clutter group that this doubtful clutter is rolled into a ball; Otherwise group does not handle to doubtful clutter;
7. discern the clutter that anomalous propagation causes in the radar return on 0.5 ° of elevation angle face; Detailed process is: judge on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone in whether have the superrefraction clutter of having rejected or do not have any clutter; And there is dubiety sudden change zone; If, then confirm on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone be by the zone of superrefraction clutter pollution, and mark should the zone be the clutter that anomalous propagation causes; Otherwise, on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone do not handle;
8. the clutter that the clutter in the radar return on 0.5 ° of elevation angle face is rolled into a ball and anomalous propagation causes is rejected processing; For the clutter group in the radar return on 0.5 ° of elevation angle face, judge clutter group around whether have precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value that this clutter is rolled into a ball, to reject this clutter group; The clutter that causes for the anomalous propagation in the radar return on 0.5 ° of elevation angle face; Judge clutter that anomalous propagation causes around whether have precipitation echo; If exist, then adopt numerical approximation filling mechanism to reject, otherwise; The reflectivity factor numerical value of the clutter that direct this anomalous propagation of rejecting causes is to reject the clutter that this anomalous propagation causes.
Described step 1. in the radar return process of dividing be: reflectivity factor numerical value is confirmed as clear air echo for the echo of [0dbz, (15dbz, 25dbz)], is [0dbz with reflectivity factor numerical value; (15dbz, 25dbz)] echo confirm as weak precipitation echo, be [(20dbz with reflectivity factor numerical value; 25dbz), (40dbz, 45dbz)] echo confirm as medium precipitation echo; With reflectivity factor numerical value be [(35dbz, 45dbz), (55dbz; Echo 60dbz)] is confirmed as the precipitation echo, reflectivity factor numerical value is confirmed as the dubiety echo greater than the difference of radar return reflectivity factor numerical value in vertical direction on 70dbz or 0.5 ° and 1.5 ° of two elevation angle face greater than the echo of 60dbz, wherein; Symbol " [] " expression numerical value is interval, and (a, b) expression value is that a is to any numerical value between the b.
The region area that described step is set in 2.-1 is more than or equal to 20m
2, described low value spectrum width numerical value is spectrum width numerical value less than 1 value.
Described step 2.-1 in the decision process in low value spectrum width zone be: judge reflection spectrum width numerical value and area greater than the mean value Mdsw of all the spectrum width numerical value in the zone of the region area of setting whether less than first setting value; If; Then confirm the zone to be low value spectrum width zone; Otherwise, confirm the zone to be non-low value spectrum width zone; Described first setting value is 2m/s; Described step 2.-2 in the decision process in low value radial velocity zone be: whether the mean value Mdve that judges all radial velocity numerical value that the zone of the reflection radial velocity numerical value that low value spectrum width zone is corresponding is interior less than second setting value; If; Then confirm the zone to be low value radial velocity zone; Otherwise, confirm the zone to be non-low value radial velocity zone; Described second setting value is 1m/s.
Described step 2.-3 in strong reflection rate factor numerical value be reflectivity factor numerical value greater than the reflectivity factor numerical value of precipitation.
The first threshold that described step is set in 3.-2 is 0.65~0.85.
Described numerical approximation filling machine is made as: write down the current clutter that needs processing; All echo points in the clutter that traversal has write down; Each echo point is carried out backfill handle,, define this echo point for working as the pre-echo point for the current echo point of just handling; With in the precipitation echo around the clutter that has write down with when the reflectivity factor numerical value of the immediate echo point in position of pre-echo point as reflectivity factor numerical value when pre-echo point, the backfill processing of pre-echo point is worked as in realization.
The detailed process of the sudden change zone filter method of described step in 6. is: a, judge an echo point in this dubiety zone reflectivity factor numerical value whether greater than 30dbz; If smaller or equal to 30dbz, do not handle, if greater than 30dbz; Then calculate on 1.5 ° of elevation angle faces around the corresponding echo point corresponding in the radar return mean value of the reflectivity factor numerical value of all echo points in 7 * 7 zones with this echo point position; Be designated as Mdz, whether the vertical change Vtbz that judges corresponding echo o'clock reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face again is greater than 10db, if less than 10db; Do not handle; If greater than 10db, judge that then whether the level of any echo reflectivity factor numerical value o'clock on 0.5 ° of elevation angle face in this correspondence echo point and eight echo points around it changes Elbz greater than second threshold value of setting, when the time smaller or equal to second threshold value of setting; Do not handle; When greater than second threshold value set, the next echo point in the search dubiety zone is handled, and all the echo points in the dubiety zone dispose; B, change Elbz by the level of all reflectivity factor numerical value on 0.5 ° of elevation angle face and constitute an echo group greater than the echo point of second threshold value of setting; Whether the scale-up factor on the sudden change the border whether number of judging the echo point in the echo group rolls into a ball less than the 3rd threshold value and this echo reaches the 4th threshold value; If confirm that then this echo group is doubtful clutter group.
Described second threshold value is 15db~25db, and described the 3rd threshold value is 50~100, and described the 4th threshold value is 0.5~0.7; The scale-up factor on the sudden change border that described echo is rolled into a ball is: whether the difference of the reflectivity factor numerical value of each echo point of 8 neighborhood echo points of a, the reflectivity factor numerical value of judging a borderline echo point that is arranged in this echo group and this echo point is all above the 5th threshold value; The 5th threshold value is 8db~12db, if confirm that then there is sudden change in this echo point; Define this echo point and be sudden change echo point; And continue to carry out, otherwise, confirm that there is not sudden change in this echo point; B, repeating step a continue to handle the borderline next echo point of this echo group, and borderline all echo points of rolling into a ball until this echo dispose; C, the number of adding up all sudden change echo points are calculated the ratio of the number of borderline all echo points that number and this echo of all sudden change echo points roll into a ball, and define this ratio and be the scale-up factor on the sudden change border that echo rolls into a ball.
Described step is rejected before the processing the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces and noise clutter point in 4.; Earlier judge whether superrefraction clutter or noise clutter point are in the 0.8km~1.2km scope outside medium precipitation echo or the precipitation echo border; If then superrefraction clutter or noise clutter point are not rejected processing, to carry out edge-protected to medium precipitation echo or precipitation echo; Otherwise, superrefraction clutter or noise clutter point are rejected processing; The clutter that described step causes the clutter in the radar return on 0.5 ° of elevation angle face group and anomalous propagation in is 8. rejected before the processing; Judge earlier whether the clutter that clutter group or anomalous propagation cause is in the 0.8km~1.2km scope outside medium precipitation echo or the precipitation echo border; If; Then the clutter that clutter is rolled into a ball or anomalous propagation causes is not rejected processing; Edge-protected so that medium precipitation echo or precipitation echo are carried out, otherwise, the clutter that clutter group or anomalous propagation cause is rejected processing.
Compared with prior art; The invention has the advantages that owing to mainly contain common ground clutter and AP clutter what the radar base data polluted; The reflectivity factor numerical value of these clutter aggregate performances face at the low elevation angle totally is discrete shape and distributes, and power mingles, and vertical change is big; 0.5 the radial velocity numerical value on ° elevation angle face be mixed with near zero numerical value, spectrum width numerical value is the discrete value of unordered strength distribution; And on 1.5 ° of elevation angle faces; The value that radial velocity numerical value is corresponding is substantially equal to zero; Spectrum width numerical value is less, thus the inventive method radar return is simply classified and basis to the reflectivity factor numerical analysis on, successively with reference to radial velocity numerical value and spectrum width numerical value; Discern on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on the clutter zone that exists; After correctly identifying the clutter zone, different clutters is carried out best clutter rejected strategy, promptly directly taked to reject strategy for the clutter under clear sky and the no precipitation echo situation of periphery; For being blended in the precipitation echo or the clutter adjacent with precipitation echo taked the numerical approximation mechanism of filling, can accessing like this and reject effect preferably.The inventive method has been passed through experimental demonstration; Divide distance branch rank that the inventive method is verified through 190 example classification; Checking is the result show, the inventive method general effect is excellent, is specially adapted within the 150km; Can guarantee data complete, effectively filter the non-meteorological clutter on the basis accurately; But outside 150km, because objective factor such as the Elbz of the convective echo of small scale and the bigger and remote speed of Vtbz is second-rate causes partly not ideal to AP clutter filter effect.Because the inventive method is relatively independent to the processing procedure of various clutters, so the inventive method has extraordinary expansibility, can be extended to 200km even farther to the effective range of filtering in theory.
Description of drawings
Fig. 1 is the process flow block diagram of the inventive method.
Embodiment
Embodiment describes in further detail the present invention below in conjunction with accompanying drawing.
What Radar Data adopted is the radar base data of polar format; The radar base data of radar return comprises reflectivity factor numerical value, radial velocity numerical value and spectrum width numerical value; And have the non-meteorological clutter in the radar return usually, need these non-meteorological clutters are filtered out, but when filtering non-meteorological clutter; May normal precipitation echo mistake be filtered; For avoiding mistake to filter precipitation echo, the present invention proposes the disposal route of non-meteorological clutter in a kind of radar return, the inventive method can obtain filter effect preferably.
The treatment scheme of the inventive method is as shown in Figure 1, and concrete steps are following:
1. radar return is carried out a simple classification; I.e. poor according to radar return reflectivity factor numerical value in vertical direction on the size of the reflectivity factor numerical value of radar return on 0.5 ° of elevation angle face and 0.5 ° and 1.5 ° of two elevation angle face; Radar return is divided into clear air echo, precipitation echo and dubiety echo, and precipitation echo comprises weak precipitation echo, medium precipitation echo and precipitation echo.
In this specific embodiment, the process that radar return is divided is: with reflectivity factor numerical value is [0dbz, (15dbz; Echo 25dbz)] is confirmed as clear air echo, is [0dbz, (15dbz with reflectivity factor numerical value; Echo 25dbz)] is confirmed as weak precipitation echo, with reflectivity factor numerical value be [(20dbz, 25dbz); (40dbz, 45dbz)] echo confirm as medium precipitation echo, be [(35dbz with reflectivity factor numerical value; 45dbz), (55dbz, 60dbz)] echo confirm as the precipitation echo; Reflectivity factor numerical value is confirmed as the dubiety echo greater than the difference of radar return reflectivity factor numerical value in vertical direction on 70dbz or 0.5 ° and 1.5 ° of two elevation angle face greater than the echo of 60dbz, and wherein, symbol " [] " expression numerical value is interval; (a, b) expression value be a to any numerical value between the b, this numerical value specifically depends on local forecast person's experience.In actual application, can confirm the concrete scope of the reflectivity factor numerical value of each echo according to the forecast experience of locality, like forecast experience according to the Ningbo of Zhejiang area; Can the span of the reflectivity factor numerical value of clear air echo be made as [0; 20], can the span of the reflectivity factor numerical value of weak precipitation echo be made as [0,20]; Can the span of the reflectivity factor numerical value of medium precipitation echo be made as [22; 43], can the span of the reflectivity factor numerical value of precipitation echo be made as [40,57].
Above-mentioned assorting process is not fully but that sorting result possibly exist thinks true echo to some clutter necessary in the entire process method.Based on experience, consider the precipitation echo that does not have in history greater than 70dbz, so can directly think the dubiety echo greater than the echo of 70dbz, the dubiety echo has very big probability should belong to clutter.
2. discern the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces; Detailed process is: 2.-1, on 1.5 ° of elevation angle faces in the radar return search area greater than the region area of setting and have the zone of low value spectrum width numerical value; Whether judgement should the zone be low value spectrum width zone, if then continue to carry out; Otherwise, do not handle; 2.-2, judge whether the zone of the reflection radial velocity numerical value that low value spectrum width zone is corresponding is low value radial velocity zone, if, then continue to carry out, otherwise, do not handle; 2.-3, judge whether have on 0.5 ° of elevation angle face of low value spectrum width zone below of 1.5 ° of elevation angle faces whether the zone of zone with strong reflection rate factor numerical value and should the zone corresponding reflection radial velocity numerical value is low value radial velocity zone; If having the zone of zone with strong reflection rate factor numerical value and should the zone corresponding reflection radial velocity numerical value is low value radial velocity zone; Confirm that then this zone with strong reflection rate factor numerical value is strong reflection rate zone; The area in this strong reflection rate zone is greater than the area in low value spectrum width zone, and definite low value spectrum width zone be doubtful superrefraction clutter, continuation execution then; Otherwise, do not handle; 2.-4, judge in the strong reflection rate zone whether have the dubiety echo; If exist; The level of then calculating reflectivity factor numerical value on 0.5 ° of elevation angle face change and 0.5 ° and 1.5 ° of two elevation angle face on the vertical change of reflectivity factor numerical value; The level of remembering reflectivity factor numerical value on 0.5 ° of elevation angle face is changed to Elbz, remembers that the vertical change of reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face is Vtbz, Elbz=Y (r) * | Z
I, j-Z
I, j-1|, Vtbz=Y (r) * (Z
Up-Z
Low), wherein, Y (r) representes the relevant weight of radar radial distance with the echo point position of working as pre-treatment, and its value is radar radial distance * 0.0025 of 1-when the echo point position of pre-treatment, and the unit of radar radial distance is km, Z
I, jRepresent that coordinate is (i, the reflectivity factor numerical value of echo point j), Z on 0.5 ° of elevation angle face
I, j-1Represent that coordinate is (i, the reflectivity factor numerical value of echo point j-1), Z on 0.5 ° of elevation angle face
UpThe reflectivity factor numerical value of representing the echo point on 1.5 ° of elevation angle faces, Z
LowRepresent on 0.5 ° of elevation angle face with 1.5 ° of elevation angle faces on the reflectivity factor numerical value of the corresponding echo point in echo point position; When calculating the vertical change Vtbz of reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face; Each echo point on each echo point on 0.5 ° of elevation angle face and the 1.5 ° of elevation angle faces needs corresponding on the upright position; When the value of Elbz is that the value of 25~40db and Vtbz is when being 20~40db; Doubtful superrefraction clutter is labeled as the superrefraction clutter, otherwise, doubtful superrefraction clutter is not handled.At this, the value in the span of Elbz and the span of Vtbz is different and definite according to distance.
In this specific embodiment, the region area of setting is more than or equal to 20m
2, in the actual application, directly the region area of setting is made as 20m
2, the region area of in fact setting is that the member if the region area that should set is too little, then analyzes not obvious to the experience gained of historical example statistics according to weather report.At this, value is confirmed as low value spectrum width numerical value less than 1 spectrum width numerical value.
In this specific embodiment; Because radial velocity numerical value that the superrefraction clutter is corresponding and spectrum width numerical value all are to approach zero numerical value; Therefore step 2.-1 in the decision process in low value spectrum width zone be: judge reflection spectrum width numerical value and area greater than the mean value Mdsw of all the spectrum width numerical value in the zone of the region area of setting whether less than first setting value, if then confirm the zone to be low value spectrum width zone; Otherwise, confirm the zone to be non-low value spectrum width zone; Step 2.-2 in the decision process in low value radial velocity zone be: whether the mean value Mdve that judges all radial velocity numerical value that the zone of the reflection radial velocity numerical value that low value spectrum width zone is corresponding is interior less than second setting value; If; Then confirm the zone to be low value radial velocity zone; Otherwise, confirm the zone to be non-low value radial velocity zone.Wherein, first setting value is 2m/s, and second setting value is 1m/s.The big value spectrum width numerical value that in low value spectrum width zone, possibly contain the sparse star spot shape of part, but major part all is to belong to low value spectrum width numerical value.The big value radial velocity numerical value that in low value radial velocity zone, possibly contain the sparse star spot shape of part, but major part all is to belong to low value radial velocity numerical value.
In this specific embodiment, strong reflection rate factor numerical value is the reflectivity factor numerical value greater than the reflectivity factor numerical value of precipitation.
3. discern the noise clutter point in the radar return on 1.5 ° of elevation angle faces, this mainly is to the object (ship, aircraft, flock of birds, group of insects etc.) isolated point that causes and the noise clutter point of spot that move.3.-1, calculate the coefficient of the isolated echo point in the radar return on the 1.5 ° of elevation angle faces detailed process of discerning the noise clutter point in the radar return on 1.5 ° of elevation angle faces is:, be designated as N
x,
Wherein, Mtotal=5 * 5, M for 5 * 5 region memories that are the center with this isolated echo point the number of non-NULL reflectivity factor numerical value; 3.-2, judgement factor N
xWhether less than the first threshold of setting, if then mark should isolated echo point be a noise clutter point, otherwise, do not handle isolating echo point.
In this specific embodiment; The first threshold of setting is 0.65~0.85, in the actual application, can first threshold be set according to local experience; Like the Ningbo of Zhejiang area; According to the desirable first threshold of Ningbo local experience is 0.75, and this threshold value is very effective for identification isolated point noise clutter point, but also can produce certain smoothing effect to the echo at knuckle place on the border of precipitation echo sometimes.
4. superrefraction clutter in the radar return on 1.5 ° of elevation angle faces and noise clutter point are rejected processing; Detailed process is: for the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces, judge the superrefraction clutter around whether have a precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value of this superrefraction clutter, to reject this superrefraction clutter; For the noise clutter point in the radar return on 1.5 ° of elevation angle faces, judge noise clutter point around whether have a precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value of this noise clutter point, to reject this noise clutter point.
5. discern the noise clutter point in the radar return on 0.5 ° of elevation angle face; And reject this noise clutter point, discern with reject 0.5 ° of elevation angle face on the process of noise clutter point in the radar return adopt with discern 1.5 ° of elevation angle faces in the radar return noise clutter point and reject the identical process of this noise clutter point.
6.-1, be outside 25~40db scope and the value of Vtbz is that dubiety sudden change zone is directly confirmed as in zone outside 20~40db scope with the value of Elbz 6. discern on the 0.5 ° of elevation angle face clutter group in the dubiety sudden change zone in the radar return, detailed process is:; 6.-2, adopt the regional filter method of sudden change, judge whether have doubtful clutter group in the dubiety zone, if, then continue execution, otherwise, the dubiety zone is not handled; 6.-3 the mean value Mdsw that, judges all the spectrum width numerical value in this doubtful clutter group whether less than the mean value Mdve of first setting value and all radial velocity numerical value whether less than second setting value; If the mean value Mdsw of all the spectrum width numerical value in this doubtful clutter group less than the mean value Mdve of first setting value and all radial velocity numerical value less than second setting value; Confirm that then it is clutter group that this doubtful clutter is rolled into a ball; Otherwise group does not handle to doubtful clutter.
In this specific embodiment, the detailed process of the regional filter method that suddenlys change is: a, judge an echo point in this dubiety zone reflectivity factor numerical value whether greater than 30dbz, if smaller or equal to 30dbz; Do not handle; If greater than 30dbz, then calculate on 1.5 ° of elevation angle faces around the corresponding echo point corresponding in the radar return mean value of the reflectivity factor numerical value of all echo points in 7 * 7 zones with this echo point position, be designated as Mdz; Whether the vertical change Vtbz that judges corresponding echo o'clock reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face again is greater than 10db; If less than 10db, do not handle, if greater than 10db; The level of then judging any echo reflectivity factor numerical value o'clock on 0.5 ° of elevation angle face in this correspondence echo point and eight echo points around it changes Elbz whether greater than second threshold value of setting; When smaller or equal to second threshold value set, do not handle, when the time greater than second threshold value set; Next echo point in the search dubiety zone is handled, and all the echo points in the dubiety zone dispose; B, change Elbz by the level of all reflectivity factor numerical value on 0.5 ° of elevation angle face and constitute an echo group greater than the echo point of second threshold value of setting; Whether the scale-up factor on the sudden change the border whether number of judging the echo point in the echo group rolls into a ball less than the 3rd threshold value and this echo reaches the 4th threshold value; If confirm that then this echo group is doubtful clutter group.
At this, second threshold value is 15db~25db, and the 3rd threshold value is 50~100; The 4th threshold value is 0.5~0.7, and in actual application, second threshold value, the 3rd threshold value and the 4th threshold value generally can be obtained according to the forecast experience of locality; Like the Ningbo of Zhejiang area, can get second threshold value according to the forecast experience in Ningbo is 20db, and the 3rd threshold value can get 75; The 4th threshold value can get 0.6; If in practical application, the poor effect that obtains when realizing the threshold value that use gets then can be finely tuned this threshold value as required again; If (be fine setting failure in the 15db~25db) like the scope of second threshold value, then need reanalyse the leak that whether exists in the existing experience on the logic analysis process at given threshold range.
At this; The proportionality coefficient on the sudden change border of echo group is: a; Whether the difference of reflectivity factor numerical value of each echo point of 8 neighborhood echo points of reflectivity factor numerical value and this echo point of judging the borderline echo point be arranged in this echo group is all above the 5th threshold value; The 5th threshold value is 8db~12db; In actual applications; The 5th threshold value can be confirmed according to local experience; Like the Ningbo of Zhejiang area; Can get the 5th threshold value according to the Ningbo local experience is 10db; If; Confirm that then there is sudden change in this echo point, define this echo point and be sudden change echo point, and continue to carry out; Otherwise, confirm that there is not sudden change in this echo point; B, repeating step a continue to handle the borderline next echo point of this echo group, and borderline all echo points of rolling into a ball until this echo dispose; C, the number of adding up all sudden change echo points are calculated the ratio of the number of borderline all echo points that number and this echo of all sudden change echo points roll into a ball, and define this ratio and be the proportionality coefficient on the sudden change border that echo rolls into a ball.
7. discern the clutter that anomalous propagation causes in the radar return on 0.5 ° of elevation angle face; Detailed process is: judge on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone in whether have the superrefraction clutter of having rejected or do not have any clutter; And there is dubiety sudden change zone; If, then confirm on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone be by the zone of superrefraction clutter pollution, and mark should the zone be the clutter that anomalous propagation causes; Otherwise, on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone do not handle;
8. the clutter that the clutter in the radar return on 0.5 ° of elevation angle face is rolled into a ball and anomalous propagation causes is rejected processing; For the clutter group in the radar return on 0.5 ° of elevation angle face, judge clutter group around whether have precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value that this clutter is rolled into a ball, to reject this clutter group; The clutter that causes for the anomalous propagation in the radar return on 0.5 ° of elevation angle face; Judge clutter that anomalous propagation causes around whether have precipitation echo; If exist, then adopt numerical approximation filling mechanism to reject, otherwise; The reflectivity factor numerical value of the clutter that direct this anomalous propagation of rejecting causes is to reject the clutter that this anomalous propagation causes.
The numerical approximation filling machine that uses in the above-mentioned process that various clutters are rejected is made as: write down the i.e. superrefraction clutter on 1.5 ° of elevation angle faces of the current clutter that needs to handle; Or the noise clutter point on 1.5 ° of elevation angle faces; Or the clutter group on 0.5 ° of elevation angle face; Or the clutter that causes of the anomalous propagation on 0.5 ° of elevation angle face, all the echo points in the clutter that traversal has write down carry out backfill to each echo point and handle; For the current echo point of just handling; Define this echo point for when the pre-echo point, with in the precipitation echo around the clutter that has write down with when the reflectivity factor numerical value of the immediate echo point in position of pre-echo point as reflectivity factor numerical value when pre-echo point, the backfill processing of pre-echo point is worked as in realization.
Reject the stage at clutter; For obtaining the effect of clutter rejecting preferably; Before the clutter zone is rejected, can be earlier to possibly be present on 0.5 °, 1.5 ° and the 2.4 ° of elevation angle faces the higher strong storm (the echo point existence of reflectivity factor numerical value greater than 18dbz all arranged on the vertical direction on these three layers of elevation angle faces) of corresponding development carry out edge-protected, in the methods of the invention; Only consider 0.5 ° and 1.5 ° of elevation angle faces; Detailed process is: step is rejected before the processing the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces and noise clutter point in 4., judges whether superrefraction clutter or noise clutter point are in the 1km scope outside medium precipitation echo or the precipitation echo border earlier, if; Then superrefraction clutter or noise clutter point are not rejected processing; Promptly medium precipitation echo or precipitation echo have been carried out edge-protected, otherwise, superrefraction clutter or noise clutter point are rejected processing; The clutter that step causes the clutter in the radar return on 0.5 ° of elevation angle face group and anomalous propagation in is 8. rejected before the processing; Judge earlier whether the clutter that clutter group or anomalous propagation cause is in the 1km scope outside medium precipitation echo or the precipitation echo border; If; Then the clutter that clutter is rolled into a ball or anomalous propagation causes is not rejected processing; Promptly medium precipitation echo or precipitation echo have been carried out edge-protected, otherwise, the clutter that clutter group or anomalous propagation cause is rejected processing.
The inventive method is carried out best clutter to different clutters and is rejected strategy; For the clutter or the clutter adjacent that are blended in the real precipitation echo with real precipitation echo; Utilizing in the precipitation echo reflectivity factor numerical value of putting corresponding echo point near each echo of clutter to carry out backfill after the identification handles; And for the clutter in clear sky or the no precipitation of the periphery zone, the method for directly rejecting is adopted in the identification back.
The inventive method has been considered the conservative property principle; Each threshold ratio that promptly is provided with is more wide in range; Like this clutter of very confirming is just rejected processing or carried out backfill, the clutter that can not handle is not rejected processing as far as possible, thereby can effectively avoid the strong convection precipitation of small scale to be filtered by mistake.
Claims (6)
1. the disposal route of non-meteorological clutter in the radar return, the radar base data of radar return comprises reflectivity factor numerical value, radial velocity numerical value and spectrum width numerical value, it is characterized in that this disposal route may further comprise the steps:
1. poor based on radar return reflectivity factor numerical value in vertical direction on the size of the reflectivity factor numerical value of radar return on 0.5 ° of elevation angle face and 0.5 ° and 1.5 ° of two elevation angle face; Radar return is divided into clear air echo, precipitation echo and dubiety echo, and precipitation echo comprises weak precipitation echo, medium precipitation echo and precipitation echo;
2. discern the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces; Detailed process is: 2.-1, on 1.5 ° of elevation angle faces in the radar return search area greater than the region area of setting and have the zone of low value spectrum width numerical value; Whether judgement should the zone be low value spectrum width zone, if then continue to carry out; Otherwise, do not handle; 2.-2, judge whether the zone of the reflection radial velocity numerical value that low value spectrum width zone is corresponding is low value radial velocity zone, if, then continue to carry out, otherwise, do not handle; 2.-3, judge whether have on 0.5 ° of elevation angle face of low value spectrum width zone below of 1.5 ° of elevation angle faces whether the zone of zone with strong reflection rate factor numerical value and should the zone corresponding reflection radial velocity numerical value is low value radial velocity zone; If having the zone of zone with strong reflection rate factor numerical value and should the zone corresponding reflection radial velocity numerical value is low value radial velocity zone; Confirm that then this zone with strong reflection rate factor numerical value is strong reflection rate zone; The area in this strong reflection rate zone is greater than the area in low value spectrum width zone, and definite low value spectrum width zone be doubtful superrefraction clutter, continuation execution then; Otherwise, do not handle; 2.-4, judge in the strong reflection rate zone whether have the dubiety echo; If exist; The level of then calculating reflectivity factor numerical value on 0.5 ° of elevation angle face change and 0.5 ° and 1.5 ° of two elevation angle face on the vertical change of reflectivity factor numerical value; The level of remembering reflectivity factor numerical value on 0.5 ° of elevation angle face is changed to Elbz, remembers that the vertical change of reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face is Vtbz, Elbz=Y (r) * | Z
I, j-Z
I, j-1|, Vtbz=Y (r) * (Z
Up-Z
Low), wherein, the weight that Y (r) expression is relevant with the radar radial distance, its value is 1-radar radial distance * 0.0025, the unit of radar radial distance is km, Z
I, jRepresent that coordinate is (i, the reflectivity factor numerical value of echo point j), Z on 0.5 ° of elevation angle face
I, j-1Represent that coordinate is (i, the reflectivity factor numerical value of echo point j-1), Z on 0.5 ° of elevation angle face
UpThe reflectivity factor numerical value of representing the echo point on 1.5 ° of elevation angle faces, Z
LowRepresent on 0.5 ° of elevation angle face with 1.5 ° of elevation angle faces on the reflectivity factor numerical value of the corresponding echo point in echo point position; When the value of Elbz is that the value of 25~40db and Vtbz is when being 20~40db; Doubtful superrefraction clutter is labeled as the superrefraction clutter; Otherwise, doubtful superrefraction clutter is not handled;
Described step 2.-1 in the decision process in low value spectrum width zone be: judge reflection spectrum width numerical value and area greater than the mean value Mdsw of all the spectrum width numerical value in the zone of the region area of setting whether less than first setting value; If; Then confirm the zone to be low value spectrum width zone; Otherwise, confirm the zone to be non-low value spectrum width zone; Described first setting value is 2m/s; Described step 2.-2 in the decision process in low value radial velocity zone be: whether the mean value Mdve that judges all radial velocity numerical value that the zone of the reflection radial velocity numerical value that low value spectrum width zone is corresponding is interior less than second setting value; If; Then confirm the zone to be low value radial velocity zone; Otherwise, confirm the zone to be non-low value radial velocity zone; Described second setting value is 1m/s;
3.-1, calculate the coefficient of the isolated echo point in the radar return on 1.5 ° of elevation angle faces 3. discern the noise clutter point in the radar return on 1.5 ° of elevation angle faces, detailed process is:, be designated as N
x,
Wherein, Mtotal=5 * 5, M for 5 * 5 region memories that are the center with this isolated echo point the number of non-NULL reflectivity factor numerical value; 3.-2, judgement factor N
xWhether less than the first threshold of setting, if then mark should isolated echo point be a noise clutter point, otherwise, do not handle isolating echo point;
4. superrefraction clutter in the radar return on 1.5 ° of elevation angle faces and noise clutter point are rejected processing; Detailed process is: for the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces; Judge the superrefraction clutter around whether have a precipitation echo; If exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value of this superrefraction clutter, to reject this superrefraction clutter; For the noise clutter point in the radar return on 1.5 ° of elevation angle faces, judge noise clutter point around whether have a precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value of this noise clutter point, to reject this noise clutter point;
Described numerical approximation filling machine is made as: write down the current clutter that needs processing; All echo points in the clutter that traversal has write down; Each echo point is carried out backfill handle,, define this echo point for working as the pre-echo point for the current echo point of just handling; With in the precipitation echo around the clutter that has write down with when the reflectivity factor numerical value of the immediate echo point in position of pre-echo point as reflectivity factor numerical value when pre-echo point, the backfill processing of pre-echo point is worked as in realization;
5. discern the noise clutter point in the radar return on 0.5 ° of elevation angle face; And reject this noise clutter point, discern with reject 0.5 ° of elevation angle face on the process of noise clutter point in the radar return adopt with discern 1.5 ° of elevation angle faces in the radar return noise clutter point and reject the identical process of this noise clutter point;
6.-1, be outside 25~40db scope and the value of Vtbz is that dubiety sudden change zone is directly confirmed as in zone outside 20~40db scope with the value of Elbz 6. discern on the 0.5 ° of elevation angle face clutter group in the dubiety sudden change zone in the radar return, detailed process is:; 6.-2, adopt the regional filter method of sudden change, judge whether have doubtful clutter group in the dubiety zone, if, then continue execution, otherwise, the dubiety zone is not handled; 6.-3 the mean value Mdsw that, judges all the spectrum width numerical value in this doubtful clutter group whether less than the mean value Mdve of first setting value and all radial velocity numerical value whether less than second setting value; If the mean value Mdsw of all the spectrum width numerical value in this doubtful clutter group less than the mean value Mdve of first setting value and all radial velocity numerical value less than second setting value; Confirm that then it is clutter group that this doubtful clutter is rolled into a ball; Otherwise group does not handle to doubtful clutter;
The detailed process of the sudden change zone filter method of described step in 6. is: a, judge an echo point in this dubiety zone reflectivity factor numerical value whether greater than 30dbz; If smaller or equal to 30dbz, do not handle, if greater than 30dbz; Then calculate on 1.5 ° of elevation angle faces around the corresponding echo point corresponding in the radar return mean value of the reflectivity factor numerical value of all echo points in 7 * 7 zones with this echo point position; Be designated as Mdz, whether the vertical change Vtbz that judges corresponding echo o'clock reflectivity factor numerical value on 0.5 ° and 1.5 ° of two elevation angle face again is greater than 10db, if less than 10db; Do not handle; If greater than 10db, judge that then whether the level of any echo reflectivity factor numerical value o'clock on 0.5 ° of elevation angle face in this correspondence echo point and eight echo points around it changes Elbz greater than second threshold value of setting, when the time smaller or equal to second threshold value of setting; Do not handle; When greater than second threshold value set, the next echo point in the search dubiety zone is handled, and all the echo points in the dubiety zone dispose; B, change Elbz by the level of all reflectivity factor numerical value on 0.5 ° of elevation angle face and constitute an echo group greater than the echo point of second threshold value of setting; Whether the scale-up factor on the sudden change the border whether number of judging the echo point in the echo group rolls into a ball less than the 3rd threshold value and this echo reaches the 4th threshold value; If confirm that then this echo group is doubtful clutter group;
7. discern the clutter that anomalous propagation causes in the radar return on 0.5 ° of elevation angle face; Detailed process is: judge on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone in whether have the superrefraction clutter of having rejected or do not have any clutter; And there is dubiety sudden change zone; If, then confirm on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone be by the zone of superrefraction clutter pollution, and mark should the zone be the clutter that anomalous propagation causes; Otherwise, on 1.5 ° of elevation angle faces with 0.5 ° of elevation angle face on corresponding zone, strong reflection rate zone do not handle;
8. the clutter that the clutter in the radar return on 0.5 ° of elevation angle face is rolled into a ball and anomalous propagation causes is rejected processing; For the clutter group in the radar return on 0.5 ° of elevation angle face, judge clutter group around whether have precipitation echo, if exist; Then adopt numerical approximation filling mechanism to reject; Otherwise, directly reject the reflectivity factor numerical value that this clutter is rolled into a ball, to reject this clutter group; The clutter that causes for the anomalous propagation in the radar return on 0.5 ° of elevation angle face; Judge clutter that anomalous propagation causes around whether have precipitation echo; If exist, then adopt numerical approximation filling mechanism to reject, otherwise; The reflectivity factor numerical value of the clutter that direct this anomalous propagation of rejecting causes is to reject the clutter that this anomalous propagation causes.
2. the disposal route of non-meteorological clutter in a kind of radar return according to claim 1 is characterized in that the region area of setting during described step 2.-1 is more than or equal to 20m
2, described low value spectrum width numerical value is spectrum width numerical value less than 1 value.
3. the disposal route of non-meteorological clutter in a kind of radar return according to claim 1 is characterized in that strong reflection rate factor numerical value was the reflectivity factor numerical value greater than the reflectivity factor numerical value of precipitation during described step 2.-3.
4. the disposal route of non-meteorological clutter in a kind of radar return according to claim 3 is characterized in that the first threshold of setting during described step 3.-2 is 0.65~0.85.
5. the disposal route of non-meteorological clutter in a kind of radar return according to claim 1; It is characterized in that described second threshold value is 15db~25db, described the 3rd threshold value is 50~100, and described the 4th threshold value is 0.5~0.7; The scale-up factor on the sudden change border that described echo is rolled into a ball is: whether the difference of the reflectivity factor numerical value of each echo point of 8 neighborhood echo points of a, the reflectivity factor numerical value of judging a borderline echo point that is arranged in this echo group and this echo point is all above the 5th threshold value; The 5th threshold value is 8db~12db, if confirm that then there is sudden change in this echo point; Define this echo point and be sudden change echo point; And continue to carry out, otherwise, confirm that there is not sudden change in this echo point; B, repeating step a continue to handle the borderline next echo point of this echo group, and borderline all echo points of rolling into a ball until this echo dispose; C, the number of adding up all sudden change echo points are calculated the ratio of the number of borderline all echo points that number and this echo of all sudden change echo points roll into a ball, and define this ratio and be the scale-up factor on the sudden change border that echo rolls into a ball.
6. the disposal route of non-meteorological clutter in a kind of radar return according to claim 1; It is characterized in that during described step 4. the superrefraction clutter in the radar return on 1.5 ° of elevation angle faces and noise clutter point are rejected before the processing; Earlier judge whether superrefraction clutter or noise clutter point are in the 0.8km~1.2km scope outside medium precipitation echo or the precipitation echo border; If then superrefraction clutter or noise clutter point are not rejected processing, to carry out edge-protected to medium precipitation echo or precipitation echo; Otherwise, superrefraction clutter or noise clutter point are rejected processing; The clutter that described step causes the clutter in the radar return on 0.5 ° of elevation angle face group and anomalous propagation in is 8. rejected before the processing; Judge earlier whether the clutter that clutter group or anomalous propagation cause is in the 0.8km~1.2km scope outside medium precipitation echo or the precipitation echo border; If; Then the clutter that clutter is rolled into a ball or anomalous propagation causes is not rejected processing; Edge-protected so that medium precipitation echo or precipitation echo are carried out, otherwise, the clutter that clutter group or anomalous propagation cause is rejected processing.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910153740A CN101706571B (en) | 2009-11-05 | 2009-11-05 | Method for processing non-meteorological noise in radar echoes |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910153740A CN101706571B (en) | 2009-11-05 | 2009-11-05 | Method for processing non-meteorological noise in radar echoes |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101706571A CN101706571A (en) | 2010-05-12 |
CN101706571B true CN101706571B (en) | 2012-08-29 |
Family
ID=42376809
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910153740A Expired - Fee Related CN101706571B (en) | 2009-11-05 | 2009-11-05 | Method for processing non-meteorological noise in radar echoes |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101706571B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732550A (en) * | 2018-08-01 | 2018-11-02 | 北京百度网讯科技有限公司 | Method and apparatus for predicting radar return |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102608608B (en) * | 2012-03-15 | 2013-06-26 | 天津大学 | Method for filtering super-refraction of meteorology |
CN103149552B (en) * | 2013-02-04 | 2015-11-18 | 南京信息工程大学 | A kind of Doppler radar radial velocity field move back blur method |
CN103499813B (en) * | 2013-10-11 | 2015-06-17 | 邓勇 | Method for distinguishing cloud and precipitation by radar precipitation factor |
CN103675788B (en) * | 2013-12-05 | 2015-12-30 | 中国科学院空间科学与应用研究中心 | The separation method of precipitation echo and backscattering from ocean surfaces in scatterometer echoed signal |
CN104007428B (en) * | 2014-05-20 | 2016-04-20 | 宁波市气象台 | A kind of method identifying the non-meteorological echo that running car causes in radar return |
CN105487076B (en) * | 2016-01-06 | 2018-09-11 | 北京无线电测量研究所 | A kind of millimeter wave cloud radar dense fog visibility inversion method and system |
CN106199604B (en) * | 2016-06-16 | 2019-03-12 | 杭州电子科技大学 | A kind of typhoon motion method for tracing based on correlation analysis |
CN106646411A (en) * | 2016-10-31 | 2017-05-10 | 北京市人工影响天气办公室 | Ground object detection method and device |
CN108665486A (en) * | 2018-03-31 | 2018-10-16 | 天津大学 | A kind of labeling method of strong convection monomer sample |
CN109031237B (en) * | 2018-08-09 | 2022-06-24 | 南京信息工程大学 | Method for rapidly acquiring radar combined reflectivity contour line data |
CN109459751B (en) * | 2018-08-27 | 2022-10-14 | 北京理工大学 | Migratory flight biological information monitoring method based on weather radar data |
CN109949615B (en) * | 2019-04-08 | 2021-03-12 | 中国长江电力股份有限公司 | Ship low-speed early warning detection system and method |
CN110501760A (en) * | 2019-07-29 | 2019-11-26 | 中国气象局广州热带海洋气象研究所(广东省气象科学研究所) | A kind of hail identification and nowcasting method based on weather radar |
CN110824479A (en) * | 2019-10-28 | 2020-02-21 | 兰州大方电子有限责任公司 | Radar data processing method for short-term forecasting |
CN111624606B (en) * | 2020-05-27 | 2022-06-21 | 哈尔滨工程大学 | Radar image rainfall identification method |
CN111638496A (en) * | 2020-06-08 | 2020-09-08 | 上海眼控科技股份有限公司 | Radar echo data processing method, computer device, and medium |
CN112666527B (en) * | 2020-12-03 | 2023-03-07 | 象辑科技股份有限公司 | Weather radar filtering method fusing live precipitation data |
CN114280572B (en) * | 2022-03-04 | 2022-08-09 | 浙江省气象台 | Single radar echo quality control method, system and terminal for removing signal interference clutter |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1774645A (en) * | 2003-04-17 | 2006-05-17 | 英国国防部 | Correction of troposphere induced errors in global positioning systems |
-
2009
- 2009-11-05 CN CN200910153740A patent/CN101706571B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1774645A (en) * | 2003-04-17 | 2006-05-17 | 英国国防部 | Correction of troposphere induced errors in global positioning systems |
Non-Patent Citations (1)
Title |
---|
何彩芬等.宁波非气象雷达回波的人工智能识别及滤波.《应用气象学报》.2007,第18卷(第6期),第856-864页. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732550A (en) * | 2018-08-01 | 2018-11-02 | 北京百度网讯科技有限公司 | Method and apparatus for predicting radar return |
Also Published As
Publication number | Publication date |
---|---|
CN101706571A (en) | 2010-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101706571B (en) | Method for processing non-meteorological noise in radar echoes | |
CN109871599B (en) | Unmanned vehicle operation safety index evaluation method | |
CN107766945B (en) | Method for calculating bearing capacity of urban road network | |
CN104240030B (en) | Track traffic network dynamic security risk evaluation method | |
CN102393299B (en) | Method for quantitatively calculating operational reliability of rolling bearing | |
CN105741551A (en) | Method for evaluating traffic safety of unsignalized intersection | |
CN105701862A (en) | Ground object key point extraction method based on point cloud | |
CN106296462B (en) | A kind of existing wind power plant roughness value based on double anemometer tower data determines method | |
CN112132371A (en) | Urban flood risk assessment method based on coupled entropy weight-fuzzy clustering algorithm | |
CN104318780B (en) | Consider the freeway incident detection method of meteorological factor, road alignment factor | |
CN102628944A (en) | Stratus cloud and convective cloud automatic recognition method based on Doppler radar data | |
CN106779320A (en) | A kind of gas pipeline damage from third-party methods of risk assessment based on fuzzy mathematics | |
CN112363168B (en) | Assimilation fusion method based on radar extrapolation and mode prediction | |
CN115953024A (en) | Urban inland inundation disaster early warning grade evaluation method | |
Marsh | The CONCAWE model for calculating the propagation of noise from open-air industrial plants | |
CN106872973A (en) | The method monitored using CINRAD/SA weather radars and track forest fires echo | |
CN105787959A (en) | Method for multi-agent network object tracking based on improved adaptive particle filtering | |
CN104268415A (en) | Travel congestion probability estimation method | |
CN105719510A (en) | Efficiency evaluation method of road traffic accident chain blocking system under IOV (Internet of Vehicles) environment | |
CN110459053A (en) | Highway Speed of Mixed Traffic Flow estimation method based on intelligent network connection vehicle information | |
CN109902971A (en) | A kind of road side security risk evaluations and stage division | |
CN104007428B (en) | A kind of method identifying the non-meteorological echo that running car causes in radar return | |
Xu et al. | Risk evaluation system of navigation security based on coupled wind and wave model: a case of study of Qiongzhou strait | |
CN109977569B (en) | Multi-factor fused MOS storm surge process disaster simulation method | |
Sharma et al. | Analysis of wind characteristics parameters with the application of lidar and mast |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB02 | Change of applicant information |
Address after: 315012 No. 118 meteorological Road, Zhejiang, Ningbo Applicant after: NINGBO METEOROLOGICAL SERVICE CENTER Address before: 315012 No. 118 meteorological Road, Zhejiang, Ningbo Applicant before: Ningbo Meteorological Information Center |
|
COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: NINGBO METEOROLOGICAL INFORMATION CENTER TO: NINGBO METEOROLOGICAL SERVICE CENTER |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120829 Termination date: 20211105 |