我最近开始在python中使用Geopandas来完成我的一些空间工作,并对此感到非常满意 - 我目前正在尝试阅读PostGIS功能,并且不太了解如何参数化数据库连接,这似乎并不清楚在文档中:
GeoDataFrame.from_postgis(sql, con, geom_col='geom', crs=None, index_col=None,
coerce_float=True, params=None)
Run Code Online (Sandbox Code Playgroud)
这可能是一个非常简单的问题,所有我想知道的是需要进入'con' - 我假设一个包含数据库连接信息的字符串?但是以什么格式?设置'sql'似乎很简单.任何帮助非常感谢 - 谢谢!
我想使用 R 中的 sf 包在 R 中的一个简单特征数据框中找到与另一个简单特征数据框中一组点最近的多边形。我一直在“st_join”语句中使用“st_is_within_distance”,但是这将返回给定距离内的所有内容,而不仅仅是最近的特征。
以前,我使用了 'rgeos' 包中的 'gDistance' 和这样的 'sp' 功能:
m = gDistance(a, b, byid = TRUE)
row = apply(m, 2, function(x) which(x == min(x)))
labels = unlist(b@data[row, ]$NAME)
a$NAME <- labels
Run Code Online (Sandbox Code Playgroud)
我想将这种使用 rgeos 和 sp 为一组点查找最近特征的方法转换为使用 sf。任何意见或建议非常感谢。
我习惯在 R ioslides_presentation 文档中制作两列幻灯片,但我找不到任何可以使用 R slidy_presentations 将幻灯片分成两列的方法,这是我在使用 ioslides、beamer 后第一次尝试,和转速。一定有办法吗?
我正在尝试在 R leaflet 包中添加 WMS 切片 - 使用此示例 geoserver WMS 不是问题:
leaflet() %>% addTiles() %>% setView(-93.65, 42.0285, zoom = 4) %>% addWMSTiles(
"http://sedac.ciesin.columbia.edu/geoserver/wms",
layers = "energy:energy-pop-exposure-nuclear-plants-locations_plants",
options = WMSTileOptions(format = "image/png", transparent = TRUE),
tileOptions(tms = TRUE),
attribution = "")
Run Code Online (Sandbox Code Playgroud)
然而,当我尝试使用 National Map中的 WMS 时,尽管多次尝试正确设置 url 和图层的参数,但我仍然得到空的传单结果:
leaflet() %>% addTiles() %>% setView(-93.65, 42.0285, zoom = 4) %>%addWMSTiles(
"http://basemap.nationalmap.gov/arcgis/services/USGSHydroNHD/MapServer/WmsServer?",
layers = "0",
options = WMSTileOptions(format = "image/png", transparent = TRUE),
attribution = "")
Run Code Online (Sandbox Code Playgroud)
我之前没有在 R leaflet 包之外使用过 leaflet,所以在使用这种类型的 WMS 在 …