geojson，一个有趣的 Python 库！

大家好，我是涛哥，本文内容来自 涛哥聊Python ，转载请标原创。

今天为大家分享一个有趣的 Python 库 - geojson。

Github地址：https://github.com/jazzband/geojson

---

GeoJSON是一种基于JSON的地理空间数据交换格式，广泛用于多种网络应用中，尤其是地理数据的可视化。Python的GeoJSON库提供了一个简单的接口来构造、解析和处理GeoJSON数据，非常适合用于地理信息系统（GIS）、地图创建和位置数据分析等领域。

安装

安装Python GeoJSON库非常简单，可以通过Python的包管理器pip进行安装：

pip install geojson

这条命令将安装GeoJSON库及其所有依赖。

特性

• 简洁的API：易于理解和使用的API，允许开发者快速地创建和操作GeoJSON对象。
• 广泛的格式支持：支持所有GeoJSON格式，包括点（Point）、线（LineString）、多边形（Polygon）、多点（MultiPoint）、多线（MultiLineString）和多多边形（MultiPolygon）等。
• 与其他Python库集成：可以与Pandas、Shapely等数据处理和地理处理库无缝集成。

基本功能

创建GeoJSON对象

创建基本的GeoJSON对象非常直接：

from geojson import Point, Feature, FeatureCollection

# 创建一个点
point = Point((-115.81, 37.24))

# 创建一个特征
feature = Feature(geometry=point, properties={"city": "Las Vegas"})

# 创建一个特征集合
feature_collection = FeatureCollection([feature])

print(feature_collection)

这段代码创建了一个点对象，然后将其封装在一个特征中，并最终放入一个特征集合中。输出将是GeoJSON格式的数据。

解析GeoJSON数据

GeoJSON库还可以解析GeoJSON字符串，将其转换为Python对象，以便进一步处理：

import geojson

data = '''
{
  "type": "Feature",
  "geometry": {
    "type": "Point",
    "coordinates": [-115.81, 37.24]
  },
  "properties": {
    "city": "Las Vegas"
  }
}
'''

feature = geojson.loads(data)
print(feature.geometry)
print(feature['properties']['city'])

高级功能

GeoJSON库提供了一系列高级功能，使得开发者可以执行更为复杂的地理空间数据操作和转换。

自定义属性和批量生成

GeoJSON库允许开发者自定义地理特征的属性，以及批量生成复杂的地理数据结构。

from geojson import Feature, Point, FeatureCollection

# 批量创建特征
features = [Feature(geometry=Point((lon, lat)), properties={"city": city})
            for lon, lat, city in [(-74.00597, 40.71427, "New York"),
                                   (-0.127758, 51.507351, "London"),
                                   (139.6917, 35.6895, "Tokyo")]]

feature_collection = FeatureCollection(features)
print(feature_collection)

这个例子中，创建了一个特征集合，包括世界上几个著名城市的点位。这种方式非常适合处理大量的地理位置数据。

地理数据转换

GeoJSON库可以与其他Python库如Shapely一起使用，以执行地理数据之间的转换和操作。

from geojson import Feature, Polygon
from shapely.geometry import shape, mapping

# 创建一个GeoJSON多边形
polygon = Polygon([[(2.0, 2.0), (2.0, 8.0), (8.0, 8.0), (8.0, 2.0), (2.0, 2.0)]])
feature = Feature(geometry=polygon)

# 使用Shapely转换GeoJSON多边形
shapely_polygon = shape(feature.geometry)
buffered_polygon = shapely_polygon.buffer(0.5)  # 创建缓冲区
buffered_geojson = mapping(buffered_polygon)  # 转换回GeoJSON

print(buffered_geojson)

在这个例子中，首先创建了一个GeoJSON多边形，然后利用Shapely库对其进行了缓冲区处理，并将结果转换回GeoJSON格式。

动态属性计算

GeoJSON库可以用于动态计算和更新地理特征的属性，这在地理数据分析中非常有用。

from geojson import Feature, Point

def calculate_population_density(population, area):
    return population / area if area != 0 else 0

# 创建一个带有动态计算属性的特征
feature = Feature(geometry=Point((13.4050, 52.5200)),
                  properties={"city": "Berlin",
                              "population": 3748148,
                              "area": 891.8})

# 计算并更新人口密度属性
feature['properties']['population_density'] = calculate_population_density(
    feature['properties']['population'], feature['properties']['area'])

print(feature)

这个例子中，为柏林创建了一个地理特征，并计算了其人口密度属性。

实际应用场景

GeoJSON库的功能使其非常适用于多种实际应用，从地理信息系统（GIS）到地图创建和动态数据可视化等领域。

地图可视化

GeoJSON是网页地图服务如Leaflet或Mapbox的标准数据格式之一，用于网页地图的创建和数据可视化。

from geojson import Feature, Point, FeatureCollection

# 创建多个地理位置特征
features = [
    Feature(geometry=Point((-0.1276, 51.5074)), properties={"name": "London"}),
    Feature(geometry=Point((-73.935242, 40.730610)), properties={"name": "New York"}),
    Feature(geometry=Point((139.6917, 35.6895)), properties={"name": "Tokyo"})
]

# 生成特征集合
feature_collection = FeatureCollection(features)

# 这个特征集合可以直接被地图服务读取并显示在网页地图上
print(feature_collection)

这个示例创建了一个包含几个世界大城市的特征集合，可以直接被网页地图库使用来在地图上展示这些城市。

GIS数据处理

GeoJSON库可用于执行复杂的GIS数据处理任务，如过滤、转换和地理空间分析。

from geojson import Feature, Polygon, FeatureCollection
import shapely.geometry
import shapely.ops

# 定义两个多边形
polygon1 = Polygon([[(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]])
polygon2 = Polygon([[(0.5, 0.5), (1.5, 0.5), (1.5, 1.5), (0.5, 1.5), (0.5, 0.5)]])

# 使用Shapely计算两个多边形的交集
shapely_poly1 = shapely.geometry.asShape(polygon1)
shapely_poly2 = shapely.geometry.asShape(polygon2)
intersection = shapely.ops.unary_union([shapely_poly1, shapely_poly2])

# 将计算结果转换回GeoJSON格式
intersection_geojson = shapely.geometry.mapping(intersection)
intersection_feature = Feature(geometry=intersection_geojson)

print(intersection_feature)

这个示例展示了如何计算两个多边形的交集，并将结果转换为GeoJSON格式。

动态数据跟踪和分析

在需要跟踪动态地理数据（如车辆、船只或人员位置）的应用中，GeoJSON可用于存储和分析这些数据。

from geojson import Point, Feature, FeatureCollection
import datetime

# 假设这是从设备接收到的实时位置数据
positions = [
    {"longitude": -0.1276, "latitude": 51.5074, "timestamp": datetime.datetime.now()},
    {"longitude": -73.935242, "latitude": 40.730610, "timestamp": datetime.datetime.now()}
]

# 创建实时位置特征
features = [
    Feature(geometry=Point((pos["longitude"], pos["latitude"])),
            properties={"timestamp": pos["timestamp"].isoformat()})
    for pos in positions
]

# 生成特征集合
feature_collection = FeatureCollection(features)

print(feature_collection)

这个示例处理实时位置数据，每个位置数据都被封装成一个特征，并带有时间戳，适合进行进一步的时间序列分析和可视化。

总结

Python的GeoJSON库是处理地理空间数据的强大工具，它提供了创建、解析和处理GeoJSON数据的简洁接口。通过这个库，开发者可以轻松实现地理数据的各种操作，如地图可视化、GIS数据处理和动态数据分析等。支持所有GeoJSON对象类型，包括点、线、多边形以及它们的多重形式，GeoJSON库使得地理数据的集成和应用变得更加方便快捷。无论是开发复杂的地理信息系统、实时位置跟踪应用还是简单的数据可视化项目，GeoJSON库都能提供必要的功能支持，帮助开发者有效地处理和展示地理空间信息。