diff --git a/Bericht/vtk/result_angles.vtk b/Bericht/vtk/result_angles.vtk index f0c5312a51d8a7645bef29873bb3f05a153e0734..0d7fdc3f3f39d1c8cbb7b083ec79bc56c07da473 100644 Binary files a/Bericht/vtk/result_angles.vtk and b/Bericht/vtk/result_angles.vtk differ diff --git a/UNet/visualize_vtk.ipynb b/UNet/visualize_vtk.ipynb index 708cbf71b122d020eaa813493d1ab4f54a0917a4..87a06a107015b321055262c07312989bb071a5f7 100644 --- a/UNet/visualize_vtk.ipynb +++ b/UNet/visualize_vtk.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 31, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -20,7 +20,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -106,7 +106,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 8, "metadata": {}, "outputs": [ { @@ -122,8 +122,8 @@ } ], "source": [ - "path_to_UNET = 'F:/RWTH/HiWi_IEHK/DAMASK3'\n", - "export_path = 'F:/RWTH/HiWi_IEHK/DAMASK3/Bericht/vtk/result_angles'\n", + "path_to_UNET = 'E:/Data/damask3'\n", + "export_path = 'E:/Data/damask3/Bericht/vtk/result_angles'\n", "UNet = UNet\n", "Training_data_32 = torch.load(f'{path_to_UNET}/UNet/Input/TD_norm_32_angles.pt')\n", "normalization_32 = np.load(f'{path_to_UNET}/UNet/Input/Norm_min_max_32_phase.npy', allow_pickle=True)\n", @@ -142,16 +142,67 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": 11, "metadata": {}, "outputs": [], "source": [ - "grid= pv.read('F:/RWTH/HiWi_IEHK/DAMASK3/Bericht/vtk/result_angles.vtk')\n", + "grid= pv.read('E:/Data/damask3/Bericht/vtk/result_angles.vtk')\n", "\n", "slice= grid.slice()\n", "pv.global_theme.font.color = 'black'\n", "pv.global_theme.font.size = 50 \n", - "slice['stress'] = slice['stress']/1000000000\n" + "slice['stress'] = slice['stress']/1000000000\n", + "slice['label'] = slice['label']/1000000000\n" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "<table><tr><th>Header</th><th>Data Arrays</th></tr><tr><td>\n", + "<table>\n", + "<tr><th>UniformGrid</th><th>Information</th></tr>\n", + "<tr><td>N Cells</td><td>32768</td></tr>\n", + "<tr><td>N Points</td><td>35937</td></tr>\n", + "<tr><td>X Bounds</td><td>0.000e+00, 3.200e+01</td></tr>\n", + "<tr><td>Y Bounds</td><td>0.000e+00, 3.200e+01</td></tr>\n", + "<tr><td>Z Bounds</td><td>0.000e+00, 3.200e+01</td></tr>\n", + "<tr><td>Dimensions</td><td>33, 33, 33</td></tr>\n", + "<tr><td>Spacing</td><td>1.000e+00, 1.000e+00, 1.000e+00</td></tr>\n", + "<tr><td>N Arrays</td><td>1</td></tr>\n", + "</table>\n", + "\n", + "</td><td>\n", + "<table>\n", + "<tr><th>Name</th><th>Field</th><th>Type</th><th>N Comp</th><th>Min</th><th>Max</th></tr>\n", + "<tr><td><b>label</b></td><td>Cells</td><td>float64</td><td>1</td><td>-1.215e+308</td><td>3.347e+305</td></tr>\n", + "</table>\n", + "\n", + "</td></tr> </table>" + ], + "text/plain": [ + "UniformGrid (0x20607430d60)\n", + " N Cells:\t32768\n", + " N Points:\t35937\n", + " X Bounds:\t0.000e+00, 3.200e+01\n", + " Y Bounds:\t0.000e+00, 3.200e+01\n", + " Z Bounds:\t0.000e+00, 3.200e+01\n", + " Dimensions:\t33, 33, 33\n", + " Spacing:\t1.000e+00, 1.000e+00, 1.000e+00\n", + " N Arrays:\t1" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "grid" ] }, { @@ -220,19 +271,20 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], "source": [ - "%%capture\n", - "p = pv.Plotter(off_screen=True)\n", - "sargs = dict(\n", - " title='Stress DAMASK (GPa)',\n", - " n_labels=2,\n", - ")\n", - "p.add_mesh(slice,scalars=\"label\",scalar_bar_args=sargs,name='Stress DAMASK (GPa)')\n", - "p.view_zy()\n", - "p.screenshot('F:/RWTH/HiWi_IEHK/DAMASK3/Bericht/vtk/stress_DAMASK.png',transparent_background=True)" + "plt.figure(figsize=(12,12))\n", + "data = slice['label'].reshape(32,32)\n", + "\n", + "data = np.flip(data,1)\n", + "data = np.rot90(data,1)\n", + "plt.imshow(data,cmap='viridis',origin='lower',vmin=0.9,vmax=3.2)\n", + "\n", + "plt.colorbar(orientation=\"horizontal\")\n", + "plt.show()\n", + "#plt.savefig('F:/RWTH/HiWi_IEHK/DAMASK3/Bericht/vtk/difference.png',transparent_background=True)" ] }, { @@ -272,28 +324,17 @@ "grid.spacing = (1,1,1)\n", "grid.cell_data[\"difference\"] = difference.flatten(order = \"F\")\n", "slice= grid.slice()\n", - "#data = slice['difference'].reshape(32,32)\n", - "#annotations_difference = {\n", - "# 0: '0%'\n", - "#}\n", "\n", - "#p = pv.Plotter(off_screen=False)\n", - "#sargs_difference= dict(\n", - "# title='Difference',\n", - "# n_labels = 2\n", - "# )\n", - "#grid.plot(scalars=\"error\",notebook=False,cmap=colormap_error)\n", - "#p.add_mesh(slice,scalars=\"difference\",scalar_bar_args = sargs_difference, annotations=annotations_difference,cmap='coolwarm', norm=divnorm)\n", "plt.figure(figsize=(12,12))\n", "data = slice['difference'].reshape(32,32)\n", "data = np.flip(data,1)\n", "data = np.rot90(data,1)\n", "plt.imshow(data,norm=divnorm,cmap='coolwarm',origin='lower')\n", - "#plt.tick_params(axis=’x’, which=’both’, bottom=False, top=False, labelbottom=False)\n", + "\n", "plt.colorbar(orientation=\"horizontal\")\n", "#plt.show()\n", "plt.savefig('F:/RWTH/HiWi_IEHK/DAMASK3/Bericht/vtk/difference.png',transparent_background=True)\n", - "#p.screenshot('F:/RWTH/HiWi_IEHK/DAMASK3/Bericht/vtk/difference.png',transparent_background=True)\n" + "\n" ] }, { @@ -350,7 +391,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.9.10" + "version": "3.9.5" }, "orig_nbformat": 4 },