Hello, I am new to yt and I am interesting in plotting a volume rendering with output file from enzo cosmology simulation. For the time being, I get to plot a box ; here's the result : http://imgur.com/dPwSQEh I would like now to plot borders of this box too like on this image : http://imgur.com/bJ25I0l Here's the yt script I use : ------------------------------------------------------------------------------------------------------------ from yt.mods import * from yt.utilities.amr_kdtree.api import AMRKDTree from time import time import matplotlib.colorbar as cb # Load up your dataset num='70' pf = load('RD00'+num+'/RedshiftOutput00'+num) c = [0.5]*3 # Center L = [1.0,1.0,1.0] # Viewpoint W = na.sqrt(3) # Width N = 768 # Pixels (512^2) up = [0.,0.,1.] # Get density min, max # These might take a long time, so I'd suggest hand-coding mi, ma so # that you don't have to find the maxima and minima each time you # test. Note that mi, ma, should be a log value since we are # rendering the log of density. mi, ma = pf.h.all_data().quantities['Extrema']('Density')[0] mi, ma = na.log10(mi), na.log10(ma) print mi, ma #mi, ma = -31.0, -26.0 # Construct transfer function, pad the TF space by a bit so that # gaussians sampling the data range don't hit the edge. tf = ColorTransferFunction((mi-10, ma+10), nbins=1024) # Sample transfer function with 20 gaussians. Use col_bounds keyword # to restrict color mapping to true data range. tf.add_layers(20, w=0.001, col_bounds=(mi,ma),colormap='jet') # Create the camera object cam = pf.h.camera(c, L, W, (N,N), transfer_function=tf, north_vector=up, no_ghost=True) cam.snapshot("volume_rendered%s.png" % num, clip_ratio=8.0) ------------------------------------------------------------------------------------------------------------ Is there an option for 'camera" or "snapshot" which allows to add the borders for this box ? Any help would be really appreciated.
participants (2)
-
henry wasker
-
Sam Skillman