Keywords

1 Introduction

Laser technology is very important in generating plasma waves. Plasma generated by passing a laser through a gas or a screen is utilized in various applications such as energy, new materials, semiconductor device manufacturing, and environmental fields required in the 21st century. The generation of high-efficiency plasmas is closely related to the development of high-intensity lasers. Therefore, in the study of high intensity plasma generation, visualization of intensity data is a very meaningful task [1,2,3,4].

In general, volume visualization is performed using two approaches: direct volume rendering and surface rendering. Direct volume rendering algorithms render every voxel in the volume raster directly, without conversion to geometric primitives or first converting to a different domain and then rendering from that domain. On the other hand, surface rendering algorithms first fit geometric primitives to values in the data, and then render these primitives. Therefore, we have visualized the intensity data through iso-surface methods, one of the surface rendering methods that is easy to analyze the geometric characteristics (Fig. 1).

Fig. 1.
figure 1

Direct volume rendering (left) vs surface rendering (right)

2 Problem Description

The iso-surface is a method of generating surface by connecting values having certain value, called iso-value, in volume data. The iso-value may be chosen in a purely interactive manner without any guidance. And the selected value may be inappropriate. The previous way to assist in determining iso-value is to use the histogram and previews. As a first level of support, the histogram of the dataset is displayed and labeled to enable the user to select a certain value. Previews, as a further supporting method, may be provided by rendering a down-sampled images and enable the user to select from them (Table 1).

Table 1. Data properties

The laser intensity data is characterized by scattered volumetric data. Figure 2 shows that the distribution of laser data is too broad (min: 9.0e−15, max: 2.5e + 03) to be represented through the histogram. In addition, there is a problem that big data requires a lot of time to process the preview each time. Therefore, programs such as MatLab take a built-in algorithm to compute and provide the appropriate iso-value. However, this is an advantageous method if the volume data has a clear surface, and is not a suitable method for scattered data with unclear surfaces such as laser intensity data. In addition, significant iso-values for users studying through visualization of laser intensity data may be different each time. Therefore, we propose an improved method to propose a moderate iso-value to the user [5,6,7].

Fig. 2.
figure 2

Distribution of laser intensity values

3 Camera Metering System

There is a metering system in the camera that we can see easily around us. The metering system measures and calculates the amount of light that can be confirmed by the viewfinder and sets the appropriate exposure value. Considering that the laser intensity data is light intensity, we are adopting the concept of metering system in that it is the same as iso-value selection and camera exposure value selection. In Fig. 3, each camera is slightly different, but the basic metering system is spot, center-weighted, and matrix. We analyze these three methods and consider appropriate methods for analyzing laser intensity data [8, 9].

Fig. 3.
figure 3

Three kinds of metering systems on nikon camera: (a) spot, (b) center-weighted, and (c) matrix

3.1 Spot Metering

Spot metering is a method of measuring the center of the area corresponding to about 3.5% of the area in the viewfinder (Fig. 4). This method is suitable for cases where the illumination difference between the background and the subject is large because only a specific portion is precisely measured irrespective of the overall brightness.

Fig. 4.
figure 4

Spot metering in the viewfinder (left) and output (right)

3.2 Center-Weighted Metering

Center-weighted metering focuses on the center of the viewfinder and then measures the overall average brightness. As shown in Fig. 5, the measured value of the central area corresponding to 20~30% of the viewfinder is weighted by 60~80%, and the measured value of the other area is weighted by 20~40%. Therefore, it is the most suitable measurement method when the subject is in the middle.

Fig. 5.
figure 5

Center-weighted metering in the viewfinder (left) and output (right)

3.3 Matrix Metering

Matrix metering divides the entire screen into zones and then analyzes the exposure patterns of each zone through the built-in algorithm of each manufacturer. It then determines the exposure value of the average of all the zones. Different manufacturers have different exposure patterns, and Nikon has a database that contains more than 30,000 referenceable exposure information (Fig. 6).

Fig. 6.
figure 6

Matrix metering in the viewfinder (left) and output (right)

4 Laser Intensity Metering System

The algorithm used in our proposed laser intensity metering system is based on algorithm embedded in MatLab. First, we convert the volume data to a histogram of N-bins, the range of values. Next, we average the remaining bins by removing mode value lower than the median, which are most likely to be a garbage. We used this filtering algorithm to determine the iso-values of the following measurement methods; entire volume metering, spot metering, center-weighted metering [10, 11].

4.1 Entire Volume Metering

Figure 7 shows the selected region in the entire volume data and the histogram of its region. As can be seen from Fig. 7, the bin count to the minimum value occupies the largest portion. Figure 8 shows the average values before and after removing the mode value. By eliminating the bin that is close to zero and occupies the most amount, we can see that the result is a more meaningful average.

Fig. 7.
figure 7

Histogram of the entire volume data

Fig. 8.
figure 8

Comparison of iso-value before and after removal of mode data

4.2 Spot Metering

Figure 9 shows histograms in three randomly selected spots, which are 3.5% of the total volume. As you can see from the obvious difference between Spot 2 and the other two spots, the distribution of the values depends on the selected area. Furthermore, the distribution of the values of spots 1 and 3 is similar at first glance, but the magnitude of the maximum intensity value present in the area is different. Thus, the filtered iso-value is also different for each zone (Fig. 10).

Fig. 9.
figure 9

Histograms of the selected spots

Fig. 10.
figure 10

Comparison of iso-value in each spot

4.3 Center-Weighted Metering

Figure 11 depicts about 25% of the entire volume and the histogram of selected area. The unweighted central area is equivalent to widening the area of spot 2 in Fig. 9. The central area has a relatively stable form of the distribution value as compared to the remaining area, but has relatively low intensity values.

Fig. 11.
figure 11

Histogram of the selected central data

In the center-weighted metering method, which is a mixed method of entire metering and spot metering, we use weighting method to reduce the difference in results depending on the selected position. The weighted average equation we used is:

$$ i_{avg} = \frac{{\mathop \sum \nolimits_{i = 1}^{n} w_{i} \times i_{i} }}{{\mathop \sum \nolimits_{i = 1}^{n} w_{i} }} $$

In this equation \( w_{i} \) is the weight given in each area, and \( i_{i} \) is the iso-value of each area. We used this equation to assign a weight of 75% to the central area and 25% to the other areas. The result is shown in Fig. 12.

Fig. 12.
figure 12

Comparison of iso-value before and after weighting

5 Visualization

In an experiment to visualize the iso-surface of the laser intensity data, we used VisIt. VisIt was developed by LLNL (Lawrence Livermore National Laboratory) and is an open source visualization tool for effectively representing science and engineering data by BSD licensing. In addition, VisIt is capable of visualizing data from over 120 different scientific data formats including HDF5 [12].

Table 2 shows the iso-values used in the experiment. The iso-values required to perform iso-surface were obtained through the methods described above; entire volume metering, spot metering and center-weighted metering; and those obtained through the built-in algorithm in MatLab.

Table 2. Iso-value of each algorithm

Figure 13 is the result of surface rendering the laser intensity data through the iso-surface. Figure 13(e) and (f) show larger surfaces because they have lower iso-values than Fig. 13(a). In addition, the results of Fig. 13(b), (c), and (d) show that the surfaces of the laser intensity data are broken away in the form of islands.

Fig. 13.
figure 13

The result of surface rendering by iso-values of Table 2

6 Conclusions

In this paper, we used iso-surface as a representative surface visualization method for geometric analysis to visualize laser intensity data. However, due to the properties of the scattered volumetric data of the laser intensity data, it is difficult to select appropriate iso-values for the iso-surface. Therefore, we adopted the camera metering system; A method of providing an average value by removing a value close to garbage in the entire region or selected spot, and a method of providing center-weighted average value. Through these methods, it is possible to suggest a more meaningful iso-value to the user. For further research, we are trying to provide a more appropriate iso-value by matrix pattern metering that reflects the user’s interest through machine learning.

.