(57)【特許請求の範囲】
【請求項１】
手術中に複数のカメラが撮像した被写体の画像を合成して一の画像を出力する画像合成装置であって、
第１カメラが出力する第１画像及び第２カメラが出力する第２画像が入力される画像入力部と、
前記第１カメラ及び前記第２カメラの中から最適カメラを決定する最適カメラ決定部と、
前記第１画像及び前記第２画像に基づいて合成画像を生成する画像統合部と、を備え、
前記最適カメラ決定部は、
各画像について、人の頭又は手を含む遮蔽物によって前記被写体が遮られている遮蔽領域を、前記頭又は手に着用されている物の特性を用いて検出し、
各画像における前記遮蔽領域による遮蔽率を算出し、
前記第１カメラ及び前記第２カメラのうち、より低い遮蔽率を有する画像を出力したカメラを最適カメラとして決定し、
前記画像統合部は、
前記最適カメラが出力した画像に遮蔽領域が発生しているか否かを判断し、
前記最適カメラが出力した画像に前記遮蔽領域が発生している場合には、前記最適カメラとは異なる他のカメラが出力した画像によって前記遮蔽領域を補完して合成画像を生成する、ことを特徴とする画像合成装置。
【請求項２】
最適カメラ決定部は、各画像における前記遮蔽領域の中心偏移率をさらに算出し、前記画像統合部は、前記遮蔽率と前記中心偏移率とに基づいて前記最適カメラを決定する、請求項１に記載の画像合成装置。
【請求項３】
前記画像統合部は、前記最適カメラが出力した画像を、前記最適カメラの視点とは異なる他の視点に変換する、請求項１又は２に記載の画像合成装置。
【請求項４】
前記第１カメラが有する第１レンズの第１光軸と、前記第２カメラが有する第２レンズの第２光軸とが重ならないように、前記第１カメラ及び前記第２カメラが配置されており、
前記画像統合部は、
前記第２画像の各画素が、対応する前記第１画像の各画素に重なるように座標変換を実行する、ことを特徴とする請求項１から３のいずれか一項に記載の画像合成装置。
【請求項５】
手術中に複数のカメラが撮像した被写体の画像を合成して一の画像を出力する画像合成方法であって、
第１カメラが出力する第１画像及び第２カメラが出力する第２画像を受け付けるステップと、
各画像について、人の頭又は手を含む遮蔽物によって前記被写体が遮られている遮蔽領域を、前記頭又は手に着用されている物の特性を用いて検出するステップと、
各画像における前記遮蔽領域による遮蔽率を算出するステップと、
前記第１カメラ及び前記第２カメラのうち、より低い遮蔽率を有する画像を出力したカメラを最適カメラとして決定するステップと、
前記最適カメラが出力した画像に遮蔽領域が発生しているか否かを判断するステップと、
前記最適カメラが出力した画像に前記遮蔽領域が発生している場合には、前記最適カメラとは異なる他のカメラが出力した画像によって前記遮蔽領域を補完して合成画像を生成するステップと、
を含むことを特徴とする画像合成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像合成装置及び画像合成方法に関し、具体的には、複数のカメラで撮像した手術画像を合成する画像装置及び画像方法に関する。
【背景技術】
【０００２】
手術映像を記録して残すことは、病院の資産として、学会発表の資料として、また、後進育成の資料としてなど、様々な場面での活用が期待される。さらに、手術中に医師、看護師が情報共有する上で、重要な役割を担っている。
【０００３】
特許文献１には、手術中の術部の画像情報を取得し、取得した情報を画像処理することにより、手術時の動画から患部に関する情報を得ることを可能とした医療画像記録装置が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１１−０３６３７１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記特許文献１には撮像するためのカメラの配置について詳細な記載がないが、いかなる位置に配置したとしても、手術中は医師や看護師が動くため被写体である患部とカメラとの間に遮蔽物（医師の頭等）が介在する場合がある。
【０００６】
このように患部とカメラとの間に遮蔽物が入ると撮像した画像の一部が見えなくなるので、複数のカメラで患部を撮像することが考えられる。
【０００７】
しかしながら、複数のカメラで撮像した複数の画像を同時に表示しても、どのタイミングでどの画像のどの部分が見えなくなるかは予測困難であるため、資料としての画像の使い勝手が悪くなる。
【０００８】
本発明はこのような課題に鑑みてなされたものであり、その目的は、複数のカメラで撮像した複数の画像をより見やすく合成することができる画像合成装置を提供することにある。
【課題を解決するための手段】
【０００９】
上記課題を解決するために、本発明に係る画像合成装置は、複数のカメラが撮像した被写体の画像を合成して一の画像を出力する画像合成装置であって、第１カメラが出力する第１画像及び第２カメラが出力する第２画像が入力される画像入力部と、前記第１画像及び前記第２画像に基づいて合成画像を生成する画像統合部と、を備え、前記画像統合部は、前記第１カメラと前記被写体との間に介在する遮蔽物によって前記第１画像に遮蔽領域が発生しているか否かを判断し、前記第１画像に前記遮蔽領域が発生している場合には、前記第２画像によって前記遮蔽領域を補完して合成画像を生成する、ことを特徴とする。
【００１０】
本発明によれば、第１カメラと被写体との間に介在する遮蔽物によって第１カメラが出力する第１画像に遮蔽領域が発生している場合に、第２カメラが出力する第２画像によって遮蔽領域を補完することで、遮蔽領域のない手術映像を提供することができる。遮蔽領域のない手術映像により、医師や看護師は、手術中、最良の手術映像を共有しながら手術を行うことができる。また、煩雑な確認・編集作業を必要としない、使い勝手の良い資料として活用することができる。
【発明の効果】
【００１１】
本発明によれば、複数のカメラで撮像した複数の画像をより見やすく合成することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態に係る画像合成装置を有する手術映像記録システムの概略図である。
【図２】本実施形態における画像合成の概念図である。
【図３】本実施形態における画像合成の概念図である。
【図４】本発明の実施形態に係る画像合成装置の構成を示すブロック図である。
【図５】本発明の実施形態に係る画像合成処理を示すフローチャートである。
【発明を実施するための形態】
【００１３】
以下、添付図面を参照しながら本発明の実施形態について説明する。説明の理解を容易にするため、各図面において同一の構成要素に対しては可能な限り同一の符号を付して、重複する説明は省略する。
【００１４】
図１は、本発明の実施形態に係る画像合成装置を有する手術映像記録システムの概略図である。図１では、手術映像記録システムを設置した手術室内の状況を模式的に示している。
【００１５】
画像合成装置１０は、被写体の画像を撮像する複数のカメラＣ０、Ｃ１、Ｃ２と通信可能に接続される。複数のカメラＣ０、Ｃ１、Ｃ２は、各カメラが有するレンズの光軸が重ならないように、手術室内のライトＬに取り付けられ、撮像した画像を画像合成装置１０に出力する。術者の視点からの画像を撮像するカメラＣ０を、基準カメラと称する。本実施形態では、カメラＣ０が本発明の「第１カメラ」に相当し、カメラＣ１，Ｃ２が本発明の「第２カメラ」に相当する。
【００１６】
前述したように、被写体とカメラＣ０との間には遮蔽物として医師または看護師を含む人Ｍが介在する場合がある。画像合成装置１０は、カメラＣ０と被写体との間に介在する人ＭによってカメラＣ０が出力する画像に遮蔽領域が発生している場合に、他のカメラＣ１，Ｃ２が出力する画像によって遮蔽領域を補完して合成画像を生成する。
【００１７】
画像合成装置１０が生成した合成画像は、画像表示装置３０に出力される。画像表示装置３０は、ディスプレイ３０１を有している。画像合成装置１０から出力された合成画像は、ディスプレイ３０１に表示される。
【００１８】
図２及び図３は、本実施形態における画像合成の概念図である。図２においては、被写体Ｂ１とカメラＣ０との間に、遮蔽物Ｂ２が介在している。カメラＣ０が撮像した画像Ｍ０においては、被写体Ｂ１の一部が遮蔽物Ｂ２によって遮られていて遮蔽領域となっている。一方、カメラＣ１及びカメラＣ２が撮像した画像Ｍ１及び画像Ｍ２においては、被写体Ｂ１は遮蔽されずに撮像されている。
【００１９】
図２においては、空間上の点Ｐｉは、画像Ｍ０においては点ｐｉ０として映っており、画像Ｍ１においては点ｐｉ１として映っており、画像Ｍ２においては点ｐｉ２として写っている。本実施形態では、点ｐｉ１及び点ｐｉ２を点ｐｉ０の位置に変換することで合成画像を生成している。
【００２０】
図３の（Ａ）に示される画像Ｍ１及び図３の（Ｃ）に示される画像Ｍ２を変換し、図３の（Ｂ）に示される画像Ｍ０を補完することで、図３の（Ｄ）に示される合成画像ＭＸが生成される。合成画像ＭＸにおいては、遮蔽物Ｂ２は画像から消えており、被写体Ｂ１の全体を映すことができる。
【００２１】
図４は、本発明の実施形態に係る画像合成装置１０の構成を示すブロック図である。画像合成装置１０は、カメラ画像入力部１０１と、最適カメラチャンネル判断部１０２と、画像補正・統合部１０３と、合成画像出力部１０４と、を備えている。
【００２２】
カメラ画像入力部１０１は、カメラＣ０、Ｃ１、Ｃ２から出力される画像Ｍ０、画像Ｍ１、画像Ｍ２を受け付ける。
【００２３】
カメラ画像入力部１０１から出力された画像データは、最適カメラチャンネル判断部１０２に入力される。最適カメラチャンネル判断部１０２は、最適カメラがどのカメラなのかを判断する。最適カメラチャンネル判断部１０２は、画像Ｍ０、画像Ｍ１、画像Ｍ２の遮蔽領域を検出する。次に、最適カメラチャンネル判断部１０２は、検出した遮蔽領域に基づく遮蔽率と中心偏移率とに基づいて、各カメラについて画像の適合率を算出する。
【００２４】
ここで、中心偏移率とは、遮蔽領域が画像Ｍの中心部にあるか端部にあるかに関する指標であり、例えば、画面対角距離の１／２を画素数でＲとし、画像中心から遮蔽領域の重心までの距離を画素数でＣとしたとき、中心偏移率＝Ｃ／Ｒ（０〜１までの値をとり、０のときに中心偏移率なし。）と定義することができる。
【００２５】
また、適合率については、例えば、適合率＝１−（遮蔽率×（１−中心偏移率））（０〜１までの値をとり、１のときに適合率が最良。）と定義することができる。
【００２６】
最後に、最適カメラチャンネル判断部１０２は、算出した適合率に基づいて、最適カメラを決定する。遮蔽率に加えて中心偏移率を考慮することで、画像の中心に撮像されると想定される被写体に関する遮蔽に対して、遮蔽率の重みづけをすることができる。
【００２７】
最適カメラチャンネル判断部１０２から出力された画像データは、画像補正・統合部１０３に入力される。画像補正・統合部１０３は、それらの画像データに基づいて合成画像を生成する。具体的には、最適カメラが基準カメラでない場合、画像補正・統合部１０３は、最適カメラの画像を基準カメラの視点に変換する。さらに、画像補正・統合部１０３は、最適カメラの画像に遮蔽領域が発生しているか否かを判断し、遮蔽領域が発生している場合には、最適カメラ以外の他カメラの画像を基準カメラの視点に変換して、変換した他カメラの画像によって最適カメラの遮蔽領域を補完して合成画像ＭＸを生成する。
【００２８】
画像補正・統合部１０３から出力された合成画像ＭＸの画像データは、合成画像出力部１０４に入力される。合成画像出力部１０４は、合成画像ＭＸの画像データを画像表示装置３０に出力する。
【００２９】
（画像合成処理）
続いて、図５を参照して、本発明の実施形態に係る画像合成処理を詳細に説明する。
【００３０】
ステップＳ３０１において、画像合成装置１０のカメラ画像入力部１０１が、カメラＣ０、Ｃ１、Ｃ２から出力される画像Ｍ０、画像Ｍ１、画像Ｍ２を受け付ける。
【００３１】
ステップＳ３０１に続くステップＳ３０２において、最適カメラチャンネル判断部１０２が、画像Ｍ０、画像Ｍ１、画像Ｍ２の遮蔽領域を検出する。遮蔽領域は、頭及び手とすることもできるし、手については手技を確認するための必要な画像であるとして、頭のみを遮蔽領域と定義することもできる。また、遮蔽領域の検出手法としては、手術用キャップ、及び手術用手袋の特性（色、形等）を用いることも考えられるし、他の手法を用いて遮蔽領域を検出することもできる。
【００３２】
ステップＳ３０２に続くステップＳ３０３において、最適カメラチャンネル判断部１０２は、検出した遮蔽領域に基づく遮蔽率と中心偏移率とに基づいて、各カメラについて画像の適合率を算出する。
【００３３】
ステップＳ３０４に続くステップＳ３０４において、最適カメラチャンネル判断部１０２は、算出した適合率に基づいて、最適カメラを決定する。適合率が最大のカメラを最適カメラとして決定する。尚、ステップＳ３０３及びステップＳ３０４を省略し、予め定められた基準カメラを最適カメラとみなして画像合成を実行してもよい。
【００３４】
ステップＳ３０４に続くステップＳ３０５において、画像補正・統合部１０３は、後述する座標変換処理によって、最適カメラ以外のカメラが撮像した画像を、基準カメラの視点に変換する。最適カメラが基準カメラでない場合、最適カメラの画像を基準カメラの視点に変換する。
【００３５】
ステップＳ３０５に続くステップＳ３０６においては、画像補正・統合部１０３が、ステップＳ３０２において検出した遮蔽領域が無くなるように、ステップＳ３０５において変換した画像を合成する。
【００３６】
ステップＳ３０６に続くステップＳ３０７において、合成画像出力部１０４は、画像補正・統合部１０３が生成した合成画像を出力する。
【００３７】
（事前処理）
本発明の実施形態に係る画像合成処理を行う前に、画像補正・統合部１０３は、カメラＣ１の画像Ｍ１から基準カメラであるカメラＣ０の画像Ｍ０へ変換する射影パラメータＨ１、及びカメラＣ２の画像Ｍ２からカメラＣ０の画像Ｍ０へ変換する射影パラメータＨ２を算出するために、各カメラＣのカメラ特性を表すカメラパラメータを決定する。
【００３８】
レンズの焦点距離、レンズ歪み、画像中心（光軸と画像面の交点）の座標等を含むカメラパラメータは、カメラＣで撮像した画像の特徴点と、３次元座標系内の特徴点の位置の対応関係から、計算によって求めることができる。例えば、画像補正・統合部１０３は、キャリブレーションパターンを撮像したカメラＣの画像から格子点の位置を検出して、検出した格子点の位置と３次元座標系の格子点の位置との対応関係から、カメラパラメータを計算することができる。
【００３９】
次に、画像補正・統合部１０３は、計算した各カメラＣのカメラパラメータに基づいて、射影パラメータＨ１、及び射影パラメータＨ２を算出する。最後に、画像補正・統合部１０３は、算出した射影パラメータＨ１、及び射影パラメータＨ２を用いて、カメラＣ１の画像Ｍ１、及びカメラＣ２の画像Ｍ２をカメラＣ０の視点に変換し、変換した画像をカメラＣ０の画像Ｍ０に重ねて、重畳した画像のずれが所定の閾値以内であるか否かを判断する。画像のずれが所定の閾値以内である場合には、画像補正・統合部１０３は、カメラパラメータ及び射影パラメータＨを記憶部（図示せず）に格納して事前処理を終了する。
【００４０】
（射影パラメータＨの算出方法）
図４に３次元空間上の点Ｐi（i：点の番号）と対応するカメラＣｋ（ｋ＝０、１、２）の画像の２次元点ｐｉｋの関係を示す。３次元空間上の点Ｐiと画像上の点ｐｉｋの間の関係について同次座標を用いてベクトル式で記述すると式（ｆ０１）となる。
【００４１】
【数１】
0006529345.tif 000002
【００４２】
ここで、行列ＭｋＩ、行列ＭｋＥはそれぞれカメラＣｋの内部及び外部パラメータを示し、λkは奥行情報を示す不定値である。
【００４３】
式（ｆ０１）を要素に転換すると、式（ｆ０２）となる。
【００４４】
【数２】
0006529345.tif 000003
【００４５】
ここで、ｆｘ、ｆｙは焦点距離、ｃｘ、ｃｙは光軸の画像中心からのずれを示す座標値、Ｒ、ｔは３次元座標系からカメラ座標系への回転行列と平行移動ベクトルとをそれぞれ示す。
【００４６】
（射影パラメータＨの算出方法１：全ての対象点が平面上にある場合（弱透視投影））
このとき、弱透視投影による制限から点ＰiのＺ値はＺi＝0となり、式（ｆ０２）は式（ｆ０３）のように簡略化される。
【００４７】
【数３】
0006529345.tif 000004
【００４８】
両辺をλkで割り、左辺のｘｉｋ／λｋ及びｙｉｋ/λkを簡単のため再度ｘｉｋ及びｙｉｋと表記し、右辺の1/λkは行列ＭｋＩの中に入れ全体を射影変換行列Ｈkとまとめると、式（ｆ０４）となり、最終的に式（ｆ０５）となる。
【００４９】
【数４】
0006529345.tif 000005
【００５０】
【数５】
0006529345.tif 000006
【００５１】
また、カメラＣｍ（ｋ≠ｍ）に関しては、式（ｆ０６）より、カメラＣｍの画像とカメラＣｋの画像の画像間の対応は、式（ｆ０７）となる。
【００５２】
【数６】
0006529345.tif 000007
【００５３】
【数７】
0006529345.tif 000008
【００５４】
行列Ｈmは３行３列のマトリクスであるのでｄｅｔＨm≠0であれば逆行列が存在する。これが、カメラＣｍの画像からカメラＣｋの画像へ変換する射影変換行列Ｈm2kである。
【００５５】
内部パラメータ、外部パラメータとの関係で射影変換行列Ｈm2kを記述すると、式（ｆ０８）となる。
【００５６】
【数８】
0006529345.tif 000009
【００５７】
すなわち、カメラＣｋとカメラＣｍの内部パラメータと外部パラメータとが既知であれば、各カメラ画像間の射影変換行列を算出することができる。
【００５８】
（射影パラメータＨの算出方法２：一般の場合（すなわち制限のない透視投影））
【００５９】
【数９】
0006529345.tif 000010
【００６０】
両辺をλkで割り、左辺のｘｉｋ／λｋ及びｙｉｋ/λkを簡単のため再度ｘｉｋ及びｙｉｋと表記し、右辺の1/λkは[r1,r2,r3,t]の中に入れ全体を射影変換行列Ｈkとまとめると、式（ｆ１０）となり、最終的に式（ｆ１１）となる。
【００６１】
【数１０】
0006529345.tif 000011
【００６２】
【数１１】
0006529345.tif 000012
【００６３】
また、カメラＣｍ（ｋ≠ｍ）に関しては、式（ｆ１２）となる。
【００６４】
【数１２】
0006529345.tif 000013
【００６５】
しかし、この場合、行列Ｈmは４行３列のマトリクスであるので、通常の逆行列の存在が保証されない。したがって、一般化逆行列である式（ｆ１３）を用いる。
【００６６】
【数１３】
0006529345.tif 000014
【００６７】
これにより、カメラＣｍの画像とカメラＣｋの画像の画像間の対応は、式（ｆ１４）となる。
【００６８】
【数１４】
0006529345.tif 000015
【００６９】
これが、カメラＣｍの画像からカメラＣｋの画像へ変換する射影変換行列Ｈm2kである。
【００７０】
内部パラメータ、外部パラメータとの関係で射影変換行列Ｈm2kを記述すると、式（ｆ１５）となる。
【００７１】
【数１５】
0006529345.tif 000016
【００７２】
すなわち、カメラＣｋとカメラＣｍの内部パラメータと外部パラメータならびに不定項λk が既知であれば、各カメラ画像間の射影変換行列を算出することができる。
【００７３】
（座標変換処理）
画像補正・統合部１０３は、前述した合成画像を生成するために、カメラＣ１の画像Ｍ１、及びカメラＣ２の画像Ｍ２をカメラＣ０の視点に変換する。具体的には、オフライン時に記憶部に格納した射影パラメータＨ、及びカメラ画像入力部１０１からの画像を用いて、式（ｆ１６）、式（ｆ１７）に基づいてカメラＣ１（またはカメラＣ２）の画像Ｍ１（または画像Ｍ２）上の点（p,q）をカメラＣ０の視点上の点（p’,q’）に座標変換する。
【００７４】
【数１６】
0006529345.tif 000017
【００７５】
【数１７】
0006529345.tif 000018
【００７６】
ここで、式（ｆ１６）、式（ｆ１７）に記載の記号ｚは、座標系のｚ成分を示すものではなく、不定項λkの逆数をｚと便宜上表記したものである。なお、この不定項λkはそれぞれの画素（p,q）に付随するものであり、空間上でのカメラ中心点と画素（p,q）を結ぶ直線を延長したとき最初に交差する物体までの実距離を示す。
【００７７】
さらに、画像補正・統合部１０３は、点（p’,q’）に座標変換される点（p,q）のＲＧＢ値（r,g,b）、及び画像の重なり情報αを用いて、式（ｆ１８）、式（ｆ１９）、式（ｆ２０）に基づいて、座標変換後の点（p’,q’）におけるＲＧＢ値（r’,g’,b’）を算出する。ここで、最適カメラの遮蔽領域を補完する他カメラが１つである場合、αは１となり、補完する他カメラが２つである場合、αは１／２となる。
【００７８】
【数１８】
0006529345.tif 000019
【００７９】
【数１９】
0006529345.tif 000020
【００８０】
【数２０】
0006529345.tif 000021
【００８１】
本発明によれば、第１カメラと被写体との間に介在する遮蔽物によって第１カメラが出力する第１画像に遮蔽領域が発生している場合に、第２カメラが出力する第２画像によって遮蔽領域を補完することで、遮蔽領域のない手術映像を提供することができる。遮蔽領域のない手術映像により、医師や看護師は、手術中、最良の手術映像を共有しながら手術を行うことができる。また、煩雑な確認・編集作業を必要としない、使い勝手の良い資料として活用することができる。
【符号の説明】
【００８２】
Ｃ０、Ｃ１、Ｃ２：カメラ
Ｌ：ライト
Ｍ：人
１０：画像合成装置
１０１：カメラ画像入力部
１０２：最適カメラチャンネル判断部
１０３：画像補正・統合部
１０４：合成画像出力部
３０：画像表示装置
３０１：ディスプレイ
Ｂ１：被写体
Ｂ２：遮蔽物
Ｍ０、Ｍ１、Ｍ２：画像
ＭＸ：合成画像
【図１】
0006529345.tif 000022
【図２】
0006529345.tif 000023
【図３】
0006529345.tif 000024
【図４】
0006529345.tif 000025
【図５】
0006529345.tif 000026