A Madara index representing the non-uniformity of the ultrasound image, which can help to properly diagnose the contrast findings, is calculated. For this purpose, a grid is set for the region of interest set for the ultrasonic image (steps S1 to S6), and each cell divided by the grid for each frame of the ultrasonic image is set to a luminance level of three gradations. Classify (step S8). Each cell is set as the target cell in order, and the target cell is compared with the luminance level of the target cell and the luminance level of the target cell for comparison. When there are two gradations, a mottle is determined and a madara index representing the number of mottles in one frame is obtained (step S9). The average of the madara index of each frame is calculated as an average madara index (step S12), and this is displayed on the display device and the madara index is displayed in a graph in time series of frames (step S13).造影所見の診断を適切に行なうための一助となり得る、超音波画像の不均一性を表すマダラ指数を演算する。そのために、超音波画像に対して設定された関心領域に対してグリッドを設定し（ステップＳ１～Ｓ６）、超音波画像のフレーム毎にグリッドで区分される各セルを３階調の輝度レベルに分類する（ステップＳ８）。各セルを順に対象セルとし、対象セルの上下左右のセルのうち３階調のいずれかの輝度レベルが設定されているセルを比較対象とし、対象セルの輝度レベルと比較対象のセルの輝度レベルとの差が２階調あるとき、まだらと判断し、１フレームにおけるまだらの数を表すマダラ指数を得る（ステップＳ９）。各フレームのマダラ指数の平均を平均マダラ指数として演算し（ステップＳ１２）、これを表示装置に表示すると共にフレームの時系列にマダラ指数をグラフ表示する（ステップＳ１３）。
