Python中shapefile转换geojson的示例

下面为你详细讲解“Python中shapefile转换geojson的示例”的完整攻略：

1. 安装依赖

首先，你需要安装以下两个Python库以进行shapefile和geojson的转换操作：

• pyshp： 用于读取和写入shapefile文件
• geojson： 用于读取和写入geojson文件

你可以通过以下命令来安装这两个库：

pip install pyshp
pip install geojson

2. 加载shapefile文件

接下来，你需要使用pyshp库中的Reader对象来加载shapefile文件。以下代码展示了如何加载一个名为myshapefile的shapefile文件：

import shapefile

sf = shapefile.Reader("myshapefile.shp")

3. 将shapefile转换为geojson

使用pyshp库读取shapefile文件后，我们可以使用以下代码将其转换为geojson格式：

import shapefile
import geojson

sf = shapefile.Reader("myshapefile.shp")

# 创建一个空的特征集合
feature_collection = geojson.FeatureCollection([])

# 遍历shapefile中的每个记录
for shape in sf.iterShapeRecords():

    # 将记录中的几何对象转换为geojson格式
    geometry = geojson.loads(shape.__geo_interface__)

    # 创建一个包含geojson几何对象的特征
    feature = geojson.Feature(geometry=geometry["geometry"], properties={"id": shape.record[0], "name": shape.record[1]})

    # 将特征添加到特征集合中
    feature_collection["features"].append(feature)

# 将特征集合输出为geojson格式文件
with open("mygeojsonfile.geojson", "w") as f:
    f.write(geojson.dumps(feature_collection))

以上代码将加载的shapefile文件转换为一个包含多个geojson特征的特征集合，其中每个记录都将转换为一个包含几何对象和属性信息的geojson特征。

在本代码示例中，每个记录的几何对象是shape.__geo_interface__属性的值，包含geometry和type属性，而属性信息是通过shape.record属性获取的。

示例1: 将一个多边形的shapefile转换为一个geojson文件

假设在当前工作目录下有一个名为mypolygon.shp的多边形shapefile，我们希望将其转换为一个geojson文件。下面是一个示例代码：

import shapefile
import geojson

#加载多边形shapefile
sf = shapefile.Reader("mypolygon.shp")

#创建空特征集
feature_collection = geojson.FeatureCollection([])

#循环遍历每个记录
for shape in sf.iterShapeRecords():
    #将几何对象转换为geojson格式
    geometry = geojson.loads(shape.__geo_interface__)

    #创建geojson特征
    feature = geojson.Feature(geometry=geometry["geometry"],properties={"id": shape.record[0], "name": shape.record[1]})

    #将特征添加到集合
    feature_collection["features"].append(feature)

#将特征集合输出到geojson文件
with open("mypolygon.geojson", "w") as f:
    f.write(geojson.dumps(feature_collection))

上述代码将转换一个多边形shapefile为一个GeoJSON文件，并将其保存在mypolygon.geojson文件中。

示例2: 将一个点的shapefile转换为一个geojson文件

假设在当前工作目录下有一个名为mypoint.shp的点的shapefile，我们希望将其转换为一个geojson文件。下面是一个示例代码：

import shapefile
import geojson

#加载点shapefile
sf = shapefile.Reader("mypoint.shp")

#创建空特征集
feature_collection = geojson.FeatureCollection([])

#循环遍历每个记录
for shape in sf.iterShapeRecords():
    #将几何对象转换为geojson格式
    geometry = geojson.loads(shape.__geo_interface__)

    #创建geojson特征
    feature = geojson.Feature(geometry=geometry["geometry"],properties={"id": shape.record[0], "name": shape.record[1]})

    #将特征添加到集合
    feature_collection["features"].append(feature)

#将特征集合输出到geojson文件
with open("mypoint.geojson", "w") as f:
    f.write(geojson.dumps(feature_collection))

上述代码将转换一个点的shapefile为一个GeoJSON文件，并将其保存在mypoint.geojson文件中。

总之，以上就是Python中shapefile转换geojson的完整攻略，包含安装依赖库、加载shapefile文件和将shapefile转换为GeoJSON文件的过程，以及两个实际操作示例的详细步骤说明。