CN104267403A - Rapid dynamic focusing method for shallow-water multi-beam echo sounder - Google Patents
Rapid dynamic focusing method for shallow-water multi-beam echo sounder Download PDFInfo
- Publication number
- CN104267403A CN104267403A CN201410452835.0A CN201410452835A CN104267403A CN 104267403 A CN104267403 A CN 104267403A CN 201410452835 A CN201410452835 A CN 201410452835A CN 104267403 A CN104267403 A CN 104267403A
- Authority
- CN
- China
- Prior art keywords
- array element
- phase shift
- dynamic focusing
- plural form
- result
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
- G01C13/008—Surveying specially adapted to open water, e.g. sea, lake, river or canal measuring depth of open water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a rapid dynamic focusing method for a shallow-water multi-beam echo sounder. The method comprises the steps that underwater acoustic transducer array elements collect echo sound signals of the multi-beam echo sounder, and the sound signals are converted into electric signals and then processed to obtain digital signals; the digital signals are converted to obtain the plural form V<m> of the echo sound signals; a preset angle theta <k> is selected, and a phase shift parameter tau<mki> of the ith array element is rapidly figured out according to the sampling point number N and the array element number m; the plural form of the phase shift parameters of the array elements and the plural form V<m> of the echo sound signals are correspondingly multiply-accumulated according to the channel number and the sampling point number, and a single-beam dynamic focusing result is obtained; whether a preset angle which is not selected exists or not is judged, if yes, the third step and the fourth step are executed again, and if not, all single-beam dynamic focusing results are output. The rapid dynamic focusing method has the rapid calculating ability with no reduction of the precision of dynamic focusing carried out during multi-beam sounding, and is high in precision, simple in structure and easy to implement.
Description
Technical field
The invention belongs to multibeam echosounding field, particularly relate to and can solve near-field effect, a kind of shallow water multibeam echosounder quick dynamic focusing method.
Background technology
China has large-area marine site and inland river, rivers and lakes, and major part belongs to shallow-water environment, and for understanding water-bed landform fully, shallow water multibeam echosounder plays irreplaceable effect.Along with the continuous progress of science and technology, people have more and more higher requirement to the precision obtaining water-bed landform, and calculate for simplifying in conventional multibeam echosounder, the normal far-field approximation model that adopts, and are meeting r > > D
2can ensure high precision during/λ, wherein r is the distance of target from array element center, and D is the pore size of basic matrix, and λ is the wavelength of Received signal strength, and for not meeting the near field situation of this condition, total accuracy of sounding sharply declines.In whole investigative range, all ensure total accuracy of sounding in order to ensure multibeam echosounder, near field situation does not allow ignorance, and particularly in occasions such as the river courses that shallow water is main, this problem is particularly outstanding.
For solving Near-field Problems In Civil Engineering, domestic and international experts and scholars have carried out many theoretical researches, but are limited to theoretical simulation or segmentation focus approach more.Wherein theoretical simulation is only limitted to total accuracy of sounding analysis, cannot be applied in the Quick Measurement of shallow water multibeam echosounder; Although segmentation focus approach can realize but fast to lose measuring accuracy for cost.The people such as Chen Ying propose in its article delivered the quick acoustic imaging algorithm of focused beamforming " adopt classification ", Chen Ying, Ye Qinghua, the Huanghai Sea is peaceful. adopt the quick acoustic imaging algorithm [J] of classification focused beamforming, applied acoustics, 2008,27 (3): 207-210.Although computing velocity improves, dynamic beam focusing error is larger.Therefore, a kind of new quick high accuracy dynamic beam focus method is badly in need of.
Summary of the invention
The object of this invention is to provide have high-precision, a kind of quick dynamic focusing method of shallow water multibeam echosounder.
The present invention is achieved by the following technical solutions:
A kind of quick dynamic focusing method of shallow water multibeam echosounder, comprises following step:
Step one: underwater acoustic transducer array element gathers the echo acoustical signal of multibeam echosounder, acoustical signal is converted to electric signal, then obtains digital signal through process;
Step 2: plural form V digital signal being obtained echo acoustical signal through orthogonal transformation
m;
Step 3: select a predetermined angle θ
k, according to sampling period N and array element m, calculate the phase shift parameters τ of i-th array element
mki,
Wherein angle item factor of influence
distance terms factor of influence
c is the velocity of sound, f
sfor sampling rate, π is circular constant, and d is array element distance, and λ is wavelength, and k is numbers of beams;
By the phase shift parameters τ of i-th array element
mkicalculated by Coordinate Rotation Digital, obtain the phase shift parameters τ of i-th array element
mkiplural form;
Step 4: by the plural form V of the plural form of the phase shift parameters of each array element and echo acoustical signal
m, by channel number and the corresponding multiply accumulating of sampling period, obtain the dynamic focusing result of single beam
wherein M is array element sum;
Step 5: judge whether to there is non-selected predetermined angle, if existed, repeats step 3 ~ step 4, otherwise, export the dynamic focusing result of all single beams.
Beneficial effect of the present invention:
Provide one dynamic focus technology fast, according to the difference of goal-selling from the difference of basic matrix distance and array element number, calculate the phase shift parameters of corresponding each array element, then plural form phase shift parameters being converted to plural form and echoed signal, by the corresponding multiply accumulating of channel number, obtains quick dynamic focusing result.Introducing area and speed exchange principle simultaneously, solve the quick computational problem of dynamic focusing.The effect that the present invention is useful is: when not losing dynamic focusing precision in multibeam echosounding, has the ability calculated fast, and structure is simple, is easy to realize, can be widely used in the near field Quick Measurement precision improving shallow water multibeam echosounder.
Accompanying drawing explanation
Fig. 1 method flow diagram of the present invention,
The connected system block diagram of Fig. 2 specific embodiment of the invention.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further details.
1, the quick dynamic focusing method of shallow water multibeam echosounder, comprises step as follows:
Step one: gather the echoed signal of multibeam echosounder with underwater acoustic transducer, carry out acoustic-electric conversion, then through amplification, filtering, time-varying gain control with analog to digital conversion is down-sampled obtains digital signal;
Step 2: the plural form V result of step one being obtained echoed signal through orthogonal transformation
m;
Step 3: according to the predetermined angle θ of dynamic focusing Wave beam forming
k, sampling period N and array element m, calculate the phase shift parameters of each array element fast
wherein angle item factor of influence
distance terms factor of influence
c is the velocity of sound, f
sfor sampling rate, π be circular constant, d is array element distance, λ is wavelength, k is numbers of beams;
By phase shift parameters result τ
mkcalculated by quick Coordinate Rotation Digital, produce the sinusoidal result sin (τ of phase shift parameters respectively
mk) and cosine result cos (τ
mk);
Step 4: by the predetermined angle θ of dynamic focusing Wave beam forming
kphase shift parameters be converted to the plural form V of plural form and the multibeam echosounding echoed signal after orthogonal transformation
m, press corresponding θ by channel number and sampling period
km multiplication result of angle adds up, and obtains the dynamic focusing result of single wave beam
namely be the focus beam output of k θ angle of Wave beam forming, wherein M is array element sum, V
mrepresent the plural form of the echoed signal of array element m, j represents the imaginary part of plural number;
Step 5: by different θ
kangle repeats step 4, can obtain the result of the quick dynamic focusing of arranging by isogonism or equidistant pattern.
2, aforesaid step 3, comprises step as follows:
The first step: utilize and produce array element m from adding device circulation, its value is from 0 to (M-1), and wherein M is element number of array, and resets when launching multi-beam echo sounding signal;
Second step: previous step result is input to squarer and produces m
2, and by result and parameter
the result correspondence that storer exports is multiplied, wherein parameter
c is the velocity of sound, f
sfor sampling rate, π be circular constant, d is array element distance, λ is wavelength, θ
kfor a kth focus wave beam angle;
3rd step: by first step result and parameter γ
kthe result correspondence that storer exports is multiplied, wherein
π is circular constant, d is array element distance, λ is wavelength, θ
kfor a kth focus wave beam angle;
4th step: to launch multi-beam echo sounding signal for starting point, take sampling period as interval, produces sampling period count value N with from adding device, and produces (2 with divider
p)/N, wherein 16≤p≤64, and p is integer;
5th step: the result of step 2 is multiplied by step 4 result and divided by (2
p) or p position that fixed point results is moved to right, producing near-field effect affects item;
6) result of step 3 and step 5 is added up by totalizer, obtain phase shift parameters result
wherein τ
mkfor a kth focus wave beam angle is in the phase shift of array element m.
The quick dynamic focusing method of a kind of shallow water multibeam echosounder of the present invention, comprise fast dynamically phase shift computing module 1 and quick phase shift to add up summation module 2 two parts, be connected to the quick phase shift summation module 2 that adds up by electric signal after wherein dynamically phase shift computing module 1 dynamically calculates corresponding phase shift fast according to the change of sampling period and array element number fast and calculate the result exporting dynamic focusing.At this with implementation example explanation in programmable gate array.
Dynamic phase shift computing module 1, comprises array element m counter 4, array element m squarer 3, parameter fast
with parameter γ
kstorer 5,10, sampling period N counter 7, divider 8, two multipliers 6,9, a totalizer 11 and phase shift result memory 12; For calculating dynamic phase shift fast, need internal logic Resource Design array element m counter 4 be used, realize array element m squarer 3 by multiplier resources, by internal storage Resource Design parameter
with parameter γ
kstorer 5,10, with internal logic Resource Design sampling period N counter 7, with internal multiplier Resource Design two multipliers 6,9, with internal logic Resource Design divider 8, use internal logic Resource Design totalizer 11 again, finally with memory resource design phase shift result memory 12, wherein parameter
with parameter γ
kmeet
wherein c is the velocity of sound, f
sfor sampling rate, π be circular constant, d is array element distance, λ is wavelength, θ
kfor a kth focus wave beam angle.Wherein array element m counter 4 circulates and produces channel number, by array element m squarer 3 and parameter
the pre-stored data of storer is multiplied, and then divided by sampling period N counter, can obtain
the output of array element m counter 4 simultaneously, direct and parameter γ
kthe pre-stored data of storer is multiplied and can obtains m γ
kpart, two parts are added can obtain dynamic phase shift
concrete structure is as shown in Fig. 21:
Wherein m is array element number, N is sampling period,
for custom parameter, the γ of a kth focus wave beam angle
kfor the custom parameter of a kth focus wave beam angle.Often receive one group of hyperchannel raw data, sampling period N counter 7, from increasing 1, then delivers to divider 8 as electronic signals; Array element m counter 4 Zi increasing to array number M, is then transferred to array element m squarer 3 and multiplier 9 from 1 in electrical signal form;
Array element m squarer 3 completes square operation after receiving array element m value, then delivers to multiplier 6;
Multiplier 9 is address lookup parameter γ with array element number after receiving the change of array element number
kstorer 10, and be multiplied number corresponding with Query Result for array element, result is outputted to totalizer 11 in electrical signal form; Multiplier 6 is address lookup parameter with array element number after receiving the change of array element number square
storer 5, and be multiplied square corresponding with Query Result for array element number, result is outputted to divider 8 in electrical signal form;
Divider 8 uses the result of multiplier 6 divided by sampling period N counter 7 result, and outputs to totalizer 11 in electrical signal form; The result correspondence of divider 8 and multiplier 9 is added by totalizer 11, can obtain fast dynamically phase shift result of calculation.
Quick phase shift adds up summation module 2, comprises Coordinate Rotation Digital computing module 13, raw storage 14, complex multiplier 15, totalizer 16 and dynamic beam focusing results storer 17.
Result V (θ is formed for calculating dynamic near field focused beam fast
k), wherein
In formula: V (θ
k) represent k θ angle of Wave beam forming focus beam export, V
mrepresent the raw data input signal of array element m, M represents the overall channel number of raw data, and j represents the imaginary part of plural number.Internal logic Resource Design Coordinate Rotation Digital computing module 13 need be used, with internal storage Resource Design raw storage 14, with internal multiplier and totalizer Resource Design complex multiplier 15, with internal additions device Resource Design totalizer 16 with internal storage Resource Design dynamic beam focusing results storer 17.Wherein Coordinate Rotation Digital computing module 13 is considering fast and the problem of precision, and the Coordinate Rotation Digital of 18 iteration can be adopted to calculate.By dynamic phase shift τ
mkresult is sent in Coordinate Rotation Digital computing module 13, corresponding sine and cosine is produced through interative computation, then complex multiplier 15 is passed through with the raw data of raw storage 14, finally obtaining dynamic near field focused beam by totalizer 16 forms result, and concrete structure is as shown in Fig. 22.
Wherein Coordinate Rotation Digital computing module 13 receives the phase shift result of calculation from quick dynamically phase shift computing module 1, and calculates output sine and cosine result; Coordinate Rotation Digital computing module 13 and raw storage 14 export electric signal respectively and are connected to the complex multiplication that complex multiplier 15 completes raw data and phase shift parameters, reach the object of phase shift; Complex multiplier 15, totalizer 16 and dynamic beam focusing results storer 17, press electric signal successively and connect, by by the raw data of each passage after the phase shift of correspondence is cumulative, dynamic focusing result is stored in storer.Coordinate Rotation Digital computing module 13 receives the phase-shift value from quick dynamically phase shift computing module 1 in electrical signal form, calculated by Coordinate Rotation Digital, obtain sine and the cosine result of corresponding phase shift, and be combined into the plural form of phase shift, and be sent to complex multiplier 15 in electrical signal form; Corresponding for data in the result of calculation of Coordinate Rotation Digital computing module 13 and raw storage 14 complex multiplication is obtained complex result by complex multiplier 15, transfers to totalizer 16 in electrical signal form; Totalizer 16 obtains dynamic focusing result by cumulative for complex multiplier 15 result of corresponding array element number 1 to M, and is stored in dynamic beam focusing results storer 17.
The invention provides one dynamic focus technology fast, according to the difference of goal-selling from the difference of basic matrix distance and array element number, calculate the phase shift parameters of corresponding each array element, then plural form phase shift parameters being converted to plural form and echoed signal, by the corresponding multiply accumulating of channel number, obtains quick dynamic focusing result.Introducing area and speed exchange principle simultaneously, approximate due to what do not adopt classification focusing and segmentation to focus on, relative classification focused beamforming and segmentation focused beamforming have the high feature of focus beam precision, solve the problem of traditional shallow water multibeam echosounder near field high precision depth measurement, reach the object calculating near field dynamic beam focused beamforming fast.
Claims (1)
1. the quick dynamic focusing method of shallow water multibeam echosounder, is characterized in that, comprise following step:
Step one: underwater acoustic transducer array element gathers the echo acoustical signal of multibeam echosounder, acoustical signal is converted to electric signal, then obtains digital signal through process;
Step 2: plural form V digital signal being obtained echo acoustical signal through orthogonal transformation
m;
Step 3: select a predetermined angle θ
k, according to sampling period N and array element m, calculate the phase shift parameters τ of i-th array element
mki,
Wherein angle item factor of influence
distance terms factor of influence
c is the velocity of sound, f
sfor sampling rate, π is circular constant, and d is array element distance, and λ is wavelength, and k is numbers of beams;
By the phase shift parameters τ of i-th array element
mkicalculated by Coordinate Rotation Digital, obtain the phase shift parameters τ of i-th array element
mkiplural form;
Step 4: by the plural form V of the plural form of the phase shift parameters of each array element and echo acoustical signal
m, by channel number and the corresponding multiply accumulating of sampling period, obtain the dynamic focusing result of single beam
wherein M is array element sum;
Step 5: judge whether to there is non-selected predetermined angle, if existed, repeats step 3 ~ step 4, otherwise, export the dynamic focusing result of all single beams.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410452835.0A CN104267403A (en) | 2014-09-05 | 2014-09-05 | Rapid dynamic focusing method for shallow-water multi-beam echo sounder |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410452835.0A CN104267403A (en) | 2014-09-05 | 2014-09-05 | Rapid dynamic focusing method for shallow-water multi-beam echo sounder |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104267403A true CN104267403A (en) | 2015-01-07 |
Family
ID=52158943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410452835.0A Pending CN104267403A (en) | 2014-09-05 | 2014-09-05 | Rapid dynamic focusing method for shallow-water multi-beam echo sounder |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104267403A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891835A (en) * | 2016-06-16 | 2016-08-24 | 北京海卓同创科技有限公司 | Real-time dynamic focusing wave beam forming method and system |
CN109283511A (en) * | 2018-09-01 | 2019-01-29 | 哈尔滨工程大学 | A kind of wide covering multi-beam reception basic matrix calibration method |
CN113534161A (en) * | 2021-05-31 | 2021-10-22 | 中国船舶重工集团公司第七一五研究所 | Beam mirror image focusing method for remotely positioning underwater sound source |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080137482A1 (en) * | 2006-12-07 | 2008-06-12 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Time multiplexed method and device for beamforming |
CN101644765A (en) * | 2009-04-23 | 2010-02-10 | 中国科学院声学研究所 | Amplitude and phase error correction method used for linear array of underwater acoustic transducer |
CN101936756A (en) * | 2010-08-31 | 2011-01-05 | 华南理工大学 | Multifrequency phased array ultrasonic Doppler flow detection system and method |
-
2014
- 2014-09-05 CN CN201410452835.0A patent/CN104267403A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080137482A1 (en) * | 2006-12-07 | 2008-06-12 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Time multiplexed method and device for beamforming |
CN101644765A (en) * | 2009-04-23 | 2010-02-10 | 中国科学院声学研究所 | Amplitude and phase error correction method used for linear array of underwater acoustic transducer |
CN101936756A (en) * | 2010-08-31 | 2011-01-05 | 华南理工大学 | Multifrequency phased array ultrasonic Doppler flow detection system and method |
Non-Patent Citations (2)
Title |
---|
周泽民 等: "《基于信号相位匹配原理的聚焦波束形成算法》", 《声学技术》 * |
李海森 等: "《基于FPGA的多波束实时动态聚焦波束形成方法》", 《振动与冲击》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891835A (en) * | 2016-06-16 | 2016-08-24 | 北京海卓同创科技有限公司 | Real-time dynamic focusing wave beam forming method and system |
CN109283511A (en) * | 2018-09-01 | 2019-01-29 | 哈尔滨工程大学 | A kind of wide covering multi-beam reception basic matrix calibration method |
CN109283511B (en) * | 2018-09-01 | 2022-07-29 | 哈尔滨工程大学 | Wide-coverage multi-beam receiving array calibration method |
CN113534161A (en) * | 2021-05-31 | 2021-10-22 | 中国船舶重工集团公司第七一五研究所 | Beam mirror image focusing method for remotely positioning underwater sound source |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105891835A (en) | Real-time dynamic focusing wave beam forming method and system | |
CN103197317B (en) | Synthetic aperture radar (SAR) imaging method based on field programmable gate array (FPGA) | |
USRE45823E1 (en) | System and method of acoustic doppler beamforming | |
CN106772326A (en) | A kind of multiple submatrixes synthetic aperture sonar phase error analysis method | |
CN101644773A (en) | Real-time frequency domain super-resolution direction estimation method and device | |
CN104502904B (en) | Torpedo homing beam sharpening method | |
CN102508251B (en) | Method for rapidly implementing sector conversion in multi-beam image sonar | |
CN110836981A (en) | Layered water flow high-resolution radial acoustic Doppler frequency measurement method | |
CN110133580B (en) | Stochastic resonance enhanced acoustic vector signal orientation method | |
CN204203461U (en) | A kind of 3-D scanning acoustics imaging device | |
CN104267403A (en) | Rapid dynamic focusing method for shallow-water multi-beam echo sounder | |
CN108519602A (en) | Cold seepage fast imaging system and its array signal processing method | |
CN103728464B (en) | A kind of assembled pulse speed-measuring method for acoustic Doppler fluid velocity profile instrument | |
CN102841343A (en) | Echo sounding apparatus calibration system based on industrial computer and calibration method | |
CN102928619B (en) | Signal processing method of broadband beam-control acoustic Doppler flow testing system | |
CN101504458A (en) | Phase filtering based beam forming method | |
CN103267965B (en) | Multi-target micro-variation measurement data processing system and method | |
CN109375197B (en) | Small-size vector array low-frequency scattering correction method | |
CN102087357A (en) | Method for echo direction estimation of sensor array and for multibeam echo depth sounding and bottom detection | |
CN110824193A (en) | Non-uniform water velocity estimation method based on multi-beam radial flow velocity measurement | |
Dajun et al. | Design of high accuracy ultra short baseline underwater acoustic position system | |
CN107255804A (en) | Moveing target indication filter filtering method based on clutter sketch map | |
Huang et al. | Research on analyzing and processing methods of ocean sonar signals | |
CN101598788A (en) | The rapid simulation method of synthetic aperture sonar signal | |
CN104111349A (en) | Doppler velometer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150107 |
|
WD01 | Invention patent application deemed withdrawn after publication |