Version 38 (modified by 9 years ago) ( diff ) | ,
---|
Under construction
Visualizing External Data in Shape
As we know very well, AstroBEAR generates simulations in HDF5 (Hierarchical Data Format), i.e. Chombo, files. As far as we know Shape has not yet been developed to import HDF5. Currently it accepts ASCII files. Therefore the simulations you want to visualize will need to be in ASCII form.
Shape can read and visualize 3D simulations. However make sure that the columns in your ASCII file are converted to a 3D Cartesian coordinate system. You will need a minimum of 7 columns, of the form: px, py, pz, vx, vy, vz, n
. Where all pi
are the positions in (x, y, z)
, all vi
are the respective velocities, and n
is the density. One can also have a column for temperature T
, and pressure P
. Unlike VisIt, where you visualize the time evolution of a simulation frame-by-frame, Shape will only require one frame to be visualized. In Shape, you can rotate it and study the morphology at that instant. It may be best to visualize with VisIt first, pick your frame, and then convert it to ASCII. See the Before Starting section for more on this.
Remember that Shape's goal is to be a software that reduces restrictions on physical assumptions so that astronomers can make a realistic reconstruction of an astrophysical object's morphology. Thus the software is useful for defining 3D structural elements to create a model that can be optimized to fit an observation. Here, by importing external data, we are approaching from the opposite angle. We are using a code to generate a simulation of an object, or phenomena, that is comparable to actual astrophysical objects. In Shape we can then visualize, and compare these simulations more easily with observation. Expect the animations you make with Shape to not look as defined as those one could make in VisIt with the same simulation, as they are meant to look like emission maps.
Links
In this tutorial in visualizing external data in Shape, we will be visualizing a sphere generated by the following script:
Pseudo ASCII data generator (python)
import random from math import pi,sin,cos #Creating Sphere Dataz def createSphere(r=5, N=100): lst = [] for phi in [(pi*i)/(N-1) for i in range(N)]: M = int(sin(phi)*(N-1))+1 for theta in [(2*pi*i)/M for i in range(M)]: x = r * sin(phi) * cos(theta) y = r * sin(phi) * sin(theta) z = r * cos(phi) lst.append((x, y, z)) return lst #Opens/creates new ascii file outfile = open("test_sphere_7col.dat", "w") #Writes the data to the file for x,y,z in createSphere(): rho = random.random()*1000000 vx = random.random()*10 vy = random.random()*100 vz = random.random() print("{0}\t{1}\t{2}\t{3}\t{4}\t{5}\t{6}".format(x,y,z,vx,vy,vz,rho), file=outfile) #Closes the file outfile.close()
a file generated by the script can be found in the attachments, titled test_sphere_7col.dat. In the following sections you'll observe how we come to visualize its "emission map." We suggest downloading the attachment and following along with the tutorial.
Tutorial by Example
Before Starting:
Before Starting
- Make sure your simulation is in ASCII form. Hopefully you will find this guide on how to convert HDF5 to ASCII helpful.
- (Currently Subjected to Testing) The basic format you will want for your ASCII file is a series of 7 columns minimum. They will be of the form:
px,py,pz,vx,vy,vz,n
, the Cartesian coordinates of position, velocity and the density.
- Import the data into Excel. Take the averages of all of the columns. You can open the file in terminal and scroll to the bottom to determine how many lines there are. You will use these values later to center your data.
Figure 1. A screen shot of the Excel spreadsheet, where the pseudo sphere ASCII data was imported (test_sphere_7col.dat). The averages are bolded. Headers are above the averages, and the raw data below.
- Once you have looked at your data to see what format it is in (number of columns, delimiters), and gotten the averages for each column, you can open Shape.
- Create a new project: File → New.
- Click on the 3D Module, it looks like a 2x2x2 cube.
3D Module:
3D Module
When you first open the 3D Module, it will look like this:
Figure 2. A screen shot of the 3D Module once a new project is opened.
Now following the numbers in Figure 3:
Figure 3. A screen shot of the 3D Module after creating a mesh sphere.
In order to import external data, we will need to use one of the mesh objects listed at the top of the 3D Module (Figures 2 & 3) (Sphere, Torus, Cone, Cube, Plane, Cylinder, etc.) to carry it into the software. This is all done in this module, and then we visualize it in the Render Module. We start in the General tab.
- In this case we chose the sphere. Simply click the icon, then click in any of the four windows. Then click again to choose the size of the sphere. The size does not matter, as we will get rid of the sphere later on. We just need it to carry the data into Shape.
- Note that under the general tab we can change the name of the object. We do this here in Figure 4 (below), to External Sphere. You don't have to change it. But it might be useful in case you want to import other simulations in the same session.
- Choose the ASCII format, as that is the format our simulation is in. Then click the square button next to the drop down tab, to import the data in the Input Parameters dialog. This is what it looks like before we direct Shape to where our ASCII file is.
|
Figures 4 & 5. (left) The 3D module after opening the ASCII renderer, and having changed the name of the mesh object. (right) The input parameters for the ASCII simulation we are aiming to import into Shape.
On Input Properties:
Input
- Filename: The path used to reach your ASCII file is input here.
- dNum: Keep this as the default (
1E4
). - Format: As seen in Figure 4 the default is
px,py,pz,n,vx,vy,vz,T,P
. These are all the variable columns you may input to shape: position, density, velocity, temperature, and pressure. However in this case we are simply going to use the position, velocity coordinates, and density. Therefore the format will bepx,py,pz,vx,vy,vz,n
. - Comment: How comments are denoted in your ASCII file. Keep the default of
#
. - n: Keep this as the default (
1
). - Lines per entry: How many rows are in your ASCII file? Open the file in whatever text editor your prefer, and get the number for the last line. In this example there were 6284 rows, or lines per entry.
Column Separator
How is your file delimited? Typically when you input the filename Shape will adjust the column separator automatically. However check that it is correct.
Center
For the x,y,z
inputs, enter the averages from your Excel file for the position coordinates. If you click to detect extents after entering the averages, they should be 0,0,0
.
Preview
Allows you to preview your data file.
Now you have imported your data. However you will not be able to see it in the 3D Module yet. Now move onto the Modifiers tab.
Figure 6. Under the Modifiers tab in the 3D Module.
After importing the data under the General tab, we now move to the Modifiers tab. There you will see two modifiers (temperature, density). Delete these using the X button, while they are highlighted. A new window will appear asking if you want to delete these.
Now we will want to move into the Physics Module:
Image 7. In 3D Module directing you to the Physics Module for the purpose of creating a species for our emission map to be viewed in the Render Module.
Click the icon at the top of the Shape GUI that looks like a chalkboard. Now move on to the next section of this tutorial.
Physics Module:
Physics Module
In this section of the tutorial, we will be making a species for our emission map. To do this we will need to create a new custom species, and edit the emission coefficient so that it uses a Gaussian template. Depending on what sort of simulation you've made, you may want to do something different.
Image 8. Screen capture of the default physics module. Note the red circle, indicating where to create a new species.
- Click the Add button at the top of the Physics Module.
- After doing this, you will be prompted to add a certain type of species. There are many options (
Quick, Custom, Atomic, Molecule, Dust, Scale, Multiple, Plasma, Black Body
), however here we will choose Custom.
Image 9. Prompt for adding a species. Here we have selected a Custom species.
- Once you have clicked Ok the Physics Module will now give you the option to edit a Custom species. The module will look similar to what is presented in Image 10 (sans red text, of course).
Image 11. The Physics Module after creating a new Custom species. Note the only user change here is the Name of the species, which is "Line."
Now that you have created a new Custom species we need to edit it.
- Here the Name of the species has been changed to "Line." However one may choose whatever name they prefer. The same name is used in the youtube tutorial provided above.
- Note that once you choose a Custom species, the main view in the Physics Module may still be that of the Default. Simply switch from Default to Custom in the left-hand Species list by clicking the Custom option you have made.
In this tutorial we will simply keep the default options illustrated in Image 11. Here the Emission contribution is sufficient. Now we go on to edit the Emission Coefficient, which is indicated in Image 11.
- Click the Edit button.
- Also note that if you scroll down the main options page in Shape for your new Custom species, you'll see that we have
n
, or Density, as our variable mapping. - Recall that our format string for our imported data is in terms of the columns with headers
px,py,pz,vx,vy,vz,n
.
|
Image(s) 12, 13, & 14. The Emission Coefficient dialog.
Once you click Edit for your emission coefficient, you'll see a new dialog pop up as illustrated in Image 12. You are given the analytic function:
f(x) = 1E-15*n^2
We do not want to use this analytic function.
- So as illustrated in Image 12, click the Templates option and choose Gaussian. There are the options of (
Planck (m), Gaussian, Lorentzian, Wein
). - After clicking Apply, you will notice that the analytic function, shown under Functions has changed to:
[Analytic] 1.0/(sqrt(2*pi) * c) * e^(-(x-b)^2/(2*c^2))
- Notice that this equation does not involve our variable mapping term,
n}}, otherwise known as density. As illustrated in '''Image 13''' click back to the ''Function'' property, and delete the coefficient on our Gaussian function, {{{ 1.0/(sqrt(2*pi) * c)
and replace it withn
. - Then hit enter in the white box, so that it shows up under your Functions.
- Also note that there are other variables in our Gaussian function, namely,
x,c,b
. Clearlyx
is the variable of the function, and not a constant, where asc, b
are constants. As illustrated in Figure 14 input the following numbers:
x = 9.8 c = 5.0E-7 b = 2.0E-11
Leave n
alone. Note:
n
, the coefficient, is the height of the Gaussian's peak.c
is the standard deviation, and controls the spread of the peak. Here we have it spread at 5.0E-7 Angstroms, which is essentially corresponds to the energy of a photon in a vacuum (gamma rays). This is similarly the case withb
b
is the position of the center of the Gaussian peak.- Thus our emission is within the realm of visible light in the spectrum.
- We simply choose
x
based on what is recommended in the video tutorial.
- When you are done with the Emission Coefficient dialog, simply click the X-out button at the top of the dialog.
You are done. Now we can go back to the 3D Module, so proceed to the next section of this tutorial.
Image 15. The 3D Module is the 3x3x3 cube near the left-hand corner of the Shape GUI. See the red circle.
3D Module:
3D Module
Importing Species from Physics Module
Now that are we are back in the 3D Module we will need the species we have just made. So we need to import it into this module.
Image 16. The Species tab in the 3D Module.
- Go to the Species tab in the 3D Module.
- You should see no species are listed.
- Go to the drop down, and choose the new species you have made. In this tutorial we have named it "Line."
- Click Add.
Importing Particles
Render Module:
Render Module
Attachments (39)
-
SHAPE00000.dat
(5.9 MB
) - added by 10 years ago.
ascii data — frame 0, initial conditions
-
SHAPE00022.dat
(647.3 KB
) - added by 10 years ago.
ascii data — frame 22
- Screen Shot 2015-06-16 at 11.47.22.png (111.1 KB ) - added by 10 years ago.
-
test_sphere_7col.dat
(804.1 KB
) - added by 10 years ago.
ascii data — a sphere
- Screen Shot 2015-06-16 at 13.20.02.png (133.0 KB ) - added by 10 years ago.
- Screen Shot 2015-06-16 at 13.20.28.png (254.6 KB ) - added by 10 years ago.
- Screen Shot 2015-06-16 at 13.20.57.png (245.5 KB ) - added by 10 years ago.
- Screen Shot 2015-06-16 at 13.23.32.png (72.5 KB ) - added by 10 years ago.
- Screen Shot 2015-06-16 at 13.23.50.png (242.8 KB ) - added by 10 years ago.
- Screen Shot 2015-06-16 at 13.24.19.png (250.8 KB ) - added by 10 years ago.
- Screen Shot 2015-06-16 at 13.24.31.png (82.6 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.24.45.png (10.4 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.25.03.png (82.8 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.25.24.png (39.0 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.26.17.png (39.3 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.26.28.png (36.6 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.26.42.png (22.9 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.28.58.png (308.5 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 14.57.28.png (99.5 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.26.56-new.png (27.4 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 15.02.22.png (160.3 KB ) - added by 9 years ago.
- Screen Shot 2015-06-16 at 13.28.02-new.png (249.5 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 15.09.42.png (249.8 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 15.10.06.png (23.6 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 15.10.18.png (24.0 KB ) - added by 9 years ago.
-
3d.png
(5.1 KB
) - added by 9 years ago.
3d module icon
-
physics.png
(5.0 KB
) - added by 9 years ago.
physics module icon
- render.png (4.2 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 16.40.11.png (230.1 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 16.43.06.png (120.4 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.28.25.png (36.5 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.28.40.png (35.9 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.28.53.png (28.3 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.29.02.png (23.6 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.29.10.png (34.1 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.29.37.png (22.8 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.29.45.png (24.7 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.29.54.png (6.6 KB ) - added by 9 years ago.
- Screen Shot 2015-07-13 at 17.30.03.png (153.5 KB ) - added by 9 years ago.