The scan signal obtained by scanning the target area with the ultrasonic signal is received, and the time required for the ultrasonic signal to pass through the target area varies depending on the temperature of the target area based on the scan signal. An echo shift calculating step for calculating an echo shift that is a change amount, and a moving distance when the ultrasonic signal passes through the target region apparently changes depending on the temperature of the target region based on the calculated echo shift. A strain calculation step for calculating a strain that is a change rate, a strain rate calculation step for calculating a strain rate that is a temporal change rate of the strain based on the calculated strain, a predetermined strain and a strain rate, Corresponds to the calculated strain and strain rate using the relationship with temperature Comprising a temperature estimation step of estimating the temperature of the elephant region.超音波信号によって対象領域をスキャンすることにより得られるスキャン信号を受信し、スキャン信号に基づいて、超音波信号が前記対象領域を通過するのに要する時間が対象領域の温度に依存して変化する変化量であるエコーシフトを算出するエコーシフト算出ステップと、算出されたエコーシフトに基づいて、超音波信号が対象領域を通過する際の移動距離が対象領域の温度に依存して見かけ上変化する変化割合であるストレインを算出するストレイン算出ステップと、算出されたストレインに基づいて、ストレインの時間的な変化割合であるストレインレートを算出するストレインレート算出ステップと、予め定められたストレインとストレインレートと温度との関係を用いて、算出されたストレイン及びストレインレートに対応する対象領域の温度を推定する温度推定ステップと、を含む。
