用文件中的数据（纬度，经度，海拔）画ECEF散点图

要在ECEF（地球中心惯性参考系）坐标系下使用文件中的数据（纬度、经度、海拔）创建散点图，你需要先将这些地理坐标转换为ECEF坐标，然后才能在三维空间中绘制。以下是一般的步骤：

1. **导入必要的库**：

   import numpy as np
   from pyproj import Transformer
   import matplotlib.pyplot as plt

2. **设置坐标转换**：
使用`pyproj.Transformer`，创建一个从地理坐标到ECEF坐标的转换器：

   transformer = Transformer.from_crs("epsg:4326", "crs:ECEF", always_xy=True)

其中，"epsg:4326"是WGS-84坐标系统，"crs:ECEF"是ECEF坐标系统。

3. **加载并转换数据**：
读取文件数据，然后逐行转换为ECEF坐标：

   df = pd.read_csv('your_data_file.csv')  # 假设数据在CSV文件中
   latitudes, longitudes, altitudes = df['latitude'].values, df['longitude'].values, df['altitude'].values
   ecef_x, ecef_y, ecef_z = [], [], []

   for lat, lon, alt in zip(latitudes, longitudes, altitudes):
       x, y, z = transformer.transform(lon, lat, alt)
       ecef_x.append(x)
       ecef_y.append(y)
       ecef_z.append(z)

4. **创建并显示散点图**：

   fig = plt.figure()
   ax = fig.add_subplot(projection='3d')
   ax.scatter(ecef_x, ecef_y, ecef_z, s=5)  # s控制点的大小
   ax.set_xlabel('X (m)')
   ax.set_ylabel('Y (m)')
   ax.set_zlabel('Z (m)')
   plt.title('ECEF Coordinate System Scatter Plot')
   plt.show()

记得替换`your_data_file.csv`为你实际的数据文件名。