你能指出蒙得维的亚的shapefile吗？

Question

我想将我的分析范围限定在乌拉圭首都蒙得维的亚市。

如何对这个名为乌拉圭的主要对象进行子集化？

有蒙得维的亚 shapefile 的链接吗？

uruguay <- st_read("ine_seccen.shp")
tm_shape(uruguay) + tm_polygons()`

乌拉圭

Answer 1

让我们读入数据

uruguay <- sf::st_read("uruguay/ine_seccen.shp")
#> Reading layer `ine_seccen' from data source 
#>   `uruguay/ine_seccen.shp' using driver `ESRI Shapefile'
#> Simple feature collection with 232 features and 8 fields
#> Geometry type: MULTIPOLYGON
#> Dimension:     XY
#> Bounding box:  xmin: 366582.2 ymin: 6127919 xmax: 858252.1 ymax: 6671738
#> Projected CRS: WGS 84 / UTM zone 21S

让我们检查数据本身并绘制它：

uruguay
#> Simple feature collection with 232 features and 8 fields
#> Geometry type: MULTIPOLYGON
#> Dimension:     XY
#> Bounding box:  xmin: 366582.2 ymin: 6127919 xmax: 858252.1 ymax: 6671738
#> Projected CRS: WGS 84 / UTM zone 21S
#> First 10 features:
#>    ogc_fid        area perimeter depto seccion codsec  nombdepto cdepto_iso
#> 1        1    590829.7  3320.439     1       1    101 MONTEVIDEO       UYMO
#> 2        2    482430.9  4943.781     1       2    102 MONTEVIDEO       UYMO
#> 3        3    421108.0  3805.106     1       3    103 MONTEVIDEO       UYMO
#> 4        4    565182.3  3851.338     1       4    104 MONTEVIDEO       UYMO
#> 5        5    840787.2  4471.698     1       5    105 MONTEVIDEO       UYMO
#> 6        6    542797.2  3061.679     1       6    106 MONTEVIDEO       UYMO
#> 7        7   1222503.8  5869.490     1       7    107 MONTEVIDEO       UYMO
#> 8        8   2080318.8  6619.309     1       8    108 MONTEVIDEO       UYMO
#> 9        9 124589920.4 75168.524     1       9    109 MONTEVIDEO       UYMO
#> 10      10  27682916.7 31454.642     1      10    110 MONTEVIDEO       UYMO
#>                          geometry
#> 1  MULTIPOLYGON (((572996.2 61...
#> 2  MULTIPOLYGON (((572148.3 61...
#> 3  MULTIPOLYGON (((572708.3 61...
#> 4  MULTIPOLYGON (((573700.7 61...
#> 5  MULTIPOLYGON (((574375.3 61...
#> 6  MULTIPOLYGON (((574448.7 61...
#> 7  MULTIPOLYGON (((574382.2 61...
#> 8  MULTIPOLYGON (((574342.4 61...
#> 9  MULTIPOLYGON (((570567.2 61...
#> 10 MULTIPOLYGON (((578956.6 61...

tmap::tm_shape(uruguay) + 
  tmap::tm_polygons()

从上面的数据片段中我们可以看到

nombdepto

列包含名称。让我们只过滤掉那些

MONTEVIDEO

边界的多边形。

mv <- uruguay |>
  subset(nombdepto == "MONTEVIDEO")

tmap::tm_shape(mv) + 
  tmap::tm_polygons()

或者，直接在

tm_shape()

函数中过滤它，提供我们的

sf

对象的子集：

tmap::tm_shape(uruguay[uruguay$nombdepto == "MONTEVIDEO", ]) +
  tmap::tm_polygons()

^{创建于 2023-12-27，使用 reprex v2.0.2}

PS。在将来/其他问题中，请提供代表性数据或链接所使用的原始数据文件。上面示例中使用的

ine_seccen.shp

取自https://visualizador.ide.uy/geonetwork/srv/api/records/c14e6a59-756c-44f1-bc4e-426e2573122b（Secciones Censales del Uruguay 2011 WFS 服务器层） ).