解凍が終わるとxmlファイルが現れます。 数値標高モデルの変換 まずは、標高データをGeoTiffに変換します。 すでに国土地理院の数値標高モデルデータをラスタとしてRで扱うで紹介されている手法です。 あらかじめデータを格納するmapディレクトリを作成し、解凍したxmlファイルの入ったディレクトリをその直下に置いています。 #1 ダウンロードして解答したxmlファイルのリストを作成します。 xml <- list.files(path="map/PackDLMap_DEM/", pattern=".*\\.xml", full.names=T) #2 最初のラスタを読み込む dem <- read_fgd_dem(xml[1], resolution = 10 , return_class = "raster") #3 for文でxmlファイルをラスタオブジェクトに変換して結合 for(i
![国土地理院基盤地図情報をRを使ってGISデータに変換する - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/41900dce3eb9456d86cde9ec786439efa3af016b/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-412672c5f0600ab9a64263b751f1bc81.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTk3MiZoPTM3OCZ0eHQ9JUU1JTlCJUJEJUU1JTlDJTlGJUU1JTlDJUIwJUU3JTkwJTg2JUU5JTk5JUEyJUU1JTlGJUJBJUU3JTlCJUE0JUU1JTlDJUIwJUU1JTlCJUIzJUU2JTgzJTg1JUU1JUEwJUIxJUUzJTgyJTkyUiVFMyU4MiU5MiVFNCVCRCVCRiVFMyU4MSVBMyVFMyU4MSVBNkdJUyVFMyU4MyU4NyVFMyU4MyVCQyVFMyU4MiVCRiVFMyU4MSVBQiVFNSVBNCU4OSVFNiU4RiU5QiVFMyU4MSU5OSVFMyU4MiU4QiZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTU2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9MWFmMTY2YjFhMGRiZGI4ZWMyY2Y3Njg4ZWRlZTNhMTc%26mark-x%3D142%26mark-y%3D57%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZoPTc2Jnc9NzcwJnR4dD0lNDBpc2hpaWp1bnBlaSZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTM2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9YThkOTRiODQ2NDBkN2QwZTExNjk1ODA3ZmE5ZGFmMzc%26blend-x%3D142%26blend-y%3D486%26blend-mode%3Dnormal%26s%3Dbc49427d7e4390f9d00febaaa77cc4c2)