CN110811660B - 一种校正ct射线束硬化伪影的方法 - Google Patents
一种校正ct射线束硬化伪影的方法 Download PDFInfo
- Publication number
- CN110811660B CN110811660B CN201911022549.XA CN201911022549A CN110811660B CN 110811660 B CN110811660 B CN 110811660B CN 201911022549 A CN201911022549 A CN 201911022549A CN 110811660 B CN110811660 B CN 110811660B
- Authority
- CN
- China
- Prior art keywords
- projection
- energy
- image
- ray
- soft tissue
- 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
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种校正CT射线束硬化伪影的方法,该方法基于X射线能谱分析对原图像进行补偿,消除X射线束的硬化伪影;即,根据已知的扫描条件,获得X射线的能谱分布,计算等效能量;然后,进行扫描,根据建像条件重建原始图像,将原始图像分割成软组织、骨组织和空气部分,分别对不同部分进行能谱投影和等效能量投影;计算等效能量投影与能谱投影的差值;对投影差值进行反投影得到补偿图像,将补偿图像与原始图像相加得到射线束硬化伪影校正后的目标图像。
Description
技术领域
本发明涉及一种校正CT射线束硬化伪影的方法,本发明属于图像处理技术领域。
背景技术
CT是1969年由英国电子工程师亨斯菲尔德.G.N.发明,1972年问世的。CT不同于普通X线成像,它是用X射线束对人体层面进行扫描,取得信息,经计算机处理而获得的重建图像,是数字成像而不是模拟成像。CT所显示的断层解剖图像,其密度分辨力明显优于X线成像,使X线成像不能显示的解剖结构及其病变得以显影,从而显著扩大了人体的检查范围,提高了病变检出率和诊断的准确率。
CT成像基本原理是用X射线束对人体检查部位一定厚度的层面进行扫描,由探测器接收透过该层面的X射线束,转变为可见光后,由光电转换器转变为电信号,再经模拟/数字转换器转为数字信号,输入计算机进行图像成像处理。
在X射线束穿透人体器官或组织时,由于人体器官或组织是由多种不同成分、不同密度的物质构成的,所以,实际上各点对X射线束的吸收系数是不同的,同一种物质对不同能量的光子的衰减系数也不同,而X射线束的能谱是多色的,在X射线穿过人体器官或组织的过程中,低能量的光子通常容易被吸收衰减的快,高能量的光子吸收的慢衰减的慢,因此在射线穿过物质过程中射线的能谱发生变化,射线随着穿过的路径变长,越来越不易吸收,出现了X射线束硬化现象。
由于X射线的硬化现象,导致重建的CT图像中出现伪影,这种伪影称为射线束硬化伪影。通常,射线束硬化伪影分两类,一类是:在扫描均匀的组织时,产生的呈“杯状”或“盆状”的伪影;另一类是:扫描多种不同密度物体时,在致密组织之间形成的暗色的条带。这些伪影都极易引起医生误诊,影响医生诊断。
通常,为了消除这两类伪影需要进行两步处理,首先,通过采集水模数据模拟人体的软组织,通过将水模图像校正成均匀图像从而确定多项式的系数,利用多项式拟合的方法校正射线束穿过均匀软组织产生的射线束硬化伪影;第二步,将图像中的骨头分割出来,对骨头部分的图像通过再次多项式拟合的方式消除骨头引起的射线束硬化伪影。为消除射线束硬化引起的伪影,需要进行两步校正操作。另外,通过多项式拟合的方式确定多项式系数来源于经验,缺少理论依据,因此,这种校正伪影的方法的鲁棒性比较差。
发明内容
为了更好地校正CT射线束硬化伪影,本发明的目的是提供一种基于X射线能谱分析补偿的方式校正CT射线束硬化伪影的方法。
为实现上述目的,本发明采用以下技术方案:一种校正CT射线束硬化伪影的方法,具体步骤如下:
S1:获取X光源能谱分布,计算等效能量;
S2:根据预设的扫描部位、扫描强度对人体进行X射线束扫描,图像重建,得到待校正的带有射线束硬化伪影的原始图像I0;
S3:将原始图像I0分割出软组织部分、骨组织部分和空气部分,得到分割图像I1;
S4:计算分割图像I1中X射线穿过软组织和骨组织的路径长度;
S5:根据S4计算出的X射线穿过软组织和骨组织的路径长度,进行等效能量投影得到等效能量投影P1;
等效能量投影P1计算公式如下:
P1=μb(e)Lb+μs(e)Ls
其中:P1为等效能量对应的投影;μb(e)为等效能量值对应的骨头衰减系数,该衰减系数通过公开资料可以查得;Lb为射线穿过骨头的路径长度;μs(e)为等效能量对应的软组织的衰减系数,该衰减系数通过公开资料可以查得;Ls为射线穿过软组织的路径长度;
S6:根据S4计算出的射线穿过软组织和骨组织的路径长度,进行X射线束投影得到能谱投影P2;
能谱投影P2计算公式如下:
其中:P2为射线能谱投影的结果;S(i)e为第i个能量值对应的谱分布;μf(i)为第i个能量值对应的滤过材料的衰减系数,该衰减系数通过公开资料可以查得;Lf为射线穿过滤过材料的路径长度;μb(i)为第i个能量值对应的骨头衰减系数,该衰减系数通过公开资料可以查得;Lb为射线穿过骨头的路径长度;μs(i)为第i个能量值对应的软组织的衰减系数,该衰减系数通过公开资料可以查得;Ls为射线穿过软组织的路径长度;
S7:计算等效能量投影P1与能谱投影P2的投影差值:P=P1-P2;
S8:将投影差值进行滤波反投影得到补偿图像I2;
S9:将原始图像I0与补偿图像I2相加得到校正后的消除射线束硬化伪影的目标图像I。
附图说明
图1是本发明校正CT射线束硬化伪影的方法流程图;
图2是120KV管电压对应的X光源能谱分布图;
图3是本发明具体实施例中的原始图像;
图4是本发明具体实施例中的分割图像;
图5是本发明具体实施例中的补偿图像;
图6是本发明具体实施例中校正后的消除X射线束硬化伪影后的图像。
具体实施方式
下面结合附图和实施例对本发明的结构及特征进行详细说明。需要说明的是,可以对此处公开的实施例做出各种修改,因此,说明书中公开的实施例不应该视为对本发明的限制,而仅是作为实施例的范例,其目的是使本发明的特征显而易见。
本发明校正CT射线束硬化伪影的方法是基于X射线能谱分析对原图像进行补偿,消除X射线束的硬化伪影;即,首先,根据已知的扫描条件,获得X射线的能谱分布,计算等效能量;然后,进行扫描,根据建像条件重建原始图像,将原始图像分割成软组织、骨组织和空气部分,分别对不同部分进行能谱投影和等效能量投影;计算等效能量投影与能谱投影的投影差值;对投影差值进行反投影得到补偿图像,将补偿图像与原始图像相加得到射线束硬化伪影校正后的目标图像。
如图1所示,本发明校正CT射线束硬化伪影的方法具体步骤如下:
S1:获取X光源能谱分布,计算等效能量
对于任何一台CT机,其出厂后即可确定其发射X射线的球管的管电压,根据球管管电压即可获取该CT机X光源能谱分布。如图2所示,为120KV管电压对应的X光源能谱分布图。
获取CT机X光源能谱后,以摆放在扫描视野中心的直径为20cm的水模作为模拟扫描对象,进行投影计算水模的能谱投影值,将该投影值作为目标投影;再以X光能谱中最高能量的40%-60%为能量,以0.1为步长计算单能量投影,计算X光谱中每个能量对应的水模投影结果;与目标投影均方差最接近的单能量投影对应的能量值即为等效能量。
S2:根据预设的扫描部位、扫描强度对人体进行X射线束扫描,图像重建,得到待校正的带有射线束硬化伪影的原始图像I0。
图3为对人体头部进行扫描,然后进行图像重建得到的待校正的带有射线束硬化伪影的原始图像。
S3:将原始图像I0分割出软组织部分、骨组织部分和空气部分,得到分割图像I1。
通过CT值阈值将原始图像I0分割出软组织部分、骨组织部分和空气部分。在本发明具体实施例中,如图4所示,将CT值在[-200,100]之间的图像定义为软组织;将CT值大于100的图像定义为骨组织;将CT值小于-200的图像定义为空气。
S4:计算分割图像I1中X射线穿过软组织和骨组织的路径长度。
根据CT机的几何结构,模拟计算从球管发射出的X光射线到探测器检测单元穿过被扫描患者软组织和骨组织的路径长度,为了更接近射线实际传射的路径,本发明推荐通过扇束投影的方式计算射线穿过各组织路径,扇束投影方式计算射线穿过组织路径的方法是本领域技术人员惯用技术手段,在此不再赘述。
S5:根据S4计算出的射线穿过软组织和骨组织的路径长度,进行等效能量投影得到等效能量投影P1。
等效能量投影P1计算公式如下:
P1=μb(e)Lb+μs(e)Ls
其中:P1为等效能量对应的投影;μb(e)为等效能量值对应的骨头衰减系数,该衰减系数通过公开资料可以查得;Lb为射线穿过骨头的路径长度;μs(e)为等效能量对应的软组织的衰减系数,该衰减系数通过公开资料可以查得;Ls为射线穿过软组织的路径长度。
S6:根据S4计算出的射线穿过软组织和骨组织的路径长度,进行X射线束投影得到能谱投影P2。
能谱投影P2计算公式如下:
其中:P2为射线能谱投影的结果;S(i)e为第i个能量值对应的谱分布;μf(i)为第i个能量值对应的滤过材料(例如铜或铝或聚四氟乙烯)的衰减系数,该衰减系数通过公开资料可以查得;Lf为射线穿过滤过材料(例如铜或铝或聚四氟乙烯)的路径长度;μb(i)为第i个能量值对应的骨头衰减系数,该衰减系数通过公开资料可以查得;Lb为射线穿过骨头的路径长度;μs(i)为第i个能量值对应的软组织的衰减系数,该衰减系数通过公开资料可以查得;Ls为射线穿过软组织的路径长度。
S7:计算等效能量投影P1与能谱投影P2的投影差值:P=P1-P2;
S8:将投影差值进行滤波反投影得到补偿图像I2(如图5所示),其建像条件与方法需要和原始图像I0重建时的条件方法保持一致;
S9:将原始图像I0与补偿图像I2相加得到校正后的消除射线束硬化伪影的目标图像I(如图6所示)。
通过对比本发明具体实施例的原始图像(如图3所示)和校正后的目标图像(如图6所示)可知,本发明具有以下优点:
本发明基于射线束谱分析方法校正硬化束伪影,是从射线束硬化产生的原因出发,计算扫描物体单能X光与多能X光之间成像的差异,将差异补偿到目标图像上,达到消除硬化束伪影。该方法根据硬化产生的机理通过严格计算得到的补偿项,因此,得到的补偿项更加精准,硬化束伪影校正的更加精确,消除硬化束伪影的鲁棒性更好。另外,该方法可以同时计算出软组织和骨组织利用单能X光与多能X光成像之间的差异,也就是说通过一步处理同时消除两类射线束硬化伪影,使得校正流程更加简便。
最后应说明的是:以上所述的各实施例仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或全部技术特征进行等同替换;而这些修改或替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (4)
1.一种校正CT射线束硬化伪影的方法,具体步骤如下:
S1:获取X光源能谱分布,计算等效能量;
S2:根据预设的扫描部位、扫描强度对人体进行X射线束扫描,图像重建,得到待校正的带有射线束硬化伪影的原始图像I0;
S3:将原始图像I0分割出软组织部分、骨组织部分和空气部分,得到分割图像I1;
S4:计算分割图像I1中X射线穿过软组织和骨组织的路径长度;
S5:根据S4计算出的X射线穿过软组织和骨组织的路径长度,进行等效能量投影得到等效能量投影P1;
等效能量投影P1计算公式如下:
P1=μb(e)Lb+μs(e)Ls
其中:P1为等效能量对应的投影;μb(e)为等效能量值对应的骨头衰减系数,该衰减系数通过公开资料可以查得;Lb为射线穿过骨头的路径长度;μs(e)为等效能量对应的软组织的衰减系数,该衰减系数通过公开资料可以查得;Ls为射线穿过软组织的路径长度;
S6:根据S4计算出的射线穿过软组织和骨组织的路径长度,进行X射线束投影得到能谱投影P2;
能谱投影P2计算公式如下:
其中:P2为射线能谱投影的结果;S(i)e为第i个能量值对应的谱分布;μf(i)为第i个能量值对应的滤过材料的衰减系数,该衰减系数通过公开资料可以查得;Lf为射线穿过滤过材料的路径长度;μb(i)为第i个能量值对应的骨头衰减系数,该衰减系数通过公开资料可以查得;Lb为射线穿过骨头的路径长度;μs(i)为第i个能量值对应的软组织的衰减系数,该衰减系数通过公开资料可以查得;Ls为射线穿过软组织的路径长度;
S7:计算等效能量投影P1与能谱投影P2的投影差值:P=P1-P2;
S8:将投影差值进行滤波反投影得到补偿图像I2;
S9:将原始图像I0与补偿图像I2相加得到校正后的消除射线束硬化伪影的目标图像I。
2.根据权利要求1所述的校正CT射线束硬化伪影的方法,其特征在于:所述步骤S1中的等效能量获取方法为:以摆放在扫描视野中心的直径为20cm的水模作为模拟扫描对象,进行投影,计算水模的能谱投影值,将该投影值作为目标投影;再以X光能谱中最高能量的40%-60%为能量,以0.1为步长计算单能量投影,计算X光谱中每个能量对应的水模投影结果;与目标投影均方差最接近的单能量投影对应的能量值即为等效能量。
3.根据权利要求2所述的校正CT射线束硬化伪影的方法,其特征在于:所述步骤S3根据CT值阈值将原始图像分割出软组织部分、骨组织部分和空气部分。
4.根据权利要求3所述的校正CT射线束硬化伪影的方法,其特征在于:所述步骤S6中的滤过材料为铜、铝、聚四氟乙烯之一。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911022549.XA CN110811660B (zh) | 2019-10-25 | 2019-10-25 | 一种校正ct射线束硬化伪影的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911022549.XA CN110811660B (zh) | 2019-10-25 | 2019-10-25 | 一种校正ct射线束硬化伪影的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110811660A CN110811660A (zh) | 2020-02-21 |
CN110811660B true CN110811660B (zh) | 2023-05-12 |
Family
ID=69550518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911022549.XA Active CN110811660B (zh) | 2019-10-25 | 2019-10-25 | 一种校正ct射线束硬化伪影的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110811660B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111476856B (zh) * | 2020-04-08 | 2023-06-06 | 中北大学 | 一种多谱ct成像方法 |
CN112666194B (zh) * | 2020-12-22 | 2022-12-20 | 上海培云教育科技有限公司 | 一种虚拟数字dr图像的生成方法及dr虚拟仿真仪器 |
CN114638907B (zh) * | 2022-02-23 | 2023-04-14 | 赛诺威盛科技(北京)股份有限公司 | 基于深度学习的骨硬化伪影校正方法及装置 |
CN115270075B (zh) * | 2022-08-03 | 2023-04-25 | 北京朗视仪器股份有限公司 | 基于钨靶x射线球管的双能ct系统的能谱数据校正方法 |
CN117330592B (zh) * | 2023-11-30 | 2024-03-15 | 赛诺威盛科技(北京)股份有限公司 | 能谱及散射校正方法、装置和设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05261092A (ja) * | 1992-03-18 | 1993-10-12 | Hitachi Medical Corp | ビームハードニング補正機能を持つx線ct装置 |
CN103961125A (zh) * | 2013-01-31 | 2014-08-06 | 东北大学 | 一种用于锥束ct的ct值校正方法 |
CN106618619A (zh) * | 2016-01-30 | 2017-05-10 | 上海联影医疗科技有限公司 | 计算机断层成像设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6904120B2 (en) * | 2003-07-01 | 2005-06-07 | General Electric Company | Method and apparatus for correcting bone induced spectral artifacts |
WO2017128892A1 (zh) * | 2016-01-30 | 2017-08-03 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及系统 |
-
2019
- 2019-10-25 CN CN201911022549.XA patent/CN110811660B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05261092A (ja) * | 1992-03-18 | 1993-10-12 | Hitachi Medical Corp | ビームハードニング補正機能を持つx線ct装置 |
CN103961125A (zh) * | 2013-01-31 | 2014-08-06 | 东北大学 | 一种用于锥束ct的ct值校正方法 |
CN106618619A (zh) * | 2016-01-30 | 2017-05-10 | 上海联影医疗科技有限公司 | 计算机断层成像设备 |
CN106725569A (zh) * | 2016-01-30 | 2017-05-31 | 上海联影医疗科技有限公司 | 骨硬化伪影校正方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110811660A (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110811660B (zh) | 一种校正ct射线束硬化伪影的方法 | |
US9025815B2 (en) | System and method for multi-material correction of image data | |
JP6492005B2 (ja) | X線ct装置、再構成演算装置、及び再構成演算方法 | |
US9498179B1 (en) | Methods and systems for metal artifact reduction in spectral CT imaging | |
EP2433566B1 (en) | System for blood vessel stenosis visualization and quantification using spectral CT analysis | |
US7391844B2 (en) | Method and apparatus for correcting for beam hardening in CT images | |
US6891918B2 (en) | Methods and apparatus for acquiring perfusion data | |
US7272429B2 (en) | Methods and apparatus for facilitating a reduction in artifacts | |
US8023767B1 (en) | Method and apparatus for 3D metal and high-density artifact correction for cone-beam and fan-beam CT imaging | |
US7822172B2 (en) | Method for hardening correction in medical imaging | |
CN107004284B (zh) | 用于光谱ct成像的方法和系统 | |
Jaju et al. | Artefacts in cone beam CT | |
CN111656405A (zh) | 使用深度学习来减少金属伪影 | |
US11419566B2 (en) | Systems and methods for improving image quality with three-dimensional scout | |
US9374565B2 (en) | System and method for multi-material correction of image data | |
JP2014061274A (ja) | 医用画像処理装置及びx線コンピュータ断層撮影装置 | |
EP2313865A2 (en) | System and method for spectral x-ray imaging | |
CN110458913B (zh) | 一种多阈值分割ct图像校正图像重建中骨硬化伪影的方法 | |
CN110458762B (zh) | 一种基于可调双因子的ct图像射束硬化伪影校正系统 | |
WO2019040930A1 (en) | LOCAL X-RAY BONE DENSITY TOMOGRAPHY METHOD | |
Molloi et al. | Scatter-glare estimation for digital radiographic systems: comparison of digital filtration and sampling techniques | |
Dinten et al. | Dual-energy x-ray absorptiometry using 2D digital radiography detector: application to bone densitometry | |
Xiu et al. | An Innovative Beam Hardening Correction Method for Computed Tomography Systems. | |
CN115067981A (zh) | 双能ct物质分解方法、装置、系统、电子设备及存储介质 | |
Njiti et al. | Potential of Metal Artifact Reduction (MAR) and Deep Learning-based Reconstruction (DLR) integration in CT Metal Artifact Correction: A review |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100176 1st floor, building 8, 11 Kangding street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Applicant after: Sinovision Technology (Beijing) Co.,Ltd. Address before: 100176 1st floor, building 8, 11 Kangding street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Applicant before: SAINUO WEISHENG TECHNOLOGY (BEIJING) Co.,Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |