CN103369466B - A kind of map match assists indoor orientation method - Google Patents

A kind of map match assists indoor orientation method Download PDF

Info

Publication number
CN103369466B
CN103369466B CN201310289054.XA CN201310289054A CN103369466B CN 103369466 B CN103369466 B CN 103369466B CN 201310289054 A CN201310289054 A CN 201310289054A CN 103369466 B CN103369466 B CN 103369466B
Authority
CN
China
Prior art keywords
pixel
point
match
indoor
line model
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
CN201310289054.XA
Other languages
Chinese (zh)
Other versions
CN103369466A (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.)
Harbin University of Technology Robot Group Co., Ltd.
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201310289054.XA priority Critical patent/CN103369466B/en
Publication of CN103369466A publication Critical patent/CN103369466A/en
Application granted granted Critical
Publication of CN103369466B publication Critical patent/CN103369466B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

A kind of map match assists indoor orientation method, relates to map match assisted location method, and solving existing indoor map matching algorithm needs to come the extra measurement data of sensor and the process of iterative algorithm search, causes the problem of computing redundancy.Arrange access point, reference point at indoor objects localizing environment, set up location fingerprint database; The CAD plane graph of floor to be positioned is utilized to generate the image array I of the floor level image of load(ing) point line model match; Compute location result for the elements of a fix of test point; The elements of a fix are converted to pixel, judge whether have building body to stop the mid point of the n-th pixel and (n-1) individual pixel between the n-th pixel and (n-1) individual pixel then the pixel that this point of distance is nearest on Point and Line Model is searched for and by nearest pixel as revised location pixel; The location pixel of acquisition is converted to actual position coordinates, completes indoor positioning.The present invention can be widely used in indoor positioning.

Description

A kind of map match assists indoor orientation method
Technical field
The present invention relates to map match assisted location method.
Background technology
Map match assistant positioning system has become an important component part of land vehicle navigation system.High-precision numerical map can not only provide the positional information of vehicle for driver, can also be used for improving the positioning performance of integrated navigation system.At present, this technology has been widely used in outdoor vehicle navigator fix.In indoor positioning field, user adds in the randomness that diversified indoor environment is moved the difficulty that map matching technology applies in indoor environment.Meanwhile, in the complicated radio environment in indoor, due to the impact of various error, the positioning result that the fingerprint matching indoor orientation method based on WLAN obtains may depart from the normal zone of action of people.Therefore, adopt map-matching algorithm to utilize cartographic information, by the positioning result using software engineering effectively to revise WLAN, improve positioning precision.
Current, scientific research personnel has proposed a variety of indoor map matching algorithm.But, existing indoor map matching algorithm usually and some transducer conbined usage, such as, inertial sensor, gyroscope and electronic compass etc.In addition, indoor GoogleMap is utilized to adopt iterative algorithm to search for the indoor map matching algorithm of most matching result in addition.
Summary of the invention
The present invention needs to come the extra measurement data of sensor and the process of iterative algorithm search to solve existing indoor map matching algorithm, causes the problem of computing redundancy, thus provides a kind of map match to assist indoor orientation method.
A kind of map match assists indoor orientation method, and it comprises the steps:
Step one: arrange N number of access point AP at indoor objects localizing environment i, i=1,2 ..., N, creates two-dimentional cartesian coordinate system under described localizing environment, chooses M reference point RP in described localizing environment j, j=1,2 ..., M, the coordinate recording each reference point and W the received signal strength RSS sample gathered at each reference point place, according to M reference point RP jcoordinate and the signal strength signal intensity RSS Sample Establishing location fingerprint database of each reference point;
Step 2: utilize the CAD plane graph of floor to be positioned to generate the image array I of the floor level image of load(ing) point line model match;
Described step 2: utilize the CAD plane graph of floor to be positioned to generate the image array I of the floor level image of load(ing) point line model matchprocess be:
Point and Line Model is loaded into after in the CAD plane graph of floor, the CAD plane graph of a load(ing) point line model can be obtained; CAD plane graph after load(ing) point line model and original CAD plane graph are saved as picture format respectively; Two width pictures are converted to respectively two bianry image matrix I modeland I no_model; In these two bianry image matrixes, the pixel in people's daily routines region is set to numerical value+1, and the pixel of building body structure and Point and Line Model is set to numerical value 0;
Obtain indoor map image array I match:
I no_model-2×(I model⊕I no_model)=I match
Wherein, ⊕ is xor operation.
Step 3: test point is set in described localizing environment, after receiving H signal strength signal intensity RSS sample at a test point place, ask for signal strength signal intensity RSS sample average and utilize fingerprint matching positioning mode k nearest neighbor algorithm KNN and weighting k nearest neighbor algorithm WKNN compute location result for the elements of a fix of test point;
Described step 3: after receiving N number of signal strength signal intensity RSS sample at a test point place, asks for signal strength signal intensity RSS sample average and utilizes fingerprint matching positioning mode k nearest neighbor algorithm KNN and weighting k nearest neighbor algorithm WKNN compute location result for the process of the elements of a fix is:
Formula is obtained according to fingerprint matching positioning mode k nearest neighbor algorithm KNN:
D i = ( Σ j = 1 N | rss j - rss i , j RP | q ) 1 / q , i = 1 , · · · , M L r = ( x r RP , y r RP ) , D r ∈ Min _ K ( D i ) , r = 1 , · · · , K ( x ‾ , y ‾ ) = 1 K Σ r = 1 K ( x r RP , y r RP ) = 1 K Σ r = 1 K L r
Wherein, rss jaP isignal strength signal intensity RSS data mean value, represent at reference point RP in fingerprint database jplace measure from AP isignal strength signal intensity RSS data mean value; D iit is signal strength signal intensity RSS distance; The type of q representation signal intensity RSS distance; Min_K (D i) be the set of minimum K signal strength signal intensity RSS distance; D r, r=1 ..., K is K minimum signal strength signal intensity RSS distance; belong to set Min_K (D i) middle distance D rthe coordinate of corresponding reference point;
When being assigned with weights by the position coordinates of K RP selected, described in be assigned with weights be that the inverse of signal strength signal intensity RSS distance then carries out WKNN algorithm:
( x ‾ , y ‾ ) = Σ r = 1 K ( ( x r RP , y r RP ) / D r ) Σ r = 1 K ( 1 / D r ) = Σ r = 1 K ( L r / D r ) Σ r = 1 K ( 1 / D r ) .
Step 4: the n-th positioning result obtained according to step 3, is converted to image array I matchin pixel, n>=2;
Step 5: whether whether pixel described in determining step four, in building body structure, be namely I match(x n, y n)=0; If be 0, this pixel, in building body structure, forwards step 7 to, otherwise forwards step 6 to;
Step 6: judge whether have building body to stop between the n-th pixel and (n-1) individual pixel, namely whether two pixels are stopped by 0 value pixel; If there is building body to stop, then forward step 7 to, otherwise forward step 8 to;
Step 7: the mid point calculating the n-th pixel and (n-1) individual pixel then the pixel that this point of distance is nearest on Point and Line Model is searched for and by nearest pixel as revised location pixel;
Step 8: the location pixel of acquisition is converted to actual position coordinates, completes indoor positioning.
The present invention adopts a kind of map match to assist indoor orientation method to realize indoor hi-Fix.The present invention makes full use of indoor map information, the positioning result of calculating is transformed to the relevant position of the image array for map match, by identifying the building body structure in this image array, irrational positioning result being adapted on Point and Line Model, more accurate positioning result can be obtained.Simultaneously relative to traditional indoor map matching algorithm, the inventive method can independently use does not need sensor measurement data, does not need iterative process.
Accompanying drawing explanation
Fig. 1 is the flow chart that a kind of map match of the present invention assists indoor orientation method;
Fig. 2 is the floor schematic diagram of specific embodiment;
Fig. 3 works as K=9 for described in specific embodiment, during H=2, is respectively k nearest neighbor algorithm KNN and positioning result schematic diagram of the present invention;
Fig. 4 works as K=9 for described in specific embodiment, during H=2, is respectively weighting k nearest neighbor algorithm WKNN and positioning result schematic diagram of the present invention.
Embodiment
Embodiment one, composition graphs 1 illustrate this embodiment.A kind of map match assists indoor orientation method, and it comprises the steps:
Step one: arrange N number of access point AP at indoor objects localizing environment i, i=1,2 ..., N, creates two-dimentional cartesian coordinate system under described localizing environment, chooses M reference point RP in described localizing environment j, j=1,2 ..., M, the coordinate recording each reference point and W the received signal strength RSS sample gathered at each reference point place, according to M reference point RP jcoordinate and the signal strength signal intensity RSS Sample Establishing location fingerprint database of each reference point;
Step 2: utilize the CAD plane graph of floor to be positioned to generate the image array I of the floor level image of load(ing) point line model match;
Step 3: test point is set in described localizing environment, after receiving H signal strength signal intensity RSS sample at a test point place, ask for signal strength signal intensity RSS sample average and utilize fingerprint matching positioning mode k nearest neighbor algorithm KNN and weighting k nearest neighbor algorithm WKNN compute location result for the elements of a fix of test point;
Step 4: the n-th positioning result obtained according to step 3, is converted to image array I matchin pixel, n>=2;
Step 5: whether whether pixel described in determining step four, in building body structure, be namely I match(x n, y n)=0; If be 0, this pixel, in building body structure, forwards step 7 to, otherwise forwards step 6 to;
Step 6: judge whether have building body to stop between the n-th pixel and (n-1) individual pixel, namely whether two pixels are stopped by 0 value pixel; If there is building body to stop, then forward step 7 to, otherwise forward step 8 to;
Step 7: the mid point calculating the n-th pixel and (n-1) individual pixel then the pixel that this point of distance is nearest on Point and Line Model is searched for and by nearest pixel as revised location pixel;
Step 8: the location pixel of acquisition is converted to actual position coordinates, completes indoor positioning.
Embodiment two, this embodiment and embodiment one is unlike described step 2: utilize the CAD plane graph of floor to be positioned to generate the image array I of the floor level image of load(ing) point line model matchprocess be:
Point and Line Model is loaded into after in the CAD plane graph of floor, the CAD plane graph of a load(ing) point line model can be obtained; CAD plane graph after load(ing) point line model and original CAD plane graph are saved as picture format respectively; Two width pictures are converted to respectively two bianry image matrix I modeland I no_model; In these two bianry image matrixes, the pixel in people's daily routines region is set to numerical value+1, the pixel of building body structure and Point and Line Model is set to numerical value 0;
Obtain indoor map image array I match:
I no_model-2×(I model⊕I no_model)=I match
Wherein, ⊕ is xor operation.So, at the image array I that floor plan is new matchin, the pixel of Point and Line Model is replaced by numerical value-1, and the pixel of people zone of action is replaced by numerical value+1, and the pixel of building body structure is replaced by numerical value 0.
Embodiment three, this embodiment and embodiment one or two are unlike described step 3: in described localizing environment, arrange test point, after receiving H signal strength signal intensity RSS sample at a test point place, ask for signal strength signal intensity RSS sample average and utilize fingerprint matching positioning mode k nearest neighbor algorithm KNN and weighting k nearest neighbor algorithm WKNN compute location result for the process of the elements of a fix of test point is:
Formula is obtained according to fingerprint matching positioning mode k nearest neighbor algorithm KNN:
D i = ( Σ j = 1 N | rss j - rss i , j RP | q ) 1 / q , i = 1 , · · · , M L r = ( x r RP , y r RP ) , D r ∈ Min _ K ( D i ) , r = 1 , · · · , K ( x ‾ , y ‾ ) = 1 K Σ r = 1 K ( x r RP , y r RP ) = 1 K Σ r = 1 K L r
Wherein, rss jaP isignal strength signal intensity RSS data mean value, represent at reference point RP in fingerprint database jplace measure from AP isignal strength signal intensity RSS data mean value; D iit is signal strength signal intensity RSS distance; The type of q representation signal intensity RSS distance; Min_K (D i) be the set of minimum K signal strength signal intensity RSS distance; D r, r=1 ..., K is K minimum signal strength signal intensity RSS distance; r=1 ..., K belongs to set Min_K (D i) middle distance D rthe coordinate of corresponding reference point;
When being assigned with weights by the position coordinates of K RP selected, described in be assigned with weights be that the inverse of signal strength signal intensity RSS distance then carries out WKNN algorithm:
( x ‾ , y ‾ ) = Σ r = 1 K ( ( x r RP , y r RP ) / D r ) Σ r = 1 K ( 1 / D r ) = Σ r = 1 K ( L r / D r ) Σ r = 1 K ( 1 / D r ) .
Embodiment four, this embodiment and embodiment three, unlike described q=2, namely chooses Euclidean distance.
Embodiment five, this embodiment and embodiment one or four are positive integer unlike described N, M, and W, H are positive integer.
Specific embodiment:
Composition graphs 2-4 illustrates this specific embodiment.The validity of the method is tested under the laboratory experiment environment shown in Fig. 2.Wherein, 9 LinksysWAP54GAP are arranged in the laboratory experiment environment of 24.9m × 66.4m.Experiment path is from the A point room to the B point in the corridor of 3 meters wide.Experiment utilizes an Asus A8F notebook computer image data.It is equipped with IntelPRO/Wireless3945ABG wireless network card and RSS data acquisition software NetStumbler, sampling rate 2 RSS data samples per second.In off-line phase, select in a room to select 67 reference points in the corridor of 24 reference points and 3 meters wide, with collections in 150 seconds totally 300 RSS samples in each reference point.At on-line stage, choose 65 test points along experiment path A altogether to B, spacing is 0.6m.Consider quantity and the layout of reference point, arrange K=9, when H=2, q=2, calculate traditional KNN and WKNN algorithm and assist the positioning performance after indoor orientation method with application map match, the mean error of experimental result is more as shown in table 1.
Table 1 algorithm performance compares
Algorithm Mean error (m)
KNN 2.73
WKNN 2.70
KNN+Map 2.45
WKNN+Map 2.47
Shown in Fig. 3 and Fig. 4, at turning, corridor, because Structure in Complex Structure is larger on radio propagation impact, so the positioning result of KNN algorithm and WKNN algorithm is poor, have a lot of positioning result to have passed through body of wall, this is impossible in practice.Apply after map match proposed by the invention assists indoor orientation method, positioning performance obviously improves, and has been adapted on the Point and Line Model proposed by irrational positioning result, obtains the elements of a fix more accurately.The mean error of positioning result also drops to 2.45 meters and 2.47 meters respectively from 2.73 of KNN and WKNN algorithm meters and 2.70 meters.Therefore, above-mentioned experiment has convincingly demonstrated value and the validity that map match proposed by the invention assists indoor orientation method.

Claims (2)

1. map match assists an indoor orientation method, it is characterized in that it comprises the steps:
Step one: arrange N number of access point AP at indoor objects localizing environment i, i=1,2 ..., N, creates two-dimentional cartesian coordinate system under described localizing environment, chooses M reference point RP in described localizing environment j, j=1,2 ..., M, the coordinate recording each reference point and W the received signal strength RSS sample gathered at each reference point place, according to M reference point RP jcoordinate and the signal strength signal intensity RSS Sample Establishing location fingerprint database of each reference point;
Step 2: utilize the CAD plane graph of floor to be positioned to generate the image array I of the floor level image of load(ing) point line model match;
Step 3: test point is set in described localizing environment, after receiving H signal strength signal intensity RSS sample at a test point place, ask for signal strength signal intensity RSS sample average and utilize fingerprint matching positioning mode k nearest neighbor algorithm KNN or weighting k nearest neighbor algorithm WKNN compute location result for the elements of a fix of test point;
Step 4: the n-th positioning result obtained according to step 3, is converted to image array I matchin pixel, n>=2;
Step 5: whether whether pixel described in determining step four, in building body structure, be namely I match(x n, y n)=0; If be 0, this pixel, in building body structure, forwards step 7 to, otherwise forwards step 6 to;
Step 6: judge whether have building body to stop between the n-th pixel and (n-1) individual pixel, namely whether two pixels are stopped by 0 value pixel; If there is building body to stop, then forward step 7 to, otherwise forward step 8 to;
Step 7: the mid point calculating the n-th pixel and (n-1) individual pixel then the pixel P that this point of distance is nearest on Point and Line Model is searched for n'=(x n', y n'), and by nearest pixel P n' as revised location pixel;
Step 8: the revised location pixel obtained is converted to actual position coordinates, completes indoor positioning.
2. a kind of map match according to claim 1 assists indoor orientation method, it is characterized in that described step 2: utilize the CAD plane graph of floor to be positioned to generate the image array I of the floor level image of load(ing) point line model matchprocess be:
Point and Line Model is loaded into after in the CAD plane graph of floor, the CAD plane graph of a load(ing) point line model can be obtained; CAD plane graph after load(ing) point line model and original CAD plane graph are saved as picture format respectively; Two width pictures are converted to respectively two bianry image matrix I modeland I no_model; In these two bianry image matrixes, the pixel in people's daily routines region is set to numerical value+1, the pixel of building body structure and Point and Line Model is set to numerical value 0; Obtain indoor map image array I match:
I no_model-2×(I model⊕I no_model)=I match
Wherein, ⊕ is xor operation.
CN201310289054.XA 2013-07-10 2013-07-10 A kind of map match assists indoor orientation method Active CN103369466B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310289054.XA CN103369466B (en) 2013-07-10 2013-07-10 A kind of map match assists indoor orientation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310289054.XA CN103369466B (en) 2013-07-10 2013-07-10 A kind of map match assists indoor orientation method

Publications (2)

Publication Number Publication Date
CN103369466A CN103369466A (en) 2013-10-23
CN103369466B true CN103369466B (en) 2016-03-02

Family

ID=49369806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310289054.XA Active CN103369466B (en) 2013-07-10 2013-07-10 A kind of map match assists indoor orientation method

Country Status (1)

Country Link
CN (1) CN103369466B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103822626B (en) * 2014-02-17 2018-04-10 惠州Tcl移动通信有限公司 Mobile terminal and its generation numerical map or air navigation aid, device
CN103926611B (en) * 2014-05-07 2016-08-17 中科院成都信息技术股份有限公司 A kind of indoor positioning data optimization methods in real time
CN104202818B (en) * 2014-09-03 2015-10-07 创业软件股份有限公司 A kind of floor recognition methods distance weighted based on building open edge
CN104573651B (en) * 2014-12-31 2018-02-13 北京天诚盛业科技有限公司 Fingerprint identification method and device
CN105547299A (en) * 2015-12-29 2016-05-04 哈尔滨工业大学 WLAN (Wireless Local Area Network) indoor positioning method with automatic generation of radio map corner matrix based on track matching
CN105676173B (en) * 2016-01-14 2018-05-15 广州市万联信息科技有限公司 Indoor locating system and indoor orientation method
CN106802411A (en) * 2016-12-26 2017-06-06 纵横壹旅游科技(成都)有限公司 A kind of local area localization method and system
CN107529145A (en) * 2017-09-15 2017-12-29 南京轩世琪源软件科技有限公司 The localization method of handheld terminal in a kind of high-precision office building
CN110345935B (en) * 2019-06-04 2020-12-29 中国地质大学(武汉) Indoor map matching and positioning method
CN111009036B (en) * 2019-12-10 2023-11-21 北京歌尔泰克科技有限公司 Grid map correction method and device in synchronous positioning and map construction
US11698467B2 (en) * 2021-08-30 2023-07-11 Nanning Fulian Fugui Precision Industrial Co., Ltd. Indoor positioning method based on image visual features and electronic device
TWI799969B (en) * 2021-08-30 2023-04-21 新加坡商鴻運科股份有限公司 Indoor positioning method based on image visual features and electronic device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101639527A (en) * 2009-09-03 2010-02-03 哈尔滨工业大学 K nearest fuzzy clustering WLAN indoor locating method based on REE-P
CN102802260A (en) * 2012-08-15 2012-11-28 哈尔滨工业大学 WLAN indoor positioning method based on matrix correlation
CN103167606A (en) * 2013-03-12 2013-06-19 钱钢 Wireless local area network (WLAN) indoor positioning method based on sparse representation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101639527A (en) * 2009-09-03 2010-02-03 哈尔滨工业大学 K nearest fuzzy clustering WLAN indoor locating method based on REE-P
CN102802260A (en) * 2012-08-15 2012-11-28 哈尔滨工业大学 WLAN indoor positioning method based on matrix correlation
CN103167606A (en) * 2013-03-12 2013-06-19 钱钢 Wireless local area network (WLAN) indoor positioning method based on sparse representation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A KNN-based two-step fuzzy clustering weighted algorithm for WLAN indoor positioning;Xu Yubin, Sun Yongliang, Ma Lin;《HIGH TECHNOLOGY LETTERS》;20110930;全文 *

Also Published As

Publication number Publication date
CN103369466A (en) 2013-10-23

Similar Documents

Publication Publication Date Title
CN103369466B (en) A kind of map match assists indoor orientation method
CN102802260B (en) WLAN indoor positioning method based on matrix correlation
CN107318084A (en) A kind of fingerprint positioning method and device based on optimal similarity
US20150262481A1 (en) System and method to determine occurrence of platoon
CN106597363A (en) Pedestrian location method in indoor WLAN environment
CN102012522B (en) Method and device for evaluating seismic acquisition data
CN109029429B (en) WiFi and geomagnetic fingerprint based multi-classifier global dynamic fusion positioning method
CN103220777A (en) Mobile device positioning system
CN101105524A (en) Quick cadastral survey measurement and data processing method based on real-time dynamic positioning technology
CN107607122B (en) The building of location fingerprint library and dynamic updating method towards indoor positioning
CN102928858A (en) GNSS (Global Navigation Satellite System) single-point dynamic positioning method based on improved expanded Kalman filtering
CN105046046B (en) A kind of Ensemble Kalman Filter localization method
CN108592903B (en) Vehicle geomagnetic matching positioning method based on road network geomagnetic reference library
CN102656581B (en) Messaging device, information processing method and storage medium
CN106291609A (en) A kind of RTK precision assessment method
CN104792321A (en) Auxiliary-positioning-based land information acquisition system and method
Zhuang et al. Autonomous WLAN heading and position for smartphones
US6473690B1 (en) Three-dimensional space curve comparison using spatial angle variance metric and applications thereof
CN109764876A (en) The multi-modal fusion localization method of unmanned platform
CN105717517A (en) Vehicle-mounted Beidou multi-mode GNSS high-precision road basic data collection method
CN113465616B (en) Track abnormal point detection method and device, electronic equipment and storage medium
CN103616027A (en) Gravity matching method based on improved mean-square deviation (MSD)
CN114710742A (en) Indoor positioning method for constructing fingerprint map based on multi-chain interpolation
CN114488247A (en) Method for analyzing mobility of equipment based on high-precision Beidou differential positioning
CN113484843A (en) Method and device for determining external parameters between laser radar and integrated navigation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20190611

Address after: 150000 Heilongjiang Harbin Dalian economic and Trade Zone, the North Road and Xingkai Road intersection

Patentee after: Harbin University of Technology Robot Group Co., Ltd.

Address before: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin

Patentee before: Harbin Institute of Technology