使用Geopandas将NYC shapefile替换为Manhattan

Question

我有一个NYC的shapefile，我想使用geopandas缩小为仅覆盖曼哈顿。 （来源：https://www1.nyc.gov/site/tlc/about/tlc-trip-record-data.page）。

谢谢您的建议！

已经尝试过https://mygeodata.cloud/，但无法使其正常工作...

我的总体目标是稍后检查某些经纬度点是否在Manhattan的shapefile中，但目前，它将对NYC进行此操作，从而包括实际上在曼哈顿以外的点。

Answer 1

假设您正在使用以下shapefile：

除了GeoDataFrame中的“几何”列外，它的行为类似于熊猫DataFrame。因此，您可以像对Pandas一样对GeoDataFrame进行子集化。

import geopandas as gpd

gdf = gpd.read_file('taxi_zones.shp')

# only keep Manhattan
gdf = gdf[gdf.borough == 'Manhattan']

# you can check it:
gdf.borough.value_counts()