获取在线地图上行政边界

要批量使用Python在百度地图上获取矢量行政边界并将坐标系转换为WGS84，可以使用百度地图开放平台提供的API和GDAL库。以下是一些基本步骤：

1. 在百度地图开放平台上注册并获取开发者密钥。
2. 使用Python的requests库向百度地图API发送请求以获取矢量行政边界的数据。以下是一个代码片段：

   import requests
   
   def get_admin_boundary(city):
       url = f"http://api.map.baidu.com/boundary/v1/boundary?city={city}&key=your_key"
       response = requests.get(url)
       if response.status_code == 200:
           data = response.json()
           return data
       else:
           return None
 

   请将 `your_key` 替换为您在百度地图开放平台上获取的开发者密钥。

3. 使用GDAL库将矢量数据转换为WGS84坐标系。以下是一个代码片段：

   
   from osgeo import ogr, osr
   
   def transform_to_wgs84(input_file, output_file):
       # 打开输入文件
       input_ds = ogr.Open(input_file)
       input_layer = input_ds.GetLayer()
   
       # 定义输出坐标系
       output_srs = osr.SpatialReference()
       output_srs.ImportFromEPSG(4326)  # WGS84坐标系的EPSG代码是4326
   
       # 创建输出文件
       output_driver = ogr.GetDriverByName("GeoJSON")
       output_ds = output_driver.CreateDataSource(output_file)
       output_layer = output_ds.CreateLayer("admin_boundary", srs=output_srs, geom_type=ogr.wkbMultiPolygon)
   
       # 添加字段
       field_defn = ogr.FieldDefn("name", ogr.OFTString)
       field_defn.SetWidth(50)
       output_layer.CreateField(field_defn)
   
       # 转换坐标系并写入输出文件
       for feature in input_layer:
           geometry = feature.GetGeometryRef()
           geometry.TransformTo(output_srs)
   
           output_feature = ogr.Feature(output_layer.GetLayerDefn())
           output_feature.SetGeometry(geometry)
           output_feature.SetField("name", feature.GetField("name"))
           output_layer.CreateFeature(output_feature)
   
       # 关闭文件
       input_ds = None
       output_ds = None
 

4. 最后，可以编写一个简单的脚本来循环处理每个城市的矢量行政边界。以下是一个Python脚本：

   import os
   
   input_folder = "path/to/input/folder"
   output_folder = "path/to/output/folder"
   
   for city in ["北京市", "上海市", "广州市"]:
       input_file = os.path.join(input_folder, f"{city}.json")
       output_file = os.path.join(output_folder, f"{city}_wgs84.geojson")
   
       data = get_admin_boundary(city)
       with open(input_file, "w") as f:
           f.write(data["geo"])
   
       transform_to_wgs84(input_file, output_file)
 

   请将 `input_folder` 和 `output_folder` 替换为输入文件夹和输出文件夹的路径。在此示例中，假设您已经在输入文件夹中下载了包含每个城市名称的JSON文件，并将其命名为城市名称。例如，北京市的JSON文件名为 `北京市.json`。，仅处理了三个城市，可以根据需要添加或删除城市名称。