CN111950639B - 一种同步实时显示超声与断层解剖图像的成像方法 - Google Patents
一种同步实时显示超声与断层解剖图像的成像方法 Download PDFInfo
- Publication number
- CN111950639B CN111950639B CN202010819759.8A CN202010819759A CN111950639B CN 111950639 B CN111950639 B CN 111950639B CN 202010819759 A CN202010819759 A CN 202010819759A CN 111950639 B CN111950639 B CN 111950639B
- Authority
- CN
- China
- Prior art keywords
- human body
- digital
- image
- imaging method
- coordinate system
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 17
- 230000006698 induction Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 210000000988 bone and bone Anatomy 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 238000002604 ultrasonography Methods 0.000 claims description 3
- 230000005674 electromagnetic induction Effects 0.000 claims description 2
- 238000013499 data model Methods 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Image Processing (AREA)
Abstract
本发明公开一种同步实时显示超声与断层解剖图像的成像方法,通过在普通的超声探头上安装3D定位装置,并在被采集人体上粘贴定位贴片,从而将被采集人体的3D位置,与数字人图像的3D位置对应起来,从而达到同时查看超声图像和对应数字人切片图像的目的。本发明能够迅速将人体超声切片与电子数据模型进行匹配,使用方便,适用装置简单。
Description
技术领域
本发明涉及数据模型领域,尤其涉及一种同步实时显示超声与断层解剖图像的成像方法。
背景技术
在目前的医疗成像技术中,具有多种将人体成像与数据模型匹配、结合的技术,进而使成像结果更加形象,方便为患者讲解病情,广泛应用于教育、研究等领域。
但是现有技术中对于人体成像与数据模型的匹配较为复杂,通常需要对图像多次处理,坐标换算较为复杂,需要设置专业的成像设备,成本高,使用复杂。
发明内容
本发明旨在提供一种同步实时显示超声与断层解剖图像的成像方法,迅速将人体超声切片与电子数据模型进行匹配,使用方便,适用装置简单。
为达到上述目的,本发明是采用以下技术方案实现的:
本发明公开一种同步实时显示超声与断层解剖图像的成像方法,包括以下步骤:
S1、构建数字人体模型,数字人体模型位于数字3D坐标系内;
S2、在被采集数据的人体表面特定位置设定感应点,收集感应点在定位装置坐标系下的3D坐标组P,记做p1,p2,...,pn;
S3、寻找数字人体模型表面对应特定位置的点的坐标组P′,记做p′1,p′2,...,p′n;
S4、使用最小二乘法,计算转换矩阵T,使得p′1=T*p1,p′2=T*p2,...,p′n=T*pn;
S5、通过定位装置坐标系中的位置s,计算得到数字3D坐标系中的位置s′,s′=T*s:
S6、获取采集数据的人体的超声图像,将超声图像与切片图像匹配,即根据数字3D坐标系中的位置s′,从数字人体模型中计算该位置的切片图像。
优选的,步骤S2中,通过在人体表面设置贴片,使用自由度为6的感应装置确定贴片坐标,感应装置通过电磁感应或者光学感应定位贴片。
优选的,被采集数据的人体所在环境设有多个坐标基准点。
优选的,步骤S2中,特定位置包括骨骼端部、肚脐,特定位置数量≥3个。
优选的,步骤S1中,数字人体模型为卧位,数字人体模型为沿身高方向设置的多个横截面图像组成,横截面图像之间的距离为d,每张横截面图中相邻像素之间的距离为q。
优选的,步骤S6中,超声图像中第i个像素的坐标为(xi,yi,zi),其对应第m张切片图像,
其中z0为数字人体模型在Z轴方向上的第一张图片,
第m张切片图像中对应像素的坐标为(x′i,y′i),
其中,(x0,y0)为切片图像中的原点坐标。
优选的,步骤S2中,人体表面为手臂部表面、手掌部表面、腿部表面、足部表面、躯干部表面中的任一。
优选的,使用计算机视觉技术,对被采集数据的人体进行人体关节位置的识别,确定人体轮廓。
本发明的有益效果:
1、本发明能够迅速将人体超声切片与电子数据模型进行匹配。
2、本发明适配方法简单,适配精度损失较小。
3、本发明在针对局部人体适配时,使用方便,速度快。
附图说明
图1为本发明数字人体模型中切片图像排列的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。
本发明包括以下步骤:
S1、构建数字人体模型,数字人体模型位于数字3D坐标系内;
S2、在被采集数据的人体表面特定位置设定感应点,收集感应点在定位装置坐标系下的3D坐标组P,记做p1,p2,...,pn;
S3、寻找数字人体模型表面对应特定位置的点的坐标组P′,记做p′1,p′2,...,p′n;
S4、使用最小二乘法,计算转换矩阵T,使得p′1=T*p1,p′2=T*p2,...,p′n=T*pn;
S5、通过定位装置坐标系中的位置s,计算得到数字3D坐标系中的位置s′,s′=T*s;
S6、获取采集数据的人体的超声图像,将超声图像与切片图像匹配,即根据数字3D坐标系中的位置s′,从数字人体模型中计算该位置的切片图像。
步骤S2中,通过在人体表面设置贴片,在超声探头上,安装3D定位装置,3D定位装置为电磁定位装置或光学定位装置,6自由度。在采集超声图像的过程中,能够实时的拿到每张图像在3D空间中的位置和角度,记作s,该位置和角度,都是相对于3D定位装置的坐标系。
步骤S2中,特定位置包括骨骼端部、肚脐,特定位置数量≥3个,不在同一条直线上的3个点可以实现标定。标定是指尽量将数字人模型与真实人体的位置进行匹配。比如,我们选择的三个点是头顶上一个点,两只手的中止指尖各一个点,那么我们所谓的“标定”,就是将真人的头顶这个点的坐标,与数字人的头顶点坐标尽量重合,真人的两个指尖的点与数字人两个指尖的点尽量重合,由于人体形状具有差异,不可能完全重合,完成标定后,3个点就已经固定了人体的位置了。当然,选取的点越多,通过取最小方差的形式,使得标定的平均误差越小。
如图1所示,步骤S1中,数字人体模型为卧位,数字人体模型为沿身高方向设置的多个横截面图像组成,横截面图像之间的距离为d,每张横截面图中相邻像素之间的距离为q。
步骤S6中,超声图像中第i个像素的坐标为(xi,yi,zi),其对应第m张切片图像,
其中z0为数字人体模型在Z轴方向上的第一张图片,
第m张切片图像中对应像素的坐标为(x′i,y′i),
其中,(x0,y0)为切片图像中的原点坐标。
步骤S2中,人体表面为手臂部表面、手掌部表面、腿部表面、足部表面、躯干部表面中的任一,如果要做全身的匹配,我们需要要求用户的姿势与数字人尽量保持一致。我们的数字人是平躺,那么用户的姿势最好也是平躺,如果用户是侧躺,则匹配度会不好,如果是蜷缩,则完全不能做到全身匹配。但是,在实际应用中,我们往往不需要全身匹配,只需对人体的某个部位进行匹配就可以,那么我们只需要保证所要匹配的部位,与数字人的对应部位保持姿态一样就可以了。比如我们只关注小腿,那么用户平躺或侧躺都可以进行匹配。
使用计算机视觉技术,对被采集数据的人体进行人体关节位置的识别,确定人体轮廓。也就是说,病人(或模特)躺在床上后,摄像头立刻可以找到人体各个主要关节的位置,从而对人体有了一个大致的定位,而后,再通过贴片+3D定位器的方式,进行一次较细致的定位,例如intel skeleton tracking sdk就有这样的功能,其介绍可参考https://www.intelrealsense.com/skeleton-tracking/。
当然,本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (6)
1.一种同步实时显示超声与断层解剖图像的成像方法,其特征在于包括以下步骤:
S1、构建数字人体模型,数字人体模型位于数字3D坐标系内;数字人体模型为卧位,数字人体模型为沿身高方向设置的多个横截面图像组成,横截面图像之间的距离为d,每张横截面图中相邻像素之间的距离为q;
S2、在被采集数据的人体表面特定位置设定感应点,收集感应点在定位装置坐标系下的3D坐标组P,记做p1,p2,…,pn;
S3、寻找数字人体模型表面对应特定位置的点的坐标组P′,记做p′1,p′2,…,p′n;
S4、使用最小二乘法,计算转换矩阵T,使得p′1=T*p1,p′2=T*p2,…,p′n=T*pn;
S5、通过定位装置坐标系中的位置s,计算得到数字3D坐标系中的位置s′,s′=T*s;
S6、获取采集数据的人体的超声图像,将超声图像与切片图像匹配,即根据数字3D坐标系中的位置s′,从数字人体模型中计算该位置的切片图像;超声图像中第i个像素的坐标为(xi,yi,zi),其对应第m张切片图像,
其中z0为数字人体模型在Z轴方向上的第一张图片,
第m张切片图像中对应像素的坐标为(x′i,y′i),
其中,(x0,y0)为切片图像中的原点坐标。
2.根据权利要求1所述的成像方法,其特征在于:步骤S2中,通过在人体表面设置贴片,使用自由度为6的感应装置确定贴片坐标,感应装置通过电磁感应或者光学感应定位贴片。
3.根据权利要求2所述的成像方法,其特征在于:被采集数据的人体所在环境设有多个坐标基准点。
4.根据权利要求1所述的成像方法,其特征在于:步骤S2中,特定位置包括骨骼端部、肚脐,特定位置数量≥3个。
5.根据权利要求1所述的成像方法,其特征在于:步骤S2中,人体表面为手臂部表面、手掌部表面、腿部表面、足部表面、躯干部表面中的任一。
6.根据权利要求1所述的成像方法,其特征在于:进行步骤S2前,使用计算机视觉技术,对被采集数据的人体进行人体关节位置的识别,确定人体轮廓。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010819759.8A CN111950639B (zh) | 2020-08-14 | 2020-08-14 | 一种同步实时显示超声与断层解剖图像的成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010819759.8A CN111950639B (zh) | 2020-08-14 | 2020-08-14 | 一种同步实时显示超声与断层解剖图像的成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111950639A CN111950639A (zh) | 2020-11-17 |
CN111950639B true CN111950639B (zh) | 2024-03-19 |
Family
ID=73342409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010819759.8A Active CN111950639B (zh) | 2020-08-14 | 2020-08-14 | 一种同步实时显示超声与断层解剖图像的成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111950639B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112991854A (zh) * | 2021-02-05 | 2021-06-18 | 四川大学华西医院 | 一种超声教学方法、装置及系统、电子设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2600887Y (zh) * | 2003-01-30 | 2004-01-21 | 张导华 | 一种人体仿真模型 |
CN103040525A (zh) * | 2012-12-27 | 2013-04-17 | 深圳先进技术研究院 | 一种多模医学影像手术导航方法及系统 |
CN104027106A (zh) * | 2014-05-20 | 2014-09-10 | 武汉培威医学科技有限公司 | 一种心电断层扫描成像系统及方法 |
CN107714082A (zh) * | 2017-09-04 | 2018-02-23 | 北京航空航天大学 | 一种基于光学定位的超声探头标定装置及方法 |
CN107909624A (zh) * | 2017-12-05 | 2018-04-13 | 南京大学 | 一种从三维断层成像中提取及融合二维图像的方法 |
CN109770943A (zh) * | 2019-01-28 | 2019-05-21 | 电子科技大学 | 一种利用计算机视觉定位的超声自动优化方法 |
CN109938768A (zh) * | 2019-03-11 | 2019-06-28 | 深圳市比邻星精密技术有限公司 | 超声成像方法、装置、计算机设备及存储介质 |
CN110148208A (zh) * | 2019-04-03 | 2019-08-20 | 中国人民解放军陆军军医大学 | 一种基于中国数字人的鼻咽部放疗教学模型构建方法 |
CN110251231A (zh) * | 2019-06-13 | 2019-09-20 | 艾瑞迈迪科技石家庄有限公司 | 超声三维重建的方法及装置 |
CN110400253A (zh) * | 2019-07-02 | 2019-11-01 | 西安工业大学 | 一种基于双线性插值原理确定发射层析权重矩阵的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005055665B4 (de) * | 2005-11-22 | 2013-11-07 | Siemens Aktiengesellschaft | Verfahren zur Interpolation von Volumendaten |
-
2020
- 2020-08-14 CN CN202010819759.8A patent/CN111950639B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2600887Y (zh) * | 2003-01-30 | 2004-01-21 | 张导华 | 一种人体仿真模型 |
CN103040525A (zh) * | 2012-12-27 | 2013-04-17 | 深圳先进技术研究院 | 一种多模医学影像手术导航方法及系统 |
CN104027106A (zh) * | 2014-05-20 | 2014-09-10 | 武汉培威医学科技有限公司 | 一种心电断层扫描成像系统及方法 |
CN107714082A (zh) * | 2017-09-04 | 2018-02-23 | 北京航空航天大学 | 一种基于光学定位的超声探头标定装置及方法 |
CN107909624A (zh) * | 2017-12-05 | 2018-04-13 | 南京大学 | 一种从三维断层成像中提取及融合二维图像的方法 |
CN109770943A (zh) * | 2019-01-28 | 2019-05-21 | 电子科技大学 | 一种利用计算机视觉定位的超声自动优化方法 |
CN109938768A (zh) * | 2019-03-11 | 2019-06-28 | 深圳市比邻星精密技术有限公司 | 超声成像方法、装置、计算机设备及存储介质 |
CN110148208A (zh) * | 2019-04-03 | 2019-08-20 | 中国人民解放军陆军军医大学 | 一种基于中国数字人的鼻咽部放疗教学模型构建方法 |
CN110251231A (zh) * | 2019-06-13 | 2019-09-20 | 艾瑞迈迪科技石家庄有限公司 | 超声三维重建的方法及装置 |
CN110400253A (zh) * | 2019-07-02 | 2019-11-01 | 西安工业大学 | 一种基于双线性插值原理确定发射层析权重矩阵的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111950639A (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Robotic arm based automatic ultrasound scanning for three-dimensional imaging | |
US10507002B2 (en) | X-ray system and method for standing subject | |
Mundermann et al. | Accurately measuring human movement using articulated ICP with soft-joint constraints and a repository of articulated models | |
US10881353B2 (en) | Machine-guided imaging techniques | |
EP2953569B1 (en) | Tracking apparatus for tracking an object with respect to a body | |
ES2550455T3 (es) | Método y aparato para corregir un error en el registro conjunto de sistemas de coordenadas usados para representar objetos visualizados durante la estimulación guiada del cerebro | |
Yang et al. | Automatic 3-D imaging and measurement of human spines with a robotic ultrasound system | |
US20100319100A1 (en) | Simple techniques for three-dimensional modeling | |
US20050182341A1 (en) | Posture diagnosis equipment and program therefor | |
US9715759B2 (en) | Reference object for three-dimensional modeling | |
Suligoj et al. | RobUSt–an autonomous robotic ultrasound system for medical imaging | |
CN111588467B (zh) | 基于医学影像的三维空间坐标转换为影像二维坐标的方法 | |
CN104771189B (zh) | 三维头颅图像对正方法及装置 | |
CN111028341A (zh) | 一种三维模型生成方法 | |
CN112183316B (zh) | 一种运动员人体姿态测量方法 | |
Bumacod et al. | Image-processing-based digital goniometer using OpenCV | |
CN111950639B (zh) | 一种同步实时显示超声与断层解剖图像的成像方法 | |
US9633433B1 (en) | Scanning system and display for aligning 3D images with each other and/or for detecting and quantifying similarities or differences between scanned images | |
EP3982287A1 (en) | Object keypoint detection | |
KR20160057024A (ko) | 마커리스 3차원 객체추적 장치 및 그 방법 | |
US20220309690A1 (en) | Method and system for non-contact patient registration in image-guided surgery | |
US20230074362A1 (en) | Method and system for non-contact patient registration in image-guided surgery | |
CN116485850A (zh) | 基于深度学习的手术导航影像的实时非刚性配准方法及系统 | |
CN114712835B (zh) | 一种基于双目人体位姿识别的辅助训练系统 | |
CN112991437B (zh) | 基于图像展缩技术的全自动穴位定位方法 |
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 |