PROBLEM TO BE SOLVED: To provide technique for imaging an object with higher accuracy as compared to before.SOLUTION: An optical ultrasound imaging apparatus 10 is provided for performing imaging processing in which light is made to irradiate an object and the ultrasound generated from the object is detected with multiple sensors 22, and thereby the object is imaged based on sound pressure data comprising multiple detection signals to show the sound pressure spatial distribution. Multiplication results obtained when an initial sound pressure vector is multiplied to the projection matrix of sound speed data based on the predetermined initial sound speed data are assumed as the sound pressure data. Then, the sound speed vector comprising the sound speed of multiple areas having different sound speed is changed so that an objective function value being an objective function concerning the initial sound pressure vector and representing the sparse characteristics of the image data of the object or the sparse characteristics of base data of conversion obtained when a predetermined conversion is made to the image data may become equal to or lower than a predetermined threshold, and the sound pressure vector and the sound speed vector are estimated at the same time, and thereby the image data being imaged based on the sound speed vector is calculated.SELECTED DRAWING: Figure 1COPYRIGHT: (C)2019,JPO&INPIT【課題】従来に比較して高精度で対象物を画像化する技術を提供する。【解決手段】対象物に対して光を照射して対象物から発生する超音波を複数のセンサ22で検出し、音圧空間分布を示す複数の検出信号からなる音圧データに基づいて対象物を画像化する画像化処理を行う光超音波画像化装置10において、所定の初期音速データに基づいて、音速データのプロジェクション行列に初期の音圧ベクトルを乗算したときの乗算結果を音圧データとするときに、初期の音圧ベクトルに関する目的関数であって、対象物の画像データのスパース性、もしくは、画像データに対して所定の変換を行ったときの変換の基底のデータのスパース性を表す目的関数の値が所定のしきい値以下になるように、音速が異なる複数の領域の音速を含む音速ベクトルを変化させて、音圧ベクトルと音速ベクトルを同時に推定することで、音速ベクトルに基づく画像化された画像データを計算する。【選択図】図1