CN105046644A - Ultrasonic and CT image registration method and system based on linear dependence - Google Patents

Ultrasonic and CT image registration method and system based on linear dependence Download PDF

Info

Publication number
CN105046644A
CN105046644A CN201510393418.8A CN201510393418A CN105046644A CN 105046644 A CN105046644 A CN 105046644A CN 201510393418 A CN201510393418 A CN 201510393418A CN 105046644 A CN105046644 A CN 105046644A
Authority
CN
China
Prior art keywords
image
ultrasound
ultrasonic
linear correlation
transformation matrix
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.)
Granted
Application number
CN201510393418.8A
Other languages
Chinese (zh)
Other versions
CN105046644B (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.)
Ka Heng Medical Technology (shanghai) Co Ltd
Original Assignee
Ka Heng Medical Technology (shanghai) Co Ltd
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 Ka Heng Medical Technology (shanghai) Co Ltd filed Critical Ka Heng Medical Technology (shanghai) Co Ltd
Priority to CN201510393418.8A priority Critical patent/CN105046644B/en
Publication of CN105046644A publication Critical patent/CN105046644A/en
Application granted granted Critical
Publication of CN105046644B publication Critical patent/CN105046644B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

The invention discloses an ultrasonic and CT image registration method and system based on linear dependence. The method comprises steps: linear dependence measurement is established; an ultrasonic and CT image registration transformation matrix is preset according to ultrasonic image texture information; the transformation matrix is optimized, and a linear dependence measurement optimal value is calculated; according to parameters of the transformation matrix, an ultrasonic image spatial position is changed, and an ultrasonic image and a CT image fuse. The system comprises a similarity measurement construction module, an original state setting module, a transformation matrix optimization module and a fusion display module. Ultrasonic and CT image registration can be completed rapidly and precisely, and theoretical guidance is provided for clinic focus diagnosis and treatment.

Description

Ultrasonic and CT image registration method and system based on linear correlation
Technical Field
The invention relates to the technical field of medicine, in particular to an ultrasonic and CT image registration method and system based on linear correlation.
Background
The image-guided surgery navigation utilizes medical imaging calculation and computer image technology to perform segmentation recognition, three-dimensional reconstruction and visualization of human tissues, organs and focuses on a scanned medical image of a patient before surgery, and a clinician can perform optimal surgery path planning and clinical surgery simulation through the acquired image; the method comprises the steps of registering a medical image acquired before an operation and a real physical space of a patient in the operation, converting the medical image before the operation, the actual body position of the patient in the operation and a surgical instrument into the same space three-dimensional coordinate system, positioning the relative position of the surgical instrument in human tissues and a focus in real time according to positioning equipment, and fusing and displaying the relationship between the medical image before the operation and the real physical space of the patient in real time. The clinician can observe the pose and various parameters (such as angle, depth and the like) of the surgical instrument in the tissues and the focus in real time from various angles through the fusion displayed image, thereby avoiding important tissues and organs (such as large blood vessels and the like) of the human body to the maximum extent and reaching the focus from the optimal surgical path in the shortest time to carry out accurate surgical treatment. The ultrasound image is a common intraoperative imaging device in clinical use at present by virtue of the advantages of real-time performance, no wound, no ionizing radiation and the like, and the ultrasound guided surgery navigation technology is rapidly developed. Due to the limitation of the ultrasonic imaging mechanism, the ultrasonic two-dimensional image has low resolution, narrow imaging field of view and contains special speckle noise, so that the position of the focus in the human body is difficult to accurately identify and judge through the ultrasonic two-dimensional image. Compared with an ultrasonic two-dimensional image, a Computed Tomography (CT) image is a three-dimensional image, has high imaging resolution and wide imaging field of view, can accurately identify and segment most focuses of human tissues and organs from the CT image, but has radiation injury to a patient on one hand and cannot be imaged in real time in an operation on the other hand when the CT image is acquired, so that the CT image cannot be used as a guide image of an intraoperative surgical navigation system. Therefore, in combination with the imaging characteristics of the ultrasonic two-dimensional image and the CT image, in the clinical application of the image-guided surgery navigation system, the ultrasonic two-dimensional image is used as an intraoperative guide image, a clinician can position a surgical instrument in real time through the ultrasonic two-dimensional image, and the intraoperative real-time ultrasonic image and the preoperative CT image are registered, so that fusion display of the ultrasonic image and the CT image is completed, and the clinician is assisted in positioning the positions of human tissue organs and focuses.
Disclosure of Invention
The technical problem to be solved by the embodiments of the present invention is to provide a method and a system for registering ultrasound and CT images based on linear correlation, which can complete accurate registration of a two-dimensional ultrasound image and a three-dimensional CT image.
In order to solve the above technical problem, an embodiment of the present invention provides a linear correlation-based ultrasound and CT image registration method, including:
establishing linear correlation measurement, and calculating the similarity between the ultrasonic image and the CT image;
presetting an ultrasonic and CT image registration transformation matrix according to ultrasonic image texture information;
optimizing the transformation matrix, and calculating to obtain an optimal value of linear correlation measure;
and changing the spatial pose of the ultrasonic image according to the transformation matrix parameters, and fusing the ultrasonic image and the CT image.
In other aspects of the invention, the linear correlation measure introduces an ultrasound reflection coefficient into the ultrasound and CT image similarity measure.
In other aspects of the present invention, the similarity between the ultrasound image and the CT image is calculated as the difference between corresponding regions of interest in the ultrasound and CT images.
In other aspects of the present invention, the ultrasound image texture information includes boundaries of human tissue and organs and vessel bifurcations in the ultrasound image.
In other aspects of the invention, the optimization of the transformation matrix is to change 7 parameter values in the transformation matrix.
In other aspects of the present invention, the optimal value of the linear correlation measure is the minimum difference value between the ultrasound image and the region of interest in the CT image after the spatial position of the ultrasound image is changed according to the transformation matrix.
In other aspects of the present invention, the fusion of the ultrasound image and the CT image means that the ultrasound image and the CT image are displayed in an overlapping manner according to a preset transparency.
The embodiment of the invention also provides an ultrasonic and CT image registration system based on linear correlation, which comprises:
the similarity measure building module is used for building linear correlation measure and calculating the similarity between the ultrasonic image and the CT image;
the initial state setting module is used for presetting an ultrasound and CT image registration transformation matrix according to the ultrasound image texture information;
the transformation matrix optimization module optimizes the transformation matrix and calculates to obtain an optimal value of linear correlation measure;
and the fusion display module changes the spatial pose of the ultrasonic image according to the transformation matrix parameters, and the ultrasonic image is fused with the CT image.
By using the invention, the precise registration of the two-dimensional ultrasonic image and the three-dimensional CT image can be completed.
Drawings
FIG. 1 is a flowchart of the present invention for a linear correlation-based ultrasound and CT image registration method;
FIG. 2 is a schematic diagram of ultrasound and CT image registration of the ultrasound and CT image registration method based on linear correlation according to the present invention;
fig. 3 is a schematic diagram of a fusion display of an ultrasound image and a CT image according to the present invention. Wherein, (a) is an ultrasonic image, (b) is a CT image, and (c) is a fusion display of the ultrasonic image and the CT image;
fig. 4 is a schematic structural diagram of an ultrasound and CT image registration system based on linear correlation according to an embodiment of the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and the examples, but without limiting the invention.
Fig. 1 is a schematic flow chart of an ultrasound and CT image registration method based on linear correlation according to an embodiment of the present invention, which includes the following specific steps:
in the step of S1,
and establishing linear correlation measurement, and calculating the similarity between the ultrasonic image and the CT image.
Defining the similarity of the ultrasonic image and the CT image according to the related ratio as follows:
<math> <mrow> <mi>C</mi> <mi>R</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <msub> <mi>&Sigma;</mi> <mrow> <mi>x</mi> <mo>&Element;</mo> <mi>&Omega;</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mi>U</mi> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>-</mo> <mi>f</mi> <mo>(</mo> <mi>&mu;</mi> <mo>(</mo> <mi>T</mi> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>)</mo> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mo>|</mo> <mi>&Omega;</mi> <mo>|</mo> <mi>V</mi> <mi>a</mi> <mi>r</mi> <mrow> <mo>(</mo> <mi>U</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> </math>
wherein, U represents an ultrasound image, μ represents a CT image, T represents an ultrasound/CT image registration matrix, f represents a mapping function between the ultrasound two-dimensional image U and the CT slice μ, and Ω represents a size of the image. Assuming that the mapping function f is a linear function, it is expressed as:
f(μ)=αμ+β
there are many large range of reflections and echoes of tissue organs in the ultrasound image, while there is no reflection and echo information in the CT image, introducing the ultrasound reflection coefficient r into the similarity evaluation criteria for ultrasound/CT image registration. Definition of
Representing the gray value of a pixel, the pixel gray value calculation function can be written as:
fxi=αpi+βri
ensuring minimization of unknown parameters alpha, beta and gamma
<math> <mrow> <mo>|</mo> <mo>|</mo> <mi>M</mi> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <mi>&alpha;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&beta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&gamma;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </math>
Wherein,
M = p 1 r 1 1 . . . . . . . . . p n r n 1
in the above, the similarity between the ultrasound image and the CT image may be to calculate a difference between corresponding regions of interest in the ultrasound and CT images. Of course, other criteria may also be based.
In the step of S2,
and presetting an ultrasound and CT image registration transformation matrix according to the ultrasound image texture information.
The ultrasound image texture information may include the boundaries of human tissues and organs and vessel bifurcation points in the ultrasound image, and is obtained by an image segmentation method.
In the step of S3,
and optimizing the transformation matrix, and calculating to obtain the optimal value of the linear correlation measure.
Defining an initial registration matrix of the ultrasonic/CT image as G0In general, G0If the evaluation standard measure is not the matrix which enables the evaluation standard measure to reach the optimal value, then the direction and the step length of the registration matrix optimization search are set, and the similarity evaluation standard of the ultrasonic/CT image is continuously changed G when the registration matrix is continuously optimized in an iteration mode1,G2,...,Gn. The registration matrix iterative optimization formula can be expressed as:
Gk=Gk+1+aksk
wherein s represents a direction vector of the registration matrix optimization search, and a represents a step size.
Decomposing the process of optimizing and calculating the registration matrix of the ultrasonic/CT image into a registration matrix iterative change process, representing each iteration of the registration matrix as a one-dimensional vector containing n +1 parameters, and in the registration matrix iterative process, using an initial registration matrix G0Optimizing a registration matrix in n directions as an origin, calculating an ultrasonic/CT image similarity evaluation standard function, selecting a registration matrix G with the optimal similarity evaluation standard, and calculating a secondary G0And continuing to iteratively change the registration matrix by taking G 'as a new origin until the registration matrix G' with the optimal similarity evaluation standard in the change of the G registration matrix. As shown in fig. 2, the result of the registration of the ultrasound and CT images is schematically shown.
The optimal value of the linear correlation measure may be that the difference between the ultrasound image and the region of interest in the CT image is the minimum after the spatial position of the ultrasound image is changed according to the transformation matrix.
In the step of S4,
and changing the spatial pose of the ultrasonic image according to the transformation matrix parameters, and fusing the ultrasonic image and the CT image.
And after the ultrasonic and CT image registration transformation matrix is obtained through optimization calculation, the ultrasonic image space pose for registration is changed according to the transformation matrix parameters. And simultaneously displaying the ultrasonic image and the CT image according to the preset transparency of the ultrasonic image and the CT image to complete the fusion display of the ultrasonic image and the CT image. As shown in fig. 3, the fusion result of the ultrasound image and the CT image is schematically shown. The fusion of the ultrasound image and the CT image may be displaying the ultrasound image and the CT image in an overlapping manner according to a preset transparency.
Fig. 4 is a schematic structural diagram of an ultrasound and CT image registration system based on linear correlation according to an embodiment of the present invention. As shown in fig. 4, the linear correlation-based ultrasound and CT image registration system includes a similarity measure construction module, an initial state setting module, a transformation matrix optimization module, and a fusion display module.
1. Similarity measure construction module
And establishing linear correlation measurement, and calculating the similarity between the ultrasonic image and the CT image.
Defining the similarity of the ultrasonic image and the CT image according to the related ratio as follows:
<math> <mrow> <mi>C</mi> <mi>R</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <msub> <mi>&Sigma;</mi> <mrow> <mi>x</mi> <mo>&Element;</mo> <mi>&Omega;</mi> </mrow> </msub> <msup> <mrow> <mo>(</mo> <mi>U</mi> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>-</mo> <mi>f</mi> <mo>(</mo> <mi>&mu;</mi> <mo>(</mo> <mi>T</mi> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>)</mo> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mo>|</mo> <mi>&Omega;</mi> <mo>|</mo> <mi>V</mi> <mi>a</mi> <mi>r</mi> <mrow> <mo>(</mo> <mi>U</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> </math>
wherein, U represents an ultrasound image, μ represents a CT image, T represents an ultrasound/CT image registration matrix, f represents a mapping function between the ultrasound two-dimensional image U and the CT slice μ, and Ω represents a size of the image. Assuming that the mapping function f is a linear function, it is expressed as:
f(μ)=αμ+β
there are many large range of reflections and echoes of tissue organs in the ultrasound image, while there is no reflection and echo information in the CT image, introducing the ultrasound reflection coefficient r into the similarity evaluation criteria for ultrasound/CT image registration. Definition of
Representing the gray value of a pixel, the pixel gray value calculation function can be written as:
fxi=αpi+βri
ensuring minimization of unknown parameters alpha, beta and gamma
<math> <mrow> <mo>|</mo> <mo>|</mo> <mi>M</mi> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <mi>&alpha;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&beta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&gamma;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </math>
Wherein,
M = p 1 r 1 1 . . . . . . . . . p n r n 1
2. initial state setting module
And presetting an ultrasound and CT image registration transformation matrix according to the ultrasound image texture information.
The ultrasonic image texture information comprises human tissue organ boundaries and blood vessel bifurcation points in the ultrasonic image, and the ultrasonic image texture information is obtained by an image segmentation method.
3. Transformation matrix optimization module
And optimizing the transformation matrix, and calculating to obtain the optimal value of the linear correlation measure.
Defining an initial registration matrix of the ultrasonic/CT image as G0In general, G0If the evaluation standard measure is not the matrix which enables the evaluation standard measure to reach the optimal value, then the direction and the step length of the registration matrix optimization search are set, and the similarity evaluation standard of the ultrasonic/CT image is continuously changed G when the registration matrix is continuously optimized in an iteration mode1,G2,...,Gn. The registration matrix iterative optimization formula can be expressed as:
Gk=Gk+1+aksk
wherein s represents a direction vector of the registration matrix optimization search, and a represents a step size.
Optimized calculation of ultrasound/CT image registration matrixThe process is decomposed into a registration matrix iterative change process, each iteration of the registration matrix is represented as a one-dimensional vector containing n +1 parameters, and an initial registration matrix G is used in the registration matrix iterative process0Optimizing a registration matrix in n directions as an origin, calculating an ultrasonic/CT image similarity evaluation standard function, selecting a registration matrix G with the optimal similarity evaluation standard, and calculating a secondary G0And continuing to iteratively change the registration matrix by taking G 'as a new origin until the registration matrix G' with the optimal similarity evaluation standard in the change of the G registration matrix. As shown in fig. 2, the result of the registration of the ultrasound and CT images is schematically shown.
4. Fusion display module
And changing the spatial pose of the ultrasonic image according to the transformation matrix parameters, and fusing the ultrasonic image and the CT image.
And after the ultrasonic and CT image registration transformation matrix is obtained through optimization calculation, the ultrasonic image space pose for registration is changed according to the transformation matrix parameters. And simultaneously displaying the ultrasonic image and the CT image according to the preset transparency of the ultrasonic image and the CT image to complete the fusion display of the ultrasonic image and the CT image. As shown in fig. 3, the fusion result of the ultrasound image and the CT image is schematically shown.
Of course, the foregoing is the preferred embodiment of the present invention. For convenience of illustration, the sequence numbers of steps S1, S2, etc. are used, but it should be appreciated that the steps themselves may include other processes, and there may be other steps between the steps, which are also within the scope of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principle of the present invention, and such improvements and modifications are also considered to be within the scope of the present invention.

Claims (10)

1. A linear correlation-based ultrasound and CT image registration method is characterized by comprising the following steps:
establishing linear correlation measurement, and calculating the similarity between the ultrasonic image and the CT image;
presetting an ultrasonic and CT image registration transformation matrix according to ultrasonic image texture information;
optimizing the transformation matrix, and calculating to obtain an optimal value of linear correlation measure;
and changing the spatial pose of the ultrasonic image according to the transformation matrix parameters, and fusing the ultrasonic image and the CT image.
2. The linear correlation based ultrasound and CT image registration method of claim 1, wherein the linear correlation measure introduces ultrasound reflection coefficients into the ultrasound and CT image similarity measure.
3. The linear correlation based ultrasound and CT image registration method according to claim 1, wherein the similarity between the ultrasound image and the CT image is to calculate a difference between corresponding regions of interest in the ultrasound and CT images.
4. The linear correlation based ultrasound and CT image registration method of claim 1, wherein the ultrasound image texture information includes human tissue organ boundaries, vessel bifurcation points in the ultrasound image.
5. The linear correlation based ultrasound and CT image registration method of claim 1, wherein the optimization of the transformation matrix is to change 7 parameter values in the transformation matrix.
6. The linear correlation based ultrasound and CT image registration method of claim 1, wherein the optimal value of the linear correlation measure is the smallest difference value between the ultrasound image and the region of interest in the CT image after the spatial position of the ultrasound image is changed according to the transformation matrix.
7. The linear correlation-based ultrasound and CT image registration method of claim 1, wherein the ultrasound image and CT image fusion means that the ultrasound image and the CT image are displayed in an overlapping manner according to a preset transparency.
8. A linear correlation based ultrasound and CT image registration system, comprising:
the similarity measure building module is used for building linear correlation measure and calculating the similarity between the ultrasonic image and the CT image;
the initial state setting module is used for presetting an ultrasound and CT image registration transformation matrix according to the ultrasound image texture information;
the transformation matrix optimization module optimizes the transformation matrix and calculates to obtain an optimal value of linear correlation measure;
and the fusion display module changes the spatial pose of the ultrasonic image according to the transformation matrix parameters, and the ultrasonic image is fused with the CT image.
9. The linear correlation based ultrasound and CT image registration system according to claim 8, wherein the linear correlation measure introduces ultrasound reflection coefficients into the ultrasound and CT image similarity measure.
10. The linear correlation based ultrasound and CT image registration system according to claim 8, wherein the similarity between the ultrasound image and the CT image is a calculation of a difference between corresponding regions of interest in the ultrasound and CT images.
CN201510393418.8A 2015-07-06 2015-07-06 Ultrasonic and CT image registration method and system based on linear correlation Active CN105046644B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510393418.8A CN105046644B (en) 2015-07-06 2015-07-06 Ultrasonic and CT image registration method and system based on linear correlation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510393418.8A CN105046644B (en) 2015-07-06 2015-07-06 Ultrasonic and CT image registration method and system based on linear correlation

Publications (2)

Publication Number Publication Date
CN105046644A true CN105046644A (en) 2015-11-11
CN105046644B CN105046644B (en) 2021-08-13

Family

ID=54453164

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510393418.8A Active CN105046644B (en) 2015-07-06 2015-07-06 Ultrasonic and CT image registration method and system based on linear correlation

Country Status (1)

Country Link
CN (1) CN105046644B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858479A (en) * 2018-07-10 2019-06-07 上海其高电子科技有限公司 The illegal whistle capturing system of motor vehicle based on image registration
CN111292248A (en) * 2018-12-10 2020-06-16 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic fusion imaging method and ultrasonic fusion navigation system
CN111481839A (en) * 2020-05-13 2020-08-04 戴建荣 Vertical adaptive radiotherapy system and method
CN112043377A (en) * 2020-10-12 2020-12-08 中国人民解放军总医院第五医学中心 Ultrasound visual field simulation auxiliary ablation path planning method and system for any section of CT
CN112155733A (en) * 2020-09-29 2021-01-01 苏州微创畅行机器人有限公司 Readable storage medium, bone modeling and registering system and bone surgery system
CN113056770A (en) * 2018-08-29 2021-06-29 新加坡科技研究局 Lesion localization in organs
CN113160221A (en) * 2021-05-14 2021-07-23 深圳市奥昇医疗科技有限责任公司 Image processing method, image processing device, computer equipment and storage medium
CN113393505A (en) * 2021-06-25 2021-09-14 浙江商汤科技开发有限公司 Image registration method, visual positioning method, related device and equipment
CN114219889A (en) * 2021-11-04 2022-03-22 华南理工大学 Spine in-vivo quasi-static three-dimensional measurement method based on CT and ultrasonic fusion
CN116580033A (en) * 2023-07-14 2023-08-11 卡本(深圳)医疗器械有限公司 Multi-mode medical image registration method based on image block similarity matching

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007100145A1 (en) * 2006-02-28 2007-09-07 Sony Corporation Registration device, correlation device, extraction method, and program
CN101305395A (en) * 2005-11-10 2008-11-12 皇家飞利浦电子股份有限公司 Adaptive point-based elastic image registration
CN103190078A (en) * 2010-09-03 2013-07-03 数字标记公司 Signal processors and methods for estimating transformations between signals
CN103295455A (en) * 2013-06-19 2013-09-11 北京理工大学 Ultrasonic training system based on CT (Computed Tomography) image simulation and positioning
CN103310458A (en) * 2013-06-19 2013-09-18 北京理工大学 Method for elastically registering medical images by aid of combined convex hull matching and multi-scale classification strategy
CN103455989A (en) * 2013-09-24 2013-12-18 南京大学 Method for improving imaging quality of limited-angle CT through combination of ultrasonic image
CN103503025A (en) * 2011-02-25 2014-01-08 弗劳恩霍夫应用研究促进协会 Determining model parameters based on transforming a model of an object

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101305395A (en) * 2005-11-10 2008-11-12 皇家飞利浦电子股份有限公司 Adaptive point-based elastic image registration
WO2007100145A1 (en) * 2006-02-28 2007-09-07 Sony Corporation Registration device, correlation device, extraction method, and program
CN103190078A (en) * 2010-09-03 2013-07-03 数字标记公司 Signal processors and methods for estimating transformations between signals
CN103503025A (en) * 2011-02-25 2014-01-08 弗劳恩霍夫应用研究促进协会 Determining model parameters based on transforming a model of an object
CN103295455A (en) * 2013-06-19 2013-09-11 北京理工大学 Ultrasonic training system based on CT (Computed Tomography) image simulation and positioning
CN103310458A (en) * 2013-06-19 2013-09-18 北京理工大学 Method for elastically registering medical images by aid of combined convex hull matching and multi-scale classification strategy
CN103455989A (en) * 2013-09-24 2013-12-18 南京大学 Method for improving imaging quality of limited-angle CT through combination of ultrasonic image

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JONGHYE WOO 等: "Multimodal Registration via Mutual Information", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
冯兆美: "基于互信息的医学图像非刚性配准研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858479A (en) * 2018-07-10 2019-06-07 上海其高电子科技有限公司 The illegal whistle capturing system of motor vehicle based on image registration
CN113056770A (en) * 2018-08-29 2021-06-29 新加坡科技研究局 Lesion localization in organs
CN111292248A (en) * 2018-12-10 2020-06-16 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic fusion imaging method and ultrasonic fusion navigation system
CN111292248B (en) * 2018-12-10 2023-12-19 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic fusion imaging method and ultrasonic fusion navigation system
CN111481839A (en) * 2020-05-13 2020-08-04 戴建荣 Vertical adaptive radiotherapy system and method
CN112155733B (en) * 2020-09-29 2022-01-28 苏州微创畅行机器人有限公司 Readable storage medium, bone modeling and registering system and bone surgery system
CN112155733A (en) * 2020-09-29 2021-01-01 苏州微创畅行机器人有限公司 Readable storage medium, bone modeling and registering system and bone surgery system
CN112043377A (en) * 2020-10-12 2020-12-08 中国人民解放军总医院第五医学中心 Ultrasound visual field simulation auxiliary ablation path planning method and system for any section of CT
CN113160221A (en) * 2021-05-14 2021-07-23 深圳市奥昇医疗科技有限责任公司 Image processing method, image processing device, computer equipment and storage medium
CN113393505B (en) * 2021-06-25 2023-11-03 浙江商汤科技开发有限公司 Image registration method, visual positioning method, related device and equipment
CN113393505A (en) * 2021-06-25 2021-09-14 浙江商汤科技开发有限公司 Image registration method, visual positioning method, related device and equipment
CN114219889A (en) * 2021-11-04 2022-03-22 华南理工大学 Spine in-vivo quasi-static three-dimensional measurement method based on CT and ultrasonic fusion
CN116580033A (en) * 2023-07-14 2023-08-11 卡本(深圳)医疗器械有限公司 Multi-mode medical image registration method based on image block similarity matching
CN116580033B (en) * 2023-07-14 2023-10-31 卡本(深圳)医疗器械有限公司 Multi-mode medical image registration method based on image block similarity matching

Also Published As

Publication number Publication date
CN105046644B (en) 2021-08-13

Similar Documents

Publication Publication Date Title
CN105046644B (en) Ultrasonic and CT image registration method and system based on linear correlation
Huang et al. Fully automatic three-dimensional ultrasound imaging based on conventional B-scan
US10660613B2 (en) Measurement point determination in medical diagnostic imaging
CN106456125B (en) System for linking features in medical images to anatomical model and method of operation thereof
US8831708B2 (en) Multi-modal medical imaging
US9554772B2 (en) Non-invasive imager for medical applications
EP2701607B1 (en) Bone surface image reconstruction using ultrasound
JP5335280B2 (en) Alignment processing apparatus, alignment method, program, and storage medium
CN112384146B (en) Identifying optimal images from several ultrasound images
JP7253560B2 (en) Guided Transcranial Ultrasound Imaging Using Neural Networks and Related Apparatuses, Systems and Methods
EP2807978A1 (en) Method and system for 3D acquisition of ultrasound images
Zhang et al. Synthetic tracked aperture ultrasound imaging: design, simulation, and experimental evaluation
WO2010064348A1 (en) Information processing apparatus, information processing method, and program for positioning medical image
EP3463097B1 (en) Correcting probe induced deformation in an ultrasound fusing imaging system
Jiang et al. Deformation-aware robotic 3D ultrasound
WO2009136461A1 (en) Ultrasonograph
US10980509B2 (en) Deformable registration of preoperative volumes and intraoperative ultrasound images from a tracked transducer
Herickhoff et al. Low-cost volumetric ultrasound by augmentation of 2D systems: Design and prototype
JP7075854B2 (en) Ultrasonic diagnostic equipment and display method
US20200051257A1 (en) Scan alignment based on patient-based surface in medical diagnostic ultrasound imaging
JP2008183063A (en) Medical image diagnostic apparatus, medical image display device and program
BR112020014733A2 (en) COMPUTER IMPLEMENTED METHOD FOR OBTAINING ANATOMICAL MEASUREMENTS ON AN ULTRASOUND IMAGE, COMPUTER PROGRAM MEDIA, IMAGE ANALYSIS DEVICE AND ULTRASOUND IMAGING METHOD
Bø et al. Versatile robotic probe calibration for position tracking in ultrasound imaging
JP2022545219A (en) Ultrasonic guidance dynamic mode switching
JP2012071138A (en) Slice image display ultrasonic diagnostic apparatus of target object and method thereof

Legal Events

Date Code Title Description
C06 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