GIS小技术分享：Python中JSON数据转GeoJSON或SHP格式

在地理信息系统（GIS）中，数据通常以不同的格式存储和交换，如Shapefile（SHP）、GeoJSON等。有时，我们可能会遇到JSON格式的数据，并需要将其转换为GIS可用的格式。在Python中，我们可以使用一些库来轻松完成这个任务。
首先，确保你已经安装了必要的Python库。你可以使用pip来安装Fiona和GeoPandas：

pip install fiona
pip install geopandas

下面是一个简单的示例，演示如何将JSON数据转换为GeoJSON格式：

import json
import fiona
from shapely.geometry import shape
# 加载JSON数据
with open('data.json', 'r') as f:
data = json.load(f)
# 创建空的GeoJSON对象
geojson = {'type': 'FeatureCollection', 'features': []}
# 遍历JSON数据中的每个要素
for feature in data['features']:
# 将每个要素的geometry字段转换为Shapely几何对象
geometry = shape(feature['geometry'])
# 创建一个新的GeoJSON要素，并将几何对象添加到其中
geojson['features'].append({
'type': 'Feature',
'geometry': geometry.__geo_interface__,
'properties': feature['properties']
})
# 将GeoJSON对象写入文件
with open('output.geojson', 'w') as f:
json.dump(geojson, f)

在这个示例中，我们首先打开并读取JSON文件中的数据。然后，我们创建一个空的GeoJSON对象，遍历JSON数据中的每个要素，并将其几何对象转换为Shapely几何对象。最后，我们将每个要素添加到GeoJSON对象中，并将其写入文件。
要将GeoJSON转换为SHP格式，可以使用Fiona库。首先，确保你已经安装了Fiona：

pip install fiona

然后，你可以使用以下代码将GeoJSON转换为SHP：

import fiona
from fiona.crs import from_epsg
from shapely.geometry import shape, mapping, asShape, Polygon, mapping, Point, LineString, MultiLineString, MultiPoint, MultiPolygon, GeometryCollection, mapping, box2d, box3d, from_shape, from_shapefile, shapefile_box2d, shapefile_box3d, shapefile_points_box2d, shapefile_points_box3d, shapefile_multipoints_box2d, shapefile_multipoints_box3d, shapefile_multilines_box2d, shapefile_multilines_box3d, shapefile_multipolygons_box2d, shapefile_multipolygons_box3d, shapefile_geometrycollection_box2d, shapefile_geometrycollection_box3d, shapefile_multipolygons_box2d, shapefile_multipolygons_box3d, shapefile_multilines_box2d, shapefile_multilines_box3d, shapefile_geometrycollection_box2d, shapefile_geometrycollection_box3d