function box_image, PROJ=proj
    alpha = proj.slice.alpha
    beta = proj.slice.beta
    w = proj.slice.image_size
    axes = proj.slice.axes
    ; draw box lines
    set_plot, 'z'
    DEVICE, SET_RESOLUTION=[w,w]
    erase, 0
    save_x = !X & save_y = !Y & save_p = !P
    !X.margin = 0. & !Y.margin = 0. & !X.thick = 0.  & !Y.thick = 0.
    !X.ticks = 1   & !Y.ticks = 1 & !X.range = [0.,1.] & !Y.range = [0.,1.]
    !X.style=4 & !Y.style=4 & !P.ticklen = 0 & !P.T3D = 0
    ca = cos(alpha) & sa = sin(alpha) & cb = cos(beta) & sb = sin(beta)
    T = [[1,0,0,.5], [0,1,0,.5], [0,0,1,0], [0,0,0,1]] ## $
        [[1,0,0,0], [0,-cb,sb,0], [0,-sb,-cb,0], [0,0,0,1]] ## $
        [[sa,-ca,0,0], [ca,sa,0,0], [0,0,1,0], [0,0,0,1]]
    X = .5*[0,0,0,0,1,1,1,1] - .25
    Y = .5*[0,0,1,1,0,0,1,1] - .25
    Z = .5*[0,1,0,1,0,1,0,1] - .25
    O = [1,1,1,1,1,1,1,1]
    XYZ = T ## [[X],[Y],[Z],[O]]
    order = [0,1,5,4,0,2,6,4,5,7,3,1,3,2,6,7]
    PLOT, XYZ[order,0], XYZ[order,1]
    ; draw info text
    xyouts, 0, 0, STRING(FORMAT='(I0,I0,I0,F5.1,F5.1)', axes[0], axes[1], axes[2], alpha, beta)
    image = TVRD(TRUE=0)
    !X = save_x & !Y = save_y & !P = save_p
    set_plot, 'x'
    return, image / max(image)
end