Geopandas：带有不同ESPG代码的GeodataFrame的数据框

Question

我有一个dataframe（DF2）：其中x，y在RD新epsg中指定：28992坐标。

            x           y      z     batch_nr  batch_description
0  117298.377  560406.392  0.612         5800   PRF Grasland (l)
1  117297.803  560411.756  1.015                                
2  117296.327  560419.840  1.580                                
3  117295.470  560425.716  2.490                                
4  117296.875  560429.976  4.529

更多CRS信息：

# def CRS, used in geopandas
from pyproj import CRS
crs_rd = CRS.from_user_input(28992)
crs_rd

<Derived Projected CRS: EPSG:28992>
Name: Amersfoort / RD New
Axis Info [cartesian]:
- X[east]: Easting (metre)
- Y[north]: Northing (metre)
Area of Use:
- name: Netherlands - onshore, including Waddenzee, Dutch Wadden Islands and 12-mile offshore coastal zone.
- bounds: (3.2, 50.75, 7.22, 53.7)
Coordinate Operation:
- name: RD New
- method: Oblique Stereographic
Datum: Amersfoort
- Ellipsoid: Bessel 1841
- Prime Meridian: Greenwich

如何将DF2转换为GeodataFame，其中将几何图形设置为CRS：EPSG 28992？

Answer 1

这是使用 geopandas crs参数和points_from_xy（）

import geopandas as gpd
import pandas as pd
import io

df2 = pd.read_csv(io.StringIO("""            x           y      z     batch_nr  batch_description
0  117298.377  560406.392  0.612         5800   PRF Grasland (l)
1  117297.803  560411.756  1.015                                
2  117296.327  560419.840  1.580                                
3  117295.470  560425.716  2.490                                
4  117296.875  560429.976  4.529"""), sep="\s\s+", engine="python")

gdf = gpd.GeoDataFrame(df2, geometry=gpd.points_from_xy(df2["x"], df2["y"], df2["z"]), crs="epsg:28992")

gdf

输出

	x	y	z	batch_nr	batch_description	几何
0	117298	560406	0.612	5800	PRF GRASLAND GRASLAND 的简单情况。 （L）	Z点（117298.377 560406.392 0.612）
560412	1.015	NAN	点	Z		117297.803 560411.756 1.015） 26
Point	Z	（	2.49	（		117295.47
117298	1	2.49	560425.716	4.529		）
4	117297	560430	Nan	Nan		Point Z（117296.875 560429.976 4.529）