Presently I am attempting sequential gaussian simulation using SGeMS software. I have followed the steps as shown in various videos available online and also the user guide (Remy et. al., 2006). I have set a grid of 50x20x100 as per procedure and finally derived 50 realisations. Now, I want to retrieve the simulated data along with the grid point locations so that I can plot them in some different software.
I followed the steps of "object -> save project -> grid to save -> property to save -> csv format" and saved the file. However, when i open it, I get the simulated values of variable but not its XYZ coordinates. I have attempted multiple times unsuccessfully.
Can any one point out where am I missing the step?
TIA.
Greetings,
You haven't missed any steps. The simulation results are delivered by SGEMS in a file with as many rows as there are cells in your grid. In your case, I believe you have a file with 1,000,000 rows of different values.
If you want to plot it in another software, reshape the result array to the shape of your computation grid, i.e., from (1 x 1,000,000) to (50, 20, 100).
It's simple to do in Python with Numpy's 'reshape' method.
Make sure the reshape operation doesn't flip your simulation.
You should compare your plot to the SGEMS's built-in visualization tool.
If you notice that the results are flipped, you can use Numpy's 'flipud' method to fix it.
If you truly require the X, Y, and Z coordinates, you must define them yourself. It is simple to do because the grid has the same dimensions (dx, dy, dz) in each direction. If you're familiar with Python, some guidelines and code snippets can be found at https://github.com/robinthibaut/pysgems.
For example, in Python, you could define an empty Numpy array of shape (50, 20, 100, 3). The `3` in parentheses represents the three spatial dimensions. You can then replace each point coordinate by your `true` spatial coordinates. The function `blocks_from_rc` in the pysgems package does this for you.
I hope this helps!
Robin
Thank you Dr Robin for your prompt answer. However, my problem still persists. The reasons are:
1. My cell sizes are 2mx25mx5m in XYZ direction, hence even if I reshape my computation grid of 1x100000 to 50x20x100, I wont be able to assign the right XYZ coordinates to each centroid unless I know the sequence of numbering of centroids.
2. I don't know use of Python.
I wonder whether/why the software which takes input of cordinates in proper format with clearly defined X,Y,Z coordinates, does not deliver the output in similar fashion.
Thank you for your help anyway.
Best wishes
Ajoy Kumar Padhi
Thank you very much. I will, God willing, follow the steps and thoroughly search for them. Greetings.
- Badges
- Science topic
- Similar topics
- Geoscience
- Cartography
- Mapping