【R】simplevizでGIS

こちらでsimplevizパッケージを見てみましたが、それをGISへ応用してみます。

こちらのサイトをなぞってみました。

データの読み込み

library(simplevis)
library(tidyverse)

example_sf_point

> example_sf_point
Simple feature collection with 112 features and 3 fields
Geometry type: POINT
Dimension:     XY
Bounding box:  xmin: 1175354 ymin: 4853914 xmax: 2025939 ymax: 6096100
Projected CRS: NZGD2000 / New Zealand Transverse Mercator 2000
First 10 features:
     site_id median trend_category                geometry
1  ARC-00001 0.0140      Improving POINT (1735609 5916179)
2  ARC-00008 0.0610      Improving POINT (1753479 5976281)
3  ARC-00013 0.1310      Improving POINT (1742066 5915382)
4  ARC-00014 0.9900      Improving POINT (1764285 5907017)
5  ARC-00015 1.0300      Improving POINT (1767401 5907336)
6  ARC-00016 0.2980      Improving POINT (1768314 5908177)
7  ARC-00017 0.3550      Improving POINT (1751305 5933319)
8  ARC-00018 0.7350  Indeterminate POINT (1769952 5912814)
9  ARC-00019 0.5000      Improving POINT (1769452 5910614)
10 ARC-00026 0.1295      Improving POINT (1748608 5953465)

ggplotのwrapperとして使ってみます。

gg_sf_col(example_sf_point,  
          col_var = trend_category, 
          borders = nz)

leafletでも使えます。インタラクティブなマップを利用できます。

leaflet_sf_col(example_sf_point, 
               col_var = median, 
               col_method = "quantile",
               col_cuts = seq(0, 1, 0.25),
               title = "Monitored medians, 2008-17") %>% 
  leaflet::addPolygons(data = nz, 
                       color = "green", 
                       weight = 2, 
                       fillOpacity = 0, 
                       opacity = 1)