CN109308366A - Halation emulation mode in low-light level television imaging - Google Patents

Halation emulation mode in low-light level television imaging Download PDF

Info

Publication number
CN109308366A
CN109308366A CN201710673245.4A CN201710673245A CN109308366A CN 109308366 A CN109308366 A CN 109308366A CN 201710673245 A CN201710673245 A CN 201710673245A CN 109308366 A CN109308366 A CN 109308366A
Authority
CN
China
Prior art keywords
low
indicate
microchannel plate
fluorescent screen
electron
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.)
Pending
Application number
CN201710673245.4A
Other languages
Chinese (zh)
Inventor
颜馨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201710673245.4A priority Critical patent/CN109308366A/en
Publication of CN109308366A publication Critical patent/CN109308366A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Image-Pickup Tubes, Image-Amplification Tubes, And Storage Tubes (AREA)

Abstract

The invention discloses the halation emulation mode in a kind of imaging of low-light level television, mainly solve the problems, such as that traditional experiment research method cost height and effect are emerald green low.Implementation step is: (1) on the three-dimensional scenic emulation platform based on OGRE, importing the threedimensional model of intense light source and low-light level television image device, generate three-dimensional low-light scene;(2) in three-dimensional low-light scene, according to the night television system internal electron characteristics of motion, photoelectricity transformation principle and microchannel plate scattering theory, statistics reaches the amount of electrons of fluorescent screen;(3) using the amount of electrons for reaching fluorescent screen as foundation, by system voltage principle of signal conversion and grey level quantization principle, the simulation of halation intensity profile is realized.Precision of the present invention is high, adapts to wide, strong real-time, is able to achieve the accurate simulation that intense light source is imaged in night television system.

Description

Halation emulation mode in low-light level television imaging
Technical field
Halation emulation side the invention belongs to computer simulation technique field, in particular in a kind of low-light level television imaging Method can be used for the low-light scene objects identification field under intense light source.
Background technique
Work as intense light source, such as when street lamp, car light, flare are appeared in night television system visual field, system exports image It will form halation around the upper luminous point.The strong light in part on halation can cover surrounding weak signal, directly affect the imaging point of system Resolution and signal-to-noise ratio;And under the constant same observation condition of the factors such as dimension of light source, observed range, different model imaging system The size of halation is different on system output image.Therefore, the quantitatively characterizing and its emulation mode for studying halation are under intense light source The raising etc. of the identification of target, image resolution ratio and signal-to-noise ratio has important practical value.
In recent years, the research in terms of domestic and foreign scholars have been substantially carried out following two to the halation in low-light level television imaging: (1) it by replacing micro-channel tubes with hollow pipe, is tested under different light sources, analyzing influence system exports halation in image Principal element;(2) by changing light source, change the conditions such as the parameter of image intensifier and tested, research influences halation size Factor.In those references, it is all characterized without the accurate quantification of halation in further research system output image, to cause micro- Photoelectricity regard imaging simulation output image accuracy is low and the small problem of the light source scope of application.
Summary of the invention
It is an object of the invention to be directed to the deficiency of above-mentioned prior art, propose that the halation in a kind of imaging of low-light level television is imitative True method improves the accuracy of low-light level television imaging simulation output image, expands micro- to provide the quantitatively characterizing of halation in image The scope of application of photoelectricity view imaging simulation light source.
Realizing the technical solution of the object of the invention is: according to low-light level television imaging system photoelectricity transformation principle, micro- logical Guidance tape scattering theory come analyze by intense light source generate electron motion to microchannel plate the characteristics of motion;It is fixed according to the conservation of energy Rule calculates electronics under the different motion state of microchannel plate surface, and the offset distance of impacting electron and each location of pixels receive Electron number;In conjunction with system voltage principle of signal conversion and grey level quantization principle, calculates these electronics and generated on output image Gray value, obtain the space annular spread of these gray values, as halation image.Implementation step includes the following:
(1) intense light source and optical lens, photocathode, electron lens, microchannel plate, fluorescent screen are inputted these constitutes low-lights The calibrating parameters of television imaging device, using the threedimensional model of 3Dmax Software Create intense light source and low-light level television image device, and The threedimensional model is imported in the three-dimensional scenic simulated program based on OGRE, low-light three-dimensional scenic is generated;
(2) according to the illuminance E of Facing material reflectivity ρ and bias light in low-light three-dimensional scenic and intense light source Brightness L, the brightness L of calculating optical lens surfacea(i, j), wherein the two dimension that (i, j) is optical lens surface face element is sat Mark;
(3) according to the brightness L of optical lens surfacea(i, j) calculates photocathode surface by light by photoelectric conversion According to the electron number n (i, j) of generation;
(4) inclined when calculating photocathode surface electronic by electron lens arrival microchannel plate according to law of conservation of energy Move distance s1(i, j);
(5) according to microchannel plate scattering theory, the electronics for reaching microchannel plate is calculated after microchannel plate Multiple Scattering The offset distance s (i, j) of generation;According to the Two dimensional Distribution of image on fluorescent screen vegetarian refreshments, calculates the offset distance and correspond to the every of fluorescent screen A pixel position (x, y), and count each received electron number n of pixel position institutee(x, y), wherein (x, y) is fluorescent screen The two-dimensional coordinate of surface all pixels point;
(6) according to the received electron number n of each pixel position of corresponding fluorescent screen institutee(x, y) passes through the electricity of fluorescent screen Light conversion calculates the gray value Gray (x, y) of all pixels point on fluorescent screen, these gray values are sat in the two dimension of fluorescence screen surfaces Mark central point is formed about a kind of annular distribution, and this annular distribution is halation image, and is exported.
Compared with prior art, the present invention having the following obvious advantages:
1. the present invention has fully considered that electronics reaches the different motion state of microchannel plate, light is realized by voltage quantization Dizzy emulation, it is thus possible to accurately emulate the halation phenomenon in low-light level television imaging;
2. the present invention is that low-light three-dimensional scenic is generated in the environment of demarcating intense light source and low-light level television image device parameter, Therefore it is directed to different intense light sources, can realize that its is corresponding by changing calibrating parameters to match low-light level television image-forming condition Emulation, expands the scope of application of low-light level television imaging simulation light source.
Invention Detailed description of the invention
Fig. 1 is general flow chart of the invention;
Fig. 2 is movement schematic diagram of step 4 electronics between electron lens and microchannel plate in the present invention.
Specific embodiment
Referring to Fig.1, the detailed implementation process of the halation emulation mode in low-light level television imaging of the present invention is as follows:
Step 1. generates low-light three-dimensional scenic.
It inputs intense light source and these composition low-lights of optical lens, photocathode, electron lens, microchannel plate, fluorescent screen is electric Depending on the calibrating parameters of image device, using the threedimensional model of 3Dmax Software Create intense light source and low-light level television image device, and will The threedimensional model imports in the three-dimensional scenic simulated program based on OGRE, generates low-light three-dimensional scenic.
The brightness of step 2. calculating optical lens surface.
(2a) calculates intense light source at optical lens surface face element according to light measurement formula according to the brightness L of intense light source Illuminance Ei(i, j)
Ei (i, j)=LAcos θ icos θ j12]] >
Wherein (i, j) is the two-dimensional coordinate of optical lens surface face element, and A indicates intense light source area, θiIndicate that low-light is three-dimensional The angle of surface face element normal in scene and intense light source between the incident ray at optical lens surface face element, θjIndicate strong The angle of light source surface panel method line and intense light source between the incident ray at optical lens surface face element, l indicate intense light source The distance between optical lens surface;
(2b) is according to the illuminance E and intense light source of Facing material reflectivity ρ, bias light in low-light three-dimensional scenic in optics Illuminance E at lens surface1, pass through the brightness L of reflected light law calculating optical lens surfacea
La=ρ (Ei+ E)/π,
(2c) combines (2a) and (2b), the brightness L of calculating optical lens surfacea(i, j)
La (i, j)=ρ (LAcos θ icos θ jl2+E)/π]] >
Step 3. calculates the electron number that photocathode surface generates.
According to the brightness L of optical lens surfacea(i, j), according to photoelectric conversion formula calculate photocathode surface by The electron number n (i, j) that illumination generates:
N (i, j)=π 4 (Dff) 2La (i, j) τ ·Sk ,]] >
Wherein DfIndicating System Optics effective aperture, f indicates optical system focal length, and τ indicates optical lens transmitance, SkIndicate photocathode sensitivity.
Step 4. calculates offset distance when electronics reaches microchannel plate.
Referring to Fig. 2, this step is implemented as follows:
(4a), according to law of conservation of energy, when calculating photocathode surface electronic by electron lens arrival microchannel plate Speed V1The angle theta of (i, j) and the speed and electron lens normal1(i, j):
V1 (i, j)=V02+2eume]] >
θ 1 (i, j)=arctanV0sin θ 0 (V0cos θ 0) 2+2eume ,]] >
Wherein V0、θ0Initial velocity and initial angle that photocathode surface generates electronics are respectively indicated, e indicates single electricity The quantity of electric charge of son, meIndicate the quality of Single Electron, u indicates the operating voltage between electron lens and microchannel plate;
(4b) reaches speed V when microchannel plate according to electronics1The angle theta of (i, j) and the speed and electron lens normal1 (i, j) calculates offset distance s when it reaches microchannel plate by law of conservation of energy1(i, j):
S1 (i, j)=(1 (i, j)-V0cos θ 0 of V1 (i, j) cos θ) meV0sin θ 0Deu ,]] >
Wherein, D indicates the distance between electron lens and microchannel plate.
Step 5. counts the received electron number of each pixel position institute.
(5a) reaches the speed and the angle theta of electron lens normal when microchannel plate according to electronics1(i, j), by micro- logical Guidance tape scattering theory calculates electronics and enters the offset distance s that primary collision occurs for microchannel plate2:
S2=2meV02+2eumesin θ 1 (i, j) eu ,]] >
Wherein meIndicate the quality of Single Electron, V0Indicate that photocathode surface generates the initial velocity of electronics, e indicates single The quantity of electric charge of a electronics, u indicate the operating voltage between electron lens and microchannel plate;
(5b) reaches offset distance s when microchannel plate according to electronics1(i, j), primary collision occurs into microchannel plate Offset distance s2With collision frequency N, the offset distance s (i, j) that electronics generates after microchannel plate Multiple Scattering is calculated:
S (i, j)=s1(i, j)+Ns2,
(5c) is by s in (4b)1S in (i, j) and (5a)2Expression formula substitute into the expression formula of (5b), calculate electronics by micro- logical The offset distance s (i, j) generated after guidance tape Multiple Scattering:
S (i, j)=(V1 (i, j) cos θ (i, j)-V0cos θ 0) meV0sin θ 0Deu+2NmeV02+2eumesin θ 1 (i, J) eu ,]] >
Wherein, V0、θ0Respectively indicate initial velocity and initial angle that photocathode surface generates electronics, meIndicate single The quality of electronics, e indicate that the quantity of electric charge of Single Electron, u indicate the operating voltage between electron lens and microchannel plate, and D indicates electricity Distance between sub-lens and microchannel plate, N indicate the number that electronics and the non-perforated wall of microchannel plate collide, and value is nature Number.
(5d) calculates the electronics process for reaching microchannel plate by following formula according to the Two dimensional Distribution of image on fluorescent screen vegetarian refreshments Each pixel position (x, y) of the corresponding fluorescent screen of the offset distance s (i, j) generated after microchannel plate Multiple Scattering:
| (xdw) 2+ (ydh) 2-s (i, j) |<δ ,]]>
Wherein, dwIndicate the width of fluorescent screen single pixel, dhIndicate the height of fluorescent screen single pixel, δ=min (dw, dh), s (i, j) indicates offset distance;
(5e) counts each received electron number n of pixel position institutee(x, y).
Step 6. calculates the gray value of all pixels point on fluorescent screen.
(6a) is according to the received electron number n of each pixel position institute for corresponding to fluorescent screene(x, y) passes through voltage signal Conversion formula calculates the quantization voltage V (x, y) of each pixel:
V (x, y)=π 4 τ e (Defe) 2R&CenterDot;As&CenterDot;Gv&CenterDot;v&CenterDot; G&CenterDot;k&CenterDot;Ne (x, y) &CenterDot;E/a2 ,]] >
Wherein τeIndicate the transmitance of coupled lens, DeIndicate coupled lens effective aperture, feIndicate coupled lens focal length, R indicates detector response rate, AsIndicate effectively photosensitive elemental area, GvIndicate that vision signal amplification factor, v indicate electron lens Operating voltage, G indicate the gain of microchannel plate, and k indicates fluorescent screen luminous efficacy, and a indicates electron-optical amplification factor, ne(x, Y) indicate that the received electron number of each location of pixels institute, e indicate the quantity of electric charge of Single Electron;
(6b) calculates the gray value Gray (x, y) of each pixel according to the quantization voltage V (x, y) of each pixel:
Gray (x, y)=255VMax-VMinV (x, y) ,]] >
Wherein, VMaxAnd VMinRespectively indicate the maximum value and minimum value of the quantization voltage of all pixels point;
The gray value of (6c) all pixels point is formed about a kind of cyclic annular minute in the two-dimensional coordinate central point of fluorescence screen surfaces Cloth, this annular distribution are the halation image exported.
Above description is only example of the present invention, does not constitute any limitation of the invention.Obviously for this It, all may be without departing substantially from the principle of the invention, structure after having understood the content of present invention and principle for the professional in field In the case of, various modifications and variations in form and details are carried out, but these modifications and variations based on inventive concept are still Within the scope of the claims of the present invention.

Claims (6)

1. the halation emulation mode in a kind of low-light level television imaging, includes the following steps:
(1) intense light source and optical lens, photocathode, electron lens, microchannel plate, fluorescent screen are inputted these constitutes low-light level televisions The calibrating parameters of image device, using the threedimensional model of 3Dmax Software Create intense light source and low-light level television image device, and should Threedimensional model imports in the three-dimensional scenic simulated program based on OGRE, generates low-light three-dimensional scenic;
(2) according to the illuminance E of Facing material reflectivity ρ and bias light in low-light three-dimensional scenic and the light of intense light source Spend L, the brightness L of calculating optical lens surfacea(i, j), wherein (i, j) is the two-dimensional coordinate of optical lens surface face element;
(3) according to the brightness L of optical lens surfacea(i, j) calculates photocathode surface by photoelectric conversion and is produced by illumination Raw electron number n (i, j);
(4) according to law of conservation of energy, calculate photocathode surface electronic by electron lens reach microchannel plate hour offset away from From s1(i, j)
(4a) calculates speed when photocathode surface electronic reaches microchannel plate by electron lens according to law of conservation of energy V1The angle theta of (i, j) and the speed and electron lens normal1(i, j) is calculated by following formula:
V1 (i, j)=V02+2eume]] >
θ 1 (i, j)=arctanV0sin θ 0 (V0cos θ 0) 2+2eume]] >
Wherein V0、θ0Indicate that photocathode surface generates initial velocity, the initial angle of electronics, e indicates the charge of Single Electron Amount, m0Indicate the quality of Single Electron, u indicates the operating voltage between electron lens and microchannel plate;
(4b) reaches speed V when microchannel plate according to electronics1The angle theta of (i, j) and the speed and electron lens normal1(i, J), offset distance s when it reaches microchannel plate is calculated by law of conservation of energyi(i, j)
S1 (i, j)=(1 (i, j)-V0cos θ 0 of V1 (i, j) cos θ) meV0sin θ 0Deu]] >
Wherein D indicates distance between electron lens and microchannel plate;
(5) according to microchannel plate scattering theory, the electronics for calculating arrival microchannel plate generates after microchannel plate Multiple Scattering Offset distance s (i, j);According to the Two dimensional Distribution of image on fluorescent screen vegetarian refreshments, the corresponding fluorescent screen of offset distance s (i, j) is calculated Each pixel position (x, y), and count each received electron number n of pixel position institutee(x, y), wherein (x, y) is fluorescence The two-dimensional coordinate of screen surfaces all pixels point;
(6) according to the received electron number n of each pixel position of corresponding fluorescent screen institutee(x, y) is turned by the electric light of fluorescent screen The gray value Gray (x, y) for calculating all pixels point on fluorescent screen is changed, these gray values are in the two-dimensional coordinate of fluorescence screen surfaces Heart point is formed about a kind of annular distribution, and this annular distribution is halation image, and is exported.
2. the halation emulation mode in a kind of low-light level television imaging according to claim 1, wherein meter described in step (2) Calculate the brightness L of optical lens surfacea(i, j) is calculated by following formula:
La (i, j)=ρ (LAcos θ icos θ j12+E)/π]] >
Wherein A indicates intense light source area, θiSurface face element normal and intense light source in expression low-light three-dimensional scenic is in optical lens table The angle between incident ray at the face element of face, θjIndicate intense light source surface face element normal and intense light source in optical lens surface face The angle between incident ray at member, 1 indicates the distance between intense light source and optical lens surface.
3. the halation emulation mode in a kind of low-light level television imaging according to claim 1, wherein leads to described in step (3) It crosses photoelectric conversion and calculates the electron number n (i, j) that photocathode surface is generated by illumination, calculation formula is as follows:
N (i, j)=π 4 (Dff) 2La (i, j) τ &CenterDot;Sk]]>
Wherein DfIndicate System Optics effective aperture, f indicates optical system focal length, and τ indicates optical lens transmitance, SkTable Show photocathode sensitivity.
4. the halation emulation mode in a kind of low-light level television imaging according to claim 1, wherein root described in step (5) According to microchannel plate scattering theory, calculates and reach the offset distance that the electronics of microchannel plate generates after microchannel plate Multiple Scattering S (i, j), is calculated as follows:
S (i, j)=(1 (i, j)-V0cos θ 0 of V1 (i, j) cos θ) meV0sin θ 0Deu+2NmeV02+2eumesin θ 1 (i, j) eu]]>
Wherein, N indicates the number that electronics and the non-perforated wall of microchannel plate collide, and value is natural number.
5. the halation emulation mode in a kind of low-light level television imaging according to claim 1, wherein root described in step (5) According to the Two dimensional Distribution of image on fluorescent screen vegetarian refreshments, each pixel position (x, y) of the corresponding fluorescent screen of offset distance s (i, j) is calculated, It is calculated by following formula:
| (xdw) 2+ (ydh) 2-s (i, j) |<δ]]>
Wherein, dwIndicate the width of fluorescent screen single pixel, dhIndicate the height of fluorescent screen single pixel, δ=min (dw, dh)。
6. the halation emulation mode in a kind of low-light level television imaging according to claim 1, wherein meter described in step (6) The gray value Gray (x, y) of all pixels point on fluorescent screen is calculated, is calculated as follows:
(6a) calculates the quantization voltage V (x, y) of each pixel according to voltage signal conversion formula
V (x, y)=π 4 τ e (Defe) 2R&CenterDot;As&CenterDot;Gv&CenterDot;v&CenterDot;G& CenterDot;k&CenterDot;Ne (x, y) &CenterDot;e/a2]]>
Wherein τeIndicate the transmitance of coupled lens, DeIndicate coupled lens effective aperture, feIndicate that coupled lens focal length, R indicate Detector response rate, AsIndicate effectively photosensitive elemental area, GvIndicate that vision signal amplification factor, v indicate the work electricity of electron lens Pressure, G indicate the gain of microchannel plate, and k indicates fluorescent screen luminous efficacy, and a indicates electron-optical amplification factor;
(6b) calculates the gray value Gray (x, y) of each pixel according to the quantization voltage V (x, y) of each pixel
Gray (x, y)=255VMax-VMinV (x, y)]] >
Wherein VMaxAnd VMinRespectively indicate the maximum value and minimum value of the quantization voltage of all pixels point.
CN201710673245.4A 2017-07-28 2017-07-28 Halation emulation mode in low-light level television imaging Pending CN109308366A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710673245.4A CN109308366A (en) 2017-07-28 2017-07-28 Halation emulation mode in low-light level television imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710673245.4A CN109308366A (en) 2017-07-28 2017-07-28 Halation emulation mode in low-light level television imaging

Publications (1)

Publication Number Publication Date
CN109308366A true CN109308366A (en) 2019-02-05

Family

ID=65207595

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710673245.4A Pending CN109308366A (en) 2017-07-28 2017-07-28 Halation emulation mode in low-light level television imaging

Country Status (1)

Country Link
CN (1) CN109308366A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113432838A (en) * 2021-06-09 2021-09-24 北方夜视技术股份有限公司 Automatic testing system and method for signal-to-noise ratio and halo of low-light-level image intensifier

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113432838A (en) * 2021-06-09 2021-09-24 北方夜视技术股份有限公司 Automatic testing system and method for signal-to-noise ratio and halo of low-light-level image intensifier
CN113432838B (en) * 2021-06-09 2022-08-09 北方夜视技术股份有限公司 Automatic testing system and testing method for signal-to-noise ratio and halo of low-light-level image intensifier

Similar Documents

Publication Publication Date Title
CN104976960A (en) Raindrop physical property observation method and device
CN108279115B (en) Pulse discharge-based high ultra-thin flow field space-time evolution diagnosis device and method
US10247750B2 (en) Online measuring method of particle velocity in multiphase system
CN104280670A (en) Corona detection method based on solar-blind ultraviolet imager
CN104111485A (en) Stereo imaging based observation method for raindrop size distribution and other rainfall micro physical characteristics
CN103400005A (en) Quantifying method for intense light source to interfere imaging features of glimmer system
CN102252678B (en) High dynamic and high update rate star sensor and implementation method thereof
CN110763600A (en) Real-time online detection device for suspended particles
CN115639571A (en) Streak tube imaging laser radar image coordinate correction method and device
CN109308366A (en) Halation emulation mode in low-light level television imaging
CN107036710A (en) Using the light field luminous intensity distribution measurement method of multi-detector
CN106405566A (en) High-measurement-precision laser radar distance measurement method
CN103606185B (en) Halation emulation mode in low-light level television imaging
CN101886951B (en) System and method for analyzing movable target aiming performance of high repetition frequency laser
CN109709460A (en) A kind of blind type ultraviolet imagery system Photoncounting methods
JP5349257B2 (en) Visual status measurement device
Chen et al. Depth resolution improvement of streak tube imaging lidar system using three laser beams
Marshall et al. High-speed measurements of small-scale features in sprites: Sizes and lifetimes
CN104458192A (en) Measurement system and measurement method for wind tunnel flow field vortex motion frequency
Cui et al. Test and analysis of the halo in low-light-level image intensifiers
Ma et al. The research on the reconstruction of intensity image based on streak tube imaging lidar
CN102052905B (en) Device and method for testing effective diameter of cathode of micro-optic device with wide spectrum
Ye et al. Effect of time bin size on accuracy of streak tube imaging lidar
CN106383352A (en) Laser radar range finding method
CN112763189A (en) Measuring device for EBCMOS resolution parameter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190205