pyshp读取矢量图层，用pyproj获取图层原始坐标系，并投影

要使用pyshp读取矢量图层，并使用pyproj获取图层的原始坐标系并进行投影，你可以按照以下步骤操作：

1. 安装所需的库：
- `pyshp`：用于读取矢量数据的库
- `pyproj`：用于进行坐标系投影的库
- `shapely`：用于处理几何对象的库

你可以使用以下命令安装这些库：

   pip install pyshp pyproj shapely

2. 导入所需的库和函数：

   import shapefile
   from pyproj import Proj, Transformer
   from shapely.geometry import shape

3. 读取矢量图层数据：

   # 打开矢量图层文件
   sf = shapefile.Reader('your_shapefile.shp')

   # 获取矢量图层的几何对象和属性数据
   shapes = sf.shapes()
   records = sf.records()

4. 获取图层的原始坐标系：

   # 获取矢量图层的原始坐标系信息
   original_crs = sf.crs

   # 打印原始坐标系信息
   print(original_crs)

5. 定义目标坐标系并创建投影转换对象：

   # 定义目标坐标系（例如 EPSG:3857）
   target_crs = 'EPSG:3857'

   # 创建投影转换对象
   transformer = Transformer.from_crs(original_crs, target_crs, always_xy=True)

6. 遍历矢量图层的几何对象，并进行投影转换：

   for shape_obj, record in zip(shapes, records):
       # 将几何对象转换为Shapely对象
       shapely_obj = shape(shape_obj)

       # 对几何对象进行投影转换
       projected_obj = shapely_obj.transform(transformer.transform)

       # 打印投影后的几何对象坐标
       print(projected_obj)

在这个示例中，我们首先使用`shapefile`模块的`Reader`函数打开矢量图层文件，并使用`shapes`和`records`属性获取几何对象和属性数据。然后，我们使用`crs`属性获取图层的原始坐标系信息，并打印输出。

接下来，我们定义了目标坐标系（例如EPSG:3857），并使用`pyproj`库的`Transformer`类创建了一个投影转换对象。

最后，我们遍历了矢量图层的几何对象和对应的属性数据，将每个几何对象转换为Shapely对象，并使用投影转换对象对其进行投影转换。你可以根据需要修改代码中的文件路径和目标坐标系。

请确保你已经安装了所需的库，并替换代码中的`your_shapefile.shp`为你自己的矢量图层文件路径。