1. 項目概述當AI代理遇見地理空間最近在AI和地理信息科學(xué)GIS的交叉領(lǐng)域一個名為“OpenEarthAgent”的項目引起了我的注意。簡單來說它試圖解決一個核心痛點如何讓AI智能體Agent像人類專家一樣去理解和操作復(fù)雜的地理空間數(shù)據(jù)與工具。我們常說的AI代理比如基于大語言模型LLM的助手已經(jīng)能很好地處理文本、代碼甚至圖像但一旦涉及到帶有坐標、投影、拓撲關(guān)系的地理數(shù)據(jù)以及需要調(diào)用ArcGIS、QGIS插件或?qū)I(yè)遙感分析庫時它們往往就“抓瞎”了。OpenEarthAgent正是瞄準了這個空白它不是一個單一的工具而是一個統(tǒng)一的框架旨在為地理空間任務(wù)構(gòu)建“工具增強型”的AI代理。想象一下你是一個城市規(guī)劃師想分析某個區(qū)域過去十年的土地利用變化。傳統(tǒng)流程是打開專業(yè)GIS軟件加載多期遙感影像進行預(yù)處理、分類、變化檢測最后出圖制表。每一步都需要專業(yè)知識。而OpenEarthAgent框架的目標是讓你能用自然語言告訴AI代理“幫我分析一下A市B區(qū)2013年到2023年的土地利用變化并生成一份簡報?！?這個代理就能自動理解你的意圖規(guī)劃任務(wù)步驟并調(diào)用背后集成的各種地理空間工具如GDAL進行數(shù)據(jù)讀取scikit-learn進行分類GeoPandas進行空間運算來完成任務(wù)。它把分散的、專業(yè)門檻高的工具鏈通過一個智能的“大腦”統(tǒng)一調(diào)度起來這就是“工具增強”和“統(tǒng)一框架”的核心價值。這個框架適合誰呢我認為有三類人最需要關(guān)注一是地理信息領(lǐng)域的科研人員和工程師他們可以基于此框架快速構(gòu)建面向特定場景如災(zāi)害評估、環(huán)境監(jiān)測的自動化分析流水線二是希望將AI能力融入現(xiàn)有地理信息產(chǎn)品的開發(fā)者框架提供了標準的集成接口三是對空間智能Spatial AI感興趣的研究者這提供了一個絕佳的實驗平臺。接下來我將深入拆解這個框架的設(shè)計思路、核心模塊并分享如何從零開始構(gòu)建一個簡易版的地理空間代理以及其中會遇到的各種“坑”。2. 框架核心設(shè)計思路與架構(gòu)拆解要理解OpenEarthAgent我們不能只把它看作一堆代碼的集合而應(yīng)該從它要解決的問題和設(shè)計哲學(xué)入手。其核心思路是構(gòu)建一個感知-規(guī)劃-執(zhí)行-學(xué)習(xí)的閉環(huán)智能體并專門針對地理空間數(shù)據(jù)的特性進行優(yōu)化。2.1 為什么需要“工具增強”而非“端到端”模型這是第一個關(guān)鍵設(shè)計決策。目前雖然有像Segment Anything Model (SAM)這樣的通用視覺模型但要處理多樣、復(fù)雜且專業(yè)的地理空間分析任務(wù)如計算NDVI、進行水文分析、處理投影變換訓(xùn)練一個“全能”的端到端模型成本極高且難以保證精度和可解釋性。相反地理信息領(lǐng)域經(jīng)過數(shù)十年發(fā)展已經(jīng)積累了如GDAL、PROJ、PostGIS、WhiteboxTools等大量成熟、穩(wěn)定、高效的開源工具庫。OpenEarthAgent框架采用了“工具增強”Tool-Augmented的策略即讓AI代理作為協(xié)調(diào)者和規(guī)劃者而具體的“重活累活”交給這些經(jīng)過驗證的專業(yè)工具去執(zhí)行。這帶來了幾個顯著優(yōu)勢可靠性直接利用業(yè)界標準工具結(jié)果可信度高?？蓴U展性新增一個分析能力只需為代理“裝配”一個新的工具函數(shù)無需重新訓(xùn)練核心模型?？山忉屝源淼拿恳徊讲僮髡{(diào)用哪個工具、輸入什么參數(shù)都是清晰可追溯的符合地理分析對過程嚴謹性的要求。資源效率避免為每一個細分任務(wù)都訓(xùn)練大模型節(jié)省計算資源。2.2 統(tǒng)一框架的層次化架構(gòu)根據(jù)我對類似系統(tǒng)如LangChain的Agent架構(gòu)和地理信息處理流程的理解一個合理的OpenEarthAgent框架應(yīng)包含以下層次交互層Interface Layer負責(zé)與用戶或上游系統(tǒng)對接。支持多種輸入方式如自然語言指令、圖形界面操作或API調(diào)用。這一層的關(guān)鍵是將模糊的用戶需求轉(zhuǎn)化為結(jié)構(gòu)化的“任務(wù)意圖”。例如將“找出所有受洪水影響的居民區(qū)”解析為包含關(guān)鍵要素災(zāi)害類型洪水分析對象居民區(qū)操作空間篩選的機器可讀格式。代理核心層Agent Core Layer這是框架的“大腦”通常由一個或一組大語言模型驅(qū)動。它的核心功能是任務(wù)規(guī)劃與工具調(diào)度。任務(wù)分解將復(fù)雜的宏觀任務(wù)如“評估風(fēng)災(zāi)損失”分解為一系列有序的原子操作如“獲取風(fēng)場數(shù)據(jù)”、“獲取人口分布數(shù)據(jù)”、“進行疊加分析”、“計算經(jīng)濟損失”。工具匹配為每一個原子操作從工具庫中選擇最合適的工具。這需要模型對工具的功能描述通常用結(jié)構(gòu)化文本定義有深刻理解。參數(shù)推理根據(jù)任務(wù)上下文自動推斷工具調(diào)用所需的參數(shù)。例如當任務(wù)是需要計算NDVI時代理應(yīng)能自動識別輸入數(shù)據(jù)中的近紅外波段和紅光波段索引。工具抽象層Tool Abstraction Layer這是框架的“手”和“腳”。它將底層的、異構(gòu)的地理空間工具可能是命令行工具、Python庫函數(shù)、REST API進行統(tǒng)一封裝向上提供標準化的調(diào)用接口通常是函數(shù)。一個工具定義通常包括工具名稱、功能描述、所需參數(shù)列表名稱、類型、描述、返回結(jié)果類型。例如一個buffer_analysis工具會被描述為“對輸入的矢量要素進行緩沖區(qū)分析。參數(shù)input_geojsonGeoJSON格式的要素distance緩沖距離單位米。返回緩沖后的GeoJSON。”地理空間運行時Geospatial Runtime這是所有工具實際執(zhí)行的環(huán)境。它需要管理地理數(shù)據(jù)的I/O、坐標參考系統(tǒng)CRS的統(tǒng)一與轉(zhuǎn)換、內(nèi)存與計算資源。一個健壯的運行時必須能優(yōu)雅地處理地理空間數(shù)據(jù)特有的問題比如不同數(shù)據(jù)源之間的CRS不一致、大規(guī)模柵格數(shù)據(jù)的分塊處理等。記憶與反饋層Memory Feedback Layer使代理具備持續(xù)學(xué)習(xí)能力。短期記憶保存當前會話的上下文確保在多輪對話中理解用戶意圖。長期記憶則可以存儲歷史任務(wù)的成功模式或失敗教訓(xùn)用于優(yōu)化未來的規(guī)劃決策。用戶對結(jié)果的反饋“這個范圍不對再大一點”也可以被用來調(diào)整代理的行為。注意這種分層架構(gòu)的關(guān)鍵在于“松耦合”。代理核心層不需要知道GDAL庫是用C寫的它只需要調(diào)用read_raster這個工具函數(shù)。這極大提高了框架的靈活性和可維護性。3. 核心模塊深度解析與實操要點理解了宏觀架構(gòu)我們深入到幾個核心模塊看看它們具體如何工作以及在實現(xiàn)時需要注意什么。3.1 工具的定義與封裝讓AI“懂”地理工具工具封裝的質(zhì)量直接決定了代理的能力上限。封裝不僅僅是寫一個Python函數(shù)包裝器那么簡單。一個完整的工具定義示例from typing import Dict, Any import geopandas as gpd from shapely.geometry import shape def spatial_join( target_features: Dict, # GeoJSON格式 join_features: Dict, # GeoJSON格式 op: str intersects # 空間關(guān)系如 ‘within’ ‘contains’ ) - Dict: 對兩個矢量數(shù)據(jù)集進行空間連接Spatial Join。 將join_features的屬性連接到與其有空間關(guān)系的target_features上。 參數(shù): target_features: 目標要素GeoJSON格式的字典。 join_features: 連接要素GeoJSON格式的字典。 op: 空間關(guān)系操作默認為‘intersects’相交。 返回: 空間連接后的GeoJSON格式字典。 # 1. 將GeoJSON轉(zhuǎn)換為GeoDataFrame gdf_target gpd.GeoDataFrame.from_features(target_features[features]) gdf_join gpd.GeoDataFrame.from_features(join_features[features]) # 2. 確保坐標系一致關(guān)鍵步驟 if gdf_target.crs ! gdf_join.crs: gdf_join gdf_join.to_crs(gdf_target.crs) # 3. 執(zhí)行空間連接 gdf_result gdf_target.sjoin(gdf_join, howleft, predicateop) # 4. 轉(zhuǎn)換回GeoJSON result_geojson gdf_result.to_json() return result_geojson # 提供給Agent框架的工具描述 TOOL_DESCRIPTION { name: spatial_join, description: Performs a spatial join between two vector datasets. It transfers attributes from the join features to the target features based on a spatial relationship (e.g., intersects, within)., parameters: { target_features: {type: object, description: GeoJSON object for target features.}, join_features: {type: object, description: GeoJSON object for join features.}, op: {type: string, description: Spatial predicate: intersects, within, contains, etc., default: intersects} }, returns: {type: object, description: GeoJSON object of the joined features.} }實操要點與避坑指南描述要精準且面向任務(wù)工具描述description是AI理解工具的“說明書”。避免使用“處理空間數(shù)據(jù)”這種模糊描述而要用“計算多邊形要素的質(zhì)心坐標”這樣具體、可操作的句子。描述應(yīng)包含輸入、輸出和核心功能。參數(shù)類型和默認值至關(guān)重要LLM對強類型和默認值很敏感。明確參數(shù)類型string,number,object能減少調(diào)用錯誤。合理的默認值可以簡化用戶指令用戶說“緩沖一下這個區(qū)域”代理可以自動使用默認距離。內(nèi)部必須處理CRS這是地理空間編程中最常見的坑。任何涉及多個數(shù)據(jù)源的工具在操作前必須檢查并統(tǒng)一坐標系。如上例所示忽略CRS轉(zhuǎn)換會導(dǎo)致空間分析結(jié)果完全錯誤。一個最佳實踐是在框架層面約定一個內(nèi)部統(tǒng)一坐標系如WGS84 Web墨卡托EPSG:3857或WGS84經(jīng)緯度EPSG:4326所有工具在輸入輸出時都明確遵循此約定或在工具內(nèi)部進行強制轉(zhuǎn)換。錯誤處理與友好反饋工具函數(shù)內(nèi)部必須有完善的try-except并將地理庫可能拋出的專業(yè)錯誤如CRSError,TopologicalError轉(zhuǎn)化為代理和用戶能理解的友好信息例如“無法執(zhí)行疊加分析可能是因為兩個數(shù)據(jù)層的坐標系不一致”。3.2 任務(wù)規(guī)劃與工具鏈編排AI的“思考”過程代理核心如何將“分析城市公園的服務(wù)盲區(qū)”變成一序列工具調(diào)用這依賴于提示工程Prompt Engineering和規(guī)劃算法。一個簡化的規(guī)劃流程可能是指令解析用戶輸入 - LLM - 結(jié)構(gòu)化任務(wù)描述。提示詞示例“請將以下用戶指令解析為包含‘目標’、‘約束條件’和‘期望輸出’的JSON格式。指令{用戶指令}”任務(wù)分解結(jié)構(gòu)化任務(wù) - LLM - 步驟列表。提示詞示例“給定任務(wù)目標‘{目標}’請將其分解為不超過5個順序執(zhí)行的子步驟。每個步驟應(yīng)是一個可執(zhí)行的動作描述例如‘加載道路網(wǎng)絡(luò)數(shù)據(jù)’、‘計算公園可達性’?！惫ぞ咂ヅ渑c參數(shù)填充對每一個子步驟LLM需要從工具庫中選擇工具并填充參數(shù)。提示詞示例“現(xiàn)有工具庫{工具描述列表}。當前步驟是‘{步驟描述}’。請選擇最合適的工具名稱并基于常識推斷其所需的參數(shù)值。如果信息不足請指出需要用戶澄清什么。”我的實操心得思維鏈Chain-of-Thought提示至關(guān)重要直接讓LLM輸出最終工具調(diào)用序列成功率低。更好的方法是引導(dǎo)它“一步一步思考”。例如在提示詞中要求“首先要計算服務(wù)盲區(qū)我需要知道公園位置和居民區(qū)位置。然后我需要計算每個居民區(qū)到最近公園的距離。最后根據(jù)一個閾值比如500米篩選出距離過遠的居民區(qū)?！?LLM輸出這樣的推理過程后再將其映射到工具上會準確得多。給LLM提供“范例”Few-Shot Learning在提示詞中提供2-3個從自然語言到工具調(diào)用鏈的完整示例能極大提升模型在陌生任務(wù)上的表現(xiàn)。這相當于給AI看了幾個“標準作業(yè)流程”。設(shè)計“驗證-重試”機制代理調(diào)用工具后應(yīng)對結(jié)果進行簡單驗證如檢查返回數(shù)據(jù)是否為空、幾何是否有效。如果失敗應(yīng)將錯誤信息反饋給LLM讓它重新規(guī)劃或調(diào)整參數(shù)。這構(gòu)成了一個簡單的自我修正循環(huán)。3.3 地理空間數(shù)據(jù)與模型的中間表示AI模型尤其是LLM通常處理文本而地理數(shù)據(jù)是復(fù)雜的二進制或結(jié)構(gòu)化數(shù)值。如何讓兩者溝通需要一個高效的中間表示。矢量數(shù)據(jù)GeoJSON是目前最通用的文本化表示格式。它用JSON描述地理要素可讀性好與Web技術(shù)棧兼容性極佳。代理內(nèi)部流轉(zhuǎn)可以使用GeoJSON在調(diào)用底層工具如Shapely, GeoPandas時再臨時轉(zhuǎn)換為幾何對象。對于特別大的矢量數(shù)據(jù)可以傳遞數(shù)據(jù)索引或URI而非數(shù)據(jù)本身由專門的數(shù)據(jù)加載工具按需讀取。柵格數(shù)據(jù)直接傳遞TIFF/IMG文件的二進制數(shù)據(jù)不現(xiàn)實。通常采用以下策略元數(shù)據(jù)摘要傳遞柵格文件的元數(shù)據(jù)范圍、分辨率、波段數(shù)、CRS給代理用于規(guī)劃?？s略圖或統(tǒng)計信息對于需要視覺判斷或概要分析的任務(wù)可以生成一個小尺寸的PNG預(yù)覽圖或各波段的統(tǒng)計值最小值、最大值、均值、標準差供LLM參考。切片或區(qū)域查詢當代理確定需要處理某塊具體區(qū)域時再調(diào)用工具讀取該區(qū)域的像素數(shù)據(jù)。任務(wù)上下文表示除了數(shù)據(jù)本身任務(wù)執(zhí)行中的中間狀態(tài)如“已加載A市邊界”、“已計算NDVI指數(shù)圖”也需要一種方式在代理的“記憶”中留存。這可以通過維護一個鍵值對狀態(tài)字典來實現(xiàn)其中值可以是數(shù)據(jù)的引用URI、結(jié)果的簡短文本描述或關(guān)鍵數(shù)值。4. 從零構(gòu)建一個簡易地理空間代理實戰(zhàn)理論說了這么多我們來動手實現(xiàn)一個簡化版的“微縮OpenEarthAgent”。這個實戰(zhàn)將聚焦于一個具體任務(wù)“給定一個點的坐標找出其周圍10公里內(nèi)所有的醫(yī)院并計算到每個醫(yī)院的直線距離。”4.1 環(huán)境準備與工具庫搭建我們選擇Python作為實現(xiàn)語言因為它在地理空間和數(shù)據(jù)科學(xué)領(lǐng)域有最豐富的生態(tài)。1. 核心依賴安裝# 地理數(shù)據(jù)處理核心庫 pip install geopandas shapely pyproj folium # 用于在線獲取數(shù)據(jù)的庫演示用 pip install requests # LLM交互這里以O(shè)penAI API為例也可用本地模型如Ollama pip install openai # 代理框架基礎(chǔ)我們這里簡化不直接用LangChain但借鑒其思想 # pip install langchain langchain-openai2. 定義我們的微型工具庫創(chuàng)建一個geo_tools.py文件。# geo_tools.py import math import requests import geopandas as gpd from shapely.geometry import Point, shape from typing import List, Dict, Any def create_buffer(center_point: Dict, radius_km: float) - Dict: 根據(jù)中心點和半徑創(chuàng)建緩沖區(qū)圓形。 參數(shù): center_point: 包含‘lon’和‘lat’鍵的字典單位度。 radius_km: 緩沖區(qū)半徑單位公里。 返回: 緩沖區(qū)多邊形的GeoJSON字典。 lon, lat center_point[lon], center_point[lat] # 簡化的地理緩沖區(qū)計算更精確需使用投影 # 將公里轉(zhuǎn)換為近似的度數(shù)1度約111公里 radius_deg radius_km / 111.0 center Point(lon, lat) buffer_polygon center.buffer(radius_deg) return gpd.GeoSeries([buffer_polygon]).__geo_interface__ def fetch_pois_within_bounds(bounds: Dict, poi_type: str hospital) - Dict: 模擬從外部API獲取興趣點POI數(shù)據(jù)。 實際中可替換為Overpass API、Nominatim或商業(yè)POI數(shù)據(jù)源。 參數(shù): bounds: 包含‘west’, ‘south’, ‘east’, ‘north’鍵的字典表示地理范圍。 poi_type: 興趣點類型。 返回: 模擬的醫(yī)院POI的GeoJSON字典。 # 這里是模擬數(shù)據(jù)真實場景需要調(diào)用真實API。 # 例如使用OpenStreetMap的Overpass API # url fhttps://overpass-api.de/api/interpreter?data[out:json];node[amenity{poi_type}]({bounds[south]},{bounds[west]},{bounds[north]},{bounds[east]});out body; # response requests.get(url).json() # 為演示我們生成一些模擬點 import random west, south, east, north bounds[west], bounds[south], bounds[east], bounds[north] features [] for i in range(5): sim_lon west random.random() * (east - west) sim_lat south random.random() * (north - south) point Point(sim_lon, sim_lat) feature { type: Feature, geometry: point.__geo_interface__, properties: {name: f{poi_type.capitalize()} {i1}, id: i} } features.append(feature) return {type: FeatureCollection, features: features} def calculate_distance(point1: Dict, point2: Dict) - float: 計算兩點之間的近似大圓距離Haversine公式。 參數(shù): point1/point2: 包含‘lon’和‘lat’鍵的字典。 返回: 距離單位公里。 lon1, lat1 math.radians(point1[lon]), math.radians(point1[lat]) lon2, lat2 math.radians(point2[lon]), math.radians(point2[lat]) dlon, dlat lon2 - lon1, lat2 - lat1 a math.sin(dlat/2)**2 math.cos(lat1) * math.cos(lat2) * math.sin(dlon/2)**2 c 2 * math.asin(math.sqrt(a)) radius_earth_km 6371.0 return c * radius_earth_km def spatial_analysis_main(center_lon: float, center_lat: float, radius_km: float) - List[Dict]: 主分析函數(shù)串聯(lián)以上工具完成完整任務(wù)。 這是我們的“手工編排”版本后續(xù)將由Agent自動完成。 # 1. 創(chuàng)建緩沖區(qū) center {lon: center_lon, lat: center_lat} buffer_geojson create_buffer(center, radius_km) # 從GeoJSON中提取邊界框用于查詢POI buffer_geom shape(buffer_geojson[features][0][geometry]) bounds buffer_geom.bounds # (minx, miny, maxx, maxy) bounds_dict {west: bounds[0], south: bounds[1], east: bounds[2], north: bounds[3]} # 2. 獲取緩沖區(qū)內(nèi)的醫(yī)院 hospitals_geojson fetch_pois_within_bounds(bounds_dict, hospital) # 3. 計算每個醫(yī)院到中心點的距離 results [] for feature in hospitals_geojson[features]: hosp_coords feature[geometry][coordinates] hosp_point {lon: hosp_coords[0], lat: hosp_coords[1]} distance calculate_distance(center, hosp_point) results.append({ name: feature[properties][name], distance_km: round(distance, 2), coordinates: hosp_coords }) # 按距離排序 results.sort(keylambda x: x[distance_km]) return results # 工具描述列表用于提供給LLM TOOLS_FOR_AGENT [ { name: create_buffer, description: Creates a circular buffer polygon around a given center point with a specified radius in kilometers., parameters: { center_point: {type: object, description: Dict with keys lon and lat representing longitude and latitude in degrees.}, radius_km: {type: number, description: Radius of the buffer in kilometers.} } }, { name: fetch_pois_within_bounds, description: Fetches points of interest (POIs) of a specified type (e.g., hospital) within a given geographic bounding box., parameters: { bounds: {type: object, description: Dict with keys west, south, east, north defining the bounding box.}, poi_type: {type: string, description: Type of POI to fetch, e.g., hospital, school., default: hospital} } }, { name: calculate_distance, description: Calculates the great-circle distance between two geographic points (in degrees) using the Haversine formula, returning distance in kilometers., parameters: { point1: {type: object, description: First point with lon and lat.}, point2: {type: object, description: Second point with lon and lat.} } } ]4.2 構(gòu)建代理核心與任務(wù)執(zhí)行引擎接下來我們創(chuàng)建一個簡單的代理它使用LLM來規(guī)劃任務(wù)并調(diào)用我們定義的工具。這里我們使用OpenAI的GPT-4 API作為“大腦”。# agent_core.py import openai import json from geo_tools import create_buffer, fetch_pois_within_bounds, calculate_distance, TOOLS_FOR_AGENT class SimpleGeoAgent: def __init__(self, api_key): openai.api_key api_key self.client openai.OpenAI() self.tools TOOLS_FOR_AGENT # 一個簡單的上下文記憶存儲中間結(jié)果 self.context {} def _call_llm(self, prompt, system_messageYou are a helpful geospatial assistant.): 調(diào)用LLM的通用函數(shù) try: response self.client.chat.completions.create( modelgpt-4, # 或 gpt-3.5-turbo messages[ {role: system, content: system_message}, {role: user, content: prompt} ], temperature0.1, # 低溫度保證輸出穩(wěn)定性 ) return response.choices[0].message.content except Exception as e: return fError calling LLM: {e} def _parse_llm_plan(self, plan_text: str): 解析LLM返回的規(guī)劃文本。 期望格式一個JSON字符串包含步驟列表每個步驟有‘tool’和‘inputs’。 try: # 嘗試從文本中提取JSON部分 import re json_match re.search(rjson\n(.*?)\n, plan_text, re.DOTALL) if json_match: plan_text json_match.group(1) plan json.loads(plan_text) return plan except json.JSONDecodeError: # 如果解析失敗返回一個簡單回退計劃 print(fFailed to parse LLM plan. Raw output:\n{plan_text}) return {steps: []} def plan_task(self, user_query: str): 讓LLM根據(jù)用戶查詢和可用工具制定計劃 tools_description json.dumps(self.tools, indent2) system_msg You are an expert in geospatial analysis. Your job is to break down a users request into a sequence of tool calls. Available tools are described below. prompt f User Request: {user_query} Available Tools (in JSON format): {tools_description} Based on the users request and the available tools, generate a step-by-step execution plan. The plan should be a JSON object with a key steps, which is a list. Each step in the list should be an object with: - step_number: integer - tool: the exact name of the tool to use (must match one of the available tool names) - inputs: an object containing the input parameters for the tool. You must infer reasonable values from the user request or use defaults. If a value cannot be inferred, set it to null. Output only the JSON plan, no other text. Example plan for find hospitals within 5km of point (116.4, 39.9): {{ steps: [ {{ step_number: 1, tool: create_buffer, inputs: {{ center_point: {{lon: 116.4, lat: 39.9}}, radius_km: 5 }} }}, {{ step_number: 2, tool: fetch_pois_within_bounds, inputs: {{ bounds: This should be the OUTPUT from step 1s buffer geometry bounds., poi_type: hospital }} }} ] }} llm_response self._call_llm(prompt, system_msg) return self._parse_llm_plan(llm_response) def execute_plan(self, plan): 執(zhí)行LLM生成的計劃 results {} for step in plan.get(steps, []): tool_name step[tool] inputs step[inputs] print(fExecuting Step {step[step_number]}: {tool_name} with inputs {inputs}) # 處理動態(tài)輸入將上一步的結(jié)果作為下一步的輸入 # 這里簡單實現(xiàn)如果輸入值是字符串且以step_開頭則從results中獲取 resolved_inputs {} for key, value in inputs.items(): if isinstance(value, str) and value.startswith(step_): prev_step_num int(value.split(_)[1]) # 這里需要更復(fù)雜的邏輯來映射上一步的哪個輸出作為當前輸入 # 為簡化我們假設(shè)每個工具只有一個主要輸出存儲在results中 if prev_step_num in results: resolved_inputs[key] results[prev_step_num] else: resolved_inputs[key] value else: resolved_inputs[key] value # 調(diào)用對應(yīng)的工具函數(shù) try: if tool_name create_buffer: output create_buffer(**resolved_inputs) elif tool_name fetch_pois_within_bounds: output fetch_pois_within_bounds(**resolved_inputs) elif tool_name calculate_distance: output calculate_distance(**resolved_inputs) else: output fError: Unknown tool {tool_name} results[step[step_number]] output print(f Result: {str(output)[:100]}...) # 打印前100字符 except Exception as e: error_msg fTool execution error: {e} results[step[step_number]] error_msg print(f Error: {error_msg}) break # 或根據(jù)策略決定是否繼續(xù) return results def run(self, user_query): 運行代理規(guī)劃并執(zhí)行 print(fProcessing query: {user_query}) plan self.plan_task(user_query) print(fGenerated Plan:\n{json.dumps(plan, indent2)}) final_results self.execute_plan(plan) return final_results # 使用示例 if __name__ __main__: # 注意你需要設(shè)置自己的OPENAI_API_KEY環(huán)境變量 import os api_key os.getenv(OPENAI_API_KEY) if not api_key: print(Please set OPENAI_API_KEY environment variable.) else: agent SimpleGeoAgent(api_key) # 測試查詢 query Find all hospitals within 10 kilometers of longitude 116.4074 and latitude 39.9042, and tell me their distances. results agent.run(query) # 后續(xù)可以添加一個“結(jié)果總結(jié)”步驟讓LLM解析final_results并生成自然語言回答。4.3 執(zhí)行流程與結(jié)果解析運行上述代碼你會看到類似以下的輸出具體內(nèi)容因LLM輸出而異Processing query: Find all hospitals within 10 kilometers of longitude 116.4074 and latitude 39.9042, and tell me their distances. Generated Plan: { steps: [ { step_number: 1, tool: create_buffer, inputs: { center_point: {lon: 116.4074, lat: 39.9042}, radius_km: 10 } }, { step_number: 2, tool: fetch_pois_within_bounds, inputs: { bounds: step_1, // LLM可能會聰明地引用上一步的結(jié)果 poi_type: hospital } }, { step_number: 3, tool: calculate_distance, inputs: { point1: {lon: 116.4074, lat: 39.9042}, point2: step_2 // 這里需要更精細的設(shè)計LLM可能無法準確表達對每個POI循環(huán)計算 } } ] } Executing Step 1: create_buffer with inputs {center_point: {lon: 116.4074, lat: 39.9042}, radius_km: 10} Result: {type: FeatureCollection, features: [{id: 0, type: Feature, properties: {}, geometry: {type: Polygon, coordinates:... Executing Step 2: fetch_pois_within_bounds with inputs {bounds: {west: 116.3174, south: 39.8142, east: 116.4974, north: 39.9942}, poi_type: hospital} Result: {type: FeatureCollection, features: [{type: Feature, geometry: {type: Point, coordinates: [116.367..., 39.934...]}, properti... Executing Step 3: calculate_distance with inputs {point1: {lon: 116.4074, lat: 39.9042}, point2: step_2} Error: Tool execution error: calculate_distance() argument after ** must be a mapping, not str問題暴露了我們的簡易代理在第三步失敗了因為LLM生成的計劃不夠精確它無法自動處理“對fetch_pois_within_bounds返回的每一個醫(yī)院點調(diào)用calculate_distance”這種循環(huán)邏輯。這引出了下一個關(guān)鍵章節(jié)常見問題與系統(tǒng)優(yōu)化。5. 常見問題、挑戰(zhàn)與優(yōu)化策略實錄在構(gòu)建和調(diào)試這樣一個地理空間代理框架時你會遇到一系列典型問題。以下是我在實踐中總結(jié)的“坑”和應(yīng)對策略。5.1 LLM規(guī)劃能力的局限性及應(yīng)對問題1LLM無法生成復(fù)雜的控制流如循環(huán)、條件判斷。如上例所示LLM擅長將任務(wù)分解為線性步驟但難以生成“對列表中的每個元素執(zhí)行某操作”這樣的代碼邏輯。它更傾向于輸出靜態(tài)的參數(shù)。解決方案A工具層面抽象創(chuàng)建更高級的復(fù)合工具。例如創(chuàng)建一個新工具calculate_distances_to_points(center_point, points_geojson)它內(nèi)部處理循環(huán)一次性返回所有距離。這樣LLM只需要調(diào)用這一個工具。解決方案B代理層面增強引入“子代理”或“遞歸規(guī)劃”。當主代理發(fā)現(xiàn)需要處理一個集合時它可以生成一個新的規(guī)劃子任務(wù)專門處理集合中的第一個元素并將模式應(yīng)用于其余元素。這需要更復(fù)雜的框架設(shè)計。解決方案C提示工程引導(dǎo)在工具描述中明確說明其處理集合的能力?；蛘咴谝?guī)劃提示詞中明確要求“如果某一步驟需要對一個列表中的每個項目進行操作請將該步驟命名為‘for_each_XXX’并說明循環(huán)的內(nèi)部操作。”問題2參數(shù)推斷不準確或模糊。用戶說“找附近的學(xué)?！盠LM可能無法推斷“附近”是多少米。解決方案在工具定義中提供合理的默認值和清晰的單位。同時框架應(yīng)支持多輪對話澄清。當代理無法確定關(guān)鍵參數(shù)時應(yīng)主動向用戶提問例如“您所說的‘附近’大概是指多少米范圍內(nèi)請?zhí)峁┮粋€具體數(shù)值例如500米、1公里或5公里。” 將用戶的回答補充到上下文中再繼續(xù)執(zhí)行。問題3工具選擇錯誤或順序混亂。LLM可能先調(diào)用需要A數(shù)據(jù)作為輸入的工具B但卻沒有先調(diào)用獲取A數(shù)據(jù)的工具A。解決方案在提供給LLM的工具描述中顯式聲明工具的前置條件和產(chǎn)出。例如fetch_pois_within_bounds的前置條件是“需要一個邊界框bounds”產(chǎn)出是“一個POI集合”。在規(guī)劃時可以要求LLM檢查每一步的輸入是否已被之前的步驟產(chǎn)出。更高級的框架會使用圖規(guī)劃算法來保證順序的正確性。5.2 地理空間數(shù)據(jù)處理的特殊挑戰(zhàn)問題4坐標系CRS混亂導(dǎo)致空間關(guān)系錯誤。這是地理分析中最常見、最隱蔽的錯誤。不同來源的數(shù)據(jù)WGS84經(jīng)緯度、Web墨卡托、各種地方坐標系混在一起計算結(jié)果毫無意義。解決方案建立嚴格的內(nèi)部CRS規(guī)范。我強烈建議在框架層面設(shè)定一個默認的、統(tǒng)一的工作坐標系例如EPSG:4326用于全球數(shù)據(jù)EPSG:3857用于Web地圖。所有工具在接收外部輸入時第一件事就是檢查并轉(zhuǎn)換到內(nèi)部CRS所有輸出時再根據(jù)需求轉(zhuǎn)換回目標CRS。在工具描述中必須注明“本工具要求輸入數(shù)據(jù)的CRS為EPSG:4326WGS84”。問題5大規(guī)模數(shù)據(jù)處理性能瓶頸。讓代理直接處理一個10GB的全國遙感影像是不現(xiàn)實的。解決方案采用懶加載和分塊處理策略。工具不直接處理原始大數(shù)據(jù)文件而是處理數(shù)據(jù)的引用如文件路徑、數(shù)據(jù)庫查詢、切片URL?？蚣軕?yīng)提供“數(shù)據(jù)加載”工具它可以根據(jù)分析范圍只讀取需要的那部分數(shù)據(jù)。對于必須全量處理的任務(wù)框架應(yīng)能生成可在高性能計算環(huán)境如Spark集群上運行的腳本而不是在交互式代理中直接運行。問題6地理操作的多樣性和復(fù)雜性。一個“疊加分析”可能意味著相交intersection、聯(lián)合union、擦除difference等不同操作。解決方案工具設(shè)計要粒度適中功能單一。不要設(shè)計一個萬能的spatial_analysis工具而是設(shè)計intersect_features,union_features,erase_features等具體工具。這樣LLM更容易理解和匹配。同時提供詳盡的工具描述和示例。5.3 系統(tǒng)魯棒性與用戶體驗問題7工具執(zhí)行失敗后的處理。網(wǎng)絡(luò)超時、數(shù)據(jù)不存在、參數(shù)越界等都可能導(dǎo)致單個工具調(diào)用失敗。解決方案實現(xiàn)重試機制和備選方案。例如獲取POI的API失敗后可以嘗試換一個備用數(shù)據(jù)源。更重要的是框架需要將詳細的錯誤信息包括堆棧跟蹤進行摘要反饋給LLM讓它有機會重新規(guī)劃Replan。例如錯誤是“坐標超出數(shù)據(jù)范圍”LLM可能會推斷出需要先獲取一個更大范圍的基礎(chǔ)數(shù)據(jù)。問題8如何呈現(xiàn)復(fù)雜的地理結(jié)果最終輸出可能是一個GeoJSON文件、一張統(tǒng)計圖表、或一段文字報告。解決方案提供多樣化的結(jié)果渲染工具。除了返回原始數(shù)據(jù)框架應(yīng)集成如plot_map用Folium/Matplotlib生成交互地圖、generate_summary_statistics生成文本摘要、export_to_geojson_file導(dǎo)出文件等工具。讓代理在最后一步根據(jù)用戶指令“把結(jié)果在地圖上標出來”選擇合適的渲染方式。問題9代理的“幻覺”問題。LLM可能會編造一個不存在的工具或參數(shù)。解決方案實施嚴格的工具驗證。在代理執(zhí)行計劃前對每一步的tool_name進行校驗確保其在注冊的工具列表中。對輸入?yún)?shù)的類型和范圍進行基礎(chǔ)校驗。這屬于“護欄”設(shè)計是生產(chǎn)級系統(tǒng)必不可少的。構(gòu)建OpenEarthAgent這樣的框架是一個在AI的靈活性與地理計算的嚴謹性之間尋找平衡的藝術(shù)。它不是一個能解決所有問題的“銀彈”而是一個強大的“力量倍增器”能將地理空間專家的知識沉淀為可復(fù)用的工具并通過自然語言界面釋放給更廣泛的用戶。從我們上面的簡易實現(xiàn)可以看出核心難點不在于單個工具的實現(xiàn)而在于如何讓AI可靠地、準確地組合它們。這需要精心的工具設(shè)計、巧妙的提示工程以及健壯的框架邏輯。隨著多模態(tài)大模型和代碼生成能力的進步我相信這類框架會越來越成熟最終讓每個人都能像專家一樣進行空間思考與分析。