CN106952285B - 基于先验统计运动模型及自适应配准的肺部运动估计方法 - Google Patents

基于先验统计运动模型及自适应配准的肺部运动估计方法 Download PDF

Info

Publication number
CN106952285B
CN106952285B CN201710080490.4A CN201710080490A CN106952285B CN 106952285 B CN106952285 B CN 106952285B CN 201710080490 A CN201710080490 A CN 201710080490A CN 106952285 B CN106952285 B CN 106952285B
Authority
CN
China
Prior art keywords
motion
lung
group
samples
motion information
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
CN201710080490.4A
Other languages
English (en)
Other versions
CN106952285A (zh
Inventor
顾力栩
陈栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Inner Mongolia Jingmai Medical Technology Co.,Ltd.
Jingmai Medical Technology Nantong Co ltd
Shanghai Jingmai Medical Technology Co ltd
Shanghai Robo Robot Technology Co ltd
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201710080490.4A priority Critical patent/CN106952285B/zh
Publication of CN106952285A publication Critical patent/CN106952285A/zh
Application granted granted Critical
Publication of CN106952285B publication Critical patent/CN106952285B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30061—Lung

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明涉及一种基于先验统计运动模型及自适应配准的肺部运动估计方法，包括：S1：获取多组病人呼吸始末两个相位的胸部CT图像；S2：对其进行肺部组织的图像分割；S3：将每组病人两个不同呼吸相位的胸部CT图像进行配准；S4：提取肺实质的运动信息，并以此建立肺部呼吸运动样本库，每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像做配准，获取样本间的转换矩阵；S5：采集一组新增肺部运动模型EI时刻的全局肺区CT及一组EE时刻的肺部病灶区域CT；S6：在新增的两种CT中选取十组的特征点对；S7：稀疏先验运动模型的生成；S8：病灶区域的精确运动信息获取。与现有技术相比，本发明具有计算速度快、实现方便、应用灵活等优点。

Description

基于先验统计运动模型及自适应配准的肺部运动估计方法
技术领域
本发明涉及肺部呼吸运动信息获取领域，尤其是涉及一种基于先验统计运动模型及自适应配准的肺部运动估计方法。
背景技术
精确治疗是当今医学的重要发展方向，是高速发展的计算技术、影像学技术与医学交叉发展的产物，可以为病人提供有效并且附加伤害最小的治疗，最大程度地降低风险，改善病人预后。肺癌的诊断率也随之有了一定的提高。目前，大多数的穿刺手术是基于传统三维CT图像引导下进行的，也是诊断及鉴别肺癌的一种新兴的有效方法，其可以实现微创下直接获取病变标本，并且将现代医学影像学技术与病理学相结合，做出组织病理学诊断，具有诊断率高，创伤小等特点。但由于人体解剖结构的动态特性，例如呼吸运动、胃肠道蠕动的变化等会导致病人器官和胸部、腹部等部位的肿瘤位置在治疗过程中发生运动。而三维CT图像是在呼吸周期中瞬时扫描所得，是静态的，基本不包含运动信息。为此穿刺活检往往需要重复多次扫描引导穿刺针进入并观察有无并发症，辐射剂量大，一定程度上限制了应用。因此，建立一种安全，高效，精确可靠并且实时的肺部肿瘤运动获取技术具有重要意义。
经过对现有文献的检索发现，现有的肺部运动模型按其建立方法主要分为三大类：基于图像灰度的数学模型，生物力学模型和统计学模型。基于图像的方法主要是运用图像配准的技术，通过使两幅或多幅图像之间达到最大相似度，从而得到一个肺部运动的流场。假定序列图像中相邻图像间的时间间隔很短，同时相邻图像间的差异也很小，通过构建一个包含时间和空间梯度的信息函数达到配准的目的。这类方法在建立图像与运动数学关系时，需要多个时刻肺部的CT数据，甚至需要采集患者的4DCT才可以完成上述目的。这势必对患者产生了大量的辐射量。并且由于图像配准要处理很大的计算量，因此，它的实时性就会收到制约。基于生物力学模型是把两个时刻肺部的模型从CT数据中重建出，然后将其中一个时刻模型上的顶点和三角面移动到另一个时刻模型上对应的位置实现呼吸模拟。基于生物力学模型的方法通过建立能量方程的形式代替相似度作为评判肺部呼吸运动的依据，然而复杂的组织间挤压情况，软组织的各向异性和非均一性都使得这类方法在获取生物属性参数及边界条件上存在一定困难，无法快速生成各异性的力学模型。
本发明要解决的关键技术问题有：
1、利用不同时刻CT图像信息，结合配准方法得到多组样本的运动信息，以此建立肺部呼吸运动模型的数据库。
2、针对当前运动统计模型无法体现个体差异或保留局部呼吸运动细节的不足，本研究将对统计模型的建立进行探索。
3、针对个体化统计学模型生后存在的估计误差，结合二次局部CT图像信息对感兴趣区域的运动进行进一步修正，以提高统计学模型在肺部运动获取中的精度。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于先验统计运动模型及自适应配准的肺部运动估计方法。
本发明的目的可以通过以下技术方案来实现：
一种基于先验统计运动模型及自适应配准的肺部运动估计方法，包括依次执行的以下步骤：
S1：事先获取多组病人呼吸始末两个相位的胸部CT图像，其中两个相位分别为吸气末EI、呼气末EE；
S2：处理从步骤S1中获取的胸部CT图像，对其进行肺部组织的图像分割，得到肺实质；
S3：将每组病人两个不同呼吸相位的胸部CT图像进行配准；
S4：根据配准后的序列图像信息，提取肺实质的运动信息，并以此建立肺部呼吸运动样本库，每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像做配准，获取样本间的转换矩阵；
S5：采集一组新增肺部运动模型EI时刻的全局肺区CT及一组EE时刻的肺部病灶区域CT；
S6：在新增的两种CT中选取十组的特征点对，获取这些特征点对的运动信息作为步骤S7的输入；
S7：稀疏先验运动模型的生成：将步骤S5中EI时刻的CT与运动样本库中的第一组CT数据进行配准，得到两者之间的转换矩阵，结合步骤S6中的特征点对运动信息获取其在运动样本库中各个样本中的各异性表达，计算获得步骤S6提供的运动信息在运动样本库中的稀疏线性表达，并以此生成肺实质的运动信息；
S8：病灶区域的精确运动信息获取：统计运动模型提供的先验运动信息作为依据，结合步骤S5中的两组图像利用配准算法得到病灶区域的精确运动信息。
优选地，所述的步骤S3中配准公式为：
式中，Ω为图像区域；Ip,EE和Ip,EI为样本库中每个病例在呼气末和吸气末的CT图像；xi为图像像素的位置；d为像素的运动向量；|Δd|为运动向量的梯度，α为正则化系数。
优选地，所述的步骤S4中样本间的对应关系获取描述为：
采用S3中的配准方法对每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像进行配准，获取样本间的转换矩阵得到运动样本库中每组样本Ip,EI对应第一组样本EI时刻CT图像I1,EI的关系为：
优选地，所述的步骤S7中稀疏先验运动模型的生成表示为：
其中yj为步骤S6中特征点对的运动信息；M为这些特征点在样本库每组样本中对应位置的运动信息，其中k为特征点对数量，n为样本数量；x为一组样本线性表达的系数；e为输入yj的误差；λ1和λ2控制x和e的稀疏量。
优选地，所述的步骤S8中病灶区域的精确运动信息获取表示为：
其中，Ω为图像区域；为B样条配准中控制点影响系数；Ij,EE和Ij,EI为当前病例在呼气末和吸气末的CT图像；gprior为步骤S7提供的先验统计运动信息；g(X)为B样条控制点的运动量，gn+1(X)和gn(X)为第n+1次和第n次配准时控制点的运动量；λ为控制点运动量正则化系数。
优选地，所述的步骤S8中B样条配准中控制点影响系数的表达公式为：
其中，Ri,j,k为第ijk个控制点在EI时刻图像上所影响的区域；为其在EE时刻图像上所影响的区域，为ξi,j,k的集合。
与现有技术相比，本发明具有以下优点：
1、采用配准多组病例不同时刻CT图像的方法获取肺部呼吸运动信息，从而建立运动数据库。该运动数据库能够较好的描述不同呼吸模式下肺部运动信息，解决了肺部运动的各向异性及非均一性问题所引起的问题；
2、稀疏先验运动模型有效地反映了新增病例的肺部局部呼吸运动细节；
3、针对稀疏先验运动模型不能提供足够精确的运动信息，本发明基于先验统计运动模型及自适应配准的肺部运动估计方法实现了对肺实质病灶组织运动的精确计算方法。
附图说明
图1为本发明的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应属于本发明保护的范围。
在一个具体实施中，采用CPU为Xeon E5405 2.0GHZ，显卡为NVIDIA Quadro VX200,内存为4.0GB的计算机配置，实现方式为MATLAB，
一种基于先验统计运动模型及自适应配准的肺部运动估计方法，包括依次执行的以下步骤：
S1：事先获取多组病人呼吸始末两个相位(吸气末EI、呼气末EE)的胸部CT图像；
S2：处理从步骤S1中获取的胸部CT图像，对其进行肺部组织的图像分割，得到肺实质；
S3：将每组病人两个不同呼吸相位的胸部CT图像进行配准；
S4：根据配准后的序列图像信息，提取肺实质的运动信息，并以此建立肺部呼吸运动样本库。每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像做配准，获取样本间的转换矩阵；
S5：采集一组新增肺部运动模型EI时刻的全局肺区CT及一组EE时刻的肺部病灶区域CT；
S6：在新增的两种CT中选取约十组的特征点对，获取这些特征点对的运动信息作为步骤S7的输入；
S7：稀疏先验运动模型的生成：将步骤S5中EI时刻的CT与运动样本库中的第一组CT数据进行配准，得到两者之间的转换矩阵。结合步骤S6中的特征点对运动信息获取其在运动样本库中各个样本中的各异性表达。计算获得步骤S6提供的运动信息在运动样本库中的稀疏线性表达，并以此生成肺实质的运动信息；
S8：病灶区域的精确运动信息获取：统计运动模型提供的先验运动信息作为依据，结合第五步中的两组图像利用配准算法得到病灶区域的精确运动信息。
在一个优选实施例中，步骤S3中配准公式为：
式中，Ω为图像区域；Ip,EE和Ip,EI为样本库中每个病例在呼气末和吸气末的CT图像；xi为图像像素的位置；d为像素的运动向量；|Δd|为运动向量的梯度，α为正则化系数。
在一个优选实施例中，步骤S4中样本间的对应关系获取可描述为：
采用S3中的配准方法对每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像进行配准，获取样本间的转换矩阵得到运动样本库中每组样本对应第一组样本EI时刻CT图像的关系为：
在一个优选实施例中，步骤S7中稀疏先验运动模型的生成可表示为：
其中,yj为步骤S6中特征点对的运动信息；M为这些特征点在样本库每组样本中对应位置的运动信息，其中k为特征点对数量，n为样本数量；x为一组样本线性表达的系数；e为输入yj的误差；λ1和λ2控制x和e的稀疏量。
在一个优选实施例中，步骤S8中病灶区域的精确运动信息获取可表示为：
其中，Ω为图像区域；为B样条配准中控制点影响系数；Ij,EE和Ij,EI为当前病例在呼气末和吸气末的CT图像；gprior为步骤S7提供的先验统计运动信息；g(X)为B样条控制点的运动量，gn+1(X)和gn(X)为第n+1次和第n次配准时控制点的运动量；λ为控制点运动量正则化系数。
在一个优选实施例中，步骤S8中病灶区域的精确运动信息获取可表示为：步骤S8中B样条配准中控制点影响系数的表达公式为：
其中，Ri，j，k为第ijk个控制点在EI时刻图像上所影响的区域；为其在EE时刻图像上所影响的区域。为ξi,j,k的集合。
上述所有的公式或方程仅属于本发明实施中的一部分例子，本领域技术人员也可以采用其他现有的公式或方程进行计算，但均应落入本发明保护的范围之内。
上述列举的各种实施例，在不矛盾的前提下，可以相互组合实施，本领域技术人员可结合附图和上文对实施例的解释，作为对不同实施例中的技术特征进行组合的依据。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (4)

1.一种肺部运动信息的计算机获取方法，其特征在于，该方法采用计算机，实现方式为MATLAB，所述的方法包括依次执行的以下步骤：
S1：事先获取多组病人呼吸始末两个相位的胸部CT图像，其中两个相位分别为吸气末EI、呼气末EE；
S2：处理从步骤S1中获取的胸部CT图像，对其进行肺部组织的图像分割，得到肺实质；
S3：将每组病人两个不同呼吸相位的胸部CT图像进行配准；
S4：根据配准后的序列图像信息，提取肺实质的运动信息，并以此建立肺部呼吸运动样本库，每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像做配准，获取样本间的转换矩阵；
S5：采集一组新增肺部运动模型EI时刻的全局肺区CT及一组EE时刻的肺部病灶区域CT；
S6：在新增的两种CT中选取十组的特征点对，获取这些特征点对的运动信息作为步骤S7的输入；
S7：稀疏先验运动模型的生成：将步骤S5中EI时刻的CT与运动样本库中的第一组CT数据进行配准，得到两者之间的转换矩阵，结合步骤S6中的特征点对运动信息获取其在运动样本库中各个样本中的各异性表达，计算获得步骤S6提供的运动信息在运动样本库中的稀疏线性表达，并以此生成肺实质的运动信息；
S8：病灶区域的精确运动信息获取：稀疏先验运动模型提供的先验运动信息作为依据，结合步骤S5中的两组图像利用配准算法得到病灶区域的精确运动信息；
所述的步骤S7中稀疏先验运动模型的生成表示为：
其中yj为步骤S6中特征点对的运动信息；M为这些特征点在样本库每组样本中对应位置的运动信息，其中k为特征点对数量，n为样本数量；x为一组样本线性表达的系数；e为输入yj的误差；λ1和λ2控制x和e的稀疏量；
所述的步骤S8中病灶区域的精确运动信息获取表示为：
其中，Ω为图像区域；为B样条配准中控制点影响系数；Ij，EE和Ij，EI为当前病例在呼气末和吸气末的CT图像；gprior为稀疏先验运动模型提供的先验统计运动信息；g(X)为B样条控制点的运动量，gn+1(X)和gn(X)为第n+1次和第n次配准时控制点的运动量；λ为控制点运动量正则化系数。
2.根据权利要求1所述的一种肺部运动信息的计算机获取方法，其特征在于，所述的步骤S3中配准所需的能量公式E(d)为：
式中，Ω为图像区域；Ip，EE和Ip，EI为样本库中每个病例在呼气末和吸气末的CT图像；xi为图像像素的位置；d为像素的运动向量；|Δd|为运动向量的梯度，α为正则化系数。
3.根据权利要求1所述的一种肺部运动信息的计算机获取方法，其特征在于，所述的步骤S4中样本间的对应关系获取描述为：
采用S3中的配准方法对每组样本的EI时刻CT图像都与第一组样本的EI时刻CT图像进行配准，获取样本间的转换矩阵得到运动样本库中每组样本Ip，EI对应第一组样本EI时刻CT图像I1，EI的关系为：
4.根据权利要求1所述的一种肺部运动信息的计算机获取方法，其特征在于，所述的步骤S8中B样条配准中控制点影响系数的表达公式为：
其中，Ri，j，k为第ijk个控制点在EI时刻图像上所影响的区域；为其在EE时刻图像上所影响的区域，为ξi，j，k的集合。
CN201710080490.4A 2017-02-15 2017-02-15 基于先验统计运动模型及自适应配准的肺部运动估计方法 Active CN106952285B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710080490.4A CN106952285B (zh) 2017-02-15 2017-02-15 基于先验统计运动模型及自适应配准的肺部运动估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710080490.4A CN106952285B (zh) 2017-02-15 2017-02-15 基于先验统计运动模型及自适应配准的肺部运动估计方法

Publications (2)

Publication Number Publication Date
CN106952285A CN106952285A (zh) 2017-07-14
CN106952285B true CN106952285B (zh) 2019-12-10

Family

ID=59467340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710080490.4A Active CN106952285B (zh) 2017-02-15 2017-02-15 基于先验统计运动模型及自适应配准的肺部运动估计方法

Country Status (1)

Country Link
CN (1) CN106952285B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110335358B (zh) * 2019-06-18 2021-11-05 大连理工大学 可变形数字人解剖学模型的个性化变形方法
CN111724364B (zh) * 2020-06-12 2023-04-14 深圳技术大学 基于肺叶及气管树的方法及装置、电子设备和存储介质
CN112785588B (zh) * 2021-02-05 2021-09-14 南京钺曦医疗科技有限公司 一种ct与mr脑灌注数据的运动幅度自动估计方法
CN113538481B (zh) * 2021-06-18 2023-06-23 深圳技术大学 肺轮廓特征点的确定方法及装置、电子设备和存储介质
CN114176726B (zh) * 2021-12-10 2023-08-04 佗道医疗科技有限公司 一种基于相位配准的穿刺方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8467497B2 (en) * 2007-10-25 2013-06-18 Tomotherapy Incorporated System and method for motion adaptive optimization for radiation therapy delivery
CN101623198A (zh) * 2008-07-08 2010-01-13 深圳市海博科技有限公司 动态肿瘤实时跟踪方法
CN102292743B (zh) * 2009-01-22 2015-01-14 皇家飞利浦电子股份有限公司 核图像重建
US8923590B2 (en) * 2011-01-20 2014-12-30 Siemens Aktiengesellschaft Method and system for 3D cardiac motion estimation from single scan of C-arm angiography
US9020229B2 (en) * 2011-05-13 2015-04-28 Broncus Medical, Inc. Surgical assistance planning method using lung motion analysis
CN102411781B (zh) * 2011-09-09 2013-11-13 华南理工大学 一种双能减影胸部x射线图像的运动校正系统
CN103761745B (zh) * 2013-07-31 2017-04-12 深圳大学 一种肺部运动模型估计方法及系统
JP6237353B2 (ja) * 2014-03-07 2017-11-29 コニカミノルタ株式会社 画像処理装置、画像撮影システム及び画像処理プログラム
CN104361568A (zh) * 2014-09-18 2015-02-18 南方医科大学 基于配准的肺4d-ct图像呼气过程中间相位图像的重建方法
CN104361606B (zh) * 2014-09-29 2017-07-11 北京理工大学 X射线造影图像序列中的心动周期恢复方法
CN105574882B (zh) * 2015-12-30 2018-08-10 中国科学院深圳先进技术研究院 基于胸部横断面ct图像的肺部分割提取方法及系统

Also Published As

Publication number Publication date
CN106952285A (zh) 2017-07-14

Similar Documents

Publication Publication Date Title
CN106952285B (zh) 基于先验统计运动模型及自适应配准的肺部运动估计方法
CN112885453B (zh) 用于标识后续医学图像中的病理变化的方法和系统
US8666128B2 (en) Methods, systems, and computer readable media for mapping regions in a model of an object comprising an anatomical structure from one image data set to images used in a diagnostic or therapeutic intervention
US8787648B2 (en) CT surrogate by auto-segmentation of magnetic resonance images
CN109598722B (zh) 基于递归神经网络的图像分析方法
RU2711140C2 (ru) Редактирование медицинских изображений
Tian et al. Medical image processing and analysis
JP2005078176A (ja) 複数画像間の非剛体レジストレーション方法
CN116797612B (zh) 基于弱监督深度活动轮廓模型的超声图像分割方法及装置
Ho et al. An unsupervised image registration method employing chest computed tomography images and deep neural networks
CN115830016A (zh) 医学图像配准模型训练方法及设备
Li et al. Pulmonary CT image registration and warping for tracking tissue deformation during the respiratory cycle through 3D consistent image registration
Liang et al. ORRN: An ODE-based recursive registration network for deformable respiratory motion estimation with lung 4DCT images
Xue et al. Lung respiratory motion estimation based on fast Kalman filtering and 4D CT image registration
Cai et al. Registration on DCE-MRI images via multi-domain image-to-image translation
CN106446572B (zh) 基于边界元模型和局部区域修正的肺部呼吸运动获取方法
Zhao et al. Target organ non-rigid registration on abdominal CT images via deep-learning based detection
Wu et al. Reconstruction of 4D-CT from a single free-breathing 3D-CT by spatial-temporal image registration
Liu et al. Real-time estimation of lung deformation from body surface using a general CoordConv CNN
JP5068334B2 (ja) 医用画像変換装置および方法並びにプログラム
JP2022052210A (ja) 情報処理装置、情報処理方法及びプログラム
CN111127636A (zh) 一种智能的复杂关节内骨折桌面级三维诊断系统
Prabu et al. A statistical deformation model (SDM) based regularizer for non-rigid image registration: application to registration of multimodal prostate MRI and histology
Chambon et al. CT-PET landmark-based lung registration using a dynamic breathing model
Chen et al. A statistical weighted sparse-based local lung motion modelling approach for model-driven lung biopsy

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200827

Address after: 226014 Jiangsu province Nantong City Xinsheng Road No. 158 Building 2 Room 502

Patentee after: JINGMAI MEDICAL TECHNOLOGY NANTONG Co.,Ltd.

Address before: 200240 Dongchuan Road, Shanghai, No. 800, No.

Patentee before: SHANGHAI JIAO TONG University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221101

Address after: Room 409, South Building, building 1, 951 Jianchuan Road, Minhang District, Shanghai 201100

Patentee after: SHANGHAI JINGMAI MEDICAL TECHNOLOGY CO.,LTD.

Patentee after: JINGMAI MEDICAL TECHNOLOGY NANTONG CO.,LTD.

Patentee after: Shanghai Robo Robot Technology Co.,Ltd.

Patentee after: Inner Mongolia Jingmai Medical Technology Co.,Ltd.

Address before: Room 502, Building 158, Xinsheng Road, Nantong City, Jiangsu Province

Patentee before: JINGMAI MEDICAL TECHNOLOGY NANTONG CO.,LTD.