(57)【要約】【課題】測定装置固有のノイズと、粗度形状とが同一領域に混在する場合にも、粗度断面曲線から、測定装置固有のノイズのみを除去する粗度断面曲線評価方法、粗度解析装置及び粗度解析プログラム、並びに、このような粗度断面曲線評価方法を採用する粗度計測方法及び粗度計測装置を提供する。【解決手段】測定対象の測定断面曲線をフーリエ変換することによって、測定対象の粗度スペクトルを生成し、予め測定した基準面の測定断面曲線をフーリエ変換することによって、基準面の粗度スペクトルを生成し、測定対象の粗度スペクトルから、基準面の粗度スペクトルを減算することによって、測定対象の粗度スペクトルからノイズ除去を行い、ノイズ除去後の測定対象の粗度スペクトルを逆フーリエ変換することによって、測定対象の断面曲線を生成する。【選択図】図１
2020134356.tif 000002
【特許請求の範囲】
【請求項１】
測定対象の測定断面曲線から該測定対象の表面粗度評価を行う粗度断面曲線評価方法であって、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを生成し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを生成し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を生成する、粗度断面曲線評価方法。
【請求項２】
前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出する、請求項１に記載の粗度断面曲線評価方法。
【請求項３】
請求項１または２に記載の粗度断面曲線評価方法をコンピュータにより実行させるための粗度解析プログラム。
【請求項４】
請求項３に記載の粗度解析プログラムを記録したコンピュータ可読記録媒体。
【請求項５】
測定対象の測定断面曲線から該測定対象の表面粗度評価を行う粗度解析装置であって、
前記測定対象の測定断面曲線を特定可能な測定断面曲線特定データをデータ入力可能な入力手段と、
前記測定断面曲線特定データに基づき、前記測定対象の測定断面曲線を生成するとともに、前記測定対象の測定断面曲線から表面粗度パラメーターを算出する演算装置と、
を備え、
前記演算装置が、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを取得し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを取得し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を取得するように構成される、粗度解析装置。
【請求項６】
前記演算装置が、
前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出するように構成される、請求項５に記載の粗度解析装置。
【請求項７】
測定対象の表面粗度を計測する粗度計測方法であって、
変位計によって、Ｘ軸方向の一定間隔毎に変位を測定することにより変位データを取得することによって、前記測定対象の測定断面曲線を生成し、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを生成し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを生成し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を生成する、粗度計測方法。
【請求項８】
前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出する、請求項７に記載の粗度計測方法。
【請求項９】
変位計と、
前記変位計をＸ軸方向に移動させる移動機構と、
前記変位計のＸ軸方向の移動距離を読み取る移動距離計測装置と、
前記変位計により取得された変位データと、前記移動距離計測装置により取得された移動距離データとに基づき、測定対象の測定断面曲線を生成するとともに、前記測定対象の測定断面曲線から表面粗度パラメーターを算出する演算装置と、
を備える粗度計測装置であって、
前記演算装置が、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを取得し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを取得し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を取得するように構成される、粗度計測装置。
【請求項１０】
前記演算装置が、
前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出するように構成される、請求項９に記載の粗度計測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、物体や構造物などの表面粗度を評価するために測定装置によって取得される粗度断面曲線から、測定装置固有のノイズを除去する粗度断面曲線評価方法、粗度解析装置及び粗度解析プログラム、並びに、このような粗度断面曲線評価方法を採用する粗度計測方法及び粗度計測装置に関する。
【背景技術】
【０００２】
物体の表面粗度形状は、その表面の摩擦特性を把握する上で重要である。近年では、非特許文献１に示すように、物体表面の水流摩擦特性の把握には、粗度高さだけではなく、波長等を含めた形状パラメーターの測定の重要性が示されている。
【０００３】
非特許文献１では、船体の表面粗度（塗膜表面粗度）に関する評価について述べられており、船体塗膜表面粗度の評価には、波長の評価が重要であることが示されている。また、粗度のパラメーターを評価する場合、その粗度の波長に対して十分な水平方向の解像度を有する測定が必要となる。
【０００４】
しかしながら、水平方向の解像度を上げた場合、測定装置固有のノイズが粗度断面曲線に現れる場合がある。測定装置固有のノイズの原因としては、例えば、触針や光プローブ（レーザー）の機械的振動、電気信号に乗る熱雑音などの電気的ノイズ、レーザー光の散乱などによるものが考えられる。
【０００５】
このような測定装置固有のノイズを除去する方法として、例えば、非特許文献２のように短波長λSを非通過とするフィルタを用いる方法がある。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】三重野  紘央、増田  宏、「船底塗料の塗膜表面粗度による抵抗増加について-船底外板塗料による船体抵抗低減」、日本マリンエンジニアリング学会誌、第48巻、第 3 号（2013）、p300-307
【非特許文献２】ISO 4287、 「Geometrical Product specifications (GPS)-Surface texture : Profile method: Terms, definitions and surface texture parameters.」、 (1997)
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、非特許文献２に開示されている手法では、測定装置固有のノイズと、粗度形状の波長とが、同一領域にあり混在している場合には、短波長λSを非通過とするフィルタの適用により、粗度形状も同時に除去してしまうことになり、正確な粗度評価が行えないという問題がある。
【０００８】
本発明では、このような現状に鑑み、測定装置固有のノイズと、粗度形状とが同一領域に混在する場合にも、粗度断面曲線から、測定装置固有のノイズのみを除去する粗度断面曲線評価方法、粗度解析装置及び粗度解析プログラム、並びに、このような粗度断面曲線評価方法を採用する粗度計測方法及び粗度計測装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、前述するような従来技術における課題を解決するために発明されたものであって、本発明は、例えば、以下の態様を含む。
【００１０】
［１］  測定対象の測定断面曲線から該測定対象の表面粗度評価を行う粗度断面曲線評価方法であって、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを生成し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを生成し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を生成する、粗度断面曲線評価方法。
【００１１】
［２］  前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出する、項［１］に記載の粗度断面曲線評価方法。
【００１２】
［３］  項［１］または［２］に記載の粗度断面曲線評価方法をコンピュータにより実行させるための粗度解析プログラム。
【００１３】
［４］  項［３］に記載の粗度解析プログラムを記録したコンピュータ可読記録媒体。
【００１４】
［５］  測定対象の測定断面曲線から該測定対象の表面粗度評価を行う粗度解析装置であって、
前記測定対象の測定断面曲線を特定可能な測定断面曲線特定データをデータ入力可能な入力手段と、
前記測定断面曲線特定データに基づき、前記測定対象の測定断面曲線を生成するとともに、前記測定対象の測定断面曲線から表面粗度パラメーターを算出する演算装置と、
を備え、
前記演算装置が、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを取得し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを取得し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を取得するように構成される、粗度解析装置。
【００１５】
［６］  前記演算装置が、
前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出するように構成される、項［５］に記載の粗度解析装置。
【００１６】
［７］  測定対象の表面粗度を計測する粗度計測方法であって、
変位計によって、Ｘ軸方向の一定間隔毎に変位を測定することにより変位データを取得することによって、前記測定対象の測定断面曲線を生成し、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを生成し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを生成し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を生成する、粗度計測方法。
【００１７】
［８］  前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出する、項［７］に記載の粗度計測方法。
【００１８】
［９］  変位計と、
前記変位計をＸ軸方向に移動させる移動機構と、
前記変位計のＸ軸方向の移動距離を読み取る移動距離計測装置と、
前記変位計により取得された変位データと、前記移動距離計測装置により取得された移動距離データとに基づき、測定対象の測定断面曲線を生成するとともに、前記測定対象の測定断面曲線から表面粗度パラメーターを算出する演算装置と、
を備える粗度計測装置であって、
前記演算装置が、
前記測定対象の測定断面曲線をフーリエ変換することによって、前記測定対象の粗度スペクトルを取得し、
予め測定した基準面の測定断面曲線をフーリエ変換することによって、前記基準面の粗度スペクトルを取得し、
前記測定対象の粗度スペクトルから、前記基準面の粗度スペクトルを減算することによって、前記測定対象の粗度スペクトルからノイズ除去を行い、
ノイズ除去後の前記測定対象の粗度スペクトルを逆フーリエ変換することによって、前記測定対象の断面曲線を取得するように構成される、粗度計測装置。
【００１９】
［１０］  前記演算装置が、
前記測定対象の断面曲線について、パラメーター解析することによって、前記測定対象の表面粗度パラメーターを算出するように構成される、項［９］に記載の粗度計測装置。
【発明の効果】
【００２０】
本発明によれば、測定された断面曲線をフーリエ変換することによって、ノイズの特性が現れやすい周波数領域においてノイズ除去を行っているため、測定対象である粗面特有の要素を除去することなく、正確な粗度評価を行うことができる。
【００２１】
これにより、例えば、移動機構の偏心によって生じる長波長側の影響など、粗度測定装置固有のノイズを除去することができ、粗度測定装置の個体差によらず、正確な粗度評価を行うことができる。
【図面の簡単な説明】
【００２２】
【図１】図１は、本発明の粗度解析装置を含む粗度計測装置の一例を説明するための概略構成図であり、図１（ａ）は概略側面図、図１（ｂ）は概略正面図である。
【図２】図２は、本実施形態の粗度計測装置を用いて粗度断面曲線を評価する流れを説明するフローチャートである。
【図３】図３は、実施例における粗面Ａ〜Ｄ及び基準面について取得された測定断面曲線における粗度スペクトルを示すグラフである。
【図４】図４は、実施例における粗面Ａ〜Ｄのノイズ除去後の粗度スペクトルを示すグラフである。
【図５】図５は、比較例における粗面Ａ〜Ｄの短波長ガウシアンフィルタ適用後の粗度スペクトルを示すグラフである。
【発明を実施するための形態】
【００２３】
以下、本発明の実施の形態（実施例）を図面に基づいて、より詳細に説明する。
図１は、本発明の粗度解析装置を含む粗度計測装置の一例を説明するための概略構成図であり、図１（ａ）は概略側面図、図１（ｂ）は概略正面図である。
【００２４】
図１に示すように、本実施形態の粗度計測装置１０は、二次元レーザー変位計１２と、移動機構１４と、移動距離計測装置１６と、演算装置１８と、二次元レーザー変位計１２，移動距離計測装置１６，演算装置１８を動作させるための可搬型バッテリー２０と、を備えている。
【００２５】
二次元レーザー変位計１２としては、特に限定されるものではないが、例えば、キーエンス社製ＬＪ−Ｖ７０８０（基準距離におけるレーザー幅３２ｍｍ）などを用いることができる。なお、本実施形態では、変位計として二次元レーザー変位計１２を用いているが、本発明における変位計はこれに限定されることはなく、例えば、触針型変位計であってもよく、また、二次元レーザー変位計１２以外のレーザー型変位計などであってもよい。
【００２６】
また、二次元レーザー変位計１２は、一定間隔毎にＹ軸方向の座標の変位データを一度に読み取れるように、二次元レーザー変位計１２の幅方向が、移動機構１４の移動方向（Ｘ軸方向）と垂直、すなわち、Ｙ軸方向と一致するように配置されている。
【００２７】
また、本実施例において、移動機構１４は、所定方向（Ｘ軸方向）のみに移動可能な二対の車輪（すなわち、四輪）としているが、所定方向のみに移動可能な機構であれば特に限定されるものではなく、例えば、クローラーなどであってもよい。
【００２８】
また、本実施例では、移動距離計測装置１６として、ロータリーエンコーダーを用いており、移動機構１４である車輪の回転数を読み取ることにより、車輪の周長と回転数とに基づいて移動距離を算出している。
【００２９】
なお、移動距離計測装置１６は、これに限らず、例えば、移動機構１４と連動した円形スケールにより移動距離をスケールセンサーで読み取るような装置とすることもできる。
【００３０】
また、演算装置１８は、後述するように、二次元レーザー変位計１２により取得される変位データや、移動距離計測装置１６により取得される移動距離データなどの測定断面曲線を特定可能なデータ（本明細書においては、「測定断面曲線特定データ」と言う。）を入力可能なデータ入力手段１９を備えるとともに、このような測定断面曲線特定データに基づき、表面粗度パラメーターなどを算出可能なものであれば特に限定されるものではない。例えば、タブレットコンピューターやスマートフォン、パーソナルコンピューターなどの外部端末に変位データと移動距離データとを送信して、この外部端末を演算装置１８として演算処理を実施し、演算結果を外部端末に表示するように構成することができる。
【００３１】
なお、演算装置１８は、粗度計測装置１０に内蔵することができる。この場合、粗度計測装置１０に演算結果を表示する表示手段を設けてもよいし、データ出力手段を設けて、演算結果を上述するような外部端末に表示するように構成することもできる。
【００３２】
演算装置１８には、後述するように動作する粗度解析プログラムがインストールされており、粗度解析装置として動作する。
【００３３】
また、可搬型バッテリー２０としては、特に限定されるものではなく、例えば、鉛蓄電池やニッケル水素二次電池、リチウムイオン二次電池など、既存の二次電池を用いることもできる。
【００３４】
なお、本実施形態において、可搬型バッテリー２０を搭載することにより、持ち運び容易な粗度計測装置１０としているが、二次元レーザー変位計１２や演算装置１８などの電源を、例えば、外部バッテリーや商用電源などから供給するように構成することも可能である。
【００３５】
このように構成される本実施形態の粗度計測装置１０は、以下のようにして、測定対象である粗面の断面曲線、及び、この断面曲線の解析によって得られる粗度パラメーターを算出することができる。
【００３６】
図２は、本実施形態の粗度計測装置を用いて粗度断面曲線を評価する流れを説明するフローチャートである。
【００３７】
まず、粗度計測装置１０を移動機構１４により基準面の表面上で移動させることにより、一定間隔毎に移動距離計測装置１６により得られる移動距離に基づく移動距離データと、二次元レーザー変位計１２により得られる変位データとを演算装置１８に記録する。これにより、基準面の測定断面曲線が得られる（Ｓ１）。
【００３８】
本明細書において、「基準面」とは、理想的には、一切の凹凸がない平滑面のことをいい、本実施形態においては、例えば、鏡面仕上げがなされた樹脂板や金属板など、表面の平滑度が高い平板を用いることができる。このような基準面を用いて取得された測定断面曲線には、粗度計測装置１０固有のノイズのみが含まれている。
【００３９】
また、基準面の測定断面曲線は、測定対象２２の測定の都度、取得する必要はなく、予め取得しておき、演算装置１８に記憶しておけばよい。なお、二次元レーザー変位計１２の経時変化などを考慮すると、基準面の測定断面曲線を定期的に取得し、更新することが好ましい。
【００４０】
次いで、粗度計測装置１０を移動機構１４により測定対象２２の表面上で移動させることにより、一定間隔毎に移動距離計測装置１６により得られる移動距離に基づく移動距離データと、二次元レーザー変位計１２により得られる変位データとを演算装置１８に記録する。これにより、測定対象２２の測定断面曲線が得られる（Ｓ２）。
【００４１】
演算装置１８は、このようにして得られた基準面の測定断面曲線及び測定対象２２の測定断面曲線について、それぞれフーリエ変換し、波長成分に分解する。これにより、各断面曲線を、波長λに対する振幅Ｈの関数として表す粗度スペクトルとすることができる（Ｓ３）。なお、演算装置１８によるフーリエ変換としては、離散フーリエ変換（ＤＦＴ）や高速フーリエ変換（ＦＦＴ）などを用いることができる。
【００４２】
次いで、演算装置１８は、測定対象２２の粗度スペクトルの各波長での振幅から、基準面の粗度スペクトルの各波長での振幅を減算する。
これにより、測定対象２２の粗度スペクトルから、粗度計測装置１０固有のノイズが除去されることになる（Ｓ４）。
【００４３】
次いで、演算装置１８は、ノイズ除去後の測定対象２２の粗度スペクトルについて、逆フーリエ変換することによって、断面曲線を得ることができる（Ｓ５）。なお、演算装置１８による逆フーリエ変換としては、逆離散フーリエ変換（ＩＤＦＴ）や逆高速フーリエ変換（ＩＦＦＴ）などを用いることができる。
【００４４】
そして、演算装置１８は、この断面曲線をパラメーター解析することによって、例えば、Ｒｚ（最大高さ粗さ）、ＲＳｍ（要素の平均長さ）、ＲＺＪＩＳ（十点平均粗さ）、Ｒａ（算術平均粗さ）、Ｒｑ（二乗平均平方根粗さ）、Ｒｃ（粗さ曲線の平均高さ）、Ｒｓｋ（スキューネス）、Ｒｋｕ（クルトシス）などの表面粗度パラメーターを算出する（Ｓ６）。
【００４５】
断面曲線から表面粗度パラメーターを算出する方法は、例えば、ＩＳＯ４２８７：１９９７やＩＳＯ４２８８：１９９６などで規定された方法を用いることができる。
【００４６】
なお、本発明の粗度解析装置を含む粗度計測装置、及び粗度断面曲線評価方法において、変位データ及び移動距離データの測定間隔は、いずれも制限されるものではないが、計測精度と計測時間のバランスの観点から、１０μｍ〜１００μｍとすることが好ましく、３０μｍ〜８０μｍとすることがより好ましい。
【００４７】
［実施例］
このように構成される本実施形態の粗度計測装置１０による測定対象２２の表面粗度評価の具体例を以下に示す。
【００４８】
本実施例においては、二次元レーザー変位計１２のレーザーの幅方向（Ｙ軸方向）に５０μｍ間隔、進行方向（Ｘ軸方向）に６２．５μｍ間隔で、幅方向に３０ｍｍ、進行方向に３２ｍｍの範囲（Ｘ軸方向に５１３点、Ｙ軸方向に６０１点）の変位データを各測定対象について測定し、断面曲線を取得した。
【００４９】
測定対象として、粗面Ａ，粗面Ｂ，粗面Ｃ，粗面Ｄの４つの粗面を測定した。また、同様の条件で、平滑面についても測定し、これを基準面とした。
【００５０】
粗面Ａ〜Ｄ及び基準面について取得された測定断面曲線をそれぞれフーリエ変換し、波長成分に分解し、粗度スペクトルを取得した。図３は、粗面Ａ〜Ｄ及び基準面について取得された測定断面曲線における粗度スペクトルを示すグラフである。なお、図３では、粗度スペクトルのピークを詳細に確認するため、縦軸をＨ／λとしている。
【００５１】
図３に示すように、粗面Ａ〜Ｄの粗度スペクトルには、各粗面特有のピークが現れているが、短波長側に粗面だけでなく基準面にも一様に現れるピークが確認できる。これは、全測定対象に一様に現れていることから、粗度計測装置１０固有のノイズに起因するものであると考えられる。
【００５２】
短波長側において、粗面のピークと粗度計測装置１０固有のノイズとが混ざり合っているものもあり、この状態のまま、各粗面の測定断面曲線のパラメーター解析を行い、表面粗度パラメーターを算出したとしても、粗面の断面曲線を正確に評価することができない。
【００５３】
表１は、粗面Ａ〜Ｄ及び基準面について取得された測定断面曲線について、パラメーター解析を行い算出された表面粗度パラメーターである。この表面粗度パラメーターには、粗度計測装置１０固有のノイズを含むため、例えば、粗度高さに関するパラメーター（ＲｚやＲｃなど）を過大評価し、粗度長さに関するパラメーター（ＲＳｍなど）を過小評価する恐れがある。
【００５４】
【表１】
2020134356.tif 000003
【００５５】
そこで、各粗面Ａ〜Ｄの粗度スペクトルから、基準面の粗度スペクトルを減算し、各粗面Ａ〜Ｄの粗度スペクトルから粗度計測装置１０固有のノイズを除去した。
図４は、粗面Ａ〜Ｄのノイズ除去後の粗度スペクトルを示すグラフである。なお、図４でも、図３と同様に、粗度スペクトルのピークを詳細に確認するため、縦軸をＨ／λとしている。
【００５６】
図４に示すように、図３では短波長側に現れていた粗度計測装置１０固有のノイズが除去され、各粗面の特徴のみを取得することが可能となっている。
【００５７】
このようにして得られた粗面Ａ〜Ｄのノイズ除去後の粗度スペクトルについて、逆フーリエ変換を行い、粗面Ａ〜Ｄのノイズ除去後の断面曲線を得た。この断面曲線について、パラメーター解析を行い、表面粗度パラメーターを算出した。その結果を表２に示す。
【００５８】
【表２】
2020134356.tif 000004
【００５９】
［比較例］
実施例１と同一の条件で、粗面Ａ〜Ｄの測定断面曲線を取得した。
このように取得した測定断面曲線に対して、ノイズ除去の方法としてＩＳＯ４２８７：１９９７で規定されるように、カットオフ値λｓ（＝１０００μｍ）以下の短波長ガウシアンフィルタを適用した。
【００６０】
短波長ガウシアンフィルタ適用後の測定断面曲線をフーリエ変換し、粗面Ａ〜Ｄの粗度スペクトルを取得した。図５は、粗面Ａ〜Ｄの短波長ガウシアンフィルタ適用後の粗度スペクトルを示すグラフである。
【００６１】
カットオフ値λｓの短波長ガウシアンフィルタを適用した結果、短波長側に現れる粗度計測装置１０固有のノイズが取り除かれているものの、粗面特有のピークの一部も同時に除去してしまっている。
【００６２】
表３は、短波長ガウシアンフィルタ適用後の測定断面曲線について、パラメーター解析を行い算出された表面粗度パラメーターである。このように短波長ガウシアンフィルタを適用することによって、粗度計測装置１０固有のノイズとともに、粗度の要素も同時に除去されてしまっているため、粗度高さに関するパラメーターを過小評価し、粗度長さに関するパラメーターを過大評価してしまっていることが明らかである。
【００６３】
【表３】
2020134356.tif 000005
【００６４】
以上、本発明の好ましい実施形態を説明したが、本発明はこれに限定されることはなく、本発明の目的を逸脱しない範囲で種々の変更が可能である。
【００６５】
また、上述する粗度解析プログラムを記録した、例えば、磁気テープ（デジタルデータストレージ（ＤＤＳ）など）、磁気ディスク（ハードディスクドライブ（ＨＤＤ）、フレキシブルディスク（ＦＤ）など）、光ディスク（コンパクトディスク（ＣＤ）、デジタルバーサタイルディスク（ＤＶＤ）、ブルーレイディスク（ＢＤ）など）、光磁気ディスク（ＭＯ）、フラッシュメモリ（ＳＳＤ（Solid State Drive）、メモリーカード、ＵＳＢメモリなど）などのコンピュータ可読記録媒体も本発明の一態様として含まれる。
【符号の説明】
【００６６】
１０      粗度計測装置
１２      二次元レーザー変位計
１４      移動機構
１６      移動距離計測装置
１８      演算装置
１９      データ入力手段
２０      可搬型バッテリー
２２      測定対象
【図１】
2020134356.tif 000006
【図２】
2020134356.tif 000007
【図３】
2020134356.tif 000008
【図４】
2020134356.tif 000009
【図５】
2020134356.tif 000010