"所定の物理量が実質的に所定の間隔毎に離散した等間隔信号及び前記所定の間隔から外れるよう離散した不等間隔信号が混在する信号群から前記等間隔信号を抽出する等間隔信号の抽出方法であって、前記信号群に含まれる各信号に対応する前記所定の物理量を前記所定の間隔で除した値を整数部と小数部に分けて、整数部の値及び小数部の値を各軸とした直交座標系にプロットし、前記直交座標系において近似直線推定処理によって求めた直線上及び前記直線に近接してプロットされた信号を前記等間隔信号であるとして抽出することを特徴とする等間隔信号の抽出方法。前記近似直線推定処理がハフ変換であることを特徴とする請求項1に記載の等間隔信号の抽出方法。前記等間隔信号は鉄道路線に沿って等距離間隔に配置された目標物を鉄道車両が通過時に生成されるものであり、前記所定の物理量は、車両の走行速度検出値に基づいて求められる鉄道車両の推定走行距離であることを特徴とする請求項1又は請求項2に記載の等間隔信号の抽出方法。しかし、この方法は、ノイズなどの除去すべきパルスが、抽出したいパルスと同じ間隔で発生した場合に、そのパルスが抽出されてしまうといった問題がある。また、一般に測定におけるサンプリングは一定時間間隔であるのに対し、上述のような方法をキロポストマーカーの抽出に適用するには、パルス信号を一定距離間隔でサンプリングする必要がある。上述した問題に鑑み、本発明の課題は、実質的に所定の間隔に離散した等間隔信号及び前記所定の間隔から外れるよう離散した不等間隔信号が混在する信号群から等間隔信号を適切に抽出する等間隔信号の抽出方法を提供することである。上述した課題を解決するため、本発明の等間隔信号の抽出方法は、所定の物理量が実質的に所定の間隔毎に離散した等間隔信号及び前記所定の間隔から外れるよう離散した不等間隔信号が混在する信号群から前記等間隔信号を抽出する等間隔信号の抽出方法であって、前記信号群に含まれる各信号に対応する前記所定の物理量を前記所定の間隔で除した値を整数部と小数部に分けて、整数部の値及び小数部の値を各軸とした直交座標系にプロットし、前記直交座標系において近似直線推定処理によって求めた直線上及び前記直線に近接してプロットされた信号を前記等間隔信号であるとして抽出することを特徴とする。本発明によれば、パルス間隔に誤差が生じた場合であっても、傾きやばらつきは生ずるものの各パルスのプロット位置は直線的に配列されることになるため、この直線を検出してその近傍の点を抽出することによって、適切に等間隔信号を抽出することができる。本発明において、前記近似直線推定処理がハフHough変換である構成とすることができる。ハフ変換は、直線の近傍の点が最も多くなる直線を推定する直線検出方法であり、これによれば、等間隔信号以外の外れ点の影響を低減しつつ直線的に並んだ等間隔信号のプロット位置近傍を通る直線を推定することができ、このような直線の近傍の点を抽出することによって、より適切に等間隔信号を抽出することができる。本発明において、前記等間隔信号は鉄道路線に沿って等距離間隔に配置された目標物を鉄道車両が通過時に生成されるものであり、前記所定の物理量は、車両の走行速度検出値に基づいて求められる鉄道車両の推定走行距離である構成とすることができる。これによれば、走行速度検出値に基づいて求められる鉄道車両の推定走行距離を適切に補正し、鉄道車両の走行位置を正確に検出することが可能となる。以上のように、本発明によれば、実質的に所定の間隔に離散した等間隔信号及び前記所定の間隔から外れるよう離散した不等間隔信号が混在する信号群から等間隔信号を適切に検出することができる。以下、本発明を適用した等間隔信号の抽出方法の実施形態について、図面を参照して詳細に説明する。実施形態の等間隔信号の抽出方法は、鉄道車両の走行試験において、例えば1km間隔に設置されたキロポスト通過時、及び、その他の構造物等通過時に生成されるパルス信号utから、キロポスト通過時の等距離間隔パルスを抽出するものである。先ず、実施形態の等間隔信号の抽出方法を実行するのに先立って、鉄道車両を実際に走行させながら、上述したパルス信号を記録する。このとき、鉄道車両の走行速度履歴も併せて記録する。鉄道車両の走行速度検出手法として、例えば輪軸の回転速度に応じた速度信号vtを出力する車速発電機を用いることができる。次に、速度信号のオフセット誤差を除去する。具体的には速度信号vtがあらかじめ設定した閾値以下の場合には、鉄道車両が停止中であるものとみなして、その停止中の速度データのメディアン中間値を、速度データ全体から差し引いて、これを速度信号vtとして再定義する。なお、速度信号の記録前に予めオフセット誤差を除去していれば、この手順は不要である。次に、上述した速度信号vtを積分することにより、出発記録開始地点からの累積距離xtを算出する。本実施形態の等間隔信号の抽出方法は、この累積距離xtを用いて、上述したパルス信号utの中から等距離間隔パルスを抽出するものである。この抽出手法については、後に詳しく説明する。抽出したパルスが立ち上がる時の累積距離をxi=x1x2xnとする。次に、抽出した各パルス間の所定の距離と、累積距離との比から、補正係数αiを求める。パルス間の所定の距離は、パルス間隔λの整数倍となるように丸めることで算出する。このとき、最初のパルスx1より前、及び、最後のパルスxnより後の補正係数は、各パルス間の補正係数のメディアン値を用いる。ここで、平均値を使うことも可能であるが、特定のパルス間において、比較的大きな空転や滑走が発生した場合の影響を受けないようにするためには、上述したようにメディアン値を用いることが好ましい。αii=1,2,,n-1は、以下の式1により表される。また、α0,αnは、以下の式2により表される。α0=αn=medianαii=1,2,,n-1式2そして、式3、式4のように、各パルス間ごとの距離に、それぞれ上述した補正係数αiを乗じ、それを補正累積距離x’tとする。x’i1=x’iαixi1-xi,i=1,2,,n-1式3次に、補正累積距離x’tを、当該路線の起点からの位置距離程dtに変換する。先ず、検出したいずれかの等距離間隔パルスの距離程を求める。例えば、出発地点の距離程dstartから、最初のパルスの距離程d1を、以下の式5を用いて推定する。等間隔距離パルスの距離程がλの整数倍であるとして、推定距離程を丸めることによって、dstartの誤差が許容できる。ここで、γは進行方向を示しており、起点から遠ざかるときはγ=1、起点に近づくときはγ=-1である。d1=λround[dstartγx’1/λ]式5そして、以下の式6に示すように、上述したd1を用いて補正累積距離x’tをシフトし、全体の距離程dtを求める。その際に、進行方向γをかける。dt=γ[x’t-x’1]d1式6以下、本実施形態の等間隔信号の抽出方法を用いた計算例について説明する。図1は、鉄道車両の走行試験時の速度信号vtの時間履歴の一例を示すグラフである。図2は、この走行試験時に入力されるマーカーパルス信号utの履歴を示すグラフである。このマーカーパルス信号には、距離程1km毎に等間隔に設置されたキロポスト通過時に入力される信号等間隔信号、及び、その他の目標通過時に入力される信号不等間隔信号が含まれている。図3は、図1に示す速度信号を積分して求められる累積距離xtを示すグラフである。ただし、この累積距離には、車速発電機が設けられた輪軸の車輪の磨耗や、空転、滑走等の影響により、実際の走行距離に対して誤差が生じている場合がある。図4は、図2のマーカーパルス信号から抽出されたキロポストマーカパルスを示すグラフである。このキロポストマーカパルスの抽出方法については、後に詳しく説明する。図5は、図4に示す抽出されたキロポストマーカパルスを累積距離上にプロットしたグラフである。図5に示すように、累積距離は誤差を含むため、抽出されたキロポストマーカパルスの間隔は正確に1kmになっていない。図6は、上述した式1乃至式6を用いて算出された距離程dtの一例を示すグラフである。図6に示すように、抽出されたキロポストマーカパルスは、距離程1km毎に端数なく配列されていることがわかる。以下、本実施形態における等間隔信号等間隔パルスの抽出方法について、詳細に説明する。なお、この等間隔信号の抽出方法は、鉄道車両の走行位置検出に限らず、他の用途にも適用可能なものである。ここでは、一般化した記述によって方法を説明し、後にキロポストマーカパルス検出へ適用した応用例を示す。先ず、検出したいパルスが等間隔となる物理量をxとする。ここで、xは、例えば距離や時間である。閾値処理等により、パルス信号uxの中のパルスを検出し、各パルスの立ち上がり時のxを並べた数列xk={x1,x2,,xm}を用意する。検出したいパルス間隔をλとして、式7のようにxkをλで除して正規化し、それをXk={X1,X2,,Xm}とする。Xk=xk/λ式7そして、そのXkを、式8のように整数部pkと小数部qkとに分ける。Xk=pkqk式8仮にuxが完全なλ間隔のパルス信号であれば、Xkは1ずつ増える等差数列になるので、その小数部qkは一定の値をとることになる。これを、横軸p、縦軸qのp-q平面の1点としてプロットすると、横軸pと平行な直線となる。もし、これにλ間隔以外のパルスが加われば、このパルスはその直線から離れた点として表される。そこで、p-q平面で直線を構成する点を抽出することによって、等間隔パルスを抽出することが可能となる。実際は、パルス間隔に少なからず誤差があるので、その直線は傾きやばらつきを持つことが想定される。その誤差に対応するため、本実施形態においては、近似直線を求めてそれの近傍の点を抽出している。ここで、近似直線を求める一般的な方法として、最小二乗法がある。最小二乗法による近似直線の推定は、直線と各点との差残差の二乗和が最小となるような直線を算出するものである。全ての点が等間隔パルスである場合には、全ての点が直線的に並ぶため、最小二乗法はそれらの点の近傍を通る直線を推定するが、等間隔パルス以外の外れ点が含まれる場合、最小二乗法ではその外れ点を含めた残差の二乗和を最小にするため、直線的に並んだ等間隔パルスから外れた近似直線が推定される問題がある。これを解決するために、本実施形態においては、画像から直線を検出する手法として用いられるハフHough変換を用いて近似直線を算出している。ハフ変換は、直線の近傍の点が最も多くなる直線を推定するもので、一般にノイズや直線の途切れに強い特徴がある。等間隔パルス以外の外れ点が多く含まれていても、それらは直線的に並ばないので、推定される直線は、最小二乗法と異なり外れ点の影響を受けにくく、直線的に並んだ等間隔パルス近傍を通る直線を推定できる。その直線の直近の点を抽出することで、適切に等間隔パルスを抽出できる。具体的には、求めた近似直線に近い点を選び出すために、まず近似直線上に乗る模範数列Xjを用意する。近似直線方程式は、横軸が整数でなければならないので、式9のように表すことができる。qj=AjB,j=,-1,0,1,2,式9式8のように、横軸と縦軸の加算がもとの数列であるので、模範数列は式10のように表すことができる。Xj=jqj式10式10に式9を代入すると、模範数列は式11で求めることができる。Xj=A1jB式11そして、模範数列Xj={,X0,X1,}の各要素に対して直近のXkを抽出していく。ただし、Xjと直近のXkとの差が予め設定した許容範囲の外であればそれは抽出しない。抽出したものをxi={x1,x2,,xn}とする。最後に、抽出したXiに、式12の通りパルス間隔λを乗じてもとの物理量に戻す。xi=λXi式12以下、計算例として、パルス間隔λ=1[km]間隔のキロポストマーカを抽出した結果を示す。図7は、p-q平面上にプロットされたXk、ハフ変換により求めた近似直線、及び、抽出したパルスXiの一例を示すグラフである。なお、参考のため、最小二乗法によって求めた近似直線も記載している。図8は、図7の近似直線を用いて抽出されたキロポストマーカパルスを示すグラフである。図7に示すように、最小二乗法では、キロポストマーカ以外の外れ点を含めて残差の二乗和を最小化しようとするため、等間隔のキロポストマーカの並びに対して、得られる直線のずれが大きいことがわかる。これに対し、ハフ変換においては、キロポストマーカの直線的な並びに沿って近似直線を算出することが可能であり、図8に示すように、適切に等間隔のキロポストマーカを抽出し、鉄道車両の走行位置を正確に検出することが可能となる。他の実施形態なお、本発明は上述した実施形態のみに限定されるものではなく、種々の応用や変形が考えられる。上述した実施形態では、速度信号を積分して累積距離を算出しているが、これに限らず、例えば車軸の回転に応じて速度発電機が出力する回転パルスの発生数をカウントして累積距離を算出してもよい。また、本発明の等間隔信号の抽出方法は、実施形態のような鉄道車両の走行位置検出に限らず、他の用途にも適用することが可能である。例えば、高速道路等の道路の距離標を用いて自動車等の位置検出を行なうこともできる。また、等間隔で離散する物理量は走行距離に限らず、例えば時間やその他の物理量であってもよい。"