"ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成する画像推定部と、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出し、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算部を有する画像処理装置。前記画像推定部は、前記ジャギー発生画像から生成した輝度画像を入力してジャギー発生前推定輝度画像を生成し、前記加重加算部は、前記ジャギー発生前推定輝度画像に基づいて前記重みを算出し、前記ジャギー発生画像から生成した輝度画像を適用して、前記重みを適用した画素値の加重加算処理を実行する請求項1に記載の画像処理装置。前記画像推定部は、事前学習データを適用したベイズ推定処理により、ジャギー発生前推定画像を生成する請求項1に記載の画像処理装置。前記画像推定部は、前記ジャギー発生画像の平滑化処理によってジャギー発生前推定画像を生成する請求項1に記載の画像処理装置。前記画像推定部は、エッジを表現する特徴量の統計データを事前学習データとして適用したベイズ推定処理により、ジャギー発生前推定画像を生成する請求項1に記載の画像処理装置。前記加重加算部は、NL-MeansNonLocal-Means法を適用した加重加算処理を実行する請求項1に記載の画像処理装置。前記加重加算部は、JointBirateralFilter法を適用した加重加算処理を実行する請求項1に記載の画像処理装置。前記画像処理装置は、さらに、前記加重加算部の算出した補正画素値を持つジャギー低減輝度画像を参照画像として、JointBilateralFilterJBFを適用した相関処理を実行して補正色差情報を生成するフィルタ処理部を有する請求項1に記載の画像処理装置。前記画像処理装置は、さらに、前記加重加算部の算出した補正画素値を持つジャギー低減輝度画像と、補正前のジャギー発生画像に基づいて生成されるジャギー発生輝度画像とを、画像のエッジ領域に応じた重みに従ってブレンドするブレンド処理部を有する請求項1に記載の画像処理装置。画像処理装置において画像処理を実行する画像処理方法であり、画像推定部が、ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成する画像推定ステップと、加重加算部が、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出し、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算ステップを実行する画像処理方法。画像処理装置において画像処理を実行させるプログラムであり、画像推定部に、ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成させる画像推定ステップと、加重加算部に、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出させ、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算ステップを実行させるプログラム。本開示は、例えば上記問題点に鑑みてなされたものであり、ジャギーの削減を簡易な構成で実現する画像処理装置、および画像処理方法、並びにプログラムを提供することを目的とする。本開示の第1の側面は、ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成する画像推定部と、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出し、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算部を有する画像処理装置にある。さらに、本開示の画像処理装置の一実施態様において、前記画像推定部は、前記ジャギー発生画像から生成した輝度画像を入力してジャギー発生前推定輝度画像を生成し、前記加重加算部は、前記ジャギー発生前推定輝度画像に基づいて前記重みを算出し、前記ジャギー発生画像から生成した輝度画像を適用して、前記重みを適用した画素値の加重加算処理を実行する。さらに、本開示の画像処理装置の一実施態様において、前記画像推定部は、事前学習データを適用したベイズ推定処理により、ジャギー発生前推定画像を生成する。さらに、本開示の画像処理装置の一実施態様において、前記画像推定部は、前記ジャギー発生画像の平滑化処理によってジャギー発生前推定画像を生成する。さらに、本開示の画像処理装置の一実施態様において、前記画像推定部は、エッジを表現する特徴量の統計データを事前学習データとして適用したベイズ推定処理により、ジャギー発生前推定画像を生成する。さらに、本開示の画像処理装置の一実施態様において、前記加重加算部は、NL-MeansNonLocal-Means法を適用した加重加算処理を実行する。さらに、本開示の画像処理装置の一実施態様において、前記加重加算部は、JointBirateralFilter法を適用した加重加算処理を実行する。さらに、本開示の画像処理装置の一実施態様において、前記画像処理装置は、さらに、前記加重加算部の算出した補正画素値を持つジャギー低減輝度画像を参照画像として、JointBilateralFilterJBFを適用した相関処理を実行して補正色差情報を生成するフィルタ処理部を有する。さらに、本開示の画像処理装置の一実施態様において、前記画像処理装置は、さらに、前記加重加算部の算出した補正画素値を持つジャギー低減輝度画像と、補正前のジャギー発生画像に基づいて生成されるジャギー発生輝度画像とを、画像のエッジ領域に応じた重みに従ってブレンドするブレンド処理部を有する。さらに、本開示の第2の側面は、画像処理装置において画像処理を実行する画像処理方法であり、画像推定部が、ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成する画像推定ステップと、加重加算部が、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出し、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算ステップを実行する画像処理方法にある。さらに、本開示の第3の側面は、画像処理装置において画像処理を実行させるプログラムであり、画像推定部に、ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成させる画像推定ステップと、加重加算部に、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出させ、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算ステップを実行させるプログラムにある。なお、本開示のプログラムは、例えば、様々なプログラムコードを実行可能な情報処理装置やコンピュータシステムに対して、コンピュータ可読な形式で提供する記憶媒体、通信媒体によって提供可能なプログラムである。このようなプログラムをコンピュータ可読な形式で提供することにより、情報処理装置やコンピュータシステム上でプログラムに応じた処理が実現される。本開示のさらに他の目的、特徴や利点は、後述する本開示の実施例や添付する図面に基づくより詳細な説明によって明らかになるであろう。なお、本明細書においてシステムとは、複数の装置の論理的集合構成であり、各構成の装置が同一筐体内にあるものには限らない。本開示の一実施例の構成によれば、画像に発生するジャギーを低減させた高品質な画像を生成することが可能となる。具体的には、画像推定部において、例えばベイズ推定によりジャギー発生画像からジャギー発生前画像を生成し、生成したジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出し、処理領域と類似領域に対応するジャギー発生画像の対応領域を検出し、検出した対応領域の画素値に対して、重みを適用した加重加算処理を実行して、ジャギー発生画像の処理領域の補正画素値を算出する。これらの処理により、ジャギーが低減された画像が生成される。以下、図面を参照しながら本開示の画像処理装置、および画像処理方法、並びにプログラムの詳細について説明する。なお、説明は以下の項目に従って行う。1.画像におけるジャギー発生について2.画像処理装置の構成と処理について2-1.処理1RGB画像から輝度画像とその色差情報への色変換処理2-2.処理2ジャギーの発生している輝度画像からベイズ推定によってジャギー発生前の輝度画像を推定する処理2-3.処理3推定された輝度画像から類似度を求め、NL-Meansによってジャギー低減を行う処理2-4.処理4処理3によるジャギー低減画像のエッジ領域と、原画像のエッジ領域以外の部分をブレンドする処理2-5.処理5色差画像に対してジャギー低減済み輝度画像を参照するJBFによる相関処理に基づく補正色差情報の生成2-6.処理6輝度色差合成処理3.画像処理装置の構成例について4.本開示の構成のまとめ[1.画像におけるジャギー発生について]先に説明したように、カメラで撮影した画像を印刷出力や表示出力した場合、出力画像に含まれるオブジェクトの画像輪郭が本来の輪郭と異なり、階段状のギザギザを持つ場合がある。このギザギザは、一般にジャギーと呼ばれる。このジャギーは、例えば画素数が多い高解像度の撮像素子を利用した撮影画像に対して、画素の間引きや合成処理を行って、画素数を削減した出力画像を生成する場合にも発生する。昨今の撮像装置は、高解像度の画像撮影のために、数100万~数1000万画素の極めて多数の画素を持つ撮像素子を備えたものが多くなってきている。しかし、撮像素子の有する画素数に対応する高解像度画像の出力可能な表示装置が利用されることは少なく、また、撮像素子から出力される高解像度画像をメモリにそのまま記録すると記録に必要となるメモリ容量が増大してしまい、記録可能な画像枚数が減少してしまうといった事態が発生する。このような状況を考慮し、高画素数の撮像素子を備えた撮像装置では、撮像素子からの出力画素信号をそのままメモリに記録せず、出力画素数を間引く処理や、複数画素の加算演算などを含む合成処理によって、総画素数を削減してメモリに記録する処理が多く行われている。しかし、このような画素数変換処理によって、出力画像にジャギーが発生するという問題が顕在化してきている。この画素数ジャギーの削減に基づくジャギー発生について、図を参照して説明する。図1には、a撮像素子の画素配列ベイヤーBayer配列、b出力画素の画素重心、これらの図を示している。図1aに示す画素配列は、撮像素子の高画素数の画素配列であり、撮影画像は、高画素数の画素情報を有している。この画素数を削減してメモリ格納画像を生成する。例えば、図1aの撮像素子の同じ色の画素4画素に基づいて、出力画像の1画素の画素値を設定し出力する。すなわち、4画素を1画素に集約して出力して総画素数を削減する。図1b出力画素の画素重心は、画素数削減後の出力画素各々についての元の撮像素子における画素重心を示した図である。例えば、図1bのGb画素31は、図1aの3×3画素ブロック21の4隅のGb画素の画素値を均等に用いて決定される画素値である。3×3画素ブロック21の中心位置に画素重心が設定される。この重心位置を示しているのが、図1bのGb画素31である図1に示す例は、aに示す8×8画素の64画素を、bに示す4×4=16画素として、出力画像の画素数を撮像素子の画素数の1/4に削減する処理を行っている例である。この処理のために、例えば図1aに示す3×3画素のブロック21の4隅のGb画素の加算平均処理を実行して、出力画像における1つのGb画素の画素値を算出している。すなわち、図1aに示す3×3画素のブロック21に含まれる4つのGb画素に基づいて、図1bに示すGb画素31の画素値を算出している。この場合、出力画像におけるGb画素31の重心は、水平右方向をx、垂直下方向をyとした座標軸において、x,y=2,2の位置、すなわち、図1bに示すGb画素31の位置となる。また、図1aに示す3×3画素のブロック22に含まれる4つのB画素に基づいて、図1bに示すB画素32の画素値が算出される。この場合、出力画像におけるB画素32の重心は、x,y=3,2の位置、すなわち、図1bに示すB画素32の位置となる。同様に、図1aに示す3×3画素のブロック23に含まれる4つのGb画素に基づいて、図1bに示すGB画素33の画素値が算出される。この場合、出力画像におけるGb画素33の重心は、x,y=6,2の位置、すなわち、図1bに示すB画素33の位置となる。また、図1aに示す3×3画素のブロック24に含まれる4つのB画素に基づいて、図1bに示すB画素34の画素値が算出される。この場合、出力画像におけるB画素34の重心は、x,y=7,2の位置、すなわち、図1bに示すB画素34の位置となる。図1bに示す計16個の画素は、出力画像とする場合、4×4画素の画像として出力される。すなわち、図2cに示すような4×4画素の画像70として出力される。図2には、b出力画素の画素重心図1bと同じc出力画像の画素位置これらを示している。図2cでは左上の2×2画素、すなわちGb画素31、B画素32を含む2×2画素のブロックを固定して考えている。この設定とした場合、その他の3つの2×2の画素プロックは、いずれも図2cに示す矢印α、β、γに従って移動させて図2cに示す4×4画素の画像70の構成画素として出力されることになる。このシフト処理によって、以下のような問題が発生する。例えば、画素重心の位置がx,y=6,2のGb画素33は、出力画像では、画素位置がx,y=3,2のGb画素53として設定される。また、画素重心の位置がx,y=7,2のB画素34は、出力画像では、画素位置がx,y=3,3のB画素54として設定される。ここで、縮尺率を算出する。画素位置x,y=2,2のGB画素31を固定位置とした基準画素と仮定して考える。図2bに示す画素位置x,y=6,2のGb画素33は、基準画素であるGb画素31から4画素離れている。これが、出力画像では、画素位置がx,y=3,2のGb画素53として設定されるので、その基準画素:Gb画素31からの距離は2画素となる。すなわち、縮尺率は、2画素/4画素=1/2となる。一方、図2bに示す画素位置x,y=7,2のB画素34は、基準画素であるGb画素31から5画素離れている。これが、出力画像では、画素位置がx,y=4,2のB画素54として設定されるので、その基準画素:Gb画素31からの距離は3画素となる。すなわち、縮尺率は、3画素/5画素=3/5となる。このように、画素間の縮尺率がばらばらとなり、撮像素子の画素配列に対応する撮影画像の各画素の相対位置と異なる相対位置を有する出力画像が生成されてしまうことになる。すなわち、撮像素子において撮影されたオリジナル画像の各画素の間隔が不均一に縮小されて出力画像が生成されてしまう。このような画素間隔の不均等は画質劣化を引き起こす。具体的には、例えば図3に示すようなジャギーの拡大などの劣化を発生させる。図3Aのオリジナル画像は、撮像素子の撮影画像に相当する画素数の多い高解像度の画像であり、この画像はジャギーが小さい。この図3Aのオリジナル画像に基づいて、図1~図2を参照して説明した相対的な画素位置をばらばらにしてしまう画素数削減処理を行うと、図3Bに示すようなジャギーの拡大した画像が生成されてしまう。なお、ジャギーは、折り返し雑音の一種である。加算後の画素重心の間隔の不均等により、ジャギー劣化は大きくなる。なお、この他にも、RGBの各配色の位置関係の乱れなどにより、出力画像の色もオリジナル画像と差が発生するなど、様々な画質劣化を引き起こすことになる。本開示に係る画像処理装置では、例えばこのような画素数削減に伴って発生するジャギーの発生を抑制する処理を実現する。なお、ジャギーの発生要因は、上述の画素数削減処理に限らず、様々な要因がある。本開示の画像処理装置は、様々な要因に基づいて発生したジャギーを低減させる処理を行う。具体的には、例えば画像の輝度成分を抽出し、輝度成分に基づいてジャギー発生前の形状の推定を行い、この推定した形状情報を利用して、画像内から類似する画像領域を検出し、検出した類似領域の画素同士を加重平均することでジャギーの低減を図る。この処理は、ジャギー発生要因である例えば画素加算処理や画素間引き処理の態様に依存せず、任意のジャギー発生パターンに対応可能である。また、特定の方向に依存した処理ではないため、任意角度のエッジに発生するジャギーに対して有効な処理となる。[2.画像処理装置の構成と処理について]本開示に係る画像処理装置の構成と処理について、図4以下を参照して説明する。図4は、本開示に係る画像処理装置100の主要部の構成例を示す図である。画像処理装置100は、例えばRGB画像等の入力画像80を入力して、ジャギーを低減した画像である出力画像90を生成して出力する。画像処理装置100は、図4に示すように、輝度色差成分分離部101、画像推定部ベイズ推定部102、加重加算部103、エッジ抽出部104、ブレンド処理部105、フィルタ処理部106、輝度色差成分合成部107を有する。画像処理装置100は、入力画像80に含まれるジャギー、例えば前述した画素加算間引き処理によって既に発生してしまったジャギーを低減する処理を実行する。具体的には、入力画像80の輝度成分に基づいて、一度ジャギー発生前の形状を推定し、その形状情報を利用して画像内の類似領域を探索し、類似領域の画素同士を加重平均することでジャギーの低減を図る。この処理は画素加算間引き処理に依存せず、任意のジャギー発生パターンに対応することが可能である。また、特定の方向に依存した処理を行っていないため、任意の角度のエッジに発生するジャギーに対して有効である。なお、図4に示す画像処理装置100は、画像処理装置の主要部の構成を示しており、画像処理装置100は、この他、例えば全体の処理制御を実行する制御部、例えば所定の処理プログラムの実行機能持つCPUを備えた制御部や、プログラム、画像、各種のパラメータを記憶するためのメモリ等を備えている。画像処理装置100は、具体的には、例えば、PCや、撮像装置カメラ等として実現可能であり、これらの装置においては、図4に示す構成の他、それぞれの装置に応じた構成が具備されたものとなる。画像処理装置100の実行する処理は大きく分けて以下の6つの処理に分類される。処理1RGB画像から輝度画像とその色差へと色変換を行う。またこの逆変換を行う。[輝度色差分離部101、輝度色差合成部107の処理]処理2ジャギーの発生している輝度画像からベイズ推定によってジャギー発生前の輝度画像を推定する。[画像推定部ベイズ推定部102の処理]処理3推定された輝度画像から類似度を求め、NL-Meansによってジャギー低減を行う。[主に加重加算部103の処理]処理4処理3によるジャギー低減画像のエッジ領域と、原画像のエッジ領域以外の部分をブレンドする。[エッジ抽出部104と、ブレンド処理部105の処理]処理5色差画像に対してジャギー低減済み輝度画像を参照するJBFによっての相関処理に基づく補正色差情報の生成。[フィルタ処理部106の処理]処理6輝度色差合成処理。[輝度色差合成部107の処理]以下、この6つの処理の詳細について順に説明する。2-1.処理1RGB画像から輝度画像とその色差情報への色変換処理まず処理1RGB画像から輝度画像とその色差情報への色変換処理について説明する。この処理は、図4に示す画像処理装置100の輝度色差分離部101において実行する処理である。本手法ではジャギー低減処理を輝度成分に対してのみ行う。これは輝度成分がジャギーの影響を特に強く受けるためである。図4に示す輝度色差成分分離部101は、ジャギーが含まれる入力画像例えばRGB画像10を入力し、例えば以下の式1に従って、輝度色差変換を行う。なお、上記式1で示す輝度色差変換処理は、一例であり、この他の様々な既存の処理を適用してもよい。例えば、一般的なYCbCr,Lab方式など、輝度と色差の成分に分離できる方式であれば任意の色変換方式を利用することができる。輝度色差成分分離部101は、入力画像80の各構成画素のRGB値を入力し、上記式1に従って、各画素単位の、輝度:L=R2GB/4色差:Cr=L-R色差:Cb=L-Bこれらを算出する。図4に示すように、輝度色差成分分離部101の算出した輝度[L=R2GB/4]は、画像推定部ベイズ推定部102に入力される。また、色差[Cr=L-R,Cb=L-B]は、フィルタ処理部106に入力される。なお、最終処理として輝度色差合成部107において実行される輝度色差成分からRGB成分への復元処理は、上記の式1をR,G,Bについて解くことで求めることができる。2-2.処理2ジャギーの発生している輝度画像からベイズ推定によってジャギー発生前の輝度画像を推定する処理次にベイズ推定による輝度画像推定方法について説明する。この処理は、図1に示す画像推定部ベイズ推定部102において実行する処理である。画像推定部ベイズ推定部102は、輝度色差成分分離部101の算出した輝度[L=R2GB/4]を入力して、この輝度情報からなる輝度画像に基づいてベイズ推定によってジャギー発生前の輝度画像を推定する。今回の処理対象としているジャギーは、例えば画素加算間引き処理によって発生するジャギーであるため、事前のシミュレーションでジャギー発生パターン、つまり、ある角度、コントラストのエッジ近辺でのジャギーがどのような形状で、どの程度の振幅強度で分布するかを事前に知ることが可能である。画像推定部ベイズ推定部102ではこの条件を利用して、事前にあらゆるパターンのエッジに対して擬似的にジャギーを発生させ、ジャギー発生前後でエッジがどのように変化するかを統計的に観測しておき、その統計量を利用して確率的にジャギー発生前のエッジ形状の推定を行う。すなわち、図4に示す事前学習データ82を予め生成して、この事前学習データを利用して処理対象画像のジャギー発生前のエッジ形状の推定を行う。本手法ではジャギー発生前後の輝度画像の画素値を観測し、ジャギー発生後の画素値分布から尤もらしい画素値を推定する方法を採用している。具体的な手順について説明する。まず、事前学習データ82の生成処理について説明する。図5Aに示すようなエッジを含む学習用画像を用意する。図5Aジャギー発生前の学習用画像である。この画像Aに対してジャギーを発生させて、Bジャギー発生後の学習用画像を生成する。なお、このジャギー発生処理は、実際の処理対象とする画像に対する画素加算間引き処理と同様の処理を実行して行うことが好ましい。次に、Aジャギー発生前の画像内の注目位置i,jの画素値xと、Bジャギー発生後の画像内で位置i,jのn個の近傍画素の画素値y1,y2,...,ynをそれぞれ観測する。図5に示す例では、n=9とした例を示している。学習用画像に含まれる画素について、同様の処理を繰り返していくと、ある画素値xに対してn個のジャギー発生後の画像における近傍画素のヒストグラムを観測することができる。ヒストグラムというのは画素値の出現頻度に他ならないため、観測したヒストグラムは確率分布と考えられる。よって、この事前観測により以下の確率が求められる。Px:ジャギー発生前の画像で画素値xが出現する確率=ジャギー発生前の画像のヒストグラム、Py1,y2,...,yn|x:ジャギー発生前の画像において、画素値xが出現した時に、ジャギー発生後の画像の対応画素位置のn個の近傍画素の画素値が順にy1,y2,...,ynとなる確率、これらの確率が求められる。上記の、Py1,y2,...,yn|xは、あるエッジ領域において発生するジャギーパターンの発生確率に相当する。この算出確率に基づいて、さらに、ベイズの定理より以下の式2に示す確率を求めることができる。式2ただし、Px:ジャギー発生前の画像で画素値xが出現する確率=ジャギー発生前の画像のヒストグラム、Py1,y2,...,yn:ジャギー発生後の画像で画素値y1,y2,...,ynのn個の画素ブロックが出現する確率、Py1,y2,...,yn|x:ジャギー発生前の画像において、画素値xが出現した時に、ジャギー発生後の画像の対応画素位置のn個の近傍画素の画素値が順にy1,y2,...,ynとなる確率、Px|y1,y2,...,yn:ジャギー発生後の画像におけるn個の画素ブロックの画素値が順にy1,y2,...,ynとなる場合に、ブロック中心の画素位置に対応するジャギー発生前の画像の画素位置の画素値がxとなる確率、である。つまり、上記式2に従って、ジャギー発生後画像から選択されるn個の画素ブロックで、画素値:y1,y2,...,ynを観測した時に、この画素ブロックの中央の画素に対応するジャギー発生前画像の画素に画素値xが発生する確率を算出できる。従って、ジャギー発生後の画像における画素ブロックの画素値を観測した時に確率Px|y1,y2,...,ynが最大となるxを求めれば、ジャギー前の画素値を推定することが可能となる。具体的には、以下に示す式3に従って、ジャギー前の画素値xを推定することが可能となる。なお、ここで画素値は、輝度画像の画素値であり輝度値に相当する。式3ただし、x:推定されたジャギー発生前の画素値輝度値ここで、ジャギー発生後の近傍画素値、y1,y2,...,ynを観測する確率は1つの観測領域で一意に決定されるため、上記の式3は以下の式4に示すように簡略化できる。式4画像推定部ベイズ推定部102は、輝度色差分離部101から入力する入力画像ジャギー発生画像10に基づいて生成された輝度画像から画素ブロックを選択し、上記式4を適用して、ジャギー発生前の画像における画素値輝度値:xを推定する。この推定処理によって、ジャギー発生前推定画像、すなわちジャギー発生前の輝度画像が生成される。なお、ここで推定した輝度画像は階調数が少なく、人工的な見た目となる傾向が強い。これは確率的に尤もらしい画素値が、全階調のうちの一部に偏っているためである。そこで、本手法の一例ではベイズ推定によって求めた輝度情報はエッジの形状情報としてのみ利用することとする。上述の画像推定部ベイズ推定部102におけるベイズ推定処理では、注目画素とその近傍画素の画素値を利用して、ジャギー発生後画像の画素値から、ジャギー発生前の画素値を推定する処理を行っている。この処理例以外にも、例えばエッジを表現する特徴量を用いたベイズ推定により、ジャギー発生後画像の画素値から、ジャギー発生前の画素値を推定する処理を行うことも可能である。この場合、事前学習データとして、事前にジャギー発生前後の画像からエッジを表現する特徴量を抽出して統計量を保持しておく。ジャギー低減処理の対象となるジャギー発生画像が入力された時に、ジャギー発生画像からエッジを表現する特徴量を抽出してベイズ推定によってジャギー発生前のエッジ特徴量を求める。この特徴量を利用するとジャギー発生前のエッジ間の類似度を求めることができるため、これに続く平均化処理で必要な類似度として利用することができる。画像推定部ベイズ推定部102では、ジャギー発生前のエッジ形状を正確に求めるために確率的なアプローチを採用しているが、ジャギー発生画像に対して、例えば平滑化フィルタを施して平滑化した画像を生成し、平滑化画像をジャギー発生前の画像とみなす構成としてもよい。ジャギーが見かけ上すべて低減される程度の平滑化フィルタを画像全体にかけると画像全体の解像度が劣化するが、ジャギーが特に顕著に現れるコントラストの強いエッジの構造自体は維持される。この構造を利用してエッジ近辺の類似度を算出することでベイズ推定を用いる方法の代替手法となる。具体的にはジャギー発生後の画像に対して、例えば、次の式5に示すような線形フィルタを適用して平滑化を行う。式5上記式5に示すフィルタのカットオフ周波数はジャギーの原因となる折り返し発生周波数でもある。すでに折り返しが発生している画像に対する平滑化のため、折り返し自身を完全に取り除くことは困難であるが、おおよそジャギー低減効果を得られるため、このようなフィルタを利用する。2-3.処理3推定された輝度画像から類似度を求め、NL-Meansによってジャギー低減を行う処理次に、画像推定部ベイズ推定部102において推定されたジャギー発生前推定輝度画像から類似度を求め、NL-MeansNonLocal-Means:ノンローカル手法によってジャギー低減を行う処理について説明する。この処理は、主に図4に示す画像処理装置100中の加重加算部NL-Means103において実行される。画像処理装置100において低減対象としているジャギーは、多くの場合、エッジに沿って規則的なパターンとなって出現する。そのため、類似エッジ周辺画素を平均化することでジャギーをキャンセルすることができる。ただし、ジャギーが残ったままのエッジ形状ではエッジに沿って平均化することが非常に困難である。従って、上述した処理2で説明したエッジ形状推定結果を利用することにより、ジャギーの影響を無視して類似エッジ周辺画素の平均化を行う。本手法では平均化処理にNL-Means法を用いる。NL-Means法による平均化処理は、以下の式6に従った処理として実行される。式6なお、各パラメータは以下のパラメータである。x:対象画素値のx座標y:参照画素値のy座標I:画像全体の画素の集合gx:NL-Means適用後の画素値fx:NL-Means適用前の画素値fx:座標xの近傍画素値をベクトルで表現したものwx,y:座標x,yにある画素値から得られる重みh:NL-Meansの強度調整パラメータCx:重みの総和上記の式6に示すように、NL-Meansは、対象画素と、対象画素の近傍領域と類似している領域の中心画素値を、その類似度に応じて決まる重みを乗算して加算する、という手法である。加重加算部NL-Means103は、画像推定部ベイズ推定部102において生成した推定輝度画像83を入力して、この推定輝度画像83に基づいて類似度を判定し、加重平均の重み算出に用いる。この処理によって、推定輝度画像83として生成されるジャギー発生前画像において、類似するエッジ領域には大きめの重みをかけて加重平均を行うことができる。これにより規則的なパターンとして出現するジャギーの影響を低減することが可能となる。図6を参照して、加重加算部NL-Means103の実行する処理の具体例について説明する。図6aは、画像推定部ベイズ推定部102の生成したジャギー発生前の推定輝度画像である。例えば、加重加算部NL-Means103は、この推定輝度画像から、ジャギー低減対象領域としての処理領域201を選択し、この処理領域201に対する類似領域を画像内から検出する。例えば図6bに示すように3つの類似領域211~213が検出されたものとする。次に、加重加算部NL-Means103は、これらの3つの類似領域211~213と、処理領域の類似度を判定する。例えば画素値の変化パターン等に応じた類似度を各類似領域211~213に設定する。さらに、各類似領域211~213に設定した類似度に応じた重みを設定する。処理領域201との類似度の高い類似領域には大きめの重みを設定し、類似度の低い類似領域には小さめの重みを設定する。次に、加重加算部NL-Means103は、図4に示す輝度色差分離部101の出力するジャギー低減前の輝度画像81を入力して、この輝度画像に対する加重加算処理を行う。図4に示す輝度色差分離部101の出力するジャギー低減前の輝度画像81は、図6cに示す画像である。加重加算部NL-Means103は、図6bに示す画像推定部ベイズ推定部102の生成したジャギー発生前の推定輝度画像に基づいて設定した処理領域201と検出した類似領域211~213に対応する位置の領域を図6cに示すジャギー低減前の輝度画像81から選択して、各対応する類似領域の位置の画素に対して、算出済みの重みを用いて加重平均を行い、処理領域の画素値輝度値の補正画素値補正輝度値を算出する。この処理によりジャギーの低減された輝度画像を生成することができる。この結果画像が、図6dに示すジャギー低減輝度画像である。なお、上記処理例では、加重加算部NL-Means103は、NL-Meansノンローカル手法によってジャギー低減を行う構成として説明したが、NL-Meansの代わりにJointBilateralFilterを利用することも可能である。NL-Meansはブロック単位で領域間の類似度を算出するが、ブロックの代わりに画素単位の類似度を利用することも可能である。JointBilateralFilterを適用する場合、画素単位の類似度算出に基づく処理も可能である。2-4.処理4処理3によるジャギー低減画像のエッジ領域と、原画像のエッジ領域以外の部分をブレンドする処理次に、処理3によるジャギー低減画像のエッジ領域と、原画像のエッジ領域以外の部分をブレンドする処理について説明する。この処理は、図1に示す画像処理装置100のエッジ抽出部104と、ブレンド処理部105の処理として実効される。前述の処理3、すなわち、加重加算部103の実行する推定輝度画像から類似度を求め、NL-Meansによってジャギー低減を行う処理によってエッジ近辺のジャギーを低減することができる。しかし、画像内のすべての領域でNL-Meansを適用するとテクスチャ領域などでどうしても画像の解像度劣化が顕著となってしまう。これはテクスチャのようなランダム性の強い信号を多く含む領域では、テクスチャそのものがノイズのように扱われてNL-Meansによって信号レベルが低減されてしまうためであり、NL-Meansの原理的な問題である。そこで、処理対象画像内のジャギーの多く発生するエッジ近辺にはNL-Meansの結果を積極的に利用し、それ以外の領域は原画像を積極的に用いるようなブレンド処理マスク処理を行う。すなわち、NL-Means適用結果を利用したいのはエッジ近辺のみであり、通常のエッジ検出フィルタで反応するようなテクスチャ成分は原画像を利用するという画像特徴に応じて、それぞれの出力をブレンドする処理を行う。このブレンド処理を実行するのが、図1に示す画像処理装置100のブレンド処理部105である。エッジ検出処理は、エッジ抽出部LaplacianofGaussian104において実行される。まず、エッジ抽出部104の処理について説明する。図4に示す画像処理装置100のエッジ抽出部104は、LaplacianofGaussianLoGフィルタによってエッジ領域を抽出する。LoGフィルタとはガウシアンによって平滑化をかけた画像に対して、ラプラシアンによるエッジ検出を行うフィルタであり、次の式7によって示される。式7ただし、i,j:フィルタの中心からの水平方向垂直方向の距離σ:ガウシアンの分散値である。上記式7に従って、テクスチャ領域は最初のガウシアンフィルタによって高周波成分が除去された状態でラプラシアンフィルタをかけることになるため、そのままラプラシアンフィルタをかける場合よりも検出感度を抑えることができる。一方、ジャギーが特に目立つようなコントラストの強いエッジでは、ガウシアンによる平滑化がかかっても高周波成分を依然として保持しているため、ラプラシアンによる反応感度を高いレベルで維持することができる。エッジ抽出部104は、上記式7に示すフィルタによって得られた検出値LoGフィルタ検出値を、図4に示す輝度高周波成分情報84としてブレンド処理部105に出力する。ブレンド処理部105は、エッジ抽出部104から入力する輝度高周波成分情報84に対して、予め既定した上限閾値thupperと下限閾値thlowerを利用して閾値間を直線的に変化するような重みブレンド率分布関数に基づいて、各画素単位の輝度高周波成分情報に応じたブレンド率を算出する。ブレンド処理部105は、図4に示す輝度色差分離部101の出力するジャギー補正前の輝度画像L,81と、加重加算部NL-Means103の出力するNL-Meansの適用結果であるジャギー補正輝度画像Lp,85とのブレンド処理を実行してブレンド済みジャギー補正輝度画像Lq,86を生成する。上記のブレンド率は、ブレンド処理部105の生成するブレンド済みジャギー補正輝度画像Lq,86の各画素における、加重加算部NL-Means103の出力するジャギー補正輝度画像Lp,85のブレンド率に相当する。次式式8はブレンド率分布関数αlを示し、図7にその例を示す。式8なお、上記式8において、各パラメータは以下のパラメータである。l:LoGフィルタ検出値αl:LoGフィルタ検出値がlの時のNL-Means結果ブレンド率thupper,thlower:それぞれ上限、下限閾値図4に示すブレンド処理部105は、上記のブレンド率分布関数によって得られたブレンド率を適用して、図4に示す輝度色差分離部101の出力するジャギー補正前の輝度画像L,81=foと、加重加算部NL-Means103の出力するNL-Meansの適用結果であるジャギー補正輝度画像Lp,85=fNとのブレンド処理を実行してブレンド済みジャギー補正輝度画像Lq,86を生成する。ブレンド処理部105におけるブレンド処理は、以下の式9に従って実行する。式9ただし、fx:ブレンド後の輝度画像=図4に示すブレンド済みジャギー補正輝度画像Lq,86αx:座標xでのブレンド率fNx:ジャギー補正輝度画像Lp,85の画素値輝度値fox:ジャギー補正前の輝度画像L,81の画素値輝度値である。以上により、最終出力となるジャギー補正画像の輝度成分からなる輝度画像が生成される。すなわち、図4に示すブレンド済みジャギー補正輝度画像Lq,86である。ブレンド処理部105は、生成したジャギー補正輝度画像Lq,86をフィルタ処理部106と、輝度色差合成部107に出力する。2-5.処理5色差画像に対してジャギー低減済み輝度画像を参照するJBFによる相関処理に基づく補正色差情報の生成次に、色差画像に対してジャギー低減済み輝度画像を参照するJBFによる相関処理に基づく補正色差情報の生成処理について説明する。この処理は、図4に示す画像処理装置100のフィルタ処理部106において実行される。フィルタ処理部106は、ブレンド処理部105の生成したブレンド済みジャギー補正輝度画像Lq,86を入力し、入力したブレンド済みジャギー補正輝度画像Lq,86を参照画像として利用したJointBilateralFilterJBFを適用した相関処理により補正色差情報Cr,Cbを生成する。JBFは以下の式10によって表される。式10ただし、上記式10において、各パラメータは以下のパラメータである。x:対象画素値の座標y:参照画素値の座標N:対象画素xの近傍画素の集合gx:JBF適用後の画素値fx:JBF適用前の画素値rx:JBFで類似度を算出する参照画像wx,y:座標x,yにある画素値から得られる重みh:JBFの強度調整パラメータCx:重みの総和なお、輝度成分を参照画像として色差成分にJBFを適用する処理は、輝度と色差の相関を取る処理に他ならない。既に輝度成分は前述した各種の処理によりエッジ近辺のジャギーが既に低減されている状態である。この状態の輝度との相関を取ることで色差成分の補正を行う。この結果、補正された輝度画素に応じて各画素の色差成分が補正される。フィルタ処理部106の生成する補正色差Cr,Cbは、輝度色差合成部107に出力される。2-6.処理6輝度色差合成処理次に、図4に示す画像処理装置100の輝度色差合成部107の実行する輝度色差合成処理について説明する。輝度色差合成部107は、ブレンド処理部105の出力するブレンド済みジャギー補正輝度画像Lq,86、フィルタ処理部106の出力する補正色差Cr,Cb、これらを入力する。輝度色差合成部107は、これらの輝度色差成分からRGB変換への逆変換処理を実行してジャギー低減画像を出力画像90として生成する。なお、輝度色差成分からRGB変換への逆変換処理は、先に処理1RGB画像から輝度画像とその色差情報への色変換処理において説明した式1に、ブレンド処理部105の出力するブレンド済みジャギー補正輝度画像Lq,86、の輝度値と、フィルタ処理部106の出力する補正色差Cr,Cbを入力して、R,G,Bを求める処理として実効することができる。輝度色差合成部107は、この逆変換処理によって、ジャギー低減画像としての出力画像90を生成する。なお、上述した実施例において、色差成分に前述したジャギー低減を行わないのは主に次の2つの理由からである。第一に自然画において人間が解像度劣化を強く認知するのは輝度成分の劣化によるものであり、色差成分の高周波成分劣化は輝度成分ほどの影響を与えない。そのため、輝度成分には計算時間をかけてでも解像度劣化を極力抑えたジャギー低減処理を行い、色差成分には多少の解像度劣化が発生してしまうが比較的計算量が少ないJBFを利用している。第二に、画素加算間引きによって発生するアーチファクトは輝度成分のジャギーだけでなく、偽色も同様に観測される。これは画素加算間引きによるリサンプリングにより、R,G,B成分の高周波成分にそれぞれの位相で折り返しが発生してしまうためである。偽色の低減で有効なのは、すべての色成分を輝度成分と相関を持つように補正することである。輝度成分を参照画像としたJBFには強い相関を与える効果があるため、色差にJBFを適用することで輝度のジャギーだけでなく偽色対策も行うことが可能である。以上、本開示の画像処理装置の実行する一実施例について説明してきた。なお、本開示の画像処理装置の実行するジャギー低減処理には、例えば以下の特徴がある。任意のジャギーパターンに対応可能解像度劣化が小さい線形フィルタで表現できないような角度のエッジでもジャギー低減可能偽色低減を行える[3.画像処理装置の構成例について]先に、図4を参照して画像処理装置100の構成について説明したが、前述したように、図4に示す画像処理装置100は、画像処理装置の主要部の構成であり、画像処理装置100は、この他、例えば全体の処理制御を実行する制御部、例えば所定の処理プログラムの実行機能持つCPUを備えた制御部や、プログラム、画像、各種のパラメータを記憶するためのメモリを備えている。画像処理装置100は、具体的には、例えば、PC、撮像装置カメラとして構成することが可能である。図4に示す画像処理構成を撮像装置に組み込んだは場合の構成例を図8に示す。図8は、撮像装置の構成例を示すブロック図である。光学レンズ301を介して入射される光は例えばCMOSイメージセンサなどによって構成される撮像素子302に入射し、光電変換による画像データを出力する。出力画像データは信号処理部303に入力される。信号処理部303は、例えばホワイトバランスWB調整、ガンマ補正、デモザイク処理等、一般的なカメラにおける信号処理を行う。さらに、信号処理部303は、例えば、画素加算処理や画素間引き処理等の画素数の削減処理と、上述したジャギー低減処理を実行して出力画像320を生成する。出力画像320は、例えば図示しない記憶部に格納される。あるいは表示部に出力される。制御部305は、例えば図示しないメモリに格納されたプログラムに従って各部に制御信号を出力し、各種の処理の制御を行う。具体的には、例えば図4~図7を参照して説明したジャギー低減処理、すなわち、処理1RGB画像から輝度画像とその色差へと色変換を行う。またこの逆変換を行う。[輝度色差分離部101、輝度色差合成部107の処理]処理2ジャギーの発生している輝度画像からベイズ推定によってジャギー発生前の輝度画像を推定する。[画像推定部ベイズ推定部102の処理]処理3推定された輝度画像から類似度を求め、NL-Meansによってジャギー低減を行う。[主に加重加算部103の処理]処理4処理3によるジャギー低減画像のエッジ領域と、原画像のエッジ領域以外の部分をブレンドする。[エッジ抽出部104と、ブレンド処理部105の処理]処理5色差画像に対してジャギー低減済み輝度画像を参照するJBFによっての相関処理に基づく補正色差情報の生成。[フィルタ処理部106の処理]処理6輝度色差合成処理。[輝度色差合成部107の処理]これらの処理に対する制御等を行う。[4.本開示の構成のまとめ]以上、特定の実施例を参照しながら、本開示の実施例について詳解してきた。しかしながら、本開示の要旨を逸脱しない範囲で当業者が実施例の修正や代用を成し得ることは自明である。すなわち、例示という形態で本発明を開示してきたのであり、限定的に解釈されるべきではない。本開示の要旨を判断するためには、特許請求の範囲の欄を参酌すべきである。なお、本明細書において開示した技術は、以下のような構成をとることができる。1ジャギー発生画像からジャギー発生前の画像を推定して、ジャギー発生前推定画像を生成する画像推定部と、前記ジャギー発生前推定画像におけるシャギー発生領域を処理領域として選択し、該処理領域に類似する画素領域である類似領域を検出し、検出した類似領域各々に対して処理領域との類似度に応じた重みを算出し、前記処理領域と類似領域に対応する前記ジャギー発生画像の対応領域を検出し、検出した対応領域の画素値の前記重みを適用した加重加算処理により、前記ジャギー発生画像の処理領域の補正画素値を算出する加重加算部を有する画像処理装置。2前記画像推定部は、前記ジャギー発生画像から生成した輝度画像を入力してジャギー発生前推定輝度画像を生成し、前記加重加算部は、前記ジャギー発生前推定輝度画像に基づいて前記重みを算出し、前記ジャギー発生画像から生成した輝度画像を適用して、前記重みを適用した画素値の加重加算処理を実行する前記1に記載の画像処理装置。3前記画像推定部は、事前学習データを適用したベイズ推定処理により、ジャギー発生前推定画像を生成する前記1または2に記載の画像処理装置。4前記画像推定部は、前記ジャギー発生画像の平滑化処理によってジャギー発生前推定画像を生成する前記1~3いずれかに記載の画像処理装置。5前記画像推定部は、エッジを表現する特徴量の統計データを事前学習データとして適用したベイズ推定処理により、ジャギー発生前推定画像を生成する前記1~4いずれかに記載の画像処理装置。6前記加重加算部は、NL-MeansNonLocal-Means法を適用した加重加算処理を実行する前記1~5いずれかに記載の画像処理装置。7前記加重加算部は、JointBirateralFilter法を適用した加重加算処理を実行する前記1~5いずれかに記載の画像処理装置。8前記画像処理装置は、さらに、前記加重加算部の算出した補正画素値を持つジャギー低減輝度画像を参照画像として、JointBilateralFilterJBFを適用した相関処理を実行して補正色差情報を生成するフィルタ処理部を有する前記1~7いずれかに記載の画像処理装置。9前記画像処理装置は、さらに、前記加重加算部の算出した補正画素値を持つジャギー低減輝度画像と、補正前のジャギー発生画像に基づいて生成されるジャギー発生輝度画像とを、画像のエッジ領域に応じた重みに従ってブレンドするブレンド処理部を有する前記1~8いずれかに記載の画像処理装置。また、明細書中において説明した一連の処理はハードウェア、またはソフトウェア、あるいは両者の複合構成によって実行することが可能である。ソフトウェアによる処理を実行する場合は、処理シーケンスを記録したプログラムを、専用のハードウェアに組み込まれたコンピュータ内のメモリにインストールして実行させるか、あるいは、各種処理が実行可能な汎用コンピュータにプログラムをインストールして実行させることが可能である。例えば、プログラムは記録媒体に予め記録しておくことができる。記録媒体からコンピュータにインストールする他、LANLocalAreaNetwork、インターネットといったネットワークを介してプログラムを受信し、内蔵するハードディスク等の記録媒体にインストールすることができる。なお、明細書に記載された各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。また、本明細書においてシステムとは、複数の装置の論理的集合構成であり、各構成の装置が同一筐体内にあるものには限らない。"