小编and*_*bri的帖子

我想用leaflet()in 实现一个特殊的WMS层R.但我对我喜欢使用的WMS层没有成功.

library("leaflet")
library("sp")

# this example from the tutorial works
leaflet() %>% addTiles() %>% setView(-93.65, 42.028, zoom = 4) %>%
  addWMSTiles(
    "http://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r.cgi",
    layers = "nexrad-n0r",
    options = WMSTileOptions(format = "image/png", transparent = TRUE) 
)

# but for my implementation it don't work
leaflet() %>% addTiles() %>% setView(lat = 53.8, lng = 12.6, zoom = 8) %>% 
  addWMSTiles("http://www.umweltkarten.mv-regierung.de/script/mv_a3_wasserschutz_wms.php?", 
              layers = "t2_bib_p", 
              options = WMSTileOptions(format = "image/png")
)

我认为坐标参考系统存在问题.因为我喜欢使用的WMS层,没有标准的CRS.所以我尝试用这个来实验WMSTileOptions():

crs=CRS("+init=epsg:5650")

我们的想法是使用WMS-Layers作为我的矢量数据的背景图像.矢量数据主要来自形状文件和GPS数据(文本文件).目的是创建我可以在报告中使用的地图.如果你知道leaflet在R中使用WMS的其他可能性,也欢迎吗？

我喜欢用它来准备情节,添加水平线,最后绘制条形图.当我像下面这样做时,R将每个东西都绘制两次,这看起来不太好.对于正常情节,我搜索的等价物是:plot(NULL, xlim=c(1,2), ylim=c(3,5), axes=F)

data = c(1,4,2) 
barplot(data)
abline(h=seq(0,5,1), col="red")
barplot(data, add=T)

我更喜欢有一个基础解决方案.

我想知道0周围的人数1.但如果有多个0不被中断,1它只计为一个.

string <- "1101000010101111001110"

这是我能做的最接近的:

length(gregexpr(pattern ="101",string)[[1]])

预期产量:

5