CN103471586B - The terminal combinations localization method that a kind of sensor is auxiliary and device - Google Patents

The terminal combinations localization method that a kind of sensor is auxiliary and device Download PDF

Info

Publication number
CN103471586B
CN103471586B CN201310431932.7A CN201310431932A CN103471586B CN 103471586 B CN103471586 B CN 103471586B CN 201310431932 A CN201310431932 A CN 201310431932A CN 103471586 B CN103471586 B CN 103471586B
Authority
CN
China
Prior art keywords
location
personnel
time
positioning
localization method
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.)
Active
Application number
CN201310431932.7A
Other languages
Chinese (zh)
Other versions
CN103471586A (en
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.)
Chongqing Institute of Green and Intelligent Technology of CAS
Original Assignee
Chongqing Institute of Green and Intelligent Technology of CAS
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 Chongqing Institute of Green and Intelligent Technology of CAS filed Critical Chongqing Institute of Green and Intelligent Technology of CAS
Priority to CN201310431932.7A priority Critical patent/CN103471586B/en
Publication of CN103471586A publication Critical patent/CN103471586A/en
Application granted granted Critical
Publication of CN103471586B publication Critical patent/CN103471586B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

The invention discloses the auxiliary terminal combinations localization method of a kind of sensor and device, belong to space orientation technique field.This method obtains the elements of a fix by WLAN, WSN, GPS wireless location and satellite positioning tech, utilize accelerometer to treat the walking step number of locating personnel between this is located and positioning time last time is poor to measure, and according to estimating the translational speed of step size computation personnel to be positioned; By filtering algorithm, filtering is carried out to this positioning result, obtain the accurate coordinates of this location and the translational speed in all directions, translational speed is combined with above-mentioned location technology and realizes locating next time.Positioning precision is promoted to 3 ~ 5m by 5 ~ 10m when not increasing hardware facility by the method for the invention, and then significantly can promote indoor and outdoor seamless location service quality, expanding location service new business, apply in fields such as public safety, emergency management and rescue, production safety, Internet of Things, smart city, obtain important economic results in society.

Description

The terminal combinations localization method that a kind of sensor is auxiliary and device
Technical field
The invention belongs to space orientation technique field, relate to the auxiliary terminal combinations localization method of a kind of sensor and device.
Background technology
According to statistics, have more than 80% relevant with " position " in the information used people.Along with the development of mobile interchange technology, location-based service (LocationBasedSevers, LBS) technology and market will obtain development at full speed.Increasing along with the flourish of LBS and heavy construction, the demand of people to the service of indoor and outdoor seamless location constantly increases, and the fields such as public safety, emergency management and rescue, anti-terrorism are anti-probably, production safety, Internet of Things, special population love, smart city also all need to use positional information accurately.
GPS, the satellite navigation systems such as the Big Dipper achieve outdoor positioning, but cannot in-door covering be realized, in recent years for the further investigation of general fit calculation and Distributed Communication Technology, indoor wireless communication and network technology is made to enter the develop rapidly stage, and then derived based on WLAN(WirelessLocalAreaNetworks, WLAN (wireless local area network)), Bluetooth(bluetooth), WSN(wirelesssensornetwork, wireless sensor network) indoor positioning technologies, that location-based service enters the new stage of indoor and outdoor seamless service by outdoor service, wherein, due to WLAN popularizing on mobile phone and the application of other intelligent electronic device, location technology based on WLAN is the most universal.
At present, indoor positioning technologies emphasis place based on networks such as WLAN, WSN can reach 2 ~ 3m, but input cost with expend time in all larger, precision in most of commercial situation is generally 5 ~ 10m, in outdoor positioning, the civilian precision of mobile phone of the location technologies such as GPS is also generally 5 ~ 10m in urban environment.How to utilize the hardware facility of existing mobile phone, intelligent terminal, the positioning precision promoting the systems such as WLAN, WSN, GPS is further the target that location-based service manufacturer is constantly pursued.Current smart mobile phone, intelligent terminal etc. all possess the sensor such as accelerometer, electronic compass, all obvious acceleration change can be produced because pedestrian often transports to move a step, meter step can be carried out to the motion of personnel by accelerometer, by the electronic compass then direction run to personnel of cocoa, both combine, and can try to achieve the movable information of personnel, thus adopt the filtering techniques such as conventional Kalman, the velocity information that the positional information obtain systems such as WLAN, WSN, GPS and sensor obtain carries out filtering jointly, realize integrated positioning, improve positioning precision.
But at present the electronic compass of a lot of mobile phone points to and is inaccurate, and error reaches more than 40 °, and the direction of motion of people is pointed to electronic compass and inconsistent (as adept machine, retrogressing is walked) sometimes.Above-mentioned situation will cause pedestrian movement's walking direction to occur comparatively serious mistake, now, adopt the sensor auxiliary positioning of classic method not only can not promote positioning precision, positioning precision can be caused even on the contrary to decline.
Summary of the invention
In view of this, the terminal combinations localization method that the object of the present invention is to provide a kind of sensor auxiliary and device, the method utilizes walking step number and the direction situation of change of accelerometer in mobile phone or other intelligent terminal and electronic compass survey crew, the imperfect velocity information that cannot be obtained by sensor in exactly determined situation in personnel's direction of motion (only has velocity magnitude, lack velocity reversal) combine with wireless/satellite positioning tech, improve positioning precision.
For achieving the above object, the invention provides following technical scheme:
The terminal combinations localization method that a kind of sensor is auxiliary, the elements of a fix are obtained by WLAN, WSN, GPS wireless location and satellite positioning tech, utilize accelerometer to treat the walking step number of locating personnel between this is located and positioning time last time is poor to measure, and according to estimating the translational speed of step size computation personnel to be positioned; By filtering algorithm, filtering is carried out to this positioning result, obtain the accurate coordinates of this location and the translational speed in all directions, translational speed is combined with above-mentioned location technology and realizes locating next time.
Further, described filtering algorithm adopts Kalman filtering, and personnel to be positioned obtain the coordinate (x of kth time location by location technology k', y k'), the walking step number N during utilizing accelerometer measures kth time location to locate for-1 time with kth k, obtain the velocity magnitude of current persons | v k|=N kl k/ T, wherein l kfor the personnel's step-length estimated, T is the interval time of twice location; By (x k', y k') carry out Kalman filtering as positioning result, obtain result (x after filtering k, y k, vx k, vy k), wherein x k, y k, vx k, vy kbe respectively the x coordinate finally obtaining kth time location, y coordinate, x direction movement velocity size, y direction movement velocity size; On carrying out once Kalman filtering iterative process before, to vx k, vy kcarry out correction as follows:
vx k = vx k · | v k | / vx k 2 + vy k 2 vy k = vy k · | v k | / vx k 2 + vy k 2 , Thus will | v k| information and wireless, satellite positioning information combine, improve positioning precision.
Further, when personnel's direction of motion keeps stablizing, current positioning result (x is utilized k, y k) with the N-1 time front positioning result (x k-N+1, y k-N+1) line direction calculating personnel motion velocity reversal.
Further, described personnel's direction of motion keeps stable condition to be: think that direction of motion keeps stable when the standard deviation of electronic compass measured value is less than pre-determined threshold delt in N continuous time location.
Further, according to the velocity reversal obtained and velocity magnitude, speed amount is decomposed x, y-axis, obtain measured value vx k', vy k', obtain the secondary coordinate (x located of kth with wireless, satellite positioning tech k', y k') common as input quantity, carry out Kalman filtering calculating.
Further, when direction of motion keeps stablizing, current positioning result (x is utilized k, y k), positioning result (x before the N-1 time k-N+1, y k-N+1), accumulative walking number Δ N between twice location k, k+1-N, step-length l is estimated in adjustment kfor l k = ( x k - x k - N + 1 ) 2 + ( y k - y k - N + 1 ) 2 Δ N k , k + 1 - N .
Further, when the standard deviation of electronic compass measured value is greater than pre-determined threshold delt2 in continuous N time location, think that personnel's direction of motion changes, now need to recalculate velocity reversal and realize location.
Present invention also offers the terminal combinations locating device that a kind of sensor is auxiliary, this application of installation combinations thereof localization method, adopt walking step number and the direction situation of change of accelerometer and the electronic compass survey crew carried in current mobile phone or other intelligent terminal, thus be combined with conventional alignment systems, realize accurately locating.
Beneficial effect of the present invention is: utilize the accelerometer in existing mobile phone or intelligent terminal and electronic compass, / satellite positioning tech wireless with existing WLAN, WSN, GPS etc. is combined, when not increasing hardware facility, general positioning precision is promoted to 3 ~ 5m by 5 ~ 10m, and then significantly can promote indoor and outdoor seamless location service quality, expanding location service new business, the fields such as, production safety, Internet of Things, special population love anti-probably in public safety, emergency management and rescue, anti-terrorism, smart city are applied, and obtain important economic results in society.
Accompanying drawing explanation
In order to make object of the present invention, technical scheme and beneficial effect clearly, the invention provides following accompanying drawing and being described:
Fig. 1 is the schematic flow sheet of the method for the invention.
Embodiment
Below in conjunction with accompanying drawing, the preferred embodiments of the present invention are described in detail.
Fig. 1 is the schematic flow sheet of the method for the invention, the invention provides a kind of terminal combinations localization method possessing higher accuracy, specific as follows:
Utilize WLAN, WSN, GPS etc. wireless/satellite positioning tech obtains the coordinate (x of kth time location k', y k'), utilize accelerometer measures personnel kth time location and kth locate for-1 time between walking step number N k, obtain the velocity magnitude of current persons | v k|=N kl k/ T, wherein l kfor the personnel's step-length estimated, T is the interval time of twice location.Make (x k', y k') carry out Kalman filtering as positioning result, obtain result (x after filtering k, y k, vx k, vy k), wherein x k, y k, vx k, vy kbe respectively the x coordinate finally obtaining kth time location, y coordinate, x direction movement velocity size, y direction movement velocity size.On carrying out once Kalman filtering iterative process before, to vx k, vy kcarry out correction as follows:
vx k = vx k · | v k | / vx k 2 + vy k 2 vy k = vy k · | v k | / vx k 2 + vy k 2 - - - ( 1 )
Thus will | v k| information and wireless/satellite positioning information combine, improve positioning precision.
In the present embodiment, Kalman filtering process is specific as follows:
Output vector is made to be m k=[x k, y k, vx k, vy k] t, initial value m 1can obtain according to front twice wireless/satnav result, state-transition matrix is A = 1 0 T 0 0 1 0 T 0 0 1 0 0 0 0 1 , Process noise covariance matrix Q = q 1 0 0 0 0 q 2 0 0 0 0 q 3 0 0 0 0 q 4 , Wherein q 1, q 2, q 3, q 4for preset value, representative value can get 0.2,0.2,0.05,0.05 respectively, calculation matrix H = 1 0 0 0 0 1 0 0 , Measurement noises covariance matrix R = r 1 0 0 r 2 , R 1, r 2for preset value, determine according to accuracy of wireless location system, representative value can get 49,49, measurement result Z respectively k=[x k', y k'] t.Matrix P is the matrix of a 4*4, can upgrade in Kalman's iterative process always, and initial value can rule of thumb not grudged, as P = 9 0 2 0 0 9 0 2 2 0 1 0 0 2 0 1 . Following kth time output vector is m kcomputation process as follows:
S k=A·m k-1(2)
P k=A·P k-1·A T+A T·Q·A(3)
K k=P kH T/(HPH T+R)(4)
P k=P k-K(HPH T+R)K T(5)
m k=S k+K k(Z-HS k-1)(6)
m k ( 3,1 ) = m k ( 3,1 ) · | v k | / vx k 2 + vy k 2 m k ( 4,1 ) = m k ( 4,1 ) · | v k | / vx k 2 + vy k 2 - - - ( 7 )
When personnel's direction of motion keeps stablizing, utilize current positioning result (x k, y k) with the N-1 time front positioning result (x k-N+1, y k-N+1) line direction calculating personnel motion velocity reversal (such as can be set to also can be according to different coordinates definition and adopt different formulas).Above-mentioned judgement personnel direction of motion keeps stable condition to be: when the standard deviation of electronic compass measured value in N continuous time location (can be set to 15 in the present embodiment) is less than pre-determined threshold delt (can be set to 10 ° in the present embodiment).That is: although do not obtain accurate personnel's direction of motion information by electronic compass, when personnel direction can be utilized constant compass can not illustrate change this characteristic judge whether personnel's direction of motion is stablized.
The velocity reversal obtained according to said process again and velocity magnitude, decompose x by speed amount, y-axis, obtain measured value vx k', vy k', obtain the secondary coordinate (x located of kth with wireless/satellite positioning tech k', y k') jointly carry out Kalman filtering calculating.Specific implementation process for: calculation matrix described in above-mentioned steps is rewritten as H = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 , Measurement noises covariance matrix is rewritten as R = r 1 0 0 0 0 r 2 0 0 0 0 r 3 0 0 0 0 r 4 , Measurement result is rewritten as Z k=[x k', y k', vx k', vy k'] t, calculate kth time positioning result by formula (2) to formula (6), positioning result can be promoted further.
Meanwhile, when described personnel's direction of motion keeps stablizing, current positioning result (x is utilized k, y k), positioning result (x before the N-1 time k-N+1, y k-N+1), accumulative walking number Δ N between twice location k, k+1-N, in set-up procedure 1, estimate step-length l kfor l k = ( x k - x k - N + 1 ) 2 + ( y k - y k - N + 1 ) 2 Δ N k , k + 1 - N , Further lifting positioning precision.
In above process, when the standard deviation of electronic compass measured value in continuous N time location (can be set to 5 in the present embodiment) is greater than pre-determined threshold delt2 (30 ° can be set in the present embodiment), think that personnel move to change, then velocity reversal and velocity magnitude are redefined.
What finally illustrate is, above preferred embodiment is only in order to illustrate technical scheme of the present invention and unrestricted, although by above preferred embodiment to invention has been detailed description, but those skilled in the art are to be understood that, various change can be made to it in the form and details, and not depart from claims of the present invention limited range.

Claims (7)

1. the terminal combinations localization method that a sensor is auxiliary, the elements of a fix are obtained by WLAN, WSN wireless location and satellite positioning tech, it is characterized in that: utilize accelerometer to treat the walking step number of locating personnel between this is located and positioning time last time is poor and measure, and according to estimating the translational speed of step size computation personnel to be positioned; By filtering algorithm, filtering is carried out to this positioning result, obtain the accurate coordinates of this location and the translational speed in all directions, translational speed is combined with above-mentioned location technology and realizes locating next time;
Described filtering algorithm adopts Kalman filtering, personnel to be positioned by location technology obtain kth time location coordinate (x ' k, y ' k), the walking step number N during utilizing accelerometer measures kth time location to locate for-1 time with kth k, obtain the velocity magnitude of current persons | v k|=N kl k/ T, wherein l kfor the personnel's step-length estimated, T is the interval time of twice location; Will (x ' k, y ' k) carry out Kalman filtering as positioning result, obtain result (x after filtering k, y k, vx k, vy k), wherein x k, y k, vx k, vy kbe respectively the x coordinate finally obtaining kth time location, y coordinate, x direction movement velocity size, y direction movement velocity size; On carrying out once Kalman filtering iterative process before, to vx k, vy kcarry out correction as follows:
v x k = v x k · | v k | / vx k 2 + vy k 2 vy k = vy k · | v k | / vx k 2 + vy k 2 , Thus will | v k| information and wireless, satellite positioning information combine, improve positioning precision.
2. the terminal combinations localization method that a kind of sensor according to claim 1 is auxiliary, is characterized in that: when personnel's direction of motion keeps stablizing, utilize current positioning result (x k, y k) with the N-1 time front positioning result (x k-N+1, y k-N+1) line direction calculating personnel motion velocity reversal.
3. the terminal combinations localization method that a kind of sensor according to claim 2 is auxiliary, is characterized in that: described personnel's direction of motion keeps stable condition to be: think that direction of motion keeps stable when the standard deviation of electronic compass measured value is less than pre-determined threshold delt in N continuous time location.
4. the terminal combinations localization method that a kind of sensor according to claim 3 is auxiliary, is characterized in that: according to the velocity reversal obtained and velocity magnitude, speed amount is decomposed x, y-axis, obtain measured value vx ' k, vy ' k, with wireless, satellite positioning tech obtain coordinate that kth time locates (x ' k, y ' k) common as input quantity, carry out Kalman filtering calculating.
5. the terminal combinations localization method that a kind of sensor according to claim 4 is auxiliary, is characterized in that: when direction of motion keeps stablizing, utilize current positioning result (x k, y k), positioning result (x before the N-1 time k-N+1, y k-N+1), accumulative walking number Δ N between twice location k, k+1-N, step-length l is estimated in adjustment kfor
6. the terminal combinations localization method that a kind of sensor according to claim 5 is auxiliary, it is characterized in that: when in continuous N time location, the standard deviation of electronic compass measured value is greater than pre-determined threshold delt2, think that personnel's direction of motion changes, now need to recalculate velocity reversal and component realization location thereof.
7. the device of the terminal combinations localization method using sensor according to any one of claim 1 to 6 to assist.
CN201310431932.7A 2013-09-22 2013-09-22 The terminal combinations localization method that a kind of sensor is auxiliary and device Active CN103471586B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310431932.7A CN103471586B (en) 2013-09-22 2013-09-22 The terminal combinations localization method that a kind of sensor is auxiliary and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310431932.7A CN103471586B (en) 2013-09-22 2013-09-22 The terminal combinations localization method that a kind of sensor is auxiliary and device

Publications (2)

Publication Number Publication Date
CN103471586A CN103471586A (en) 2013-12-25
CN103471586B true CN103471586B (en) 2016-01-06

Family

ID=49796538

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310431932.7A Active CN103471586B (en) 2013-09-22 2013-09-22 The terminal combinations localization method that a kind of sensor is auxiliary and device

Country Status (1)

Country Link
CN (1) CN103471586B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104023398B (en) * 2014-06-23 2017-03-22 哈尔滨工业大学 Method for quickly establishing positioning Radio Map based on WiFi and accelerometer
CN105548956B (en) * 2014-10-28 2017-11-21 ***通信集团公司 Pedestrian's alignment system, method and relevant device in car networking system
CN105788349B (en) * 2014-12-22 2019-01-15 联想(北京)有限公司 A kind of information processing method and electronic equipment
US9781569B2 (en) * 2015-03-12 2017-10-03 GM Global Technology Operations LLC Systems and methods for resolving positional ambiguities using access point information
CN113766636A (en) * 2015-05-12 2021-12-07 三星电子株式会社 Apparatus and method for estimating position in wireless communication system
CN105975753A (en) * 2016-04-28 2016-09-28 努比亚技术有限公司 Method for calculating moving speed of user and mobile terminal
EP3454614B1 (en) * 2016-05-20 2021-05-05 Huawei Technologies Co., Ltd. Positioning method and mobile device
CN107576330A (en) * 2017-09-07 2018-01-12 西北大学 A kind of localization method of the indoor dynamic sensing strategy based on WLAN fingerprints
CN111586580A (en) * 2020-04-29 2020-08-25 杭州十域科技有限公司 Position event capturing method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841334A (en) * 2012-08-24 2012-12-26 北京邮电大学 Method and device for acquiring locating point
CN103033835A (en) * 2011-09-30 2013-04-10 卡西欧计算机株式会社 Positioning apparatus and positioning method
CN103298105A (en) * 2012-02-24 2013-09-11 美国博通公司 Wireless communication device capable of accurately performing position estimations

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8831909B2 (en) * 2011-09-22 2014-09-09 Microsoft Corporation Step detection and step length estimation
JP5949061B2 (en) * 2012-03-30 2016-07-06 ソニー株式会社 Information processing apparatus, information processing method, and program

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103033835A (en) * 2011-09-30 2013-04-10 卡西欧计算机株式会社 Positioning apparatus and positioning method
CN103298105A (en) * 2012-02-24 2013-09-11 美国博通公司 Wireless communication device capable of accurately performing position estimations
CN102841334A (en) * 2012-08-24 2012-12-26 北京邮电大学 Method and device for acquiring locating point

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于Android的室内导航定位***设计与实现;胡天琨;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20130715;第29-46页 *

Also Published As

Publication number Publication date
CN103471586A (en) 2013-12-25

Similar Documents

Publication Publication Date Title
CN103471586B (en) The terminal combinations localization method that a kind of sensor is auxiliary and device
Ban et al. Indoor positioning method integrating pedestrian Dead Reckoning with magnetic field and WiFi fingerprints
Zhou et al. Activity sequence-based indoor pedestrian localization using smartphones
CN109470238B (en) Positioning method and device and mobile terminal
CN105589064A (en) Rapid establishing and dynamic updating system and method for WLAN position fingerprint database
CN104655137A (en) Wi-Fi signal fingerprint positioning algorithm for assisting in speculating flight tracks of pedestrians
CN102322858B (en) Geomagnetic matching navigation method for geomagnetic-strapdown inertial navigation integrated navigation system
CN104197935B (en) Indoor localization method based on mobile intelligent terminal
CN109974694B (en) Indoor pedestrian 3D positioning method based on UWB/IMU/barometer
CN105021198B (en) A kind of location estimation method navigated based on MULTISENSOR INTEGRATION
CN102944889A (en) Positioning system and positioning method for sensor auxiliary positioning terminal
CN108519084B (en) Pedestrian geomagnetic positioning method and system assisted by dead reckoning
CN104121905A (en) Course angle obtaining method based on inertial sensor
CN105704652A (en) Method for building and optimizing fingerprint database in WLAN/Bluetooth positioning processes
CN102281629A (en) Positioning method and positioning device
Xuan et al. Making indoor maps with portable accelerometer and magnetometer
CN113923596A (en) Indoor positioning method, device, equipment and medium
CN105512722A (en) BP neural network-based Kalman gain correction method
CN105357754A (en) Wireless-network-based mobile node combined positioning method
CN104640204A (en) Wireless sensor network node positioning method in indirect wave environment
CN109708631A (en) A kind of pedestrian navigation navigational calibration method based on mobile phone modal idenlification
CN114111802A (en) Pedestrian dead reckoning assisted UWB positioning method
CN106211318A (en) A kind of path loss localization method based on WiFi and system
CN103596265A (en) Multiple-user indoor positioning method based on voice distance measuring and movement vector
CN110530359A (en) It is a kind of to be corrected based on bluetooth and map match seeks an indoor orientation method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: CHONGQING INSTITUTE OF GREEN AND INTELLIGENT TECHN

Free format text: FORMER OWNER: CHONGQING INSTITUTE OF GREEN AND INTELLIGENT TECHNOLOGY

Effective date: 20140912

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 401122 YUBEI, CHONGQING TO: 400714 BEIPEI, CHONGQING

TA01 Transfer of patent application right

Effective date of registration: 20140912

Address after: 400714 Chongqing Road, Beibei District, No. 266

Applicant after: Chongqing Institute of Green and Intelligent Technology of the Chinese Academy of Sciences

Address before: 401122, B, block 9, Han Chinese Center, 85 Jin Yu Road, North New District, Chongqing

Applicant before: Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences

C14 Grant of patent or utility model
GR01 Patent grant