CN113065454A - High-altitude parabolic target identification and comparison method and device - Google Patents
High-altitude parabolic target identification and comparison method and device Download PDFInfo
- Publication number
- CN113065454A CN113065454A CN202110338987.8A CN202110338987A CN113065454A CN 113065454 A CN113065454 A CN 113065454A CN 202110338987 A CN202110338987 A CN 202110338987A CN 113065454 A CN113065454 A CN 113065454A
- Authority
- CN
- China
- Prior art keywords
- parabolic
- target
- images
- area
- image
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000009467 reduction Effects 0.000 claims abstract description 27
- 238000012549 training Methods 0.000 claims description 49
- 238000012545 processing Methods 0.000 claims description 34
- 238000013145 classification model Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000013135 deep learning Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 241000238631 Hexapoda Species 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 241000256626 Pterygota <winged insects> Species 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/50—Context or environment of the image
- G06V20/52—Surveillance or monitoring of activities, e.g. for recognising suspicious objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
- G06V10/443—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by matching or filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/56—Extraction of image or video features relating to colour
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30232—Surveillance
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a method and a device for identifying and comparing high-altitude parabolic targets, wherein the method comprises the steps of obtaining multi-frame images of a video to be identified, identifying two continuous frame images in the multi-frame images to obtain target points, removing noise points in the two continuous frame images by using a single-classification noise reduction model, removing non-parabolic objects in the two continuous frame images by using the single-classification non-parabolic model, carrying out area comparison and color comparison on the targets in the two continuous frame images from which the noise points and the non-parabolic objects are removed, and determining the parabolic targets in the two continuous frame images. Because the single classification noise reduction model and the single classification non-parabolic model are used for removing noise points and non-parabolic objects in the image, the interference given in the image can be reduced, and the accuracy of parabolic target identification is improved.
Description
Technical Field
The invention relates to the technical field of high-altitude parabolas, in particular to a method and a device for identifying and comparing high-altitude parabolas.
Background
High altitude parabolas are generally sporadic events, and in video detection, a great amount of detected objects are noise (such as noise of a camera and noise caused by light and shadow changes) and non-parabolic objects (such as flying birds, flying insects, fallen leaves and the like. Thus, in high altitude parabola detection scenarios, identification and validation of parabolas is critical to the problem. If the problem cannot be solved, a large number of false detections occur, and noise points and non-parabolic objects are identified as parabolas, so that false alarms are caused.
General processing methods of image processing include expansion, erosion, noise reduction, and the like. These processing methods are also mentioned in numerous patents for high altitude parabolic detection. However, for high altitude parabolic detection, the parabolic target is small in the picture when the parabolic target is located on a high floor (as shown in fig. 1). In this case, the erosion or noise reduction process is performed, which causes the parabola to disappear from the picture.
In addition, even if the above-described image processing method is used, a large noise point (for example, due to a change in light and shadow, fig. 2) cannot be removed. The probability of the continuous occurrence of the noise point is very small, but still larger than the probability of the occurrence of the true parabola, and if not processed, tens or even tens of false alarms are given every day.
For non-parabolic interferences such as birds and insects, the method using the above-mentioned general image processing is difficult to solve.
Disclosure of Invention
The embodiment of the invention provides a high-altitude parabolic target identification and comparison method and device, and solves the problem that noise and non-parabolic interference cannot be removed in image processing in the prior art.
In a first aspect, an embodiment of the present invention provides a method for identifying and comparing high-altitude parabolic targets, including:
acquiring multi-frame images of a video to be identified, and identifying two continuous frames of images in the multi-frame images to obtain a target point;
removing noise points in the two continuous frames of images by using a single classification noise reduction model;
removing the non-parabolic object in the two continuous frames of images by using a single-classification non-parabolic model;
and carrying out area comparison and color comparison on the targets in the two continuous frames of images without the noise points and the non-parabolic objects, and determining the parabolic targets in the two continuous frames of images.
In the technical scheme, the single-classification noise reduction model and the single-classification non-parabolic model are used for removing noise points and non-parabolic objects in the image, so that the interference given in the image can be reduced, and the accuracy of parabolic target identification is improved.
Optionally, the identifying two consecutive images of the multiple frames of images to obtain the target point includes:
carrying out graying processing on the multi-frame image;
and performing difference calculation on two continuous frames of images in the grayed multi-frame images to determine a target point.
Optionally, the single classification noise reduction model is determined according to the following steps:
acquiring a training sample marked with a noise point;
extracting features of each frame of image in the training sample marked with the noise points;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification noise reduction model.
Optionally, determining the single classification non-parabolic model according to the following steps includes:
acquiring a training sample marked with a non-parabolic object;
extracting features from each frame of image in the training sample marked with the non-parabolic object;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification non-parabolic model.
Optionally, the extracting features for each frame of image includes:
determining the area of the noise points or the non-parabolic areas in each frame of image according to the number of pixel points covered by the noise points or the non-parabolic connected domains;
determining the ratio of the area to the perimeter of each noise point or non-parabolic object in each frame of image;
determining the hue difference between each noise point or non-parabolic object and the background in each frame of image according to the color space;
determining the area of each noise point or non-parabolic object, the ratio of the area of each noise point or non-parabolic object to the perimeter and the hue difference of each noise point or non-parabolic object and the background in each frame image as the characteristics of each frame image.
Optionally, the area comparing the noise-removed point with the target in two consecutive frames of images of the non-parabolic object includes:
for the two continuous frames of images, if the area of the target in the previous frame of image is smaller than or equal to the area threshold, determining that the target in the previous frame of image and the target in the next frame of image are continuous parabolic targets when the ratio of the area of the target in the next frame of image to the area of the target in the previous frame of image is determined to be within a first area range;
if the area of the target in the previous frame image is larger than the area threshold, determining that the target in the previous frame image and the target in the next frame image are continuous parabolic targets when the ratio of the area of the target in the next frame image to the area of the target in the previous frame image is determined to be in a second area range;
the area threshold, the first area range, and the second area range are determined from parabolic experimental data.
Optionally, the color comparing the target in the two consecutive frames of images without the noise point and the non-parabolic object includes:
extracting the outline of the target by the frame difference of the gray level images aiming at the two continuous frames of images;
determining the outline of the target in the color map from the color map corresponding to the gray map according to the outline of the target in the gray map;
and calculating the picture similarity of the outlines of the targets in the two continuous frames of images, and if the picture similarity meets a similarity threshold, determining that the targets in the two continuous frames of images are continuous parabolic targets.
In a second aspect, an embodiment of the present invention provides an apparatus for identifying and comparing high-altitude parabolic targets, including:
the device comprises an acquisition unit, a processing unit and a processing unit, wherein the acquisition unit is used for acquiring multi-frame images of a video to be identified;
the processing unit is used for identifying two continuous frames of images in the multi-frame images to obtain a target point; removing noise points in the two continuous frames of images by using a single classification noise reduction model; removing the non-parabolic object in the two continuous frames of images by using a single-classification non-parabolic model; and carrying out area comparison and color comparison on the targets in the two continuous frames of images without the noise points and the non-parabolic objects, and determining the parabolic targets in the two continuous frames of images.
Optionally, the processing unit is specifically configured to:
carrying out graying processing on the multi-frame image;
and performing difference calculation on two continuous frames of images in the grayed multi-frame images to determine a target point.
Optionally, the processing unit is specifically configured to:
determining the single classification noise reduction model according to the following steps:
acquiring a training sample marked with a noise point;
extracting features of each frame of image in the training sample marked with the noise points;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification noise reduction model.
Optionally, the processing unit is specifically configured to:
determining the single classification non-parabolic model according to the following steps comprising:
acquiring a training sample marked with a non-parabolic object;
extracting features from each frame of image in the training sample marked with the non-parabolic object;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification non-parabolic model.
Optionally, the processing unit is specifically configured to:
determining the area of the noise points or the non-parabolic areas in each frame of image according to the number of pixel points covered by the noise points or the non-parabolic connected domains;
determining the ratio of the area to the perimeter of each noise point or non-parabolic object in each frame of image;
determining the hue difference between each noise point or non-parabolic object and the background in each frame of image according to the color space;
determining the area of each noise point or non-parabolic object, the ratio of the area of each noise point or non-parabolic object to the perimeter and the hue difference of each noise point or non-parabolic object and the background in each frame image as the characteristics of each frame image.
Optionally, the processing unit is specifically configured to:
for the two continuous frames of images, if the area of the target in the previous frame of image is smaller than or equal to the area threshold, determining that the target in the previous frame of image and the target in the next frame of image are continuous parabolic targets when the ratio of the area of the target in the next frame of image to the area of the target in the previous frame of image is determined to be within a first area range;
if the area of the target in the previous frame image is larger than the area threshold, determining that the target in the previous frame image and the target in the next frame image are continuous parabolic targets when the ratio of the area of the target in the next frame image to the area of the target in the previous frame image is determined to be in a second area range;
the area threshold, the first area range, and the second area range are determined from parabolic experimental data.
Optionally, the processing unit is specifically configured to:
extracting the outline of the target by the frame difference of the gray level images aiming at the two continuous frames of images;
determining the outline of the target in the color map from the color map corresponding to the gray map according to the outline of the target in the gray map;
and calculating the picture similarity of the outlines of the targets in the two continuous frames of images, and if the picture similarity meets a similarity threshold, determining that the targets in the two continuous frames of images are continuous parabolic targets.
In a third aspect, an embodiment of the present invention further provides a computing device, including:
a memory for storing program instructions;
and the processor is used for calling the program instructions stored in the memory and executing the high-altitude parabolic target identification and comparison method according to the obtained program.
In a fourth aspect, an embodiment of the present invention further provides a computer-readable non-volatile storage medium, which includes computer-readable instructions, and when the computer-readable instructions are read and executed by a computer, the computer is caused to execute the above-mentioned method for identifying and comparing high-altitude parabolic targets.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic diagram of a parabolic target according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating a noise point according to an embodiment of the present invention;
FIG. 3 is a diagram illustrating a system architecture according to an embodiment of the present invention;
FIG. 4 is a schematic flow chart illustrating a method for identifying and comparing high altitude parabolic targets according to an embodiment of the present invention;
FIG. 5 is a schematic flow chart illustrating a method for identifying and comparing high altitude parabolic targets according to an embodiment of the present invention;
FIG. 6 is a schematic flow chart of a single classification noise reduction model training according to an embodiment of the present invention;
FIG. 7 is a schematic flow chart of a single classification non-parabolic model training process according to an embodiment of the present invention;
FIG. 8 is a schematic flow chart of continuous target area comparison provided by an embodiment of the present invention;
FIG. 9 is a schematic flow chart of continuous target color distribution comparison provided by an embodiment of the present invention;
fig. 10 is a schematic structural diagram of an apparatus for identifying and comparing high-altitude parabolic targets according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the present invention will be described in further detail with reference to the accompanying drawings, and it is apparent that the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 3 is a system architecture according to an embodiment of the present invention. As shown in fig. 3, the system architecture may be a server 100, and the server 100 may include a processor 110, a communication interface 120, and a memory 130.
The communication interface 120 is used for communicating with other terminal devices, and transceiving information transmitted by the other terminal devices to implement communication.
The processor 110 is a control center of the server 100, connects various parts of the entire server 100 using various interfaces and lines, performs various functions of the server 100 and processes data by running or executing software programs and/or modules stored in the memory 130 and calling data stored in the memory 130. Alternatively, processor 110 may include one or more processing units.
The memory 130 may be used to store software programs and modules, and the processor 110 executes various functional applications and data processing by operating the software programs and modules stored in the memory 130. The memory 130 may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function, and the like; the storage data area may store data created according to a business process, and the like. Further, the memory 130 may include high speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other volatile solid state storage device.
It should be noted that the structure shown in fig. 3 is only an example, and the embodiment of the present invention is not limited thereto.
Based on the above description, fig. 4 shows in detail a flow of a method for high-altitude parabolic target recognition and comparison according to an embodiment of the present invention, where the flow may be executed by an apparatus for high-altitude parabolic target recognition and comparison.
As shown in fig. 4, the process specifically includes:
In the embodiment of the invention, the monitoring video can be obtained, and each frame image of the video is obtained. When the target point is identified, the multi-frame image is subjected to graying processing, and then the difference value of two continuous frames of images in the grayed multi-frame image is calculated to determine the target point.
And 402, removing noise points in the two continuous frames of images by using a single classification noise reduction model.
In embodiments of the present invention, noise is common in high-altitude parabolic video, which is often not noticeable if the video is viewed by the naked eye, but is easily detected if frame differencing is performed, especially when there is a change in light and shadow. In contrast, the parabolic events are sporadic events, the video of the actually occurring parabolic events cannot easily collect a large amount of data, and the collection of data by artificially and intentionally parabolic events is difficult to really cover most parabolic varieties. For this case, a single classification method can be used to learn the characteristics of the noise, even if it is not known what the parabola is, at least what the noise is. This allows the noise to be removed.
Single classification methods that may be used include isolated forests, single class support vector machines (one-class-SVMs), and the like. The method comprises the steps of collecting high-altitude parabolic prevention monitoring videos under various conditions, and extracting noise point data through a frame difference method. Model training is performed using the area size, shape, color, and the like of the noise point as its features.
The area of the noise point has various methods, such as calculating the area by using a circumscribed rectangle of the noise point. In order to distinguish noise from a true parabola, the number of pixels covered by a noise point connected domain is used for measurement. The shape feature is characterized by the ratio of area to perimeter. The shape of the noise point is generally irregular, and hollowness (such as noise caused by camera shaking) sometimes occurs; these two features are therefore good features for defining noise points.
The noise caused by the change in shade generally fluctuates in brightness without a change in color. That is, in the HSV color space (hue (H), saturation (S), and lightness (V)), the hue of a noise point is close to the hue of the background. Thus, in principle, noise can be eliminated by comparing the detected hue of the target with the hue at the position where it is normal. Therefore, the tonal difference between the noise point and the background can be used as the characteristic to participate in the single classification machine learning training.
The above mentioned features are mainly trained using conventional machine learning. Of course, a deep learning method (convolutional neural network) can also be used to perform single class learning on the noise. However, deep learning generally requires more computing resources and is more costly. If the cost is reduced, the model is required to be lightened. Consistent with the foregoing from a single classification concept.
Specifically, when a single classification noise reduction model is trained, a training sample marked with noise points needs to be obtained; then extracting features of each frame of image in the training sample marked with the noise points; and finally, inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining a single classification noise reduction model.
And 403, removing the non-parabolic object in the two continuous frames of images by using a single-classification non-parabolic model.
Based on the technical concept of the single-classification noise reduction model, the single-classification non-parabolic model can be trained, and a training sample marked with a non-parabolic model is obtained; extracting features of each frame of image in the training sample marked with the non-parabolic object; and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining a single classification non-parabolic model.
The process of extracting the features may specifically be: determining the area of the noise points or the non-parabolic areas in each frame of image according to the number of pixel points covered by the noise points or the non-parabolic connected domains; determining the ratio of the area to the perimeter of each noise point or non-parabolic object in each frame of image; determining the hue difference between each noise point or non-parabolic object and the background in each frame of image according to the color space; determining the area of each noise point or non-parabolic object, the ratio of the area of each noise point or non-parabolic object to the perimeter and the hue difference of each noise point or non-parabolic object and the background in each frame image as the characteristics of each frame image.
Under the condition that enough training samples of non-parabolic objects (flying birds, flying insects and the like) are collected, the noise can be removed in the same way, namely the characteristics of various non-parabolic objects are learned by using a single classification method.
In addition, the movement mode of the non-parabolic object is generally different from the parabolic falling mode, and in the case of failure of identification by using a single classification model, the non-parabolic object can be distinguished according to the characteristics of free falling body movement. See in particular the prior art high altitude parabolic starting point determination scheme.
And step 404, performing area comparison and color comparison on the targets in the two continuous frames of images without the noise points and the non-parabolic objects, and determining parabolic targets in the two continuous frames of images.
The area comparison mainly aims at two continuous frames of images, and if the area of the target in the previous frame of image is smaller than or equal to an area threshold value, when the ratio of the area of the target in the next frame of image to the area of the target in the previous frame of image is determined to be within a first area range, the target in the previous frame of image and the target in the next frame of image are determined to be continuous parabolic targets; and if the area of the target in the previous frame image is larger than the area threshold, determining that the target in the previous frame image and the target in the next frame image are continuous parabolic targets when the ratio of the area of the target in the next frame image to the area of the target in the previous frame image is determined to be in a second area range.
Wherein the area threshold, the first area range and the second area range are determined from parabolic experimental data.
When color comparison is carried out, extracting the outline of a target through the frame difference of the gray level images aiming at two continuous frames of images; determining the outline of the target in the color map from the color map corresponding to the gray map according to the outline of the target in the gray map; and calculating the image similarity of the outlines of the targets in the two continuous frames of images, and if the image similarity meets a similarity threshold, determining that the targets in the two continuous frames of images are continuous parabolic targets. Wherein, the similarity threshold value can be set according to experience.
In the specific implementation process, the targets detected at the two moments before and after need to be confirmed to be the same target, so that target comparison is needed. The comparison method comprises target area and color distribution.
First consider the change in parabolic area.
Normally, in the falling process of the parabola, the distance from the camera is closer and closer, and the area in the image is larger and larger. The target area at the previous time is atThe target area at the latter time is at+1Then a ist+1>at. However, when the parabola is small or is far from the camera, the detection error is large, and the area cannot be strictly restricted. It is reasonable to consider the area ratio of parabolic recognition at the front and back time instants (typically two pictures before and after) to be within a certain range. For example, the parameter is set to r (r)<1) Then the reasonable area ratio range is [ r,1/r ]](first area range). That is, the target area ratio of the current frame and the next frame is in the range, and the current frame and the next frame are regarded as the same target.
When the target is large enough (a)t>athre,athreArea threshold), the target area error has a small influence, and generally, a can be considered ast+1/at>1, and at+1/atK. K is a ratio calculated from image range (mapping conversion of pixel distance to true distance) without considering an error.
Actual length L represented by pixel points in horizontal directionxThe coordinate y, which depends on the vertical direction of the pixel, can be expressed as:
the coordinate y of the pixel point in the vertical direction of the plane and used for calculating the actual distance′Can be expressed as:
y′=γLx(y)……………………………………(2)
the three parameters in equations (1) and (2) can be obtained by multi-point calibration of the pixel length and the actual length in the picture.
For convenience of explanation, it is assumed that a rectangular object with an actual height H and a width W translates on the floor, and appears as a pixel with a height H on the camera picturepixWidth of wpix. Its center coordinate on the screen is (x, y). Then, one can obtain:
H=γLx(y-0.5hpix)-γLx(y+0.5hpix)…………………(4)
from (3) can be obtained:
solving the quadratic equation by (4) can obtain:
then, the area ratio at time t +1 and t is:
taking into account the presence of errors, at+1/atShould be in the range centered on K, i.e., [ K- ε K, K + ε K](second area range). ε is a positive number less than 1.
The three parameters (r, a)threEpsilon) need to be determined by parabolic experiments.
Next, the color distribution difference of the two frames of objects before and after the comparison is also needed. The target contour can be extracted through the frame difference of the gray-scale image, and then the target at the same position in the original color image is extracted. Next, in essence, comparison of picture similarity is performed, and various algorithms can be implemented, such as three-channel color histogram comparison, grayscale histogram comparison, and the like, and a deep learning method can also be used, and the model needs to be lightened in consideration of consumption of computing resources.
In order to better explain the embodiment of the present invention, the following describes the above-mentioned process of high-altitude parabolic target recognition in a specific implementation scenario, specifically as follows:
as shown in fig. 5, the specific steps are as follows:
And 502, performing difference between the two grayed frames to obtain a target point.
And 503, removing noise points by using a single classification noise reduction model.
As shown in fig. 6, the training process of the single-classification noise reduction model is specifically as follows:
And step 604, calculating the hue difference between the noise point and the background in the HSV color space.
And 606, testing on the test set to ensure that sufficient accuracy and recall rate are achieved, and obtaining the trained single-classification noise reduction model.
The specific process of the training process of the single-classification non-parabolic model shown in fig. 7 is as follows:
And step 702, calculating the non-parabolic area according to the number of pixel points covered by the non-parabolic connected domain.
And 505, performing area comparison of adjacent time points on the target subjected to noise removal and non-parabolic target, wherein the target meeting the area comparison condition is a parabolic target.
A specific process of continuous target area comparison may be as shown in fig. 8, and specifically includes:
The area of the target at the previous time is atThe area of the target at the subsequent time is at+1. If area atA is less than or equal tothre(area threshold value), then judge at+1/atWhether or not in [ r,1/r ]]In the range (first area range), 0<r<1. If so, the target is considered to be a continuous parabolic target.
If object atArea greater than athreUsing the formula (7) to obtain at+1/atThe standard ratio K of (A) is determined by considering the error coefficientt+1/atWhether or not [ K- ε K, K + ε K]Within a range (second area range). If so, the target is considered to be a continuous parabolic target.
And step 506, performing color comparison of adjacent time points on the target with the noise and the non-parabolic target removed, wherein the target meeting the color comparison condition is the parabolic target.
The specific process of continuous target color contrast may be as shown in fig. 9, specifically as follows:
and step 901, extracting a target contour through a frame difference of the gray-scale image.
And step 902, taking out the target at the same position in the original color image.
The target obtained through the filtering in step 507 from step 503 to step 506 is considered to be a parabolic target.
In the embodiment of the invention, a multi-frame image of a video to be identified is obtained, two continuous frames of images in the multi-frame image are identified to obtain a target point, a single-classification noise reduction model is used for removing noise points in the two continuous frames of images, a single-classification non-parabolic model is used for removing non-parabolic objects in the two continuous frames of images, the noise points and the targets in the two continuous frames of images without the parabolic objects are subjected to area comparison and color comparison, and parabolic targets in the two continuous frames of images are determined. Because the single classification noise reduction model and the single classification non-parabolic model are used for removing noise points and non-parabolic objects in the image, the interference given in the image can be reduced, and the accuracy of parabolic target identification is improved.
Based on the same technical concept, fig. 10 exemplarily shows the structure of an apparatus for high-altitude parabolic target recognition and comparison provided by an embodiment of the present invention, and the apparatus can perform a flow of high-altitude parabolic target recognition and comparison.
As shown in fig. 10, the apparatus specifically includes:
an acquisition unit 1001 configured to acquire a multi-frame image of a video to be identified;
the processing unit 1002 is configured to identify two consecutive frames of images in the multiple frames of images to obtain a target point; removing noise points in the two continuous frames of images by using a single classification noise reduction model; removing the non-parabolic object in the two continuous frames of images by using a single-classification non-parabolic model; and carrying out area comparison and color comparison on the targets in the two continuous frames of images without the noise points and the non-parabolic objects, and determining the parabolic targets in the two continuous frames of images.
Optionally, the processing unit 1002 is specifically configured to:
carrying out graying processing on the multi-frame image;
and performing difference calculation on two continuous frames of images in the grayed multi-frame images to determine a target point.
Optionally, the processing unit 1002 is specifically configured to:
determining the single classification noise reduction model according to the following steps:
acquiring a training sample marked with a noise point;
extracting features of each frame of image in the training sample marked with the noise points;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification noise reduction model.
Optionally, the processing unit 1002 is specifically configured to:
determining the single classification non-parabolic model according to the following steps comprising:
acquiring a training sample marked with a non-parabolic object;
extracting features from each frame of image in the training sample marked with the non-parabolic object;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification non-parabolic model.
Optionally, the processing unit 1002 is specifically configured to:
determining the area of the noise points or the non-parabolic areas in each frame of image according to the number of pixel points covered by the noise points or the non-parabolic connected domains;
determining the ratio of the area to the perimeter of each noise point or non-parabolic object in each frame of image;
determining the hue difference between each noise point or non-parabolic object and the background in each frame of image according to the color space;
determining the area of each noise point or non-parabolic object, the ratio of the area of each noise point or non-parabolic object to the perimeter and the hue difference of each noise point or non-parabolic object and the background in each frame image as the characteristics of each frame image.
Optionally, the processing unit 1002 is specifically configured to:
for the two continuous frames of images, if the area of the target in the previous frame of image is smaller than or equal to the area threshold, determining that the target in the previous frame of image and the target in the next frame of image are continuous parabolic targets when the ratio of the area of the target in the next frame of image to the area of the target in the previous frame of image is determined to be within a first area range;
if the area of the target in the previous frame image is larger than the area threshold, determining that the target in the previous frame image and the target in the next frame image are continuous parabolic targets when the ratio of the area of the target in the next frame image to the area of the target in the previous frame image is determined to be in a second area range;
the area threshold, the first area range, and the second area range are determined from parabolic experimental data.
Optionally, the processing unit 1002 is specifically configured to:
extracting the outline of the target by the frame difference of the gray level images aiming at the two continuous frames of images;
determining the outline of the target in the color map from the color map corresponding to the gray map according to the outline of the target in the gray map;
and calculating the picture similarity of the outlines of the targets in the two continuous frames of images, and if the picture similarity meets a similarity threshold, determining that the targets in the two continuous frames of images are continuous parabolic targets.
Based on the same technical concept, an embodiment of the present invention further provides a computing device, including:
a memory for storing program instructions;
and the processor is used for calling the program instructions stored in the memory and executing the high-altitude parabolic target identification and comparison method according to the obtained program.
Based on the same technical concept, embodiments of the present invention further provide a computer-readable non-volatile storage medium, which includes computer-readable instructions, and when the computer reads and executes the computer-readable instructions, the computer is caused to execute the above-mentioned method for identifying and comparing high-altitude parabolic targets.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including preferred embodiments and all such alterations and modifications as fall within the scope of the invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (10)
1. A method for identifying and comparing high-altitude parabolic targets is characterized by comprising the following steps:
acquiring multi-frame images of a video to be identified, and identifying two continuous frames of images in the multi-frame images to obtain a target point;
removing noise points in the two continuous frames of images by using a single classification noise reduction model;
removing the non-parabolic object in the two continuous frames of images by using a single-classification non-parabolic model;
and carrying out area comparison and color comparison on the targets in the two continuous frames of images without the noise points and the non-parabolic objects, and determining the parabolic targets in the two continuous frames of images.
2. The method of claim 1, wherein said identifying a target point for two consecutive images of said plurality of images comprises:
carrying out graying processing on the multi-frame image;
and performing difference calculation on two continuous frames of images in the grayed multi-frame images to determine a target point.
3. The method of claim 1, wherein the single classification noise reduction model is determined according to the following steps:
acquiring a training sample marked with a noise point;
extracting features of each frame of image in the training sample marked with the noise points;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification noise reduction model.
4. The method of claim 3, wherein determining the single classification non-parabolic model is based on the following steps comprising:
acquiring a training sample marked with a non-parabolic object;
extracting features from each frame of image in the training sample marked with the non-parabolic object;
and inputting the extracted features into a preset single classification model for training until a preset training target is met, and obtaining the single classification non-parabolic model.
5. The method of claim 4, wherein said extracting features for each frame of image comprises:
determining the area of the noise points or the non-parabolic areas in each frame of image according to the number of pixel points covered by the noise points or the non-parabolic connected domains;
determining the ratio of the area to the perimeter of each noise point or non-parabolic object in each frame of image;
determining the hue difference between each noise point or non-parabolic object and the background in each frame of image according to the color space;
determining the area of each noise point or non-parabolic object, the ratio of the area of each noise point or non-parabolic object to the perimeter and the hue difference of each noise point or non-parabolic object and the background in each frame image as the characteristics of each frame image.
6. The method of any one of claims 1 to 5, wherein the area comparing the noise point removed and the target in the two consecutive images of the non-parabolic object comprises:
for the two continuous frames of images, if the area of the target in the previous frame of image is smaller than or equal to the area threshold, determining that the target in the previous frame of image and the target in the next frame of image are continuous parabolic targets when the ratio of the area of the target in the next frame of image to the area of the target in the previous frame of image is determined to be within a first area range;
if the area of the target in the previous frame image is larger than the area threshold, determining that the target in the previous frame image and the target in the next frame image are continuous parabolic targets when the ratio of the area of the target in the next frame image to the area of the target in the previous frame image is determined to be in a second area range;
the area threshold, the first area range, and the second area range are determined from parabolic experimental data.
7. The method of any one of claims 1 to 5, wherein the color comparing the noise point removed and the target in the two consecutive images of the non-parabolic object comprises:
extracting the outline of the target by the frame difference of the gray level images aiming at the two continuous frames of images;
determining the outline of the target in the color map from the color map corresponding to the gray map according to the outline of the target in the gray map;
and calculating the picture similarity of the outlines of the targets in the two continuous frames of images, and if the picture similarity meets a similarity threshold, determining that the targets in the two continuous frames of images are continuous parabolic targets.
8. An apparatus for identifying and comparing high altitude parabolic target, comprising:
the device comprises an acquisition unit, a processing unit and a processing unit, wherein the acquisition unit is used for acquiring multi-frame images of a video to be identified;
the processing unit is used for identifying two continuous frames of images in the multi-frame images to obtain a target point; removing noise points in the two continuous frames of images by using a single classification noise reduction model; removing the non-parabolic object in the two continuous frames of images by using a single-classification non-parabolic model; and carrying out area comparison and color comparison on the targets in the two continuous frames of images without the noise points and the non-parabolic objects, and determining the parabolic targets in the two continuous frames of images.
9. A computing device, comprising:
a memory for storing program instructions;
a processor for calling program instructions stored in said memory to perform the method of any of claims 1 to 7 in accordance with the obtained program.
10. A computer-readable non-transitory storage medium including computer-readable instructions which, when read and executed by a computer, cause the computer to perform the method of any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110338987.8A CN113065454B (en) | 2021-03-30 | 2021-03-30 | High-altitude parabolic target identification and comparison method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110338987.8A CN113065454B (en) | 2021-03-30 | 2021-03-30 | High-altitude parabolic target identification and comparison method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113065454A true CN113065454A (en) | 2021-07-02 |
CN113065454B CN113065454B (en) | 2023-01-17 |
Family
ID=76564441
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110338987.8A Active CN113065454B (en) | 2021-03-30 | 2021-03-30 | High-altitude parabolic target identification and comparison method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113065454B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114998794A (en) * | 2022-05-31 | 2022-09-02 | 天翼爱音乐文化科技有限公司 | High-altitude parabolic recognition method, system, device and storage medium |
CN115294744A (en) * | 2022-07-29 | 2022-11-04 | 杭州海康威视数字技术股份有限公司 | Image display system, method, device and equipment |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060245002A1 (en) * | 2005-04-21 | 2006-11-02 | Koji Kita | Image processing method and program for restricting granular noise and module for implementing the method |
CN101493400A (en) * | 2008-01-25 | 2009-07-29 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic classification correcting method based on shape characteristic |
CN101795400A (en) * | 2010-03-16 | 2010-08-04 | 上海复控华龙微***技术有限公司 | Method for actively tracking and monitoring infants and realization system thereof |
CN102298709A (en) * | 2011-09-07 | 2011-12-28 | 江西财经大学 | Energy-saving intelligent identification digital signage fused with multiple characteristics in complicated environment |
CN103975342A (en) * | 2012-01-12 | 2014-08-06 | 柯法克斯公司 | Systems and methods for mobile image capture and processing |
CN109101944A (en) * | 2018-08-27 | 2018-12-28 | 四创科技有限公司 | A kind of real-time video monitoring algorithm identifying rubbish of jettisoninging into river |
CN110751044A (en) * | 2019-09-19 | 2020-02-04 | 杭州电子科技大学 | Urban noise identification method based on deep network migration characteristics and augmented self-coding |
CN111260693A (en) * | 2020-01-20 | 2020-06-09 | 北京中科晶上科技股份有限公司 | Detection method of high-altitude object throwing |
CN111768431A (en) * | 2020-06-28 | 2020-10-13 | 熵康(深圳)科技有限公司 | High-altitude parabolic moving target detection method, detection equipment and detection system |
CN111868260A (en) * | 2017-08-07 | 2020-10-30 | 约翰斯霍普金斯大学 | Methods and materials for assessing and treating cancer |
CN111931599A (en) * | 2020-07-20 | 2020-11-13 | 浙江大华技术股份有限公司 | High altitude parabolic detection method, equipment and storage medium |
CN112016414A (en) * | 2020-08-14 | 2020-12-01 | 熵康(深圳)科技有限公司 | Method and device for detecting high-altitude parabolic event and intelligent floor monitoring system |
CN112037266A (en) * | 2020-11-05 | 2020-12-04 | 北京软通智慧城市科技有限公司 | Falling object identification method and device, terminal equipment and storage medium |
CN112213707A (en) * | 2020-10-14 | 2021-01-12 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Method for predicting action range of ship-borne pulse radar in evaporation waveguide environment |
CN112257557A (en) * | 2020-10-20 | 2021-01-22 | 中国电子科技集团公司第五十八研究所 | High-altitude parabolic detection and identification method and system based on machine vision |
-
2021
- 2021-03-30 CN CN202110338987.8A patent/CN113065454B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060245002A1 (en) * | 2005-04-21 | 2006-11-02 | Koji Kita | Image processing method and program for restricting granular noise and module for implementing the method |
CN101493400A (en) * | 2008-01-25 | 2009-07-29 | 深圳迈瑞生物医疗电子股份有限公司 | Automatic classification correcting method based on shape characteristic |
CN101795400A (en) * | 2010-03-16 | 2010-08-04 | 上海复控华龙微***技术有限公司 | Method for actively tracking and monitoring infants and realization system thereof |
CN102298709A (en) * | 2011-09-07 | 2011-12-28 | 江西财经大学 | Energy-saving intelligent identification digital signage fused with multiple characteristics in complicated environment |
CN103975342A (en) * | 2012-01-12 | 2014-08-06 | 柯法克斯公司 | Systems and methods for mobile image capture and processing |
CN111868260A (en) * | 2017-08-07 | 2020-10-30 | 约翰斯霍普金斯大学 | Methods and materials for assessing and treating cancer |
CN109101944A (en) * | 2018-08-27 | 2018-12-28 | 四创科技有限公司 | A kind of real-time video monitoring algorithm identifying rubbish of jettisoninging into river |
CN110751044A (en) * | 2019-09-19 | 2020-02-04 | 杭州电子科技大学 | Urban noise identification method based on deep network migration characteristics and augmented self-coding |
CN111260693A (en) * | 2020-01-20 | 2020-06-09 | 北京中科晶上科技股份有限公司 | Detection method of high-altitude object throwing |
CN111768431A (en) * | 2020-06-28 | 2020-10-13 | 熵康(深圳)科技有限公司 | High-altitude parabolic moving target detection method, detection equipment and detection system |
CN111931599A (en) * | 2020-07-20 | 2020-11-13 | 浙江大华技术股份有限公司 | High altitude parabolic detection method, equipment and storage medium |
CN112016414A (en) * | 2020-08-14 | 2020-12-01 | 熵康(深圳)科技有限公司 | Method and device for detecting high-altitude parabolic event and intelligent floor monitoring system |
CN112213707A (en) * | 2020-10-14 | 2021-01-12 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Method for predicting action range of ship-borne pulse radar in evaporation waveguide environment |
CN112257557A (en) * | 2020-10-20 | 2021-01-22 | 中国电子科技集团公司第五十八研究所 | High-altitude parabolic detection and identification method and system based on machine vision |
CN112037266A (en) * | 2020-11-05 | 2020-12-04 | 北京软通智慧城市科技有限公司 | Falling object identification method and device, terminal equipment and storage medium |
Non-Patent Citations (2)
Title |
---|
MOHAMMAD ALI SAGHAFI 等: "Review of person re-identification techniques", 《IET COMPUTER VISION》 * |
胡燕: "基于分类器融合的SAR图像目标识别", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114998794A (en) * | 2022-05-31 | 2022-09-02 | 天翼爱音乐文化科技有限公司 | High-altitude parabolic recognition method, system, device and storage medium |
CN115294744A (en) * | 2022-07-29 | 2022-11-04 | 杭州海康威视数字技术股份有限公司 | Image display system, method, device and equipment |
CN115294744B (en) * | 2022-07-29 | 2024-03-22 | 杭州海康威视数字技术股份有限公司 | Image display system, method, device and equipment |
Also Published As
Publication number | Publication date |
---|---|
CN113065454B (en) | 2023-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9230175B2 (en) | System and method for motion detection in a surveillance video | |
EP3246874B1 (en) | Method and apparatus for updating a background model used for background subtraction of an image | |
CN109636771B (en) | Flight target detection method and system based on image processing | |
CN109284673B (en) | Object tracking method and device, electronic equipment and storage medium | |
CN106910204B (en) | A kind of method and system to the automatic Tracking Recognition of sea ship | |
CN113065454B (en) | High-altitude parabolic target identification and comparison method and device | |
CN111027546B (en) | Character segmentation method, device and computer readable storage medium | |
CN111325769A (en) | Target object detection method and device | |
WO2006115676A2 (en) | Video ghost detection by outline | |
CN111223129A (en) | Detection method, detection device, monitoring equipment and computer readable storage medium | |
CN116485779B (en) | Adaptive wafer defect detection method and device, electronic equipment and storage medium | |
CN111274964B (en) | Detection method for analyzing water surface pollutants based on visual saliency of unmanned aerial vehicle | |
CN117456371B (en) | Group string hot spot detection method, device, equipment and medium | |
CN112261402B (en) | Image detection method and system and camera shielding monitoring method and system | |
CN113962900A (en) | Method, device, equipment and medium for detecting infrared dim target under complex background | |
CN113743378A (en) | Fire monitoring method and device based on video | |
CN111402185B (en) | Image detection method and device | |
CN110334703B (en) | Ship detection and identification method in day and night image | |
CN109977965B (en) | Method and device for determining detection target in remote sensing airport image | |
CN115661475A (en) | Image foreign matter identification method, device, equipment and storage medium | |
CN115909151A (en) | Method for identifying serial number of motion container under complex working condition | |
CN112801072B (en) | Elevator non-flat-layer door opening fault recognition device and method based on computer vision | |
CN112307916A (en) | Alarm monitoring method based on visible light camera | |
CN111985334A (en) | Gun detection method, system and device and storage medium | |
CN117690096B (en) | Contact net safety inspection system adapting to different scenes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |