- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有一组点,带有一个 LINEID 和一个 ID_SEQ。 LINEID 确定唯一的 Line-IDS,而 ID_SEQ 确定 Line-ID 中点的顺序。
我要转型指向线,按 ID_SEQ 排序并按 LINEID 分组 .
随包sp
,我能够达到预期的结果,但我想用 sf
做到这一点包裹。我在这里缺少什么?
这是一些虚拟数据,用 sp
说明了所需的结果-functions 和两次尝试对 sf
做同样的事情,但绘图显示不同的结果。
library(sp)
library(sf)
sfpoints <- {
structure(list(LINEID = c(4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 10L, 10L,
10L), ID_SEQ = c(10L, 11L, 12L, 13L, 16L, 30L, 31L, 32L, 33L,
34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L,
1L, 2L, 3L), geometry = structure(list(structure(c(15.423568,
47.06248), class = c("XY", "POINT", "sfg")), structure(c(15.423644,
47.062523), class = c("XY", "POINT", "sfg")), structure(c(15.423691,
47.062553), class = c("XY", "POINT", "sfg")), structure(c(15.423712,
47.06257), class = c("XY", "POINT", "sfg")), structure(c(15.423716,
47.062576), class = c("XY", "POINT", "sfg")), structure(c(15.423712,
47.062588), class = c("XY", "POINT", "sfg")), structure(c(15.423731,
47.062595), class = c("XY", "POINT", "sfg")), structure(c(15.423779,
47.062626), class = c("XY", "POINT", "sfg")), structure(c(15.423835,
47.062664), class = c("XY", "POINT", "sfg")), structure(c(15.423879,
47.062714), class = c("XY", "POINT", "sfg")), structure(c(15.423897,
47.062767), class = c("XY", "POINT", "sfg")), structure(c(15.423862,
47.062828), class = c("XY", "POINT", "sfg")), structure(c(15.423783,
47.062897), class = c("XY", "POINT", "sfg")), structure(c(15.423681,
47.062973), class = c("XY", "POINT", "sfg")), structure(c(15.423564,
47.06306), class = c("XY", "POINT", "sfg")), structure(c(15.423437,
47.063164), class = c("XY", "POINT", "sfg")), structure(c(15.42331,
47.06327), class = c("XY", "POINT", "sfg")), structure(c(15.423186,
47.063385), class = c("XY", "POINT", "sfg")), structure(c(15.423062,
47.063496), class = c("XY", "POINT", "sfg")), structure(c(15.422941,
47.063602), class = c("XY", "POINT", "sfg")), structure(c(15.422821,
47.063717), class = c("XY", "POINT", "sfg")), structure(c(15.422699,
47.063824), class = c("XY", "POINT", "sfg")), structure(c(15.422518,
47.061687), class = c("XY", "POINT", "sfg")), structure(c(15.422617,
47.06179), class = c("XY", "POINT", "sfg")), structure(c(15.422717,
47.061893), class = c("XY", "POINT", "sfg"))), class = c("sfc_POINT",
"sfc"), precision = 0, bbox = structure(c(15.422518, 47.061687,
15.423897, 47.063824), .Names = c("xmin", "ymin", "xmax", "ymax"
), class = "bbox"), crs = structure(list(epsg = NA_integer_,
proj4string = NA_character_), .Names = c("epsg", "proj4string"
), class = "crs"), n_empty = 0L)), .Names = c("LINEID", "ID_SEQ",
"geometry"), row.names = c(NA, -25L), class = c("sf", "data.frame"
), sf_column = "geometry", agr = structure(c(NA_integer_, NA_integer_
), .Names = c("LINEID", "ID_SEQ"), .Label = c("constant", "aggregate",
"identity"), class = "factor"))}
par(mfrow=c(1,3))
## SP - way
tstssp <- as(sfpoints, "Spatial")
tstssp <- SpatialLines(lapply(split(tstssp, tstssp$LINEID), function(x)
Lines(list(Line(coordinates(x))), x$LINEID[1L])))
plot(tstssp, col=1:2, lwd=3, main="SP-Desired Result")
## SF - way ???
tst <- sfpoints %>%
group_by(LINEID) %>%
st_coordinates() %>%
st_linestring()
plot(st_geometry(tst), col=1:2 , main="SF-Wrong Result")
tst <- sfpoints %>%
group_by(LINEID) %>%
summarise() %>%
st_cast("LINESTRING")
plot(st_geometry(tst), col=1:2, main="SF-Wrong Result")
最佳答案
正如 Edzer 在 this issue 中解释的那样,您必须提供summarise
带参数的函数 do_union = FALSE
:
tst <- sfpoints %>%
group_by(LINEID) %>%
summarise(do_union = FALSE) %>%
st_cast("LINESTRING")
plot(st_geometry(tst), col=1:2)
关于r - 从点创建多线,用 sf 包按 ID 分组,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/50908771/
前言: 有时候,一个数据库有多个帐号,包括数据库管理员,开发人员,运维支撑人员等,可能有很多帐号都有比较大的权限,例如DDL操作权限(创建,修改,删除存储过程,创建,修改,删除表等),账户多了,管理
所以我用 Create React App 创建并设置了一个大型 React 应用程序。最近我们开始使用 Storybook 来处理和创建组件。它很棒。但是,当我们尝试运行或构建应用程序时,我们不断遇
遵循我正在创建的控件的代码片段。这个控件用在不同的地方,变量也不同。 我正在尝试编写指令来清理代码,但在 {{}} 附近插入值时出现解析错误。 刚接触 Angular ,无法确定我错过了什么。请帮忙。
我正在尝试创建一个 image/jpeg jax-rs 提供程序类,它为我的基于 post rest 的 Web 服务创建一个图像。我无法制定请求来测试以下内容,最简单的测试方法是什么? @POST
我一直在 Windows 10 的模拟器中练习 c。后来我改用dev C++ IDE。当我在 C 中使用 FILE 时。创建的文件的名称为 test.txt ,而我给出了其他名称。请帮助解决它。 下面
当我们创建自定义 View 时,我们将 View 文件的所有者设置为自定义类,并使用 initWithFrame 或 initWithCode 对其进行实例化。 当我们创建 customUITable
我正在尝试为函数 * Producer 创建一个线程,但用于创建线程的行显示错误。我为这句话加了星标,但我无法弄清楚它出了什么问题...... #include #include #include
今天在做项目时,遇到了需要创建JavaScript对象的情况。所以Bing了一篇老外写的关于3种创建JavaScript对象的文章,看后跟着打了一遍代码。感觉方法挺好的,在这里与大家分享一下。 &
我正在阅读将查询字符串传递给 Amazon 的 S3 以进行身份验证的文档,但似乎无法理解 StringToSign 的创建和使用方式。我正在寻找一个具体示例来说明 (1) 如何构造 String
前言:我对 C# 中任务的底层实现不太了解,只了解它们的用法。为我在下面屠宰的任何东西道歉: 对于“我怎样才能开始一项任务但不等待它?”这个问题,我找不到一个好的答案。在 C# 中。更具体地说,即使任
我有一个由一些复杂的表达式生成的 ILookup。假设这是按姓氏查找人。 (在我们简单的世界模型中,姓氏在家庭中是唯一的) ILookup families; 现在我有两个对如何构建感兴趣的查询。 首
我试图创建一个 MSI,其中包含 和 exe。在 WIX 中使用了捆绑选项。这样做时出错。有人可以帮我解决这个问题。下面是代码: 错误 error LGH
在 Yii 中,Create 和 Update 通常使用相同的形式。因此,如果我在创建期间有电子邮件、密码、...other_fields...等字段,但我不想在更新期间专门显示电子邮件和密码字段,但
上周我一直在努力创建一个给定一行和一列的 QModelIndex。 或者,我会满足于在已经存在的 QModelIndex 中更改 row() 的值。 任何帮助,将不胜感激。 编辑: QModelInd
出于某种原因,这不起作用: const char * str_reset_command = "\r\nReset"; const char * str_config_command = "\r\nC
现在,我有以下由 original.df %.% group_by(Category) %.% tally() %.% arrange(desc(n)) 创建的 data.frame。 DF 5),
在今天之前,我使用/etc/vim/vimrc来配置我的vim设置。今天,我想到了创建.vimrc文件。所以,我用 touch .vimrc cat /etc/vim/vimrc > .vimrc 所
我可以创建一个 MKAnnotation,还是只读的?我有坐标,但我发现使用 setCooperative 手动创建 MKAnnotation 并不容易。 想法? 最佳答案 MKAnnotation
在以下代码中,第一个日志语句按预期显示小数,但第二个日志语句记录 NULL。我做错了什么? NSDictionary *entry = [[NSDictionary alloc] initWithOb
我正在使用与此类似的代码动态添加到数组; $arrayF[$f+1][$y][$x+1] = $value+1; 但是我在错误报告中收到了这个: undefined offset :1 问题:尝试创
我是一名优秀的程序员,十分优秀!