はじめに
いまは公開停止となっているiOSアプリfeildaceess風に、写真を上部に並べて、写真と写真撮影位置と線で結んで表示するQGISスタイルとモデルを作ったので紹介します。
feildaccess2（iOSアプリ）のスクリーンショット（ちなみに場所は愛媛県の滑床渓谷です）
QGISスタイル、モデルファイルの置き場所
使い方の流れ
ジオタグ付き写真の写真を用意する
GPS付カメラで写真を撮るなどします
モデルの読み込み
QGISを起動し、プロセッシングツールボックスの歯車のアイコンから「ツールボックスにモデル」を追加を選択し、上記に置いている2つの.model3ファイルを指定します。
そのあと、PhotoView_Inport を右クリック→「出力のためのレンダリングスタイルを編集」→「PhotoView.qml」を設定します。
プロジェクトのCRSに投影座標系（平面直角座標）を設定
メニューのプロジェクト→プロパティ→座標参照系 を開き、写真を撮影した地域の投影座標系（平面直角座標）を設定します。
https://www.gsi.go.jp/sokuchikijun/jpc.html
残念ながら、地理座標系では動きません。また、地図を回転すると動きません。
ジオタグ付き写真からレイヤを作成
プロセッシングツールボックス→「PhotoView_Inport」を開き、Folderにジオタグ付き写真を保存したフォルダを指定し、実行します。
マップキャンバスの上部に写真が表示されます。
このモデルは、下記を行い、ポイントレイヤを出力します。
- プロセッシング「ジオタグ（位置情報）付きの写真」
- プロセッシング「レイヤの再投影」。上記をプロジェクトCRS（投影座標が期待される）に再投影
- プロセッシング「属性のリファクタリング」2つのフィールド追加
- invisible(真偽値） 真の場合、写真の表示対象から除外
- label（STRING) 値が入っている場合、ラベルとして使用。ない場合filenameフィールドを使用
表示のパラメータの設定・更新
「PhotoView_Parameter」のモデルを使ってプロジェクト変数を設定・更新すると、写真の枚数や余白などのパラメータをを変えることができます。以下は写真枚数を変更した例
プリントコンポーザ
プリントコンポーザでも同様に表示されます。
内部の処理について
座標の計算など
写真はラスター画像マーカで、写真の上中央をポイントジオメトリで指定しており、この位置の計算がこのスタイルのキモになります。
画面中央から近い画像を表示することを実現するには、他のポイントとの比較をする必要があり、表示範囲のポイントを抽出して、画面中央との距離の配列を作って並べ替えた後、指定数を抽出し、さらにX座標で並べ替えて位置決めをしています。
以下は、ラスター画像マーカのアンカー位置を計算する式です。
/*
ラスタ画像マーカーのアンカー位置、ラベルのアンカー位置 共通
*/
with_variable('Number', to_int( coalesce(@PhotoView_Number , 5) ) ,
with_variable('TopMargin', to_real(coalesce(@PhotoView_TopMargin ,8.0))*@map_scale/1000,
with_variable('PhotoBySideMargin',to_real(coalesce(@PhotoView_PhotoBySideMargin , 15.0)),
with_variable('SideMargin', @map_extent_width /(1+(@PhotoBySideMargin+1)*@Number),
with_variable('Extent',intersection( @map_extent , translate( @map_extent ,0,-@PhotoBySideMargin*@SideMargin*0.800-@TopMargin)),
with_variable('i', --左から何番目か
array_find(
array_foreach(
array_sort(
array_foreach(
array_slice(
array_sort(
array_agg( array( distance( @map_extent_center , $geometry ),$id ,$x ) ,filter:=if( attribute( 'invisible' ),false,true) and within( $geometry , @Extent ))
)
,0,@Number-1)
, array_reverse(@element))
)
,@element[1])
,$id), -- i の定義終わり
with_variable('p0',translate(@map_extent_center,(1+@PhotoBySideMargin/2)*@SideMargin-@map_extent_width/2,
@map_extent_height/2-@TopMargin),
translate(@p0,(1+@PhotoBySideMargin)*@SideMargin*@i,0)
)))))))
下図は、座標計算のメモです。
with_variable()は関数内で使えるローカル変数を定義できる関数で、これなしではこれは作れません。
構文 with_variable(name,value,expression)
引数
name 設定する変数の名前
value 設定する値
expression 変数が使用できる式
曲線の引き出し線
参考にしたfeildaceessと違い、写真の引き出し線は曲線になっています。
これは、写真の縦横比が取得できないため、写真の下端が設定できず、直線だと適切な位置に写真側の端点を設定できなかったため、写真からいったん下げて、写真と結んだラインをスムージングするという処理にしたためです。
処理速度
個別の地物ごとに表示範囲の地物のジオメトリの配列化→画面中心との距離計算→並べ替えを行っているので、表示範囲の地物の数が増えると極端に重くなります。下記は地物が500個の例で、ノートPCのi7で12秒かかりました。100個くらいなら、そんなに処理は重くありません。
おわりに
旅行で写真を撮ったら、試しに使ってみてください。