### 使用Shapely和Fiona判断Polygon和MultiPolygon中的岛屿 在处理GeoJSON或Shapefile文件时，可以通过分析`Polygon`和`MultiPolygon`对象来识别哪些几何体代表岛屿。通常情况下，“岛屿”是指位于主陆地区域内部的小型封闭区域，在GIS术语中也称为内环或多边形内的孔。 对于`Polygon`对象而言，如果存在内部的空洞，则这些空洞可能被视作湖泊或其他非陆地特征；而对于`MultiPolygon`来说，独立存在的较小多边形可能是岛屿。具体实现方式如下： #### 判断Polygon中的岛屿 当创建一个带有空洞的`Polygon`时，可以利用其构造函数接收两个列表作为参数：一个是外轮廓点集，另一个是包含零个或多个内环（即空洞）点集合的列表[^4]。因此，可以直接通过访问该属性来获取潜在的“岛屿”。 ```python from shapely.geometry import Polygon # 定义一个多边形及其内部的一个空洞 pol = Polygon( [(0, 0), (4, 1), (5, 3), (0, 7)], [[(1, 1), (1, 2), (2, 2), (2, 1)]] ) islands = pol.interiors # 获取所有的内环（假设为空洞） for island in islands: print(f"Island coordinates: {list(island.coords)}") ``` #### 处理MultiPolygon中的岛屿 针对`MultiPolygon`结构，由于它是由若干个不相连的`Polygon`组成，所以需要遍历每一个子多边形，并依据面积大小或者其他逻辑条件判定是否属于岛屿。这里提供了一种基于面积阈值的方法来进行初步筛选[^1]。 ```python from shapely.ops import transform import pyproj def is_island(polygon, threshold_area=1e6): # 面积单位取决于输入坐标系，默认采用平方米 project = partial(pyproj.transform, pyproj.Proj(init='epsg:4326'), # WGS84经纬度投影 pyproj.Proj('+proj=merc')) # Web墨卡托投影 projected_polygon = transform(project, polygon) return projected_polygon.area < threshold_area from shapely.geometry import MultiPolygon polygons = [ Polygon([(0, 0), (1, 1), (1, 0)]), Polygon([(1, 1.1), (1, 2), (2, 2)]) ] multi_poly = MultiPolygon(polygons) possible_islands = [] for geom in multi_poly.geoms: if isinstance(geom, Polygon) and not geom.is_empty: if is_island(geom): possible_islands.append(geom) print("Possible Islands:") for i, isl in enumerate(possible_islands, start=1): print(f"{i}: Coordinates={list(isl.exterior.coords)}, Area={isl.area}") ``` 上述代码片段展示了如何使用Shapely库结合简单的面积计算方法来检测可能的岛屿。需要注意的是，实际应用中还需要考虑更多因素，比如地理位置、行政划分等信息才能更精确地标记岛屿。