经纬度和不同高度消光系数数据绘制三维路径图请教

Linda.W

王老师，请问有经纬度数据和对应的不同高度消光系数数据怎么用MI绘制三维路径图，试着写了程序但是没有出图也没有报错，
附上数据。

MeteoInfo

你说的三位路径图具体指什么样的图？

Linda.W

王老师，我想绘制成示例图片这种形式，这是个示意图，想叠加进三维坐标系中，目前出的图可能数据读取的问题只能构建出坐标系。

MeteoInfo

Linda.W 发表于 2023-4-4 13:35
王老师，我想绘制成示例图片这种形式，这是个示意图，想叠加进三维坐标系中，目前出的图可能数据读取的问题 ...

参考这个脚本，把excel文件另存为csv文件，注意先要把excel文件中右边多余的列删除。

1. fn = 'D:/Temp/ascii/ws-aod.csv'
   
2. df = DataFrame.read_table(fn, delimiter=',', format='%422f')
   
3. lon = df['Longitude'].values
   
4. lat = df['Latitude'].values
   
5. alt = array(df.columns.names[2:], 'float')
   
6. n_alt = len(alt)
   
7. alt, lon = meshgrid(alt, lon)
   
8. lat = np.expand_dims(lat, axis=1)
   
9. lat = np.repeat(lat, n_alt, axis=1)
   
10. ec = df[:,2:].values
    
11. 
    
12. axes3d()
    
13. surf(lon, lat, alt, ec, edgecolor=None, facecolor='interp')
    
14. colorbar

fn = 'D:/Temp/ascii/ws-aod.csv'<br /> df = DataFrame.read_table(fn, delimiter=',', format='%422f')<br /> lon = df['Longitude'].values<br /> lat = df['Latitude'].values<br /> alt = array(df.columns.names[2:], 'float')<br /> n_alt = len(alt)<br /> alt, lon = meshgrid(alt, lon)<br /> lat = np.expand_dims(lat, axis=1)<br /> lat = np.repeat(lat, n_alt, axis=1)<br /> ec = df[:,2:].values<br /> <br /> axes3d()<br /> surf(lon, lat, alt, ec, edgecolor=None, facecolor='interp')<br /> colorbar

Linda.W

王老师，试了多个版本，运行到surf函数的时候总是报错，希望您能再指导指导
fn =  r'D:\test\data-3D\ws6.2-1aod.csv'
df = DataFrame.read_table(fn, delimiter=',', format='%422f')
lon = df['Longtitude'].values
lat = df['Latitude'].values
alt = array(df.columns.names[2:533], 'float')
n_alt = len(alt)
print n_alt
alt, lon = meshgrid(alt, lon)
lat = np.expand_dims(lat, axis=1)
lat = np.repeat(lat, n_alt, axis=1)
ec = df[:,2:533].values

ax=axes3d()
surf(lon, lat, alt, ec, edgecolor=None, facecolor='interp')

xlim(114.2,114.4)
ylim(34.2,34.6)
zlim(0,3000)
colorbar()
savefig(r'D:\test\ws6.2aod-3D.jpg',dpi = 300,width=800, height=600)

Linda.W

王老师，我这个运行还是会报错，可以帮我看看是什么原因吗？
fn =  r'D:\test\data-3D\ws6.2-1aod.csv'
df = DataFrame.read_table(fn, delimiter=',', format='%422f')
lon = df['Longtitude'].values
lat = df['Latitude'].values
alt = array(df.columns.names[2:533], 'float')
n_alt = len(alt)
print n_alt
alt, lon = meshgrid(alt, lon)
lat = np.expand_dims(lat, axis=1)
lat = np.repeat(lat, n_alt, axis=1)
ec = df[:,2:533].values

ax=axes3d()
surf(lon, lat, alt, ec, edgecolor=None, facecolor='interp')

xlim(114.2,114.4)
ylim(34.2,34.6)
zlim(0,3000)
colorbar()
savefig(r'D:\test\ws6.2aod-3D.jpg',dpi = 300,width=800, height=600)

MeteoInfo

Linda.W 发表于 2023-4-11 08:55
王老师，我这个运行还是会报错，可以帮我看看是什么原因吗？
fn =  r'D:\test\data-3D\ws6.2-1aod.csv'
d ...

你的数据有多少列？422还是533？

Linda.W

源文件数据放了600列，读取的时候只用了前533列数据

Linda.W

王老师，是不是format='%422f'这是限定了长度？我之前没懂这个代码的意思，我试试改一下

Linda.W

好像不是列数的问题，因为我改成surf(lon, lat, alt, c=ec, edgecolor=None, facecolor='interp')的时候，能出高度图