0

NetCDF ファイルからの降水データに基づいて、地域の降水レベルのマップを作成しました。降水量が 800mm 未満の場合は 1 つの色、800 ～ 1000mm の場合は別の色になるように、カスタム スケールを追加したいと思います。http: //www.metmalawi.com/climate/climate. php

現在、グラデーション スケールを使用していますが、必要な詳細が表示されません。これは、現時点でのプロットのコードです (「平均」は、既にフォーマットしたデータです)。

    #load color palette
    colourA = mpl_cm.get_cmap('BuPu')
    
    #plot map with physical features 
    ax = plt.axes(projection=cartopy.crs.PlateCarree())
    ax.add_feature(cartopy.feature.COASTLINE)   
    ax.add_feature(cartopy.feature.BORDERS)
    ax.add_feature(cartopy.feature.LAKES, alpha=0.5)
    ax.add_feature(cartopy.feature.RIVERS)

    #set map boundary
    ax.set_extent([32.5, 36., -9, -17]) 
    
    #set axis tick marks
    ax.set_xticks([33, 34, 35]) 
    ax.set_yticks([-10, -12, -14, -16]) 
    lon_formatter = LongitudeFormatter(zero_direction_label=True)
    lat_formatter = LatitudeFormatter()
    ax.xaxis.set_major_formatter(lon_formatter)
    ax.yaxis.set_major_formatter(lat_formatter)
    
    #plot data and set colour range
    plot = iplt.contourf(Average, cmap=colourA, levels=np.arange(0,15500,500), extend='both')

    #add colour bar index and a label
    plt.colorbar(plot, label='mm per year')

    #give map a title
    plt.title('Pr 1990-2008 - Average_ERAINT ', fontsize=10)

    #save the image of the graph and include full legend
    plt.savefig('ERAINT_Average_Pr_MAP_Annual', bbox_inches='tight')

    plt.show()

どうすればこれができるか知っている人はいますか？

ありがとうございました！

4

1 に答える 1