ggplot指定经度/纬度轴刻度值

huangapple
go评论110阅读模式
英文:

ggplot specify longitude/latitude axis breaks

问题

我正在使用以下代码来创建一个地图，该地图基于一个shapefile中的一个'trial_area'，并在一个第二个shapefile中叠加了一个海岸线和'prod_areas'。然后，我使用coord_sf来缩放地图以适应trial_area多边形的st_bbox。然而，对于某些区域（参见下面的示例），轴刻度文本标签最终会重叠在一起，是否有一种方法可以指定轴刻度标记的间隔以避免这种情况（例如，纬度为0.1，经度为0.5）？

  1. poly <- trial_areas %>%
  2.   filter(Id==5) 
  3. ext <- st_bbox(poly)
  5. plot_SoundOfSleat <- ggplot() +
  6. theme(panel.background = element_rect(fill = 'light blue'), element_line()) +
  7. geom_sf(data=poly) +
  8. geom_sf(data=prod_areas, fill=mycol) +
  9. geom_sf(data = Scot, aes(),
  10.         fill = "lightgreen", col="darkgreen") +
  11. coord_sf(xlim = c(ext[1], ext[3]), ylim = c(ext[2], ext[4]))  +
  12. ggtitle("Sound of Sleat Trial Area 5") +
  13. geom_sf_text(aes(label = Producti_1), data=prod_areas, size=3, hjust=0, vjust=0) +
  14. labs(x = "Longitude", y= "Latitude")
  16. plot_SoundOfSleat

ggplot指定经度/纬度轴刻度值

  2. <details>
  3. <summary>英文:</summary>
  5. I am using the following code to create a map of one &#39;trial_area&#39; from a shapefile and overlaying a coastline and &#39;prod_areas&#39; from a second shapefile. I&#39;m then using the coord_sf to zoom the map to the st_bbox of the trial_area polygon. However for some areas (see the example below) the axis tick text labels end up overlapping, is there a way that i can specify the axis tick mark intervals to avoid this (e.g. .1 for latitude and .5 for longitude)? 
  7.     poly &lt;- trial_areas %&gt;% 
  8.     filter(Id==5) 
  9.     ext &lt;- st_bbox(poly)
  10.   
  11.     plot_SoundOfSleat &lt;- ggplot() +
  12.     theme(panel.background = element_rect(fill = &#39;light blue&#39;),element_line()) +
  13.     geom_sf(data=poly)+
  14.     geom_sf(data=prod_areas,fill=mycol) +
  15.     geom_sf(data = Scot, aes(),
  16.             fill = &quot;lightgreen&quot;,col=&quot;darkgreen&quot;) +
  17.     coord_sf(xlim = c(ext[1], ext[3]), ylim = c(ext[2], ext[4]))  +
  18.     ggtitle(&quot;Sound of Sleat Trial Area 5&quot;) +
  19.     geom_sf_text(aes(label = Producti_1), data=prod_areas,size=3,hjust=0, vjust=0) +
  20.     labs(x = &quot;Longitude&quot;, y= &quot;Latitude&quot;)
  22.     plot_SoundOfSleat
  24. [![enter image description here][1]][1]
  27.   [1]: https://i.stack.imgur.com/iZeme.png
  29. </details>
  32. # 答案1
  33. **得分**: 6
  35. `geom_sf()` 应该能够与 `ggplot2::scale_*_continuous()` 无缝配合使用，您可以使用 `breaks = ` 参数。请注意西经度，因为它们在数据中是负数，但在标签中是正数。
  37. 以下是一些示例：
  39. ```R
  40. library(sf)
  42. # 样本数据
  43. nc <- st_read(system.file("shape/nc.shp", package="sf"))
  45. # 不编辑格网线
  46. nc_1 <- ggplot(nc) + 
  47.           geom_sf() +
  48.           ggtitle('原始')
  50. nc_2 <- ggplot(nc) +
  51.           geom_sf() + 
  52.           scale_y_continuous(breaks = c(34, 35, 36)) + 
  53.           scale_x_continuous(breaks = seq(-84, -76, by = 1)) +
  54.           ggtitle('更少的纬度，更多的经度')
  56. nc_3 <- ggplot(nc) +
  57.           geom_sf(data = st_graticule(nc, 
  58.                                       lat = seq(34, 36, by = 1),
  59.                                       lon = seq(-84, -76, by = 4)), 
  60.                   color = 'orange') +
  61.           geom_sf() +
  62.           coord_sf(datum = NA) +
  63.           ggtitle('使用 st_graticule')
  65. # 使用 cowplot 输出单个图像，而不是3个
  66. cowplot::plot_grid(nc_1, nc_2, nc_3, ncol = 1)

您应该能够使用 st_bbox 的输出来自动化合理数量的格网线（网格线），如果有必要。

英文:

geom_sf() should work seamlessly with ggplot2::scale_*_continuous(), where you can use the breaks = argument. Be careful with west longitudes, as they're negative numbers in the data but positive in the labels.

I've included some examples below:

  1. library(sf)
  3. # sample data
  4. nc &lt;-  st_read(system.file(&quot;shape/nc.shp&quot;, package=&quot;sf&quot;))
  6. # No edits to graticules
  7. nc_1 &lt;- ggplot(nc) + 
  8.           geom_sf() +
  9.           ggtitle(&#39;original&#39;)
  11. nc_2 &lt;- ggplot(nc) +
  12.           geom_sf() + 
  13.           scale_y_continuous(breaks = c(34, 35, 36)) + 
  14.           scale_x_continuous(breaks = seq(-84, -76, by = 1)) +
  15.           ggtitle(&#39;fewer lat, more lon&#39;)
  17. nc_3 &lt;- ggplot(nc) +
  18.           geom_sf(data = st_graticule(nc, 
  19.                                       lat = seq(34, 36, by = 1),
  20.                                       lon = seq(-84, -76, by = 4)), 
  21.                   color = &#39;orange&#39;) +
  22.           geom_sf() +
  23.           coord_sf(datum = NA) +
  24.           ggtitle(&#39;using st_graticule&#39;)
  26. # using cowplot to output single image, rather than 3
  27. cowplot::plot_grid(nc_1, nc_2, nc_3, ncol = 1)

ggplot指定经度/纬度轴刻度值

You should be able to use the output from st_bbox to automate a sensible number of graticules (grid lines) if necessary.

huangapple
  • 本文由 发表于 2020年1月6日 18:26:51
  • 转载请务必保留本文链接：https://go.coder-hub.com/59610405.html
  • r
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定