python基于站点经纬度绘制降水空间散点分布图
导入库并读取数据import numpy as npimport pandas as pdimport cartopy.crs as ccrsimport cartopy.feature as cfeatureimport matplotlib.pyplot as pltfrom matplotlib import colors, cmimport cmapsimport geocat.viz as gvfrom cartopy.io.shapereader import Readerdata=pd.read_csv('20190722.csv',dtype=np.float64,header=None,delimiter=',',encoding='gbk')lat=np.array(data)lon=np.array(data)rain=np.array(data)设置colorbar刻度及区间色调scales = cmap = cmaps.rainbowboundaries = norm = colors.BoundaryNorm(boundaries, cmap.N)mappable = cm.ScalarMappable(norm=norm, cmap=cmap)
设置散点标记的颜色区间marker_colors = mappable.to_rgba(boundaries)sizes = np.geomspace(10, 250, len(boundaries))plt.figure(figsize=(9, 6))projection = ccrs.PlateCarree()ax = plt.axes(projection=projection)ax.set_extent(, crs=projection)
添加四川地图shap=Reader('SCmap.shp').geometries()sichuan = cfeature.ShapelyFeature(shap,crs=ccrs.PlateCarree(),edgecolor='k', facecolor='none')ax.add_feature(sichuan)
设置x、y轴经纬度刻度gv.set_axes_limits_and_ticks(ax,xticks=np.linspace(97, 109, 5),yticks=np.linspace(26, 34, 5))gv.add_lat_lon_ticklabels(ax)gv.add_major_minor_ticks(ax,x_minor_per_major=1,y_minor_per_major=1,labelsize=12)# Remove ticks on the top and right sides of the plotax.tick_params(axis='both', which='both', top=False, right=False)
绘制不同降水区间散点图masked_lon = np.where(rain < scales, lon, np.nan)masked_lat = np.where(rain < scales, lat, np.nan)plt.scatter(masked_lon,masked_lat,s=sizes,color=marker_colors,zorder=1)for x in range(1, len(scales)): masked_lon = np.where(rain >= scales, lon, np.nan) masked_lon = np.where(rain < scales, masked_lon, np.nan) masked_lat = np.where(rain >= scales, lat, np.nan) masked_lat = np.where(rain < scales, masked_lat, np.nan) plt.scatter(masked_lon,masked_lat,s=sizes,color=marker_colors,zorder=1)masked_lon = np.where(rain >= scales[-1], lon, np.nan)masked_lat = np.where(rain >= scales[-1], lat, np.nan)plt.scatter(masked_lon,masked_lat,s=sizes[-1],color=marker_colors[-1],zorder=1)
考标记出某一站点plt.colorbar(mappable=mappable,ax=ax,orientation='horizontal',label='Rainfall Amount(mm)', drawedges=True,format='%.2f',ticks=scales)plt.scatter(103.12,30.08,s=20)plt.annotate(r'$mingshan$', xy=(103.12,30.08),xytext=(4,-100),xycoords='data',textcoords='offset points', fontsize=16,arrowprops=dict(arrowstyle='->',connectionstyle='arc3'))plt.savefig('test.png')
页:
[1]