"被検体の診断対象を含む所定領域内に超音波を送信し、前記所定領域内からの反射波を受信し、前記反射波に基づいて超音波画像データを取得する画像データ取得ユニットと、前記超音波画像データの階調補正を実行する補正ユニットと、を具備し、前記補正ユニットは、前記階調補正において、前記画像データの輝度に関するヒストグラムを算出し、前記ヒストグラムを用いて前記診断対象に対応する輝度分布範囲と階調補正関数とを算出する超音波診断装置。前記補正ユニットは、前記ヒストグラムの累積和を算出し、前記輝度分布範囲における前記累積和の傾きを利用して前記階調補正関数の傾きを算出する請求項1記載の超音波診断装置。前記補正ユニットは、前記ヒストグラムの二階微分から得られる特徴点を用いて、前記輝度分布範囲と前記階調補正関数とを算出する請求項1記載の超音波診断装置。前記補正ユニットは、前記ヒストグラムの変曲点を前記特徴点として、前記輝度分布範囲と前記階調補正関数とを算出する請求項3記載の超音波診断装置。前記補正ユニットは、前記ヒストグラム二階微分の最低値、最大値、第2零点、第3零点の少なくともいずれかを前記特徴点として、前記階調補正関数を算出する請求項3記載の超音波診断装置。前記補正ユニットは、入力装置を介して指定される任意の点を前記特徴点として、前記階調補正関数を算出する請求項3記載の超音波診断装置。前記補正ユニットは、前記ヒストグラムのn階微分ただし、nは自然数から得られる特徴点を用いて、前記輝度分布範囲と前記階調補正関数とを算出する請求項1乃至3のうちいずれか一項記載の超音波診断装置。前記補正ユニットは、前記ヒストグラムの変曲点を前記特徴点として、前記輝度分布範囲と前記階調補正関数とを算出する請求項7記載の超音波診断装置。前記補正ユニットは、前記ヒストグラムn階微分の第i番目の極値、第j番目の零点ただし、n,i,jはそれぞれ自然数の少なくともいずれかを前記特徴点として、前記階調補正関数を算出する請求項7記載の超音波診断装置。前記補正ユニットは、入力装置を介して指定される任意の点を前記特徴点として、前記階調補正関数を算出する請求項7記載の超音波診断装置。前記補正ユニットは、ゲイン調整がなされた場合には、ゲイン調整後の前記画像データの輝度に関する前記ヒストグラムを算出する請求項1乃至10のうちいずれか一項記載の超音波診断装置。前記診断対象は心臓であり、前記補正ユニットは、前記ヒストグラムの心筋に対応する領域を前記輝度分布範囲として決定する請求項1乃至11のうちいずれか一項記載の超音波診断装置。前記階調補正の開始指示を入力するための入力ユニットをさらに具備し、前記補正ユニットは、前記入力ユニットを介して入力される前記開始指示に応答して、前記階調補正を実行する請求項1乃至12のうちいずれか一項記載の超音波診断装置。前記補正ユニットは、前記画像データ取得ユニットが超音波画像データを取得する毎に、前記階調補正を実行する請求項1乃至13のうちいずれか一項記載の超音波診断装置。被検体の診断対象を含む所定領域内に超音波を送信し、前記所定領域内からの反射波を受信し、前記反射波に基づいて取得された超音波画像データを記憶する記憶ユニットと、前記超音波画像データの階調補正を実行する補正ユニットと、を具備し、前記補正ユニットは、前記階調補正において、前記画像データの輝度に関するヒストグラムを算出し、前記ヒストグラムを用いて前記診断対象に対応する輝度分布範囲と階調補正関数とを算出する超音波画像処理装置。被検体の診断対象を含む所定領域に関する画像データを取得する画像データ取得ユニットと、前記波画像データの階調補正を実行する補正ユニットと、を具備し、前記補正ユニットは、前記階調補正において、前記画像データの輝度に関するヒストグラムを算出し、前記ヒストグラムを用いて前記診断対象に対応する輝度分布範囲と階調補正関数とを算出する医用画像診断装置。被検体の診断対象を含む所定領域に関する画像データを記憶する記憶ユニットと、前記画像データの階調補正を実行する補正ユニットと、を具備し、前記補正ユニットは、前記階調補正において、前記画像データの輝度に関するヒストグラムを算出し、前記ヒストグラムを用いて前記診断対象に対応する輝度分布範囲と階調補正関数とを算出する医用画像処理装置。コンピュータに、被検体の診断対象を含む所定領域に関する画像データの階調補正を実行させるものであって、前記階調補正において、前記画像データの輝度に関するヒストグラムを算出させるヒストグラム算出機能と、前記ヒストグラムを用いて前記診断対象に対応する輝度分布範囲と階調補正関数とを算出させる関数算出機能と、を実現させる医用画像処理プログラム。ところで、生体内の構造、部位によっては、診断に不要なノイズ等が多くなり、超音波画像のコントラストを好適にするための調整が困難な場合や、ユーザ自身が調整可能な範囲内では適切に調整することができない場合等がある。しかしながら、従来の超音波診断装置は操作性が低く、例えば上記各場合において、画像の階調の設定調整を適切に支援することができない。上記事情を鑑みてなされたもので、従来の超音波診断装置に比して、画像の階調を適切に設定調整することができる超音波診断装置、超音波画像処理装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラムを提供することを目的としている。一実施形態に係る超音波診断装置は、被検体の診断対象を含む所定領域内に超音波を送信し、前記所定領域内からの反射波を受信し、前記反射波に基づいて超音波画像データを取得する画像データ取得ユニットと、前記超音波画像データの階調補正を実行する補正ユニットと、を具備し、前記補正ユニットは、前記階調補正において、前記画像データの輝度に関するヒストグラムを算出し、前記ヒストグラムを用いて前記診断対象に対応する輝度分布範囲と階調補正関数とを算出すること、を特徴とするものである。以下、本実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。第1実施形態図1は、本実施形態に係る超音波診断装置1のブロック構成図を示している。同図に示すように、本超音波診断装置1は、超音波プローブ12、入力装置13、モニター14、超音波送信ユニット21、超音波受信ユニット22、Bモード処理ユニット23、血流検出ユニット24、RAWデータメモリ25、画像処理ユニット26、階調最適化ユニット280を有する表示処理ユニット28、制御プロセッサCPU29、記憶ユニット30、インタフェースユニット31を具備している。以下、個々の構成要素の機能について説明する。超音波プローブ12は、被検体に対して超音波を送信し、当該送信した超音波に基づく被検体からの反射波を受信するデバイス探触子であり、その先端に複数に配列された圧電振動子、整合層バッキング材等を有している。圧電振動子は、超音波プローブ12は、超音波送信ユニット21からの駆動信号に基づきスキャン領域内の所望の方向に超音波を送信し、当該被検体からの反射波を電気信号に変換する。整合層は、当該圧電振動子に設けられ、超音波エネルギーを効率良く伝播させるための中間層である。バッキング材は、当該圧電振動子から後方への超音波の伝播を防止する。当該超音波プローブ12から被検体Pに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ12に受信される。このエコー信号の振幅は、反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分に依存して、周波数偏移を受ける。なお、ボリュームデータを取得する場合には、超音波プローブ12として、例えば二次元アレイプローブ複数の超音波振動子が二次元マトリックス状に配列されたプローブ、又はメカニカル4Dプローブ超音波振動子列をその配列方向と直交する方向に機械的に煽りながら超音波走査を実行可能なプローブを採用するようにしてもよい。しかしながら、当該例に拘泥されず、例えば一次元アレイプローブを超音波プローブ12として採用し、これを手動によって揺動させながら超音波走査をすることでも、ボリュームデータを取得することは可能である。入力装置13は、装置本体11に接続され、オペレータからの各種指示、条件、関心領域ROIの設定指示、種々の画質条件設定指示等を装置本体11にとりこむための各種スイッチ、ボタン、トラックボール、マウス、キーボード等を有している。また、入力装置13は、後述する階調最適化機能において、その開始指示や対象とする超音波画像の選択指示等を有している。モニター14は、画像処理ユニット28からのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。超音波送信ユニット21は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。トリガ発生回路では、所定のレート周波数frHz周期;1/fr秒で、送信超音波を形成するためのトリガパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各トリガパルスに与えられる。パルサ回路は、このトリガパルスに基づくタイミングで、プローブ12に駆動パルスを印加する。なお、超音波送信ユニット21は、制御プロセッサ28の指示に従って所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に送信駆動電圧の変更については、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は複数の電源ユニットを電気的に切り替える機構によって実現される。超音波受信ユニット22は、図示していないアンプ回路、A/D変換器、加算器等を有している。アンプ回路では、プローブ12を介して取り込まれたエコー信号をチャンネル毎に増幅する。A/D変換器では、増幅されたエコー信号に対し受信指向性を決定し、受信ダイナミックフォーカスを行うのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。Bモード処理ユニット23は、受信ユニット22からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるデータを生成する。血流検出ユニット24は、受信ユニット22から受け取ったエコー信号から血流信号を検出し、血流データを生成する。血流信号の検出は、通常CFMColorFlowMappingで行われる。この場合、血流信号を解析し、血流データとして平均速度、分散、パワー等の血流情報を多点について求める。RAWデータメモリ25は、Bモード処理ユニット23から受け取ったBモードデータ、血流検出ユニット24から受け取った血流データを用いて、フレーム毎のBモードRAWデータ、血流RAWデータをそれぞれ生成する。また、RAWデータメモリ25は、必要に応じてRAW-ボクセル変換を実行することにより、RAWデータからボリュームデータを生成する。画像処理ユニット28は、RAWデータメモリ25から受け取るRAWデータに対して、スキャンコンバージョン処理等を実行する。また、画像処理ユニット28は、RAWデータメモリ25から受け取るボリュームデータに対して、ボリュームレンダリング、多断面変換表示MPR:multiplanarreconstruction、最大値投影表示MIP:maximumintensityprojection等の所定の画像処理を行う。なお、ノイズ低減や画像の繋がりを良くすることを目的として、画像処理ユニット28の後に二次元的なフィルタを挿入し、空間的なスムージングを行うようにしてもよい。表示処理ユニット28は、画像処理ユニット28において生成処理された各種画像データに対し、ダイナミックレンジ、輝度ブライトネス、コントラスト、γカーブ補正、RGB変換等の各種を実行する。さらに、表示処理ユニット28が有する階調最適化ユニット280は、後述する階調最適化機能に従う処理を、制御プロセッサ29からの制御に従って実行する。制御プロセッサ29は、情報処理装置計算機としての機能を持ち、本超音波診断装置本体の動作を制御する。制御プロセッサ29は、記憶ユニット30から後述する階調最適化機能を実現するための専用プログラムを読み出して自身が有するメモリ上に展開し、各種処理に関する演算制御等を実行する。記憶ユニット30は、後述する階調最適化機能を実現するための専用プログラムや、診断情報患者ID、医師の所見等、診断プロトコル、送受信条件、スペックル除去機能を実現するためのプログラム、ボディマーク生成プログラム、その他のデータ群が保管されている。また、必要に応じて、RAWデータメモリ中の画像の保管などにも使用される。記憶ユニット30のデータは、インタフェースユニット31を経由して外部周辺装置へ転送することも可能となっている。インタフェースユニット31は、入力装置13、ネットワーク、新たな外部記憶装置図示せずに関するインタフェースである。当該装置によって得られた超音波画像等のデータや解析結果等は、インタフェースユニット31よって、ネットワークを介して他の装置に転送可能である。階調最適化機能次に、本超音波診断装置1が有する、階調最適化機能について説明する。この機能は、取得された個々の超音波画像データを用いて生成されたヒストグラム、当該ヒストグラムの一階微分、二階微分から特徴点を算出し、診断対象に対応する輝度分布範囲を設定する。また、算出された特徴点を用いて制御条件を設定し、設定された制御条件を用いて、画像毎の階調補正関数を算出するものである。図2、図3は、本階調最適化機能に従う処理階調最適化処理において実行される各処理の流れを示したフローチャートである。以下、各ステップにおける処理の内容について説明する。なお、以下の説明においては、説明を具体的にするため、Bモードによって心筋を含む二次元領域を超音波走査し所得する超音波画像に対して階調最適化処理を行う場合を例とする。また、超音波画像処理装置において本階調最適化処理を実現する場合には、例えばステップS2において生成される超音波画像データを予め記憶しておき、ステップS3以降の処理を実行することになる。係る場合の超音波画像処理装置は、図1の破線内の構成を具備するものとなる。[患者情報送受信条件を入力受:ステップS1]入力装置13を介して患者情報の入力、送受信条件被走査領域の大きさ、焦点位置、送信電圧等、被検体の所定領域を超音波走査するための撮像モード、スキャンシーケンス等の選択が実行されるステップS1。入力、選択された各種情報条件等は、自動的に記憶ユニット30に記憶される。[超音波走査超音波画像の生成:ステップS2、S3]超音波プローブ12が被検体表面の所望の位置に当接され、診断部位今の場合、心臓を含む二次元領域を被走査領域として、Bモードによる超音波走査が実行されるステップS2。Bモードによる超音波走査によって取得されたエコー信号は、逐次超音波受信ユニット22を経由してBモード処理ユニット23に送られる。Bモード処理ユニット23は、対数増幅処理、包絡線検波処理等を実行し複数のBモードデータを生成する。RAWデータメモリ25は、Bモード処理ユニット23から受け取った複数のBモードデータを用いてBモードRAWデータを生成する。画像処理ユニット26は、生成されたBモードRAWデータに対してスキャンコンバージョンを実行することで、フレーム毎の超音波画像を生成するステップS3。[階調最適化処理:ステップS4]次に、階調最適化ユニット280は、画像処理ユニット26において生成された超音波画像に対して、図3に示すような階調最適化処理を実行する。すなわち、階調最適化ユニット280は、まず、画像処理ユニット26から受け取った超音波画像の輝度に関するヒストグラムを生成し、得られたヒストグラムに対して平滑化処理を実行するステップS41。その後、階調最適化ユニット280は、得られたヒストグラムを用いて、図4に示すようにヒストグラムの累積和を算出し、平滑化処理を実行するステップS42。次に、階調最適化ユニット280は、ヒストグラムの微分或いは差分を算出した後平滑化処理を実行し、図5に示すような結果を得るステップS43。さらに、階調最適化ユニット280は、ヒストグラムの二階微分或いは二階差分を算出した後平滑化処理を実行し、図6に示すような結果を得るステップS44。次に、階調最適化ユニット280は、算出されたヒストグラムの二階微分を用いて、ヒストグラムの特徴点を算出するステップS45。本実施形態では、ヒストグラムの特徴点の例として、ヒストグラム二階微分の負の最低値ヒストグラムの正のピーク位置最大値或いは極大値に相当、ヒストグラム二階微分の第二零点ヒストグラムの裾領域における変曲点:ヒストグラム微分の負のピーク位置に相当、ヒストグラムの裾領域における変曲点、ヒストグラム二階微分の第三零点例えば、心筋領域からの信号とそれ以外の領域からの信号との分岐点に相当の4点を採用するものとする。次に、階調最適化ユニット280は、算出された特徴点を用いて輝度分布範囲を設定するステップS46。例えば、階調最適化ユニット280は、図7に示すように、ヒストグラム二階微分の負の最低値極大値P1を下限としヒストグラム二階微分の第三零点P4を上限とする輝度分布範囲Lを設定する。次に、階調最適化ユニット280は、算出された特徴点と設定される制御条件と輝度分布範囲Lとを用いて、階調特性関数を算出するステップS47。一例を示せば、階調最適化ユニット280は、以下の制御条件1~5によって階調特性関数を算出する。1ヒストグラム二階微分の負の最低値P1における平滑化累積値の傾きを正規化して、階調特性関数の傾きとする。2入力最小値は出力最小値に対応させる図7下の点B1。3ノイズレベル上限は、変曲点二階微分の第2零点とし、所望の出力値に対応させる図7下の点B2。4対象物領域今の場合心筋領域は、ノイズレベル条件変曲点から第3零点とし、第3零点における出力値は、所望の傾きで求めた値に対応させる図7下の点B3。5入力最大値は出力最大値に対応させる図7下の点B4。次に、階調最適化ユニット280は、ステップS47において得られた階調特性関数を補間例えばスプライン補間して、例えば図7下の階調補正関数C1を算出しステップS48、算出した階調補正関数C1を用いて、超音波画像に対して階調最適化処理を実行するステップS49。[超音波画像表示:ステップS5]次に、モニター14は、階調最適化処理が施された超音波画像を、所定の形態で表示するステップS5。図8は、本階調最適化処理が施されていない超音波画像の一例を示した図であり、図9は、図8の画像に対して本階調最適化処理を施した超音波画像の一例を示した図である。図8と図9とを比較して解るように、本階調最適化処理によって階調調整された超音波画像の方が、コントラストが明瞭であり見易い画像となっている。なお、階調最適化処理の内容は、上記例に拘泥されない。以下、本階調最適化処理の変形例について説明する。変形例1例えば、図10に示すように、第2零点P2における平滑化累積値の傾きを正規化して算出したA2の傾きを採用し、階調補正関数C2を算出することも可能である。また、対象物領域に対応する範囲例えば、輝度分布範囲P2-P3の傾きを平均して算出したものを採用するようにしてもよい。変形例2上記1~5に対し、更なる制御条件を追加して階調特性関数を算出するようにしてもよい。例えば、高輝度領域の傾きを抑制することを目的とした、次の制御条件を加えることができる。6第3零点から入力最大値例えば階調256までの間に存在する入力値を選択し、当該入力値が所望の出力値に対応させる例えば、図11下の点B5。なお、この様な制御条件6の追加は、階調特性関数の算出において制御条件1~5に対し並列的に追加してもよい。また、制御条件1~5を用いて算出された階調特性関数の結果を踏まえ必要に応じて制御条件6を追加し、改めて階調特性関数を算出するようにしてもよい。また、追加した制御条件に対応する位置の階調特性が例えば255を下回る場合には、特性の乱れを抑制するために、階調特性関数が単調増加となるように補正することが好ましい。変形例3必要に応じて階調伸張処理を行うようにしてもよい。階調伸張処理を行うか否かの判定例としては、例えば、入力画像のヒストグラム累積和の所定量以上になる位置が、平滑化の端部処理領域より前において到達する場合には、画像取得時のダイナミックレンジDRや階調特性が適当でなかったとして、例えばステップS47の算出において、図12から図13に示すような階調特性の伸張を行うようにする。変形例4輝度分布範囲は任意に変更することができる。例えば、上述した図7の輝度分布範囲L下限:ヒストグラム二階微分の負の最低値。上限:第三零点を、任意のタイミングで図14に示す輝度分布範囲L1下限:第2零点。上限:第3零点に変更することができる。これにより、対象物をより好適に抽出できる場合がある。変形例5本階調最適化処理は、動画像を取得し表示する場合においても適用可能である。図15は、本階調最適化処理を動画像表示に適用する場合のフローチャートの一例である。図3に示したフローチャートと比較した場合、ステップS40のヒストグラム解析画像の選択処理が追加される点が異なる。ヒストグラム解析画像の選択処理ステップS40においては、取得された複数の超音波画像例えば、一心拍や複数心拍に亘る心臓のついての超音波画像のうち、いずれかの画像を選択してステップS41以降の階調最適化処理が実行されることになる。ヒストグラム解析画像は、例えば、予め設定された時相やタイミングに対応する画像、所定期間の最初、中央或いは最後の時相に対応する画像、画像の輝度平均が最大である画像、画像の輝度平均が複数画像のうちで平均的である画像等を選択することができる。どの様な画像をヒストグラム解析画像とするかは自由であるが、診断対象を心臓とする場合、ヒストグラム累積和は収縮拡張に合わせて周期的に変動する等の理由から、真っ暗な画像、対象物が適切に映像化されていない画像等は、選択しないようにすべきである。変形例6本階調最適化処理は、例えば超音波画像データを取得する毎に実行することができる。また、全ての画像に対して常に実行する必要がない場合には、例えばユーザから入力装置13を介して入力される指示に応答して、任意のタイミングで実行するようにしてもよい。応用例以上述べた実施形態では、超音波画像の輝度値に関するヒストグラムの第二次導関数の第1番目の極値、第2、第3番目の零点等を特徴点として算出する例を示した。しかしながら、本階調最適化機能は、当該例に拘泥されず、例えば第n次導関数の、第i番目の極値、第j番目の零点ただし、n,i,jはそれぞれ任意の自然数を特徴点として算出し、得られた特徴点を用いて輝度分布範囲、階調補正曲線を算出するようにしてよい。以上述べた構成によれば、以下の効果を得ることができる。本超音波診断装置によれば、取得された超音波画像データを用いて生成されたヒストグラム、当該ヒストグラムの一階微分、二階微分等から特徴点を算出し、診断対象に対応する輝度分布範囲を設定する。また、算出された特徴点を用いて制御条件を設定し、設定された制御条件を用いて、階調補正関数を算出する。従って、個々の画像の階調特性に応じた適切な輝度分布範囲、階調補正関数を定義することができる。その結果、診断対象や個体差、撮像状況に応じて、迅速且つ簡単に階調補正を最適化することができる。第2実施形態次に、第2実施形態に係る超音波診断装置について説明する。本階調最適化処理の対象となる超音波画像のゲイン調整を行った場合、ゲイン調整後の輝度値に関するヒストグラムの形状は変化し、その結果ヒストグラムの累積和、第n次階導関数ただし、nは自然数の形状も変化することになる。本実施形態では、このようにゲイン調整によってヒストグラムの形状が変化する場合の階調最適化処理について説明する。所定の超音波画像が取得され、当該画像のゲイン調整前の輝度値に関するヒストグラムが図16の様に得られた場合を想定する。係る場合、ヒストグラムの第一次導関数、第二次導関数、第三次導関数は、それぞれ図17、図18、図19に示す様になる。現状の超音波画像図16のヒストグラムに対応する超音波画像に対してゲインを上げる処理を施した場合、ゲインを上げた後の輝度値に関するヒストグラム、ヒストグラムの第一次導関数、第二次導関数、第三次導関数は、それぞれ図20、図21、図22、図23に示す様になる。ここで、例えば図16と図20とを比較すると、ゲインを上げた結果ヒストグラムのピーク位置が移動右方向にシフトしていることがわかる。また、図17と図21、図18と図22、図19と図23とをそれぞれ比較すると、ゲインを上げた結果、各階の導関数において各極値と各零点の位置も移動右方向にシフトしていることがわかる。特に、図19と図23とを見ると、図19において最初に出現する極値の極性は負であるのに対し、図23において最初に出現する極値の極性は正となっている。従って、例えばゲイン調整前における図19に示す第三次導関数の第1、第2極値を特徴点として採用するという条件は、ゲイン調整後における図23に示すゲイン調整前の第三次導関数の第2、第3極値を特徴点として採用するという条件に対応することになる。また、現状の超音波画像に対してゲインを下げる処理を施した場合、ゲインを下げた後の輝度値に関するヒストグラム、ヒストグラムの第一次導関数、第二次導関数、第三次導関数は、それぞれ図24、図25、図26、図27に示す様になる。ここで、例えば図16と図24とを比較すると、ゲインを下げた結果ヒストグラムのピーク位置が左方向にシフトしていることがわかる。また、図17と図25、図18と図26、図19と図27とをそれぞれ比較すると、ゲインを下げた結果、各次各階の導関数において各極値と各零点の位置も左方向にシフトしていることがわかる。特に、ゲイン上昇前後において第一次導関数、第二次導関数、第三次導関数のそれぞれを比較すると、それぞれ同じ順番の極値の極性が逆転していることがわかる例えば、第一次導関数において、ゲイン上昇前の最初の極値の極性は正であるのに対し、ゲイン上昇後の最初の極値の極性は負となっている。従って、例えばゲイン調整前における図18に示す第二次導関数の第2、第3零点を特徴点として採用するという条件は、ゲイン調整後における図26に示すゲイン調整前の第二次導関数の第1、第2零点を特徴点として採用するという条件に対応することになる。すなわち、ヒストグラムのピーク位置、及びヒストグラムの第n次階導関数の各極値の極性、各零点の位置は、ゲイン調整によって変わる。従って、階調最適化処理における特徴点の算出に用いる第i番目の極値、第j番目の零点i,jは自然数は、ゲイン補正後のヒストグラム、第n次導関数の位置、形状に応じて個別に適切に選択されることが好ましい。なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。1本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスクフロッピー登録商標ディスク、ハードディスクなど、光ディスクCD-ROM、DVDなど、半導体メモリなどの記録媒体に格納して頒布することも可能である。2上記実施形態においては、Bモードによって心筋を含む二次元領域を超音波走査し所得する超音波画像に対して階調最適化処理を行う場合を例とした。しかしながら、本階調最適化機能は、当該例に拘泥されない。例えば、心臓以外の部位を診断対象とする場合、Bモードではなく血流検出モードCFMモード等によって撮像する場合、三次元領域を超音波走査する場合等にも適用可能である。特に、三次元領域を超音波走査する場合においては、ボリュームデータを構成する各二次元画像、MPR画像等のそれぞれに対して、上記階調最適化処理を適用することで、好適な超音波画像を取得することができる。3上記各実施形態においては、超音波診断装置によって取得された超音波画像データに対して階調最適化処理を実行する場合を例に説明した。しかしながら、本階調最適化処理に対象は超音波画像データに限定されず、X線診断装置、X線コンピュータ断層撮像装置、磁気共鳴イメージング装置等の医用画像診断装置によって取得された各医用画像に対しても適用可能である。また、各種医用画像診断装置によって取得された画像を、医用画像処理装置を用いて、階調最適化処理を事後的に行うようにしてもよい。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。"