gis-distance-tool 0.1.0__tar.gz

gis_distance_tool-0.1.0/PKG-INFO

@@ -0,0 +1,43 @@
+Metadata-Version: 2.4
+Name: gis-distance-tool
+Version: 0.1.0
+Summary: GIS距离分析工具-基于Haversine公式计算大圆距离
+Author-email: belywe-spec <belywe79@gmail.com>
+License: MIT
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: GIS
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+Provides-Extra: numpy
+Requires-Dist: numpy>=1.20.0; extra == "numpy"
+Provides-Extra: scipy
+Requires-Dist: scipy>=1.7.0; extra == "scipy"
+Provides-Extra: full
+Requires-Dist: numpy>=1.20.0; extra == "full"
+Requires-Dist: scipy>=1.7.0; extra == "full"
+Provides-Extra: dev
+Requires-Dist: pytest>=7.0.0; extra == "dev"
+Requires-Dist: black>=22.0.0; extra == "dev"
+Requires-Dist: flake8>=5.0.0; extra == "dev"
+Requires-Dist: mypy>=0.990; extra == "dev"
+
+# GIS Distance Tool
+
+一个基于 Haversine 公式的 GIS 距离计算工具，支持单点和批量计算。
+
+# 安装
+
+# 基础安装（无依赖）
+pip install .
+
+# 安装 NumPy 支持
+pip install .[numpy]
+

gis_distance_tool-0.1.0/README.md

@@ -0,0 +1,12 @@
+# GIS Distance Tool
+
+一个基于 Haversine 公式的 GIS 距离计算工具，支持单点和批量计算。
+
+# 安装
+
+# 基础安装（无依赖）
+pip install .
+
+# 安装 NumPy 支持
+pip install .[numpy]
+

gis_distance_tool-0.1.0/pyproject.toml

@@ -0,0 +1,80 @@
+[build-system]
+# 指定构建工具
+requires = ["setuptools>=61.0", "wheel"]
+build-backend = "setuptools.build_meta"
+
+[project]
+# 项目基本信息
+name = "gis-distance-tool"
+version = "0.1.0"
+description = "GIS距离分析工具-基于Haversine公式计算大圆距离"
+readme = "README.md"
+authors = [
+    {name = "belywe-spec", email = "belywe79@gmail.com"}
+]
+license = {text = "MIT"}
+requires-python = ">=3.8"
+classifiers = [
+    "Development Status :: 3 - Alpha",
+    "Intended Audience :: Science/Research",
+    "License :: OSI Approved :: MIT License",
+    "Programming Language :: Python :: 3",
+    "Programming Language :: Python :: 3.8",
+    "Programming Language :: Python :: 3.9",
+    "Programming Language :: Python :: 3.10",
+    "Programming Language :: Python :: 3.11",
+    "Programming Language :: Python :: 3.12",
+    "Topic :: Scientific/Engineering :: GIS",
+]
+
+# 核心依赖
+dependencies = [
+    # 基础版无需第三方库
+]
+
+# 可选依赖（按功能分组，用户用 pip install .[extra] 安装）
+[project.optional-dependencies]
+numpy = ["numpy>=1.20.0",]
+scipy = ["scipy>=1.7.0",]
+full = [
+    "numpy>=1.20.0",
+    "scipy>=1.7.0",
+]
+dev = [
+    "pytest>=7.0.0",
+    "black>=22.0.0",
+    "flake8>=5.0.0",
+    "mypy>=0.990",
+]
+
+# 命令行入口点
+[project.scripts]
+gis-distance = "gis_tool.cli:main"
+
+# 包目录结构配置
+[tool.setuptools.packages.find]
+where = ["src"]
+include = ["gis_tool*"]
+exclude = ["tests*", "docs*"]
+
+# 工具配置（可选，用于代码质量检查）
+[tool.black]
+line-length = 100
+target-version = ['py38', 'py39', 'py310', 'py311']
+
+[tool.isort]
+profile = "black"
+line_length = 100
+
+[tool.mypy]
+python_version = "3.8"
+warn_return_any = true
+warn_unused_configs = true
+ignore_missing_imports = true
+
+[tool.pytest.ini_options]
+minversion = "7.0"
+testpaths = ["tests"]
+python_files = "test_*.py"
+python_classes = "Test*"
+python_functions = "test_*"

gis_distance_tool-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

gis_distance_tool-0.1.0/src/gis_distance_tool.egg-info/PKG-INFO

@@ -0,0 +1,43 @@
+Metadata-Version: 2.4
+Name: gis-distance-tool
+Version: 0.1.0
+Summary: GIS距离分析工具-基于Haversine公式计算大圆距离
+Author-email: belywe-spec <belywe79@gmail.com>
+License: MIT
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: GIS
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+Provides-Extra: numpy
+Requires-Dist: numpy>=1.20.0; extra == "numpy"
+Provides-Extra: scipy
+Requires-Dist: scipy>=1.7.0; extra == "scipy"
+Provides-Extra: full
+Requires-Dist: numpy>=1.20.0; extra == "full"
+Requires-Dist: scipy>=1.7.0; extra == "full"
+Provides-Extra: dev
+Requires-Dist: pytest>=7.0.0; extra == "dev"
+Requires-Dist: black>=22.0.0; extra == "dev"
+Requires-Dist: flake8>=5.0.0; extra == "dev"
+Requires-Dist: mypy>=0.990; extra == "dev"
+
+# GIS Distance Tool
+
+一个基于 Haversine 公式的 GIS 距离计算工具，支持单点和批量计算。
+
+# 安装
+
+# 基础安装（无依赖）
+pip install .
+
+# 安装 NumPy 支持
+pip install .[numpy]
+

gis_distance_tool-0.1.0/src/gis_distance_tool.egg-info/SOURCES.txt

@@ -0,0 +1,13 @@
+README.md
+pyproject.toml
+src/gis_distance_tool.egg-info/PKG-INFO
+src/gis_distance_tool.egg-info/SOURCES.txt
+src/gis_distance_tool.egg-info/dependency_links.txt
+src/gis_distance_tool.egg-info/entry_points.txt
+src/gis_distance_tool.egg-info/requires.txt
+src/gis_distance_tool.egg-info/top_level.txt
+src/gis_tool/__init__.py
+src/gis_tool/cli.py
+src/gis_tool/core/__init__.py
+src/gis_tool/core/processor.py
+src/gis_tool/core/utils.py

gis_distance_tool-0.1.0/src/gis_distance_tool.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

gis_distance_tool-0.1.0/src/gis_distance_tool.egg-info/entry_points.txt

@@ -0,0 +1,2 @@
+[console_scripts]
+gis-distance = gis_tool.cli:main

gis_distance_tool-0.1.0/src/gis_distance_tool.egg-info/requires.txt

@@ -0,0 +1,16 @@
+
+[dev]
+pytest>=7.0.0
+black>=22.0.0
+flake8>=5.0.0
+mypy>=0.990
+
+[full]
+numpy>=1.20.0
+scipy>=1.7.0
+
+[numpy]
+numpy>=1.20.0
+
+[scipy]
+scipy>=1.7.0

gis_distance_tool-0.1.0/src/gis_distance_tool.egg-info/top_level.txt

@@ -0,0 +1 @@
+gis_tool

gis_distance_tool-0.1.0/src/gis_tool/cli.py

@@ -0,0 +1,142 @@
+#!/usr/bin/env python3
+"""
+GIS 距离分析工具 - 命令行入口
+使用示例:
+    # 单点距离计算
+    python cli.py single --lat1 31.23 --lon1 121.47 --lat2 39.90 --lon2 116.40
+    
+    # 批量距离计算（需安装 NumPy）
+    python cli.py batch --coords1 31.23 121.47 39.90 116.40 --coords2 31.23 121.47 35.68 139.76
+"""
+
+import argparse
+import sys
+from gis_tool.core.processor import haversine_distance, haversine_distance_batch
+
+
+def parse_coordinates(flat_list):
+    """
+    将扁平的坐标列表转为 [(lat, lon), ...] 格式
+    例如: [31.23, 121.47, 39.90, 116.40] -> [(31.23, 121.47), (39.90, 116.40)]
+    """
+    if len(flat_list) % 2 != 0:
+        raise ValueError("坐标个数必须为偶数（每两个数为一组经纬度）")
+    
+    coords = []
+    for i in range(0, len(flat_list), 2):
+        coords.append((flat_list[i], flat_list[i+1]))
+    return coords
+
+
+def main():
+    # 创建主解析器
+    parser = argparse.ArgumentParser(
+        description="GIS 距离分析工具 - 基于 Haversine 公式计算大圆距离",
+        formatter_class=argparse.RawDescriptionHelpFormatter,
+        epilog="""
+示例:
+  # 单点距离（无依赖）
+  python cli.py single --lat1 31.2304 --lon1 121.4737 --lat2 39.9042 --lon2 116.4074
+  
+  # 单点距离（输出米）
+  python cli.py single --lat1 31.2304 --lon1 121.4737 --lat2 39.9042 --lon2 116.4074 --unit m
+  
+  # 批量距离（需安装 NumPy）
+  python cli.py batch --coords1 31.23 121.47 39.90 116.40 --coords2 31.23 121.47 35.68 139.76
+        """
+    )
+    
+    # 添加子命令
+    subparsers = parser.add_subparsers(dest="command", required=True, help="计算模式")
+    
+    # ---------- 子命令1: single（单点距离） ----------
+    parser_single = subparsers.add_parser(
+        "single",
+        help="计算两个单点之间的大圆距离（无第三方依赖）"
+    )
+    parser_single.add_argument("--lat1", type=float, required=True, help="起点纬度 (-90 ~ 90)")
+    parser_single.add_argument("--lon1", type=float, required=True, help="起点经度 (-180 ~ 180)")
+    parser_single.add_argument("--lat2", type=float, required=True, help="终点纬度 (-90 ~ 90)")
+    parser_single.add_argument("--lon2", type=float, required=True, help="终点经度 (-180 ~ 180)")
+    parser_single.add_argument(
+        "--unit",
+        type=str,
+        default="km",
+        choices=["km", "m"],
+        help="距离单位，默认 km"
+    )
+    
+    # ---------- 子命令2: batch（批量距离） ----------
+    parser_batch = subparsers.add_parser(
+        "batch",
+        help="批量计算两组坐标间的距离（需安装 NumPy）"
+    )
+    parser_batch.add_argument(
+        "--coords1",
+        type=float,
+        nargs="+",
+        required=True,
+        help="第一组坐标，格式: lat1 lon1 lat2 lon2 ... (必须为偶数个)"
+    )
+    parser_batch.add_argument(
+        "--coords2",
+        type=float,
+        nargs="+",
+        required=True,
+        help="第二组坐标，格式: lat1 lon1 lat2 lon2 ... (必须为偶数个)"
+    )
+    parser_batch.add_argument(
+        "--unit",
+        type=str,
+        default="km",
+        choices=["km", "m"],
+        help="距离单位，默认 km"
+    )
+    
+    # 解析参数
+    args = parser.parse_args()
+    
+    try:
+        if args.command == "single":
+            # 调用基础版距离计算
+            result = haversine_distance(
+                args.lat1, args.lon1,
+                args.lat2, args.lon2,
+                args.unit
+            )
+            print(f"✅ 两点之间的大圆距离: {result:.4f} {args.unit}")
+            
+        elif args.command == "batch":
+            # 解析坐标列表
+            coords1 = parse_coordinates(args.coords1)
+            coords2 = parse_coordinates(args.coords2)
+            
+            # 调用 NumPy 版批量计算
+            result = haversine_distance_batch(coords1, coords2, args.unit)
+            
+            # 打印结果（格式化输出）
+            import numpy as np
+            if isinstance(result, np.ndarray):
+                if result.ndim == 1:
+                    print(f"✅ 逐点距离向量 (长度={len(result)}):")
+                    for i, d in enumerate(result):
+                        print(f"   点{i}: {d:.4f} {args.unit}")
+                else:
+                    print(f"✅ 距离矩阵 (形状={result.shape}):")
+                    print(result)
+            else:
+                print(f"✅ 距离计算结果: {result}")
+            
+    except ImportError as e:
+        print(f"❌ 缺少依赖库: {e}", file=sys.stderr)
+        sys.exit(1)
+    except ValueError as e:
+        print(f"❌ 输入错误: {e}", file=sys.stderr)
+        sys.exit(1)
+    except Exception as e:
+        print(f"❌ 运行错误: {e}", file=sys.stderr)
+        sys.exit(1)
+
+
+if __name__ == "__main__":
+    main()

gis_distance_tool-0.1.0/src/gis_tool/core/processor.py

@@ -0,0 +1,101 @@
+import math
+from gis_tool.core.utils import validate_coordinates
+
+# 基础版：纯 Python 实现 Haversine 距离计算
+def haversine_distance(lat1: float, lon1: float, lat2: float, lon2: float, unit: str = "km") -> float:
+    """
+    使用 Haversine 公式计算地球表面两点间的大圆距离
+    
+    参数:
+        lat1, lon1: 起点坐标（度）
+        lat2, lon2: 终点坐标（度）
+        unit: 输出单位，'km' 或 'm'
+    
+    返回:
+        距离值（浮点数）
+    """
+    # 1. 校验坐标合法性
+    validate_coordinates(lat1, lon1)
+    validate_coordinates(lat2, lon2)
+    
+    # 2. 地球平均半径（公里）
+    R = 6371.0
+    
+    # 3. 角度转弧度
+    phi1 = math.radians(lat1)
+    phi2 = math.radians(lat2)
+    delta_phi = math.radians(lat2 - lat1)
+    delta_lambda = math.radians(lon2 - lon1)
+    
+    # 4. Haversine 公式
+    a = math.sin(delta_phi / 2) ** 2 + \
+        math.cos(phi1) * math.cos(phi2) * math.sin(delta_lambda / 2) ** 2
+    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))
+    distance_km = R * c
+    
+    # 5. 按单位返回
+    return distance_km if unit == "km" else distance_km * 1000.0
+
+
+# 进阶版：使用 NumPy 实现批量距离计算
+def haversine_distance_batch(coords1, coords2, unit: str = "km"):
+    """
+    使用 NumPy 批量计算两组坐标间的距离
+    
+    参数:
+        coords1: 第一组坐标，形状为 (N, 2)，每行 [lat, lon]
+        coords2: 第二组坐标，形状为 (M, 2)，每行 [lat, lon]
+        unit: 输出单位，'km' 或 'm'
+    
+    返回:
+        如果 N == M: 返回长度为 N 的距离向量
+        如果 N != M: 返回 N x M 的距离矩阵
+    
+    注意:
+        此函数需要安装 NumPy，否则会抛出 ImportError
+    """
+    try:
+        import numpy as np
+    except ImportError:
+        raise ImportError(
+            "NumPy 未安装，请运行 'pip install numpy' 后使用此功能"
+        )
+    
+    # 转换为 numpy 数组
+    arr1 = np.asarray(coords1, dtype=np.float64)
+    arr2 = np.asarray(coords2, dtype=np.float64)
+    
+    # 检查形状
+    if arr1.shape[1] != 2 or arr2.shape[1] != 2:
+        raise ValueError("坐标数组必须为 (N, 2) 形状，每行 [lat, lon]")
+    
+    # 校验所有坐标
+    for lat, lon in arr1:
+        validate_coordinates(lat, lon)
+    for lat, lon in arr2:
+        validate_coordinates(lat, lon)
+    
+    # 角度转弧度
+    lat1 = np.radians(arr1[:, 0])
+    lon1 = np.radians(arr1[:, 1])
+    lat2 = np.radians(arr2[:, 0])
+    lon2 = np.radians(arr2[:, 1])
+    
+    # 根据坐标数量决定计算方式
+    if len(arr1) == len(arr2):
+        # 逐点计算（一对一）
+        dlat = lat2 - lat1
+        dlon = lon2 - lon1
+        a = np.sin(dlat / 2) ** 2 + np.cos(lat1) * np.cos(lat2) * np.sin(dlon / 2) ** 2
+        c = 2 * np.arcsin(np.sqrt(np.clip(a, 0, 1)))  # 防止浮点误差导致 >1
+        distance_km = 6371.0 * c
+    else:
+        # 计算距离矩阵（广播机制）
+        # 扩展维度: (N, 1) 和 (1, M) 进行广播
+        dlat = lat2[np.newaxis, :] - lat1[:, np.newaxis]
+        dlon = lon2[np.newaxis, :] - lon1[:, np.newaxis]
+        a = np.sin(dlat / 2) ** 2 + np.cos(lat1[:, np.newaxis]) * np.cos(lat2[np.newaxis, :]) * np.sin(dlon / 2) ** 2
+        c = 2 * np.arcsin(np.sqrt(np.clip(a, 0, 1)))
+        distance_km = 6371.0 * c
+    
+    return distance_km if unit == "km" else distance_km * 1000.0

gis_distance_tool-0.1.0/src/gis_tool/core/utils.py

@@ -0,0 +1,8 @@
+def validate_coordinates(lat: float, lon: float) -> bool:
+    
+    # 校验经纬度是否在合法范围内
+    if not (-90 <= lat <= 90):
+        raise ValueError(f"纬度 {lat} 超出有效范围 [-90, 90]")
+    if not (-180 <= lon <= 180):
+        raise ValueError(f"经度 {lon} 超出有效范围 [-180, 180]")
+    return True