JP6448142B2 - 放射線照射による皮膚変化予測装置と検証装置 - Google Patents

放射線照射による皮膚変化予測装置と検証装置 Download PDF

Info

Publication number
JP6448142B2
JP6448142B2 JP2016510041A JP2016510041A JP6448142B2 JP 6448142 B2 JP6448142 B2 JP 6448142B2 JP 2016510041 A JP2016510041 A JP 2016510041A JP 2016510041 A JP2016510041 A JP 2016510041A JP 6448142 B2 JP6448142 B2 JP 6448142B2
Authority
JP
Japan
Prior art keywords
skin
change
image
unit
skin image
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.)
Expired - Fee Related
Application number
JP2016510041A
Other languages
English (en)
Other versions
JPWO2015146164A1 (ja
Inventor
礼明 松原
礼明 松原
成弘 松藤
成弘 松藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NATIONAL INSTITUTES FOR QUANTUM AND RADIOLOGICALSCIENCE AND TECHNOLOGY
Original Assignee
NATIONAL INSTITUTES FOR QUANTUM AND RADIOLOGICALSCIENCE AND TECHNOLOGY
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 INSTITUTES FOR QUANTUM AND RADIOLOGICALSCIENCE AND TECHNOLOGY filed Critical NATIONAL INSTITUTES FOR QUANTUM AND RADIOLOGICALSCIENCE AND TECHNOLOGY
Publication of JPWO2015146164A1 publication Critical patent/JPWO2015146164A1/ja
Application granted granted Critical
Publication of JP6448142B2 publication Critical patent/JP6448142B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/1032—Determining colour for diagnostic purposes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1079—Measuring physical dimensions, e.g. size of the entire body or parts thereof using optical or photographic means
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/44—Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
    • A61B5/441—Skin evaluation, e.g. for skin disorder diagnosis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/44—Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
    • A61B5/441—Skin evaluation, e.g. for skin disorder diagnosis
    • A61B5/445—Evaluating skin irritation or skin trauma, e.g. rash, eczema, wound, bed sore
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4848—Monitoring or testing the effects of treatment, e.g. of medication
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103—Treatment planning systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/90—Determination of colour characteristics
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10024—Color image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30088—Skin; Dermal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Dermatology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Psychiatry (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Urology & Nephrology (AREA)
  • Quality & Reliability (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Radiation-Therapy Devices (AREA)

Description

この発明は、例えば、放射線照射による皮膚変化を事前に予測するような放射線照射による皮膚変化予測装置と皮膚変化予測プログラム、および、予測結果と実際の皮膚変化を検証するような検証装置と検証プログラムに関する。
人や動物等の生体の外から放射線を照射する放射線治療において、治療の対象となる腫瘍は、多くの場合に体内に存在する。このため、放射線治療は、放射線の通り道となる皮膚への被曝という副作用を常に伴う。
放射線治療後の皮膚は、放射線の照射領域において、急性の皮膚反応が生じることによって赤く炎症している。従って、放射線治療によって腫瘍に投与する放射線量は、この皮膚反応によって制限される。また、この皮膚反応は、治療後の患者の生活の質（クオリティオブライフ：ＱＯＬ）とも密接に関係している。このため、放射線治療によって生じる皮膚反応を事前に予測することができれば、治療で照射する放射線量の決定を容易にでき、また、治療後の患者の生活の質を向上することができる。
ここで、特許文献を調査したところ、放射線治療用放射線モニタシステムが提案されていた（特許文献１参照）。このシステムは、放射線治療を受けている患者に浴びせられる放射線の吸収線量を推定することを目的としている。
しかし、このシステムは実際に照射された放射線の吸収線量を測定しようとするものであり、放射線照射による皮膚反応を事前に予測するようなものではなかった。
特表２００５−５１２０２８号公報
この発明は、上述の問題に鑑みて、放射線照射による皮膚反応を事前に精度よく予測できる放射線照射による皮膚変化予測装置、検証装置、放射線照射による皮膚変化予測プログラム、および検証プログラムを提供することを目的とする。
この発明は、照射予定の放射線の放射線情報の入力を受け付ける放射線情報受付部と、生体の皮膚を撮影した皮膚画像を取得する皮膚画像取得部と、前記放射線情報により定まる放射線を照射したことによる前記皮膚の変化を演算して前記皮膚画像から変化後の予測皮膚画像を得る変化演算部と、前記予測皮膚画像を出力する出力部とを備え、前記変化演算部は、前記皮膚画像を、元の色成分による表現形式である皮膚画像色成分データから生体要素成分による表現形式である皮膚画像生体要素成分データに変換する表現形式変換部と、一部の前記皮膚画像生体要素成分データが前記放射線情報で定められる放射線照射によって変化する変化量を決定する変化量決定部と、前記変化量に応じて前記一部の皮膚画像生体要素成分データを変化させる変化部と、変化後の画像を元のカラー画像の表現形式に復元する表現形式復元部とを有する放射線照射による皮膚変化予測装置、これを利用した検証装置、これらのプログラム、および皮膚変化予測方法であることを特徴とする。
この発明により、放射線照射による皮膚反応を事前に精度よく予測できる放射線照射による皮膚変化予測装置、検証装置、放射線照射による皮膚変化予測プログラム、および検証プログラムを提供できる。
皮膚変化予測システムの構成を示すブロック図。 放射線照射によるグレード決定の説明図。 皮膚変化の説明図。 グレード判定画面の画面構成図。 皮膚変化検証システムの構成を示すブロック図。 照射領域の検証用画像を示す説明図。 皮膚領域抽出効果を説明する説明図。
光子線治療および陽子線・重粒子線治療等の放射線治療において、放射線照射による皮膚反応は、一般に次の４段階のグレードで表される。
グレード１：軽い紅斑
グレード２：強い紅斑
グレード３：水泡、糜爛（ビラン）
グレード４：潰瘍
このグレード表現は、非常に荒いものであり、医師が皮膚の状態を目視確認してグレードを判別しているのが現状である。このため、このグレードを定量的に精度よく把握できるようにし、医師間でも共通尺度で意思疎通を図れるようにする必要があった。
発明者らは、放射線照射による皮膚反応を放射線治療前に精度よく予測できる方法を研究した。当初は、皮膚反応によって肌が赤くなることに着目し、カメラで撮影したＲＧＢ表現形式の画像からＲ成分（赤色成分）を取り出すことで、皮膚が被曝する放射線量（以下、皮膚線量という。）と皮膚反応の関係を検証して放射線照射後の皮膚画像を作成することを検討した。
しかし、実際に試験を行い、照射による皮膚線量と、照射後の皮膚の撮影画像から得られる皮膚反応の程度との相関性を確認すると、皮膚線量とＲ成分の変化量とが必ずしも対応しなかった。このため、放射線照射前に皮膚を撮影した皮膚画像から成分を分離して、放射線照射による皮膚反応後の皮膚画像を精度良く作成することは、非常に困難であった。
発明者らは、鋭意研究したものの、研究を開始してから約１０年の間、皮膚線量から皮膚反応を精度よく予測する方法は発見されなかった。
ここで、皮膚が赤く炎症する反応は、ダメージを受けた細胞を修復するために皮膚表面付近の毛細血管が拡張して血流量を増やそうとする生体反応の結果である。このため、発明者らは、血液中の色素であるヘモグロビン流量と、赤い炎症として現れる皮膚反応の程度が相関している可能性があることに着目した。また、ヘモグロビン流量は、治療に要する放射線量とも相関している可能性にも着目した。
そうすると、皮膚線量とヘモグロビン流量の相関性を臨床データから把握できれば、皮膚線量と皮膚反応の関係を定量的に把握できる可能性がある。
しかし、ヘモグロビン流量という尺度のみでは、医師間の意思疎通においても感覚的に理解し難く、目で見てわかる皮膚反応（軽い紅斑、強い紅斑、水泡・糜爛、潰瘍）、特に変化後の皮膚画像で示すことが望まれる。
一方、化粧効果や健康状態の定量的評価のために顔の肌色を構成する主成分であるメラニンとヘモグロビンの相対的色素含量をデジタル画像のＲＧＢ値から抽出する技術が開発された（特許文献：特開２００２−２０００５０、非特許文献：Ｔｓｕｍｕｒａ Ｎ， ｅｔ ａｌ．： Ｉｍａｇｅ−ｂａｓｅｄ ｓｋｉｎ ｃｏｌｏｒ ａｎｄ ｔｅｘｔｕｒｅ ａｎａｌｙｓｉｓ／ｓｙｎｔｈｅｓｉｓ ｂｙ ｅｘｔｒａｃｔｉｎｇ ｈｅｍｏｇｌｏｂｉｎ ａｎｄ ｍｅｌａｎｉｎ ｉｎｆｏｒｍａｔｉｏｎ ｉｎ ｔｈｅ ｓｋｉｎ， ＡＣＭ Ｔｒａｎｓ． Ｇｒａｐｈｉｃｓ ２２， ７７０， ２００３．）。この技術を適用すると、撮影した皮膚画像からヘモグロビン量だけを意図的に操作したような合成画像を得ることが出来、この操作を画像の中の特定の領域にだけ選択的に行うことも可能である。
しかし、放射線治療による皮膚反応は、グレード１（軽い紅斑）からグレード４（潰瘍）まであるために赤味だけで表現できるわけではなく、過去の研究で撮影画像のＲ成分の変化と皮膚線量との相関が取れなかったことからも、この技術を応用することには困難性があると考えられた。
発明者らは、鋭意研究した結果、放射線治療においては、グレード１〜グレード２の間で放射線治療計画を立てられることから、グレード１〜グレード２を重視すれば良いことを着想した。そして、グレード１〜グレード２であれば、ヘモグロビン流量と皮膚線量との間に相関がとれ、かつ、ヘモグロビン流量と皮膚反応との間にも相関がとれるのではないかということを着想した。
そして、発明者らは、皮膚を撮影したＲＧＢ表現形式の皮膚画像（反射光）を、ヘモグロビン成分等の生体要素（吸光成分）で表現する生体要素表現形式に変換し、そのうちのヘモグロビン成分に対して、照射予定の放射線による皮膚線量に対応するヘモグロビン流量変化量を適用し、生体要素表現形式からＲＧＢ表現形式に皮膚画像を戻すことで、皮膚線量に応じた皮膚反応を予測した予測皮膚画像の作成に成功した。
また、これとともに、精度の高いグレード判定を実現し、医師間の意思疎通の尺度を定量化するとともに細かくすることを実現した。
以下、この発明の一実施形態を図面と共に説明する。
図１は、皮膚変化予測システム１の構成を示すブロック図であり、図２は、放射線照射によるグレード決定の説明図であり、図３は、皮膚変化の説明図である。
図１に示すように、皮膚変化予測システム１は、撮影装置であるカメラ２と、パーソナルコンピュータである表示入力装置３と、パーソナルコンピュータである皮膚変化予測装置４とで構成されている。なお、図３の各画像は、元はカラー画像であるが、モノクロ画像となる特許図面上で違いをわかりやすくするために、実際のものよりも全て一定濃度で濃度を濃くしたものである。
カメラ２は、生体の皮膚を撮影した撮影画像データ（皮膚画像色成分データ）を表示入力装置３へ送信する。この撮影画像データは、図３（Ａ）に示すＲＧＢ表現形式（図示はモノクロで表示されるが、実際にはカラー画像である）の皮膚画像のデータである。ここで取得される撮影画像データは、撮影環境として撮影対象の皮膚に当たっている照明量がある程度の一様性を有していることが必要である。この一様性は厳密である必要はなく、不自然な不均一さがない程度でよい。
なお、撮影画像の色成分は、光の三原色（ＲＧＢ表現形式）であるＲ成分（赤）とＧ成分（緑）とＢ成分（青）、または、色の三原色（ＳＭＹ表現形式）であるＣ成分（シアン）とＭ成分（マゼンタ）とＹ成分（イエロー）、あるいは、印刷の四原色（ＣＭＹＫ表現形式）であるＣ成分（シアン）とＭ成分（マゼンタ）とＹ成分（イエロー）とＫ成分（ブラック）等、色を表現する適宜の成分で構成することができる。この実施例では、好ましい色成分として光の三原色を用いている。
表示入力装置３は、操作入力を受け付けるキーボードとマウスまたはタッチパネルなどの入力部と、画像表示を行う液晶ディスプレイまたはＣＲＴモニタ等の表示部と、データおよびプログラムを記憶するハードディスク等の記憶部と、前記プログラムに従って各種動作や演算を行うＣＰＵとＲＯＭとＲＡＭによる制御部と、カメラ２や皮膚変化予測装置４等の外部機器と接続されてデータ送受信を行うＵＳＢ等の接続インターフェースとを備えている。
表示入力装置３は、カメラ２から受信した撮影画像データを皮膚変化予測装置４に送信する処理、撮影画像と入力部で構成される入力画面を表示して皮膚線量、照射領域、およびテスト領域を入力させる処理、入力された皮膚線量、照射領域、およびテスト領域を皮膚変化予測装置４に送信する処理を実行する。
皮膚変化予測装置４は、少なくともプログラムおよびデータを記憶する記憶部と、前記プログラムに従って各種動作や演算を行うＣＰＵと、表示入力装置３等の外部機器を接続するＵＳＢ等の接続インターフェースとを備えている。記憶部には、記録媒体９から皮膚変化予測プログラムがインストールされている。なお、皮膚変化予測装置４は、表示入力装置３と同様に、入力部および表示部等をさらに備えても良い。この皮膚変化予測装置４は、画像処理を行う画像処理装置でもある。
皮膚変化予測装置４は、記憶部のプログラムに従って制御部が動作する処理機能部として、放射線情報受付部１１（要因情報受付部）、変化量対応データ１２、変化量決定部１３、グレード判定部１４、グレード対応データ１５、グレード出力部１６、照射領域受付部１７、照射前皮膚画像取得部２１（皮膚画像取得部）、皮膚領域抽出部２１１、表現形式変換部２２、画像変化部２３、表現形式復元部２４、予測皮膚画像出力部２５、テスト領域画像取得部３１、皮膚領域抽出部３１１、主成分分析部３２、独立成分分析部３３、色素ベクトル適切性判定部３４、色素ベクトル決定部３５、及び色素ベクトル標準データ３６を備えている。このうちの表現形式変換部２２と、変化量決定部１３と、画像変化部２３と、表現形式復元部２４は、放射線照射による皮膚変化を演算する変化演算部として機能する。
放射線情報受付部１１は、放射線治療計画における放射線情報（要因情報）の入力を表示入力装置３により受け付ける。この放射線情報の受付は、手入力によって入力される、ＵＳＢ等の情報媒体から入力される、あるいは治療計画装置から通信手段を介して受信するなど、適宜の方法によって実行するとよい。このとき受け付ける放射線情報は、少なくとも放射線種別と線量を必須とし、放射線種別によってはさらに線質を必須とする。具体的には、放射線種別がＸ線であれば、線量の入力を受け付ける。放射線種別が荷電粒子線（陽子、炭素を含む重イオン粒子）であれば、線量と線質の入力を受け付ける。線量は物理的に測定可能な吸収線量（単位 グレイ（Ｇｙ）、または臨床効果の指標として生物学的効果比（ＲＢＥ）を乗じた臨床線量（単位 Ｇｙ （ＲＢＥ））で表現され、線質は粒子の単位長さあたりのエネルギー付与量（ＬＥＴまたは線エネルギー、単位 ｋｅＶ／ｍｉｃｒｏｍｅｔｅｒ）で表現される。また、放射線情報には、放射線を照射するタイミングに関して、基準日からの経過時間（経過日数）も受け付ける。なお、本明細書において、「皮膚線量」の用語は、Ｘ線であれば線量を指すもの、荷電粒子線であれば線量と線質を指すものとして用いる。
変化量対応データ１２は、皮膚線量とヘモグロビン流量の変化量を対応付けたデータである。このデータは、皮膚線量とヘモグロビン流量の数値を対応付けた表データとする、あるいは、皮膚線量を入力するとヘモグロビン流量が算出される計算式データとするなど、適宜の形式のデータとすることができる。この変化量対応データ１２は、皮膚線量照射前後の照射領域のヘモグロビン量がレーザー血流計等の測定器によって測定され、皮膚線量と対応づけられた臨床データを用いて、事前に作成され登録されたものである。
ここで、本実施例における変化量対応データ１２を詳述する。
図２（Ａ）は、ヘモグロビン色素画像の画素値であるＨ（ｘ，ｙ）とヘモグロビン流量（皮膚血流量）の変化量の測定結果をプロットしたグラフである。全データに対して原点を通る直線でフィットすると、傾きは０．０１３９（２０）（ｍｌ／ｍｉｎ１００ｇ）−１であった。相関係数は０．７４であり、両者の間に相関があると言える。なお図中に点線で図示するエラーバンドに、照明環境に起因する１５％の画素値の誤差、血流量測定に起因する０．２ｍｌ／ｍｉｎ／１００ｇの誤差、そしてフィットによる誤差を二乗和して算出されている。血流量測定に起因する誤差とは主には心拍による生体的な効果が原因である。
図２（Ｂ）は、皮膚線量とヘモグロビン流量の関係のデータをプロットしたグラフである。この定量関係が皮膚障害における皮膚の紅斑の度合いを数値化する重要な関係となる。グラフ中の塗りつぶしシンボルはレーザー血流計によって直接測定されたデータであり、中抜きシンボルは図２（Ａ）の関係式を用いて画像解析した画素値から変換して得たヘモグロビン流量（血流量）である。その結果、０．２０５（１８）ｍｌ／ｍｉｎ／１００ｇ／Ｇｙ（ＲＢＥ）の傾きを得た。レーザーによるデータだけの場合の相関係数は０．７５、画像も含めた全データの場合は０．６０であり相関が認められる。ここで図中に点線で図示するエラーバンドは、血流量測定に起因する０．２ｍｌ／ｍｉｎ／１００ｇの誤差、そしてフィットによる誤差を二乗和して算出されている。
従って、本実施例においては、ヘモグロビン色素画像の画素値に対するヘモグロビン流量（皮膚血流量）の変化量の傾き０．０１３９（２０）（ｍｌ／ｍｉｎ１００ｇ）−１、および、ヘモグロビン流量に対する皮膚線量の傾き０．２０５（１８）ｍｌ／ｍｉｎ／１００ｇ／Ｇｙ（ＲＢＥ）を変化量対応データ１２としている。
このように、ヘモグロビン色素画像の画素値とヘモグロビン流量（皮膚血流量）に相関関係があり、かつ、ヘモグロビン流量（皮膚血流量）と皮膚線量に相関関係があることから、ヘモグロビン色素画像における画素値変化量から皮膚線量を推定する、逆に皮膚線量からヘモグロビン色素画像における画素値変化量を推定することができる。
変化量決定部１３は、放射線情報受付部１１から取得する放射線情報と変化量対応データ１２に基づいて、ヘモグロビン流量を変化させる変化量を決定する。決定した変化量は、画像変化部２３へ送信する。
グレード判定部１４は、放射線情報受付部１１から変化量決定部１３を介して受け取る放射線情報と、グレード対応データ１５に基づいて、グレードを判定する。この時判定するグレードは、例えば従来の４段階のグレードをさらに１０倍の精細度に高めた４０段階のグレードとするなど、精度よく詳細なグレードとする。グレード判定部１４は、判定したグレードをグレード出力部１６へ送信する。
グレード対応データ１５は、放射線情報とグレードを対応付けるためのデータであり、表形式のデータとする、あるいは、放射線情報の数値を代入することで解を求める演算式とするなど、適宜の形式のデータとすることができる。
このグレード対応データ１５は、臨床データによって事前に作成されたものである。詳述すると、図２（Ｃ）及び図２（Ｄ）は、画素値変化量から算出できるヘモグロビン流量（血流量）とグレード判定の関係を示すグラフであり、第一縦軸（図示左側の縦軸）がヘモグロビン流量（血流量）、横軸がグレードを示している。図２（Ｃ）における第二縦軸（図示右側の縦軸）は画素値変化量、図２（Ｄ）における第二縦軸（図示右側の縦軸）は皮膚線量を示している。
このグラフより、本実施例では、グレード対応データ１５として、ヘモグロビン流量（血流量）が６．４ｍｌ／ｍｉｎ／１００ｇ以上増加する場合にグレードは２以上と定めた。また、本実施例では、ヘモグロビン流量変化（血流量変化）との変換係数を使うと、グレード１／２の境界線に対応する画素値変化量、皮膚線量が算出でき、それぞれ０．０９，３２Ｇｙ（ＲＢＥ）程度と推測される。この場合の誤差は、＋−２０％である。
なお、グレード判定の閾値は上記数値に限られるものではなく、適宜定めることができる。重要なことは、グラフに示すように、ヘモグロビン流量（血流量）、皮膚線量、画素値変化量のそれぞれが精密なグレード値と１：１で定量的に対応付けられていることである。
グレード出力部１６は、グレード判定部１４で得られたグレードを表示入力装置３へ送信し、表示入力装置３にグレード表示画面として表示する、
照射領域受付部１７は、生体の皮膚のうち放射線治療計画において放射線を照射する照射領域の入力を表示入力装置３から受け付ける。照射領域受付部１７は、受け付けた照射領域を画像変化部２３へ送信する。
照射前皮膚画像取得部２１は、放射線照射前に皮膚を撮影した照射前皮膚画像データを表示入力装置３から取得する。この照射前皮膚画像データは、ＲＧＢ表現形式であるため、Ｒ成分とＧ成分とＢ成分に容易に分離できる。照射前皮膚画像取得部２１は、この分離可能な照射前皮膚画像データを表現形式変換部２２へ送信する。
皮膚領域抽出部２１１は、照射前皮膚画像取得部２１で取得した照射前皮膚画像データから皮膚領域（肌領域）を抽出する。この皮膚領域の抽出は、照射前皮膚画像データの各ピクセルのうち、皮膚領域条件を満たすピクセルを抽出することで実行する。
皮膚領域条件は、ＲＧＢ表現形式における次の条件を全て満たすピクセルを皮膚領域とする。
＜皮膚領域条件＞
Ｒ＞Ｒ＿ｔｈ
Ｇ／Ｒ＞Ｒ＿ｒａｔｉｏ１
Ｂ／Ｒ＞Ｒ＿ｒａｔｉｏ２
（注）ＲはＲ成分の濃度（例えば０〜２５５）を示す。
ＧはＧ成分の濃度（例えば０〜２５５）を示す。
ＢはＢ成分の濃度（例えば０〜２５５）を示す。
Ｒ＿ｔｈは所定値であり、例えば１００〜１６０とする、あるいは１０とすることができる。なお、Ｒ＿ｔｈは、１０以上とすることができ、５０以上とすることが好ましく、１００以上とすることがより好ましい。また、Ｒ＿ｔｈは、２００以下とすることが好ましく、１６０以下とすることがより好ましい。
Ｒ＿ｒａｔｉｏ１は所定値であり、例えば１．１とする、あるいは１．０とする等、適宜の値とすることができる。なお、Ｒ＿ｒａｔｉｏ１は、０．５〜１．５とすることが好ましく、０．９〜１．２とすることがより好ましく、１．０〜１．１とすることがさらに好ましい。
Ｒ＿ｒａｔｉｏ２は所定値であり、例えば１．１とする、あるいは１．０とする等、適宜の値とすることができる。なお、Ｒ＿ｒａｔｉｏ２は、０．５〜１．５とすることが好ましく、０．９〜１．２とすることがより好ましく、１．０〜１．１とすることがさらに好ましい。
なお、上記皮膚領域の抽出条件により適切な画像を得るための撮影条件は、（１）肌色が認識できる程度の明るさがあること、（２）不自然な影や白とびが生じることなく皮膚領域が照らされていること、の２点のみである。この２点の条件をクリアできていれば、上記皮膚領域の抽出が良好に行われる。なお、条件（２）においては、皮膚領域が略均一に照らされていることがより好ましい。
表現形式変換部２２は、取得した照射前皮膚画像データ（皮膚領域抽出済み）を、ＲＧＢ表現形式から色素ベクトル表現形式へ変換する。この変換は、色素ベクトル決定部３５から受け取る色素ベクトルを用いることで実行する。ここで、色素ベクトルとは、本実施形態においてヘモグロビンベクトルとメラニンベクトルを指す。色素ベクトル表現形式は、ヘモグロビン成分と、メラニン成分と、その他成分の３成分に情報が分離されている。表現形式変換部２２による表現形式の変更の詳細については後述する。
画像変化部２３は、照射前皮膚画像データから表現形式変換部２２によって抽出されたヘモグロビン成分のうち、照射領域受付部１７によって定められた照射領域に対して、変化量決定部１３で決定した変化量を適用し、当該照射領域のヘモグロビン成分を変化させる。この画像変化部２３による画像変化処理の詳細は後述する。
表現形式復元部２４は、変化量決定部１３で変化させた後の色素ベクトル表現形式の画像データ（皮膚画像生体要素成分データ）を、ＲＧＢ表現形式に変換することで表現形式の復元を実行する。表現形式復元部２４は、復元後のＲＧＢ表現形式の画像データを予測皮膚画像として予測皮膚画像出力部２５へ送信する。
予測皮膚画像出力部２５は、受信した予測皮膚画像を表示入力装置３へ送信し、表示させる。
テスト領域画像取得部３１は、表示入力装置３によって入力されたテスト領域を取得する。このテスト領域は、撮影した皮膚画像の一部を選択した領域であり、例えば１００ピクセル×１００ピクセルの小領域とするなど、適宜の領域とすることができる。
皮膚領域抽出部３１１は、テスト領域の皮膚画像に対して皮膚領域抽出部２１１と同一の操作を実施して肌画像を抽出する。なお、この皮膚領域抽出部３１１を設けず、皮膚領域抽出部２１１で抽出された皮膚領域からテスト領域画像取得部３１がテスト領域を抽出する構成としてもよい。
主成分分析部３２は、テスト領域の皮膚画像（皮膚領域抽出済み）に対して主成分分析（多変量解析の一種）を行う。主成分分析部３２は、最も相関の高い方向（最も分布の広い方向）へ主成分の座標軸（第一成分の第一軸）を取り出し、以降、それまでの座標軸に直交する平面（多次元上での平面）の上で最も相関の高い方向（最も分布の広い方向）へ第ｎ成分の第ｎ軸を取り出す作業を繰り返す。演算した結果は、最も相関のある高い第一成分（第一軸）と、この第一成分に直交する第二成分（第二軸）の二軸により、９８％以上を表現できている。このため、残りの成分をまとめて第三成分として破棄し、以降の処理で第一成分と第二成分のみを用いる。これにより、皮膚変化予測装置４は、高精度かつ高速の処理を実現している。テスト領域画像での画像データを主成分分析したベクトルで表現すると、次の（数１）に示す式となる。
Figure 0006448142
 （注）（ｒ０，ｇ０，ｂ０）は、肌色のベース色を示す。
ｒａｖｅ、ｇａｖｅ、ｂａｖｅは、画像中の（ｒ，ｇ，ｂ）の値から（ｒ０，ｇ０，ｂ０）の値を減算した後の値の平均値を示す。
ｐｃ１は第一成分、ｐｃ２は第二成分、ｐｃ３は第三成分を示す。
（ｒ１，ｇ１，ｂ１）は、第一成分を表すベクトルである。
（ｒ２，ｇ２，ｂ２）は、第二成分を表すベクトルである。
（ｒ３，ｇ３，ｂ３）は、第三成分を表すベクトルである。
なお、各ベクトルは１に規格化されており、互いに直交しているためにお互いの内積はゼロとなる。
独立成分分析部３３は、主成分分析部３２で得られた第一成分と第二成分で表現されるデータに対して独立成分分析（多変量解析の一種）を実行し、各軸が最も独立な状態となる直交しない二軸の座標軸を取る。得られた座標軸の各成分は、一方の成分がヘモグロビンベクトルとなり、他方の成分がメラニンベクトルとなる。テスト領域画像での画像データを独立成分分析したベクトルで表現すると、次の（数２）に示す式となる。
Figure 0006448142
 （注）ｒａｖｅ、ｇａｖｅ、ｂａｖｅは、画像中の（ｒ，ｇ，ｂ）の値から（ｒ０，ｇ０，ｂ０）の値を減算した後の値の平均値を示す。
（ｒｈ，ｇｈ，ｂｈ）は、ヘモグロビンベクトルを表す。
（ｒｍ，ｇｍ，ｂｍ）は、メラニンベクトルを表す。
ｈは、ヘモグロビン量（ヘモグロビン座標）を示す係数である。
ｍは、メラニン量（メラニン座標）を示す係数である。
なお、等号ではなく近似等号を使っているのは、第三成分を無視したためである。また、各色素ベクトルは１に規格化されているが、色素ベクトル同士に直交関係はない。
また、上記（数１）および（数２）では平均値ｒａｖｅ、ｇａｖｅ、ｂａｖｅを各成分に加算することを明確に表記しているが、平均値ｒａｖｅ、ｇａｖｅ、ｂａｖｅを（数１）及び（数２）の数式におけるＲ０，Ｂ０，Ｇ０に含ませて表現上は見えない形に表現してもよい。この場合でも、同じ意味の数式であるから、ＲＧＢ表現形式をヘモグロビンベクトル等が用いられた生体要素表現形式に変換することができる。
色素ベクトル適切性判定部３４は、得られた色素ベクトルが適切か否かを判定する。詳述すると、テスト領域画像取得部３１で取得されるテスト領域は、主成分分析部３２でヘモグロビンとメラニンを第一成分と第二成分に選び、独立成分分析部３３で独立成分として分離してくれるものでなければならない。このようなテスト領域が操作者によって表示入力装置３で入力されたか否かを判定するために、色素ベクトル適切性判定部３４は、得られた色素ベクトルがヘモグロビンベクトルとメラニンベクトルであるか否かを判定する。
色素ベクトル決定部３５は、色素ベクトル適切性判定部３４による判定結果が「適切」であれば、独立成分分析部３３から送信されたヘモグロビンベクトルとメラニンベクトルを色素ベクトルとして決定する。色素ベクトル適切性判定部３４による判定結果が「不適切」であれば、色素ベクトル決定部３５は、色素ベクトル標準データ３６から標準的なヘモグロビンベクトルとメラニンベクトルを取得して、これを色素ベクトルとして決定する。色素ベクトル決定部３５は、決定したヘモグロビンベクトルとメラニンベクトルを表現形式変換部２２へ送信する。なお、この実施例では、一度不適切であれば色素ベクトル標準データ３６を使用することとしているが、これに限らず、テスト領域の取得から適否の判定までを複数回繰り返し、所定回数繰り返しても全て不適切であった場合に色素ベクトル標準データ３６を用いる構成としてもよい。
色素ベクトル標準データ３６は、予め算出された標準的な色素ベクトル（ヘモグロビンベクトルとメラニンベクトル）を記憶している。
＜表現形式変換部２２による表現形式の変更の詳細＞
表現形式変換部２２は、独立成分分析部３３で得られた色素ベクトルを用いて、（ｒ，ｇ，ｂ）の三成分を、次の（数３）の式に示す（ｈ，ｍ，ｓ）の三成分に線形変換（座標変換）する。
Figure 0006448142
 （注）ｓは、照明の単位ベクトル（陰影ベクトル）の係数であり、ヘモグロビンおよびメラニンによる色素ベクトルにて表現できない情報をすべて押し込んだものに該当する。
ここで、（数３）への線形変換について説明する。（数３）をより正確に表現すると、次の（数４）で表すことができる。
Figure 0006448142
 （注）（ｘ，ｙ）の二次元の位置関数になっている変数は画像データのピクセル毎に異なる値が入る。
（ｘ、ｙ）の関数になっていない箇所は定数であることを意味している。
上記（数４）は、（ｒ（ｘ，ｙ），ｇ（ｘ，ｙ）ｂ（ｘ，ｙ））で表現されていたフルカラー画像が（Ｈ（ｘ，ｙ），Ｍ（ｘ，ｙ），Ｓ（ｘ，ｙ））の三種類の色素画像に分解されることを示している。
ここで、Ｈ（ｘ，ｙ），Ｍ（ｘ，ｙ），Ｓ（ｘ，ｙ）は、次の（数５）で表すことができる。
Figure 0006448142
 （注）ｈａｖｅ，ｍａｖｅ，ｓａｖｅは、ｈ，ｍ，ｓの各成分の平均値を示す。この平均値は、画像全体の平均値とする等、適宜の平均値とできるが、皮膚を表現する領域だけのＲＧＢ成分から算出する平均値であることが好ましい。従って、この実施例では皮膚領域のＲＧＢの平均値を用いる。
上記（数５）におけるｈ（ｘ，ｙ），ｍ（ｘ，ｙ），ｓ（ｘ，ｙ）は、上述した（数１）により示した主成分分析を行って得た第一成分と第二成分の主成分ベクトルを用いて、次の（数６）の計算式により導出する。
Figure 0006448142
ｈａｖｅ，ｍａｖｅ，ｓａｖｅについては、上記（数６）におけるｒ（ｘ，ｙ），ｇ（ｘ，ｙ），ｂ（ｘ，ｙ），ｈ（ｘ，ｙ）の箇所をｒａｖｅ，ｇａｖｅ，ｂａｖｅ，ｈａｖｅにそれぞれ置き換えた式により導出できる。
この表現形式の変換により、図３（Ａ）に示した皮膚画像が、図３（Ｂ）に示すヘモグロビン成分皮膚画像、図３（Ｃ）に示すメラニン成分皮膚画像、図３（Ｄ）に示すその他成分皮膚画像に分離される。図はモノクロ画像となっているが、実際には、図３（Ｂ）は薄い赤色系の単一色スケール画像であり、図３（Ｃ）は薄い黄土色系の単一色スケール画像であり、図３（Ｄ）は薄い灰色系の単一色スケール画像（グレースケール画像）である。例えば青い血管は、ヘモグロビンまたはメラニンの色素ベクトルでは表現できないため、その影響が影となってその他成分である図３（Ｄ）に見えている。図３（Ｂ）〜図３（Ｄ）に図示する例は、皮膚領域抽出を行わないで変換した場合の画像であるが、例えば、背景など皮膚以外が写りこんでいる画像の場合は、平均値ｒａｖｅ、ｇａｖｅ、ｂａｖｅに対して皮膚領域抽出を行うとさらに良好なヘモグロビン成分皮膚画像が得られる（後述の図７参照）。
なお、この実施例では、生体要素成分としてヘモグロビン成分とメラニン成分による色素成分を用いているが、これに限らず、他の成分を用いる構成としてもよい。この場合、ヘモグロビン成分、ユーメラニン成分、フェオメラニン成分、カロチン成分、またはこれらの複数等とすることができ、このうちのユーメラニン成分とフェオメラニン成分をまとめてメラニン成分とすることができる。どのように構成する場合でも、生体要素成分は、少なくともヘモグロビン成分を有することが好ましい。
＜画像変化部２３による画像変化処理の詳細＞
画像変化部２３は、上記（数３）の式におけるｈ成分に対して、変化量決定部１３で決定された変化量（ヘモグロビン増加量）を増加させる。
このようにして画像変化させた後に各成分を合成し元のＲＧＢ表現形式に復元すると、図３（Ｅ）に示す予測皮膚画像が得られる。図３（Ｅ）の予測皮膚画像は、特許図面としてモノクロ画像になっているが、実際にはカラー画像である。この図３（Ｅ）の例は、四角形の第１照射領域４１と、この第１照射領域４１と一部重なる四角形の第２照射領域４２にそれぞれ放射線を照射した場合の予測皮膚画像を示している。この第１照射領域４１と第２照射領域４２は、照射領域受付部１７で受け付けた領域である。第１照射領域４１は、第２照射領域４２よりも皮膚線量が多いために第２照射領域４２よりも濃く（赤く）表示されており、かつ、第１照射領域４１と第２照射領域４２が重なっている重複照射領域４３はさらに濃く（赤く）表示されている。
なお、表示入力装置３と皮膚変化予測装置４は、別の装置として構成したが、１つの装置としてまとめて構成してもよい。この場合、１つのコンピュータとして構成することができる。例えばタッチパネルモニタなどの表示入力部、あるいはキーボードやマウス等の入力部と液晶ディスプレイやＣＲＴモニタ等の表示部を、表示入力装置３とすることができる。また、例えば、制御部と記憶部を皮膚変化予測装置４とすることができる。これにより１つのコンピュータで表示入力装置３と皮膚変化予測装置４を構成できる。
図４は、グレード出力部１６が表示入力装置３に出力して表示するグレード判定画面５０の画面構成図である。
グレード判定画面５０は、予測グレード値表示部５１、予測グレード目盛表示部５２、急性期スコア表示部５６、および皮膚反応表示部５８を備えている。
予測グレード値表示部５１は、グレード判定部１４によって判定されたグレード値が表示される。このグレード値は、従来の１〜４の４段階よりもさらに細かい段階を示せる値であり、例えば０．１〜４．０までの４０段階で表示される。
なお、このグレード値の細かい段階は、どこからどこまでが従来の４段階のどこに入るかが明確な段階とすることが好ましく、この実施例では、本発明によるグレード（以下新グレードという）が０．１〜１．０までは従来グレードの１に対応し、新グレードの１．１〜２．０までが従来グレードの２に対応し、新グレードの２．１〜３．０までが従来グレードの３に対応し、新グレードの３．１〜４．０までが従来グレードの４に対応している。
予測グレード目盛表示部５２は、従来グレードを新グレードに対応させて細かく分割した目盛を示し、この目盛を指す予測値マーク５３を表示する。
予測値マーク５３は、予測値を正確に示す予測値明示部５４と、この予測値から誤差が生じ得る範囲を示す誤差範囲表示部５５とを備えている。この誤差範囲表示部５５が予測グレード目盛表示部５２の目盛の並ぶ報告に幅をもっており、その内側で重なるように予測値マーク５３が設けられていることで、予測値の中央値と誤差の範囲を直観的にわかりやすく認識できる。誤差範囲表示部５５で表示される誤差は、予め定めた一定の範囲とする、あるいは、その都度演算で求めるなど、適宜の方法によって定めることができる。予め定めた一定の範囲とする場合は、例えば、レーザー血流計によって測定するヘモグロビン量の誤差（例えば１０％〜１５％）に撮影環境に依存する誤差（例えば１５％）を加味した誤差（例えば１８％〜２１％）とすることができる。
急性期スコア表示部５６は、従来のグレード１〜４を示しており、皮膚反応表示部５８は、従来のグレード１〜４に対応する皮膚反応である「軽い紅斑」「強い紅斑」「水泡、糜爛（ビラン）」「潰瘍」をそれぞれグレードに対応させて表示している。
急性期スコア表示部５６と皮膚反応表示部５８は、予測される新グレードが位置するグレード部分を、色付きで表示する等の差別化表示５８を行う。これにより、従来グレードであればどのグレードに該当するかを容易に把握できるようにしている。
図５は、これまでに説明した皮膚変化予測システム１に検証装置６を追加した皮膚変化検証システム１０の構成を示すブロック図である。
検証装置６は、予測皮膚画像取得部６１と、照射後皮膚画像取得部６２と、比較検証部６３と、検証結果出力部６４とを備えている。検証装置６は、少なくともプログラムおよびデータを記憶する記憶部と、前記プログラムに従って各種動作や演算を行うＣＰＵと、皮膚変化予測装置４等の外部機器を接続するＵＳＢ等の接続インターフェースとを備えている。検証装置６の記憶部には、記録媒体９から検証プログラムがインストールされている。この検証装置６は、皮膚変化予測装置４と別のハードウェアとしても良いが、皮膚変化予測装置４と同じハードウェアに機能を搭載して構成してもよい。
予測皮膚画像取得部６１は、皮膚変化予測装置４の予測皮膚画像出力部２５（図１参照）から予測皮膚画像を取得する。
照射後皮膚画像取得部６２は、放射線照射後にカメラ２で皮膚が撮影された照射後皮膚画像を表示入力部３を介して取得する。
比較検証部６３は、予測皮膚画像取得部６１から受信する予測皮膚画像と、照射後皮膚画像取得部６２から受信する照射後皮膚画像を比較し、放射線照射計画に対する照射領域および皮膚変化の適切性を検証する。具体的には、放射線照射によって皮膚反応の生じている領域と放射線照射計画における照射領域の一致度を算出し、この一致度が所定範囲内であるか否かによって適否を判定する。ここで、予測皮膚画像取得部６１から受信する予測皮膚画像は、別途の治療計画装置で計画された放射線情報等に基づいて予測された画像である。このため、この時点で治療計画に関する放射線情報は無くとも良いが、放射線情報を得てより詳細に比較検証する構成としてもよい。
また、皮膚反応によって回りの皮膚と色が異なっている程度を放射線治療計画によって色が変化する程度との一致度を算出し、この一致度が所定範囲内であるか否かによって適否を判定する。この色変化の一致度の算出は、予測皮膚画像と照射後皮膚画像を直接比較して行う、あるいは、上述した表現形式変換部２２による色素ベクトル表現への変換をしてヘモグロビン成分のみを比較して行うなど、適宜の構成とすることができる。
検証結果出力部６４は、比較検証部６３で検証した結果である適否判定と一致度を表示入力装置３に出力して表示させる。このとき、検証結果出力部６４は、次に図６と共に説明する検証用画像も表示させる。なお、図６の各画像は、元はカラー画像であるが、モノクロ画像となる特許図面上で違いをわかりやすくするために、実際のものよりも全て一定濃度で濃度を濃くしたものである。
図６は、照射領域の検証用画像を示す説明図である。図６（Ａ）は、治療計画から算出した皮膚線量の分布を示す皮膚線量分布画像である。図示の例では、１門１２．５Ｇｙ（ＲＢＥ）の処方線量を二方向から照射しており、二門照射が重なった部分の皮膚線量２２Ｇｙ（ＲＢＥ）を１００％としている。図６（Ｂ）は、図６（Ａ）に示した治療計画による照射後に写真撮影した照射後皮膚画像である。この照射後皮膚画像は、図６（Ａ）における長方形の白枠部分を拡大表示している。図６（Ｃ）は、照射後皮膚画像を色素分解したヘモグロビン色素画像である。図６（Ｄ）は、図６（Ａ）〜図６（Ｃ）における照射領域を説明するための説明図である。図示するように、本実施例の照射では、第一照射領域７１と、第二照射領域７５と、第一照射領域７１および第二照射領域７５の重なり部分である重複照射領域７３が存在している。このように、図６（Ａ）〜図６（Ｃ）に示す各画像を検証結果出力部６４により出力することで、図６（Ａ）の治療計画における皮膚線量分布画像と図６（Ｃ）の治療後のヘモグロビン色素画像を比較して計画通りに適切な照射が出来たか否かを検証できる。
なお、比較検証部６３および検証結果出力部６４は、治療計画における皮膚線量分布画像（図６（Ａ）参照）から推定した予測皮膚画像と照射後皮膚画像（図６（Ｃ）のヘモグロビン色素画像）を並べて表示入力装置３に表示し、医師の目視判断によって最終決定されて表示入力装置３に入力された適否を受け取る構成としてもよい。この場合、比較検証部６３で検証した結果を表示入力装置３に表示した上で最終判断を医師にゆだねる構成とすることができるが、比較検証部６３では比較検証を行わずに適否判断を完全に医師にゆだねる構成とすることも可能である。
以上の構成及び動作により、放射線照射による皮膚反応を事前に精度よく予測することができる。
また、皮膚変化予測装置４は、放射線照射後の皮膚の変化を予測皮膚画像で明瞭に表示することができる。このため、操作者は、放射線治療を行う前に、どの程度の皮膚変化が生じるかを目視で分かり易く把握することができる。
また、皮膚変化予測装置４は、従来よりも詳細なグレード値を出力することができる。このため、操作者は、例えばグレード２であっても、グレード１に近いグレード２なのか、それともグレード３に近いグレード２なのかを、明瞭に把握できる。従って、操作者が予測されたグレード値を見て皮膚線量を治療前に最適に調節することを容易に実現できる。
また、皮膚変化予測装置４は、精度の良いグレード値を定量的に出力できる。このため、グレード値は、医師間の意思疎通の共通尺度として利用されることができる。そして、皮膚変化予測装置４は、放射線照射による皮膚変化を示す尺度としてグレード値の標準化を実現し得る。
また、皮膚変化予測装置４は、グレード判定画面５０に従来グレードと新グレードとを並べて表示するため、従来グレードに慣れている医師が新グレードを違和感なく利用することができる。特に、予測値マーク５３を表示することにより、皮膚変化予測装置４は、医師が直観的にわかり易い形で表示することができる。
また、皮膚変化予測装置４は、放射線照射による皮膚の変化を、元の皮膚画像に対して相対量の変化として予測皮膚画像を作成する。このため、操作者は、撮影環境等の厳密性が要求されず、皮膚変化予測装置４を簡易に利用することができる。
また、皮膚変化予測装置４は、偏光照明などの特殊な撮影環境を必要とせず、照明（若しくは自然光）に不自然な不均一さがなければ十分である。このため、様々な環境で様々なカメラ２を用いて簡易に撮影した皮膚画像を利用することができる。
また、検証装置６は、皮膚変化予測装置４により予測された予測皮膚画像と、実際に放射線が照射された後の照射後皮膚画像とを比較検証することができる。これにより、操作者は、計画通りの放射線治療が行われたか否かを容易に確認することができる。また、予測皮膚画像と照射後皮膚画像の誤差が大きいような場合には、その原因を追究して予測精度を高めることに活用することができる。
また、皮膚領域抽出部２１１により皮膚領域を抽出してから計算を行うことで、図７の説明図に示すように色素分解の精度を向上させることができる。すなわち、図７（Ａ）は照射後に撮影した撮影画像であり、ここから皮膚領域を抽出することなくヘモグロビン色素画像を求めると図７（Ｂ）に示すヘモグロビン色素画像となる。これに対して皮膚領域を抽出して計算することで、図７（Ｃ）に示す撮影画像から図７（Ｄ）に示す高精度なヘモグロビン色素画像を得ることができる。なお、図７（Ａ）と図７（Ｃ）は同一の撮影画像である。また、図７の各画像は、元はカラー画像であるが、モノクロ画像となる特許図面上で違いをわかりやすくするために、実際のものよりも全て一定濃度で濃度を濃くしたものである。
また、ＲＧＢ表現形式と色素ベクトル表現形式との変換に、ＲＧＢ表現形式の画素の平均値を用いることで、変換を実現できる。
また、この平均値を、画像全体の平均値ではなく皮膚領域のみのＲＧＢ画素値から算出した平均値とすることで、照明環境や撮影環境の影響を受けずに精度よくヘモグロビン量を推定することができる。すなわち、背景色や衣服の色による影響を防止することができ、ゆるやかな撮影条件で誰でも容易に利用することができる。また、過去に撮影した撮影画像で解析を意図していない撮影画像であっても、精度よく解析してヘモグロビン量を推定することができる。
この発明は、上述の実施形態の構成のみに限定されるものではなく、多くの実施の形態を得ることができる。
例えば、色素ベクトルとしてヘモグロビンベクトルとメラニンベクトルを用いたが、メラニン成分をその他成分としてまとめてヘモグロビンベクトルとその他としてもよい。この場合でも、ヘモグロビン成分のみに画像変化を行えば良いため、皮膚線量の照射による皮膚変化の予測画像を適切に作成することができる。
また、放射線の影響による皮膚の変化に限らず、他の要因によっておこる皮膚紅斑の症状の定量化、および皮膚紅斑の予測に用いても良い。この場合、予測に使うことに加えて、現状の定量化による正確な把握、および医師間の正確な情報伝達に役立てることができる。
この発明は、放射線の照射による皮膚への影響を予測する産業、予測の適否を検証する産業、皮膚紅斑の影響を定量化する産業に利用することができる。
４…皮膚変化予測装置
６…検証装置
１１…放射線情報受付部
１３…変化量決定部
１４…グレード判定部
２１…照射前皮膚画像取得部
２２…表現形式変換部
２３…画像変化部
２４…表現形式復元部
２５…予測皮膚画像出力部
６１…予測皮膚画像取得部
６２…照射後皮膚画像取得部
６３…比較検証部
２１１，３１１…皮膚領域抽出部

Claims (6)

  1. 照射予定の放射線の放射線情報の入力を受け付ける放射線情報受付部と、
    生体の皮膚を撮影した皮膚画像を取得する皮膚画像取得部と、
    前記放射線情報により定まる放射線を照射したことによる前記皮膚の変化を演算して前記皮膚画像から変化後の予測皮膚画像を得る変化演算部と、
    前記予測皮膚画像を出力する出力部とを備え、
    前記変化演算部は、
    前記皮膚画像を、元の色成分による表現形式である皮膚画像色成分データから生体要素成分による表現形式である皮膚画像生体要素成分データに変換する表現形式変換部と、
    一部の前記皮膚画像生体要素成分データが前記放射線情報で定められる放射線照射によって変化する変化量を決定する変化量決定部と、
    前記変化量に応じて前記一部の皮膚画像生体要素成分データを変化させる変化部と、
    変化後の画像を元のカラー画像の表現形式に復元する表現形式復元部とを有する
    放射線照射による皮膚変化予測装置。
  2. 前記変化量決定部により決定した変化量に基づいて、従来グレードの区分単位である４段階よりも細かい新グレードでのグレード値を算出するグレード値算出部を備えた
    請求項１記載の放射線照射による皮膚変化予測装置。
  3. 照射予定の放射線の放射線情報の入力を受け付ける放射線情報受付部と、生体の皮膚を撮影した皮膚画像を取得する皮膚画像取得部と、前記放射線情報により定まる放射線を照射したことによる前記皮膚の変化を演算して前記皮膚画像から変化後の予測皮膚画像を得る変化演算部と、前記予測皮膚画像を出力する出力部とを備えた放射線照射による皮膚変化予測装置から前記予測皮膚画像を取得する予測皮膚画像取得部と、
    放射線照射後に皮膚が撮影された照射後皮膚画像を取得する照射後皮膚画像取得部と、
    前記予測皮膚画像と前記照射後皮膚画像とに基づいて放射線照射が適切に行われたか否かを検証する検証部とを備えた
    検証装置。
  4. コンピュータを、
    照射予定の放射線の放射線情報の入力を受け付ける放射線情報受付部と、
    生体の皮膚を撮影した皮膚画像を取得する皮膚画像取得部と、
    前記放射線情報により定まる放射線を照射したことによる前記皮膚の変化を演算して前記皮膚画像から変化後の予測皮膚画像を得る変化演算部と、
    前記予測皮膚画像を出力する出力部として機能させ、
    前記変化演算部は、
    前記皮膚画像を、元の色成分による表現形式である皮膚画像色成分データから生体要素成分による表現形式である皮膚画像生体要素成分データに変換する表現形式変換部と、
    一部の前記皮膚画像生体要素成分データが前記放射線情報で定められる放射線照射によって変化する変化量を決定する変化量決定部と、
    前記変化量に応じて前記一部の皮膚画像生体要素成分データを変化させる変化部と、
    変化後の画像を元のカラー画像の表現形式に復元する表現形式復元部として機能する
    放射線照射による皮膚変化予測プログラム。
  5. コンピュータを、
    照射予定の放射線の放射線情報の入力を受け付ける放射線情報受付部と、生体の皮膚を撮影した皮膚画像を取得する皮膚画像取得部と、前記放射線情報により定まる放射線を照射したことによる前記皮膚の変化を演算して前記皮膚画像から変化後の予測皮膚画像を得る変化演算部と、前記予測皮膚画像を出力する出力部とを備えた放射線照射による皮膚変化予測装置から前記予測皮膚画像を取得する予測皮膚画像取得部と、
    放射線照射後に皮膚が撮影された照射後皮膚画像を取得する照射後皮膚画像取得部と、
    前記予測皮膚画像と前記照射後皮膚画像とに基づいて放射線照射が適切に行われたか否かを検証する検証部として機能させる
    検証プログラム。
  6. 照射予定の放射線の放射線情報の入力を受け付ける放射線情報受付部により受け付け、
    前記生体の皮膚を撮影した皮膚画像を皮膚画像取得部により取得し、
    前記放射線情報により定まる放射線を照射したことによる前記皮膚の変化を変化演算部により演算して前記皮膚画像から変化後の予測皮膚画像を得、
    前記予測皮膚画像を出力部により出力する皮膚変化予測方法であって、
    前記変化演算部は、
    前記皮膚画像を、元の色成分による表現形式である皮膚画像色成分データから生体要素成分による表現形式である皮膚画像生体要素成分データに表現形式変換部により変換し、
    一部の前記皮膚画像生体要素成分データが前記放射線情報で定められる放射線照射によって変化する変化量を変化量決定部により決定し、
    前記変化量に応じて前記一部の皮膚画像生体要素成分データを変化部により変化させ、
    変化後の画像を元のカラー画像の表現形式に表現形式復元部により復元する
    皮膚変化予測方法。
JP2016510041A 2014-03-28 2015-03-25 放射線照射による皮膚変化予測装置と検証装置 Expired - Fee Related JP6448142B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2014068506 2014-03-28
JP2014068506 2014-03-28
PCT/JP2015/001692 WO2015146164A1 (ja) 2014-03-28 2015-03-25 放射線照射による皮膚変化予測装置と検証装置

Publications (2)

Publication Number Publication Date
JPWO2015146164A1 JPWO2015146164A1 (ja) 2017-04-13
JP6448142B2 true JP6448142B2 (ja) 2019-01-09

Family

ID=54194708

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016510041A Expired - Fee Related JP6448142B2 (ja) 2014-03-28 2015-03-25 放射線照射による皮膚変化予測装置と検証装置

Country Status (4)

Country Link
US (1) US10315052B2 (ja)
EP (1) EP3124077B1 (ja)
JP (1) JP6448142B2 (ja)
WO (1) WO2015146164A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111034153B (zh) * 2017-07-31 2022-12-23 株式会社理光 通信***，分散处理***，分布式处理方法和记录介质
WO2019078327A1 (ja) * 2017-10-20 2019-04-25 パナソニック株式会社 炎症の評価システム、評価方法、プログラム、及び非一時的記録媒体
EP4151276A1 (en) * 2021-09-20 2023-03-22 Koninklijke Philips N.V. Monitoring and detection of cutaneous reactions caused by radiotherapy
CN114209288A (zh) * 2022-01-14 2022-03-22 平安普惠企业管理有限公司 皮肤状态预测方法、皮肤状态预测装置、设备及存储介质

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2891350B2 (ja) * 1997-03-04 1999-05-17 技術研究組合医療福祉機器研究所 三次元画像処理方法
JP3877959B2 (ja) 2000-12-28 2007-02-07 花王株式会社 肌色測定装置および肌画像処理装置
JP5596290B2 (ja) * 2006-01-12 2014-09-24 コーニンクレッカ フィリップス エヌ ヴェ ラジオロジーにおける患者皮膚線量の改良された表示
JP5274526B2 (ja) * 2010-09-09 2013-08-28 三菱電機株式会社 皮膚線量表示装置及び皮膚線量表示方法
CN103702716B (zh) * 2011-08-17 2016-03-30 三菱电机株式会社 皮肤剂量评估辅助装置及治疗计划装置
JP6259078B2 (ja) * 2013-06-27 2018-01-10 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 外照射療法における皮膚火傷のリアルタイム定量化

Also Published As

Publication number Publication date
US10315052B2 (en) 2019-06-11
WO2015146164A1 (ja) 2015-10-01
EP3124077A4 (en) 2017-11-22
EP3124077B1 (en) 2019-02-27
EP3124077A1 (en) 2017-02-01
US20170128748A1 (en) 2017-05-11
JPWO2015146164A1 (ja) 2017-04-13

Similar Documents

Publication Publication Date Title
US11452455B2 (en) Skin reflectance and oiliness measurement
JP6448142B2 (ja) 放射線照射による皮膚変化予測装置と検証装置
US10638968B2 (en) Skin gloss evaluation device, skin gloss evaluation method, and skin gloss evaluation program
CN110522466B (zh) 确定患者重量和/或体重指数的方法和装置
JP6040327B2 (ja) 肺の測定
KR101711228B1 (ko) 신체, 특히, 종양의 용적 변화를 수반하는 병리 진행을 정량화하기 위한 방법
Tsang et al. Interinstitutional measurements of left ventricular volumes, speckle-tracking strain, and dyssynchrony using three-dimensional echocardiography
Valentinitsch et al. Regional analysis of age-related local bone loss in the spine of a healthy population using 3D voxel-based modeling
Zhang et al. The effect of nonlinear human visual system components on performance of a channelized Hotelling observer in structured backgrounds
US11039807B2 (en) System and method for dark-field-imaging
EP3082105B1 (en) Image processing apparatus, image processing system, image processing method, and program
Nightingale et al. Frugal 3D scanning using smartphones provides an accessible framework for capturing the external ear
Nölke et al. Invertible neural networks for uncertainty quantification in photoacoustic imaging
US9465983B2 (en) Method and imaging apparatus to automatically display and/or measure bone variations in medical image data
Qaq et al. Sex estimation using lateral cephalograms: a statistical analysis
Marcal et al. Evaluation of the Menzies method potential for automatic dermoscopic image analysis.
EP2980757B1 (en) Quantification and imaging methods of the echo-texture feature
GB2577700A (en) Method and device for determining nature or extent of skin disorder
US20220301687A1 (en) Uncertainty maps for deep learning eletrical properties tomography
KR101587720B1 (ko) 의료 영상의 분산 이미지 노이즈를 활용한 골밀도 진단 장치 및 이의 골밀도 진단 정보 제공 방법
KR102335210B1 (ko) 다채널 근적외선 분광법을 이용한 근피로 측정 방법 및 가시화 장치
US20220068468A1 (en) Medical image processing apparatus, system, and method
Khider et al. Multifractal analysis by the large deviation spectrum to detect osteoporosis

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180129

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180129

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180815

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181015

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20181113

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20181203

R150 Certificate of patent or registration of utility model

Ref document number: 6448142

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees