QGIS上で衛星画像への「ピン刺し」と「描き込み」③ 直接、イラスト的に図描する方法

you1126 @Rb3qzWR5AK5nooA

Sentinel-1，2やLandsat-8，9で見たウクライナ情勢について呟いていきたいと思います。

you1126 @Rb3qzWR5AK5nooA

ネタがなくなりましたね。そして今日は欲しい画像の撮影日でもありません。という訳で、QGISを使った塹壕マッピングの方法、もう一回紹介しちゃいますか。こういうのは何回やっても良いですもんね。

2023-01-07 15:33:56

you1126 @Rb3qzWR5AK5nooA

すみません。ちょっと準備に時間をいただきます。なるべく良い画像を揃えてやりたいので。あと、方法も色々と調べます。少しでも簡単な方が良いので。

2023-01-07 17:54:57

you1126 @Rb3qzWR5AK5nooA

よし、決めた。メリトポリ周辺を舞台にやります。Landsatが2枚、Sentinel-2が4枚の計6枚です。使用するLandsatの画像IDは以下のとおりです。 LC09_L1TP_177027_20230102_20230102_02_T1 LC09_L1TP_177028_20230102_20230102_02_T1 pic.twitter.com/oOiABJBJVF

2023-01-07 21:36:33

拡大

you1126 @Rb3qzWR5AK5nooA

Sentinel-2は以下の4枚を使います。 S2A_MSIL1C_20230103T083341_N0509_R021_T36TXS_20230103T092252 S2A_MSIL1C_20230103T083341_N0509_R021_T36TXT_20230103T092252 S2A_MSIL1C_20230103T083341_N0509_R021_T36TYS_20230103T092252 S2A_MSIL1C_20230103T083341_N0509_R021_T36TYT_20230103T092252

2023-01-07 21:37:18

you1126 @Rb3qzWR5AK5nooA

遅くなりました。では始めます。今日は行けるところまで。まずはQGISを起動。 pic.twitter.com/IC8ePe72LR

2023-01-07 22:23:22

拡大

you1126 @Rb3qzWR5AK5nooA

Landsat、Sentinel-2は既にEISEIにてGeotiff形式でカラー化されているものとして解説します。開き方は簡単。ドラッグ＆ドロップでQGISに放り込みます。 pic.twitter.com/6wl5Ek39gX

2023-01-07 22:31:16

拡大

you1126 @Rb3qzWR5AK5nooA

開きました。黒いエッジ部分が邪魔なので消去しましょう。また、Landsatは解像度15m（パンシャープン処理済の状態）、Sentinel-2は解像度10mなので、解像度の高いSentinel-2を優先して使いましょう。方法を説明していきます。 pic.twitter.com/VEQUHdrfsH

2023-01-07 22:35:31

拡大

you1126 @Rb3qzWR5AK5nooA

まずはエッジを消去しましょう。画像の名称を右クリック。プロパティを選択します。 pic.twitter.com/koEYa6Ll3C

2023-01-07 22:41:07

拡大

you1126 @Rb3qzWR5AK5nooA

リストから「透明度」を選択し、「追加のnodata値」に0を入力。OKを押します。全ての画像にこの設定を適用しましょう pic.twitter.com/cC81VqyDh6

2023-01-07 22:48:48

拡大

you1126 @Rb3qzWR5AK5nooA

この設定を全ての画像に適用すると、こんな感じにエッジが消えてくれます。……なんかクリミア方面を撮ったLandsatが消えるべきでない部分も消えてるけど、まぁいいや。当然、エッジを消さずにやる方法もあります。そこはお好みで。次は、Sentinel-2を優先して使えるよう、レイヤの順番を変えます。 pic.twitter.com/zIJAmkCzeB

2023-01-07 22:56:06

拡大

you1126 @Rb3qzWR5AK5nooA

データは「レイヤ」のリストの上から順で表示されます。ですので、2枚のLandsatをリストの一番後ろに持って行きましょう。方法は簡単です。リストの一番上に来ているLandsatを左クリック長押しでドラッグし、 pic.twitter.com/dmyD1CtBpE

2023-01-07 23:07:09

拡大

you1126 @Rb3qzWR5AK5nooA

リストの一番下でドロップ。最背面で表示されるようになりました。同様にして、もう1枚のLandsatも最背面に持って行きましょう。 pic.twitter.com/fjVSHSyu8U

2023-01-07 23:09:36

拡大

you1126 @Rb3qzWR5AK5nooA

こんな感じになる筈です。Sentinel-2が前面に表示され、優先して使えるようになりました。これにて準備は完了です。ラインデータ作成機能を使って塹壕等をマッピングしていきましょう。その方法を解説して行きます。 pic.twitter.com/ZdkILO7H30

2023-01-07 23:14:27

拡大

you1126 @Rb3qzWR5AK5nooA

マッピングは緯度経度座標の状態でやるのが無難です。左下の赤四角で囲われた部分をクリック。 pic.twitter.com/o12b2cKNjj

2023-01-07 23:25:27

拡大

you1126 @Rb3qzWR5AK5nooA

こんな画面が出てきます。「座標参照系（CRS）」をクリックし、フィルタに4326と入力。WGS84座標系を選択してOKを押します。 pic.twitter.com/hELhDLFWhS

2023-01-07 23:43:14

拡大

you1126 @Rb3qzWR5AK5nooA

赤枠内の「座標」に出てくる数字のとおり、表記が緯度経度に変わりました。 pic.twitter.com/BnCcaGFN5O

2023-01-07 23:46:49

拡大

you1126 @Rb3qzWR5AK5nooA

すみません。続きは明日にします。明日の夜まで良い画像は入らなさそうなので、のんびりやらせてもらいます。ご了承ください。

2023-01-07 23:56:12

you1126 @Rb3qzWR5AK5nooA

再開します。せっかくなので、メリトポリの東にある例の長い防御線をマッピングしてみましょう。46.819902806,36.251189355の辺りからスタートします。画面下に「座標」とありますが、ここに座標値をペーストしてEnterを押せば、その位置に移動してくれます。 pic.twitter.com/54jufdPzDZ

2023-01-08 09:01:23

拡大

you1126 @Rb3qzWR5AK5nooA

なぜか「°」が入った状態ではEnterが反応しないので注意。表示では「°」が入っているんですけどね……。移動したら、これも同じく画面下の「縮尺」を調整してズームします。あんまりズームしすぎても見づらいので、1:10000スケールくらいで作業するのが良いでしょう。

2023-01-08 09:04:50

you1126 @Rb3qzWR5AK5nooA

ラインデータの作成機能を使ってマッピングしていきます。ここをクリック。こんな感じの画面が立ち上がるので、パラメータを設定していきます。 pic.twitter.com/VIUEzzQipv

2023-01-08 09:12:27

拡大

you1126 @Rb3qzWR5AK5nooA

基本設定の部分はこのように設定しましょう。これから作成するラインデータはshp形式ですが、複数のデータで構成されるため、きちんとフォルダを作って格納しておきましょう。文字コードはUTF-8、ジオメトリ型はラインストリング、座標系はEPSG:4326 WGS84に設定します。 pic.twitter.com/uub1ZlfAXb

2023-01-08 09:28:28