CN106056637A - Method for processing smog area test data of powder charging flow field - Google Patents
Method for processing smog area test data of powder charging flow field Download PDFInfo
- Publication number
- CN106056637A CN106056637A CN201610423474.6A CN201610423474A CN106056637A CN 106056637 A CN106056637 A CN 106056637A CN 201610423474 A CN201610423474 A CN 201610423474A CN 106056637 A CN106056637 A CN 106056637A
- Authority
- CN
- China
- Prior art keywords
- smog
- matrix
- image sequence
- flow field
- transmitance
- 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
Classifications
-
- 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/10048—Infrared image
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a method for processing smog area test data of a powder charging flow field. The method comprises the main steps of: firstly, converting video data into a gray-scale image sequence; secondly, converting the image sequence into a smog transmittance matrix with the same resolution; thirdly, collecting statics on the quantity and the distribution of pixels of the transmittance within a threshold range; and fourthly, calculating the corresponding smog region area via a measuring scale. The method combines the smog transmittance characteristic of the flow field with the spatial characteristic to obtain such parameters as smog region distribution and area within the specific smog transmittance range in the flow field. The method can effectively promote the development of an assessment technology of influence of propellant plume smog on the guidance performance and a clean propellant smog performance assessment technology, and is suitable for processing video data of engines, gas generators and barrel weapons in flow field smog area test.
Description
Technical field
The present invention relates to the processing method of a kind of smog area test data, it is adaptable to powder charge flow field smog area is surveyed
Have a try and test the processing method of middle flow field video data.
Background technology
Low characteristic signal propellant and clean propellant powder are important directions of China's gunpowder development.Powder charge stream
Field smog has a negative impact for the stealthy of armament systems and guidance performance.On the one hand, when electromotor and weapons with pipes work
Powder burning produces the smog in flow field can expose position and the movement locus of armament systems.On the other hand, propellant loading plumage
Stream flow field smog for the guidance signal passed through have decay, scatter, the multiple effect such as absorption, cause the guidance of weapon to be believed
Number intensity reduces, and even loses efficacy time serious.Smog during gun propellant combustion produces flow field in weapons with pipes can be to weapon operation hands
Or the secondary of video guidance system aims at and interferes.The smog annoyance level to guidance signal, not only the concentration with smog has
Close, also and the area of smog, thickness and guidance signal pass the path-dependent of smog.
At present, mostly the domestic method of testing to powder charge flow field smog is to measure propellant loading or weapons with pipes test
During the smog transmitance in flow field, clear and definite method of testing be there is no for powder charge flow field smog area.At other field
In fire prevention, the test of smoke agent smog characteristic, generally morphological characteristic and gray threshold according to smog are to video or figure
As the smoke region in data is extracted, and calculate smoke region area.Smog district determined by this data processing method
In territory, smog is indefinite for the interference effect of armament systems guidance signal, it is impossible to Efficient Characterization flow field smog is to guidance
The spatial dimension of signal specific effect degree (i.e. the optical transmittance of smog is within the specific limits).
Summary of the invention
In order to characterize associating between powder charge flow field smog transmitance and smog area, determine that in flow field, smog passes through
Rate smoke region area within the specific limits, the invention provides the place of a kind of powder charge flow field smog area test data
Reason method.This data processing method is applicable to the stream of the weapons with pipes (such as rifle, big gun etc.) of propellant loading and gun propellant charge
The video data of field smog area testing experiment processes.
The field layout of powder charge flow field smog area testing experiment, optics background board lays respectively at imageing sensor
Both sides, powder charge flow field.Optics background board is greyish white alternate striped, its edge and sides aligned parallel in image, and in process of the test
Keep static.Imageing sensor can comprise visible light wave range imageing sensor, infrared band imageing sensor, laser wavelength image
Sensor.Before test, the scale of length utilizing meter ruler that imageing sensor is obtained background board image is demarcated, it is thus achieved that test
The image scale K of video.Imageing sensor triggers before powder charge is lighted a fire, and records flow field before and after powder charge works
Video data.In the video data of test gained, the initial image frame of powder charge flow field smog area test video data
Number is no less than 50 frames.
In data processing method provided by the present invention, mainly comprise the steps that
(1) video data in powder charge flow field is processed, be converted to grayscale image sequence.
(2) grayscale image sequence is converted to the smog transmitance matrix of equal resolution.First, grayscale image sequence is taken
The average synthesis benchmark image S of front 50 two field picture gray scalesM×N;Secondly, by each two field picture B in grayscale image sequenceM×NWith base
The gray scale of quasi-image corresponding pixel points makees ratio, obtains gray level ratio matrix τk(i, j), then, at gray level ratio matrix τk(i,
J), in, formula (1) is utilized to eliminate the veiling glare impact on gray level ratio, it is thus achieved that smog transmitance matrix τ 'k(i, j);Formula
τ ' in formulak(i, j) represents the transmitance distribution matrix of kth frame image in grayscale image sequence, and wherein k represents image
Position in image sequence, i represents the position of data point horizontal direction in the picture, and j represents that data point is the most vertical
The position in direction;Gs(i j) represents benchmark image SM×NAt point (i, j) gray value at place;τk(i j) represents grayscale image sequence
The gray level ratio matrix of middle kth frame image, w represents the pixel that in image sequence, wall scroll background striped occupies in the horizontal direction
Number.
(3) the pixel quantity in threshold range of smog transmitance in every two field picture and distribution are added up.
(4) formula (2) is utilized to calculate the smoke region area S in every two field picturek, and draw smog area and change over
Curve.Formula
Sk=K2Nk (2)
S in formulakRepresenting the smoke region area of kth frame image in image sequence, unit is square metre;K is image sequence
Length ratio coefficient, unit is rice/pixel;NkIt is transmitance pixel quantity in threshold range in k-th frame image, single
Position is pixel.
(5) with the smog area-time curve of mark threshold value as final data processed result.
As smog transmitance matrix τ ' in step (3)k(i, when j) noise is bigger, in statistics smog transmitance at threshold value model
Before enclosing interior pixel quantity, can be to the smog transmitance matrix τ ' of each two field picturek(i, j) carries out linear filtering, sets
The neighborhood of filtering is m × n (m, n are odd number), the most filtered smog transmitance matrix τ "k(i j) is represented by
Threshold value in step (3) is a bounded interval, can be (τa, τb)、[τa, τb]、[τa, τb) or (τa, τb] form,
Wherein τa≤τb;τa, τbCan determine according to actual needs.
In step (3) during statistics, pixel reunion quantity N can be increasedrAs additional qualifications, the rate that will transmit through is in threshold value
In the range of reunion quantity < NrPixel get rid of outside scope of statistics.
The light transmission rate of smog, compared with the data processing method of other smog areas, is tied by the present invention with spatial character
Process altogether, it is possible to obtain there is the smoke region smog area parameters of specific smog transmitance scope.Gunpowder is made to fill
The sign of medicine plume smog space performance more has actual application value.
Accompanying drawing explanation
The field layout figure of Fig. 1 powder charge flow field smog area testing experiment, 1-optics background board, 2-visible images
Sensor, 3-infrared image sensor, 4-laser image sensor, 5-portable computer, 6-flow field regions.
The initial pictures of Fig. 2 powder charge flow field smog area test video data
The flow field smog area-time curve of Fig. 3 propellant loading
Detailed description of the invention
Below by embodiment, the present invention is further explained explanation.
The video file obtained by certain propellant loading flow field smog area testing experiment and the length of video image
Scale K=0.02 rice/pixel, ask transmitance in the range of [20%, 60%], the smoke region area of pixel reunion number >=4.
Such as Fig. 1, carrying out the field layout of powder charge flow field smog area testing experiment, optics background board 1 passes with image
Sensor lays respectively at powder charge flow field regions 6 both sides, and imageing sensor is connected with portable computer 5 respectively.Optics background
Plate is greyish white alternate striped, its edge and sides aligned parallel in image, and keeps static in process of the test.Imageing sensor comprises can
See optical image sensor 2, infrared image sensor 3, laser image sensor 4.Powder charge flow field smog area test video
Initial pictures (such as Fig. 2) frame number of data is no less than 50 frames.
The step that test data process is as follows:
(1) it is grayscale image sequence by smog field video conversion;Image sequence is carried out cutting, remove optics background board it
Outer image section.
(2) grayscale image sequence is converted to the smog transmitance matrix τ ' of equal resolutionk(i,j);First, gray scale is taken
The average synthesis benchmark image S of 50 two field picture gray scales before image sequence400×90;Secondly, by each frame figure in grayscale image sequence
As B400×90Make ratio with the gray scale of benchmark image corresponding pixel points, obtain gray level ratio matrix τk(i, j), then, at gray scale ratio
Value matrix τk(i, j) in, utilize formula (1) eliminate veiling glare to gray scale ratio
The impact of value, it is thus achieved that smog transmitance matrix τ 'k(i, j).At formula
In formula: M=400, the horizontal resolution of expression image:
W=20, is the pixel count that occupies in the horizontal direction of wall scroll background striped.
(3) formula (4) the smog transmitance matrix τ ' to each two field picture is utilizedk(i j) carries out linear filtering, formula
(4) it is represented by
(4) smog transmitance in every two field picture is added up in the range of [20%, 60%] and pixel reunion quantity >=4
Pixel quantity Nk。
(6) formula (2) is utilized to calculate the smoke region area S in every two field picturek, and draw smog area and change over
Curve.Formula
Sk=K2Nk (2)
In formula: K=0.02 rice/pixel
(7) with smog area-time curve (as shown in Figure 3) of mark threshold value as final data processed result.
Claims (4)
1. the processing method of powder charge flow field smog area test data, it is characterised in that step is as follows:
(1) video data in powder charge flow field is processed, be converted to grayscale image sequence;
(2) grayscale image sequence is converted to the smog transmitance matrix of equal resolution;First, 50 are taken before grayscale image sequence
The average synthesis benchmark image S of two field picture gray scaleM×N;Secondly, by each two field picture B in grayscale image sequenceM×NWith reference map
As the gray scale of corresponding pixel points makees ratio, obtain gray level ratio matrix τk(i, j), then, at gray level ratio matrix τk(i, j) in,
Formula (1) is utilized to eliminate the veiling glare impact on gray level ratio, it is thus achieved that smog transmitance matrix τ 'k(i, j);Formula
τ ' in formulak(i, j) represents the transmitance distribution matrix of kth frame image in grayscale image sequence, and wherein k represents that image is at figure
As the position in sequence, i represents the position of data point horizontal direction in the picture, and j represents data point vertical direction in the picture
Position;Gs(i j) represents benchmark image SM×NAt point (i, j) gray value at place;τk(i j) represents kth in grayscale image sequence
The gray level ratio matrix of two field picture, w represents the pixel count that in image sequence, wall scroll background striped occupies in the horizontal direction;
(3) the pixel quantity in threshold range of smog transmitance in every two field picture and distribution are added up;
(4) formula (2) is utilized to calculate the smoke region area S in every two field picturek, and draw the time dependent song of smog area
Line;Formula
Sk=K2Nk (2)
S in formulakRepresenting the smoke region area of kth frame image in image sequence, unit is square metre;K is the length of image sequence
Proportionality coefficient, unit is rice/pixel;NkBeing transmitance pixel quantity in threshold range in k-th frame image, unit is picture
Element;
(5) with the smog area-time curve of mark threshold value as final data processed result.
The processing method of data the most according to claim 1, it is characterised in that as smog transmitance matrix τ ' in step (3)k
(i, when j) noise is bigger, before statistics smog transmitance pixel quantity in threshold range, the cigarette to each two field picture
Mist transmitance matrix τ 'k(i, j) carries out linear filtering, sets the neighborhood of filtering as m × n (m, n are odd number), the most filtered
Smog transmitance matrix τ "k(i j) is expressed as
。
The processing method of data the most according to claim 1, it is characterised in that the threshold value in step (3) is one battery limit (BL)
Between, for (τa, τb)、[τa, τb]、[τa, τb) or (τa, τb] form, wherein τa≤τb;τa, τbDetermine according to actual needs.
The most according to claim 1, the processing method of data, it is characterised in that in step (3) during statistical pixel point, increase picture
Vegetarian refreshments reunion quantity NrAs additional qualifications, the rate that will transmit through reunion quantity < N in threshold rangerPixel get rid of exist
Outside scope of statistics.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610423474.6A CN106056637B (en) | 2016-06-15 | 2016-06-15 | A kind of processing method of powder charge flow field smog area test data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610423474.6A CN106056637B (en) | 2016-06-15 | 2016-06-15 | A kind of processing method of powder charge flow field smog area test data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106056637A true CN106056637A (en) | 2016-10-26 |
CN106056637B CN106056637B (en) | 2018-07-31 |
Family
ID=57168151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610423474.6A Active CN106056637B (en) | 2016-06-15 | 2016-06-15 | A kind of processing method of powder charge flow field smog area test data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106056637B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108174113A (en) * | 2016-12-06 | 2018-06-15 | 佳能株式会社 | Image processing apparatus |
CN110021011A (en) * | 2019-03-15 | 2019-07-16 | 横店集团东磁有限公司 | A kind of method that detection chip draws glue surface product |
CN110458825A (en) * | 2019-08-08 | 2019-11-15 | 西安近代化学研究所 | A kind of flow field smog transmitance distribution tester and data processing method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6650280B2 (en) * | 2000-12-08 | 2003-11-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Measurement system and method |
CN104809711A (en) * | 2015-05-14 | 2015-07-29 | 西安近代化学研究所 | Method for processing video image of plume smog of solid propellant |
CN104849241A (en) * | 2015-05-14 | 2015-08-19 | 西安近代化学研究所 | Calibration method of system for testing light shielding capability of smoke of propellants |
-
2016
- 2016-06-15 CN CN201610423474.6A patent/CN106056637B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6650280B2 (en) * | 2000-12-08 | 2003-11-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Measurement system and method |
CN104809711A (en) * | 2015-05-14 | 2015-07-29 | 西安近代化学研究所 | Method for processing video image of plume smog of solid propellant |
CN104849241A (en) * | 2015-05-14 | 2015-08-19 | 西安近代化学研究所 | Calibration method of system for testing light shielding capability of smoke of propellants |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108174113A (en) * | 2016-12-06 | 2018-06-15 | 佳能株式会社 | Image processing apparatus |
CN108174113B (en) * | 2016-12-06 | 2021-03-19 | 佳能株式会社 | Image processing apparatus |
CN110021011A (en) * | 2019-03-15 | 2019-07-16 | 横店集团东磁有限公司 | A kind of method that detection chip draws glue surface product |
CN110458825A (en) * | 2019-08-08 | 2019-11-15 | 西安近代化学研究所 | A kind of flow field smog transmitance distribution tester and data processing method |
CN110458825B (en) * | 2019-08-08 | 2022-05-24 | 西安近代化学研究所 | Flow field smoke transmittance distribution testing device and data processing method |
Also Published As
Publication number | Publication date |
---|---|
CN106056637B (en) | 2018-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106056637A (en) | Method for processing smog area test data of powder charging flow field | |
CN105842205A (en) | Muzzle smoke concentration test device and method | |
CN110110654B (en) | Amplitude inversion method and device for descending ocean isolated waves | |
CN203745204U (en) | Visual testing device for spraying and combustion dynamic conditions in diesel engine cylinder | |
CN103940714A (en) | Imitated artificial haze monitoring system and method | |
CN103792090A (en) | Diesel engine in-cylinder spraying and combustion visualization testing apparatus and testing method | |
CN103384895A (en) | Fast image enhancement and three-dimensional depth calculation | |
Shawal et al. | High-speed imaging of early flame growth in spark-ignited engines using different imaging systems via endoscopic and full optical access | |
CN104976960A (en) | Raindrop physical property observation method and device | |
CN104634784B (en) | atmospheric visibility monitoring method and device | |
Aptowicz et al. | Decomposition of atmospheric aerosol phase function by particle size and asphericity from measurements of single particle optical scattering patterns | |
Mann et al. | Experimental study of relative, turbulent diffusion | |
CN103617635A (en) | Transient flame detection method based on image processing | |
CN111666925B (en) | Training method and device for face recognition model | |
CN106769732A (en) | A kind of rectilinear diesel vehicle smoke intensity detection method | |
Goschütz et al. | Endoscopic imaging of early flame propagation in a near-production engine | |
CN103942786B (en) | The self adaptation block objects detection method of unmanned plane visible ray and infrared image | |
CN103185706A (en) | Laser measurement method and device for light obscuration of inorganization-emission particulate smoke plume | |
Barnie et al. | Quantification of gas and solid emissions during Strombolian explosions using simultaneous sulphur dioxide and infrared camera observations | |
CN115690211A (en) | Air explosion point three-dimensional coordinate detection device and measurement method | |
CN103761762B (en) | A kind of Simulations of Infrared Image method based on OpenGL tinter | |
CN104849241A (en) | Calibration method of system for testing light shielding capability of smoke of propellants | |
CN104809711B (en) | Solid propellant plume smog method of video image processing | |
CN204555826U (en) | A kind of electronic sighting device with transparent display device | |
CN205067360U (en) | Gaseous formation of image telemetry unit of sulfur dioxide |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |