地図データの可視化に役立つcptcityのカラーパレットをためす

[cptcity](http://soliton.vm.bytemark.co.uk/pub/cpt-city/)パッケージは膨大なカラーパレットを提供しています。この中から、地図表現で利用できるカラーパレットを試しました。

# パッケージ読み込み
```{r}
library(cptcity)
```

# 遺跡分布図の背景に適したカラーパレットをさぐる

描画する地図は厚沢部町の遺跡分布図です。

```{r}
p<-
ggplot() +
  geom_raster(data = df_dem ,   # 標高ラスタ
            aes(x, y, fill = Elevation),
            hjust = 0, vjust = 0) +
  geom_contour(data = df_dem, aes(x, y, z = Elevation), col = "gray70", size = 0.2) +  # 等高線描画
  geom_sf(data = asbWl , fill = "#B4CDD8" , colour = "#B4CDD8") +  # 河川
  geom_sf(data = asb , col = "gray40" , fill = NA, lwd = 0.8) +   # 厚沢部町域
  geom_sf(data=asbSite , aes( colour = class ) , show.legend = "point") +   # 遺跡
  theme_void() +
  scale_colour_viridis(option = "cividis" , discrete =TRUE ) +
  guides(fill=FALSE) 
```

![01.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/2b4f5ea9-d6cf-387d-0c9c-a8ee1130a6bb.png)



# 地図に特化したcptパレット

[cptcity](http://soliton.vm.bytemark.co.uk/pub/cpt-city/)のHPでは膨大なカラーパレットが紹介されています。登録されているカラーパレットは7140件！！あります。

```{r}
> length(cpt_names)
[1] 7140

```

このうち、地図に特化したカラーパレットは「Palettes for topography」で紹介されています。
DEM_posterとか、派手ですね。
![02.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/966d861c-7f7c-30db-811f-273714d8c687.png)
cptcityに登録されているパレット名を調べるためには`find_cpt()`関数を使います。含まれている文字列から検索できます。DEM_posterは、`td_DEM_poster`という名前で登録されているようです。

```{r}
> find_cpt("poster")
[1] "td_DEM_poster"
```

### td_DEM_poster

予想通り派手です・・・

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "td_DEM_poster" , n = 50)) 
```
![03tdDEMposter.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/2090fbfd-1837-1cf1-4d31-5d84a3d8597c.png)

### td_DEM_screen

おだやかになりましたが、ちょっと緑がくどい。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "td_DEM_screen" , n = 50)) 
```
![04tdDEMscreen.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/cefe1e3d-4f54-27fc-2438-34e8561c7fc4.png)

### jjg_dem_c3t1

落ち着いているので、背景図としては好ましいですね。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "jjg_dem_c3t1" , n = 50)) 
```
![05c3t1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/20070760-0edf-f9a0-e515-11d99f3d2b92.png)


### jjg_dem_c3t3

こちらも落ち着いています。標高の高い部分の表現が多少違います。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "jjg_dem_c3t3" , n = 50)) 
```
![06c3t3.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/0f19a4e3-b8be-14a8-f752-9ea96384b61e.png)

### fme_metres_neutral

標高の高い部分の表現がちょっと違和感あります。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "fme_metres_neutral" ,n=50)) 
```
![07fme_neutral.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/4865acf6-fbd4-1909-33ee-56078a5edeed.png)

### jm_tv_tv_a

ちょっと安っぽい感じかなあ・・・ビビッドすぎる感じが好ましくない。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "jm_tv_tv_a" ,n=50)) 
```
![08tv-a.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/7cf954a4-f425-a3fc-b79c-2d08e48ca39d.png)

### jm_o2_o2_a

控えめで良い感じ。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "jm_o2_o2_a" ,n=50)) 
```
![09o2-a.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/fd2569bb-2f35-23ea-8dcb-73426ab27fb3.png)

### jm_ao_ao_a 

地形がわかりやすいとはいえませんが、背景図としては良いかもしれません。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "jm_ao_ao_a" ,n=50)) 
```
![11ao-a.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/8136bb26-c068-cc15-8b54-27f0dc3f325a.png)

### jm_sd_sd_a

オーソドックスな感じ。地形は見やすいです。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "jm_sd_sd_a" ,n=50)) 
```
![12sd-a.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/179dc242-9e70-667e-e42f-af2a473ddde7.png)

### esri_hypsometry_na_palm_springs_4

なんか、カラーパレットが反転しているようにも感じる。ノースアメリカというテーマ。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hypsometry_na_palm_springs_4" ,n=50)) 
```
![13palm-springs4.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/7582cea2-1221-a9ec-0558-0beeedcd59a8.png)

### esri_hypsometry_na_georgia

同じノースアメリカでもこちらはまた別の意味で特徴的なジョージアというパレット

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hypsometry_na_georgia" ,n=50)) 
```
![14georgia.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/569d4239-21f3-e3f9-c50c-55f811627ef5.png)

### esri_hypsometry_eu_europe_2

esriのeuropaというテーマ。2〜9までありますが、europa2か3が一番しっくりきます。地形がわかりやすいパレットです。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hypsometry_eu_europe_5" ,n=50)) 
```
![15europa2.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/b4fb9cef-99de-4a1f-d33b-e087ded24490.png)

### tp_tpushum

Tom Pattersonのtopographyテーマはどれも好ましいものが多いのですが、topshumが気に入りました。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "tp_tpushum" ,n=50)) 
```
![16tp-tpushum.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/8e81a9f2-5e02-877c-37ea-470dfe9bb4bd.png)

### esri_hypsometry_eu_iceland

Icelandというパレット。うむ・・・。日本の遺跡分布図で使うのは難しいかな？

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hypsometry_eu_iceland" ,n=50)) 
```
![17iceland.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/22ad91fd-154c-3db1-1544-b42627354a58.png)

### esri_hypsometry_eu_spain

Spainというパレット。平地が茶色で標高が高いほうが緑というのがスペイン的なんでしょうか？

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hypsometry_eu_spain" ,n=50)) 
```
![18spain.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/01f63cfe-5774-a0d7-13fb-91a784107951.png)

### wkp_schwarzwald_wiki_schwarzwald_cont

wikipedeaのカラーパレット群から。色みの振れ幅が大きい印象。目が痛い。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "wkp_schwarzwald_wiki_schwarzwald_cont" ,n=50) ) 
```
![19wikiSchwald.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/ca19b940-af03-644f-604a-371d6028cb57.png)

### wkp_tubs_nrwc

遺跡分布図の背景には向いていないですね。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "wkp_tubs_nrwc" ,n=50) ) 
```
![20nrwc.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/a8d00913-e426-3bb9-b4b7-8288698e7a11.png)

### arendal_zambezi_proximity

オーソドックス、地形がわかりやすい、うるさくないと三拍子揃ったカラーパレットです。好きです。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "arendal_zambezi_focus" ,n=50) ) 
```
![21zambezi.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/f00e7a80-7df0-f9be-a0c0-70f52a374e2c.png)

### wkp_country_wiki_andorra

wikipediaから。控えめで見やすいです。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "wkp_country_wiki_andorra" ,n=50) ) 
```
![22andorra.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/88f07cb9-7a3c-7c0d-3448-9f179a87a5e3.png)

### wkp_country_wiki_zurich

wikipediaから。こちらも控えめ。andorraと比べるとほんの少しコントラストが強い。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "wkp_country_wiki_zurich" ,n=50) ) 
```
![23zurich.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/297ef2d1-7f4d-1a41-ef76-a24f795b4117.png)

###esri_hypsometry_planet_mars_2

火星の地表面をイメージしたカラーパレット・・・地球上では使わなくても良いかもしれません。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hypsometry_planet_mars_2" ,n=50) ) 
```
![25mars2.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/03410189-c458-a719-3ca2-80ce2e49cece.png)

# 陰影起伏図のカラーパレット

`cptcity`には陰影図に特化したパレットも用意されています。

### esri_hillshade_blue_green_dk

少し青みがかったパレット。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_blue_green_dk" ,n=50) ) 
```
![30blue_greenDK.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/ae573e86-be1f-aaa0-89f1-8d5897182bd0.png)

### esri_hillshade_craterlake_dk

もう真っ青なパレット。使える場面があるんでしょうか？

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_craterlake_dk" ,n=50) ) 
```
![31craterlade_dk.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/65d6e9b8-eac0-eea1-a059-0d1d4a02833d.png)

### esri_hillshade_craterlake_dki

青みが強いけれど、陰影が見やすいパレット

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_craterlake_dki" ,n=50) ) 
```
![32craterlake_dki.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/5aca3ab1-17c6-66da-50cb-c9c5a46c48fa.png)


### esri_hillshade_craterlake_shore

薄ら青いけれど等高線も見やすいパレット

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_craterlake_dk" ,n=50) ) 
```
![34craterlake_shore.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/df6e23a6-f3ed-a876-0c47-4fc061f265f1.png)


### esri_hillshade_gray_dk

オーソドックスなグレースケールのパレット

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_gray_dk" ,n=50) ) 
```
![35gray_dk.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/3fc97751-0bed-6bfd-deef-78262fb2db9b.png)


### esri_hillshade_illumination

急傾斜に明るい色が割り当てられているパレット

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_illumination" ,n=50) ) 
```
![36illumination.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/ca3f739f-b04a-72dd-29fe-05818b0534a1.png)


### esri_hillshade_pnw

茶色みがわずかに入ったパレット。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_pnw" ,n=50) ) 
```
![37pnw.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/f6f26a83-69dc-88e6-84aa-f099ab46a7d5.png)

### esri_hillshade_yellow_dk

こちらも茶色みがわずかに入ったパレット。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esri_hillshade_yellow_dk" ,n=50) ) 
```
![38yellwo_dk.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/948562d2-df66-3edb-a4c4-8dc2629707f6.png)


# 傾斜区分図

slopeに特化したパレット。

### colo_alpen_ski_slope_lunch

濃緑基調のパレット。自然な感じで好ましい。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "colo_alpen_ski_slope_lunch" ,n=50) ) 
```
![41alpen_ski_lunch.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/a1c3abbe-e29e-6c4e-1437-2509a98c48cc.png)

### esdb_slopese

濃茶基調のパレット。傾斜は激しく見やすい。

```{r}
p +  scale_fill_gradientn(colours = cptcity::cpt(pal = "esdb_slopese" ,n=50) ) 
```
![42esdb.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/c40a53b0-0af4-dbbf-8160-9b18bcdaab62.png)