JP2011062253A - Tomographic image processor, muscle cross-sectional area measurement system, and tomographic image processing method - Google Patents
Tomographic image processor, muscle cross-sectional area measurement system, and tomographic image processing method Download PDFInfo
- Publication number
- JP2011062253A JP2011062253A JP2009213565A JP2009213565A JP2011062253A JP 2011062253 A JP2011062253 A JP 2011062253A JP 2009213565 A JP2009213565 A JP 2009213565A JP 2009213565 A JP2009213565 A JP 2009213565A JP 2011062253 A JP2011062253 A JP 2011062253A
- Authority
- JP
- Japan
- Prior art keywords
- tomographic image
- muscle
- subject
- area
- psoas
- 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.)
- Pending
Links
- 210000003205 muscle Anatomy 0.000 title claims abstract description 163
- 238000005259 measurement Methods 0.000 title claims abstract description 36
- 238000003672 processing method Methods 0.000 title claims abstract description 10
- 238000003384 imaging method Methods 0.000 claims abstract description 120
- 210000002097 psoas muscle Anatomy 0.000 claims abstract description 92
- 238000004364 calculation method Methods 0.000 claims abstract description 71
- 238000012545 processing Methods 0.000 claims description 72
- 230000005484 gravity Effects 0.000 claims description 29
- 238000012937 correction Methods 0.000 abstract description 8
- 238000000034 method Methods 0.000 description 25
- 238000002595 magnetic resonance imaging Methods 0.000 description 14
- 230000003187 abdominal effect Effects 0.000 description 10
- 238000004590 computer program Methods 0.000 description 10
- 230000003387 muscular Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 210000001061 forehead Anatomy 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 206010073306 Exposure to radiation Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 235000001916 dieting Nutrition 0.000 description 1
- 230000037228 dieting effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 210000001624 hip Anatomy 0.000 description 1
- 210000004394 hip joint Anatomy 0.000 description 1
- 230000000968 intestinal effect Effects 0.000 description 1
- 210000004705 lumbosacral region Anatomy 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 210000001087 myotubule Anatomy 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 210000002027 skeletal muscle Anatomy 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Abstract
Description
本発明は、断層画像処理装置、筋断面積測定システム及び断層画像処理方法に関する。 The present invention relates to a tomographic image processing apparatus, a muscle cross-sectional area measuring system, and a tomographic image processing method.
近年の高齢化社会が到来した環境下では、一定の筋力を保持することが自立した生活を行うために重要であることが指摘されている。高齢者が一定の筋力を保持するためには、定期的に筋力を測定し、筋力水準を把握した上でその水準に適した方法・内容で筋力の維持又は向上を図ることが重要である。筋力測定方法として、筋の最大筋力を発揮させてその力を測定する方法が従来から用いられているが、この筋力測定方法は心循環器系や骨格筋系に負荷がかかるため、高齢者や疾患を有する人のように最大筋力を発揮することで負傷する可能性の高い人は筋力測定を行うことができない。したがって、筋力の維持又は向上を達成できない場合がある。 It has been pointed out that maintaining a certain muscular strength is important for living independently in an environment where an aging society has come in recent years. In order for an elderly person to maintain a certain muscular strength, it is important to regularly measure the muscular strength, grasp the muscular strength level, and maintain or improve the muscular strength by a method / content suitable for the level. As a method for measuring muscle strength, a method of measuring the strength of a muscle by exerting the maximum muscle strength of the muscle has been used conventionally, but this muscle strength measurement method places a burden on the cardiovascular system and the skeletal muscle system. A person who is likely to be injured by exerting maximum muscle strength, such as a person with a disease, cannot measure muscle strength. Therefore, the maintenance or improvement of muscle strength may not be achieved.
これに対して、近年では、筋（筋肉）が発揮可能な最大の力とその筋の横断面積（筋断面積）とがおおよそ比例関係にあることを利用して、磁気共鳴画像法（Magnetic Resonance Imaging：ＭＲＩ）やコンピュータ断像撮影法（Computed Tomography：ＣＴ）により被検者の体内の断層画像を撮像し、この断層画像から筋断面積を算出することで、その被検者の筋力水準のポテンシャルを推定する方法が検討されている。この筋力の推定方法は、ＭＲＩやＣＴによる断層画像の撮像のみで筋力を推定するため、従来の筋力測定方法のように被検者に負荷をかけることがない。 On the other hand, in recent years, magnetic resonance imaging (Magnetic Resonance) is utilized by utilizing the fact that the maximum force that a muscle (muscle) can exert and the cross-sectional area (muscle cross-sectional area) of the muscle are approximately proportional. By taking a tomographic image of the subject's body by imaging (MRI) or computerized tomography (CT) and calculating the muscle cross-sectional area from the tomographic image, the muscle strength level of the subject can be calculated. Methods for estimating the potential are being studied. Since this muscle strength estimation method estimates muscle strength only by taking a tomographic image by MRI or CT, it does not place a load on the subject unlike conventional muscle strength measurement methods.
この筋力の推定方法で用いられる筋断面積とは、測定対象となる筋の最も太い部分を筋の走行方向に対して垂直な方向に切断した際の断面積を指す。そして、この筋断面積を正確に算出する場合には、複数の方向から被検者の断層画像を撮像することで被検者の筋の走行方向を正確に求めた上で、その筋の最も太い位置で走行方向に対して垂直な断層画像を撮像する必要がある。しかし、複数の方向から被検者の断層画像を撮像する方法では、筋断面積の算出に必要な撮像データが多いために被検者の撮像時間が長時間化し、被検者の負担も大きくなる可能性がある。また、一被検者あたりに必要な撮像時間が長くなるため、一台の撮像装置で多人数の撮像を行うことが困難になる。 The muscle cross-sectional area used in this muscle strength estimation method refers to a cross-sectional area when the thickest portion of the muscle to be measured is cut in a direction perpendicular to the running direction of the muscle. When accurately calculating the muscle cross-sectional area, the tomographic images of the subject are taken from a plurality of directions to accurately determine the traveling direction of the subject's muscle, It is necessary to capture a tomographic image perpendicular to the traveling direction at a thick position. However, in the method of capturing a tomographic image of a subject from multiple directions, the amount of imaging data required for calculation of the muscle cross-sectional area is large, so the imaging time of the subject is prolonged and the burden on the subject is also large. There is a possibility. Moreover, since the imaging time required per subject becomes long, it becomes difficult to perform imaging of a large number of people with one imaging device.
これに対して、一般的に水平なベッドに横たわった被検者に対して垂直な断層画像を撮像する従来のＭＲＩやＣＴを用いて筋断面積を概算する方法が用いられることが多い。この場合、従来のＭＲＩやＣＴによる撮像によって得られる断層画像を利用するため、撮像時間は短縮される。ただし、従来のＭＲＩやＣＴによる断層画像は、測定対象となる筋の走行方向に対して垂直な方向の筋断面とは一致しない場合があるため、例えば非特許文献１では、死体を用いて筋の走行方向の平均値を算出し、この平均値を用いてＭＲＩによる被検者の断層画像における筋断面積を補正することで実際の筋断面積との誤差を小さくする方法が示されている。 On the other hand, generally, a method of roughly estimating a muscle cross-sectional area using conventional MRI or CT that captures a vertical tomographic image for a subject lying on a horizontal bed is often used. In this case, since a tomographic image obtained by conventional MRI or CT imaging is used, the imaging time is shortened. However, since a conventional tomographic image by MRI or CT may not coincide with a muscle cross-section in a direction perpendicular to the running direction of the muscle to be measured, in Non-Patent Document 1, for example, A method is shown in which an average value in the traveling direction is calculated, and an error with the actual muscle cross-sectional area is reduced by correcting the muscle cross-sectional area in the tomographic image of the subject by MRI using the average value. .
しかしながら、非特許文献１記載の方法では、実際の被検者では一人ひとりで体格が異なるため、被検者の断層画像から算出した筋断面積と実際の筋断面積との誤差を十分に補正できない場合がある。また、同じ被検者でも撮像時の体勢の違いによって誤差が大きくなる場合がある。したがって、トレーニングやダイエットの実施、サプリメントの摂取等による筋断面積の変化の測定のように微小な変化を検出したい場合に、高い精度で測定することが困難であった。 However, in the method described in Non-Patent Document 1, each physical subject has a different physique, so the error between the muscle cross-sectional area calculated from the tomographic image of the subject and the actual muscle cross-sectional area cannot be sufficiently corrected. There is a case. Even in the same subject, an error may increase due to a difference in posture during imaging. Therefore, it is difficult to measure with high accuracy when it is desired to detect minute changes such as measurement of changes in the muscle cross-sectional area due to training, dieting, supplement intake, and the like.
本発明は上記を鑑みてなされたものであり、ＭＲＩやＣＴによる被検者の断層画像からの筋断面積の算出をより高精度で行うことができる断層画像処理装置、筋断面積測定システム及び断層画像処理方法を提供することを目的とする。 The present invention has been made in view of the above, and a tomographic image processing apparatus, a muscle cross-sectional area measuring system, and a muscle cross-sectional area measuring system capable of calculating a muscle cross-sectional area from a tomographic image of a subject by MRI or CT with higher accuracy. An object is to provide a tomographic image processing method.
上記目的を達成するため、本発明に係る断層画像処理装置は、被検者の大腰筋が撮像された複数の断層画像を取得する画像取得手段と、断層画像に含まれている被検者の大腰筋に対応する部分（大腰筋領域）の面積を断層画像毎に算出する面積算出手段と、断層画像に含まれている被検者の大腰筋領域の重心位置を断層画像毎に算出する重心算出手段と、重心算出手段により算出された重心位置に基づいて、大腰筋の走行方向を算出する走行方向算出手段と、走行方向算出手段により算出された走行方向と面積算出手段により算出された断層画像毎の大腰筋領域の面積とに基づいて、走行方向に対して垂直な大腰筋の筋断面に対応する部分の面積を算出する筋断面積算出手段と、を備えることを特徴とする。 In order to achieve the above object, a tomographic image processing apparatus according to the present invention includes an image acquisition means for acquiring a plurality of tomographic images obtained by imaging the psoas major muscle of a subject, and a subject included in the tomographic image. Area calculation means for calculating for each tomographic image the area of the portion corresponding to the psoas major muscle (large psoas muscle region), and for each tomographic image the center of gravity position of the subject's psoas muscle region included in the tomographic image Based on the center of gravity calculated by the center of gravity calculating means, the traveling direction calculating means for calculating the traveling direction of the psoas major, and the traveling direction and area calculating means calculated by the traveling direction calculating means. And a muscle cross-sectional area calculating means for calculating an area of a portion corresponding to the muscle cross section of the psoas major perpendicular to the running direction based on the area of the psoas muscle region calculated for each tomographic image. It is characterized by that.
上記の断層画像処理装置によれば、被検者の複数の断層画像の大腰筋領域の面積及び重心位置が画像毎に算出され、この重心位置に基づいて被検者の大腰筋の走行方向が算出される。そして、大腰筋の走行方向と複数の断層画像毎の被検者の大腰筋領域の面積とに基づいて、被検者の大腰筋の筋断面に対応する部分の面積が算出される。このように、被検者の複数の断層画像に基づいてその被検者の大腰筋の走行方向を求め、求められた走行方向に基づいて筋断面積を算出するため、被検者の体格や測定時の体勢に応じた補正を行うことができ、より高い精度で筋断面積の算出を行うことができる。また、上記の断層画像処理装置による断層画像処理は、トレーニングやリハビリによる効果の評価対象として好適な大腰筋に対して適用することが好ましい。すなわち、上記の断層画像処理装置によって被検者の大腰筋の筋断面積の算出を行うことによって、被検者の大腰筋の微小な変化を適切に計測することができ、被検者の筋肉の状態をより適切に把握することができる。 According to the tomographic image processing apparatus described above, the area of the psoas muscle region and the center of gravity position of a plurality of tomographic images of the subject are calculated for each image, and the subject's psoas muscle travels based on the center of gravity position. A direction is calculated. Then, based on the running direction of the psoas muscle and the area of the psoas muscle region of the subject for each of a plurality of tomographic images, the area of the portion corresponding to the muscle cross section of the psoas muscle of the subject is calculated. . Thus, in order to determine the traveling direction of the psoas major of the subject based on the tomographic images of the subject and to calculate the muscle cross-sectional area based on the determined traveling direction, the physique of the subject And correction according to the posture at the time of measurement, and the muscle cross-sectional area can be calculated with higher accuracy. Moreover, it is preferable that the tomographic image processing by the above-described tomographic image processing apparatus is applied to the psoas major muscle that is suitable as an evaluation target of the effects of training and rehabilitation. That is, by calculating the muscle cross-sectional area of the large psoas muscle of the subject by the tomographic image processing device, it is possible to appropriately measure a minute change in the psoas major muscle of the subject. Can better grasp the state of muscles.
ここで、上記作用を効果的に奏する構成として、具体的には、複数の断層画像は、互いに平行な大腰筋の断面が撮像されたものであり、走行方向算出手段は、複数の断層画像における被検者の大腰筋の重心位置と当該断層画像の撮像位置とを用いて走行方向を算出する態様が挙げられる。 Here, as a configuration that effectively exhibits the above-described operation, specifically, the plurality of tomographic images are obtained by imaging cross sections of the psoas muscle parallel to each other, and the traveling direction calculation means includes a plurality of tomographic images. In which the running direction is calculated using the center of gravity position of the psoas muscle of the subject and the imaging position of the tomographic image.
また、本発明に係る筋断面積測定システムは、撮像装置と断層画像処理装置とを含んで構成される筋断面積測定システムであって、撮像装置は、被検者の大腰筋の複数の断層画像を撮像する撮像手段を備え、断層画像処理装置は、撮像装置で撮像された複数の断層画像を取得する画像取得手段と、断層画像に含まれている被検者の大腰筋領域の面積を断層画像毎に算出する面積算出手段と、断層画像に含まれている被検者の大腰筋領域の重心位置を断層画像毎に算出する重心算出手段と、重心算出手段により算出された重心位置に基づいて、大腰筋の走行方向を算出する走行方向算出手段と、走行方向算出手段により算出された走行方向と面積算出手段により算出された断層画像毎の大腰筋領域の面積とに基づいて、走行方向に対して垂直な大腰筋の筋断面に対応する部分の面積を算出する筋断面積算出手段と、を備えることを特徴とする。 The muscle cross-sectional area measurement system according to the present invention is a muscle cross-sectional area measurement system including an imaging device and a tomographic image processing device, and the imaging device includes a plurality of psoas major muscles of a subject. The tomographic image processing apparatus includes an image capturing unit configured to capture a tomographic image, and the tomographic image processing apparatus includes: an image acquiring unit configured to acquire a plurality of tomographic images captured by the imaging apparatus; and a subject's psoas muscle region included in the tomographic image. Calculated by area calculation means for calculating the area for each tomographic image, centroid calculation means for calculating the position of the center of gravity of the subject's psoas muscle included in the tomographic image for each tomographic image, and centroid calculation means Based on the center-of-gravity position, a travel direction calculation unit that calculates the travel direction of the psoas major muscle, a travel direction calculated by the travel direction calculation unit, and an area of the psoas muscle region for each tomographic image calculated by the area calculation unit Perpendicular to the direction of travel And Sujidan area calculation means for calculating the area of a portion corresponding to the muscle cross-section of the psoas muscle, characterized in that it comprises a.
なお、本発明は、上記のように断層画像処理装置の発明として記述できる他に、以下のように断層画像処理方法の発明としても記述することができる。これはカテゴリが異なるだけで、実質的に同一の発明であり、同様の作用及び効果を奏する。 In addition to being described as an invention of a tomographic image processing apparatus as described above, the present invention can also be described as an invention of a tomographic image processing method as follows. This is substantially the same invention only in different categories, and has the same operations and effects.
すなわち、本発明に係る断層画像処理方法は、断層画像処理装置による断層画像処理方法であって、断層画像処理装置が、被検者の大腰筋が撮像された複数の断層画像を取得する画像取得ステップと、断層画像処理装置が、断層画像に含まれている被検者の大腰筋領域の面積を断層画像毎に算出する面積算出ステップと、断層画像処理装置が、断層画像に含まれている被検者の大腰筋領域の重心位置を断層画像毎に算出する重心算出ステップと、断層画像処理装置が、重心算出ステップにおいて算出された重心位置に基づいて、大腰筋の走行方向を算出する走行方向算出ステップと、断層画像処理装置が、走行方向算出ステップにおいて算出された走行方向と面積算出ステップにおいて算出された断層画像毎の大腰筋領域の面積とに基づいて、走行方向に対して垂直な大腰筋の筋断面に対応する部分の面積を算出する筋断面積算出ステップと、を備えることを特徴とする。 That is, the tomographic image processing method according to the present invention is a tomographic image processing method by a tomographic image processing apparatus, in which the tomographic image processing apparatus acquires a plurality of tomographic images obtained by imaging the subject's psoas muscle. The tomographic image includes an acquisition step, an area calculation step in which the tomographic image processing device calculates the area of the psoas major region of the subject included in the tomographic image for each tomographic image, and a tomographic image processing device. A center-of-gravity calculation step for calculating the position of the center of gravity of the psoas muscle region of the subject who is tomographic image, and the tomographic image processing device based on the position of the center of gravity calculated in the center-of-gravity calculation step Based on the travel direction calculated in the travel direction calculation step and the area of the psoas muscle region for each tomographic image calculated in the area calculation step. It characterized in that it and a Sujidan area calculation step of calculating the area of a portion corresponding to the muscle cross-section of the psoas major perpendicular to the direction of travel.
本発明によれば、ＭＲＩやＣＴによる撮像データからの筋断面積の算出をより高精度で行うことができる断層画像処理装置、筋断面積測定システム及び断層画像処理方法が提供される。 According to the present invention, a tomographic image processing apparatus, a muscle cross-sectional area measuring system, and a tomographic image processing method capable of calculating a muscle cross-sectional area from imaging data by MRI or CT with higher accuracy are provided.
以下、添付図面を参照して、本発明を実施するための最良の形態を詳細に説明する。なお、図面の説明においては同一要素には同一符号を付し、重複する説明を省略する。 The best mode for carrying out the present invention will be described below in detail with reference to the accompanying drawings. In the description of the drawings, the same elements are denoted by the same reference numerals, and redundant description is omitted.
本実施形態に係る筋断面積測定システム１は、撮像装置１０と、断層画像処理装置２０とを備える。撮像装置１０は、大腰筋を含む被検者の断層画像を撮像する装置であって、撮像装置１０によって撮像された断層画像が断層画像処理装置２０に送られる。断層画像処理装置２０は、被検者の大腰筋が含まれる断層画像に基づいて被検者の大腰筋の筋断面積を算出する装置である。 The muscle cross-sectional area measurement system 1 according to this embodiment includes an imaging device 10 and a tomographic image processing device 20. The imaging device 10 is a device that captures a tomographic image of a subject including the psoas muscle, and the tomographic image captured by the imaging device 10 is sent to the tomographic image processing device 20. The tomographic image processing device 20 is a device that calculates a muscle cross-sectional area of the subject's psoas muscle based on a tomographic image including the subject's psoas muscle.
撮像装置１０は、被検者の断面を撮像する装置であり、撮像のための光線やＸ線等の電磁波を被検者に対して照射する照射装置と、被検者からの電磁波を受信する検出装置とを含んで構成される。撮像装置１０としては、具体的には、ＭＲＩ装置やＸ線ＣＴ装置等が用いられる。ＭＲＩ装置は、具体的には、被検者に高周波（ＲＦ波）を照射し、被検者中の原子核に生じる核磁気共鳴現象を観察することで体内組織の解剖学的情報、及び血液量などの機能情報を取得する装置である。また、Ｘ線ＣＴ装置は、被験者にＸ線を照射して体内を透過させ、その透過率の違いから体内組織の解剖学的情報を取得する装置である。本実施形態に係る撮像装置１０では、大腰筋の筋断面積を算出するために、被検者の大腰筋を含む断層画像が撮像される。なお、本実施形態の筋断面積測定システム１による筋断面積の測定対象となる大腰筋（psoas major muscle）とは、被検者の腰椎と大腿骨を結ぶ筋肉である腸腰筋の一部であり、股関節や脊柱の屈曲、体幹の安定（姿勢維持）に寄与する筋肉である。このように、大腰筋は被検者の歩行能力や運動能力にも大きく影響を与える筋肉であることから、トレーニングやリハビリによる効果を確認するための指標として大腰筋の筋断面積を測定し、この大腰筋の筋断面積から被検者の筋力水準を推定する方法が用いられる。また、本実施形態における筋断面積とは、測定対象となる筋のうち最も太い部分の筋の走行方向に対して垂直な方向に筋を切断した際の断面積のことをいう。 The imaging device 10 is a device that images a cross section of a subject, and receives an irradiation device that irradiates the subject with electromagnetic waves such as light rays and X-rays for imaging, and electromagnetic waves from the subject. And a detection device. Specifically, an MRI apparatus, an X-ray CT apparatus, or the like is used as the imaging apparatus 10. Specifically, the MRI apparatus irradiates a subject with a high frequency (RF wave) and observes a nuclear magnetic resonance phenomenon generated in a nucleus in the subject, thereby analyzing anatomical information of a body tissue and blood volume. It is a device that acquires functional information such as. The X-ray CT apparatus is an apparatus that irradiates a subject with X-rays and transmits the subject through the body, and acquires anatomical information of the body tissue from the difference in transmittance. In the imaging apparatus 10 according to the present embodiment, a tomographic image including the psoas major muscle of the subject is taken in order to calculate the muscle cross-sectional area of the psoas major. Note that the psoas major muscle, which is a measurement target of the muscle cross-sectional area by the muscle cross-sectional area measurement system 1 of the present embodiment, is one of the intestinal lumbar muscles, which are muscles connecting the lumbar spine and the femur of the subject. It is a muscle that contributes to flexion of the hip joint and spinal column and stability of the trunk (maintaining posture). In this way, the psoas muscle is a muscle that greatly affects the walking ability and exercise ability of the subject, so the muscle cross-sectional area of the psoas muscle is measured as an index to confirm the effects of training and rehabilitation. And the method of estimating a muscular strength level of a subject from the muscle cross-sectional area of this psoas muscle is used. In addition, the muscle cross-sectional area in the present embodiment refers to a cross-sectional area when the muscle is cut in a direction perpendicular to the running direction of the thickest muscle among the muscles to be measured.
断層画像処理装置２０は、撮像装置１０で撮像され、有線のネットワーク等を介して撮像装置１０から送信された被検者の大腰筋を含む断層画像に基づいて、被検者の大腰筋の筋断面積を算出する。断層画像処理装置２０によって算出された被検者の筋断面積はモニタやプリンタ等を介して出力することで被検者等に通知される。断層画像処理装置２０における被検者の筋断面積の算出に係る具体的な方法については後述する。 The tomographic image processing device 20 is imaged by the imaging device 10 and is transmitted from the imaging device 10 via a wired network or the like, based on the tomographic image including the subject's large psoas muscle. Calculate the muscle cross-sectional area. The muscle cross-sectional area of the subject calculated by the tomographic image processing apparatus 20 is output to the subject or the like through a monitor, a printer, or the like. A specific method for calculating the muscle cross-sectional area of the subject in the tomographic image processing apparatus 20 will be described later.
次に、図１及び図２を用いて、撮像装置１０及び断層画像処理装置２０の各装置の構成について説明する。このうち、まず撮像装置１０について説明する。撮像装置１０は、図２に示すように、物理的には、ＣＰＵ（Central Processing Unit）１０ａ、ＲＯＭ（Read Only Memory）１０ｂ、ＲＡＭ（RandomAccess Memory）１０ｃ、通信装置１０ｄ及び記憶装置１０ｅ等を有しており、これらＣＰＵ１０ａ〜記憶装置１０ｅはバス１０ｆに接続されている。ＣＰＵ１０ａは、ＲＯＭ１０ｂ等の内蔵メモリに格納された所定のコンピュータプログラムをＲＡＭ１０ｃにロードして実行することによって、撮像装置１０を統括的に制御する。通信装置１０ｄは、外部の装置（本実施形態では断層画像処理装置２０）と通信を行うためのインターフェースを有する。記憶装置１０ｅは、書き込み／読出しが自在なメモリであって、各種コンピュータプログラム（特に、図７のフローチャートに示す処理のうち撮像装置１０で行われる処理を実行するためのプログラム）を格納する。記憶装置１０ｅは、コンピュータプログラムの実行に必要な各種データ及びコンピュータプログラムの実行によって得られる各種データ（特に、被検者の断層画像の撮像に関するデータ（撮像データ））等を格納するためのデータベースが含まれる（図２では不図示）。 Next, the configuration of each of the imaging apparatus 10 and the tomographic image processing apparatus 20 will be described with reference to FIGS. 1 and 2. Among these, first, the imaging device 10 will be described. As shown in FIG. 2, the imaging device 10 physically includes a CPU (Central Processing Unit) 10a, a ROM (Read Only Memory) 10b, a RAM (Random Access Memory) 10c, a communication device 10d, a storage device 10e, and the like. These CPU 10a to storage device 10e are connected to a bus 10f. The CPU 10a centrally controls the imaging device 10 by loading a predetermined computer program stored in a built-in memory such as the ROM 10b into the RAM 10c and executing it. The communication device 10d has an interface for communicating with an external device (the tomographic image processing device 20 in the present embodiment). The storage device 10e is a readable / writable memory, and stores various computer programs (particularly, a program for executing processing performed by the imaging device 10 among the processing shown in the flowchart of FIG. 7). The storage device 10e has a database for storing various data necessary for the execution of the computer program, various data obtained by the execution of the computer program (particularly, data relating to imaging of tomographic images of the subject (imaging data)), and the like. Included (not shown in FIG. 2).
また、撮像装置１０は、図１に示すように、機能的には、撮像部１１（撮像手段）及び撮像データ生成部１２を有する。撮像部１１及び撮像データ生成部１２は、ＣＰＵ１０ａがＲＯＭ１０ｂ等の撮像装置１０の内蔵メモリに格納された上記コンピュータプログラムを実行し、図２に示す撮像装置１０の各構成部を動作させることによって（撮像部１１及び撮像データ生成部１２を用いて）、図７のフローチャートに示す処理を実行する。 Further, as illustrated in FIG. 1, the imaging device 10 functionally includes an imaging unit 11 (imaging unit) and an imaging data generation unit 12. The imaging unit 11 and the imaging data generation unit 12 execute the above-described computer program stored in the built-in memory of the imaging device 10 such as the ROM 10b by the CPU 10a and operate each component of the imaging device 10 illustrated in FIG. The processing shown in the flowchart of FIG. 7 is executed using the imaging unit 11 and the imaging data generation unit 12.
撮像部１１は、被検者の断層画像を撮像する機能を有する。撮像装置１０が例えばＭＲＩ装置である場合、撮像部１１は、照射装置と検出装置とを駆動させることによって磁場中に置かれた被検者に高周波（ＲＦ波）を照射し、被検者の内部の原子核に生じる核磁気共鳴現象を検出する。ここで、撮像部１１では、被検者に対して傾斜磁場を印加することによって撮像部１１による被検者の測定位置が特定され、傾斜磁場の印加により特定される測定位置と、当該測定位置での原子核の核磁気共鳴現象とを組み合わせることで、被検者を特定の方向に輪切りしたような断面の画像（断層画像）が得られる。被検者のどの方向に対して垂直な断面（例えば、被検者の体軸方向に対して垂直な断面）の断層画像を撮像するか、及び、断層画像の撮像枚数等の撮像のためのパラメータは予め撮像装置１０のオペレータ等によって予め指定されるが、被検者の大腰筋を含む腹部の撮像を目的に撮像を行う場合、一般的には、被検者の体軸方向に対して（すなわち被検者が横たわるベッドに対して）垂直な断面の断層画像（腹部横断像）や、被検者の前後方向に対して垂直な断面の断層画像（腹部前額断像）がそれぞれ複数枚撮像される。これらの複数の腹部横断像及び複数の腹部前額断像は、それぞれ被検者について互いに平行な断層画像である。撮像部１１によって撮像された断層画像は撮像データ生成部１２へ送られる。 The imaging unit 11 has a function of capturing a tomographic image of the subject. When the imaging apparatus 10 is, for example, an MRI apparatus, the imaging unit 11 irradiates a subject placed in a magnetic field by driving the irradiation apparatus and the detection apparatus, and irradiates the subject with high frequency (RF waves). Detects nuclear magnetic resonance phenomena occurring in internal nuclei. Here, in the imaging unit 11, the measurement position of the subject by the imaging unit 11 is specified by applying a gradient magnetic field to the subject, the measurement position specified by the application of the gradient magnetic field, and the measurement position By combining with the nuclear magnetic resonance phenomenon of nuclei at, a cross-sectional image (tomographic image) as if the subject was cut in a specific direction was obtained. To capture a tomographic image of a cross section perpendicular to which direction of the subject (for example, a cross section perpendicular to the body axis direction of the subject) and the number of tomographic images to be captured The parameters are specified in advance by the operator of the imaging apparatus 10 or the like. However, when imaging is performed for the purpose of imaging the abdomen including the psoas major muscle of the subject, generally the body axis direction of the subject is (I.e., to the bed on which the subject lies) a tomographic image (abdominal cross-sectional image) perpendicular to the cross-section and a tomographic image (abdominal forehead image) perpendicular to the longitudinal direction of the subject. Multiple images are taken. The plurality of abdominal cross-sectional images and the plurality of abdominal forehead images are tomographic images parallel to each other for the subject. The tomographic image captured by the imaging unit 11 is sent to the imaging data generation unit 12.
なお、本実施形態における筋断面積の算出に用いられる断層画像の枚数は断層画像の撮像に係る被検者の負担と筋断面積算出の精度とを考慮すると５〜８枚であることが好ましい。５枚未満である場合には、断層画像処理装置２０で算出される筋断面積の精度が低くなる可能性がある。また、８枚よりも多い枚数の断層画像を用いることで精度は高められる可能性があるが、撮像枚数の増加による被検者に対する負担が増加する可能性がある。また、筋断面積測定システム１における筋断面積の測定対象である大腰筋は、被検者の体軸方向とのなす角が比較的小さいことから、腹部前額断像よりも腹部横断像を後述の筋断面積の算出に用いることが好ましい。 The number of tomographic images used for calculation of the muscle cross-sectional area in the present embodiment is preferably 5 to 8 in consideration of the burden on the subject related to the tomographic image capturing and the accuracy of the muscle cross-sectional area calculation. . When the number is less than 5, the accuracy of the muscle cross-sectional area calculated by the tomographic image processing apparatus 20 may be lowered. In addition, the accuracy may be improved by using more than eight tomographic images, but the burden on the subject may increase due to an increase in the number of images to be captured. Further, the psoas major muscle, which is a measurement object of the muscle cross-sectional area in the muscle cross-sectional area measurement system 1, has a relatively small angle with the body axis direction of the subject, and therefore, the abdominal cross-sectional image rather than the abdominal frontal frame image. Is preferably used for calculation of the muscle cross-sectional area described later.
撮像データ生成部１２は、撮像部１１によって撮像された断層画像に基づいて断層画像処理装置２０へ送信するための撮像データを生成する。ここでいう撮像データとは、撮像部１１によって撮像された断層画像に対して、撮像位置を特定する情報を対応付けたものである。ここでいう撮像位置を特定する情報とは、例えば断層画像の中心位置をｘｙｚ軸による三次元座標を用いて表したもの等である。また、複数の断層画像が互いに平行でない場合には、当該断層画像の撮像方向を特定する情報等を断層画像毎に対応付けることもできる。このようにして断層画像に対してその撮像位置を特定する情報を対応付けた撮像データが撮像データ生成部１２から断層画像処理装置２０に対して送られる。 The imaging data generation unit 12 generates imaging data to be transmitted to the tomographic image processing device 20 based on the tomographic image captured by the imaging unit 11. The imaging data here is information in which information specifying the imaging position is associated with the tomographic image captured by the imaging unit 11. The information specifying the imaging position here is, for example, information representing the center position of the tomographic image using three-dimensional coordinates based on the xyz axis. Further, when a plurality of tomographic images are not parallel to each other, information for specifying the imaging direction of the tomographic image can be associated with each tomographic image. In this manner, imaging data in which information specifying the imaging position is associated with the tomographic image is sent from the imaging data generation unit 12 to the tomographic image processing apparatus 20.
なお、以下の実施形態では、被検者の左右方向をｘ軸とし、被検者の前後方向をｙ軸とし、被検者の体軸の延びる方向をｚ軸とする。そして、筋断面積測定システム１での筋断面積の算出では断層画像として上述の腹部横断像を用い、この断層画像の撮像位置及び断層画像に含まれる測定位置（断層画像における画素の位置）を示す位置情報をｘｙｚ座標を用いて表す。また、複数の断層画像として、ｚ軸方向に沿って互いに異なる複数の位置で被検者を輪切りにした８枚の断層画像Ｌ１〜Ｌ８（すなわち、ｘｙ軸方向の撮像範囲は互いに同一であり、ｚ軸方向における撮像位置のみが異なる被検者の断層画像）を用いて説明する。 In the following embodiments, the left-right direction of the subject is the x-axis, the front-rear direction of the subject is the y-axis, and the direction in which the body axis of the subject extends is the z-axis. In the calculation of the muscle cross-sectional area in the muscle cross-sectional area measurement system 1, the above-described abdominal cross-sectional image is used as the tomographic image, and the imaging position of the tomographic image and the measurement position (pixel position in the tomographic image) included in the tomographic image are determined. The position information shown is expressed using xyz coordinates. Further, as a plurality of tomographic images, eight tomographic images L1 to L8 obtained by cutting the subject at a plurality of different positions along the z-axis direction (that is, the imaging ranges in the xy-axis direction are the same as each other, A description will be given using a tomographic image of a subject that differs only in the imaging position in the z-axis direction.
次に、図１及び図２を用いて、断層画像処理装置２０の各装置の構成について説明する。断層画像処理装置２０は、図２に示すように、物理的には、ＣＰＵ２０ａ、ＲＯＭ２０ｂ、ＲＡＭ２０ｃ、通信装置２０ｄ及び記憶装置２０ｅ等を有しており、これらＣＰＵ２０ａ〜記憶装置２０ｅはバス２０ｆに接続されている。ＣＰＵ２０ａは、ＲＯＭ２０ｂ等の内蔵メモリに格納された所定のコンピュータプログラムをＲＡＭ２０ｃにロードして実行することによって、断層画像処理装置２０を統括的に制御する。通信装置２０ｄは、外部の装置（本実施形態では撮像装置１０）と通信を行うためのインターフェースを有する。記憶装置２０ｅは、書き込み／読出しが自在なメモリであって、各種コンピュータプログラム（特に、図７のフローチャートに示す処理のうち断層画像処理装置２０で行われる処理を実行するためのプログラム）を格納する。記憶装置２０ｅは、コンピュータプログラムの実行に必要な各種データ及びコンピュータプログラムの実行によって得られる各種データ（撮像データや断層画像処理装置２０における処理の結果生成されるデータ）等を格納するためのデータベースが含まれる（図２では不図示）。 Next, the configuration of each apparatus of the tomographic image processing apparatus 20 will be described with reference to FIGS. 1 and 2. As shown in FIG. 2, the tomographic image processing apparatus 20 physically includes a CPU 20a, a ROM 20b, a RAM 20c, a communication device 20d, a storage device 20e, and the like. These CPU 20a to storage device 20e are connected to a bus 20f. Has been. The CPU 20a centrally controls the tomographic image processing apparatus 20 by loading a predetermined computer program stored in a built-in memory such as the ROM 20b into the RAM 20c and executing it. The communication device 20d has an interface for communicating with an external device (the imaging device 10 in this embodiment). The storage device 20e is a readable / writable memory and stores various computer programs (particularly, a program for executing processing performed by the tomographic image processing device 20 in the processing shown in the flowchart of FIG. 7). . The storage device 20e has a database for storing various data necessary for execution of the computer program and various data obtained by execution of the computer program (imaging data and data generated as a result of processing in the tomographic image processing device 20). Included (not shown in FIG. 2).
また、断層画像処理装置２０は、図１に示すように、機能的には、撮像データ取得部２１（画像取得手段）、面積算出部２２（面積算出手段）、重心算出部２３（重心算出手段）、走行方向算出部２４（走行方向算出手段）、筋断面積算出部２５（筋断面積算出手段）及び出力部２６を有する。撮像データ取得部２１〜出力部２６は、ＣＰＵ２０ａがＲＯＭ２０ｂ等の断層画像処理装置２０の内蔵メモリに格納された上記コンピュータプログラムを実行し、図２に示す断層画像処理装置２０の各構成部を動作させることによって（撮像データ取得部２１〜出力部２６を用いて）、図７のフローチャートに示す処理を実行する。 Further, as shown in FIG. 1, the tomographic image processing apparatus 20 functionally has an imaging data acquisition unit 21 (image acquisition unit), an area calculation unit 22 (area calculation unit), and a centroid calculation unit 23 (centroid calculation unit). ), A traveling direction calculation unit 24 (traveling direction calculation unit), a muscle cross-sectional area calculation unit 25 (muscle cross-sectional area calculation unit), and an output unit 26. In the imaging data acquisition unit 21 to the output unit 26, the CPU 20a executes the computer program stored in the internal memory of the tomographic image processing apparatus 20 such as the ROM 20b, and operates each component of the tomographic image processing apparatus 20 shown in FIG. By doing so (using the imaging data acquisition unit 21 to the output unit 26), the processing shown in the flowchart of FIG. 7 is executed.
撮像データ取得部２１は、撮像装置１０から送信された被検者の撮像データを受信する。撮像データ取得部２１によって受信された情報は、断層画像処理装置２０の記憶装置２０ｅに格納されると共に面積算出部２２へ送られる。ここで、断層画像毎に対応付けられた撮像位置を特定する情報についても、断層画像と共に断層画像に対応付けて記憶装置２０ｅに格納されると共に面積算出部２２へ送られる。 The imaging data acquisition unit 21 receives the imaging data of the subject transmitted from the imaging device 10. The information received by the imaging data acquisition unit 21 is stored in the storage device 20e of the tomographic image processing device 20 and sent to the area calculation unit 22. Here, the information specifying the imaging position associated with each tomographic image is also stored in the storage device 20e in association with the tomographic image together with the tomographic image and is also sent to the area calculating unit 22.
面積算出部２２は、記憶装置２０ｅに格納された撮像データに含まれる複数の断層画像における被検者の大腰筋領域の面積を断層画像毎に算出する。面積算出部２２での被検者の大腰筋領域の面積の算出について図３及び図４を用いて具体的に説明する。一般的に筋断面積測定システム１のように被検者の断層画像を撮像しその撮像画像を利用することを目的としたシステムでは、ＤＩＣＯＭ（Digital Imaging and Communication in Medicine）規格に基づいてデータの送受信が行われる。そして、面積算出部２２で用いられる断層画像（すなわち、撮像装置１０から送信されて撮像データ取得部２１によって取得される撮像データに含まれる断層画像）は例えば図３（Ａ）に示すような画像である。面積算出部２２では、ＤＩＣＯＭ規格に基づいたデータ形式から例えばＢＭＰ（Bitmap）形式へのファイル形式の変換が行われた後、ＢＭＰ形式の断層画像に含まれる画素毎の輝度情報を用いてヒストグラムが作成される。このヒストグラムの一例を図４に示す。図４に示すヒストグラムでは、筋領域を示す低輝度側のピークＰ１と脂肪領域を示す高輝度側のピークＰ２とがある。面積算出部２２は、このヒストグラムにおける各輝度に対応する画素数を特定の範囲毎に平均化することにより得られた曲線において、ある閾値よりも高い輝度を有する画素（図４における矢印Ｆに対応する領域の輝度を有する画素）の領域を脂肪領域と判断する。閾値の設定には、例えばピークＰ１の画素数に対して５０％の画素数に対応したピークＰ１より高輝度側の輝度を閾値とする方法や、ピークＰ１が出現する輝度に３０（輝度が２５６階調の場合）を加えた輝度を閾値とする方法などがある。図３（Ｂ）は、この脂肪領域Ｆを白く塗りつぶした画像である。そして、この図３（Ｂ）に示す画像で脂肪領域Ｆとして白く塗りつぶされなかった領域のうち、大腰筋があると推定される位置に設けられた領域を大腰筋領域Ｍ１（図示左側），Ｍ２（図示右側）として塗りつぶした画像を図３（Ｃ）として示す。そして、図３（Ｃ）に示す画像で塗りつぶされた領域Ｍ１，Ｍ２の画素数をカウントすることにより断層画像における大腰筋領域の面積が算出される。面積算出部２２では、上記の処理を複数の断層画像毎に行われることで、断層画像毎の大腰筋領域の面積が算出される。ここで算出された大腰筋領域の面積は、各断層画像に対応付けられて記憶装置２０ｅに格納される。また、大腰筋領域の面積値及び断層画像に含まれる大腰筋領域Ｍ１に含まれる画素の位置情報が重心算出部２３及び筋断面積算出部２５へ送られる。 The area calculation unit 22 calculates the area of the large psoas muscle region of the subject in a plurality of tomographic images included in the imaging data stored in the storage device 20e for each tomographic image. The calculation of the area of the large psoas muscle region of the subject in the area calculation unit 22 will be specifically described with reference to FIGS. 3 and 4. In general, in a system aiming at taking a tomographic image of a subject and using the taken image like the muscle cross-sectional area measuring system 1, the data of the cross section is measured based on the DICOM (Digital Imaging and Communication in Medicine) standard. Transmission / reception is performed. The tomographic image used in the area calculation unit 22 (that is, the tomographic image included in the imaging data transmitted from the imaging device 10 and acquired by the imaging data acquisition unit 21) is an image as shown in FIG. It is. In the area calculation unit 22, after the file format is converted from the data format based on the DICOM standard to, for example, the BMP (Bitmap) format, a histogram is generated using luminance information for each pixel included in the tomographic image in the BMP format. Created. An example of this histogram is shown in FIG. In the histogram shown in FIG. 4, there is a low luminance side peak P1 indicating a muscle region and a high luminance side peak P2 indicating a fat region. In the curve obtained by averaging the number of pixels corresponding to each luminance in the histogram for each specific range, the area calculating unit 22 corresponds to a pixel having a luminance higher than a certain threshold (corresponding to an arrow F in FIG. 4). The pixel having the luminance of the region to be determined) is determined as a fat region. The threshold value can be set by, for example, a method in which the luminance on the higher luminance side than the peak P1 corresponding to 50% of the number of pixels of the peak P1 is used as a threshold, or 30 (luminance is 256) for the luminance at which the peak P1 appears. For example, there is a method of using a luminance with added gradation as a threshold value. FIG. 3B is an image in which the fat region F is painted white. Then, in the image shown in FIG. 3B, among the regions that are not painted white as the fat region F, the region provided at the position where it is estimated that the psoas muscle is present is the psoas region M1 (left side in the figure). , M2 (shown on the right side) is shown in FIG. 3C. Then, the area of the psoas muscle region in the tomographic image is calculated by counting the number of pixels in the regions M1 and M2 filled with the image shown in FIG. The area calculation unit 22 calculates the area of the psoas muscle region for each tomographic image by performing the above processing for each of the plurality of tomographic images. The area of the psoas major region calculated here is stored in the storage device 20e in association with each tomographic image. Further, the area value of the psoas muscle region and the position information of the pixels included in the psoas muscle region M1 included in the tomographic image are sent to the centroid calculation unit 23 and the muscle cross-sectional area calculation unit 25.
重心算出部２３は、断層画像に含まれている被検者の大腰筋領域の重心位置を断層画像毎に算出する。具体的には、重心算出部２３によって、面積算出部２２において大腰筋領域Ｍ１，Ｍ２として判断された領域に含まれる各画素の位置情報（ｘｙ座標）を用いて大腰筋領域Ｍ１，Ｍ２の重心位置を算出する。重心算出部２３によって画像毎に算出された大腰筋領域Ｍ１，Ｍ２の重心位置は、各断層画像に対応付けられて記憶装置２０ｅに格納されると共に走行方向算出部２４へ送られる。 The center-of-gravity calculation unit 23 calculates the position of the center of gravity of the subject's psoas muscle region included in the tomographic image for each tomographic image. Specifically, the psoas major regions M1, M2 using the position information (xy coordinates) of each pixel included in the region determined by the center of gravity calculating unit 23 as the psoas major regions M1, M2 by the area calculating unit 22. The position of the center of gravity is calculated. The barycentric positions of the psoas major regions M1 and M2 calculated for each image by the barycentric calculating unit 23 are stored in the storage device 20e in association with each tomographic image and sent to the traveling direction calculating unit 24.
走行方向算出部２４は、重心算出部２３により算出された被検者の大腰筋領域の重心位置に基づいて、大腰筋の走行方向を算出する。被検者の大腰筋の走行方向とは、図５に示すように、８枚の断層画像Ｌ１〜Ｌ８に含まれる被検者の断層画像をその撮像位置に応じてｚ軸に沿って配置し、断層画像Ｌ１〜Ｌ８のそれぞれにおける大腰筋領域Ｍ１の重心同士を結ぶことで得られる大腰筋の重心が延びる方向のことであり、走行方向算出部２４は、この大腰筋の重心が延びる方向と被検者の体軸方向とのなす角を求めることで走行方向を算出する。 The traveling direction calculation unit 24 calculates the traveling direction of the psoas muscle based on the position of the center of gravity of the psoas muscle region of the subject calculated by the gravity center calculation unit 23. As shown in FIG. 5, the traveling direction of the large psoas muscle of the subject is the arrangement of the tomographic images of the subject included in the eight tomographic images L1 to L8 along the z-axis according to the imaging position. In the tomographic images L1 to L8, the center of gravity of the psoas muscle obtained by connecting the centers of gravity of the psoas muscle region M1 extends, and the travel direction calculation unit 24 calculates the center of gravity of the psoas muscle. The traveling direction is calculated by obtaining an angle formed by the direction in which the body extends and the body axis direction of the subject.
具体的には、走行方向算出部２４は、図６に示すように、断層画像Ｌ１〜Ｌ８に含まれる大腰筋領域Ｍ１のみを取り出し、これを各大腰筋領域Ｍ１の３次元座標軸に対応する位置に配置する（図６に示すＬ１ａ〜Ｌ８ａは、断層画像Ｌ１〜Ｌ８のそれぞれにおける大腰筋領域Ｍ１を示す）。そして、走行方向算出部２４は、Ｌ１ａ〜Ｌ８ａの重心位置に近似する近似直線を最小二乗法等を用いて算出し、この近似直線をＬ１ａとＬ８ａと接する点でそれぞれ区切った線分Ｋをこの３次元座標軸に書き加える。このとき、図６における辺ａは線分Ｋをｘｙ平面に投影した辺であり、辺ｂは被検者の体軸方向に相当する辺であり、辺ｃは断層画像Ｌ１〜Ｌ８の大腰筋領域Ｌ１ａ〜Ｌ８ａの重心位置に基づく線分Ｋである。また、点Ｂのｘｙ座標（Ｘ１，Ｙ１）は、線分ＫのＬ１ａ側の端部のｘｙ座標と同じであり、点Ｃのｘｙ座標（Ｘ２，Ｙ２）は、線分ＫのＬ８ａ側の端部のｘｙ座標である。そして、これらの情報を用いて、頂点Ａでの辺ｂと辺ｃとのなす角αを求めて、これを大腰筋の走行方向とする。そして、上記の処理と同様の処理を断層画像Ｌ１〜Ｌ８の大腰筋領域Ｍ２についても行うことで、対側の大腰筋の走行方向も算出される。このように走行方向算出部２４によって算出された大腰筋の走行方向は、被検者の一連の断層画像と対応付けて記憶装置２０ｅに格納されると共に筋断面積算出部２５へ送られる。 Specifically, as shown in FIG. 6, the traveling direction calculation unit 24 extracts only the psoas muscle region M1 included in the tomographic images L1 to L8, and corresponds to the three-dimensional coordinate axis of each psoas muscle region M1. (L1a to L8a shown in FIG. 6 indicate the psoas muscle region M1 in each of the tomographic images L1 to L8). Then, the traveling direction calculation unit 24 calculates an approximate straight line that approximates the center of gravity position of L1a to L8a by using the least square method or the like, and the segment K that is divided by the points that touch the L1a and L8a, respectively. Add to the 3D coordinate axes. At this time, the side a in FIG. 6 is a side obtained by projecting the line segment K onto the xy plane, the side b is a side corresponding to the body axis direction of the subject, and the side c is the large waist of the tomographic images L1 to L8. This is a line segment K based on the barycentric positions of the muscle regions L1a to L8a. In addition, the xy coordinates (X 1 , Y 1 ) of the point B are the same as the xy coordinates of the end of the line segment K on the L1a side, and the xy coordinates (X 2 , Y 2 ) of the point C are the line segment K It is an xy coordinate of the edge part by the side of L8a. Then, using these pieces of information, an angle α formed by the side b and the side c at the vertex A is obtained, and this is set as the running direction of the psoas muscle. Then, the running direction of the contralateral psoas muscle is also calculated by performing the same process as described above for the psoas muscle region M2 of the tomographic images L1 to L8. Thus, the travel direction of the psoas major calculated by the travel direction calculation unit 24 is stored in the storage device 20e in association with a series of tomographic images of the subject and sent to the muscle cross-sectional area calculation unit 25.
筋断面積算出部２５は、走行方向算出部２４により算出された大腰筋の走行方向と、面積算出部２２により算出された断層画像毎の被検者の大腰筋領域の面積とに基づいて、走行方向に対して垂直な大腰筋の筋断面領域の面積を算出する。具体的には、撮像装置１０により撮像された被検者の断層画像Ｌ１〜Ｌ８は被検者の体軸方向（すなわち図６のｚ軸）に沿って撮像されたため、断層画像Ｌ１〜Ｌ８における大腰筋領域Ｌ１ａ〜Ｌ８ａはｘｙ平面に平行な断面であって、大腰筋の走行方向に対して垂直な断面ではない。したがって、面積算出部２２により算出された断層画像Ｌ１〜Ｌ８における大腰筋領域Ｌ１ａ〜Ｌ８ａの面積を、走行方向を示す角αを用いてそれぞれ補正する（具体的には、大腰筋領域Ｌ１ａ〜Ｌ８ａの面積に対してｃｏｓαを乗じる）ことによって、走行方向に対して垂直な断面における面積が断層画像毎に算出される。そして、筋断面積算出部２５により、断層画像Ｌ１〜Ｌ８のそれぞれにおける補正後の大腰筋領域の面積のうち最も大きいものが選ばれ、これが大腰筋の筋断面積（補正後の面積、実際には大腰筋の筋断面に対応する部分の面積である。）とされる。このようにして筋断面積算出部２５により算出された大腰筋の筋断面に対応する部分の面積は、被検者についての一連の断層画像に対応付けて記憶装置２０ｅに格納されると共に出力部２６へ送られる。 The muscle cross-sectional area calculator 25 is based on the travel direction of the psoas major calculated by the travel direction calculator 24 and the area of the psoas muscle region of the subject for each tomographic image calculated by the area calculator 22. Thus, the area of the muscle cross-sectional area of the psoas major perpendicular to the running direction is calculated. Specifically, since the tomographic images L1 to L8 of the subject imaged by the imaging device 10 were imaged along the body axis direction of the subject (that is, the z axis in FIG. 6), the tomographic images L1 to L8 The psoas muscle regions L1a to L8a are cross sections parallel to the xy plane and are not cross sections perpendicular to the running direction of the psoas muscle. Therefore, the areas of the psoas muscle regions L1a to L8a in the tomographic images L1 to L8 calculated by the area calculation unit 22 are respectively corrected using the angle α indicating the running direction (specifically, the psoas muscle region L1a). By multiplying the area of .about.L8a by cos.alpha.), The area in the cross section perpendicular to the traveling direction is calculated for each tomographic image. The muscle cross-sectional area calculation unit 25 selects the largest area of the corrected psoas muscle region in each of the tomographic images L1 to L8, and this is the muscle cross-sectional area (corrected area, Actually, it is the area of the part corresponding to the muscle cross section of the psoas major). Thus, the area of the portion corresponding to the muscle cross section of the psoas major calculated by the muscle cross-sectional area calculation unit 25 is stored in the storage device 20e in association with a series of tomographic images of the subject and output. Sent to the unit 26.
出力部２６は、筋断面積算出部２５により算出された大腰筋の筋断面積（補正後の面積）を被検者や断層画像処理装置２０のオペレータ等へ通知する機能を有する。具体的には、例えば、出力部２６は、断層画像処理装置２０のモニタ等に表示する方法や、プリンタ等に出力する方法等を用いて、筋断面積算出部２５により算出された大腰筋の筋断面積（補正後の面積）を出力する。このとき、出力部２６は、記憶装置２０ｅに格納された他の情報（断層画像、大腰筋の走行方向等）についても筋断面に対応する部分の面積の算出結果と対応させて出力することもできる。出力部２６によってどの情報をどのような形式で出力するかは、オペレータによって予め設定される態様としてもよいし、測定毎に変更する態様としてもよい。 The output unit 26 has a function of notifying the examinee, the operator of the tomographic image processing apparatus 20 and the like of the muscle cross-sectional area (corrected area) of the psoas major calculated by the muscle cross-sectional area calculation unit 25. Specifically, for example, the output unit 26 uses the method of displaying on the monitor or the like of the tomographic image processing apparatus 20 or the method of outputting to a printer or the like, and the psoas muscle calculated by the muscle cross-sectional area calculation unit 25. The cross-sectional area (area after correction) is output. At this time, the output unit 26 outputs other information (tomographic image, travel direction of the psoas muscle, etc.) stored in the storage device 20e in association with the calculation result of the area corresponding to the muscle cross section. You can also. Which information is output in what format by the output unit 26 may be set in advance by the operator, or may be changed every measurement.
次に、図７を参照して、筋断面積測定システム１の動作を説明する。まず、筋断面積測定システム１の撮像装置１０の撮像部１１によって被検者の断層画像が撮影される（Ｓ０１）。撮像部１１によって撮像された被検者の断層画像は、撮像データ生成部１２により撮像データとして断層画像処理装置２０に送信され、断層画像処理装置２０の撮像データ取得部２１によって受信される（画像取得ステップ）。 Next, the operation of the muscle cross-sectional area measurement system 1 will be described with reference to FIG. First, a tomographic image of a subject is taken by the imaging unit 11 of the imaging device 10 of the muscle cross-sectional area measurement system 1 (S01). The tomographic image of the subject imaged by the imaging unit 11 is transmitted as imaging data to the tomographic image processing device 20 by the imaging data generation unit 12 and received by the imaging data acquisition unit 21 of the tomographic image processing device 20 (image). Acquisition step).
次に、面積算出部２２によって、撮像データ取得部２１で取得された撮像データのうち断層画像を用いて、被検者の大腰筋領域の面積が画像毎に算出される（Ｓ０２：面積算出ステップ）。次に、重心算出部２３によって、断層画像毎に大腰筋領域の重心位置が算出される（Ｓ０３：重心算出ステップ）。次に、走行方向算出部２４によって、大腰筋領域の重心位置の情報を用いて、大腰筋の走行方向が算出される（Ｓ０４：走行方向算出ステップ）。そして、筋断面積算出部２５によって、面積算出部２２で算出された断層画像毎の大腰筋領域の面積と、走行方向算出部２４で算出された大腰筋の走行方向と、に基づいて、大腰筋の筋断面に対応する部分の面積が算出される（Ｓ０５：筋断面積算出ステップ）。そして、筋断面積算出部２５で算出された筋断面に対応する部分の面積が出力部２６によって出力されることによって、被検者や断層画像処理装置２０のオペレータがその結果（大腰筋の筋断面積の算出結果）を知ることができる（Ｓ０６：出力ステップ）。以上の処理が筋断面積測定システム１によって行われることで、被検者の大腰筋の筋断面積が算出される。 Next, the area calculation unit 22 calculates the area of the psoas major region of the subject for each image using the tomographic image of the imaging data acquired by the imaging data acquisition unit 21 (S02: area calculation). Step). Next, the barycentric position of the psoas major region is calculated for each tomographic image by the barycentric calculating unit 23 (S03: barycentric calculating step). Next, the travel direction calculation unit 24 calculates the travel direction of the psoas muscle using the information on the center of gravity of the psoas muscle region (S04: travel direction calculation step). Then, based on the area of the psoas muscle region for each tomographic image calculated by the area calculator 22 by the muscle cross-sectional area calculator 25 and the running direction of the psoas muscle calculated by the running direction calculator 24. Then, the area of the part corresponding to the muscle cross section of the psoas major is calculated (S05: muscle cross-sectional area calculating step). Then, the area of the portion corresponding to the muscle cross section calculated by the muscle cross-sectional area calculating unit 25 is output by the output unit 26, so that the examinee or the operator of the tomographic image processing apparatus 20 can obtain the result (the psoas major muscle). It is possible to know the calculation result of the muscle cross-sectional area (S06: output step). By performing the above processing by the muscle cross-sectional area measurement system 1, the muscle cross-sectional area of the large psoas muscle of the subject is calculated.
ここで、上記の筋断面積測定システム１によって算出された大腰筋の断面積の精度について、実施例を参照しながら説明する。表１は、被検者の体軸方向に対して垂直な断面で被検者の大腰筋が撮像された８枚の腹部横断像（断層画像）における大腰筋領域の面積と、この断層画像を用いて本実施形態で説明した方法により大腰筋の走行方向を求め、この走行方向を用いて各断層画像の大腰筋領域の面積を補正した面積（筋断面に対応する部分の面積：補正後の面積）と、被検者の大腰筋の断面積の実測値と、を比較したものである。なお、被検者の大腰筋の断面積の実測値とは、被検者の大腰筋近傍における複数の腹部前額断像と複数の腹部横断像とを撮像しこれらの断層画像から被検者の筋の走行方向を正確に求めた上で、横断像１〜８が撮像された位置毎に大腰筋の走行方向に対して垂直な面の断層画像を撮像し、その横断像から大腰筋に対応する部分の面積を算出したものである。そして、断層画像における大腰筋領域の面積及び補正後の面積と、実測値とを比較した結果を面積及び百分率で示す。なお、表１では左右の大腰筋のうちの一方のみについて示す。この８枚の横断像のうち、大腰筋の両端間のｘｙ平面における距離（図６に示す辺ａの長さ）は３．９９６ｃｍであり、８枚の断層画像から求められた大腰筋の走行方向（図６に示す角αの大きさ）は１７．８４度であった。 Here, the accuracy of the cross-sectional area of the psoas major calculated by the muscle cross-sectional area measurement system 1 will be described with reference to the examples. Table 1 shows the area of the psoas muscle region in the eight abdominal cross-sectional images (tomographic images) obtained by imaging the subject's psoas muscle in a cross section perpendicular to the body axis direction of the subject, Using the image, the travel direction of the psoas muscle is obtained by the method described in this embodiment, and the area of the psoas muscle region of each tomographic image is corrected using this travel direction (the area of the portion corresponding to the muscle cross section). : Area after correction) and the measured value of the cross-sectional area of the psoas major of the subject. Note that the measured value of the cross-sectional area of the patient's psoas muscle is the number of abdominal forehead images and a plurality of abdominal cross-sectional images in the vicinity of the subject's psoas muscle that are taken from these tomographic images. After accurately determining the traveling direction of the examiner's muscle, a tomographic image of a plane perpendicular to the traveling direction of the psoas major is taken at each position where the transverse images 1 to 8 are taken, and from the transverse image The area of the part corresponding to the psoas major is calculated. And the result of having compared the area of the psoas muscle area in a tomographic image, the area after correction | amendment, and an actual measurement value is shown by an area and a percentage. In Table 1, only one of the left and right psoas muscles is shown. Among these eight cross-sectional images, the distance in the xy plane between both ends of the psoas muscle (the length of the side a shown in FIG. 6) is 3.996 cm, and the psoas muscle determined from the eight tomographic images. The traveling direction (the angle α shown in FIG. 6) was 17.84 degrees.
表１に示すように、体軸方向に垂直な断面の断層画像を用いて本実施形態で説明した方法により大腰筋の走行方向を求め、この走行方向を用いて各断層画像の大腰筋領域の面積を補正した面積（補正後の面積）と、大腰筋の断面積の実測値との差は±１％程度であり、補正前（断層画像における大腰筋領域の面積）と実測値との差が４〜７％あったのと比較すると、誤差が大幅に低減され、高い精度で大腰筋の面積を補正することができ、高精度で筋断面積を算出できることが確認された。 As shown in Table 1, the running direction of the psoas muscle is obtained by the method described in the present embodiment using the tomographic image of the cross section perpendicular to the body axis direction, and the psoas muscle of each tomographic image is obtained using this running direction. The difference between the corrected area of the area (the area after correction) and the measured value of the cross-sectional area of the psoas major is about ± 1%. It is measured before and after correction (the area of the psoas major in the tomographic image). Compared with the difference of 4-7%, the error is greatly reduced, the area of the psoas muscle can be corrected with high accuracy, and the muscle cross-sectional area can be calculated with high accuracy. It was.
そして、本実施形態で説明したように、複数毎の断層画像から大腰筋の走行方向を求め、この走行方向を用いて各断層画像の大腰筋領域の面積を補正することで筋断面積を算出する方法では、筋断面積の実測値を求める場合のように複数の方向から被検者の断層画像を撮像する必要がなく、より短時間且つ低コストで断層画像の撮像から大腰筋の筋断面積の算出までの一連の処理を行うことが確認された。 Then, as described in the present embodiment, the cross-sectional area of the muscle is obtained by determining the running direction of the psoas muscle from the plurality of tomographic images and correcting the area of the psoas muscle region of each tomographic image using this running direction. In the method for calculating the tomographic image, it is not necessary to take a tomographic image of the subject from a plurality of directions as in the case of obtaining an actual measurement value of the muscle cross-sectional area. It was confirmed that a series of processes up to the calculation of the muscle cross-sectional area was performed.
このように、本実施形態に示す筋断面積測定システム１によれば、断層画像処理装置２０によって、被検者の大腰筋領域の面積及び重心位置が撮像装置１０によって撮像された被検者の断層画像毎に算出され、この大腰筋領域の重心位置に基づいて被検者の大腰筋の走行方向が算出される。そして、断層画像処理装置２０の筋断面積算出部２５では、大腰筋の走行方向と複数の断層画像のそれぞれにおける被検者の大腰筋領域の面積とに基づいて、被検者の大腰筋の筋断面に対応する部分の面積が算出される。このように、筋断面積測定システム１では、被検者の複数の断層画像に基づいてその被検者の大腰筋の走行方向が求められ、得られた走行方向に基づいて被検者の大腰筋の筋断面積が算出されるため、被検者の体格や測定時の体勢に応じた補正を行うことができる。したがって、被検者の体格や測定時の体勢等に関係なく高い精度で筋断面積の算出を行うことができる。そして、上記の筋断面積測定システム１による断層画像処理により筋断面積を算出する方法をトレーニングやリハビリによる効果の評価対象として好適に用いられる大腰筋に対して適用することで、被検者の大腰筋の微小な変化を適切に計測することができ、被検者の筋肉の状態をより適切に把握することができる。 As described above, according to the muscle cross-sectional area measurement system 1 shown in the present embodiment, the tomographic image processing device 20 captures the area and the center of gravity position of the subject's psoas muscle region by the imaging device 10. The tomographic image is calculated for each of the tomographic images, and the running direction of the subject's psoas muscle is calculated based on the position of the center of gravity of the psoas muscle region. Then, the muscle cross-sectional area calculation unit 25 of the tomographic image processing device 20 determines the large size of the subject based on the traveling direction of the psoas major and the area of the psoas major region of the subject in each of the plurality of tomographic images. The area of the portion corresponding to the muscle cross section of the lower back is calculated. As described above, in the muscle cross-sectional area measurement system 1, the traveling direction of the psoas major of the subject is obtained based on the plurality of tomographic images of the subject, and the subject's traveling direction is determined based on the obtained traveling direction. Since the muscle cross-sectional area of the psoas major is calculated, it is possible to perform correction according to the physique of the subject and the posture at the time of measurement. Therefore, the muscle cross-sectional area can be calculated with high accuracy regardless of the physique of the subject and the posture at the time of measurement. Then, by applying the method for calculating the muscle cross-sectional area by the tomographic image processing by the muscle cross-sectional area measuring system 1 described above to the psoas muscle that is preferably used as an evaluation target of the effect by training and rehabilitation, the subject The minute change of the psoas major muscle can be measured appropriately, and the state of the subject's muscle can be grasped more appropriately.
また、大腰筋の断層画像を種々の方向から撮像することで筋の走行方向を３次元的に算出することで筋断面積の実測値を得る方法と比較して、被検者の負担を減らすことができる。具体的には、長時間の撮像による被検者の疲労や、磁場（ＭＲＩの場合）又は放射線の曝露（Ｘ線ＣＴの場合）等による被検者への影響を減らすことができる。また、本実施形態の筋断面積測定システム１では、大腰筋の筋断面積の算出に必要な断層画像を短時間で得ることができるため、より多くの被検者に対して測定を行うことができる。 Compared with the method of obtaining the measured value of the muscle cross-sectional area by three-dimensionally calculating the running direction of the muscle by capturing tomographic images of the psoas muscle from various directions, the burden on the subject is reduced. Can be reduced. Specifically, it is possible to reduce the influence on the subject due to fatigue of the subject due to long-time imaging, magnetic field (in the case of MRI), or exposure to radiation (in the case of X-ray CT). Further, in the muscle cross-sectional area measurement system 1 of the present embodiment, a tomographic image necessary for calculating the muscle cross-sectional area of the psoas major muscle can be obtained in a short time, so that more subjects are measured. be able to.
以上、本発明の実施形態について説明したが、本実施形態に係る筋断面積測定システム１は種々の変更を行うことができる。例えば、上記実施形態では、撮像装置１０と断層画像処理装置２０とが別体である場合について説明したが、これらが一体化された装置であってもよい。また、上記実施形態では、撮像装置１０で撮像された断層画像を含む撮像データが撮像装置１０から送信され、断層画像処理装置２０によって受信されることで断層画像処理装置２０は撮像データを取得する構成について説明したが、例えば過去に撮像された被検者の断層画像に係る情報が書き込まれた外部ファイルをオペレータが断層画像処理装置２０に読み込ませることで、撮像データ取得部２１が撮像データを取得する態様とすることもできる。すなわち、本発明では撮像装置１０は必須の構成ではない。 As mentioned above, although embodiment of this invention was described, the muscle cross-sectional area measurement system 1 which concerns on this embodiment can perform a various change. For example, in the above-described embodiment, the case where the imaging device 10 and the tomographic image processing device 20 are separate has been described. However, an apparatus in which these are integrated may be used. In the above-described embodiment, the tomographic image processing apparatus 20 acquires the imaging data when the imaging data including the tomographic image captured by the imaging apparatus 10 is transmitted from the imaging apparatus 10 and received by the tomographic image processing apparatus 20. The configuration has been described. For example, when the operator reads an external file in which information related to a tomographic image of a subject imaged in the past is written into the tomographic image processing apparatus 20, the imaging data acquisition unit 21 captures the imaging data. It can also be set as the aspect acquired. That is, in the present invention, the imaging device 10 is not an essential configuration.
また、上記実施形態では、被検者の大腰筋が撮像された複数の断層画像として、被検者の体軸方向に対して垂直な断面の断層画像を用いて説明したが、他の方向の断層画像を用いても大腰筋の筋断面積の撮像を行うことができる。その場合であっても、断層画像の大腰筋領域の重心位置から算出される筋の走行方向と断層画像の撮像方向とのなす角を算出し、この結果に基づいて断層画像における大腰筋領域の面積を補正することで、大腰筋の筋断面積を算出することができる。また、被検者の断面の撮像方向が互いに異なる複数の断層画像を用いて筋断面積の撮像を行うことも可能であるが、この場合にはそれぞれの断層画像がどの方向から被検者の断面を撮像したものであるかを特定する情報を用いて、各断層画像の補正を行う必要があるため、互いに平行な断層画像を用いて筋断面積の算出を行う場合によりも処理量が増大する可能性がある。 Further, in the above-described embodiment, the cross-sectional tomographic image perpendicular to the body axis direction of the subject has been described as a plurality of tomographic images obtained by imaging the subject's psoas muscle. The tomographic image of the psoas major can also be imaged using the tomographic image of the above. Even in that case, the angle between the muscle running direction calculated from the barycentric position of the psoas muscle region of the tomographic image and the imaging direction of the tomographic image is calculated, and based on this result, the psoas muscle in the tomographic image is calculated. By correcting the area of the region, the muscle cross-sectional area of the psoas major can be calculated. It is also possible to capture a muscle cross-sectional area using a plurality of tomographic images having different imaging directions of the cross section of the subject. Since it is necessary to correct each tomographic image using information that identifies whether a cross-sectional image is taken, the amount of processing increases even when the muscle cross-sectional area is calculated using mutually parallel tomographic images there's a possibility that.
１…筋断面積測定システム、１０…撮像装置、１１…撮像部、１２…撮像データ生成部、２０…断層画像処理装置、２１…撮像データ取得部、２２…面積算出部、２３…重心算出部、２４…走行方向算出部、２５…筋断面積算出部、２６…出力部、Ｌ１〜Ｌ８…断層画像、Ｌ１ａ〜Ｌ８ａ，Ｍ１，Ｍ２…大腰筋領域
DESCRIPTION OF SYMBOLS 1 ... Muscle cross-sectional area measurement system, 10 ... Imaging device, 11 ... Imaging part, 12 ... Imaging data generation part, 20 ... Tomographic image processing apparatus, 21 ... Imaging data acquisition part, 22 ... Area calculation part, 23 ... Gravity center calculation part , 24 ... Traveling direction calculation unit, 25 ... Muscle cross-sectional area calculation unit, 26 ... Output unit, L1-L8 ... Tomographic image, L1a-L8a, M1, M2 ... Large psoas muscle region
Claims (4)
前記断層画像に含まれている前記被検者の大腰筋領域の面積を前記断層画像毎に算出する面積算出手段と、
前記断層画像に含まれている前記被検者の大腰筋領域の重心位置を前記断層画像毎に算出する重心算出手段と、
前記重心算出手段により算出された前記重心位置に基づいて、前記大腰筋の走行方向を算出する走行方向算出手段と、
前記走行方向算出手段により算出された前記走行方向と前記面積算出手段により算出された前記断層画像毎の前記大腰筋領域の面積とに基づいて、前記走行方向に対して垂直な前記大腰筋の筋断面に対応する部分の面積を算出する筋断面積算出手段と、
を備えることを特徴とする断層画像処理装置。 Image acquisition means for acquiring a plurality of tomographic images obtained by imaging the psoas major muscle of the subject;
An area calculating means for calculating the area of the psoas major region of the subject included in the tomographic image for each tomographic image;
Centroid calculating means for calculating, for each tomographic image, the centroid position of the psoas major region of the subject contained in the tomographic image;
A traveling direction calculating means for calculating a traveling direction of the psoas major, based on the position of the center of gravity calculated by the gravity center calculating means;
Based on the travel direction calculated by the travel direction calculation means and the area of the psoas muscle region for each tomographic image calculated by the area calculation means, the psoas major perpendicular to the travel direction Muscle cross-sectional area calculating means for calculating the area of the portion corresponding to the muscle cross-section,
A tomographic image processing apparatus comprising:
前記走行方向算出手段は、前記複数の断層画像における前記被検者の大腰筋領域の重心位置と当該断層画像の撮像位置とを用いて前記走行方向を算出する
ことを特徴とする請求項１記載の断層画像処理装置。 The plurality of tomographic images are obtained by imaging cross sections of the psoas major parallel to each other,
The travel direction calculation means calculates the travel direction using a center of gravity position of the psoas major region of the subject in the plurality of tomographic images and an imaging position of the tomographic image. The tomographic image processing apparatus described.
前記撮像装置は、
被検者の大腰筋の複数の断層画像を撮像する撮像手段を備え、
前記断層画像処理装置は、
前記撮像装置で撮像された前記複数の断層画像を取得する画像取得手段と、
前記断層画像に含まれている前記被検者の大腰筋領域の面積を前記断層画像毎に算出する面積算出手段と、
前記断層画像に含まれている前記被検者の大腰筋領域の重心位置を前記断層画像毎に算出する重心算出手段と、
前記重心算出手段により算出された前記重心位置に基づいて、前記大腰筋の走行方向を算出する走行方向算出手段と、
前記走行方向算出手段により算出された前記走行方向と前記面積算出手段により算出された前記断層画像毎の前記大腰筋領域の面積とに基づいて、前記走行方向に対して垂直な前記大腰筋の筋断面に対応する部分の面積を算出する筋断面積算出手段と、
を備えることを特徴とする筋断面積測定システム。 A muscle cross-sectional area measuring system including an imaging device and a tomographic image processing device,
The imaging device
Comprising imaging means for imaging a plurality of tomographic images of the psoas major of the subject;
The tomographic image processing apparatus comprises:
Image acquisition means for acquiring the plurality of tomographic images captured by the imaging device;
An area calculating means for calculating the area of the psoas major region of the subject included in the tomographic image for each tomographic image;
Centroid calculating means for calculating, for each tomographic image, the centroid position of the psoas major region of the subject contained in the tomographic image;
A traveling direction calculating means for calculating a traveling direction of the psoas major, based on the position of the center of gravity calculated by the gravity center calculating means;
Based on the travel direction calculated by the travel direction calculation means and the area of the psoas muscle region for each tomographic image calculated by the area calculation means, the psoas major perpendicular to the travel direction Muscle cross-sectional area calculating means for calculating the area of the portion corresponding to the muscle cross-section,
A muscle cross-sectional area measurement system comprising:
前記断層画像処理装置が、被検者の大腰筋が撮像された複数の断層画像を取得する画像取得ステップと、
前記断層画像処理装置が、前記断層画像に含まれている前記被検者の大腰筋領域の面積を前記断層画像毎に算出する面積算出ステップと、
前記断層画像処理装置が、前記断層画像に含まれている前記被検者の大腰筋領域の重心位置を前記断層画像毎に算出する重心算出ステップと、
前記断層画像処理装置が、前記重心算出ステップにおいて算出された前記重心位置に基づいて、前記大腰筋の走行方向を算出する走行方向算出ステップと、
前記断層画像処理装置が、前記走行方向算出ステップにおいて算出された前記走行方向と前記面積算出ステップにおいて算出された前記断層画像毎の前記大腰筋領域の面積とに基づいて、前記走行方向に対して垂直な前記大腰筋の筋断面に対応する部分の面積を算出する筋断面積算出ステップと、
を備えることを特徴とする断層画像処理方法。 A tomographic image processing method using a tomographic image processing apparatus,
The tomographic image processing apparatus obtains a plurality of tomographic images in which the subject's psoas major is imaged; and
An area calculating step in which the tomographic image processing apparatus calculates an area of the psoas major region of the subject included in the tomographic image for each tomographic image;
A center-of-gravity calculation step in which the tomographic image processing apparatus calculates a center-of-gravity position of the psoas major region of the subject included in the tomographic image for each tomographic image;
The tomographic image processing apparatus calculates a traveling direction of the psoas muscle based on the center of gravity position calculated in the center of gravity calculating step;
The tomographic image processing apparatus is configured for the traveling direction based on the traveling direction calculated in the traveling direction calculating step and the area of the psoas muscle region for each tomographic image calculated in the area calculating step. A muscle cross-sectional area calculating step for calculating an area of a portion corresponding to a muscle cross-section of the psoas major perpendicular to
A tomographic image processing method comprising:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009213565A JP2011062253A (en) | 2009-09-15 | 2009-09-15 | Tomographic image processor, muscle cross-sectional area measurement system, and tomographic image processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009213565A JP2011062253A (en) | 2009-09-15 | 2009-09-15 | Tomographic image processor, muscle cross-sectional area measurement system, and tomographic image processing method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2011062253A true JP2011062253A (en) | 2011-03-31 |
Family
ID=43949164
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009213565A Pending JP2011062253A (en) | 2009-09-15 | 2009-09-15 | Tomographic image processor, muscle cross-sectional area measurement system, and tomographic image processing method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2011062253A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012157406A1 (en) * | 2011-05-18 | 2012-11-22 | 株式会社 日立メディコ | Image analysis device, program, and image-capturing device |
| WO2017179256A1 (en) * | 2016-04-13 | 2017-10-19 | 日本メジフィジックス株式会社 | Automatic removal of physiological accumulations from nuclear medicine image, and automatic segmentation of ct image |
| CN108882907A (en) * | 2016-04-13 | 2018-11-23 | 日本医事物理股份有限公司 | Automatic estimating method, device and the computer program in the bone region in CT image |
-
2009
- 2009-09-15 JP JP2009213565A patent/JP2011062253A/en active Pending
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012157406A1 (en) * | 2011-05-18 | 2012-11-22 | 株式会社 日立メディコ | Image analysis device, program, and image-capturing device |
| WO2017179256A1 (en) * | 2016-04-13 | 2017-10-19 | 日本メジフィジックス株式会社 | Automatic removal of physiological accumulations from nuclear medicine image, and automatic segmentation of ct image |
| CN108882907A (en) * | 2016-04-13 | 2018-11-23 | 日本医事物理股份有限公司 | Automatic estimating method, device and the computer program in the bone region in CT image |
| JPWO2017179255A1 (en) * | 2016-04-13 | 2019-02-21 | 日本メジフィジックス株式会社 | Method, apparatus, and computer program for automatically estimating bone region in CT image |
| EP3443906A4 (en) * | 2016-04-13 | 2019-11-20 | Nihon Medi-Physics Co., Ltd. | Method, device and computer program for automatic estimation of bone region in ct image |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Livanelioglu et al. | The validity and reliability of “Spinal Mouse” assessment of spinal curvatures in the frontal plane in pediatric adolescent idiopathic thoraco-lumbar curves | |
| CN104582578B (en) | Quantitative calibration system is integrated in computer tomography scanner | |
| Fiorentino et al. | Accuracy of functional and predictive methods to calculate the hip joint center in young non-pathologic asymptomatic adults with dual fluoroscopy as a reference standard | |
| Munn et al. | Changes in face topography from supine-to-upright position—and soft tissue correction values for craniofacial identification | |
| Kim et al. | Body fat assessment method using CT images with separation mask algorithm | |
| Li et al. | Could clinical ultrasound improve the fitting of spinal orthosis for the patients with AIS? | |
| WO2009107500A1 (en) | Medical image processing device, medical image processing method, and program | |
| Girdler et al. | Emerging techniques in diagnostic imaging for idiopathic scoliosis in children and adolescents: a review of the literature | |
| Kibsgård et al. | Radiosteriometric analysis of movement in the sacroiliac joint during a single-leg stance in patients with long-lasting pelvic girdle pain | |
| Ilharreborde et al. | Spinal penetration index assessment in adolescent idiopathic scoliosis using EOS low-dose biplanar stereoradiography | |
| Muyor et al. | Test-retest reliability and validity of a motion capture (MOCAP) system for measuring thoracic and lumbar spinal curvatures and sacral inclination in the sagittal plane | |
| Jiang et al. | A fast 3-D ultrasound projection imaging method for scoliosis assessment | |
| Zhao et al. | Three-dimensional analysis of cervical spine segmental motion in rotation | |
| Hu et al. | Value and clinical application of orthopedic metal artifact reduction algorithm in CT scans after orthopedic metal implantation | |
| Xu et al. | Back shape measurement and three-dimensional reconstruction of spinal shape using one kinect sensor | |
| Bauer et al. | MRI-derived body segment parameters of children differ from age-based estimates derived using photogrammetry | |
| Courvoisier et al. | Evaluation of a three-dimensional reconstruction method of the rib cage of mild scoliotic patients | |
| Kinel et al. | Normative 3D opto-electronic stereo-photogrammetric sagittal alignment parameters in a young healthy adult population | |
| Kandasamy et al. | Posture and back shape measurement tools: A narrative | |
| Passmore et al. | Measuring femoral torsion in vivo using freehand 3-D ultrasound imaging | |
| Liu et al. | Development of a 3-dimensional back contour imaging system for monitoring scoliosis progression in children | |
| JP2011062253A (en) | Tomographic image processor, muscle cross-sectional area measurement system, and tomographic image processing method | |
| Tay et al. | A method for in-vivo kinematic analysis of the forearm | |
| KR101508178B1 (en) | Scoliosis analysis system and method the same | |
| Humbert et al. | 3D Evaluation of the acetabular coverage assessed by biplanar X-rays or single anteroposterior X-ray compared with CT-scan |