須崎純一の研究紹介

測量学を基礎に、衛星リモートセンシングや写真測量、レーザ計測等の計測、画像処理を通じて、都市を対象にした空間データや地図の生成に関する研究を行っています。

[業績表（簡易版）] [業績表（大学のWebページ）]

衛星リモートセンシング

三次元変動地盤沈下・土木インフラストラクチャの変状モニタリング

衛星搭載型のマイクロ波レーダである合成開口レーダ(Synthetic Aperture Radar: SAR)で取得された複数時期のデータから、地物の変動を推定できます。レーダでは、散乱体（主に地上の地物）から反射された信号を受信します。振幅と位相が反映された信号から複数時期のデータ間の位相差を計算し、位相差から変動量に換算します。二時期のデータからは、ある基準面からの相対的な高さが求められます（干渉SAR: Interferometric SAR または SAR interferometry (InSAR)）。また二時期のデータから推定された標高データと、別途用意された標高データとの差を調べることで、二時期間に生じた地盤変動を推定することも可能です（差分干渉SAR: Differential InSAR (DInSAR)）。

更に、多時期の画像を用いて変動速度を推定する時系列SAR（Time-series SAR）解析も広く普及しています。代表的な手法として、Permanent Scatterers Interfrometry (PSI) が挙げられます。干渉位相に含まれる大気による遅延や軌道の誤差による位相を、時間的、空間的な平滑フィルタリングを行うことで取り除き、変動速度と、標高データ (Digital Elevation Model: DEM) と散乱体とのずれを表すDEM誤差を同時に推定する技術です。

しかし、差分干渉SARや時系列SAR解析によって観測される変動は、衛星視線方向のみの一次元変動です。我々の研究グループでは、まず衛星が南から北へ向かう上昇軌道と、北から南へ向かう下降軌道という異なる軌道で取得されたSAR画像から、各々の変動ベクトルを求めます。それに加えて、例えば、国土地理院が運用しています、10 km〜20 km間隔で設置されている電子基準点のデータといった、Global Navigation Satellite System (GNSS) の三次元変動測量結果も使います。GNSSのデータは東西、南北、鉛直の三方向の変動に分解し、クリッギング (Kriging)等の空間統計学の手法を用いて内挿します。それらの観測データに対し、衛星観測時の天頂角や相対方位角等を考慮して最小二乗法で解くことで、地上での三次元変動を推定しています。

関西国際空港の東西方向（左）、南北方向（中）、鉛直方向（右）の変動成分。ALOS-2/PALSAR-2のレンジ・アジマス分解能3 mのUBSモードで、上昇軌道データ13シーン、下降軌道データ17シーンを用いた。

抜け上がり

地盤沈下が発生している地域では相対的に構造物が地盤面から浮き出てくる、抜け上がりと呼ばれる現象が発生します。抜け上がりが発生すると、地下で構造物と接続する埋設パイプの破断や構造物の支持杭の破損を引き起こします。抜け上がりが発生している地域を特定し、抜け上がり量を定量的に把握することは土木構造物の維持管理の観点からも重要ですが、多大な時間と労力を消費するため効率的な監視手法が必要とされています。

本研究では、時系列SAR解析手法の一種であるPermanent Scatterer Interferometry (PSI) を改良して、地盤面恒久散乱体 (permanent scatterer: PS) と構造物PSの変動速度の差を利用して抜け上がり速度を推定する手法を開発しました。

タイ王国バンコク市周辺の抜け上がり発生地域・非発生地域に対して、提案手法を適用しました。現地計測結果と比較した結果、RMSE=約6 mm/yearと十分な精度で推定できていることが確認できた。さらに、本手法をSAR画像のシーン全体に対しても適用した結果、バンコク中心部から郊外部へと抜け上がり速度が大きくなっている様子が確認され、郊外部において抜け上がりの影響が深刻である様子が示唆されました。

バンコク空港周辺の抜け上がり推定結果： （左上）恒久散乱体(PS)の分類結果。sPSは構造体由来のPS、gPSは地盤面由来のPS、（右上）地盤面由来のPSの変動速度、（左下）構造体由来のPSの変動速度、（右下）抜け上がり速度。20枚のTerraSAR-X画像を用いた。

都市密度推定

衛星搭載型の多偏波合成開口レーダ(Polarimetric Synthetic Aperture Radar: PolSAR)を用いて、建蔽率や容積率のような建物の密度を反映した指標を作成し、都市密度を推定する研究を行っています。通常、衛星画像や航空写真から3次元データを推定するには、同一地点を異なる視点から撮影した2枚の画像が必要になりますが、この手法では一時期の多偏波データがあれば高さを反映した都市密度を推定できます。そのため、統計データが不足する地域を含めた都市構造の比較に役立てられる可能性があります。

人工衛星を用いて推定した都市密度分布図： （左）JAXAが打ち上げたALOS衛星に搭載されたPALSARセンサの4偏波データを利用して生成した都市密度分布図、（中）ZENRINのZmap TownIIを用いて生成した建蔽率、（右）同じく容積率。生成された都市密度分布図は建蔽率、容積率の両方の中間に相当する性質を持ち合わせている。

写真測量・画像処理

建設現場におけるリアルタイム三次元地図作成

現在、建設現場では省力化、効率化を目指して、Information and Communication Technology (ICT) を活用した技術開発が進められています。我々は、民間会社との共同研究を通じて、建設機械に取り付けられたビデオで撮影された動画像から、リアルタイムに高精度な三次元地図を生成、更新することを目標に研究を進めています。

基礎的な段階ではありますが、建設機械から撮影された動画像に加えて、大学内でも基礎実験を繰り返しています。ORB SLAM2という既存のアルゴリズムを用いて、一定の精度で周辺環境の三次元地図を生成できることが分かりました。

地上撮影した動画像からORB SLAM2を用いて生成された三次元点群。青いフレームはカメラの位置と向きを、黒い点は三次元座標が確定して点群、赤い点は暫定的に生成された（更新の可能性のある）点群を表す。

桂キャンパス屋外での実験現場（左：上記の赤線、右：上記の青線に対応する視線）

3次元点群処理（（航空機・地上）レーザ・航空写真）

フィルタリング、地盤高推定

航空機LiDAR計測により得られるデータは地表面と建物や樹木等の表面を含んだ標高データ (Digital Surface Model: DSM) であり、建物等の下面の地盤面データは取得できません。森林の樹高推定や建物モデル作成のためには、DSMから地盤面データが取得されていない部分を補間した地盤面だけの標高データ (Digital Terrain Model: DTM) を推定する必要があります。この処理は一般的にフィルタリングと呼ばれています。

我々の研究グループでは、これまでは困難であった傾斜が平坦、急峻な地形が混在した地域においても、地盤面データの抽出精度、および地盤面データが存在しない地点の地盤面高さの推定精度を向上するアルゴリズムを開発しました。従来の手法とは異なり、本手法では局所的に異なる最大許容傾斜値を自動的に推定するため、精確なDTM生成アルゴリズムとして有用性が高いといえます。

フィルタリングの対象地域（左）、航空機LiDARデータ（中）、推定された地盤高 (DTM)（右）

植生の自動抽出

地上LiDARで得られた3次元点群データから、後述する緑視率を推定するために、点群を植生と非植生に分類する手法を開発しました。ボクセルを分類の最小単位として、ボクセル内の点の散布状況を利用して、主成分分析から植生の形状を反映する特徴量を計算し、分類しました。その上で、多段階の空間スケールを用いて処理を組み合わせることで点群分類精度の向上を図りました。

地上LiDARデータから自動抽出した植生。緑はTrue Positive（植生を正しく抽出した結果）、青はTrue Negative（非植生を正しく抽出した結果）、赤はFalse Positive（非植生を誤って植生として抽出した結果）、白はFalse Negative（植生を過って抽出できなかった結果）を表す。

景観指標推定

航空機LiDARデータや地上LiDARデータ、または車載搭載型LiDAR (Mobile Mapping System: MMS) データを用いて、建物や植生等の一切の地物による囲まれ感を表す囲繞度（いじょうど）を推定することができます。視野内に含まれる植生が抽出できれば、その割合を表す緑視率を計算できます。我々の研究グループでは、航空写真から植生を抽出したり、または地上LiDARやMMSデータから直接植生の分布を把握したりして、緑視率を推定する手法を確立しました。

航空機LiDARデータと航空写真を用いて推定した京都市での囲繞度（左）と緑視率（中）の推定結果。（右）MMSデータを用いて推定した京都市東山区での緑視率。

交通流画像・データ処理

車両軌跡抽出

交通流を撮影した動画像から車両の位置や軌跡、速度を自動で抽出・推定する手法が長年研究されてきました。昼間の動画像において太陽光による建物や車両の影が発生する場合、画像空間上で車両領域を抽出しようとすると影が混入してしまいますが、このような影に頑健な手法は存在しません。また昼間だけでなく夜間も運用可能な手法も限られます。そこで私達は、日照変化に頑健で車両や建物の影の影響を受けにくく、昼夜間を問わず適用可能な交通流からの四輪車の位置と軌跡の自動抽出手法を提案しました。

動画像から抽出された車両

動画像から抽出された車両の軌跡

道路ネットワーク自動更新

デジタル道路地図 (Digital Road Map: DRM) データベースを構成するノードの位置ずれを自動補正・更新することを目的とし、プローブデータを用いたノード座標更新手法を開発しました。更新する際に、ノード単位で更新する手法、閉ネットワーク単位で更新する手法、更に一括して更新する手法を構築し、最小二乗計算の計算単位の違いが更新結果に与える影響を比較・検証しました。位置情報が多く分布するリンクのみを用いて計算した場合には、3つの手法の結果には有意な差は見られませんでした。一方、位置情報の分布が少ないリンクも用いて更新計算した場合には、一括更新モデルでの更新結果が最も誤差が小さくなり、比較的良好な更新モデルであることが確認できました。

デジタル道路地図 (DRM) の対象地域（左）と自動更新結果の一例（右）