gpt4 book ai didi

根据 NPA(地区)R map 瑞士

转载 作者:行者123 更新时间:2023-12-04 07:39:20 31 4
gpt4 key购买 nike

我打算在瑞士做一个调查。将询问 NPA。

NPA (postal codes)包含 4 个数字。

  • 例如 1227 是 Carouge 的 NPA(瑞士日内瓦州的一部分)。
  • 例如,1784 年是 Courtepin 的 NPA(弗里堡州 - 瑞士的一部分)。

  • 我想知道如何在 map 上表示所有观察(大约 1500)。我正在考虑将 ggplot 用于其他图形(我认为 ggplot 是“漂亮的”)。但是,我愿意接受任何其他建议。

    以下是一些假数据:
    http://pastebin.com/HsuQnLP3

    瑞士 map 的输出应该有点像美国 map (来源: http://www.openintro.org)

    enter image description here

    更新:

    我试图创建一些代码:
    library(sp)
    test <- url("https://dl.dropboxusercontent.com/u/6421260/CHE_adm3.RData")
    print(load(test))
    close(test)

    gadm$NAME_3
    gadm$TYPE_3

    不过好像 http://gadm.org/不提供公社的 NPA...

    新更新:

    我找到了(感谢@yrochat)一个带有 NPA 的 shapefile:
    http://www.cadastre.ch/internet/cadastre/fr/home/products/plz/data.html

    我是名为:Shape LV03 的 ZIP 文件

    然后我试过了
    library("maptools")
    swissmap <- readShapeLines("C:/Users/yourName/YourPath/PLZO_SHP_LV03/PLZO_PLZ.shp")
    plot(swissmap)
    data <- data.frame(swissmap)
    data$PLZ #the row who gives the NPA

    由于我在 shapefile 上有 PLZ,我如何在 map 上为我的观察着色?
    我提供了一些关于数据的假数据 http://pastebin.com/HsuQnLP3

    enter image description here

    谢谢

    最佳答案

    好的,使用 shapefile,我们可以轻松地绘制内容。

    work.dir <- "directory_name_no_trailing slash"

    # open the shapefile
    require(rgdal)
    require(rgeos)
    require(ggplot2)
    ch <- readOGR(work.dir, layer = "PLZO_PLZ")

    # convert to data frame for plotting with ggplot - takes a while
    ch.df <- fortify(ch)

    # generate fake data and add to data frame
    ch.df$count <- round(runif(nrow(ch.df), 0, 100), 0)

    # plot with ggplot
    ggplot(ch.df, aes(x = long, y = lat, group = group, fill = count)) +
    geom_polygon(colour = "black", size = 0.3, aes(group = group)) +
    theme()

    # or you could use base R plot
    ch@data$count <- round(runif(nrow(ch@data), 0, 100), 0)
    plot(ch, col = ch@data$count)

    我个人觉得 ggplotplot 更容易使用并且默认输出更好看。

    screenshot

    ggplot使用简单的数据框,这使得子集化变得容易。
    # plot just a subset of NPAs using ggplot
    my.sub <- ch.df[ch.df$id %in% c(4,6), ]
    ggplot(my.sub, aes(x = long, y = lat, group = group, fill = count)) +
    geom_polygon(colour = "black", size = 0.3, aes(group = group)) +
    theme()

    结果:

    screenshot2

    关于根据 NPA(地区)R map 瑞士,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/17757281/

    31 4 0
    Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
    广告合作:1813099741@qq.com 6ren.com