探索GeoPandas：轻松处理地理空间数据的开源利器

最新推荐文章于 2024-09-27 10:11:28 发布

stjklkjhgffxw

最新推荐文章于 2024-09-27 10:11:28 发布

阅读量346

点赞数 3

文章标签： python

本文链接：https://blog.csdn.net/stjklkjhgffxw/article/details/142396146

版权

## 引言

在现代数据分析中，地理空间数据的处理变得越来越重要。GeoPandas 是一个强大的开源库，旨在简化 Python 中的地理空间数据操作。本篇文章将介绍如何使用 GeoPandas 处理地理数据，并结合 LLM 应用领域探讨其潜力。

## 主要内容

### 1. GeoPandas 简介

GeoPandas 扩展了 Pandas 的数据类型，使其能够对几何类型进行空间操作。它利用 Shapely 进行几何计算，Fiona 进行文件访问，以及 Matplotlib 进行可视化。

### 2. 安装必要的库

确保你已经安装了必要的依赖库：

```bash
%pip install --upgrade --quiet sodapy
%pip install --upgrade --quiet pandas
%pip install --upgrade --quiet geopandas

3. 创建 Geopandas 数据框

下面的代码示例展示了如何从 Open City Data 创建一个 GeoPandas 数据框：

import ast
import geopandas as gpd
import pandas as pd
from langchain_community.document_loaders import OpenCityDataLoader

# 加载开放城市数据
dataset = "tmnf-yvry"  # 旧金山犯罪数据
loader = OpenCityDataLoader(city_id="data.sfgov.org", dataset_id=dataset, limit=5000) 
docs = loader.load()

# 将字典列表转换为数据框
df = pd.DataFrame([ast.literal_eval(d.page_content) for d in docs])

# 提取纬度和经度
df["Latitude"] = df["location"].apply(lambda loc: loc["coordinates"][1])
df["Longitude"] = df["location"].apply(lambda loc: loc["coordinates"][0])

# 创建 GeoPandas 数据框
gdf = gpd.GeoDataFrame(
    df, geometry=gpd.points_from_xy(df.Longitude, df.Latitude), crs="EPSG:4326"
)

# 过滤旧金山有效经纬度
gdf = gdf[
    (gdf["Longitude"] >= -123.173825)
    & (gdf["Longitude"] <= -122.281780)
    & (gdf["Latitude"] >= 37.623983)
    & (gdf["Latitude"] <= 37.929824)
]

4. 数据可视化

使用 Matplotlib 和 GeoPandas 可视化数据：

import matplotlib.pyplot as plt

# 加载旧金山地图数据
sf = gpd.read_file("https://data.sfgov.org/resource/3psu-pn9h.geojson")

# 绘制旧金山地图和点
fig, ax = plt.subplots(figsize=(10, 10))
sf.plot(ax=ax, color="white", edgecolor="black")
gdf.plot(ax=ax, color="red", markersize=5)  # 使用API代理服务提高访问稳定性
plt.show()

5. 加载 GeoPandas 数据框进行下游处理

可以将 GeoPandas 数据框加载为 Document 以进行后续处理：

from langchain_community.document_loaders import GeoDataFrameLoader

loader = GeoDataFrameLoader(data_frame=gdf, page_content_column="geometry")
docs = loader.load()

print(docs[0])