gplot1D.py (2) ヒストグラム (ステップ)

• ヒストグラムを書くが、バー（棒グラフ）でなく、ステップでプロットするためのデータ変換

• 以下の手順で作成可能．

  • データレンジ（xMin, xMax）を連結（２次元に連結⇒１次元になおす）

  • y値のデータを２回繰り返す．

  • 端処理： y値=0 、x値=両端の値をコピー

グラフ

コード

import os, sys
import numpy as np
import nkUtilities.load__config   as lcf
import nkUtilities.gplot1D        as gp1

# ========================================================= #
# ===  draw__step                                       === #
# ========================================================= #


def draw__step():

    # ------------------------------------------------- #
    # --- [1] generate sample data                  --- #
    # ------------------------------------------------- #
    xAxis1 = np.linspace( 0,  1,  6 )
    xAxis2 = np.linspace( 2, 40, 20 )
    xMins  = np.concatenate( [ xAxis1[:-1] , xAxis2[:-1] ] )
    xMaxs  = np.concatenate( [ xAxis1[ 1:] , xAxis2[ 1:] ] )
    xAxis  = 0.5 * ( xMins + xMaxs )
    yAxis  = (xAxis)/np.sqrt( 2.0*np.pi*20.0 ) * np.exp( -0.5*( ( (xAxis-10.0)/20.0 )**2 ) )

    # ------------------------------------------------- #
    # --- [2] rearange for step plot                --- #
    # ------------------------------------------------- #
    xAxis_ = np.reshape( np.concatenate( [ xMins[:,np.newaxis], \
                                           xMaxs[:,np.newaxis] ], axis=1 ), (-1,) )
    yAxis_ = np.repeat( yAxis, 2 )
    xAxis_ = np.insert( xAxis_, [0,xAxis_.size], [ xAxis[0], xAxis_[-1] ] )
    yAxis_ = np.insert( yAxis_, [0,yAxis_.size], [ 0.0, 0.0        ] )
    
    # ------------------------------------------------- #
    # --- [2] plot                                  --- #
    # ------------------------------------------------- #
    config   = lcf.load__config()
    config_  = {
        "figure.size"        : [4.5,4.5],
        "figure.pngFile"     : "png/gplot1d__p1_02.png", 
        "figure.position"    : [ 0.16, 0.16, 0.94, 0.94 ],
        "ax1.y.normalize"    : 1.0e0, 
        "ax1.x.range"        : { "auto":False, "min": 0.0, "max":50.0, "num": 6 },
        "ax1.y.range"        : { "auto":False, "min": 0.0, "max": 2.0, "num": 5 },
        "ax1.x.label"        : "x",
        "ax1.y.label"        : "y",
        "ax1.x.minor.nticks" : 2, 
        "ax1.y.minor.nticks" : 5, 
        "plot.markersize"    : 3.0,
        "legend.fontsize"    : 9.0, 
    }
    config = { **config, **config_ }
    fig    = gp1.gplot1D( config=config )
    fig.add__plot( xAxis=xAxis_, yAxis=yAxis_ )
    fig.set__axis()
    fig.save__figure()

    
# ========================================================= #
# ===   Execution of Pragram                            === #
# ========================================================= #

if ( __name__=="__main__" ):
    draw__step()