R 由于未知错误,无法将多边形层和点层映射在一起

R 由于未知错误,无法将多边形层和点层映射在一起,r,ggplot2,sf,tmap,R,Ggplot2,Sf,Tmap,我尝试过使用tmap和ggplot的映射函数,还利用了sf包中的一些实用程序。我有一个地块数据集,我想将其可视化为点,并将它们分层到人口普查地块数据集上,我想将其可视化为多边形层。这项工作的重点是查看不同类型普查区域的点分布情况 la_geom <- read.delim("\\Users\\edwinysun\\Box\\Terner Center_Summer 2019\\Religious Organizations\\Raw data\\LandVision\\06037\\Par

我尝试过使用tmap和ggplot的映射函数,还利用了sf包中的一些实用程序。我有一个地块数据集,我想将其可视化为点,并将它们分层到人口普查地块数据集上,我想将其可视化为多边形层。这项工作的重点是查看不同类型普查区域的点分布情况

la_geom <- read.delim("\\Users\\edwinysun\\Box\\Terner Center_Summer 2019\\Religious Organizations\\Raw data\\LandVision\\06037\\Parcels_06037.txt", sep = ",")
la$PARCEL_APN <- as.character(la$PARCEL_APN)
接下来,我尝试使用tmap和ggplot来渲染贴图。我经常遇到这样一个错误“dirname(to)中的error:path太长”。在谷歌搜索了很多次之后,我无法找出这个错误的原因,或者它对我的代码意味着什么。无论我使用上述代码还是ggplot+geom_sf()+geom_points(),都会收到相同的错误消息。我也试着换掉tm_泡泡和tm_点,但也遇到了同样的错误


有人知道如何解决这个问题,或者至少了解错误消息的意思吗?

这听起来像是一个带有长路径名的操作系统错误,但脚本中的路径并不太长,所以我感到困惑。可能与tmap使用的临时文件有关?这很难帮上忙,因为如果不访问您正在使用的文件,就不可能重现您的错误。您能否提供一个会抛出此错误的自包含可复制示例(例如,删除脚本中的绝对路径并使数据文件可用)?此外,如果您发布sessionInfo()输出,这样人们就会知道您的R配置的详细信息,这也会有所帮助?
la_geom <- left_join(la, la_geom, by = "PARCEL_APN")

ca_o_map <- read.delim("\\Users\\edwinysun\\Box\\Terner_Center_Literature and_Project_Documents\\TCAC Opportunity Maps Analysis\\2020 TCAC Code for Review\\output\\2020\\final_2020_public.csv", sep = ",")

la_o_map <- ca_o_map %>% filter(county_name == "Los Angeles") 
la_o_map <- la_o_map %>% mutate(o_map = case_when(oppcat == "High Segregation & Poverty" ~ 1,
                                                      oppcat == "Low Resource" ~ 2,
                                                      oppcat == "Moderate Resource" ~ 3,
                                                      oppcat == "High Resource" ~ 4,
                                                      oppcat == "Highest Resource" ~ 5,
                                                      TRUE ~ 0)
                                    )
la_o_map$fips <- str_pad(la_o_map$fips, width = 11, side = "left", pad = 0)
la_shp <- st_read("\\Users\\edwinysun\\Box\\Terner Center_Summer 2019\\USPS Data\\Raw data\\Shapefiles\\CENSUS_TRACTS_2010.shp")
# LA census tract shapefiles already downloaded in the USPS project. 
la_o_map <- la_o_map %>% rename(GEOID10 = fips)
la_o_map <- left_join(la_o_map, la_shp, by = "GEOID10")
la_o_map <- st_as_sf(la_o_map)

la_geom <- la_geom %>% filter(USE_CODE_STD_LPS == "9101") 
la_geom <- st_as_sf(x = la_geom, coords = c("XCOORD", "YCOORD"))
st_crs(la_geom) <- st_crs(la_o_map)
tmap_mode("view")
tm_shape(la_o_map) + tm_fill(col = "o_map") + tm_shape(la_geom) + tm_bubbles(size = "ASSR_ACREAGE")