TWI679652B - Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging - Google Patents

Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging Download PDF

Info

Publication number
TWI679652B
TWI679652B TW108101383A TW108101383A TWI679652B TW I679652 B TWI679652 B TW I679652B TW 108101383 A TW108101383 A TW 108101383A TW 108101383 A TW108101383 A TW 108101383A TW I679652 B TWI679652 B TW I679652B
Authority
TW
Taiwan
Prior art keywords
image
gray matter
voxel
brain
template
Prior art date
Application number
TW108101383A
Other languages
Chinese (zh)
Other versions
TW202027091A (en
Inventor
周坤賢
Kun Hsien Chou
林慶波
Ching Po Lin
林偉哲
Wei Che Lin
Original Assignee
國立陽明大學
National Yang-Ming University
長庚醫療財團法人高雄長庚紀念醫院
Kaohsiung Chang Gung Memorial Hospital
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 國立陽明大學, National Yang-Ming University, 長庚醫療財團法人高雄長庚紀念醫院, Kaohsiung Chang Gung Memorial Hospital filed Critical 國立陽明大學
Priority to TW108101383A priority Critical patent/TWI679652B/en
Application granted granted Critical
Publication of TWI679652B publication Critical patent/TWI679652B/en
Publication of TW202027091A publication Critical patent/TW202027091A/en

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

一種個人化大腦影像評估之方法,其包含下列步驟:接收病人的大腦之第一T1加權影像;對第一T1加權影像進行組織分割程序,以取得第一灰質影像;將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板;依據配准至標準化正常灰質影像模板的第一灰質影像透過體素分析公式取得其具有的每一體素對應的評分值;判斷每一體素對應的評分值是否大於或小於預期值;以及依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告。因此,可在排除相關干擾因子下得知單一病人與健康受試者之間大腦影像之差異情形。A method for personalizing brain image evaluation includes the following steps: receiving a first T1-weighted image of a patient's brain; performing a tissue segmentation procedure on the first T1-weighted image to obtain a first gray matter image; and registering the first gray matter image To the normalized normal gray matter image template created by the superimposition procedure; according to the first gray matter image registered to the normalized normal gray matter image template, a voxel analysis formula is used to obtain the score value corresponding to each voxel; it is determined that each voxel corresponds to Whether the value of the score is greater than or less than the expected value; and based on those voxels greater than or less than the expected value, the brain region to which it belongs is pushed back, and an evaluation report is output. Therefore, the differences in brain imaging between a single patient and a healthy subject can be known by excluding relevant interference factors.

Description

個人化大腦影像評估之方法、非暫時性電腦可讀媒體及設備Method for personalizing brain image assessment, non-transitory computer-readable media and equipment

本發明涉及一種影像評估之方法、非暫時性電腦可讀媒體及設備,特別是個人化大腦影像評估之方法、非暫時性電腦可讀媒體及設備。The invention relates to a method for image evaluation, a non-transitory computer-readable medium and a device, and more particularly to a method for evaluating a personalized brain image, a non-transitory computer-readable medium and a device.

目前常見於影像的統計分析方法為T檢定(T-test),其為多對多的群組分析,例如:一群病人與一群健康人之間的比較,若將此分析方式套用於個人化評估,需將單一病人當成第一群組,把其他所有的健康受試者當成第二群組,然而在T檢定的統計假設中,需假設兩個群組具有相同的平均(mean)及方差(variance),因此,當應用在個人化評估時,由於第一群組只有一個人,存在不符合此統計假設的問題。此外,若仍將此分析方式套用於個人化評估時,會因第一群組只有一個人造成統計模型及資料庫樣本數不足,進而導致靈敏度下降之問題。At present, the statistical analysis method commonly used in imaging is T-test, which is a many-to-many group analysis, for example: comparison between a group of patients and a group of healthy people. If this analysis method is applied to personalized evaluation It is necessary to treat a single patient as the first group and all other healthy subjects as the second group. However, in the statistical test of the T test, it is necessary to assume that the two groups have the same mean and variance ( variance). Therefore, when applied to personalized evaluation, there is a problem that does not meet this statistical assumption because there is only one person in the first group. In addition, if this analysis method is still applied to the personalized evaluation, the statistical model and the number of database samples will be insufficient because there is only one person in the first group, which will cause the problem of reduced sensitivity.

綜上所述,可知先前技術中存在T檢定因第一群組只有一個人造成統計模型及資料庫樣本數不足進而導致靈敏度下降之問題,因此實有必要提出改進的技術手段,來解決此一問題。In summary, it can be known that the T test in the prior art has a problem that the statistical model and the number of database samples are insufficient due to only one person in the first group, which reduces the sensitivity. Therefore, it is necessary to propose improved technical methods to solve this problem .

本發明揭露一種個人化大腦影像評估之方法、非暫時性電腦可讀媒體及設備。The invention discloses a method, a non-transitory computer-readable medium and a device for personalizing brain image evaluation.

首先,本發明揭露一種個人化大腦影像評估之方法,其步驟包含:接收一病人的大腦之第一T1加權影像(T1 weighted image);對第一T1加權影像進行組織分割程序,以取得第一灰質影像;將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板;將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值,其中,體素分析公式為 ,TW為第i個體素對應的評分值,x*為配准至標準化正常灰質影像模板的第一灰質影像之第i個體素的體積大小, 為依據該病人的年齡與性別於資料庫取得預期影像之第i個體素的體積大小,s為依據資料庫取得估計殘差影像之第i個體素的體積標準差, 為T檢定的比例因子,n為建立資料庫的樣本數,i為大於或等於1的正整數; 判斷每一體素對應的評分值是否大於或小於預期值;依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告。 First, the present invention discloses a method for personalizing brain image evaluation. The steps include: receiving a first T1 weighted image of a patient's brain; and performing a tissue segmentation procedure on the first T1 weighted image to obtain a first Gray matter image; register the first gray matter image to the normalized normal gray matter image template created by the superposition procedure; obtain the first gray matter image registered to the normalized normal gray matter image template according to the voxel analysis formula to obtain each voxel it has Corresponding score value, where the voxel analysis formula is , TW is the score value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized gray matter image template, In order to obtain the volume of the i-th voxel of the expected image from the database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel of the estimated residual image from the database, Is the scale factor for the T test, n is the number of samples to build the database, and i is a positive integer greater than or equal to 1; determine whether the score value corresponding to each voxel is greater than or less than the expected value; based on those greater or less than the expected value The voxel pushes back the brain region to which it belongs and outputs an evaluation report.

另外,本發明揭露一種個人化大腦影像評估之非暫時性電腦可讀媒體,其經組態以儲存若干指令,該等指令在由一或多個處理器執行時使得該一或多個處理器執行以下操作:接收一病人的大腦之第一T1加權影像;對第一T1加權影像進行組織分割程序,以取得第一灰質影像;將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板;將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值,其中,體素分析公式為 ,TW為第i個體素對應的評分值,x*為配准至正常影像模板的第一灰質影像之第i個體素的體積大小, 為依據該病人的年齡與性別於資料庫取得預期影像之第i個體素的體積大小,s為依據資料庫取得估計殘差影像之第i個體素的體積標準差, 為T檢定的比例因子,n為建立資料庫的樣本數,i為大於或等於1的正整數; 判斷每一體素對應的評分值是否大於或小於預期值;依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告。 In addition, the present invention discloses a non-transitory computer-readable medium for personalized brain imaging assessment, which is configured to store a number of instructions that, when executed by one or more processors, cause the one or more processors to The following operations are performed: receiving a first T1-weighted image of a patient's brain; performing a tissue segmentation procedure on the first T1-weighted image to obtain a first gray matter image; registering the first gray matter image to a normalization established by the superposition procedure Normal gray matter image template; the first gray matter image registered to the normalized normal gray matter image template is used to obtain the score value corresponding to each voxel it has, where the voxel analysis formula is , TW is the score value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normal image template, In order to obtain the volume of the i-th voxel of the expected image from the database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel of the estimated residual image from the database, Is the scale factor for the T test, n is the number of samples to build the database, and i is a positive integer greater than or equal to 1; determine whether the score value corresponding to each voxel is greater than or less than the expected value; based on those greater or less than the expected value The voxel pushes back the brain region to which it belongs and outputs an evaluation report.

再者,本發明揭露一種個人化大腦影像評估之設備,此設備包含:一或多個處理器、儲存單元以及至少一程式,其中該至少一程式儲存於儲存單元中且經組態以由該一或多個處理器執行,該至少一程式整體上包含用於以下操作指令:接收一病人的大腦之第一T1加權影像;對第一T1加權影像進行組織分割程序,以取得第一灰質影像;將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板;將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值,其中,體素分析公式為 ,TW為第i個體素對應的評分值,x*為配准至標準化正常灰質影像模板的第一灰質影像之第i個體素的體積大小, 為依據該病人的年齡與性別於資料庫取得預期影像之第i個體素的體積大小,s為依據資料庫取得估計殘差影像之第i個體素的體積標準差, 為T檢定的比例因子,n為建立資料庫的樣本數,i為大於或等於1的正整數; 判斷每一體素對應的評分值是否大於或小於預期值;依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告。 Furthermore, the present invention discloses a device for personalizing brain image evaluation. The device includes: one or more processors, a storage unit, and at least one program, wherein the at least one program is stored in the storage unit and configured to be executed by the device. Executed by one or more processors, the at least one program generally includes the following operation instructions: receiving a first T1 weighted image of a patient's brain; performing a tissue segmentation procedure on the first T1 weighted image to obtain a first gray matter image ; Register the first gray matter image to the normalized normal gray matter image template created by the superposition procedure; obtain the first gray matter image registered to the normalized normal gray matter image template according to the voxel analysis formula to obtain the corresponding Scoring value, where the voxel analysis formula is , TW is the score value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized gray matter image template, In order to obtain the volume of the i-th voxel of the expected image from the database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel of the estimated residual image from the database, Is the scale factor for the T test, n is the number of samples to build the database, and i is a positive integer greater than or equal to 1; determine whether the score value corresponding to each voxel is greater than or less than the expected value; The voxel pushes back the brain region to which it belongs and outputs an evaluation report.

本發明所揭露之系統與方法如上,與先前技術的差異在於本發明是透過接收病人的大腦之第一T1加權影像;對第一T1加權影像進行組織分割程序,以取得第一灰質影像;將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板;依據配准至標準化正常灰質影像模板的第一灰質影像透過體素分析公式取得其具有的每一體素對應的評分值;判斷每一體素對應的評分值是否大於或小於預期值;以及依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告。The system and method disclosed in the present invention are as above. The difference from the prior art is that the present invention receives the first T1 weighted image of the patient's brain; performs a tissue segmentation procedure on the first T1 weighted image to obtain the first gray matter image; The first gray matter image is registered to the normalized normal gray matter image template created by the superimposition process; the first gray matter image registered to the normalized normal gray matter image template is used to obtain the corresponding score value for each voxel it has through the voxel analysis formula. ; Determine whether the score value corresponding to each voxel is greater than or less than the expected value; and push back the brain region to which it belongs according to those voxels greater than or less than the expected value, and output an evaluation report.

透過上述的技術手段,本發明可在排除相關干擾因子(年齡因子與性別因子)下得知單一病人與健康受試者之間大腦影像之差異情形,有助於大腦影像數據化以及偵測肉眼察覺不到之細微大腦結構的改變。Through the above-mentioned technical means, the present invention can know the difference between the brain image between a single patient and a healthy subject without excluding relevant interference factors (age factor and gender factor), which is helpful for data quantification of the brain image and detection of the naked eye. Undetectable subtle changes in brain structure.

以下將配合圖式及實施例來詳細說明本發明之實施方式,藉此對本發明如何應用技術手段來解決技術問題並達成技術功效的實現過程能充分理解並據以實施。需注意的是,本案所述之影像皆為三維影像,為避免圖式過於複雜,僅以二維影像呈現每一實施例的示意圖。此外,儘管這裡可以使用詞彙“第一”、“第二”等來描述各種影像,但是這些影像不應受這些詞彙的限制。In the following, the embodiments of the present invention will be described in detail with reference to the drawings and examples, so as to fully understand and implement the implementation process of how the present invention applies technical means to solve technical problems and achieve technical effects. It should be noted that the images described in this case are all three-dimensional images. In order to avoid the drawings being too complicated, only the two-dimensional images are used to present the schematic diagrams of each embodiment. In addition, although the terms "first", "second", etc. may be used herein to describe various images, these images should not be limited by these terms.

請先參閱「第1圖」,「第1圖」為本發明個人化大腦影像評估之方法的一實施例方法流程圖,其步驟包含:接收一病人的大腦之第一T1加權影像(步驟110);對第一T1加權影像進行組織分割程序,以取得第一灰質影像(步驟120);將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板(步驟130);將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值,其中,體素分析公式為 ,TW為第i個體素對應的評分值,x*為配准至標準化正常灰質影像模板的第一灰質影像之第i個體素的體積大小, 為依據該病人的年齡與性別於資料庫取得預期影像之第i個體素的體積大小,s為依據資料庫取得估計殘差影像之第i個體素的體積標準差, 為T檢定的比例因子,n為建立資料庫的樣本數,i為大於或等於1的正整數(步驟140);判斷每一體素對應的評分值是否大於或小於預期值;(步驟150);依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告(步驟160)。 Please refer to "Figure 1" first. "Figure 1" is a method flowchart of an embodiment of a method for personalizing brain image assessment according to the present invention. The steps include: receiving a first T1-weighted image of a patient's brain (step 110). ); Perform a tissue segmentation procedure on the first T1 weighted image to obtain a first gray matter image (step 120); register the first gray matter image to a normalized normal gray matter image template created by the superposition procedure (step 130); The first gray matter image registered to the normalized normal gray matter image template obtains a score value corresponding to each voxel it has according to a voxel analysis formula, where the voxel analysis formula is , TW is the score value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized gray matter image template, In order to obtain the volume of the i-th voxel of the expected image from the database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel of the estimated residual image from the database, Is the scale factor for the T test, n is the number of samples to build the database, i is a positive integer greater than or equal to 1 (step 140); determine whether the score value corresponding to each voxel is greater than or less than the expected value; (step 150); According to those voxels that are larger or smaller than the expected value, the brain region to which they belong is pushed back, and an evaluation report is output (step 160).

步驟110所述之病人可為但不限於欲進行大腦可能存在異常狀態評估之病患,該病患透過上述步驟110至步驟160後可在排除相關干擾因子下得知其與健康受試者之間大腦影像之差異情形,且有助於大腦影像數據化以及偵測肉眼察覺不到之細微大腦結構的改變。其中,該病人為成年人。The patient described in step 110 may be, but is not limited to, a patient who may be involved in the assessment of a possible abnormal state of the brain. After passing through the above steps 110 to 160, the patient can learn that he is related to healthy subjects by excluding relevant interference factors. The differences between brain images can help to digitize brain images and detect subtle changes in brain structure that cannot be detected by the naked eye. Among them, the patient is an adult.

步驟120所述之組織分割程序可為但不限於利用開源軟體FSL對第一T1加權影像進行組織分割,可取得該病人的灰質影像(即第一灰質影像)、白質影像與腦脊髓液影像,但本實施例並非用以限定本發明,可依據實際需求進行調整。舉例而言,也可透過統計參數圖譜(Statistical Parametric Mapping,SPM)的腦功能磁共振資料處理技術取得該病人的灰質影像(即第一灰質影像)。The tissue segmentation procedure described in step 120 may be, but is not limited to, using the open source software FSL to perform tissue segmentation on the first T1 weighted image to obtain the patient's gray matter image (ie, the first gray matter image), white matter image, and cerebrospinal fluid image. However, this embodiment is not intended to limit the present invention, and can be adjusted according to actual needs. For example, the gray matter image (ie, the first gray matter image) of the patient may also be obtained through a functional magnetic resonance data processing technique of statistical parametric mapping (SPM).

請參閱「第1圖」與「第2圖」,「第2圖」為「第1圖」之步驟130所述的疊合程序之一實施例方法流程圖。在本實施例中,步驟130所述之疊合程序包含以下步驟:取得n個健康受試者的大腦之第二T1加權影像,n為大於或等於三十之正整數(步驟210);對每一第二T1加權影像分別進行組織分割程序,以取得對應的第二灰質影像(步驟220);將該些第二灰質影像分別配准至MNI( Montreal Neurological Institute)標準空間模板(步驟230);以及將配准至MNI標準空間模板的該些第二灰質影像進行疊合,並以影像亮度取平均值的方式建立標準化正常灰質影像模板(步驟240)。其中,MNI標準空間模板為國際腦成像領域廣泛使用的標準腦模板。Please refer to "Fig. 1" and "Fig. 2". "Fig. 2" is a flowchart of one embodiment of the superimposing procedure described in step 130 of "Fig. 1". In this embodiment, the superposition procedure described in step 130 includes the following steps: obtaining second T1-weighted images of the brains of n healthy subjects, where n is a positive integer greater than or equal to thirty (step 210); Each second T1 weighted image is subjected to a tissue segmentation procedure to obtain corresponding second gray matter images (step 220); the second gray matter images are respectively registered to the MNI (Montreal Neurological Institute) standard space template (step 230) ; And superimposing the second gray matter images registered to the MNI standard space template, and establishing a normalized normal gray matter image template by averaging the image brightness (step 240). Among them, the MNI standard space template is a standard brain template widely used in the field of international brain imaging.

由於本實施例係用以對成年人進行個人化大腦影像評估,因此,步驟210所述之健康受試者為健康成年人。Since this embodiment is used to perform personalized brain image evaluation on adults, the healthy subjects described in step 210 are healthy adults.

步驟240所述之以影像亮度取平均值的方式建立標準化正常灰質影像模板,係透過配准至MNI標準空間模板的每一第二灰質影像之每一體素具有其亮度值,當配准至MNI標準空間模板的該些第二灰質影像進行疊合後,可藉由不同第二灰質影像之第i個體素(i為大於或等於1的正整數)的亮度值相加後取平均值,而取得正常影像模板之第i個體素的亮度值,進而建立標準化正常灰質影像模板。In step 240, the normalized normal gray matter image template is established by averaging the image brightness. Each voxel of each second gray matter image registered to the MNI standard space template has its brightness value. When registered to the MNI, After superimposing the second gray matter images of the standard space template, the luminance values of the i-th voxel (i is a positive integer greater than or equal to 1) of different second gray matter images can be added to obtain the average value, and Obtain the brightness value of the i-th voxel of the normal image template, and then establish a normalized normal gray matter image template.

由於上述n個健康受試者皆可為台灣人,故可透過步驟210至步驟240取得適用於台灣的標準空間模板(即標準化正常灰質影像模板)。因此,可藉由不同地區的健康受試者的大腦之第一T1加權影像建立適用於該地區的標準空間模板。Since the above n healthy subjects can all be Taiwanese, a standard space template (ie, a normalized normal gray matter image template) applicable to Taiwan can be obtained through steps 210 to 240. Therefore, the first T1-weighted images of the brains of healthy subjects in different regions can be used to establish a standard spatial template suitable for that region.

請參閱「第1圖」與「第3圖」,「第3圖」為「第1圖」之步驟140所述的資料庫的建立方法之一實施例方法流程圖。在本實施例中,由於在執行步驟140所述之「將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值」時,需藉由資料庫取得體素分析公式所需的 與s,因此,在執行步驟140之前要先建立該資料庫,而該資料庫的建立方法包含以下步驟:將每一第二灰質影像配准至標準化正常灰質影像模板(步驟310);對配准至標準化正常灰質影像模板的每一第二灰質影像進行平滑化處理(步驟320);將平滑化的該些第二灰質影像依據一般線性模型並整合年齡因子與性別因子,估算出年齡因子的權重係數與性別因子的權重係數,並取得平滑化的每一第二灰質影像對應的殘差影像(步驟330);以及將每一殘差影像的第i個體素的體積大小進行標準差計算,以取得估計殘差影像的第i個體素之體積標準差,進而建立資料庫(步驟340)。 Please refer to "Figure 1" and "Figure 3", and "Figure 3" is a method flowchart of an embodiment of the method for establishing a database described in step 140 of "Figure 1". In this embodiment, since “the first gray matter image registered to the normalized normal gray matter image template is obtained according to the voxel analysis formula to obtain the score value corresponding to each voxel” described in step 140, it is necessary to borrow Required for obtaining voxel analysis formula from database And s, therefore, the database must be established before executing step 140, and the method of establishing the database includes the following steps: registering each second gray matter image to a normalized normal gray matter image template (step 310); Each second gray matter image normalized to a normal gray matter image template is smoothed (step 320); the smoothed second gray matter images are based on a general linear model and integrate age factors and gender factors to estimate the age factor. A weighting coefficient and a weighting coefficient of a gender factor, and obtaining a smoothed residual image corresponding to each second gray matter image (step 330); and calculating a standard deviation of the volume size of the i-th voxel of each residual image, To obtain the volume standard deviation of the i-th voxel of the estimated residual image, and then establish a database (step 340).

在本實施例中,步驟320所述之平滑化處理可為但不限於三維高斯平滑化處理,可依據實際需求進行調整。其中,平滑化處理可減少第二灰質影像配准至標準化正常灰質影像模板時的誤差。In this embodiment, the smoothing process described in step 320 may be, but is not limited to, a three-dimensional Gaussian flattening process, and may be adjusted according to actual needs. Among them, the smoothing process can reduce the error when the second gray matter image is registered to the normalized normal gray matter image template.

由於年齡及性別係為會影響大腦變化的相關因子,因此,在步驟330中,可在使用一般線性模型對平滑化的該些第二灰質影像進行分析時,整合年齡因子與性別因子,以估算出年齡因子的權重係數與性別因子的權重係數,並取得平滑化的每一第二灰質影像對應的殘差影像。其中,估算出來的年齡因子的權重係數與性別因子的權重係數可分別為一數值範圍。接著,將每一殘差影像的第i個體素的體積大小進行標準差計算,以取得估計殘差影像的第i個體素之體積標準差(即體素分析公式的s),進而建立資料庫(即步驟340)。因此,資料庫係為利用基於年齡因子與性別因子的體素線性回歸之資料庫。在本實施例中,使用一般線性模型進行分析所整合的因子可為年齡因子與性別因子,但本實施例並非用以限定本發明。舉例而言,使用一般線性模型進行分析所整合的因子還可針對各種會影響大腦變化之因子進行整合,例如:教育程度或認知功能表現等,可依據實際需求進行調整。Since age and gender are related factors that affect brain changes, in step 330, when the smoothed second gray matter images are analyzed using a general linear model, the age factor and gender factor are integrated to estimate The weight coefficient of the age factor and the weight factor of the gender factor are obtained, and a residual image corresponding to each smoothed second gray matter image is obtained. Among them, the weight coefficients of the estimated age factor and the gender factor may be a numerical range, respectively. Next, calculate the standard deviation of the volume of the i-th voxel of each residual image to obtain the estimated standard deviation of the volume of the i-th voxel of the residual image (ie, s of the voxel analysis formula), and then establish a database (Ie step 340). Therefore, the database is a database using linear regression of voxels based on age and gender factors. In this embodiment, the factors integrated by using a general linear model for analysis may be an age factor and a gender factor, but this embodiment is not intended to limit the present invention. For example, using a general linear model to analyze the integrated factors can also integrate various factors that affect brain changes, such as education level or cognitive function performance, etc., which can be adjusted according to actual needs.

透過步驟310至步驟340建立資料庫後,執行步驟140所述之「將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值」。其中, 係為將該病人的年齡與性別導入利用基於年齡因子與性別因子的體素線性回歸之資料庫後所取得的預期影像的第i個體素的體積大小,s為依據資料庫取得估計殘差影像之第i個體素的體積標準差。 After the database is established through steps 310 to 340, the "registration of the first gray matter image registered to the normalized normal gray matter image template according to the voxel analysis formula to obtain the score value corresponding to each voxel" described in step 140 is performed. among them, This is the volume of the i-th voxel of the expected image obtained by using the age and gender voxel linear regression database based on the age and sex of the patient, and s is the estimated residual image obtained from the database The volume standard deviation of the i-th voxel.

請繼續參閱「第1圖」,在步驟150中,預期值可為但不限於0.05,可依據實際需求進行調整。Please continue to refer to "Figure 1". In step 150, the expected value can be, but not limited to, 0.05, and can be adjusted according to actual needs.

在步驟160中,大於或小於預設值的該些體素所屬的腦區可為但不限於額葉(Frontal Lobe)、頂葉(Parietal lobe)、枕葉(Occipital Lobe)、顳葉(Temporal lobe)、島葉(Insula)、海馬體(Hippocampus)或杏仁體(Amygdala)。評估報告的內容可包括透過可視化方式將評分值大於或小於預期值的該些體素標註其上之第一灰質影像。因此,透過上述步驟110至步驟160,可在排除相關干擾因子(年齡因子與性別因子)下得知單一病人與健康受試者之間大腦影像之差異情形,有助於大腦影像數據化以及偵測肉眼察覺不到之細微大腦結構的改變。In step 160, the brain regions to which these voxels larger or smaller than a preset value belong may be, but not limited to, Frontal Lobe, Parietal Lobe, Occipital Lobe, and Temporal lobe), island leaves (Insula), hippocampus (Hippocampus) or amygdala (Amygdala). The content of the evaluation report may include visually marking the first gray matter image on the voxels whose score value is greater than or less than the expected value. Therefore, through the above steps 110 to 160, the difference between the brain image between a single patient and a healthy subject can be learned under the exclusion of relevant interference factors (age factor and gender factor), which is helpful for the data visualization and detection of brain images. Changes in subtle brain structures that are not detectable to the naked eye.

接著,請參閱「第4圖」,「第4圖」為本發明個人化大腦影像評估之非暫時性電腦可讀媒體之一實施例方塊圖。個人化大腦影像評估之非暫時性電腦可讀媒體400包含可由一或多個處理器404用於執行包含上文所闡述的大腦影像評估之方法之操作指令402。該等操作指令可包含本文中所闡述之任何(一或多個)其他步驟。實施諸如本文中所闡述之操作指令402可儲存於個人化大腦影像評估之非暫時性電腦可讀媒體400上。個人化大腦影像評估之非暫時性電腦可讀媒體400可係諸如一磁碟或光碟或一磁帶之一儲存媒體,或此項技術中已知之任何其他適合之非暫時性電腦可讀媒體。Next, please refer to "Figure 4". "Figure 4" is a block diagram of an embodiment of a non-transitory computer-readable medium for personalized brain image evaluation according to the present invention. The non-transitory computer-readable medium 400 for personalized brain imaging assessment includes operating instructions 402 that can be used by one or more processors 404 to perform methods including the brain imaging assessment described above. Such operating instructions may include any one or more other steps set forth herein. Implementing operating instructions 402 such as those described herein may be stored on a non-transitory computer-readable medium 400 of a personalized brain imaging assessment. The non-transitory computer-readable medium 400 for personalized brain imaging assessment may be a storage medium such as a magnetic disk or optical disk or a magnetic tape, or any other suitable non-transitory computer-readable medium known in the art.

接著,請參閱「第5圖」,「第5圖」為本發明個人化大腦影像評估之設備之一實施例設備方塊圖。個人化大腦影像評估之設備500包含處理器502、儲存單元504以及程式506,其中程式506儲存於儲存單元504中且經組態以由處理器502執行,程式506整體上可包含上文所闡述的大腦影像評估之方法之操作指令。其中,處理器502、儲存單元504以及程式506的數量可為但不限於一,可依據實際需求進行調整。Next, please refer to "FIG. 5". "FIG. 5" is a block diagram of an embodiment of a device for personalizing brain image evaluation according to the present invention. The device 500 for personalized brain imaging evaluation includes a processor 502, a storage unit 504, and a program 506. The program 506 is stored in the storage unit 504 and configured to be executed by the processor 502. The program 506 as a whole may include the above-mentioned Instructions for the method of brain image evaluation. The number of the processor 502, the storage unit 504, and the program 506 may be, but is not limited to, one, and may be adjusted according to actual needs.

以下配合「第6圖」至「第7B圖」以實例的方式進行如下說明。The following description is made by way of example in conjunction with "Figure 6" to "Figure 7B".

(健康成人資料庫之樣本數對個人化大腦影像評估的穩定度分析)(Stability analysis of samples from a healthy adult database for personalized brain imaging assessment)

由於個人化大腦影像評估之方法係為一種評量單一病人的大腦灰質影像相較於透過正常健康受試者的大腦灰質影像所建立的資料庫之變化程度的方法,因此,用以建立資料庫之正常健康受試者的大腦灰質影像之樣本數可能會影響評估單一病人的大腦灰質影像的穩定性。請參閱「第6圖」,「第6圖」為不同樣本數所建立的資料庫對同一病人進行個人化大腦影像評估之方法後輸出的評估報告之第一灰質影像。其中,n為建立資料庫的樣本數,每一第一灰質影像係透過可視化方式將評估值大於或小於預期值的體素標註其上。在本實施例中,分別將61歲女性阿茲海默氏症患者的大腦灰質影像利用不同樣本數所建立的資料庫進行個人化大腦影像評估,由「第6圖」可知,當用以建立資料庫之正常健康受試者的大腦灰質影像之樣本數為30人時即可跟樣本數為152人所建立的資料庫達到極為相近的成果。換句話說,個人化大腦影像評估之方法可以利用小樣本數所建立的資料庫得到極為穩定的統計成果。Because the method of personalized brain image assessment is a method for evaluating the degree of change of the gray matter image of a single patient compared to a database established through the gray matter image of a normal healthy subject, it is used to establish a database The number of gray matter images of normal healthy subjects may affect the evaluation of the gray matter image stability of a single patient. Please refer to "Figure 6", "Figure 6" is the first gray matter image of the evaluation report output by the personalized brain image evaluation method for the same patient based on the database created by different sample sizes. Among them, n is the number of samples for establishing a database, and each first gray matter image is visually labeled with voxels whose evaluation value is greater than or less than the expected value. In this example, the brain gray matter images of 61-year-old female Alzheimer's patients were individually evaluated using a database created by different sample numbers. As shown in "Figure 6", when used to establish When the number of gray matter images of normal healthy subjects in the database is 30, the results can be very similar to the database established by the sample number of 152. In other words, the method of personalized brain image evaluation can use the database established by a small number of samples to obtain extremely stable statistical results.

(個人化大腦影像評估之方法的靈敏度分析)(Sensitivity analysis of personalized brain imaging assessment methods)

在本實施例中,係以偵測海馬體體積的萎縮之靈敏度進行討論。請參閱「第7A圖」與「第7B圖」,「第7A圖」與「第7B圖」分別為利用「第1圖」的個人化大腦影像評估之方法與傳統T檢定統計方式對同一病人的海馬體進行萎縮電腦模擬評估之示意圖。其中,縱軸為海馬體萎縮百分比;橫軸為對應群聚多重比較(clustering-wise multiple comparison)校正的體素數量;淺色區域為未校正的P值小於0.05時的情況;深色區域為群聚層次校正(cluster level corrected)的P值小於0.05時的情況。在本實施例中,把P值小於0.05當作統計上具有顯著意義來看。In this embodiment, the sensitivity of detecting hippocampal volume atrophy is discussed. Please refer to "Figure 7A" and "Figure 7B", "Figure 7A" and "Figure 7B" respectively use the method of personalized brain image evaluation of "Figure 1" and the traditional T test statistical method for the same patient Schematic illustration of computer simulation of hippocampal atrophy. Among them, the vertical axis is the percentage of hippocampal atrophy; the horizontal axis is the number of voxels corrected for corresponding clustering-wise multiple comparison; the light-colored area is when the uncorrected P value is less than 0.05; the dark-colored area is When the cluster level corrected P value is less than 0.05. In the present embodiment, a P value of less than 0.05 is regarded as statistically significant.

本實施例之個人化大腦影像評估之方法係利用樣本數為152人所建立的資料庫進行評估,透過「第7A圖」可知,個人化大腦影像評估之方法可在海馬體體積產生9%的萎縮時就被偵測出該病人有異常的體積下降現象;而傳統T檢定統計方式則要在海馬體體積產生12%的萎縮時才可以被偵測出該病人有異常的體積下降現象。因此,可知個人化大腦影像評估之方法在單一病人大腦灰質體積大小的偵測上可以提供更好的靈敏度。The method of personalizing brain image evaluation in this embodiment uses a database established by a sample size of 152 people for evaluation. According to "Figure 7A", the method of personalizing brain image evaluation can generate 9% of the hippocampal volume. The patient was detected to have an abnormal volume decrease during atrophy; the traditional T-test statistical method can only be detected when the hippocampal body has a 12% atrophy before the patient has detected an abnormal volume decrease. Therefore, it can be seen that the method of personalized brain imaging assessment can provide better sensitivity in detecting the gray matter volume of a single patient's brain.

綜上所述,可知本發明與先前技術之間的差異在於透過接收病人的大腦之第一T1加權影像;對第一T1加權影像進行組織分割程序,以取得第一灰質影像;將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板;依據配准至標準化正常灰質影像模板的第一灰質影像透過體素分析公式取得其具有的每一體素對應的評分值;判斷每一體素對應的評估值是否大於或小於預期值;以及依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告,藉由此一技術手段可以解決先前技術所存在的問題,進而可在排除相關干擾因子(年齡因子與性別因子)下得知單一病人與健康受試者之間大腦影像之差異情形,有助於大腦影像數據化以及偵測肉眼察覺不到之細微大腦結構的改變。In summary, the difference between the present invention and the prior art is that the first T1 weighted image of the patient's brain is received; the tissue segmentation procedure is performed on the first T1 weighted image to obtain the first gray matter image; The image is registered to the normalized normal gray matter image template created by the superimposition procedure; the first gray matter image registered to the normalized normal gray matter image template is used to obtain the corresponding score value of each voxel it has through the voxel analysis formula; Whether the evaluation value corresponding to the integrated element is greater than or less than the expected value; and based on those voxels greater than or less than the expected value, the brain region to which the voxel belongs is pushed back and an evaluation report is output. The problem of the brain image can be obtained by excluding related interference factors (age factor and gender factor), which can help to visualize the brain image and detect the invisible to the naked eye. Subtle changes in brain structure.

雖然本發明以前述之實施例揭露如上,然其並非用以限定本發明,任何熟習相像技藝者,在不脫離本發明之精神和範圍內,當可作些許之更動與潤飾,因此本發明之專利保護範圍須視本說明書所附之申請專利範圍所界定者為準。Although the present invention is disclosed in the foregoing embodiments as above, it is not intended to limit the present invention. Any person skilled in similar arts can make some modifications and retouches without departing from the spirit and scope of the present invention. The scope of patent protection shall be determined by the scope of the patent application attached to this specification.

400‧‧‧個人化大腦影像評估之非暫時性電腦可讀媒體400‧‧‧ Non-transitory computer-readable media for personalized brain imaging assessment

402‧‧‧操作指令 402‧‧‧operation instructions

404、502‧‧‧處理器 404, 502‧‧‧ processors

500‧‧‧個人化大腦影像評估之設備 500‧‧‧ Personalized Brain Imaging Evaluation Equipment

504‧‧‧儲存單元 504‧‧‧Storage unit

506‧‧‧程式 506‧‧‧program

步驟110‧‧‧接收一病人的大腦之第一T1加權影像 Step 110‧‧‧ receives a first T1-weighted image of the brain of a patient

步驟120‧‧‧對第一T1加權影像進行組織分割程序,以取得第一灰質影像 Step 120‧‧‧ Perform a tissue segmentation procedure on the first T1 weighted image to obtain a first gray matter image

步驟130‧‧‧將第一灰質影像配准至經疊合程序所建立的標準化正常灰質影像模板 Step 130‧‧‧ Register the first gray matter image to the normalized normal gray matter image template created by the superposition process

步驟140‧‧‧將配准至標準化正常灰質影像模板的第一灰質影像依據體素分析公式取得其具有的每一體素對應的評分值,其中,體素分析公式為,TW為第i個體素對應的評分值,x*為配准至標準化正常灰質影像模板的第一灰質影像之第i個體素的體積大小,為依據該病人的年齡與性別於資料庫取得預期影像之第i個體素的體積大小,s為依據資料庫取得估計殘差影像之第i個體素的體積標準差,為T檢定的比例因子,n為建立資料庫的樣本數,i為大於或等於1的正整數In step 140‧‧‧, the first gray matter image registered to the normalized normal gray matter image template is used to obtain a score value corresponding to each voxel according to the voxel analysis formula, where the voxel analysis formula is , TW is the score value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized gray matter image template, In order to obtain the volume of the i-th voxel of the expected image from the database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel of the estimated residual image from the database, Is the scale factor for the T test, n is the number of samples to build the database, and i is a positive integer greater than or equal to 1

步驟150‧‧‧判斷每一體素對應的評分值是否大於或小於預期值 Step 150‧‧‧ determines whether the score value corresponding to each voxel is greater than or less than the expected value

步驟160‧‧‧依據大於或小於預期值的該些體素回推其所屬的腦區,並輸出評估報告 Step 160‧‧‧ Push back the corresponding brain region according to those voxels that are larger or smaller than the expected value, and output an evaluation report

步驟210‧‧‧取得n個健康受試者的大腦之第二T1加權影像,n為大於或等於三十之正整數 Step 210‧‧‧ Obtain second T1-weighted images of the brains of n healthy subjects, where n is a positive integer greater than or equal to thirty

步驟220‧‧‧對每一第二T1加權影像分別進行組織分割程序,以取得對應的第二灰質影像 Step 220‧‧‧ performs a tissue segmentation procedure on each second T1 weighted image to obtain a corresponding second gray matter image

步驟230‧‧‧將該些第二灰質影像分別配准至MNI標準空間模板 Step 230‧‧‧ respectively register the second gray matter images to the MNI standard space template

步驟240‧‧‧將配准至MNI標準空間模板的該些第二灰質影像進行疊合,並以影像亮度取平均值的方式建立標準化正常灰質影像模板 Step 240‧‧‧ superimpose the second gray matter images registered to the MNI standard space template, and establish a normalized normal gray matter image template by averaging the image brightness

步驟310‧‧‧將每一第二灰質影像配准至標準化正常灰質影像模板 Step 310‧‧‧ Register each second gray matter image to a normalized normal gray matter image template

步驟320‧‧‧對配准至標準化正常灰質影像模板的每一第二灰質影像進行平滑化處理 Step 320‧‧‧ smoothing each second gray matter image registered to the normalized normal gray matter image template

步驟330‧‧‧將平滑化的該些第二灰質影像依據一般線性模型並整合年齡因子與性別因子,估算出年齡因子的權重係數與性別因子的權重係數,並取得平滑化的每一第二灰質影像對應的殘差影像 Step 330‧‧‧ The smoothed second gray matter images are based on a general linear model and are integrated with the age factor and the gender factor to estimate the weight factor of the age factor and the weight factor of the gender factor, and obtain each smoothed second Residual image corresponding to gray matter image

步驟340‧‧‧將每一殘差影像的第i個體素的體積大小進行標準差計算,以取得估計殘差影像的第i個體素之體積標準差,進而建立資料庫 Step 340‧‧‧ calculates the standard deviation of the volume size of the ith voxel of each residual image to obtain the estimated volume standard deviation of the ith voxel of the residual image, and then establish a database

第1圖為本發明個人化大腦影像評估之方法的一實施例方法流程圖。 第2圖為第1圖之步驟130所述的疊合程序之一實施例方法流程圖。 第3圖為第1圖之步驟140所述的資料庫的建立方法之一實施例方法流程圖。 第4圖為本發明個人化大腦影像評估之非暫時性電腦可讀媒體之一實施例方塊圖。 第5圖為本發明個人化大腦影像評估之設備之一實施例方塊圖。 第6圖為不同樣本數所建立的資料庫對同一病人進行個人化大腦影像評估之方法後輸出的評估報告之第一灰質影像。 第7A圖與第7B圖分別為利用第1圖的個人化大腦影像評估之方法與傳統T檢定統計方式對同一病人的海馬體進行萎縮電腦模擬評估之示意圖。FIG. 1 is a method flowchart of an embodiment of a method for personalizing brain image assessment according to the present invention. FIG. 2 is a flowchart of a method according to an embodiment of the overlapping procedure described in step 130 of FIG. 1. FIG. 3 is a method flowchart of an embodiment of the method for establishing a database described in step 140 of FIG. 1. FIG. 4 is a block diagram of an embodiment of a non-transitory computer-readable medium for personalized brain imaging assessment according to the present invention. FIG. 5 is a block diagram of an embodiment of a device for personalizing brain image evaluation according to the present invention. Figure 6 is the first gray matter image of the evaluation report output from the database created by different sample numbers for the same patient with a personalized brain image evaluation method. Figures 7A and 7B are schematic diagrams of computerized atrophy assessment of the hippocampus of the same patient using the personalized brain image evaluation method of Figure 1 and the traditional T test statistical method.

Claims (7)

一種個人化大腦影像評估之方法,該方法包含下列步驟: 接收一病人的大腦之一第一T1加權影像(T1 weighted image); 對該第一T1加權影像進行一組織分割程序,以取得一第一灰質影像; 將該第一灰質影像配准至經一疊合程序所建立的一標準化正常灰質影像模板; 將配准至該標準化正常灰質影像模板的該第一灰質影像依據一體素分析公式取得其具有的每一體素對應的一評分值,其中,該體素分析公式為 ,TW為第i個該體素對應的該評分值,x*為配准至該標準化正常灰質影像模板的該第一灰質影像之第i個該體素的體積大小, 為依據該病人的年齡與性別於一資料庫取得一預期影像之第i個體素的體積大小,s為依據該資料庫取得一估計殘差影像之第i個體素的體積標準差, 為T檢定的一比例因子,n為建立該資料庫的樣本數,i為大於或等於1的正整數; 判斷每一該體素對應的該評分值是否大於或小於一預期值;以及 依據大於或小於該預期值的該些體素回推其所屬的腦區,並輸出一評估報告。 A method for personalizing brain image evaluation, the method includes the following steps: receiving a first T1 weighted image of a patient's brain; and performing a tissue segmentation procedure on the first T1 weighted image to obtain a first A gray matter image; registering the first gray matter image to a standardized normal gray matter image template established by a superimposition procedure; obtaining the first gray matter image registered to the standardized normal gray matter image template according to an integrin analysis formula It has a score value corresponding to each voxel, wherein the voxel analysis formula is , TW is the scoring value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized normal gray matter image template, In order to obtain the volume of the i-th voxel of an expected image in a database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel in an estimated residual image obtained from the database, Is a scale factor for the T test, n is the number of samples to build the database, i is a positive integer greater than or equal to 1; determine whether the score value corresponding to each voxel is greater than or less than an expected value; and based on greater than Or the voxels smaller than the expected value are pushed back to the brain region to which they belong, and an evaluation report is output. 根據申請專利範圍第1項之個人化大腦影像評估之方法,其中,該疊合程序包含以下步驟: 取得n個健康受試者的大腦之一第二T1加權影像,n為大於或等於三十之正整數; 對每一該第二T1加權影像分別進行該組織分割程序,以取得對應的一第二灰質影像; 將該些第二灰質影像分別配准至MNI(Montreal Neurological Institute)標準空間模板;以及 將配准至該MNI標準空間模板的該些第二灰質影像進行疊合,並以影像亮度取平均值的方式建立該標準化正常灰質影像模板。The method for personalizing brain image evaluation according to item 1 of the patent application scope, wherein the superimposing procedure includes the following steps: obtaining a second T1-weighted image of the brain of one of n healthy subjects, where n is greater than or equal to thirty A positive integer; performing the tissue segmentation procedure on each of the second T1 weighted images to obtain a corresponding second gray matter image; registering the second gray matter images to the MNI (Montreal Neurological Institute) standard space template ; And superimposing the second gray matter images registered to the MNI standard space template, and establishing the normalized normal gray matter image template by averaging the image brightness. 根據申請專利範圍第2項之個人化大腦影像評估之方法,其中,該資料庫的建立方法包含以下步驟: 將每一該第二灰質影像配准至該標準化正常灰質影像模板; 對配准至該正常影像模板的每一該第二灰質影像進行一平滑化處理; 將平滑化的該些第二灰質影像依據一般線性模型並整合年齡因子與性別因子,估算出該年齡因子的權重係數與該性別因子的權重係數,並取得平滑化的每一該第二灰質影像對應的一殘差影像;以及 將每一該殘差影像的第i個體素的體積大小進行標準差計算,以取得該估計殘差影像的第i個體素之體積標準差,進而建立該資料庫。The method for personalizing brain image evaluation according to item 2 of the scope of patent application, wherein the method of establishing the database includes the following steps: registering each of the second gray matter images to the standardized normal gray matter image template; Each of the second gray matter images of the normal image template is subjected to a smoothing process; the smoothed second gray matter images are based on a general linear model and an age factor and a gender factor are integrated to estimate a weight coefficient of the age factor and the Weight coefficient of the gender factor, and obtain a residual image corresponding to each smoothed second gray matter image; and calculate the standard deviation of the volume size of the i-th voxel of each residual image to obtain the estimate The volume standard deviation of the i-th voxel in the residual image is used to establish the database. 一種個人化大腦影像評估之非暫時性電腦可讀媒體,其經組態以儲存若干操作指令,該等操作指令在由一或多個處理器執行時使得該一或多個處理器執行以下操作: 接收一病人的大腦之一第一T1加權影像; 對該第一T1加權影像進行一組織分割程序,以取得一第一灰質影像; 將該第一灰質影像配准至經一疊合程序所建立的一標準化正常灰質影像模板; 將配准至該標準化正常灰質影像模板的該第一灰質影像依據一體素分析公式取得其具有的每一體素對應的一評分值,其中,該體素分析公式為 ,TW為第i個該體素對應的該評分值,x*為配准至該標準化正常灰質影像模板的該第一灰質影像之第i個該體素的體積大小, 為依據該病人的年齡與性別於一資料庫取得一預期影像之第i個體素的體積大小,s為依據該資料庫取得一估計殘差影像之第i個體素的體積標準差, 為T檢定的一比例因子,n為建立該資料庫的樣本數,i為大於或等於1的正整數; 判斷每一該體素對應的評分值是否大於或小於一預期值;以及 依據大於或小於該預期值的該些體素回推其所屬的腦區,並輸出一評估報告。 A non-transitory computer-readable medium for personalized brain imaging assessment, configured to store a number of operating instructions that, when executed by one or more processors, cause the one or more processors to perform the following operations Receiving a first T1 weighted image of a patient's brain; performing a tissue segmentation procedure on the first T1 weighted image to obtain a first gray matter image; registering the first gray matter image to a superimposed procedure A normalized normal gray matter image template established; the first gray matter image registered to the normalized normal gray matter image template is used to obtain a score value corresponding to each voxel it has according to the integrin analysis formula, wherein the voxel analysis formula for , TW is the scoring value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized normal gray matter image template, In order to obtain the volume of the i-th voxel of an expected image in a database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel in an estimated residual image obtained from the database, Is a scale factor for the T test, n is the number of samples to build the database, i is a positive integer greater than or equal to 1; determine whether the score value corresponding to each voxel is greater than or less than an expected value; and based on greater than or Those voxels smaller than the expected value push back the brain region to which they belong and output an evaluation report. 根據申請專利範圍第4項之個人化大腦影像評估之非暫時性電腦可讀媒體,其中,該疊合程序包含以下操作: 取得n個健康受試者的大腦之一第二T1加權影像,n為大於或等於三十之正整數; 對每一該第二T1加權影像分別進行該組織分割程序,以取得對應的一第二灰質影像; 將該些第二灰質影像分別配准至一MNI標準空間模板;以及 將配准至該MNI標準空間模板的該些第二灰質影像進行疊合,並以影像亮度取平均值的方式建立該標準化正常灰質影像模板。A non-transitory computer-readable medium for personalized brain imaging assessment according to item 4 of the scope of the patent application, wherein the overlay procedure includes the following operations: obtaining a second T1-weighted image of one of the brains of n healthy subjects, n Is a positive integer greater than or equal to thirty; the tissue segmentation procedure is performed on each of the second T1 weighted images to obtain a corresponding second gray matter image; the second gray matter images are registered to an MNI standard respectively A spatial template; and superimposing the second gray matter images registered to the MNI standard spatial template, and establishing the normalized normal gray matter image template by averaging the image brightness. 根據申請專利範圍第5項之個人化大腦影像評估之非暫時性電腦可讀媒體,其中,該資料庫的建立方法包含以下操作: 將每一該第二灰質影像配准至該標準化正常灰質影像模板; 對配准至該標準化正常灰質影像模板的每一該第二灰質影像進行一平滑化處理; 將平滑化的該些第二灰質影像依據一般線性模型並整合年齡因子與性別因子,估算出該年齡因子的權重係數與該性別因子的權重係數,並取得平滑化的每一該第二灰質影像對應的一殘差影像;以及 將每一該殘差影像的第i個體素的體積大小進行標準差計算,以取得該估計殘差影像的第i個體素之體積標準差,進而建立該資料庫。The non-transitory computer-readable medium for personalized brain image assessment according to item 5 of the scope of patent application, wherein the method of establishing the database includes the following operations: registering each of the second gray matter images to the standardized normal gray matter image A template; performing a smoothing process on each of the second gray matter images registered to the standardized normal gray matter image template; estimating the smoothed second gray matter images according to a general linear model and integrating an age factor and a gender factor to estimate A weighting coefficient of the age factor and a weighting coefficient of the gender factor, and obtaining a residual image corresponding to each smoothed second gray matter image; and performing the size of the i-th voxel of each residual image Standard deviation calculation to obtain the volume standard deviation of the i-th voxel of the estimated residual image, and then establish the database. 一種個人化大腦影像評估之設備,該設備包含: 一或多個處理器; 一儲存單元;以及 至少一程式,其中該至少一程式儲存於該儲存單元中且經組態以由該一或多個處理器執行,該至少一程式整體上包含用於以下操作指令: 接收一病人的大腦之一第一T1加權影像; 對該第一T1加權影像進行一組織分割程序,以取得一第一灰質影像; 將該第一灰質影像配准至經一疊合程序所建立的一標準化正常灰質影像模板; 將配准至該標準化正常灰質影像模板的該第一灰質影像依據一體素分析公式取得其具有的每一體素對應的一評分值,其中,該體素分析公式為 ,TW為第i個該體素對應的該評分值,x*為配准至該標準化正常灰質影像模板的該第一灰質影像之第i個該體素的體積大小, 為依據該病人的年齡與性別於一資料庫取得一預期影像之第i個體素的體積大小,s為依據該資料庫取得一估計殘差影像之第i個體素的體積標準差, 為T檢定的一比例因子,n為建立該資料庫的樣本數,i為大於或等於1的正整數; 判斷每一該體素對應的該評分值是否大於或小於一預期值;以及 依據大於或小於該預期值的該些體素回推其所屬的腦區,並輸出一評估報告。 A device for personalizing brain imaging assessment, the device comprising: one or more processors; a storage unit; and at least one program, wherein the at least one program is stored in the storage unit and configured to be processed by the one or more The processors execute the at least one program as a whole including the following operation instructions: receiving a first T1 weighted image of a patient's brain; and performing a tissue segmentation procedure on the first T1 weighted image to obtain a first gray matter Image; registering the first gray matter image to a normalized normal gray matter image template created by a superimposition process; registering the first gray matter image to the normalized normal gray matter image template according to an integrin analysis formula A score corresponding to each voxel of, wherein the voxel analysis formula is , TW is the scoring value corresponding to the i-th voxel, and x * is the volume size of the i-th voxel of the first gray matter image registered to the normalized normal gray matter image template, In order to obtain the volume of the i-th voxel of an expected image in a database based on the age and sex of the patient, s is the volume standard deviation of the i-th voxel in an estimated residual image obtained from the database, Is a scale factor for the T test, n is the number of samples to build the database, i is a positive integer greater than or equal to 1; determine whether the score value corresponding to each voxel is greater than or less than an expected value; and based on greater than Or the voxels smaller than the expected value are pushed back to the brain region to which they belong, and an evaluation report is output.
TW108101383A 2019-01-14 2019-01-14 Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging TWI679652B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW108101383A TWI679652B (en) 2019-01-14 2019-01-14 Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW108101383A TWI679652B (en) 2019-01-14 2019-01-14 Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging

Publications (2)

Publication Number Publication Date
TWI679652B true TWI679652B (en) 2019-12-11
TW202027091A TW202027091A (en) 2020-07-16

Family

ID=69582400

Family Applications (1)

Application Number Title Priority Date Filing Date
TW108101383A TWI679652B (en) 2019-01-14 2019-01-14 Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging

Country Status (1)

Country Link
TW (1) TWI679652B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113393141A (en) * 2021-06-23 2021-09-14 科大讯飞股份有限公司 Reading capability evaluation method, device and equipment
CN113892936A (en) * 2021-09-24 2022-01-07 天津大学 Interpretable brain age prediction method based on full convolution neural network

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI788109B (en) * 2021-11-18 2022-12-21 國立陽明交通大學 System for integrating medical images based on spatial position of radiography and method thereof

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102088910A (en) * 2008-07-07 2011-06-08 浜松光子学株式会社 Brain disease diagnosis system
CN105046709A (en) * 2015-07-14 2015-11-11 华南理工大学 Nuclear magnetic resonance imaging based brain age analysis method
CN106156484A (en) * 2016-06-08 2016-11-23 中国科学院自动化研究所 Disease of brain individuation Forecasting Methodology based on nuclear magnetic resonance image and system
TWI614627B (en) * 2017-10-06 2018-02-11 林慶波 Method, non-transitory computer-readable media and apparatus for nerve tracing
TWI637187B (en) * 2017-09-19 2018-10-01 臺北榮民總醫院 Magnetic resonance image analysis method and method for assessing the risk of radiation therapy
CN109147941A (en) * 2018-10-17 2019-01-04 上海交通大学 Brain robustness appraisal procedure based on structure nuclear magnetic resonance image data

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102088910A (en) * 2008-07-07 2011-06-08 浜松光子学株式会社 Brain disease diagnosis system
US8676294B2 (en) * 2008-07-07 2014-03-18 Hamamatsu Photonics K.K. Brain disease diagnosis system
CN105046709A (en) * 2015-07-14 2015-11-11 华南理工大学 Nuclear magnetic resonance imaging based brain age analysis method
CN106156484A (en) * 2016-06-08 2016-11-23 中国科学院自动化研究所 Disease of brain individuation Forecasting Methodology based on nuclear magnetic resonance image and system
TWI637187B (en) * 2017-09-19 2018-10-01 臺北榮民總醫院 Magnetic resonance image analysis method and method for assessing the risk of radiation therapy
TWI614627B (en) * 2017-10-06 2018-02-11 林慶波 Method, non-transitory computer-readable media and apparatus for nerve tracing
CN109147941A (en) * 2018-10-17 2019-01-04 上海交通大学 Brain robustness appraisal procedure based on structure nuclear magnetic resonance image data

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113393141A (en) * 2021-06-23 2021-09-14 科大讯飞股份有限公司 Reading capability evaluation method, device and equipment
CN113393141B (en) * 2021-06-23 2023-12-01 科大讯飞股份有限公司 Method, device and equipment for evaluating reading capability
CN113892936A (en) * 2021-09-24 2022-01-07 天津大学 Interpretable brain age prediction method based on full convolution neural network
CN113892936B (en) * 2021-09-24 2023-08-18 天津大学 Interpretable brain age prediction method based on full convolution neural network

Also Published As

Publication number Publication date
TW202027091A (en) 2020-07-16

Similar Documents

Publication Publication Date Title
Das et al. Registration based cortical thickness measurement
Iscan et al. Test–retest reliability of freesurfer measurements within and between sites: Effects of visual approval process
Tucholka et al. An empirical comparison of surface-based and volume-based group studies in neuroimaging
Chang et al. Performance measure characterization for evaluating neuroimage segmentation algorithms
Smith et al. Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data
Xue et al. Automatic segmentation and reconstruction of the cortex from neonatal MRI
TWI679652B (en) Method, non-transitory computer-readable media and apparatus for evaluating personalized brain imaging
Tang et al. RABBIT: rapid alignment of brains by building intermediate templates
JP6036009B2 (en) Medical image processing apparatus and program
WO2013047278A1 (en) Medical image processing device, medical image processing method, program
Li et al. Cortical thickness estimation in longitudinal stroke studies: a comparison of 3 measurement methods
Brun et al. Mapping brain abnormalities in boys with autism
Gedamu et al. Automated quality control of brain MR images
TWI624807B (en) Iterative analysis of medical images
Shiee et al. Reconstruction of the human cerebral cortex robust to white matter lesions: method and validation
US10729389B2 (en) 3D assessment of conjugant eye deviation for the identificaiton of acute ischemic stroke
Hachaj et al. CAD system for automatic analysis of CT perfusion maps
Van Hecke et al. On the construction of a ground truth framework for evaluating voxel-based diffusion tensor MRI analysis methods
CN108701360B (en) Image processing system and method
Rostowsky et al. Macroscale white matter alterations due to traumatic cerebral microhemorrhages are revealed by diffusion tensor imaging
US9883817B2 (en) Management, assessment and treatment planning for inflammatory bowel disease
Damangir et al. Reproducible segmentation of white matter hyperintensities using a new statistical definition
US20210073965A1 (en) Devices, systems, and methods for anchor-point-enabled multi-scale subfield alignment
Styner et al. Automatic brain segmentation in rhesus monkeys
Valsasina et al. Improved Assessment of Longitudinal Spinal Cord Atrophy in Multiple Sclerosis Using a Registration‐Based Approach: Relevance for Clinical Studies