我需要从geojson获取坐标数组的简单方法

要从GeoJSON文件中提取坐标数组，你可以使用Python的json模块来解析GeoJSON数据，并提取所需的坐标信息。GeoJSON是一种用于表示地理数据的标准格式，通常包含几何对象（如点、线、多边形）及其坐标。

以下是一个详细的示例，展示了如何从GeoJSON文件中提取坐标数组。

示例GeoJSON文件

假设你有一个名为example.geojson的GeoJSON文件，内容如下：

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "geometry": {
        "type": "Point",
        "coordinates": [102.0, 0.5]
      },
      "properties": {
        "name": "Sample Point"
      }
    },
    {
      "type": "Feature",
      "geometry": {
        "type": "LineString",
        "coordinates": [
          [102.0, 0.0],
          [103.0, 1.0],
          [104.0, 0.0],
          [105.0, 1.0]
        ]
      },
      "properties": {
        "name": "Sample Line"
      }
    },
    {
      "type": "Feature",
      "geometry": {
        "type": "Polygon",
        "coordinates": [
          [
            [100.0, 0.0],
            [101.0, 0.0],
            [101.0, 1.0],
            [100.0, 1.0],
            [100.0, 0.0]
          ]
        ]
      },
      "properties": {
        "name": "Sample Polygon"
      }
    }
  ]
}

示例代码

以下是一个Python脚本，展示了如何从上述GeoJSON文件中提取所有几何对象的坐标数组。

import json

def extract_coordinates(geojson):
    coordinates = []

    for feature in geojson['features']:
        geometry = feature['geometry']
        if geometry['type'] == 'Point':
            coordinates.append(geometry['coordinates'])
        elif geometry['type'] in ['LineString', 'MultiPoint']:
            coordinates.extend(geometry['coordinates'])
        elif geometry['type'] in ['Polygon', 'MultiLineString']:
            for coord in geometry['coordinates']:
                coordinates.extend(coord)
        elif geometry['type'] == 'MultiPolygon':
            for polygon in geometry['coordinates']:
                for coord in polygon:
                    coordinates.extend(coord)
    
    return coordinates

# 读取GeoJSON文件
with open('example.geojson', 'r') as file:
    geojson_data = json.load(file)

# 提取坐标数组
coordinates = extract_coordinates(geojson_data)

# 打印坐标数组
print(coordinates)

详细解释

导入json模块：
- 使用import json导入Python的json模块，用于解析GeoJSON数据。
定义extract_coordinates函数：
- 该函数接受一个GeoJSON对象，并提取所有几何对象的坐标数组。
- 根据几何对象的类型（如Point、LineString、Polygon等），提取相应的坐标。
- 将所有坐标存储在一个列表中并返回。
读取GeoJSON文件：
- 使用with open('example.geojson', 'r') as file打开GeoJSON文件，并使用json.load(file)解析文件内容。
提取坐标数组：
- 调用extract_coordinates函数，传入解析后的GeoJSON数据，获取坐标数组。
打印坐标数组：
- 使用print(coordinates)打印提取的坐标数组。