如何使用pythong根据站点的经纬度值画广东省的站点图呢？

神奇小宝包

我想请问一下，如何使用pythong根据站点的经纬度值画广东省的站点图呢？

神奇小宝包

看了相关的例子还是画不出来

南瓜是我

画的40N以上的站点图，文件为csv格式。仅供参考

1. import pandas as pd
   
2. import matplotlib.pyplot as plt
   
3. import cartopy.crs as ccrs
   
4. import cartopy.feature as cfeature
   
5. import cartopy.mpl.ticker as cticker
   
6. import cartopy.io.shapereader as shpreader
   
7. import numpy as np
   
8. 
   
9. data = pd.read_csv("G:/GHCND/ghcnd_all/ghcnd-stations-60E-170W_40-90.csv",sep=',', \
   
10.                         header=None, names=['station','lat','lon','altitude','country','number'])  
    
11. 
    
12. #建立画布
    
13. fig2 = plt.figure(figsize=(10,8))
    
14. proj = ccrs.PlateCarree(central_longitude=110)
    
15. leftlon, rightlon, lowerlat, upperlat = (60,190,0,90)
    
16. #绘制地图
    
17. f2_ax1 = fig2.add_subplot(1,1,1,projection = proj)
    
18. #在画布的绝对坐标建立子图
    
19. f2_ax1.set_extent([leftlon, rightlon, lowerlat, upperlat], crs=ccrs.PlateCarree())
    
20. #海岸线，50m精度
    
21. f2_ax1.add_feature(cfeature.COASTLINE.with_scale('50m'))
    
22. 
    
23. #以下6条语句是定义地理坐标标签格式
    
24. f2_ax1.set_xticks(np.arange(leftlon,rightlon,20), crs=ccrs.PlateCarree())
    
25. f2_ax1.set_yticks(np.arange(lowerlat,upperlat,20), crs=ccrs.PlateCarree())
    
26. lon_formatter = cticker.LongitudeFormatter()
    
27. lat_formatter = cticker.LatitudeFormatter()
    
28. 
    
29. f2_ax1.xaxis.set_major_formatter(lon_formatter)
    
30. f2_ax1.yaxis.set_major_formatter(lat_formatter)
    
31. 
    
32. f2_ax1.set_title('Station',loc='center',fontsize =15)
    
33. 
    
34. s = f2_ax1.scatter(data.lon,data.lat,s = 5,c='k',transform=ccrs.PlateCarree())
    
35. 
    
36. #fig2.savefig('./pic/station40.jpg')
    
37. plt.show()
    

复制代码

神奇小宝包

谢谢大家！我已经会啦

yuyuci

楼主会了吗，最近要画华南地区的站点图，不知道怎么下手，求解惑

神奇小宝包

yuyuci 发表于 2021-4-6 23:13
楼主会了吗，最近要画华南地区的站点图，不知道怎么下手，求解惑

先抠出华南地区的shp图，然后使用scatter画就好啦！