将管道R中的经纬度坐标转换为UTM

问题描述

我具有海岸线的df空间值，但它们的纬度/经度，我需要将其位于NAD83 / UTM区域10N中。我可以通过在dplyr管道中进行更改来做到这一点吗？

testdf <- data.frame(long = c(-124.0048,-123.9844,-123.9691,-123.9604,-123.9810,-123.9612),lat = c(45.04352,45.10493,45.20530,45.29999,45.34960,45.40917))


testdf %>% mutate(utm = rgdal::project(c(lon,lat),proj = ???),# but then I guess I would need to separate the utm column into x and y still.

还是用map()？

编辑：对不起，我现在意识到这是一个不好的代表。实际数据集具有投影+proj=longlat +datum=wgs84 +no_defs。立即将其编辑到testdf中。

解决方法

您可以使用this软件包：

public interface Foo {
    public default String fooName() {
        return this.getClass().getName();
    }
}

^{由reprex package（v0.3.0）于2020-08-12创建}

关联的管道为：

int octet[8] = {0,1,0};