11、Python地理空间数据处理与数据库应用

Python地理空间数据处理与数据库应用

1. 地理空间距离计算基础

地球表面并非平面，因此在处理地理空间数据时，我们实际计算的是大地测量距离或大圆距离。Python的 pyproj 库可用于执行这些计算，还能实现投影坐标与地理坐标的相互转换。

几何图形（如线或多边形）本质上是由一系列相连的点组成。通过计算多边形中各点间的大地测量距离并求和，就能得到该几何图形的实际长度。下面我们通过一个实际问题来演示如何运用这些知识。

1.1 计算泰缅边境长度

为计算泰缅边境长度，我们使用之前创建的 common-borders/border.shp 文件。该文件包含一个表示两国边境的 LineString 要素。

import os.path
import osgeo.ogr

def getLineSegmentsFromGeometry(geometry):
    segments = []
    if geometry.GetPointCount() > 0:
        segment = []
        for i in range(geometry.GetPointCount()):
            segment.append(geometry.GetPoint_2D(i))
        segments.append(segment)
    for i in range(geometry.GetGeometryCount()):
        sub