为什么提供的tidycensus区域与sf :: _ st_area（）计算的不同。

Question

我正在使用tidycensus r软件包来吸收人口普查数据和几何形状。我希望能够计算人口密度，并使结果与我在。我注意到使用SF软件包SF :: ST_AREA（）函数与我计算的内容相比，从Tidycenus返回的地理变量之间有所不同。

library(tidyverse)
library(tidycensus)
census_api_key("my_api_key")
library(sf)
options(tigris_use_cache = TRUE)

pop_texas <-
  get_acs(geography = 'state',
      variables = "B01003_001", # Total Population
      year = 2020, 
      survey = 'acs5', 
      keep_geo_vars = TRUE,
      geometry = TRUE) %>%
filter(GEOID == '48') # Filter to Texas

由于我将keep_geo_vars参数包括为true，因此它返回了一个aland列，我认为这是返回平方米（M^2）的地理区域的正确区域。

> pop_texas$ALAND %>% format(big.mark=",")
[1] "676,680,588,914"

# Conversion to square miles
> (pop_texas$ALAND / 1000000 / 2.5899881) %>% format(big.mark=",")
[1] "261,267.8"

当我将Aland数量转换为平方英里时，我得到的数字与Censusreporter.org上所示的数字相同：

我还尝试使用sf :: st_area（）函数来计算该区域，但是我得到一个不同的结果：

> sf::st_area(pop_texas) %>% format(big.mark=",", scientific=FALSE)
[1] "688,276,954,146 [m^2]"

# Conversion to square miles
> (sf::st_area(pop_texas) / 1000000 / 2.5899881) %>%
+   as.numeric() %>%
+   format(big.mark=",", scientific=FALSE)
[1] "265,745.2"

请让我知道是否有一些我缺少来调和这些数字的东西。我希望直接通过平局直接获得相同的结果，或使用sf :: st_area（）计算区域。

现在，我已经离开了很多：

> (pop_texas$ALAND - as.numeric(st_area(pop_texas)) ) %>%
+   format(big.mark=",")
[1] "-11,596,365,232"

Answer 1

如果您想要像德克萨斯州这样的形状的“官方”区域，则应始终使用aland或已发布的区域值。 ST_AREA（）正在使用几何形状来计算多边形的面积，这将始终是德克萨斯州（或任何其他区域）的简化和不完美的表示。对于较小的形状（例如人口普查区），计算可能会非常接近；对于更大的状态（尤其是那些具有复杂沿海地区的国家，如德克萨斯州），您将越来越远。

Answer 2

这些差异通常是由于CRS（在您的sf对象上使用的投影）。一些预测扭曲区域，其他投影会扭曲形状。请参阅此信息以了解更多 http://wiki.gis.com/wiki/wiki/wiki/index.php/distortion#:~: text = text = 20个方向％2C％20型％20和％20 rea 。