TWI676024B - 用以組合彩色3d影像之系統及方法 - Google Patents
用以組合彩色3d影像之系統及方法 Download PDFInfo
- Publication number
- TWI676024B TWI676024B TW107103473A TW107103473A TWI676024B TW I676024 B TWI676024 B TW I676024B TW 107103473 A TW107103473 A TW 107103473A TW 107103473 A TW107103473 A TW 107103473A TW I676024 B TWI676024 B TW I676024B
- Authority
- TW
- Taiwan
- Prior art keywords
- array
- density array
- processed
- density
- original
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000003491 array Methods 0.000 claims abstract description 76
- 239000013598 vector Substances 0.000 claims abstract description 63
- 230000009471 action Effects 0.000 claims abstract description 24
- 230000005855 radiation Effects 0.000 claims abstract description 15
- 230000000149 penetrating effect Effects 0.000 claims abstract description 8
- 239000002872 contrast media Substances 0.000 claims description 25
- 238000005286 illumination Methods 0.000 claims description 11
- 230000035515 penetration Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 description 12
- 210000000056 organ Anatomy 0.000 description 9
- 238000002591 computed tomography Methods 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 210000001835 viscera Anatomy 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 208000032005 Spinocerebellar ataxia with axonal neuropathy type 2 Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 208000033361 autosomal recessive with axonal neuropathy 2 spinocerebellar ataxia Diseases 0.000 description 1
- 230000005266 beta plus decay Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000496 pancreas Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/50—Lighting effects
- G06T15/506—Illumination models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/60—Rotation of whole images or parts thereof
-
- 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/10024—Color image
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
-
- 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/30096—Tumor; Lesion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本文揭露一種用以將藉由以穿透性輻射來執行的三維掃描所獲得的資料視覺化之系統及方法。原始密度陣列從該掃描形成,各原始密度陣列為一組三維陣列。一組已處理密度陣列由一或多個操作來形成,諸如取兩個原始密度陣列之間的差值、旋轉該已處理密度陣列、將該已處理密度陣列乘上一組前照明陣列、以及將該已處理密度陣列投影到一平面上以形成一影像,而投影動作包括為各對應於該影像的一像素之一組向量中的各向量計算多個統計數據中的一或更多個,而該等多個統計數據包括:一向量平均值、一向量最大值、以及一向量標準差。
Description
根據本發明的實施例的一個或更多層面係有關於成像技術,以及更特別地係關於顯示出表示一對象的內在結構之影像的技術。
電腦化軸向斷層(CAT)掃描為可使用以獲得有關諸如患者之對象的內在結構之資訊的一診斷程序。CAT掃描可使用以檢測癌症的指標。但在諸如胰臟的一些器官中的癌症,因為該器官在身體中的位置以及周圍組織的同質性,而可能難以用CAT掃描來檢測。這些特性可能減少CAT掃描影像中可感知的細節。
因此,需要有用以處理來自CAT掃描的資料之改善的系統及方法。
本揭露內容之實施例之數個層面針對一種用以把藉由一或多個例如以穿透性輻射來執行的三維掃描所獲得的資料視覺化之系統及方法。原始密度陣列從該掃描形成,各原始密度陣列為一組三維陣列。一組已處理密度陣列由一或多個操作來形成,諸如取兩個原始密度陣列之間的差值、旋轉該已處理密度陣列、將該已處理密度陣列乘上一組前照明陣列、以及將該已處理密度陣列投影至一平面上以形成一影像,而投影動作包括為各對應於該影像的一像素之一組向量中的各向量計算多個統計數據中的一或更多個,而該等多個統計數據包括:一向量平均值、一向量最大值、以及一向量標準差。
根據本發明之一實施例提供了一種用以視覺化掃描資料之方法,該方法包括:從多個原始密度陣列形成一第一已處理密度陣列,該第一已處理密度陣列及該等原始密度陣列各為具有多個陣列元素的一組三維陣列,各原始密度陣列之各元素表示一對象之掃描部分的密度;將該第一已處理密度陣列投影至一平面上以形成多個像素位置處的多個第一像素值,該投影包括:為該第一已處理密度陣列之多個陣列元素的各像素,沿著垂直於該平面及通過該像素位置的一直線,形成一向量,該向量對應到多個陣列元素;為各向量計算多個統計數據;及作為該等多個統計數據之統計數據加權總合來為各向量計算該第一像素值;以及顯示包括多個像素的一影像,而各像素具有等於一對應第一像素值之第一值。
在一實施例中,該等多個統計數據包括從由一向量平均值、一向量最大值及一向量標準差構成的群組中選擇的兩個統計數據,其中該方法進一步包括:從該等多個原始密度陣列形成一第二已處理密度陣列;投影該第二已處理密度陣列至該平面上以形成該等多個像素位置處之多個第二像素值;從該等多個原始密度陣列形成一第三已處理密度陣列;及投影該第三已處理密度陣列至該平面上以形成該等多個像素位置處之多個第三像素值,其中顯示影像包括以與該對應第一像素值成比例的第一顏色強度、與對應第二像素值成比例的第二顏色強度、及與對應第三像素值成比例的第三顏色強度來顯示該影像的每個像素。
在一實施例中,從該等多個原始密度陣列形成該第一已處理密度陣列,包括:從該等多個原始密度陣列中的一第二原始密度陣列減掉該等多個原始密度陣列中的一第一原始密度陣列;從該等多個原始密度陣列形成該第二已處理密度陣列,包括:從該等多個原始密度陣列中的一第三原始密度陣列減掉該第一原始密度陣列;及從該等多個原始密度陣列形成該第三已處理密度陣列,包括:從該等多個原始密度陣列中的一第四原始密度陣列減掉該第一原始密度陣列。
在一實施例中,該第一原始密度陣列為將對比劑注入該對象前所執行的掃描之結果,該第二原始密度陣列為將對比劑注入該對象一段第一時間區間後所執行的掃描之結果,該第三原始密度陣列為將對比劑注入該對象一段第二時間區間後所執行的掃描之結果,該第四原始密度陣列為將對比劑注入該對象一段第三時間區間後所執行的掃描之結果,及該第二時間區間比該第一時間區間長,且該第三時間區間比該第二時間區間長。
在一實施例中,從該等多個原始密度陣列形成該第一已處理密度陣列包括:從該等多個原始密度陣列中的一第二原始密度陣列減掉該等多個原始密度陣列中的一第一原始密度陣列,以形成一第三已處理密度陣列;及繞著與該第一原始密度陣列的三個維度之一平行的軸來將該第三已處理密度陣列旋轉一第一角度,以形成一第四已處理密度陣列。
在一實施例中,該軸與該平面平行。
在一實施例中,該第三已處理密度陣列含有多個二維切片,各切片為垂直於該軸的一組二維陣列,以及其中旋轉該第三已處理密度陣列的動作包括:使該等切片中之各切片旋轉該第一角度,以形成該第四已處理密度陣列之對應切片。
在一實施例中,旋轉該第三已處理密度陣列之切片的動作包括:從該第三已處理密度陣列之該切片的多個元素內插該第四已處理密度陣列之一對應切片的一元素的值。
在一實施例中,該內插的動作包括使用雙線性內插法。
在一實施例中,各向量為該第一已處理密度陣列的一維子陣列,該平面垂直於該陣列之三個維度的第一維度,且該一維子陣列係沿著該第一維度。
在一實施例中,形成該第一已處理密度陣列的動作包括:藉由從一第二原始密度陣列減掉一第一原始密度陣列來形成一第三已處理密度陣列,該第三已處理密度陣列含有多個二維切片,各切片為垂直於該平面的二維陣列,以及逐個元素地將各切片乘上與該切片具有相同維度的一前照明陣列;該前照明陣列中較靠近沿著垂直於該平面的一直線而與該第三已處理密度陣列分離的一虛擬觀察點之一第一陣列元素,具有比該前照明陣列中較遠離該虛擬觀察點之一第二陣列元素的值更小的值。
在一實施例中,在該前照明陣列之一鄰近元素子集內,該前照明陣列之各元素具有與該元素之一第一坐標的平方成比例的一值,該第一坐標為在垂直於該平面的一方向上的坐標。
在一實施例中,在該前照明陣列之一第一鄰近元素子集內,該前照明陣列之各元素具有與該元素之一第一坐標的四次方成比例的一值,該第一坐標為在垂直於該平面的一方向上的坐標。
在一實施例中,本方法包括該前照明陣列之該第一鄰近元素子集,有超出一臨界值的一第一坐標之該前照明陣列的各元素具有為零之一值。
根據本發明之一實施例提供了一種用以產生一對象之一內部視圖的系統,該系統包括:一掃描器,用以利用穿透性輻射來掃描該對象,及測量穿過該對象的該穿透性輻射之穿透量;一處理電路;及一顯示器,該處理電路組配來:從多個原始密度陣列形成一第一已處理密度陣列,該第一已處理密度陣列和該等原始密度陣列各為包括多個陣列元素之三維陣列,各原始密度陣列的各元素表示掃描對象的密度;投影該第一已處理密度陣列至一平面上以形成多個像素位置處之多個第一像素值,該投影動作包括:為該第一已處理密度陣列之該等多個陣列元素的各像素,沿著垂直於該平面並穿過該像素位置的一直線,形成一向量,該向量對應於多個陣列元素;為各向量計算多個統計數據;及作為該等多個統計數據之統計數據加權總合來為各向量計算該第一像素值;以及顯示包括多個像素的一影像,而各像素具有等於一對應第一像素值之第一值。
在一實施例中,該等多個統計數據包括從由一向量平均值、一向量最大值及一向量標準差構成的群組中選擇的兩個統計數據,及其中從該等多個原始密度陣列形成該第一已處理密度陣列的動作包括:從該等多個原始密度陣列中的一第二原始密度陣列減掉該等多個原始密度陣列中的一第一原始密度陣列,以形成一第三已處理密度陣列;及繞著與該第一原始密度陣列的三個維度之一平行的軸來將該第三已處理密度陣列旋轉一第一角度,以形成一第四已處理密度陣列。
在一實施例中,該第三已處理密度陣列含有多個二維切片,各切片為垂直於該軸的一組二維陣列,以及其中旋轉該第三已處理密度陣列的動作包括:將該等切片中之各切片旋轉該第一角度,以形成該第四已處理密度陣列之對應切片。
在一實施例中,形成該第一已處理密度陣列的動作包括:藉由從一第二原始密度陣列減掉一第一原始密度陣列來形成一第三已處理密度陣列,該第三已處理密度陣列含有多個二維切片,各切片為垂直於該平面的二維陣列,以及藉由逐個元素地將各切片乘上具有與該切片相同維度的一前照明陣列;該前照明陣列中較靠近沿著垂直於該平面的一直線而與該第三已處理密度陣列分離的一虛擬觀察點之一第一陣列元素,具有比該前照明陣列中較遠離該虛擬觀察點之一第二陣列元素的值更小的值。
在一實施例中,在該前照明陣列之一鄰近元素子集內,該前照明陣列之各元素具有與該元素之一第一坐標的平方成比例的一值,該第一坐標為在垂直於該平面的一方向上的坐標。
根據本發明之一實施例提供了一種用以產生一對象之一內部視圖的系統,該系統包括:一處理電路、及一顯示器;該處理電路組配來:從多個原始密度陣列形成一第一已處理密度陣列,該第一已處理密度陣列和該等原始密度陣列各為包括多個陣列元素之三維陣列,各原始密度陣列的各元素表示掃描對象的密度;投影該第一已處理密度陣列至一平面上以形成多個像素位置處之多個第一像素值,該投影動作包括:為該第一已處理密度陣列之該等多個陣列元素的各像素,沿著垂直於該平面並穿過該像素位置的一直線,形成一向量,該向量對應於多個陣列元素;為各向量計算多個統計數據;及作為該等多個統計數據之統計數據加權總合為各向量計算該第一像素值;以及顯示包括多個像素的一影像,而各像素具有等於一對應第一像素值之第一值。
以下結合附加圖式闡述的詳細描述旨在作為對根據本發明提供的用以顯示三維成像資料的系統和方法的示例性實施例的描述,並且不欲表示可以構建或利用本發明的唯一形式。該描述結合繪示的實施例來闡述本發明的特徵。然而可以理解的是,相同或等效的功能與結構可由不同的實施例達成,該等實施例也欲包括於本發明的精神與範圍之中。如本文它處所示,相似的元件標號係欲標示相似的元件或特徵。
電腦化軸向斷層(CAT)掃描為一種程序,其中一對象(例如患者)從數個方向以來自輻射源的穿透性輻射(例如X射線)來照射,以及透射的輻射的一原始掃描影像每次由檢測器形成,以形成多個原始掃描影像,其各可表示為二維陣列。該輻射在不同種類的物質中可能會以不同速率衰減;因此,取決於該輻射從輻射源行進到檢測器所沿循的路徑上之物質的成分之衰減率,各影像中之各點可對應於一透射輻射強度。從一些原始掃描影像的組合即可形成該對象的“密度”的一個三維模型,其中當本文針對CAT掃描而用到時,該“密度”一詞是指穿透性輻射的局部衰減率。雖然一些範例在本揭露內容中係以人類患者的CAT掃描之情境來討論,但本發明不限於此,且在一些實施例中,諸如磁振造影掃描或正子發射斷層攝影掃描、或對人類患者之外的對象之掃描來提供三維密度資料的其它種類掃描,可以類似的方式來處理。在其它實施例中,例如使用另一個掃描方法,其產生在空間上顯示在三維至次像素層級的掃描;在其他實施例中,例如,使用產生空間上以三維對準到次體素層級的掃描之另一種掃描方法,並且其中在掃描中的值之間存在差異,而使其適用於三色顯示器。在其他類型掃描的情況下,密度可以被相對應地定義;例如在正子發射斷層攝影掃描的情況下,該密度可以為由貝他正放射(beta plus emission)衰減的原子核之密度。當於本文使用時,“對象”一詞包括可被掃描的任何事物,且包含但不限於人類患者、動物、植物、無生命物、及其組合。
當成像的對象為人類患者(或其它活體對象)時,可使用對比劑(例如注射進入或由患者攝入)來選擇性地改變一些組織的密度。該對比劑例如可包括一相對地不透光的物質(即對穿透性輻射相對地不易透射的)。因此含有該對比劑的該組織的密度可以提高,並且它可以增加到依該組織中對比劑的濃度而定的程度。圖1根據一實施例顯示用以進行一掃描以及處理及顯示結果的一系統之方塊圖。該系統包括一掃描器110、一處理電路115 (後文詳述)、用以採影片形式顯示影像或連續影像之一顯示器120、以及諸如鍵盤或滑鼠的一或多個輸入裝置125,此等輸入裝置可由一操作者(例如放射科人員)用來操作該系統,以及設定影響要顯示的該等影像的處理之參數。應該注意的是,該處理電路115、顯示器120、以及輸入裝置125可以為單一系統的一部分,或可以為分散式系統而讓處理電路115例如分離及通訊式耦合至該顯示器120以及輸入裝置125。在一些實施例中,伺服器儲存該等影像且用戶端調用該等影像,而影像處理在伺服器或用戶端上執行,或者二者皆有。
多個掃描可以執行及同時分析。例如第一次掃描可在注射對比劑之前執行,以及接著幾次掃描可在注射該對比劑之後隨著對比劑的濃度變化於不同時間(例如以固定的區間)執行。對比劑的濃度最初增加的速率、達到的尖峰濃度、以及對比劑的濃度隨後減少的速率,全部可取決於組織的類型。
在一些實施例中,可以利用多種方法自CAT掃描資料產生影像,來幫助CAT掃描使用作為診斷工具。繪示在圖2並進一步揭露細節於下文的一系列的步驟或“動作”可被使用來例如自一組四個掃描形成一彩色影片。在一些實施例中,則使用多或少於四個掃描。在一些實施例中,含有三維掃描資料之陣列可以於動作205中(從該掃描器110,經由轉換原始掃描影像成為原始密度陣列的處理電路115的一部分,或經由在處理電路115中執行來實現此轉換的程序)被接收,在動作210中相減形成差值,且於動作215中過濾(在選擇性加入一偏值之後)及剪輯。結果接著可以在動作220中旋轉,在動作225中遭受前照光程序,及在動作230中投影在一平面上。在動作235中可以組合幾個這種投影(例如三個投影)以形成一彩色影像(此三個投影例如對應於該彩色影像的紅色、綠色和藍色成分)。一系列的彩色影像(各別已旋轉例如至稍微不同的角度)可以藉由重複動作220-230來形成,然後可以在動作240中合併成一影片(例如顯示掃描對象旋轉的影片)。在動作245中可以對操作者顯示該影片。動作205-240可以由處理電路115執行。
列在列表1上的程式碼顯示用以產生可以作為影片來觀看的一系列影像之MATLABTM
程式碼。列表1的該程式碼接收多個CAT掃描資料檔案作為輸入。各檔案:sample_patient_1_a3_no_filt.mat、sample_patient_2_a3_no_filt.mat、 sample_patient_3_a3_no_filt.mat、以及sample_patient_4_a3_no_filt.mat 在變數a3中定義這裡稱為“原始密度陣列”的一個三維陣列(例如第一原始密度陣列、第二原始密度陣列、第三原始密度陣列、以及第四原始密度陣列)。各原始密度陣列的各元素表示患者體內與該元素的坐標相對應之實體位置處的小體積或“體素”組織的密度。
這些(三維)原始密度陣列可以從(二維)原始掃描影像,使用例如反摺積程序以推斷被掃描對象(例如患者)的三維密度來產生(然後在動作205中被接收)。第一個檔案可以包括從對比劑注入患者之前進行之掃描得到的資料,以及剩下三個檔案可以包括從對比劑注入患者之後各間隔兩分鐘來進行之多次掃描得到的資料。然而,應當理解的是,使用對比劑並非是必須的,及/或能夠使用其它的時間區間或影像的數量。例如單色投影、以及單色影片可以從單一掃描,使用下文進一步詳述的方法來形成。列表1的第34-41行具有讀取該四個檔案及分別複製個別3維陣列至四個個別變變scan1、scan2、scan3及scan4的效果。在列表1的第43-45行中,第一個掃描(無攝取對比劑)從其它三個掃描之各者減掉,以分別形成三組密度陣列、或“已處理”密度陣列於變數scan21、scan31、及scan41。當于本文中使用時,“已處理密度陣列”為從一或多個原始密度陣列直接或間接形成的任何三維密度陣列,以及係用來形成影像上的投影或藉由進一步的處理來形成另一個已處理密度陣列者。由(透過其中有某些部分可以任選的一或多個處理步驟)處理原始密度陣列或其它已處理密度陣列而產生的該等已處理密度陣列,可以稱作第一已處理密度陣列、第二已處理密度陣列、第三已處理密度陣列、第四已處理密度陣列、第五已處理密度陣列、第六已處理密度陣列等等,以彼此區分之。
如上所述,各密度陣列的每個體素可以對應到該對象(例如該患者)中的一實體位置。該位置的定位可以由有X、Y和Z軸的三維(例如笛卡兒)坐標系統的坐標來識別。因此,每個體素可以具有X、Y和Z坐標,其可以由各該體素相對於坐標系原點(可能在患者體內某處)的位置(例如以厘米或毫米為單位)來指定。該密度陣列的各元素可以由三個整數的元組(或三元組)索引,且各元素也可以由其對應的體素的坐標(為了簡潔起見,其可以被稱為陣列元素的坐標)來識別。此三元組可以寫成(x, y, z),並且可以包括範圍在1和Nx之間的x指數、範圍從1到Ny的y指數、以及範圍從1到Nz的z指數,其中Nx、Ny、Nz是列表1第64行定義的整數,其指明原始密度陣列的大小。
在列表1的第50-52行中,三個差值scan21、scan31、scan41各被一個三維中值濾波器予以過濾。濾波器輸出的各體素可以為濾波器輸入中對應體素周圍的一個m x m x m 鄰域的中位數,而變數“NMED”含有純量m。應注意的是,列表1內的函數medfilt3()並非建立在MATLABTM
中相同名稱的函數,而是具有相似效果的另一個函數(但接收NMED的純量值而非向量)。三個濾波器輸出被分配給個別的顏色(其中變數red對應到紅色、變數grn對應到綠色、變數blu對應到藍色),以便在列表1的後續部分中顯示。
在第55-57行中,值小於0的已分配顏色的濾波器輸出之元素被進位(或“截縮”)為零,以及在第60-62行中,三組新陣列被初始化(成單精度數字,顧及處理速度並避免使用過多的記憶體)。一偏值(變數偏值)可以用來調整在過濾及顏色分配步驟期間的對比度;藉由在截縮操作前加入一負偏值,例如可以將較低強度值設置為零,來增加對比度。列表1中的偏值設置為-1050;在一些實施例中,它由操作者例如使用輸入裝置125來設定,或者它可以基於所使用的對比劑或正在成像的身體部位來設定。
例如,可以將人類患者的縱軸(即高度軸)定義為Z軸、將前後軸定義為X軸、將左右軸定義為Y軸。然後,影像序列可以作為以不同量來旋轉的三維陣列在一平面上之一系列投影來形成。各投影可以近似於位於沿著垂直於該平面的直線(“模擬觀察方向”)而與患者相距一距離的“虛擬觀察點”處的“虛擬觀察者”,在恰似患者部分透明時所能見的圖像,從而患者的內部器官即可見。例如,於一旋轉平台上,表示組織密度的該陣列可以圍繞Z軸旋轉並投影到Y-Z平面上,以產生對應於該虛擬觀察者在恰似該患者部分透明、以及恰似該患者站在該虛擬觀察者前方時所能見的圖像之一系列影像,而該虛擬觀察者乃從沿著X軸的一模擬觀察方向來觀察該患者。在其它實施例中,三維陣列的旋轉軸可以為X軸或Y軸而非Z軸,以及模擬方向可以沿著Y軸或Z軸。
在列表1中,常數ROTXY、ROTXZ以及ROTYZ被使用來識別(i)表示組織密度的陣列被旋轉時所繞的軸線、以及(ii)該虛擬觀察者的觀察方向的不同組合。在一些實施例中,該模擬觀察方向被選擇為垂直於旋轉軸(即已處理密度矩陣所被投影上去的平面被選擇為平行於旋轉軸),從而於密度陣列旋轉時,一些器官在其他器官的前面或後面通過。為了便於計算,該平面可以被選擇為垂直於該系統的其中一軸(以及垂直於該已處理密度陣列的其中一維)。于此狀況下,投影可以一次對該已處理密度陣列的一向量進行操作,各向量是垂直於該平面的一維子陣列。在其它實施例中,該平面處於任意角度,使得垂直於該平面並穿過一像素位置的每條直線可能未與已處理密度陣列的三個維度中的任何維度對齊;在這種情況下,對直線所穿過的像素執行投影之前,可以使用內插來形成沿著每條這樣的直線之密度值的向量。
在各投影中,加權可用來增強例如某些器官之某些形貌體或例如腫瘤的異常部的可見性。例如靠近該虛擬觀察者的形貌體可以被增強,以及較遠的形貌體可以被抑制,使得在前視圖中,例如可能具有高密度但離虛擬觀察者較遠的脊椎,不一定會使較靠近虛擬觀察者的器官不明晰。
列表1的第70-72行對應到表示組織密度的該等陣列繞著Z軸旋轉及投影到Y-Z平面上的情況。當旋轉軸為Z軸,以及模擬觀察方向為沿著X軸時,變數front_lighting中的加權陣列可以用來強調較靠近虛擬觀察者的對象的部分。例如第71及72行於X-Y平面中形成一組二維陣列(即具有與X方向相關聯的第一坐標以及與Y方向相關聯的第二坐標),其隨著X指數之數值範圍內的X指數之四次方而變化(該範圍不包括X方向上的前120個元素和最後50個元素),以及在其它地方為零(即X方向上的前120個元素和最後50個元素為零)。以這種方式,若例如在一前視圖中,該脊椎完全落在距離該虛擬觀察者最遠的50個元素內,則它在旋轉的三維陣列乘以陣列front_lighting之後,可以完全從投影中消除。類似地,患者身上可能也具有高密度以及相當高的遮蔽內部器官的能力之肋骨的一部分,可能位於最靠近離虛擬觀察者的120個元素內,並且也可以透過使用第71和72行的前照光陣列來消除。其它相似的前照光陣列被形成來。其它類似的前照光陣列(與二次方相關地變化,而非隨四次方而變)被形成用於沿著Z方向來觀看、用於圍繞Y軸旋轉(在第67行中)或繞X軸旋轉(在第69行中)。在列表1中,對於ROTXZ和ROTYZ的情況使用二次方加權,以及對於ROTXY的情況使用四次方加權;在其它實施例中,該二者中任一類型的加權(或其它種類的加權,例如立方加權)可以搭配旋轉軸以及模擬觀察方向的任何組合使用。
第87-169行定義了一組迴圈(在本文中稱為“角度迴圈”),其用於每運作一次產生該系列影像中的一影像,每一影像針對已處理密度陣列的一不同旋轉角度(以及模擬穿過患者的不同觀察方向)。角度的範圍和角度增量定義在第75-80行中。
在各角度,三組已處理密度陣列(對應於紅色、綠色以及藍色)各在由第89行陸續到123行的“for”迴圈(在本文中稱為“切片迴圈”)中,一次一個切片地處理,由切片指數n來索引。各切片是在對應於沿著旋轉軸的一位置之該指數n的一個值時垂直於旋轉軸的二維陣列。例如就繞著Y軸旋轉(第91-94行)而言,平行於XZ平面的切片係從各個已處理密度陣列red_temp、grn_temp以及blu_temp中提取。
然後,各切片在第106-108行上對imrotate()的對應呼叫中,被旋轉當前角度。MATLABTM
函式imrotate()將一組二維陣列旋轉,該陣列為該函式的第一個引數,其中係旋轉了在第二個引數中指定的角度,使用在第三個引數中指定的內插方法以及具有由第四個引數所指定的大小。在第106-108行的範例程式碼中,指定了雙線性內插,而值'crop'指定由該函式送回的陣列與作為第一個引數傳遞給該函式的陣列為相同大小。
在第110-122行中,從該等切片組合出三組新的三維陣列,並且應用前照明以強調較靠近虛擬觀察者之此等三維陣列的部分。前照明是藉由將每個旋轉的切片乘以陣列front_lighting來達成的,其(如上文所述)可能含有在虛擬觀察者之方向上較大的值。如上文所述,這陣列也可能在一些元素裡包括值0,其可以具有將切片的區域(例如位於離觀察者最遠的一條體素內者及/或位於離虛擬觀察者最近的一條體素內者)設置為零的效果。
三組“單色”陣列red_im、grn_im及blu_im在第128-130行(針對平行於Z軸的觀察方向,即ROTXZ或ROTYZ),或在第150-152行(針對平行於X軸的觀察方向,即ROTXZ)形成。各單色陣列的各像素被設定為等於平行於觀察方向及穿過像素位置的體素向量之多達三個統計數據的加權和,這三個統計數據是向量平均值、向量最大值和向量標準差。例如,參考第128行,該陣列red_im之坐標(27, 52)處的一像素(為了例示說明之目的,坐標27和52係作為任意值在這裡被拿來使用),將會作為坐標為(27,53,n)的體素向量的三個統計數據(即向量平均值、向量最大值和向量標準差)之加權和來形成,其中指數n的範圍從1到Nz。此體素向量平行於觀察方向,以及垂直於三維陣列red1為了形成該單色陣列red_im所投影到的平面。MATLABTM
函式mean()沿著其第二個引數指出的一個維度計算其第一個引數指出之陣列的平均值,而送回比第一個引數指出的陣列少一維的陣列。對於三維陣列,例如,沿著第三維度(例如Z維度)的平均值是一個二維陣列,其為所含每個元素是三維陣列內平行於Z軸的對應向量之平均值的二維陣列。類似地,MATLABTM
函式max()和std()分別送回輸入陣列內的向量之最大值的陣列,或輸入陣列內的對應向量的標準差(例如樣本標準差)的陣列。權重c1、c2和c3可以被選擇(例如由系統設計者選擇、由系統自動選擇、或由使用者在操作期間透過合適的使用者介面來選擇)來增強感興趣的對象之部位的可見性,及抑制非感興趣的部位。例如特定器官可以比其他組織更快速地吸收對比劑,而且因此,它在注射了對比劑之後的第一次掃描中,可以具有比其它組織明顯更高的密度。這可以導致對於通過該器官的向量之標準差的數值會很大,而且如果在投影中標準差的加權很大(即若c3被設定為相當高時),則在每個投影中,此器官即能以對應於注入對比劑之後的第一次掃描之顏色呈現為明亮區域。
第132-140行和第144-147行是(當列表1的程式碼進行其它更改時)用於診斷的多行程式碼,且不會影響列表1中程式碼的執行。
在第157-159行中,形成包含多個像素的一彩色影像,每個像素具有第一像素值(例如紅色強度)、第二像素值(例如綠色強度)和第三像素值(例如藍色強度)。該彩色影像被組合在三維陣列col_im中,該三維陣列的前兩個維度對應於像素坐標,且其最後一個維度(其值取1、2或3)對應於三種顏色(例如分別為紅色、綠色和藍色)。col_im的三個顏色成分是三個“單色”陣列red_im、grn_im及blu_im。
列表1的第161-167行將當前彩色影像作為一圖框添加到影片中,在播放時,它將顯示由角度迴圈所生成的一系列影像,結果影片將顯示圍繞旋轉軸而旋轉的掃描對象。名為增益的變數會調整影像(和影片)的亮度。這個變數在第11行中設定;在其它實施例中,它由操作員例如使用輸入裝置125來設定。 列表1
本文中使用的“處理電路”一詞包括:用來處理資料或數位信號的硬體、韌體以及軟體之任何組合。處理電路硬體可以包括:例如特定應用積體電路(ASIC)、通用或特殊用途中央處理單元(CPU)、數位信號處理器(DSP)、圖形處理單元(GPU)、以及諸如可規劃邏輯閘陣列(FPGA)的可規劃邏輯裝置。如本文所用的一處理電路中,各功能由組配例如硬佈線來執行該功能的硬體來執行、或由組配來執行儲存在非暫態儲存媒體中的指令之諸如CPU的更通用硬體來執行。處理電路可以製作在單一印刷線路板(PWB)上,或分散在數個互連的PWB上。一處理電路可以含有其它處理電路;例如一個處理電路可以包括在一個PWB上互連的一個FPGA及一個CPU的兩個處理電路。一處理電路可以包括多個處理單元,該等處理單元為地理上相隔兩處、且例如由諸如網際網路的網路來連接。
將可知的是,儘管“第一”、 “第二”、 “第三”等等語詞,在本文中可能用來描述各種元件、構件、區域、陣列及/或部段,但這些元件、構件、區域、層體及/或部段不應受限於這些用語。這些用語僅使用來區分某一元件、構件、區域、層體或部段與另一元件、構件、區域、陣列或部段。因此,例如下文討論的第一元件、構件、區域、陣列或部段,在不脫離本發明概念的精神和範圍的情況下,可能改稱為第二元件、構件、區域、陣列或部段。
儘管本文已具體地描述及繪示出用以顯示三維成像資料的系統及方法之有限的實施例,但許多修改及變化將可被熟於此技者所明顯看出。因此,應該理解的是,根據本發明的原理採用來用以顯示三維成像資料的系統和方法,可以不同於本文中所具體描述的方式來體現。本發明也在後附申請專利範圍及其等效者中來定義。
110‧‧‧掃描器
115‧‧‧處理電路
120‧‧‧顯示器
125‧‧‧輸入裝置
205、210、215、220、225、230、235、240、245‧‧‧動作
下文結合附圖來描述多種特徵、層面、及實施例,其中:
圖1為一種根據本發明的一實施例用以產生一對象之內部影像的系統;及
圖2為根據本發明的一實施例用以產生一對象之一內部影像之方法的流程圖。
Claims (20)
- 一種用以將掃描資料視覺化的方法,該方法包含:從由一或多個三維掃描所形成的多個原始密度陣列形成一第一已處理密度陣列,該第一已處理密度陣列及該等原始密度陣列各為具有多個陣列元素的一組三維陣列,各該原始密度陣列之各元素表示一對象之一掃描部分的密度;將該第一已處理密度陣列投影到一平面上以形成多個像素位置處的多個第一像素值,該投影包含:為該第一已處理密度陣列之該等多個陣列元素的各像素,沿著垂直於該平面及通過該像素位置的一條線,形成一向量,該向量對應到多個陣列元素;為各向量計算多個統計數據;及作為該等多個統計數據之統計數據加權總合來為各向量計算該第一像素值;以及顯示包含多個像素的一影像,而各像素具有等於一對應第一像素值之一第一值。
- 如請求項1之方法,其中該等多個統計數據包含從由一向量平均值、一向量最大值及一向量標準差構成的群組中選擇的兩個統計數據,其中該方法進一步包含:從該等多個原始密度陣列形成一第二已處理密度陣列;投影該第二已處理密度陣列到該平面上以形成該等多個像素位置處之多個第二像素值;從該等多個原始密度陣列形成一第三已處理密度陣列;及投影該第三已處理密度陣列到該平面上以形成該等多個像素位置處之多個第三像素值,其中顯示該影像包含:以與該對應第一像素值成比例的一第一顏色強度、與對應第二像素值成比例的一第二顏色強度、及與對應第三像素值成比例的一第三顏色強度,來顯示該影像的每個像素。
- 如請求項2之方法,其中:從該等多個原始密度陣列形成該第一已處理密度陣列包含:從該等多個原始密度陣列中的一第二原始密度陣列減掉該等多個原始密度陣列中的一第一原始密度陣列;從該等多個原始密度陣列形成該第二已處理密度陣列包含:從該等多個原始密度陣列中的一第三原始密度陣列減掉該第一原始密度陣列;及從該等多個原始密度陣列形成該第三已處理密度陣列包含:從該等多個原始密度陣列中的一第四原始密度陣列減掉該第一原始密度陣列。
- 如請求項3之方法,其中:該第一原始密度陣列為將對比劑注入該對象前所執行的掃描之結果,該第二原始密度陣列為將對比劑注入該對象一段第一時間區間後所執行的掃描之結果,該第三原始密度陣列為將對比劑注入該對象一段第二時間區間後所執行的掃描之結果,該第四原始密度陣列為將對比劑注入該對象一段第三時間區間後所執行的掃描之結果,及該第二時間區間比該第一時間區間長,且該第三時間區間比該第二時間區間長。
- 如請求項1之方法,其中從該等多個原始密度陣列形成該第一已處理密度陣列包含:從該等多個原始密度陣列中的一第二原始密度陣列減掉該等多個原始密度陣列中的一第一原始密度陣列,以形成一第三已處理密度陣列;及繞著與該第一原始密度陣列的三個維度之一平行的一軸線來將該第三已處理密度陣列旋轉一第一角度,以形成一第四已處理密度陣列。
- 如請求項5之方法,其中該軸線與該平面平行。
- 如請求項5之方法,其中該第三已處理密度陣列含有多個二維切片,各切片為垂直於該軸線的一組二維陣列,以及其中旋轉該第三已處理密度陣列的動作包含:使該等切片中之各切片旋轉該第一角度,以形成該第四已處理密度陣列之一對應切片。
- 如請求項7之方法,其中旋轉該第三已處理密度陣列之一切片的動作包含:從該第三已處理密度陣列之該切片的多個元素內插該第四已處理密度陣列之一對應切片的一元素的值。
- 如請求項8之方法,其中該內插的動作包含使用雙線性內插法。
- 如請求項1之方法,其中各向量為該第一已處理密度陣列的一組一維子陣列,該平面垂直於該陣列之三個維度的第一維度,且該一維子陣列係沿著該第一維度。
- 如請求項1之方法,其中形成該第一已處理密度陣列的動作包含:藉由從一第二原始密度陣列減掉一第一原始密度陣列來形成一第三已處理密度陣列,該第三已處理密度陣列含有多個二維切片,各切片為垂直於該平面的二維陣列,以及逐個元素地將各切片乘上與該切片具有相同維度的一前照明陣列;該前照明陣列中較靠近沿著垂直於該平面的一直線而與該第三已處理密度陣列分離的一虛擬觀察點之一第一陣列元素,具有比該前照明陣列中較遠離該虛擬觀察點之一第二陣列元素的值更小的值。
- 如請求項11之方法,其中,在該前照明陣列之一鄰近元素子集內,該前照明陣列之各元素具有與該元素之一第一坐標的平方成比例的一值,該第一坐標為在垂直於該平面的一方向上的坐標。
- 如請求項11之方法,其中,在該前照明陣列之一第一鄰近元素子集內,該前照明陣列之各元素具有與該元素之一第一坐標的四次方成比例的一值,該第一坐標為在垂直於該平面的一方向上的坐標。
- 如請求項13之方法,其中,在包含該前照明陣列之該第一鄰近元素子集的該前照明陣列之一第二鄰近元素子集內,有超出一臨界值的一第一坐標之該前照明陣列的各元素具有為零之一值。
- 一種用以產生一對象之一內部之視圖的系統,該系統包含:一掃描器,用以利用穿透性輻射來掃描該對象,及測量穿過該對象的該穿透性輻射之穿透量;一處理電路;及一顯示器,該處理電路組配來:從由一或多個三維掃描所形成的多個原始密度陣列形成一第一已處理密度陣列,該第一已處理密度陣列和該等原始密度陣列各為包含多個陣列元素之三維陣列,各該原始密度陣列的各元素表示掃描對象的密度;投影該第一已處理密度陣列到一平面上以形成多個像素位置處之多個第一像素值,該投影動作包含:為該第一已處理密度陣列之該等多個陣列元素的各像素,沿著垂直於該平面並通過該像素位置的一條線,形成一向量,該向量對應於多個陣列元素;為各向量計算多個統計數據;及作為該等多個統計數據之統計數據加權總合來為各向量計算該第一像素值;以及顯示包含多個像素的一影像,而各像素具有等於一對應第一像素值之一第一值。
- 如請求項15之系統,其中該等多個統計數據包含從由一向量平均值、一向量最大值及一向量標準差構成的群組中選擇的兩個統計數據,及其中從該等多個原始密度陣列形成該第一已處理密度陣列的動作包含:從該等多個原始密度陣列中的一第二原始密度陣列減掉該等多個原始密度陣列中的一第一原始密度陣列,以形成一第三已處理密度陣列;及繞著與該第一原始密度陣列的三個維度之一平行的一軸線來將該第三已處理密度陣列旋轉一第一角度,以形成一第四已處理密度陣列。
- 如請求項16之系統,其中該第三已處理密度陣列含有多個二維切片,各切片為垂直於該軸線的一組二維陣列,以及其中旋轉該第三已處理密度陣列的動作包含:將該等切片中之各切片旋轉該第一角度以形成該第四已處理密度陣列之一對應切片。
- 如請求項15之系統,其中形成該第一已處理密度陣列的動作包含:藉由從一第二原始密度陣列減掉一第一原始密度陣列來形成一第三已處理密度陣列,該第三已處理密度陣列含有多個二維切片,各切片為垂直於該平面的一組二維陣列,以及藉由逐個元素地將各切片乘上具有與該切片相同維度的一前照明陣列;該前照明陣列中較靠近沿著垂直於該平面的一直線而與該第三已處理密度陣列分離的一虛擬觀察點之一第一陣列元素,具有比該前照明陣列中較遠離該虛擬觀察點之一第二陣列元素的值更小的值。
- 如請求項18之系統,其中,在該前照明陣列之一鄰近元素子集內,該前照明陣列之各元素具有與該元素之一第一坐標的平方成比例的一值,該第一坐標為在垂直於該平面的一方向上的坐標。
- 一種用以產生一對象之一內部之視圖的系統,該系統包含:一處理電路;及一顯示器;該處理電路組配來:從由一或多個三維掃描所形成的多個原始密度陣列形成一第一已處理密度陣列,該第一已處理密度陣列和該等原始密度陣列各為包含多個陣列元素之三維陣列,各該原始密度陣列的各元素表示掃描對象的密度;投影該第一已處理密度陣列到一平面上以形成多個像素位置處之多個第一像素值,該投影動作包含:為該第一已處理密度陣列之該等多個陣列元素的各像素,沿著垂直於該平面並通過該像素位置的一條線,形成一向量,該向量對應於多個陣列元素;為各向量計算多個統計數據;及作為該等多個統計數據之統計數據加權總合為各向量計算該第一像素值;以及顯示包含多個像素的一影像,而各像素具有等於一對應第一像素值之一第一值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/489,652 US10102682B1 (en) | 2017-04-17 | 2017-04-17 | System and method for combining 3D images in color |
US15/489,652 | 2017-04-17 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201839389A TW201839389A (zh) | 2018-11-01 |
TWI676024B true TWI676024B (zh) | 2019-11-01 |
Family
ID=61244689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW107103473A TWI676024B (zh) | 2017-04-17 | 2018-01-31 | 用以組合彩色3d影像之系統及方法 |
Country Status (9)
Country | Link |
---|---|
US (3) | US10102682B1 (zh) |
EP (2) | EP3613019B1 (zh) |
JP (1) | JP6854919B2 (zh) |
KR (1) | KR102275334B1 (zh) |
CN (1) | CN110753954B (zh) |
CA (1) | CA3060256A1 (zh) |
IL (1) | IL270076B (zh) |
TW (1) | TWI676024B (zh) |
WO (1) | WO2018194739A1 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10102682B1 (en) * | 2017-04-17 | 2018-10-16 | Raytheon Company | System and method for combining 3D images in color |
EP3669325B1 (en) * | 2017-08-14 | 2023-09-27 | Raytheon Company | Subtraction algorithm for detection of tumors |
US10265138B2 (en) * | 2017-09-18 | 2019-04-23 | MediVis, Inc. | Methods and systems for generating and using 3D images in surgical settings |
CN109345580B (zh) * | 2018-10-23 | 2020-03-24 | 北京字节跳动网络技术有限公司 | 用于处理图像的方法和装置 |
US11062512B2 (en) | 2019-08-09 | 2021-07-13 | Raytheon Company | System and method for generating 3D color representation of 2D grayscale images |
US11475558B2 (en) | 2019-11-13 | 2022-10-18 | Raytheon Company | Organ isolation in scan data |
US11282209B2 (en) | 2020-01-10 | 2022-03-22 | Raytheon Company | System and method for generating contours |
US11562512B2 (en) | 2020-12-09 | 2023-01-24 | Raytheon Company | System and method for generating and displaying contours |
US11893745B2 (en) | 2020-12-09 | 2024-02-06 | Raytheon Company | System and method for generating and displaying contours |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020136439A1 (en) * | 2001-03-09 | 2002-09-26 | Ruchala Kenneth J. | System and method for fusion-aligned reprojection of incomplete data |
US20110135173A1 (en) * | 2009-10-11 | 2011-06-09 | Axs Ingenierie | Method for 3d spine and full 3d trunk imaging fusion in standing position |
TW201123076A (en) * | 2009-12-30 | 2011-07-01 | Pou Yuen Technology Co Ltd | Three-dimensional display method of medical images |
US20120038644A1 (en) * | 2009-02-17 | 2012-02-16 | Uws Ventures Ltd. | Rendering caustics in computer graphics with photon mapping by relaxation |
TWI509564B (zh) * | 2014-10-30 | 2015-11-21 | Inst Nuclear Energy Res Atomic Energy Council | 三維成像的投影方法 |
TWI517093B (zh) * | 2013-10-11 | 2016-01-11 | Univ Nat Yang Ming | Computer tomography reconstruction method |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4680709A (en) * | 1984-11-13 | 1987-07-14 | General Electric Company | Back projection image reconstruction apparatus and method |
JPH0454682A (ja) * | 1990-06-22 | 1992-02-21 | Toshiba Corp | 立体画像処理方法及びその装置 |
JPH0515525A (ja) * | 1991-07-15 | 1993-01-26 | Toshiba Corp | 画像位置修正方法 |
JP3185547B2 (ja) * | 1994-06-28 | 2001-07-11 | 三菱電機株式会社 | 距離測定装置 |
US6002739A (en) * | 1998-04-28 | 1999-12-14 | Hewlett Packard Company | Computed tomography with iterative reconstruction of thin cross-sectional planes |
US6519354B1 (en) * | 1998-09-17 | 2003-02-11 | Brigham And Women's Hospital, Inc. | Method and apparatus for projecting MR angiographic data |
JP2001118057A (ja) * | 1999-10-21 | 2001-04-27 | Hitachi Medical Corp | 透視像構成表示装置 |
AU2001251539A1 (en) * | 2000-04-11 | 2001-10-23 | Cornell Research Foundation Inc. | System and method for three-dimensional image rendering and analysis |
US6867774B1 (en) * | 2002-12-02 | 2005-03-15 | Ngrain (Canada) Corporation | Method and apparatus for transforming polygon data to voxel data for general purpose applications |
US8045770B2 (en) * | 2003-03-24 | 2011-10-25 | Cornell Research Foundation, Inc. | System and method for three-dimensional image rendering and analysis |
WO2005106773A2 (en) | 2004-04-15 | 2005-11-10 | Edda Technology, Inc. | Spatial-temporal lesion detection, segmentation, and diagnostic information extraction system and method |
JP2006170635A (ja) * | 2004-12-13 | 2006-06-29 | Hitachi Engineering & Services Co Ltd | 画像処理による欠陥の検出装置および検出方法 |
US7653264B2 (en) * | 2005-03-04 | 2010-01-26 | The Regents Of The University Of Michigan | Method of determining alignment of images in high dimensional feature space |
JP4327171B2 (ja) * | 2006-04-06 | 2009-09-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 最大値投影方法および装置 |
US20110150309A1 (en) * | 2009-11-27 | 2011-06-23 | University Health Network | Method and system for managing imaging data, and associated devices and compounds |
JP5980776B2 (ja) * | 2011-04-20 | 2016-08-31 | 株式会社日立製作所 | 医用画像表示装置、医用画像表示方法 |
AU2012312896B2 (en) * | 2011-09-20 | 2015-07-30 | Siemens Aktiengesellschaft | Bayesian approach for gas concentration reconstruction based on tunable diode laser absorption spectroscopy |
JP2013094229A (ja) * | 2011-10-28 | 2013-05-20 | Ge Medical Systems Global Technology Co Llc | データ処理装置、医用装置、データ処理方法、およびプログラム |
JP5988772B2 (ja) | 2012-01-20 | 2016-09-07 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
CN104363834B (zh) * | 2012-06-07 | 2017-08-29 | 东芝医疗系统株式会社 | 图像处理装置以及x射线诊断装置 |
JP6061555B2 (ja) | 2012-08-27 | 2017-01-18 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
WO2014167466A1 (en) * | 2013-04-09 | 2014-10-16 | Koninklijke Philips N.V. | Layered two-dimensional projection generation and display |
JP6507615B2 (ja) | 2013-12-13 | 2019-05-08 | 株式会社ニデック | 光コヒーレンストモグラフィ装置、及びプログラム |
AU2016222299B1 (en) * | 2015-07-09 | 2016-11-03 | Atomic Oncology Pty Ltd | Atomic Therapeutic Indicator |
US10102682B1 (en) * | 2017-04-17 | 2018-10-16 | Raytheon Company | System and method for combining 3D images in color |
-
2017
- 2017-04-17 US US15/489,652 patent/US10102682B1/en active Active
-
2018
- 2018-01-29 WO PCT/US2018/015693 patent/WO2018194739A1/en unknown
- 2018-01-29 KR KR1020197033771A patent/KR102275334B1/ko active IP Right Grant
- 2018-01-29 CA CA3060256A patent/CA3060256A1/en active Pending
- 2018-01-29 EP EP18706015.7A patent/EP3613019B1/en active Active
- 2018-01-29 CN CN201880040220.3A patent/CN110753954B/zh active Active
- 2018-01-29 EP EP23167521.6A patent/EP4220568A1/en active Pending
- 2018-01-29 JP JP2019556353A patent/JP6854919B2/ja active Active
- 2018-01-31 TW TW107103473A patent/TWI676024B/zh active
- 2018-09-21 US US16/138,764 patent/US10366545B2/en active Active
-
2019
- 2019-06-07 US US16/435,275 patent/US10650612B2/en active Active
- 2019-10-22 IL IL270076A patent/IL270076B/en active IP Right Grant
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020136439A1 (en) * | 2001-03-09 | 2002-09-26 | Ruchala Kenneth J. | System and method for fusion-aligned reprojection of incomplete data |
US20120038644A1 (en) * | 2009-02-17 | 2012-02-16 | Uws Ventures Ltd. | Rendering caustics in computer graphics with photon mapping by relaxation |
US20110135173A1 (en) * | 2009-10-11 | 2011-06-09 | Axs Ingenierie | Method for 3d spine and full 3d trunk imaging fusion in standing position |
TW201123076A (en) * | 2009-12-30 | 2011-07-01 | Pou Yuen Technology Co Ltd | Three-dimensional display method of medical images |
TWI517093B (zh) * | 2013-10-11 | 2016-01-11 | Univ Nat Yang Ming | Computer tomography reconstruction method |
TWI509564B (zh) * | 2014-10-30 | 2015-11-21 | Inst Nuclear Energy Res Atomic Energy Council | 三維成像的投影方法 |
Also Published As
Publication number | Publication date |
---|---|
KR102275334B1 (ko) | 2021-07-08 |
EP3613019B1 (en) | 2023-04-19 |
US20190295329A1 (en) | 2019-09-26 |
WO2018194739A1 (en) | 2018-10-25 |
US10650612B2 (en) | 2020-05-12 |
CN110753954A (zh) | 2020-02-04 |
US20190035164A1 (en) | 2019-01-31 |
KR20190139968A (ko) | 2019-12-18 |
EP4220568A1 (en) | 2023-08-02 |
CA3060256A1 (en) | 2018-10-25 |
TW201839389A (zh) | 2018-11-01 |
US10366545B2 (en) | 2019-07-30 |
JP6854919B2 (ja) | 2021-04-07 |
IL270076B (en) | 2020-09-30 |
US10102682B1 (en) | 2018-10-16 |
JP2020516413A (ja) | 2020-06-11 |
CN110753954B (zh) | 2023-08-08 |
US20180300956A1 (en) | 2018-10-18 |
EP3613019A1 (en) | 2020-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI676024B (zh) | 用以組合彩色3d影像之系統及方法 | |
JP4421016B2 (ja) | 医用画像処理装置 | |
US8698806B2 (en) | System and method for performing volume rendering using shadow calculation | |
US9053565B2 (en) | Interactive selection of a region of interest in an image | |
US9665990B2 (en) | Image display to display 3D image and sectional images | |
JP6602574B2 (ja) | 医用画像処理装置 | |
JP6688618B2 (ja) | 医用画像処理装置及び医用画像診断装置 | |
US11062512B2 (en) | System and method for generating 3D color representation of 2D grayscale images | |
JP2006055213A (ja) | 画像処理装置、及びプログラム | |
JP2008259699A (ja) | 画像処理方法および装置並びにプログラム | |
CN109313816B (zh) | 深度增强的断层合成重建 | |
US20060214930A1 (en) | Image processing method and computer readable medium | |
JP7003635B2 (ja) | コンピュータプログラム、画像処理装置及び画像処理方法 | |
JP5065740B2 (ja) | 画像処理方法および装置並びにプログラム | |
CN103282942B (zh) | 基于距离的体积可视化 | |
JP7476403B2 (ja) | 医療ボリュームデータをレンダリングするためのコンピュータ実施方法 | |
Chatzigiannis | 3D Volume Data in Nuclear Medicine | |
She et al. | An interactive editing tool from arbitrary slices in 3D ultrasound volume data | |
US20080238932A1 (en) | Method and device for adjusting color value for a low-noise visualization of volumetric objects | |
JP2006065479A (ja) | プログラム、及び画像処理装置 |