## InSituNet: Deep Image Synthesis for Parameter Space Exploration of Ensemble Simulations

Wenbin He, Junpeng Wang, Hanqi Guo, Member, IEEE , Ko-Chih Wang, Han-Wei Shen, Member, IEEE , Mukund Raj, Youssef S. G. Nashed, and Tom Peterka, Member, IEEE

Abstract -We propose InSituNet, a deep learning based surrogate model to support parameter space exploration for ensemble simulations that are visualized in situ. In situ visualization, generating visualizations at simulation time, is becoming prevalent in handling large-scale simulations because of the I/O and storage constraints. However, in situ visualization approaches limit the flexibility of post-hoc exploration because the raw simulation data are no longer available. Although multiple image-based approaches have been proposed to mitigate this limitation, those approaches lack the ability to explore the simulation parameters. Our approach allows flexible exploration of parameter space for large-scale ensemble simulations by taking advantage of the recent advances in deep learning. Specifically, we design InSituNet as a convolutional regression model to learn the mapping from the simulation and visualization parameters to the visualization results. With the trained model, users can generate new images for different simulation parameters under various visualization settings, which enables in-depth analysis of the underlying ensemble simulations. We demonstrate the effectiveness of InSituNet in combustion, cosmology, and ocean simulations through quantitative and qualitative evaluations.

Index Terms -In situ visualization, ensemble visualization, parameter space exploration, deep learning, image synthesis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-16c5c2eb-a0f2-4e3d-9b85-1159eecd0002.jpg)

## 1 INTRODUCTION

Ensemble simulations [64] have been playing an increasingly important role in various scientific and engineering disciplines, such as computational fluid dynamics, cosmology, and weather research. As the computational power of modern supercomputers continues to grow, ensemble simulations are more often conducted with a large number of parameter settings in high spatial and/or temporal resolutions. Despite the advances in accuracy and reliability of simulation results, however, two challenges have emerged: (1) I/O bottleneck for the movement of the large-scale simulation data and (2) effective exploration and analysis of the simulation parameters. In situ visualization [5, 38], which generates visualization at simulation time and stores only the visualization results (that are much smaller than the raw simulation data [1, 2]) for post-hoc analysis, addresses the first challenge to some extent. However, it also limits the flexibility of post-hoc exploration and analysis, because the raw simulation data are no long available.

This study focuses on improving scientists' ability in exploring the in situ visualization results of ensemble simulations and extending their capability in investigating the influence of different simulation parameters. Several pioneering works have been proposed to facilitate post-hoc exploration of in situ visualization results. For example, the Cinema framework [1,2] visualized the simulation data from different viewpoints in situ and collected images to support post-hoc exploration. The volumetric depth images [21,22] stored ray segments with composited color and opacity values to enable post-hoc exploration of arbitrary viewpoints for volume rendering. However, these approaches focus more on extending the capability to explore the visual mapping parameters (e.g., transfer functions) and view parameters (e.g., view angles) and have little consideration of the simulation parameters, which are important in studying ensemble simulations.

Simulation parameter space exploration is not trivial, because the relationship between the simulation parameters and outputs is often

- · Wenbin He, Junpeng Wang, Ko-Chih Wang, and Han-Wei Shen are with the Department of Computer Science and Engineering, The Ohio State University. E-mail: { he.495, wang.7665, wang.3182, shen.94 } @osu.edu.
- · Hanqi Guo, Mukund Raj, Youssef S. G. Nashed, and Tom Peterka are with the Mathematics and Computer Science Division, Argonne National Laboratory. E-mail: { hguo, mraj, ynashed, tpeterka } @anl.gov.

Manuscript received xx xxx. 201x; accepted xx xxx. 201x. Date of Publication xx xxx. 201x; date of current version xx xxx. 201x. For information on obtaining reprints of this article, please send e-mail to: reprints@ieee.org. Digital Object Identifier: xx.xxxx/TVCG.201x.xxxxxxx highly complex. The majority of existing simulation parameter space exploration approaches [14, 65] resorted to visualizing a set of simulation parameters and outputs simultaneously and revealing the correspondence between the parameters and outputs through visual linkings. However, these approaches often depend on the raw simulation data that might not be available for large-scale ensemble simulations. Moreover, these approaches have limited ability in inferring simulation outputs with respect to new parameters. Hence, extra simulations have to be conducted for new parameters, which cost enormous computational resources for most scientific simulations.

In this work, we propose InSituNet, a deep learning based surrogate model to support parameter space exploration for ensemble simulations that are visualized in situ. Our work is based on the observation that images of high accuracy and fidelity can be generated with deep neural networks for various image synthesis applications, such as superresolution [18, 32, 36], inpainting [46, 68], texture synthesis [24, 71], and rendering [6, 19]. Specifically, we train InSituNet to learn the end-to-end mapping from the simulation, visual mapping, and view parameters to visualization images. The trained model enables scientists to interactively explore synthesized visualization images for different simulation parameters under various visualization settings without actually executing the expensive simulations. Our approach consists of three major steps.

- 1. In situ training data collection from ensemble simulations Given ensemble simulations conducted with different simulation parameters, we visualize the generated simulation data in situ with various visual mapping and view parameters. The resulting visualization images and the corresponding parameters are collected and used for the offline training of InSituNet.
- 2. Offline training of InSituNet Given the parameters and image pairs, we train InSituNet (i.e., a convolutional regression model) with cutting-edge deep learning techniques on image synthesis to map simulation, visual mapping, and view parameters to visualization images directly.
- 3. Interactive post-hoc exploration and analysis With the trained InSituNet, we build an interactive visual interface that enables scientists to explore and analyze the simulation from two perspectives: (1) inferring visualization results for arbitrary parameter settings within the parameter space with InSituNet's forward propagations and (2) analyzing the sensitivity of different parameters with InSituNet's backward propagations.

We demonstrate the effectiveness of the proposed approach in combustion, cosmology, and ocean simulations, and compare the predicted images of InSituNet with the ground truth and alternative methods. In

addition, we evaluate the influence of different hyperparameters of InSituNet (e.g., the choice of loss functions and the network architectures) and provide guidance in configuring the hyperparameters. In summary, the contributions of this paper are threefold:

- · A deep image synthesis model (i.e., InSituNet) that enables posthoc parameter space exploration of ensemble simulations
- · An interactive visual interface to explore and analyze the parameters of ensemble simulations with the trained InSituNet
- · A comprehensive study revealing the effects of different hyperparameters of InSituNet and providing guidance for applying InSituNet to other simulations

## 2 RELATED WORK

In this section, we review related work in image-based in situ visualization, parameter space exploration of ensemble simulations, and deep learning for visualization.

## 2.1 Image-Based In Situ Visualization

Based on the output, in situ visualization can be categorized into imagebased [1,2], distribution-based [20], compression-based [17,35], and feature-based [12] approaches. We regard our work as an image-based approach, which visualizes simulation data in situ and stores images for post-hoc analysis. Tikhonova et al. [59-61] generated images of multiple layers in situ to enable the adjustment of transfer functions in post-hoc analysis. Frey et al. [22] proposed volumetric depth images, a compact representation of volumetric data that can be rendered efficiently with arbitrary viewpoints. Fernandes et al. [21] later extended volumetric depth images to handle time-varying volumetric data. Biedert and Garth [8] combined topology analysis and imagebased data representation to preserve flexibility for post-hoc exploration and analysis. Ahrens et al. [1, 2] proposed Cinema, a framework that stores visualization images in situ and performs post-hoc analysis via exploration and composition of those images.

Compared with these approaches, our work supports not only the exploration of various visual mapping and view parameters but also the creation of visualizations under new simulation parameters without actually running the simulation.

## 2.2 Parameter Space Exploration

The existing parameter space exploration works for ensemble simulations can generally be reviewed from two perspectives, the adopted visualization techniques and the objective of parameter space exploration. Visualization techniques that designed for high-dimensional data are often borrowed to visualize the parameter space of ensemble simulations, as the simulation parameters are typically treated as multidimensional vectors. These techniques include but are not limited to: parallel coordinate plots [44, 65], radial plots [14-16], scatter plots [39, 45, 58], line charts [9], matrices [48], and glyphs [10]. For the objectives of parameter space exploration, we believe the six tasks sorted out by Sedlmair et al. [55] could best summarize the literature, which are optimization [62], partitioning [7, 65], filtering [47], outliers [47, 49], uncertainty [9, 11], and sensitivity [9]. We refer the interested readers to the work of Sedlmair et al. [55] for the detailed definition of each task, as well as the example visualization works.

The aforementioned parameter visualization techniques and analysis tasks mostly focus on a limited number of simulation inputs and outputs collected from ensemble runs. In this paper, we train a surrogate model to extend our study to arbitrary parameter settings within the parameter space, even if the simulations were not executed with those settings. In addition, our approach is incorporated with in situ visualization, which is widely used in large-scale ensemble simulations.

## 2.3 Deep Learning for Visualization

The visualization community has started to incorporate deep learning in visualization research. For example, Hong et al. [30] used long short-term memory [29] to estimate access pattern for parallel particle tracing. Han et al. [26] used autoencoders [52] to cluster streamlines and streamsurfaces. Xie et al. [67] used neural network embeddings to detect anomalous executions in high performance computing applications. Berger et al. [6] proposed a deep learning approach to assist transfer function design using generative adversarial networks (GANs) [25], which is closely related to our approach. Specifically, we focus on parameter space exploration of ensemble simulations instead of transfer function design for volume rendering.

Our work is related to deep learning based image synthesis, which has been used in various applications, including super-resolution [18, 32,36], denoising [68,70], inpainting [46,68], texture synthesis [24,71], text-to-image synthesis [50], style transfer [23, 32, 72], and rendering [6, 19]. We investigate and combine different state-of-the-art deep learning techniques on image synthesis (e.g., perpetual losses [32, 36] and GANs [25]) to improve the quality of our image synthesis results.

## 3 OVERVIEW

Fig. 1. Workflow of our approach. Ensemble simulations are conducted with different simulation parameters on supercomputers, and visualization images are generated in situ for different visual mapping and view parameters. The generated images and the parameters are collected into an image database. A deep image synthesis model (i.e., InSituNet) is then trained offline based on the collected data, which is later used for parameter space exploration through an interactive visual interface.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-ab579852-3df6-4a19-a34e-2d74964d3c06.jpg)

Figure 1 provides the workflow of our approach, which consists of three major components. First, given ensemble simulations conducted with different simulation parameters, we visualize the generated simulation outputs in situ with different visual mapping and view parameters on supercomputers. The three groups of parameters-simulation, visual mapping, and view parameters-along with the corresponding visualization results (i.e., images) are collected to constitute an image database (Section 4). Second, with the collected data pairs between parameters and the corresponding images, we train InSituNet to learn the end-to-end mapping from the simulation inputs to the visualization outputs (Section 5). To improve the accuracy and fidelity of the generated images, we use and combine different state-of-the-art deep learning techniques on image synthesis. Third, with the trained InSituNet, we build an interactive visual interface (Section 6) to explore and analyze the parameters from two aspects: (1) predicting visualization images interactively for arbitrary simulation, visual mapping, and view parameters within the parameter space and (2) investigating the sensitivity of different input parameters to the visualization results.

## 4 IN SITU TRAINING DATA COLLECTION

Fig. 2. Our in situ training data collection pipeline. Simulation data, generated with different simulation parameters, are visualized in situ with different visual mapping and view parameters. The in situ visualization generates a large number of images, which are collected along with the corresponding parameters for the training of InSituNet offline.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-7cbe04d9-3725-4398-be6d-8021e99bc7e3.jpg)

Figure 2 illustrates our in situ training data collection pipeline. Given ensemble simulations conducted with different simulation parameters, we perform in situ visualization with a desired set of visual mapping parameters (e.g., isosurfaces extraction with a set of isovalues) and different view parameters (e.g., viewpoints). We denote an instance of simulation, visual mapping, and view parameters as Psim , Pvis , and Pview , respectively, which corresponds to a visualization image I . The parameters (highlighted in green in Figure 2) and the corresponding visualization images (highlighted in blue in Figure 2) constitute data pairs, which will be stored and used to train InSituNet. InSituNet

learns a function F that maps the three groups of parameters to the corresponding visualization image, which can be defined as

$$\mathcal { F } ( P _ { s i m }, P _ { v i s }, P _ { v i e w } ) \to I, \quad \quad \quad ( 1 ) \quad \text{featu}$$

so that it can predict visualization images for unseen parameters. In the following, we discuss the three groups of parameters in detail.

Simulation parameters Psim are represented as a vector with one or more dimensions, and the value range of each dimension is defined by scientists. By sweeping the parameters within the defined ranges, ensemble simulations are conducted to generate the ensemble data.

Visual mapping parameters Pvis are predefined operations to visualize the generated simulation data, such as pseudo-coloring with predefined color schemes. Note that we limit the users' ability in selecting arbitrary visual mappings to produce and store fewer images.

View parameters Pview are used to control the viewpoints that the images are created from. In this work, we define the viewpoints by a camera rotating around the simulation data, which is controlled by azimuth q ∈ [ 0 360 , ] and elevation f ∈ -[ 90 90 . , ] For panning and zooming, we resort to image-based operations (i.e., panning and resizing the images) as proposed in [2]). To train a deep learning model that can predict visualization images for arbitrary viewpoints, we sample the azimuth and elevation and generate images from the sampled viewpoints. Based on our study, we found that taking 100 viewpoints for each ensemble member is sufficient to train InSituNet.

With the specified values for the three groups of parameters, we generate the corresponding visualization images. Our work uses RGB images compressed to the portable network graphics (PNG) format instead of more sophisticated image formats, such as volumetric depth images [21,22] or explorable images [59-61], for two reasons. First, the benefits of using those sophisticated image formats, such as supporting changing of viewpoints, can be achieved by InSituNet trained on the RGB images. Second, RGB images are more generally applicable for various visualizations and more easily to be handled by neural networks compared with those sophisticated image formats.

## 5 INSITUNET ARCHITECTURE AND TRAINING

Fig. 3. Overview of InSituNet, which is a convolutional regression model that predicts visualization images from input parameters. During training, the regression model is trained based on the losses computed with the assist of a pretrained feature comparator and a discriminator.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-f6aa0402-8093-43d5-ab8f-f59af4e6d3c8.jpg)

Figure 3 illustrates the training and inference pipelines of InSituNet. In the training stage, InSituNet consists of three subnetworks: a regressor, a feature comparator, and a discriminator. The regressor R w is a deep neural network (defined by a set of weights w ) modeling the function that maps input parameters to visualization images as defined in Equation 1. To train a regressor that can generate images of high fidelity and accuracy, we introduced the feature comparator F and the discriminator D u to compute losses by comparing the predicted and the ground truth images. The feature comparator is a pretrained neural network whose convolutional kernels are used to extract and compare image features (e.g., edges, shapes) between the predicted and the ground truth images to obtain a feature reconstruction loss. The discriminator D u is a deep neural network whose weights u are updated during training to estimate the divergence between the distributions of the predicted and the ground truth images. The divergence is known as the adversarial loss [25], which is combined with the feature reconstruction loss to train R w . In the inference stage, we need only the trained R w , which can predict visualization images for parameters that are not in the training data. In the following, we discuss the network architecture, the loss function, and the training process in detail.

## 5.1 Network Architecture

Three subnetworks are involved during training: the regressor R w , feature comparator F , and discriminator D u . The regressor R w and discriminator D u are two deep residual convolutional neural networks [27] parameterized by the weights w and u , respectively. The architectures of R w and D u are designed by following the network architecture proposed by [34,41], because the scale of our image synthesis problem is similar to theirs. For the feature comparator F , we use the pretrained VGG-19 model [57], which has been widely used in many deep image synthesis approaches [31,32,36].

## 5.1.1 Regressor R w

Fig. 4. Architecture of R w , which encodes input parameters into a latent vector with fully connected layers and maps the latent vector into an output image with residual blocks. The size of R w is defined by k , which controls the number of convolutional kernels in the intermediate layers.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-58ebbf63-5525-4a22-b2fa-5b01b5181199.jpg)

The architecture of R w is shown in Figure 4, which takes the Psim , Pvis , and Pview as inputs and outputs a predicted image I . The three types of parameters are first fed into three groups of fully connected layers separately, and the outputs are then concatenated and fed into another fully connected layer to encode them into a latent vector. Note that the parameters could also be concatenated first and then fed into fully connected layers. However, as each parameter is fully connected with all neurons in the next layer, more weights will be introduced in the network and the network size will increase. Next, the latent vector is reshaped into a low-resolution image, which is mapped to a high-resolution output image through residual blocks performing 2D convolutions and upsamplings. Following the commonly used architecture [34, 41], we use the rectified linear unit ( ReLU ) activation function [43] in all layers except the output layer. For the output layer, we use the tanh function to normalize each pixel into [ -1 1 . , ]

Note that we introduce a constant k in the network architecture to control the number of convolutional kernels in the intermediate layers. The constant k is used to balance the expressive power and the size and training time of R w to cope with datasets in different complexities.

Residual Blocks R w consists of several residual blocks (Figure 4a), which are proposed in [27] to improve the performance of neural networks with increasing depth. We adopted the residual blocks here because R w often needs to be very deep (i.e., more than 10 convolutional layers) to synthesize images with high-resolutions. Inside each residual block, the input image is first upsampled by using nearest neighbor upsampling. The upsampled image is then fed into two convolutional layers with kernel size 3 × 3. In the end, the original input image is added to the output, and the result is sent to the next layer. Batch normalizations are performed on the output of each convolutional layer to stabilize the training. Note that if the resolution or the channel number of the input image is not the same as the output, we perform the upsampling and convolution operations on the input image to transform it into the size of the output.

## 5.1.2 Discriminator D u

The architecture of D u is shown in Figure 5, which takes a predicted/ground truth image and the corresponding parameters as inputs and produces a likelihood value indicating how likely the input image is a ground truth image conditioning on the given parameters. With

Fig. 5. Architecture of D u . Input parameters and the predicted/ground truth image are transformed into latent vectors with fully connected layers and residual blocks, respectively. The latent vectors are then incorporated by using the projection-based method [42] to predict how likely the image is a ground truth image conditioning on the given parameters. Similar to R w , the size of D u is controlled by the constant k .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-e6a895c7-4983-4714-a96d-1aecc83908d5.jpg)

the likelihood values, an adversarial loss can be defined to update D u and R w (details in Section 5.2.2). Similar to R w , the three types of parameters are encoded into a latent vector in D u through fully connected layers. Meanwhile, the input image is fed through several residual blocks to derive its intermediate representation that is a latent vector. The latent vectors are then incorporated to obtain the likelihood value conditioning on the three groups of parameters. ReLU activations are used in all layers except the output layer, which instead uses the sigmoid function to derive a likelihood value within [ 0 1 . , ]

Residual blocks The architecture of the residual blocks in D u (Figure 5a) is similar to that in R w except that downsampling (average pooling in this work) is performed instead of upsampling to transform images into low-resolution representations and no batch normalization is performed, because it often hurts the performance of D u [34,37].

Projection-based condition incorporation We employed the projection-based method [42] to incorporate the conditional information (i.e., the three groups of parameters) with the input image. This method computes a dot product between the data to be incorporated, which in our work is the latent vector of the input parameters and the latent vector of the image (Figure 5b). Compared with other condition incorporation methods, such as vector concatenation [6, 40], the projection-based method improves the quality of conditional image synthesis results, as demonstrated by Miyato and Koyama [42].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-3fdbc89f-bf4a-4a11-8d87-8c81b5684334.jpg)

feature maps

Fig. 6. Architecture of F (i.e., VGG-19 network), where each layer is labeled with its name. Feature maps are extracted through convolutional layers (e.g., relu1 2 ) for feature-level comparisons.

To produce high quality image synthesis results, we also strive to minimize the feature-level difference between the generated and the ground truth image by employing a commonly used feature comparator F , namely the pretrained VGG-19 model [57] shown in Figure 6. F is a convolutional neural network, and the convolutional kernels on each layer have been pretrained to extract certain types of image features, such as edges and shapes. With it, we extract the features from a generated image, as well as its corresponding ground truth image, and minimize the difference between those features to improve the quality of the generated image (see details in Section 5.2.1). Specifically, we use the layer relu1 2 to extract feature maps for comparison based on two observations. First, early layers such as the layer relu1 2 of the VGG-19 network focus on low-level features such as edges and basic shapes, which commonly exist in scientific visualization images. Second, through our experiments we found that artifacts are introduced into the generated image by pooling layers (e.g., pool1 in Figure 6). Hence, we use the layer relu1 2 that is before the first pooling layer.

## 5.2 Loss Function

Given an image ˆ generated using I R w and the corresponding ground truth image I , a loss function L is defined by measuring the difference between them. Minimizing L can, therefore, be conducted by updating the parameters w of R w over an iterative training process. The most straightforward choice for L is the average of the pixel wise distance between ˆ and I I , such as the mean squared error. As shown in earlier works [31, 32, 36], however, the average pixel wise distance often produces over-smoothed images, lacking high-frequency features.

In this work, we define L by combining two advanced loss functions: a feature reconstruction loss [32] L F l , f eat and an adversarial loss [25] L adv R , namely

$$\mathcal { L } = \mathcal { L } _ { f e a t } ^ { F, l } + \lambda \mathcal { L } _ { a d v, R },$$

where l is the coefficient between them. L F l , f eat measures the difference between features extracted from the feature comparator F using its convolutional layer l , whereas L adv R quantifies how easily the discriminator D u can differentiate the generated images from real ones. As can be seen, minimizing L adv R requires training R w and D u together in an adversarial manner (i.e., the adversarial theory of GANs [25]). In order to train D u , an adversarial loss L adv D is used.

## 5.2.1 Feature Reconstruction Loss

The feature reconstruction loss between image ˆ and I I is defined by measuring the difference between their extracted features [31,32,36]. Specifically, for a given image I , our feature comparator F (the pretrained VGG-19) is applied on I and extracts a set of feature maps, denoted as F l ( ) I . Here l indicates which layer the feature maps are from (e.g., the relu1 2 of F ). The extracted feature maps can be considered as a 3D matrix of dimension h × × w c , where h , w , and c are the height, width, and number of channels, respectively, as shown in Figure 6. The feature reconstruction loss between ˆ and I I can, therefore, be defined as the pixel wise mean squared error between F l ( ) ˆ I and F l ( ) I . Extending this definition to a batch of images, the feature reconstruction loss between ˆ I 0: b -1 and I 0: b -1 ( b is the batch size) is

$$\mathcal { L } _ { f e a t } ^ { F, l } = \frac { 1 } { h w c b } \sum _ { i = 0 } ^ { b - 1 } \| F ^ { l } ( I _ { i } ) - F ^ { l } ( \hat { I } _ { i } ) \| _ { 2 } ^ { 2 }. \quad \quad ( 3 )$$

Using the feature reconstruction loss enables our regressor to produce images sharing similar feature maps with the corresponding ground truth images, which lead to images with sharper features.

## 5.2.2 Adversarial Loss

In addition to the feature reconstruction loss described above, we add an adversarial loss L adv R into the loss function. Unlike the feature reconstruction loss, which measures the difference between each pair of images, the adversarial loss focuses on identifying and minimizing the divergence between two image distributions following the adversarial theory of GANs. Specifically, our discriminator D u is trained along with the regressor R w to differentiate images generated by R w with ground truth images. As the regressor R w becomes stronger over the training, the discriminator D u is forced to identify more subtle differences between the generated images and the ground truth.

The adversarial loss can be used as complementary to the feature reconstruction loss for two reasons. First, the feature reconstruction loss focuses on the average difference between images, and the adversarial loss focuses on local features that are the most important to

differentiate the predicted and ground truth images. Second, the feature reconstruction loss compares the difference between each pair of the generated and ground truth images, and the adversarial loss measures divergence between two image distributions.

In this work, we use the standard adversarial loss presented in [25], which uses different loss functions for the generator and discriminator. For the generator (i.e., our regressor R w ), the adversarial loss is

$$\mathcal { L } _ { a d v. R } = - \frac { 1 } { b } \sum _ { i = 0 } ^ { b - 1 } \log D _ { \upsilon } ( \hat { I } _ { i } ), \quad \ \ ( 4 ) \quad \ \ 2 \colon \tilde { \mathbf F } \\. \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

which reaches the minimum when the discriminator cannot differentiate the generated images from the ground truth images. This loss is combined with the feature reconstruction loss to update our regressor (Equation 2). The adversarial loss of the discriminator is defined as

$$\mathcal { L } _ { a d v D } = - \frac { 1 } { b } \sum _ { i = 0 } ^ { b - 1 } ( \log D _ { \upsilon } ( I _ { i } ) + \log ( 1 - D _ { \upsilon } ( \hat { I } _ { i } ) ) ), \quad ( 5 ) \ \text{Algorit} \quad \text{param} \quad \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$$$

which estimates the divergence between the distribution of the generated images and the ground truth images.

## 5.3 Techniques to Stabilize Training

We use several techniques to stabilize the adversarial training of R w and D u . The instability of adversarial trainings is a well-known problem [25], especially when the resolution of synthesized images is high [6]. The previous work [6] divided the training into two stages for stabilization. In the first stage, the opacity GAN that produces 64 × 64 opacity images is trained, whereas the opacity-to-color translation GAN is trained in the second stage to produce 256 × 256 color images, conditioning on the 64 × 64 opacity images. In this work, we train a single pair of adversarial networks (i.e., R w and D u ) that directly produces 256 × 256 color images with the help of recent techniques in stabilizing the adversarial training, including the spectral normalization [41] and the two time-scale update rule (TTUR) [28].

## 5.3.1 Spectral Normalization

Spectral normalization [41] is used to mitigate the instability of the discriminator, which is a major challenge in stabilizing the adversarial training. Spectral normalization is a weight normalization technique, which outperforms other weight normalization techniques in many image synthesis tasks as shown in [34]. Spectral normalization normalizes the weight matrix of each layer based on the first singular value of the matrix. With spectral normalization, the discriminator is enforced to be Lipschitz continuous, such that the discriminator is constrained and stabilized to some extent. Spectral normalization is applied on each layer of the discriminator without changing the network architecture; hence spectral normalization is not labeled in Figure 5.

## 5.3.2 Learning Rate

The learning rates of R w and D u are critical for the stability of the adversarial training. This work uses the Adam optimizer [33] that changes the learning rate of each weight dynamically during training with respect to the momentum of the weight gradients. In detail, the learning rate in the Adam optimizer is controlled by three hyperparmeters: the initial learning rate a , the first-order momentum b 1 , and the second-order momentum b 2 . To stabilize the training, a small a is often preferred; and we found that 5 × 10 -5 stabilized the training in our cases. In addition, we found that a bigger b 1 often cripples the training and set b 1 to 0 as suggested in [13, 69]. Compared with b 1 , b 2 has less influence on the stability of the training, which is set to 0 999. .

In previous works on training GANs, we found that people often update the discriminator more frequently than the generator, because they do not want to update the generator based on a discriminator that is not strong enough. Doing so, however, leads to a longer training time. Our work uses the same update frequency for the regressor and discriminator but with different learning rates a D and a R (i.e., the TTUR technique [28]). Based on the empirical results shown in [13,69], we set the learning rate of the discriminator to be 4 times that of the regressor, that is, a D = 2 × 10 -4 and a R = 5 × 10 -5 .

## 5.4 Training Process

## Algorithm 1 Training process of InSituNet.

Input:

Training data includes parameters { Psim , Pvis , Pview } 0: N -1 and the corresponding images I 0: N -1 . Initial weights w and u of R w and D u , respectively. The feature comparator F .

Output: Optimized weights w and u

- 1: Repeat:
- 2: { Psim , Pvis , Pview } 0: b -1 , I 0: b -1 sampled from training data
- 3: ˆ I 0: b -1 ← R w ( { Psim , Pvis , Pview } 0: b -1 )
- 4: u ← Adam ( GLYPH&lt;209&gt; u L adv D ( I 0: b -1 , ˆ I 0: b -1 ; u ) , u a , D , b 1 , b 2 )
- 5: w ← Adam ( GLYPH&lt;209&gt; w L ( I 0: b -1 , ˆ I 0: b -1 ; w ) , w a , R , b 1 , b 2 )
- 6: Until exit criterion is satisfied

The process of training our regressor and discriminator is shown in Algorithm 1. Given the training data collected in situ, namely, N pairs of paramters { Psim , Pvis , Pview } 0: N -1 and the corresponding images I 0: N -1 , we first initialize the network weights w and u using the orthogonal initialization [54]. Then, the discriminator and regressor are updated alternatively by using the stochastic gradient descent until the exit criterion is satisfied. The exit criterion used in this work is the maximum number of iterations, which is set to 125,000 because the loss converged in our cases after 125,000 iterations.

In each iteration, a batch of parameters { Psim , Pvis , Pview } 0: b -1 and the corresponding images I 0: b -1 are sampled from the training data (line 2), where b is the batch size. Next, the current R w takes { Psim , Pvis , Pview } 0: b -1 as inputs and produces ˆ I 0: b -1 (line 3). According to the loss L adv D defined on I 0: b -1 and ˆ I 0: b -1 in Equation 5, the weights of the discriminator are updated (line 4). Similarly, the weights of the regressor are updated as well, according to the loss function L (defined in Equations 2, 3, and 4), which is computed using the feature comparator F and the updated discriminator D u (line 5). When updating the weights u and w , the gradients GLYPH&lt;209&gt; u and GLYPH&lt;209&gt; w of the loss functions L adv D and L are computed, respectively. With GLYPH&lt;209&gt; u and GLYPH&lt;209&gt; w , the weights u and w are updated through two Adam optimizers using the learning rates discussed in the preceding section.

## 6 PARAMETER SPACE EXPLORATION WITH INSITUNET

Fig. 7. Visual interface for parameter space exploration. (a) The three groups of parameters: simulation, visual mapping, and view parameters. (b) The predicted visualization image and the sensitivity analysis result.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-bc8c3496-93c4-4f82-8292-422c096d0a2f.jpg)

With the trained InSituNet, users can perform parameter space exploration of ensemble simulations from two perspectives. First, with InSituNet's forward propagations, users can interactively infer the visualization results for arbitrary parameters within the parameter space. Second, using InSituNet's backward propagations, users can investigate the sensitivity of different parameters and thus have better understanding on parameter selections. To support the parameter space exploration, we built an interactive visual interface as shown in Figure 7, which contains two views: Parameters View (Figure 7(a)) and Visualization View (Figure 7(b)). In the following, we explain how users can perform parameter space exploration with this visual interface.

## 6.1 Inference of Visualization Results

InSituNet is able to interactively infer the visualization results for any user-selected parameter values. As shown in Figure 7(a), the three

Table 1. Datasets and timings: k controls the size of InSituNet to cope with datasets in different complexities; diversity [63] measures how diverse the generated images are; t sim , t vis , and t tr are timings for running ensemble simulations, visualizing data in situ, and training InSituNet, respectively; t f p and t bp are timings for a forward and backward propagation of the trained InSituNet, respectively.

| Simulation    | P sim         | P sim   | P                                         | P view k   | P view k   | P view k   | Diversity   | Size (GB)   | Size (GB)   | Size (GB)   | Performance   | Performance   | Performance   | Performance   | Performance   |
|---------------|---------------|---------|-------------------------------------------|------------|------------|------------|-------------|-------------|-------------|-------------|---------------|---------------|---------------|---------------|---------------|
| Simulation    | Name          | Number  | vis                                       | Name       | Number     |            |             | Raw         | Image       | Network     | t sim (hr)    | t vis (hr)    | t tr (hr) t   | f p (s)       | t bp (s)      |
| SmallPoolFire | Ck , C        | 4,000   | pseudo-coloring with 5 color schemes      | N/A        | N/A        | 32         | 2.72        | ≈ 25.0      | 0.43        | 0.06        | 1,420.0       | 6.45          | 16.40         | 0.031         | 0.19          |
| Nyx           | OmM , OmB , h | 500     | volume rendering with a transfer function | q , f      | 100        | 48         | 1.72        | ≈ 30.0      | 3.92        | 0.12        | 537.5         | 8.47          | 18.02         | 0.033         | 0.22          |
| MPAS-Ocean    | BwsA          | 300     | isosurface visualization with 3 isovalues | q , f      | 100        | 48         | 1.75        | ≈ 300.0     | 3.46        | 0.15        | 229.5         | 10.73         | 18.13         | 0.033         | 0.23          |

groups of input parameters for the InSituNet are visualized by using different GUI widgets. For the simulation and view parameters, because their values are usually in continuous ranges, we visualize them using slider bars whose ranges are clipped to the corresponding parameters' predefined value ranges. Users are able to select arbitrary parameter values by interacting with those sliders. For the visual mapping parameters, users can switch among a set of predetermined options using the ratio buttons, for example, selecting different isovalues for isosurface visualizations, as shown in Figure 7(a).

The selected values for the three groups of parameters are fed into the trained InSituNet. Through a forward propagation of the network, which takes around 30 ms, the corresponding visualization image for the given set of parameters is generated and visualized in the Visualization View, as shown in Figure 7(b).

SmallPoolFire is a 2D combustion simulation from the OpenFOAM simulation package [66]. We used it as a test case to evaluate InSituNet by studying two parameters: a turbulence parameter Ck ∈ [ 0 0925 0 0975 . , . ] and a combustion parameter C ∈ [ 4 99 5 01 . We . , . ] sampled 4,000 parameter settings from the parameter space: 3,900 for training and 100 for testing. Images were generated for the temperature field by using pseudo-coloring with five predefined color schemes. To study how diverse the generated images are, we use the method proposed by Wang et al. [63], which measures the diversity as the reciprocal of the average structural similarity (SSIM) between every pair of images. The diversity of the images in this dataset is 2.72, which means the average SSIM is smaller than 0.4.

## 6.2 Sensitivity Analysis on Simulation Parameters

Because InSituNet is differentiable, users can perform sensitivity analysis for the simulation parameters using the network's backward propagations. Specifically, users can compute the derivative of a scalar value derived from the generated image (e.g., L 1 norm of the pixel values) with respect to a selected simulation parameter. The absolute value of the derivative can be treated as the sensitivity of the parameter, which indicates how much the generated image will change if the parameter gets changed. Note that the sensitivity analysis in this work is used to reflect the changes (with respect to the parameters) in the image space rather than the data space. Inspired by [6], our analysis includes overall sensitivity analysis and subregion sensitivity analysis.

In overall sensitivity analysis, we focus on analyzing the sensitivity of the entire image with respect to each simulation parameter across its value range. To this end, we sweep each parameter across its value range while fixing the values of other parameters. Images are then generated from the selected parameter values and aggregated into a scalar (i.e., the L 1 norm of the pixel values). The aggregated scalar values are then back propagated through the InSituNet to obtain the sensitivity of the selected parameter values. In the end, a list of sensitivity values is returned for each parameter and visualized as a line chart on top of the slider bar corresponding to the parameter (Figure 7(a1)) to indicate how sensitive the parameter is across its value range.

In subregion sensitivity analysis, we analyze the sensitivity of a selected parameter for different subregions of the generated image. This analysis is done by partitioning the visualization image into blocks and computing the sensitive of the parameter for the L 1 norm of the pixel values in each block. The computed sensitivity values are then color coded from white to red and overlaid on top of the visualization image to indicate what regions are more sensitive with respect to the selected parameter (red blocks in Figure 7(b1)).

## 7 RESULTS

We evaluated InSituNet using combustion, cosmology, and ocean simulations (Section 7.1) from four aspects: (1) providing implementation details and analyzing performance (Section 7.2); (2) evaluating the influence of different hyperparameters (Section 7.3); (3) comparing with alternative methods (Section 7.4); and (4) performing parameter space exploration and analysis with case studies (Section 7.5).

## 7.1 Ensemble Simulations

We evaluated the proposed approach using three ensemble simulations: SmallPoolFire [66], Nyx [4], and MPAS-Ocean [51]. They are summarized in Table 1 (left) and detailed below.

Nyx is a cosmological simulation developed by Lawrence Berkeley National Laboratory. Based on the scientists' suggestion, we studied three parameters: the total matter density ( OmM ∈ [ 0 12 0 155 ), the . , . ] total density of baryons ( OmB ∈ [ 0 0215 0 0235 ), and the Hubble . , . ] constant ( h ∈ [ 0 55 0 85 ). We sampled 500 parameter settings from . , . ] the parameter space: 400 for training and 100 for testing. The simulation was conducted with each parameter setting and generated a 256 × 256 × 256 volume representing the log density of the dark matters. The volume was visualized in situ by using volume rendering with a predefined transfer function of the wave colormap 1 and from 100 different viewpoints. The diversity of the generated images is 1.72.

MPAS-Ocean is a global ocean simulation developed by Los Alamos National Laboratory. Based on the domain scientists' interest, we studied the parameter that controls the bulk wind stress amplification ( BwsA ∈ [ 1 4 ). We generated 300 ensemble members with different , ] BwsA values. We used 270 of them for training and the rest for testing. The isosufaces of the temperature field (with isovalue= { 15, 20, 25 } ) were extracted and visualized from 100 different viewpoints for each ensemble member. The isosurfaces were colored based on salinity, using the colormap suggested by Samsel et al. [53]. The diversity of the generated images is 1.75.

## 7.2 Implementation and Performance

The proposed approach consists of three components: the in situ data collection, the training of InSituNet, and the visual exploration and analysis component. We discuss the implementation details and performance of the three components in the following.

The in situ visualization was implemented by using ParaView Catalyst 2 following the Cinema framework [1,2]. The simulations and in situ visualization were conducted on a supercomputer of 648 computation nodes. Each node contains an Intel Xeon E5-2680 CPU with 14 cores and 128 GB of main memory. We used 1, 28, and 128 processes, respectively, for the SmallPoolFire, Nyx, and MPAS-Ocean simulations. InSituNet was implemented in PyTorch 3 and trained with an NVIDIA DGX-1 system, which contains 8 NVIDIA V100 GPUs with NVlink. The visual interface was implemented based on a web server/client framework. The interface was implemented with D3.js on the client side, and the images were generated from a Python server (with the assist of the trained InSituNet) and sent to the client for visualization. The visual exploration and analysis were tested on a desktop with an Intel Core i7-4770 CPU and an NVIDIA 980Ti GPU.

The space and computation costs using the proposed approach for the three different datasets are listed in Table 1 (right). The size of InSituNet is less than 1% and 15% of the raw simulation data and the image data, respectively. In terms of data reduction, we also compare

1 https://sciviscolor.org

2 https://www.paraview.org/in-situ

3 https://pytorch.org

© 2019 IEEE. This is the author's version of the article that has been published in IEEE Transactions on Visualization and Computer Graphics. The final version of this record is available at: 10.1109/TVCG.2019.2934312

Fig. 8. Qualitative comparison of InSituNet trained with different loss functions. Combining L f eat and L adv R gives the results of high quality. our approach with several data compression methods and the results can be found in the supplementary material. The training of InSituNet generally takes more than 10 hours, but the time is much less than actually running the ensemble simulations with extra parameter settings. After training, a forward or backward propagation of InSituNet takes less than one second on a single NVIDIA 980Ti GPU.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-6a38f782-c1d1-485d-88cc-8c7715db41b8.jpg)

## 7.3 Model Evaluation for Different Hyperparameters

We evaluated InSituNet trained with different hyperparameters (i.e., loss functions, network architectures, and numbers of training samples) qualitatively and quantitatively using the data that were excluded from the training to study two questions: (1) Is InSituNet able to generate images that are close to the ground truth images? (2) How do the choices of hyperparameters influence the training results?

For quantitative evaluations, we used four metrics that focus on different aspects to compare the predicted images with the ground truth images, including peak signal-to-noise ratio (PSNR), SSIM, earth mover's distance (EMD) between color histograms [6], and Fr´ echet inception distance (FID) [28].

PSNR measures the pixel-level difference between two images using the aggregated mean squared error between image pixels. A higher PSNR indicates that the compared images are more similar pixel wise.

SSIM compares two images based on the regional aggregated statistical information (e.g., mean and standard deviation of small patches) between them. A higher SSIM means the compared images are more similar from a structural point of view.

EMD is used in [6] to quantify the distance between the color histograms of two images. A lower EMD means the compared images are more similar according to their color distributions.

FID approximates the distance between two distributions of images, which is widely used in recent image synthesis works [13, 69] as a complementary to other metrics. A lower FID suggests the two image collections are more similar statistically.

## 7.3.1 Loss Functions

We evaluated InSituNet trained with different loss functions including the mean squared error L mse , the feature reconstruction loss L f eat , the adversarial loss L adv R , and the combination of L f eat and L adv R .

Figure 8 compares the images generated by InSituNet trained with different loss functions with the ground truth (The enlarged figure can be found in the supplementary material). We can see that using L mse often generates over-smoothed images lacking high-frequency features, whereas using L f eat can mitigate the problem to some extent. Using L adv R can generate images that are as sharp as the ground truth, but the features are often not introduced in the desired positions. By combining L f eat and L adv R , we are able to generate images with sharp features, and those images are also similar to the ground truth.

Table 2. Quantitative evaluation of InSituNet trained with different loss functions. The model trained with the combination of L f eat and L adv R generates images with the best EMD and FID and only a slightly lower PSNR and SSIM compared with the model trained with L mse or L f eat .

|               |               | L mse                | L f eat                     | L adv R                     | L f eat + 10 - 2 L adv R   |
|---------------|---------------|----------------------|-----------------------------|-----------------------------|----------------------------|
| SmallPoolFire | PSNR SSIM EMD | 25.090 0.9333 0.0051 | 24.937 0.9390 0.0064 15.881 | 20.184 0.8163 0.0056 12.859 | 24.288 0.9006 0.0037       |
|               | FID           | 21.063               |                             |                             | 9.4747                     |
| Nyx           | PSNR          | 31.893               | 29.055                      | 24.592                      | 29.366                     |
|               | SSIM          | 0.8684               | 0.8698                      | 0.7081                      | 0.8336                     |
|               | EMD           | 0.0037               | 0.0083                      | 0.0064                      | 0.0022                     |
|               | FID           | 60.825               | 54.670                      | 24.036                      | 6.2694                     |
|               | PSNR          | 26.944               | 26.267                      |                             | 24.791                     |
| MPAS-Ocean    | SSIM          | 0.8908               | 0.8885                      | 17.099 0.7055               | 0.8655                     |
| MPAS-Ocean    | EMD           | 0.0025               | 0.0044                      | 0.0036                      | 0.0017                     |
| MPAS-Ocean    | FID           | 115.74               | 120.37                      | 28.927                      | 21.395                     |

Table 2 reports the quantitative results from using different loss functions. We found that using L mse gives the best PSNR, because the network using L mse is trained to minimize the mean squared error (i.e., maximize the PSNR). Using L f eat gives the best SSIM in some cases, because it focuses more on the structure of the images. However, using L mse or L f eat often results in poor performance regarding EMD and FID. When training InSituNet with L adv R , the FID value can be improved, but the values of PSNR and SSIM drop a lot. By combining L f eat and L adv R , both EMD and FID improved a lot, though the PSNR and SSIM got slightly worse than using L mse or L f eat .

Fig. 9. Images generated by InSituNet trained with L f eat that uses different layers after the first pooling layer of VGG-19: (a) relu2 1 and (b) relu3 1 . Checkerboard artifacts are introduced.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-b473dfab-1360-4245-9124-4686b156de2c.jpg)

For L f eat , using which layer of the pretrained VGG-19 to extract features from images can affect the image synthesis results. Through empirical studies, we found that using any layers after the first pooling layer of VGG-19 will introduce undesired checkerboard artifacts, because of the 'inhomogeneous gradient update' of the pooling layer [3], as shown in Figure 9. Hence, we use the last layer right before the first pooling layer, which is the layer relu1 2 .

Table 3. Evaluating the weight l of L adv R : l = 0 01 . provides the results that balance the PSNR, SSIM, EMD, and FID.

|      |   l = 0 . 005 |   l = 0 . 01 |   l = 0 . 02 |   l = 0 . 04 |
|------|---------------|--------------|--------------|--------------|
| PSNR |       30.043  |      29.366  |      29.04   |      27.232  |
| SSIM |        0.8619 |       0.8336 |       0.8253 |       0.768  |
| EMD  |        0.0041 |       0.0022 |       0.0023 |       0.0025 |
| FID  |       21.267  |       6.2694 |       6.6819 |       9.8992 |

We also evaluated the influence of the weight l for L adv R when combining it with L f eat (defined in Equation 2), and the results are shown in Table 3. We found that increasing l over 0.01 cannot improve the accuracy of the generated images any further. In addition, a small l (i.e., 0.005) will hurt the image accuracy in terms of EMD and FID, although the value of PSNR and SSIM can be improved slightly. We thereby set l to 0.01 to balance its effects on the four metrics.

## 7.3.2 Network Architectures

We evaluated InSituNet with different network architectures in terms of the accuracy of predicted images, the network size, and the training time. As mentioned in Section 5.1, the architecture of our network is controlled by a constant k , which controls the number of convolutional kernels in the intermediate layers. In this experiment, we evaluated four k values: 16, 32, 48, and 64.

Figure 10 shows the PSNR and EMD of images generated by InSituNet with the four k values. We can see that InSituNet with larger k values can generate more (or at least equally) accurate images, because

Fig. 10. Quantitative evaluation of different network architectures controlled by k with (a) PSNR and (b) EMD. Size and training time of different network architectures con-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-f1ba93ed-a825-4954-abe7-451d89dd3097.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-1732effc-0b03-486e-8fe1-2a4cd9969053.jpg)

Table 4. trolled by k for the Nyx dataset.

|                    |   k = 16 |   k = 32 |   k = 48 |   k = 64 |
|--------------------|----------|----------|----------|----------|
| Network Size (MB)  |    26.4  |    67.4  |   125.2  |   199.6  |
| Training Time (hr) |    13.73 |    16.42 |    18.02 |    20.17 |

a larger k gives more expressive power to the neural network. On the other hand, training InSituNet with a larger k also costs more time, and more storage will be needed to store the networks, as shown in Table 4 using the Nyx dataset as an example. Hence, to balance the accuracy of the generated images and the cost from both computation and storage, we set k to 32, 48, and 48 for the SmallPoolFire, Nyx, and MPAS-Ocean, respectively.

## 7.3.3 Number of Ensemble Runs used for Training

Table 5. Evaluation of the number of ensemble runs used for training.

| Simulation    | # Ensemble Runs    | PSNR                        | SSIM                        | EMD                  | FID                  |
|---------------|--------------------|-----------------------------|-----------------------------|----------------------|----------------------|
| SmallPoolFire | 900 1900 2900 3900 | 21.842 23.192 23.932 24.288 | 0.8714 0.9016 0.9018 0.9006 | 0.0040 0.0036 0.0037 | 14.398 11.732 9.5813 |
| Nyx           | 100 200            | 28.108 29.404               | 0.7951 0.8319               | 0.0022               | 9.8818 6.5481        |
|               |                    |                             |                             |                      | 6.4239               |
|               |                    |                             |                             | 0.0037               | 9.4747               |
|               |                    |                             |                             | 0.0025               |                      |
|               | 300                | 29.398                      | 0.8326                      | 0.0023               |                      |
|               | 400                | 29.366                      | 0.8336                      | 0.0022               | 6.2694               |
| MPAS-Ocean    | 70                 | 24.347                      | 0.8554                      | 0.0016               | 37.229               |
| MPAS-Ocean    | 140                | 24.593                      | 0.8607                      | 0.0017               | 28.380               |
| MPAS-Ocean    | 210                | 24.732                      | 0.8643                      | 0.0017               | 22.794               |
| MPAS-Ocean    | 270                | 24.791                      | 0.8655                      | 0.0017               | 21.395               |

We compared InSituNet trained using different numbers of ensemble runs (Table 5) to study how many ensemble runs will be needed to train a good model for the three simulations. We found that this number is different in different simulations, depending on the complexity of the mapping between simulation parameters and visualization results. Experiment results show that the accuracy of generated images becomes stable when the number of ensemble runs is greater than 2,900, 200, and 210 for the SmallPoolFire, Nyx, and MPAS-Ocean simulation, respectively. As a result, we used 3,900, 400, and 270 runs from the three simulations to train InSituNet for the rest of the study.

## 7.4 Comparison with Alternative Methods

We compared our method with two alternative methods including interpolating images from the training data that close to the target image and the GAN-based volume rendering method (GAN-VR) [6] using the Nyx dataset. For the interpolation method, we sample g images from the training data whose parameter settings are the top g closest to the parameter setting of the test image and interpolate the sampled images using inverse distance weighting interpolation [56]. We evaluated g from 1 to 5 and present the result of g = 3 in this section because it balances the four metrics (More results are in the supplementary material). For GAN-VR, we incorporated the simulation parameters into both the opacity GAN and the opacity-to-color translation GAN and removed the transfer function related parameters because we used a fixed transfer function for this dataset. For InSituNet, we selected a network architecture whose size is not greater than the size of GAN-VR network, for a fair comparison.

Figure 11 compares the ground truth images with the images generated by using interpolation, GAN-VR, and InSituNet. With the new network architecture (e.g., the projection-based condition incorporation method in Section 5.1), loss functions (e.g., the feature reconstruction

Fig. 11. Comparison of the images generated using interpolation, GANVR, and InSituNet with the ground truth images.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-fe416a96-87f5-4872-9757-1b4dca3c9b7c.jpg)

Table 6. Quantitative comparison of images generated with interpolation, GAN-VR, and InSituNet.

|               | Network Size   |   PSNR |   SSIM |    EMD |     FID |
|---------------|----------------|--------|--------|--------|---------|
| Interpolation | N/A            | 23.932 | 0.6985 | 0.007  | 58.571  |
| GAN-VR        | 80.5 MB        | 20.67  | 0.6274 | 0.0056 | 38.355  |
| InSituNet     | 67.5 MB        | 28.471 | 0.8034 | 0.0023 |  9.0152 |

loss in Section 5.2), and training strategies (e.g., the spectral normalization in Section 5.3), InSituNet can generate results that better preserve features compared with the other two methods. The quantitative comparisons between the three methods are shown in Table 6. InSituNet outperforms the other two methods in all four metrics.

## 7.5 Parameter Space Exploration

This section demonstrates the effectiveness of our deep image synthesis driven parameter space exploration through case studies on the Nyx and MPAS-Ocean simulations.

## 7.5.1 Case Study with the Nyx Simulation

Fig. 12. Parameter space exploration with the Nyx simulation. For the selected parameter values, the sensitivity of different parameters is estimated and visualized as line charts on the left, whereas the predicted image is visualized on the right.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-12d7c645-50e1-489b-bc1c-2965516d97e0.jpg)

Our first case study is focused on investigating the influence of different simulation parameters (i.e., OmM OmB , , and h ) on the Nyx simulation. The explorations of different visualization settings can be found in our associated video.

Figure 12 shows a selected parameter setting with the predicted visualization image. To understand the influence of each parameter, we computed the sensitivity of the three parameters with respect to the L 1 norm of the predicted image, shown as the three line charts in Figure 12. From the scale of the three charts (i.e., the values along the vertical axes), we see that parameter h is more sensitive to parameter OmM and parameter OmM is more sensitive to parameter OmB .

Focusing on the most sensitive parameter, namely, parameter h , we explored how it affects the visual appearance of the predicted images. Figure 13 shows five images predicted by using five different h values, while parameter OmM and OmB are fixed at the values shown on the two corresponding slider bars in Figure 12. We first evaluate the accuracy of the sensitivity curve (blue curve in Figure 13) computed by backpropagation with the central difference method. To this end, we first regularly sample the simulation parameters along the curve (128 samples are drawn) and then generate the visualization images

© 2019 IEEE. This is the author's version of the article that has been published in IEEE Transactions on Visualization and Computer Graphics. The final version of this record is available at: 10.1109/TVCG.2019.2934312

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-4dbbc2a8-c11c-4ddf-9d6f-aab2d6921880.jpg)

Figure 15 shows the predicted images when using different BwsA values. All images are generated from the temperature field (isovalue=15) of MPAS-Ocean and from the same viewpoint. The first row of images shows the result of forward inference with different BwsA values, whereas the second row of images overlays the subregion sensitivity maps onto the corresponding images of the first row. The labeled regions (i.e., Figure 15a, b) change the most when adjusting the value of BwsA , and the subregion sensitivity maps on the second row echo these bigger changes, as indicated by the darker red color.

Fig. 13. Comparison of the visual appearance of the predicted images using different h values to see the effect of this simulation parameter. with respect to the sampled simulation parameters. The L1 norm of the generated images is then computed and used to compute the sensitive curve using the central difference method. The result is shown as the orange curve in Figure 13. We can see that the sensitivity curves generated with the two methods are similar. From the guidance provided by the line chart in Figure 13, we see that parameter h is more sensitive in the first half of its range (i.e., the left side of the dashed line). The three images generated using h values from this range demonstrate a bigger variance compared with the two images shown on the right (which are images generated by using h values from the second half of its range).

7.5.2 Case Study with the MPAS-Ocean Simulation

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-bb56dfa8-bd1c-40c9-9f5a-df953c4cb739.jpg)

Fig. 14. Predicted images of the MPAS-Ocean dataset for different isosurfaces and viewpoints, which reasonably reflect the change of view projections and shading effects.

Fig. 15. Forward prediction (top row) and backward subregion sensitivity analysis (bottom row) for different BwsA . Regions that influenced by BwsA (i.e., regions a and b) are highlighted by the sensitivity map.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00407v3-figure-75865421-5c81-49bc-90de-6ae33ec232ce.jpg)

Our next case study explores different parameter settings for the MPAS-Ocean simulation and demonstrates the subregion sensitivity analysis for the simulation parameter BwsA , which characterizes the bulk wind stress. Note that here we focus only on exploration and analysis of new parameter settings, the comparison between the predicted and ground truth images is discussed in the previous sections.

Figure 14 shows isosurface visualizations of the temperature field with three different isovalues from six different viewpoints. The value of parameter BwsA is fixed at 1 in this study. The images reasonably reflect the change of view projections and shading effects. More exploration and visualization results can be found in our associated video.

## 8 LIMITATIONS, DISCUSSION, AND FUTURE WORK

This section discusses several directions that we would like to explore in the future: (1) improving the flexibility in exploring arbitrary visual mapping parameters; (2) increasing the accuracy of predicted images; and (3) increasing the resolution of predicted images.

One limitation of our approach is that we restricted the users' ability in exploring arbitrary visual mapping parameters, for example, exhausting all possible transfer functions for volume rendering. Instead, we allow users to switch only among several predefined visual mappings, for example, the three isovalues when exploring the MPAS-Ocean data. Theoretically, training a deep learning model to predict visualization images for arbitrary simulation and visualization parameters is feasible. However, it will require a large number of training images to cover the joint space of all possible simulation and visualization parameters. For example, in order to train a model that can predict volume rendering results of a single volume data for arbitrary transfer functions, 200,000 training images are required, as shown in [6]. Consequently, the size of the training data may even exceed the size of the raw simulation data, which offsets the benefit of in situ visualization. Considering this issue, we would like to explore deep learning techniques that do not require a large number of training samples, such as one- or zero-shot learning, to improve the flexibility of exploration.

Similar to most other machine learning techniques, generating prediction results that are exactly the same as the ground truth is extraordinary difficult. By taking advantage of recent advances in deep learning for image synthesis, the proposed approach has already outperformed other image synthesis based visualization techniques in terms of the fidelity and accuracy of the generated images (see the comparison in Section 7). However, we believe further improvement is still possible, and we would like explore other network architectures and/or other loss functions to improve our deep image synthesis model.

Our network architecture limits the resolution of output images to 256 × 256, which might not be sufficient for some high-resolution simulation data. We believe that our network architecture has the potential to generate images with higher resolutions by adding more residual blocks, and we will investigate this approach in the future.

## 9 CONCLUSION

In this work, we propose InSituNet, a deep learning based image synthesis model supporting the parameter space exploration of large-scale ensemble simulations visualized in situ. The model is trained to learn the mapping from ensemble simulation parameters to visualizations of the corresponding simulation outputs, conditioned on different visualization settings (i.e., visual mapping and view parameters). With a trained InSituNet, users can generate visualizations of simulation outputs with different simulation parameters without actually running the expensive simulation, as well as synthesize new visualizations with different visualization settings that are not used during the runs. Additionally, an interactive visual interface is developed to explore the space of different parameters and investigate their sensitivity using the trained InSituNet. Through both quantitative and qualitative evaluations, we validated the effectiveness of InSituNet in analyzing ensemble simulations that model different physical phenomena.

## ACKNOWLEDGMENTS

This work was supported in part by US Department of Energy Los Alamos National Laboratory contract 47145 and UT-Battelle LLC contract 4000159447 program manager Laura Biven.

## REFERENCES

- [1] J. Ahrens, S. Jourdain, P. O'Leary, J. Patchett, D. H. Rogers, P. Fasel, A. Bauer, M. Petersen, F. Samsel, and B. Boeckel. In situ MPAS-Ocean image-based visualization. In Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, Visualization &amp; Data Analytics Showcase , 2014.
- [2] J. Ahrens, S. Jourdain, P. O'Leary, J. Patchett, D. H. Rogers, and M. Petersen. An image-based approach to extreme scale in situ visualization and analysis. In Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis , pp. 424-434, 2014.
- [3] A. Aitken, C. Ledig, L. Theis, J. Caballero, Z. Wang, and W. Shi. Checkerboard artifact free sub-pixel convolution: A note on sub-pixel convolution, resize convolution and convolution resize. arXiv preprint arXiv:1707.02937 , 2017.
- [4] A. S. Almgren, J. B. Bell, M. J. Lijewski, Z. Luki´ c, and E. V. Andel. Nyx: A massively parallel AMR code for computational cosmology. The Astrophysical Journal , 765(1):39, 2013.
- [5] A. C. Bauer, H. Abbasi, J. Ahrens, H. Childs, B. Geveci, S. Klasky, K. Moreland, P. O'Leary, V. Vishwanath, B. Whitlock, and E. W. Bethel. In situ methods, infrastructures, and applications on high performance computing platforms. Computer Graphics Forum , 35(3):577-597, 2016.
- [6] M. Berger, J. Li, and J. A. Levine. A generative model for volume rendering. IEEE Transactions on Visualization and Computer Graphics , 25(4):1636-1650, 2019.
- [7] S. Bergner, M. Sedlmair, T. M¨ oller, S. N. Abdolyousefi, and A. Saad. ParaGlide: Interactive parameter space partitioning for computer simulations. IEEE Transactions on Visualization and Computer Graphics , 19(9):1499-1512, 2013.
- [8] T. Biedert and C. Garth. Contour tree depth images for large data visualization. In Proceedings of Eurographics Symposium on Parallel Graphics and Visualization , pp. 77-86, 2015.
- [9] A. Biswas, G. Lin, X. Liu, and H.-W. Shen. Visualization of time-varying weather ensembles across multiple resolutions. IEEE Transactions on Visualization and Computer Graphics , 23(1):841-850, 2017.
- [10] A. Bock, A. Pembroke, M. L. Mays, L. Rastaetter, T. Ropinski, and A. Ynnerman. Visual verification of space weather ensemble simulations. In Proceedings of 2015 IEEE Scientific Visualization Conference , pp. 17-24, 2015.
- [11] M. Booshehrian, T. M¨ oller, R. M. Peterman, and T. Munzner. Vismon: Facilitating analysis of trade-offs, uncertainty, and sensitivity in fisheries management decision making. In Computer Graphics Forum , vol. 31, pp. 1235-1244. Wiley Online Library, 2012.
- [12] P.-T. Bremer, G. H. Weber, J. Tierny, V. Pascucci, M. S. Day, and J. B. Bell. Interactive exploration and analysis of large-scale simulations using topology-based data segmentation. IEEE Transactions on Visualization and Computer Graphics , 17(9):1307-1324, 2011.
- [13] A. Brock, J. Donahue, and K. Simonyan. Large scale GAN training for high fidelity natural image synthesis. In Proceedings of International Conference on Learning Representations , 2019.
- [14] S. Bruckner and T. M¨ oller. Result-driven exploration of simulation parameter spaces for visual effects design. IEEE Transactions on Visualization and Computer Graphics , 16(6):1468-1476, 2010.
- [15] H. Chen, S. Zhang, W. Chen, H. Mei, J. Zhang, A. Mercer, R. Liang, and H. Qu. Uncertainty-aware multidimensional ensemble data visualization and exploration. IEEE Transactions on Visualization and Computer Graphics , 21(9):1072-1086, 2015.
- [16] D. Coffey, C.-L. Lin, A. G. Erdman, and D. F. Keefe. Design by dragging: An interface for creative forward and inverse design with simulation ensembles. IEEE Transactions on Visualization and Computer Graphics , 19(12):2783-2791, 2013.
- [17] S. Di and F. Cappello. Fast error-bounded lossy HPC data compression with SZ. In Proceedings of International Parallel and Distributed Processing Symposium , pp. 730-739, 2016.
- [18] C. Dong, C. C. Loy, K. He, and X. Tang. Image super-resolution using deep convolutional networks. IEEE Transactions on Pattern Analysis and Machine Intelligence , 38(2):295-307, 2016.
- [19] A. Dosovitskiy, J. T. Springenberg, and T. Brox. Learning to generate chairs with convolutional neural networks. In Proceedings of 2015 IEEE Conference on Computer Vision and Pattern Recognition , pp. 1538-1546, 2015.
- [20] S. Dutta, C.-M. Chen, G. Heinlein, H.-W. Shen, and J.-P. Chen. In
- situ distribution guided analysis and visualization of transonic jet engine simulations. IEEE Transactions on Visualization and Computer Graphics , 23(1):811-820, 2017.
- [21] O. Fernandes, S. Frey, F. Sadlo, and T. Ertl. Space-time volumetric depth images for in-situ visualization. In Proceedings of 2014 IEEE Symposium on Large Data Analysis and Visualization , pp. 59-65, 2014.
- [22] S. Frey, F. Sadlo, and T. Ertl. Explorable volumetric depth images from raycasting. In Proceedings of 2013 XXVI Conference on Graphics, Patterns and Images , pp. 123-130, 2013.
- [23] L. A. Gatys, A. S. Ecker, and M. Bethge. A neural algorithm of artistic style. arXiv preprint arXiv:1508.06576 , 2015.
- [24] L. A. Gatys, A. S. Ecker, and M. Bethge. Texture synthesis using convolutional neural networks. In Proceedings of Advances in Neural Information Processing Systems , pp. 262-270, 2015.
- [25] I. J. Goodfellow, J. Pouget-Abadie, M. Mirza, B. Xu, D. Warde-Farley, S. Ozair, A. Courville, and Y. Bengio. Generative adversarial nets. In Proceedings of Advances in Neural Information Processing Systems , pp. 2672-2680, 2014.
- [26] J. Han, J. Tao, and C. Wang. FlowNet: A deep learning framework for clustering and selection of streamlines and stream surfaces. IEEE Transactions on Visualization and Computer Graphics , 2018, Early Access.
- [27] K. He, X. Zhang, S. Ren, and J. Sun. Deep residual learning for image recognition. In Proceedings of 2016 IEEE Conference on Computer Vision and Pattern Recognition , pp. 770-778, 2016.
- [28] M. Heusel, H. Ramsauer, T. Unterthiner, B. Nessler, and S. Hochreiter. GANs trained by a two time-scale update rule converge to a local nash equilibrium. In Proceedings of Advances in Neural Information Processing Systems , pp. 6626-6637, 2017.
- [29] S. Hochreiter and J. Schmidhuber. Long short-term memory. Neural Computation , 9(8):1735-1780, 1997.
- [30] F. Hong, J. Zhang, and X. Yuan. Access pattern learning with long shortterm memory for parallel particle tracing. In Proceedings of 2018 IEEE Pacific Visualization Symposium , pp. 76-85, 2018.
- [31] X. Hou, L. Shen, K. Sun, and G. Qiu. Deep feature consistent variational autoencoder. In Proceedings of 2017 IEEE Winter Conference on Applications of Computer Vision , pp. 1133-1141, 2017.
- [32] J. Johnson, A. Alahi, and F.-F. Li. Perceptual losses for real-time style transfer and super-resolution. In Proceedings of European Conference on Computer Vision , pp. 694-711, 2016.
- [33] D. P. Kingma and J. Ba. Adam: A method for stochastic optimization. In Proceedings of International Conference on Learning Representations , 2015.
- [34] K. Kurach, M. Lucic, X. Zhai, M. Michalski, and S. Gelly. The GAN landscape: Losses, architectures, regularization, and normalization. arXiv preprint arXiv:1807.04720 , 2018.
- [35] S. Lakshminarasimhan, N. Shah, S. Ethier, S.-H. Ku, C. S. Chang, S. Klasky, R. Latham, R. Ross, and N. F. Samatova. ISABELA for effective in situ compression of scientific data. Concurrency and Computation: Practice and Experience , 25(4):524-540, 2013.
- [36] C. Ledig, L. Theis, F. Husz´ ar, J. Caballero, A. Cunningham, A. Acosta, A. Aitken, A. Tejani, J. Totz, Z. Wang, and W. Shi. Photo-realistic single image super-resolution using a generative adversarial network. In Proceedings of 2017 IEEE Conference on Computer Vision and Pattern Recognition , pp. 105-114, 2017.
- [37] M. Luˇ ci´, K. Kurach, M. Michalski, S. Gelly, and O. Bousquet. Are GANs c created equal? A large-scale study. In Proceedings of Advances in Neural Information Processing Systems , pp. 700-709, 2018.
- [38] K.-L. Ma. In situ visualization at extreme scale: Challenges and opportunities. IEEE Computer Graphics and Applications , 29(6):14-19, 2009.
- [39] K. Matkovi´ c, D. Graˇ canin, B. Klarin, and H. Hauser. Interactive visual analysis of complex scientific data as families of data surfaces. IEEE Transactions on Visualization and Computer Graphics , 15(6):1351-1358, 2009.
- [40] M. Mirza and S. Osindero. Conditional generative adversarial nets. arXiv preprint arXiv:1411.1784 , 2014.
- [41] T. Miyato, T. Kataoka, M. Koyama, and Y. Yoshida. Spectral normalization for generative adversarial networks. In Proceedings of International Conference on Learning Representations , 2018.
- [42] T. Miyato and M. Koyama. cGANs with projection discriminator. In Proceedings of International Conference on Learning Representations , 2018.
- [43] V. Nair and G. E. Hinton. Rectified linear units improve restricted boltz-

© 2019 IEEE. This is the author's version of the article that has been published in IEEE Transactions on Visualization and Computer Graphics. The final version of this record is available at: 10.1109/TVCG.2019.2934312

mann machines. In Proceedings of International Conference on Machine Learning , pp. 807-814, 2010.

Visualization and Computer Graphics , 24(6):1905-1917, 2018.

- [44] H. Obermaier, K. Bensema, and K. I. Joy. Visual trends analysis in timevarying ensembles. IEEE Transactions on Visualization and Computer Graphics , 22(10):2331-2342, 2016.
- [45] D. Orban, D. F. Keefe, A. Biswas, J. Ahrens, and D. Rogers. Drag and track: A direct manipulation interface for contextualizing data instances within a continuous parameter space. IEEE Transactions on Visualization and Computer Graphics , 25(1):256-266, 2019.
- [46] D. Pathak, P. Kr¨ ahenb¨hl, J. Donahue, T. Darrell, and A. A. Efros. Context u encoders: Feature learning by inpainting. In Proceedings of 2016 IEEE Conference on Computer Vision and Pattern Recognition , pp. 2536-2544, 2016.
- [47] H. Piringer, W. Berger, and J. Krasser. Hypermoval: Interactive visual validation of regression models for real-time simulation. In Computer Graphics Forum , vol. 29, pp. 983-992. Wiley Online Library, 2010.
- [48] J. Poco, A. Dasgupta, Y. Wei, W. Hargrove, C. R. Schwalm, D. N. Huntzinger, R. Cook, E. Bertini, and C. T. Silva. Visual reconciliation of alternative similarity spaces in climate modeling. IEEE Transactions on Visualization and Computer Graphics , 20(12):1923-1932, 2014.
- [49] K. Potter, A. Wilson, P.-T. Bremer, D. Williams, C. Doutriaux, V. Pascucci, and C. R. Johnson. Ensemble-vis: A framework for the statistical visualization of ensemble data. In 2009 IEEE International Conference on Data Mining Workshops , pp. 233-240. IEEE, 2009.
- [50] S. Reed, Z. Akata, X. Yan, L. Logeswaran, B. Schiele, and H. Lee. Generative adversarial text to image synthesis. In Proceedings of International Conference on Machine Learning , pp. 1060-1069, 2016.
- [51] T. Ringler, M. Petersen, R. L. Higdon, D. Jacobsen, P. W. Jones, and M. Maltrud. A multi-resolution approach to global ocean modeling. Ocean Modelling , 69:211-232, 2013.
- [52] D. E. Rumelhart, G. E. Hinton, and R. J. Williams. Learning internal representations by error propagation. In D. E. Rumelhart, J. L. McClelland, and CORPORATE PDP Research Group, eds., Parallel Distributed Processing: Explorations in the Microstructure of Cognition, Vol. 1 , pp. 318-362. MIT Press, Cambridge, MA, USA, 1986.
- [53] F. Samsel, M. Petersen, G. Abram, T. L. Turton, D. Rogers, and J. Ahrens. Visualization of ocean currents and eddies in a high-resolution global ocean-climate model. In Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, Visualization &amp; Data Analytics Showcase , 2015.
- [54] A. M. Saxe, J. L. McClelland, and S. Ganguli. Exact solutions to the nonlinear dynamics of learning in deep linear neural networks. arXiv preprint arXiv:1312.6120 , 2013.
- [55] M. Sedlmair, C. Heinzl, S. Bruckner, H. Piringer, and T. M¨ller. Visual o parameter space analysis: A conceptual framework. IEEE Transactions on Visualization and Computer Graphics , 20(12):2161-2170, 2014.
- [56] D. Shepard. A two-dimensional interpolation function for irregularlyspaced data. In Proceedings of the 1968 23rd ACM National Conference , pp. 517-524, 1968.
- [57] K. Simonyan and A. Zisserman. Very deep convolutional networks for large-scale image recognition. In Proceedings of International Conference on Learning Representations , 2015.
- [58] R. Splechtna, K. Matkovi´ c, D. Graˇ canin, M. Jelovi´, and H. Hauser. Interc active visual steering of hierarchical simulation ensembles. In Proceedings of 2015 IEEE Conference on Visual Analytics Science and Technology , pp. 89-96, 2015.
- [59] A. Tikhonova, C. D. Correa, and K.-L. Ma. Explorable images for visualizing volume data. In Proceedings of 2010 IEEE Pacific Visualization Symposium , pp. 177-184, 2010.
- [60] A. Tikhonova, C. D. Correa, and K.-L. Ma. An exploratory technique for coherent visualization of time-varying volume data. Computer Graphics Forum , 29(3):783-792, 2010.
- [61] A. Tikhonova, C. D. Correa, and K.-L. Ma. Visualization by proxy: A novel framework for deferred interaction with volume data. IEEE Transactions on Visualization and Computer Graphics , 16(6):1551-1559, 2010.
- [62] T. Torsney-Weir, A. Saad, T. M¨ oller, H.-C. Hege, B. Weber, J.-M. Verbavatz, and S. Bergner. Tuner: Principled parameter finding for image segmentation algorithms using visual response surface exploration. IEEE Transactions on Visualization and Computer Graphics , 17(12):1892-1901, 2011.
- [63] J. Wang, L. Gou, H. Yang, and H.-W. Shen. Ganviz: A visual analytics approach to understand the adversarial game. IEEE Transactions on
- [64] J. Wang, S. Hazarika, C. Li, and H.-W. Shen. Visualization and visual analysis of ensemble data: A survey. IEEE Transactions on Visualization and Computer Graphics , 2018, Early Access.
- [65] J. Wang, X. Liu, H.-W. Shen, and G. Lin. Multi-resolution climate ensemble parameter analysis with nested parallel coordinates plots. IEEE Transactions on Visualization and Computer Graphics , 23(1):81-90, 2017.
- [66] H. G. Weller, G. Tabor, H. Jasak, and C. Fureby. A tensorial approach to computational continuum mechanics using object-oriented techniques. Computers in Physics , 12(6):620-631, 1998.
- [67] C. Xie, W. Xu, and K. Mueller. A visual analytics framework for the detection of anomalous call stack trees in high performance computing applications. IEEE Transactions on Visualization and Computer Graphics , 25(1):215-224, 2019.
- [68] J. Xie, L. Xu, and E. Chen. Image denoising and inpainting with deep neural networks. In Proceedings of Advances in Neural Information Processing Systems , pp. 341-349, 2012.
- [69] H. Zhang, I. J. Goodfellow, D. N. Metaxas, and A. Odena. Self-attention generative adversarial networks. arXiv preprint arXiv:1805.08318 , 2018.
- [70] K. Zhang, W. Zuo, Y. Chen, D. Meng, and L. Zhang. Beyond a Gaussian denoiser: Residual learning of deep CNN for image denoising. IEEE Transactions on Image Processing , 26(7):3142-3155, 2017.
- [71] Y. Zhou, Z. Zhu, X. Bai, D. Lischinski, D. Cohen-Or, and H. Huang. Nonstationary texture synthesis by adversarial expansion. ACM Transactions on Graphics , 37(4):49:1-49:13, 2018.
- [72] J.-Y. Zhu, T. Park, P. Isola, and A. A. Efros. Unpaired image-to-image translation using cycle-consistent adversarial networks. In Proceedings of 2017 International Conference on Computer Vision , pp. 2242-2251, 2017.
Wenbin He Junpeng Wang Hanqi Guo Ko-Chih Wang Han-Wei Shen Mukund Raj Youssef S. G. Nashed Tom Peterka
InSituNet: Deep Image Synthesis for Parameter Space Exploration of Ensemble Simulations
We propose InSituNet, a deep learning based surrogate model to support
parameter space exploration for ensemble simulations that are visualized in
situ. In situ visualization, generating visualizations at simulation time, is
becoming prevalent in handling large-scale simulations because of the I/O and
storage constraints. However, in situ visualization approaches limit the
flexibility of post-hoc exploration because the raw simulation data are no
longer available. Although multiple image-based approaches have been proposed
to mitigate this limitation, those approaches lack the ability to explore the
simulation parameters. Our approach allows flexible exploration of parameter
space for large-scale ensemble simulations by taking advantage of the recent
advances in deep learning. Specifically, we design InSituNet as a convolutional
regression model to learn the mapping from the simulation and visualization
parameters to the visualization results. With the trained model, users can
generate new images for different simulation parameters under various
visualization settings, which enables in-depth analysis of the underlying
ensemble simulations. We demonstrate the effectiveness of InSituNet in
combustion, cosmology, and ocean simulations through quantitative and
qualitative evaluations.