Copper-based metal halides for X-ray and photodetection

Abstract

Copper-based metal halides have become important materials in the field of X-ray and photodetection due to their excellent optical properties, good environmental stability and low toxicity. This review presents the progress of research on crystal structure/morphology, photophysics/optical properties and applications of copper-based metal halides. We also discuss the challenges of copper-based metal halides with a perspective of their future research directions.

Graphical abstract

Keywords

Copper-based metal halides / X-ray detector / Photodetectors / Scintillators

Introduction

Precise calibration for a camera is fundamental in computer vision and vision metric. A camera is often modeled as an ideal pin-hole one, which images on the focus plane without distortion. However, lens distortion will make the image distorted from ideal one.

There are several distortion models depending on the type of lens. In this paper, radial, decentering and prism distortion models proposed by Brown [ 1, 2] were to be considered.

To identify each distortion model, several methods have been proposed. The first method uses the known object 3D world coordinates [ 3], with which the result could be inaccurate when both camera intrinsic and extrinsic parameters are estimated at the same time. The second method uses point correspondences in different views [ 4]. This method is based on fundament matrix without the intrinsic and extrinsic explicitly. However, it is difficult to find correct point correspondences in multi-views. The third method, which estimates the distortion without camera parameters, is based on projective invariants [ 5- 7]. Nowadays, more and more researches depend on this method. This method estimates the distortion without camera parameters. The most widely used projective invariant is straight line, which remains straight from different view if there is no lens distortion based on pin-hole camera. This method needs to know the scene in advance, as it is important to find out the straight line or other invariant objects in distortion image. The forth method is also based on point correspondences between image points and world points [ 8], and it uses planar invariant instead. Planar calibration pattern points and non-distortion image points are mapped by homography matrix. This method calculates the distortion model parameters to make the correspondence most fit.

There are two different ways to compute the distortion models, backward mapping and forward mapping. Zhang [ 4] used forward mapping in both model definition and computing, while Hartley and Kang [ 9] used backward mapping.

In this paper, a new camera calibration method is presented, we used only one chessboard pattern image to solve the distortion models, which includes three main types of lens distortion–radial, decentering and prism distortions. It has four steps. First, image feature points are detected by chessboard point detecting algorithm. Second, the non-distortion feature point coordinate is calculated from detected feature point on distortion image. Third, nonlinear optimization is done to get the distortion model. Last, iterate former the three steps until converge. This method is based on the idea that the distortion near the optic axis is tiny, and we can use iterative procedure to reduce influence of this distortion. Therefore, non-distortion feature points can be calculated using the near center feature points on distortion image.

The paper is organized as follows. Section 2 describes the camera lens distortion model. Section 3 proposes the method of distortion model calculating. In Section 4, several experimental results are reported on both real and synthetic data. The paper ends with some concluding remarks.

Lens distortion model

Lens distortion makes the actual image different from ideal pin-hole image. Three forms of lens distortion, namely radial distortion, decentering distortion and thin prism distortion are taken into account in this paper [ 8]. The distortion can be formed as follow:

(1) $$\left[\begin{array}{l}x\\ y\end{array}\right]=\left[\begin{array}{l}{x}^{*}\\ y^{*}\end{array}\right]+\left[\begin{array}{l}{\mathrm{\Delta }}_{x^{*}}\\ {\mathrm{\Delta }}_{y^{*}}\end{array}\right].$$

${[x,y]}^T$ is the ideal image without distortion, ${[x^{*},y^{*}]}^T$ is the actual image with distortion, ${[{\mathrm{\Delta }}_{x^{*}},{\mathrm{\Delta }}_{y^{*}}]}^T$ is the amount of distortion. Different forms of lens distortion result in the different part of distortion factor.

Radial distortion

Radial distortion is caused by inconsistent magnification in the field of view. This distortion is the significant part of lens distortion. Many researches only dealt with radial distortion [ 5, 7, 10].

The amount of radial distortion is circularly symmetric about optical center as shown in Fig. 1, which means there is no distortion at the optical center. In general, the farther the image point is away from the optical axis, the greater the distortion is. Radial distortion can be expressed as follows:

(2) $$\{\begin{array}{l}{\mathrm{\Delta }}_{rx}(x^{*},y^{*})=k_{\mathrm{1}}{x}^{*}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})+k_{\mathrm{2}}{x}^{*}{({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})}^{\mathrm{2}}\\ {\mathrm{\Delta }}_{ry}(x^{*},y^{*})=k_{\mathrm{1}}{y}^{*}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})+k_{\mathrm{2}}{y}^{*}{({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})}^{\mathrm{2}}\end{array},$$

where k₁, k₂ are distortion factors of radial distortion model.

Fig.1 (a) Barrel distortion image; (b) pincushion distortion image

Full size|PPT slide

Decentering distortion

Decentering distortion is due to non-strict collineation of the optical centers of lens elements, as shown in Fig. 2. The decentering distortion model is formed as Eq. (3).

(3) $$\{\begin{array}{l}{\mathrm{\Delta }}_{dx}(x^{*},y^{*})=p_{\mathrm{1}}(\mathrm{3}{(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})+\mathrm{2}{p}_{\mathrm{2}}{x}^{*}{y}^{*}\\ {\mathrm{\Delta }}_{dy}(x^{*},y^{*})=p_{\mathrm{2}}({(x^{*})}^{\mathrm{2}}+\mathrm{3}{(y^{*})}^{\mathrm{2}})+\mathrm{2}{p}_{\mathrm{1}}{x}^{*}{y}^{*}\end{array},$$

where p₁, p₂ are distortion factors of decentering distortion model.

Fig.2 Image with decentering distortion

Full size|PPT slide

Thin prism distortion

Thin prism distortion arises from the imperfection in lens design and manufacturing as well as camera, as shown in Fig. 3.The thin prism distortion model is described in Eq. (4).

(4) $$\{\begin{array}{l}{\mathrm{\Delta }}_{px}(x^{*},y^{*})=s_{\mathrm{1}}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})\\ {\mathrm{\Delta }}_{py}(x^{*},y^{*})=s_{\mathrm{2}}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})\end{array},$$

where s₁, s₂ are distortion factors of prism distortion model.

Fig.3 Image with thin prism distortion

Full size|PPT slide

In this paper, three forms of distortion are considered. Then, the total amount of distortion is the sum of three distortions.

(5) $$\{\begin{array}{l}{\mathrm{\Delta }}_{x^{*}}(x^{*},y^{*})={\mathrm{\Delta }}_{rx}(x^{*},y^{*})+{\mathrm{\Delta }}_{dx}(x^{*},y^{*})+{\mathrm{\Delta }}_{px}(x^{*},y^{*})\\ {\mathrm{\Delta }}_{y^{*}}(x^{*},y^{*})={\mathrm{\Delta }}_{ry}(x^{*},y^{*})+{\mathrm{\Delta }}_{dy}(x^{*},y^{*})+{\mathrm{\Delta }}_{py}(x^{*},y^{*})\end{array}.$$

Distortion correct with one chessboard pattern image

The proposed method consists following steps: 1) preprocess chessboard image and pick up the coordinate of corner points; 2) reconstruct non-distortion corner coordinate from distortion image; 3) initialize distortion model parameters using Levenberg-Marquardt nonlinear algorithm (LM); 4) use iterative procedure to minimize distortion residual; 5) use forward mapping method to correct image.

Corner points detection on distortion image

The performance of distortion correction is affected by the accuracy of the corner point coordinate. In this paper, a new corners detecting method of chessboard pattern image is proposed. This method has high corner detection rate.

The proposed method detects corner points according to the dissimilar distribution of pixels gray-scale between corner region and non-corner region. It defines three operators to calculate the dissimilarity distribution on different directions as follows. S is centrosymmetric operator. V is vertical symmetry operator. H is horizontal symmetry operator.

(6) $$\{\begin{array}{l}s(x_{\mathrm{0}},y_{\mathrm{0}})={\displaystyle \sum _{(x,y)\in neighborhood(x_{\mathrm{0}},y_{\mathrm{0}})}\left|p_{(x,y)}-p_{(\mathrm{2}{x}_{\mathrm{0}}-x,\mathrm{2}{y}_{\mathrm{0}}-y)}\right|}/card(neighborhood(x_{\mathrm{0}},y_{\mathrm{0}}))\\ v(x_{\mathrm{0}},y_{\mathrm{0}})={\displaystyle \sum _{(x,y)\in neighborhood(x_{\mathrm{0}},y_{\mathrm{0}})}\left|p_{(x,y)}-p_{(x,\mathrm{2}{y}_{\mathrm{0}}-y)}\right|}/card(neighborhood(x_{\mathrm{0}},y_{\mathrm{0}}))\\ h(x_{\mathrm{0}},y_{\mathrm{0}})={\displaystyle \sum _{(x,y)\in neighborhood(x_{\mathrm{0}},y_{\mathrm{0}})}\left|p_{(x,y)}-p_{(\mathrm{2}{x}_{\mathrm{0}}-x,y)}\right|}/card(neighborhood(x_{\mathrm{0}},y_{\mathrm{0}}))\end{array},$$

where $neighborhood(x_{\mathrm{0}},y_{\mathrm{0}})$ is 8 × 8rectangle region with point $(x_{\mathrm{0}},y_{\mathrm{0}})$ located at the center, card() returns the pixel number of the region.

Define M as pixels response, whose value is associated with three operators. The response value of corner point is calculated by the following steps:

1)Compute the value of pixels by the three operators;

2)If all of the three operator results are between the pre-defined ranges, the final response value is the result by operator S. Otherwise, M is set to zero.

(7) $$M=\{s|s\in [s_{\min },s_{\max }],v\in [v_{\min },v_{\max }],h\in [h_{\min },h_{\max }]\}\cup \{\mathrm{0}|s\notin [s_{\min },s_{\max }]orv\notin [v_{\min },v_{\max }]orh\notin [h_{\min },h_{\max }]\},$$

Median filtering algorithm is used to remove isolated detective point which is fault detection caused by noise. The corner coordinate can be computed by calculating the connected region centroid of response M, described as follows:

(8) $$centroid={\displaystyle \frac{{\displaystyle \sum _{connected}m*position}}{{\displaystyle \sum _{connected}m}}},$$

where $m\in M$ , connected region is detected by 8-adjacent connection, position_m is the coordinate of the point with value m.

In this paper, distortion of the four corner points around the center is tiny for distortion image. Non-distortion points are reconstructed by the geometry feature (considering the same size of each block) of chessboard, which exploits these four corner points, as shown in Fig. 4, defined in Eq. (9).

(9) $$\{\begin{array}{l}{x}_{(i,j)}=x_{(\mathrm{0},\mathrm{0})}+j*(x_{(\mathrm{0},\mathrm{1})}-x_{(\mathrm{0},\mathrm{0})})+i*(x_{(\mathrm{1},\mathrm{0})}-x_{(\mathrm{0},\mathrm{0})})\hfill \\ y_{(i,j)}=y_{(\mathrm{0},\mathrm{0})}+j*(y_{(\mathrm{0},\mathrm{1})}-y_{(\mathrm{0},\mathrm{0})})+i*(y_{(\mathrm{1},\mathrm{0})}-y_{(\mathrm{0},\mathrm{0})})\hfill \end{array},$$

where (i,j) is the points array order, i and j could be negative integer, center four points are (0,0), (0,1), (1,0), (1,1).

To minimize the influence of the distortion of four center points, this paper introduces iterative process. The distortion of the four points would reduce after image correction, and then the corrected image is treated as an input distortion image for next iteration to reconstruct non-distortion points, iterate this procedure until four points unchanged after correction.

Fig.4 (a) Distortion image; (b) the dots are detected corners in image (a) and the circles are ideal points location.

Full size|PPT slide

Distortion model calculate

Three forms of distortion indicate that the optical center region is always non-distortion.

As long as the distortion and non-distortion points are achieved, the total amount of distortion would be computed by Eq. (1).

Expansion Eq. (5) is the whole model formed by three distortions.

(10) $$\{\begin{array}{l}{\mathrm{\Delta }}_x(x^{*},y^{*})=k_{\mathrm{1}}{x}^{*}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})+k_{\mathrm{2}}{x}^{*}{({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})}^{\mathrm{2}}+[p_{\mathrm{1}}(\mathrm{3}{(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})+\mathrm{2}{p}_{\mathrm{2}}{x}^{*}{y}^{*}]+s_{\mathrm{1}}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})\\ {\mathrm{\Delta }}_y(x^{*},y^{*})=k_{\mathrm{1}}{y}^{*}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})+k_{\mathrm{2}}{y}^{*}{({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})}^{\mathrm{2}}+[p_{\mathrm{2}}({(x^{*})}^{\mathrm{2}}+\mathrm{3}{(y^{*})}^{\mathrm{2}})+\mathrm{2}{p}_{\mathrm{1}}{x}^{*}{y}^{*}]+s_{\mathrm{2}}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})\end{array}.$$

In this paper, the optical center is set to the image center. Different center only affects the distortion model values, which would not influence the result of distortion correction. Then, the linear equation is showed as follows

(11) $$X-X^{\mathit{*}}=AP,$$

where

$$\{\begin{array}{l}A=\left[\begin{array}{l}{x}^{*}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}}),x^{*}{({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})}^{\mathrm{2}},\mathrm{3}{(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}},\mathrm{2}{x}^{*}{y}^{*},{(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}},\mathrm{0}\\ y^{*}({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}}),y^{*}{({(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}})}^{\mathrm{2}},\mathrm{2}{x}^{*}{y}^{*},{(x^{*})}^{\mathrm{2}}+\mathrm{3}{(y^{*})}^{\mathrm{2}},\mathrm{0},{(x^{*})}^{\mathrm{2}}+{(y^{*})}^{\mathrm{2}}\end{array}\right]\\ P=\mathit{}{\left[\begin{array}{llllll}{k}_{\mathrm{1}}& k_{\mathrm{2}}& p_{\mathrm{1}}& p_{\mathrm{2}}& s_{\mathrm{1}}& s_{\mathrm{2}}\end{array}\right]}^T\\ X=\mathit{}{[x\hspace{1em}y]}^T\\ X^{\mathit{*}}={[x^{*}\hspace{1em}{y}^{*}]}^T\end{array}.$$

P is the required model parameters, X is the non-distortion position, X* is the distortion position.

The model P can be calculated by LM nonlinear optimization [ 11]. The objective function is defined as follows:

(12) $${\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}[{\Vert X_{\mathrm{0}}-X(P,X_{\mathrm{0}}^{*})\Vert }^{\mathrm{2}},\cdots ,{\Vert X_j-X(P,X_j^{*})\Vert }^{\mathrm{2}},\cdots ,{\Vert X_m-X(P,X_m^{*})\Vert }^{\mathrm{2}}]=\mathrm{0},$$

where $X_j$ is non-distortion pixel coordinate, $X(P,X_j^{*})$ is the pixel coordinate transferred by Eq. (11) related to distortion image.

The Jacobian matrix is used for LM algorithm. From Eq. (12), the Jacobian matrix can be easy to get.

Then, the procedure for calculating distortion parameters is:

1)use corner detection to obtain point coordinate;

2)calculate the non-distortion point location depending on near center point;

3)using LM method to optimize distortion model parameters;

4)calculate the distortion residual, repeat steps 1), 2) and 3) until distortion residual convergence;

The parameter of convergence is the best solution of proposed method.

Distortion image correction

Image is corrected by the distortion parameters obtained in last steps. This procedure contains two steps.

1)Coordinate transfer: transfer the pixel coordinate from distortion image to non-distortion image (forward-mapping, as shown in Fig. 5) or from non-distortion image to distortion image (backward-mapping). This paper uses forward-mapping method to rectify image.

Fig.5 Forward-mapping from distortion image to correction image

Full size|PPT slide

2)Gray-scale interpolation: in this paper, bilinear interpolation is used to reconstruct the correction image.

Experiments

Corner points detection algorithm

The performance of the proposed corner detective algorithm is tested by comparing with the chessboard corner detection algorithm in OpenCV. It is shown in Fig. 6 that the proposed method achieves higher detection rate and is more accurate than the OpenCV method.

Fig.6 White rectangles are the detected corners: (a) the proposed method; (b) the OpenCV method

Full size|PPT slide

Synthetic data experiment

The performance of the proposed method is tested in several experiments. First, use the proposed method to estimate parameters of artificial distortion images. We distorted the origin image (Fig. 7) by adding all the three forms of distortion, and then corrected them with the proposed method.

Fig.7 Origin image without distortion

Full size|PPT slide

Fig.8 (a) and (b) are images with different distortion parameters; (c) and (d) are corrected image by proposed method, which parameters are shown in Table 1

Full size|PPT slide

Tab.1 Compares the distortion between man-made parameters and the results of proposed method

origin distortion image	Figure 8(a)	Figure 8 (b)
distortion model	$k_{\mathrm{1}}=\mathrm{3.0}e-\mathrm{6},k_{\mathrm{2}}=-\mathrm{5.8}e-\mathrm{13}$ $p_{\mathrm{1}}=-\mathrm{2.4}e-\mathrm{5},p_{\mathrm{2}}=\mathrm{2.86}e-\mathrm{5}$ $s_{\mathrm{1}}=\mathrm{6.5}e-\mathrm{6},s_{\mathrm{2}}=-\mathrm{7.2}e-\mathrm{6}$	$k_{\mathrm{1}}=-\mathrm{2.3}e-\mathrm{6},k_{\mathrm{2}}=\mathrm{4.3}e-\mathrm{13}$ $p_{\mathrm{1}}=\mathrm{5.0}e-\mathrm{5},p_{\mathrm{2}}=-\mathrm{3.4}e-\mathrm{5}$ $s_{\mathrm{1}}=\mathrm{7.8}e-\mathrm{6},s_{\mathrm{2}}=\mathrm{5.6}e-\mathrm{6}$
correct image by proposed method	Figure 8(c)	Figure 8(d)
correct model by proposed method	$k_{\mathrm{1}}=\mathrm{3.67}e-\mathrm{6},k_{\mathrm{2}}=-\mathrm{3.47}e-\mathrm{12}$ $p_{\mathrm{1}}=-\mathrm{1.99}e-\mathrm{5},p_{\mathrm{2}}=\mathrm{2.36}e-\mathrm{5}$ $s_{\mathrm{1}}=-\mathrm{1.26}e-\mathrm{6},s_{\mathrm{2}}=\mathrm{2.87}e-\mathrm{7}$	$k_{\mathrm{1}}=-\mathrm{2.67}e-\mathrm{6},k_{\mathrm{2}}=\mathrm{1.37}e-\mathrm{13}$ $p_{\mathrm{1}}=\mathrm{4.90}e-\mathrm{5},p_{\mathrm{2}}=-\mathrm{2.79}e-\mathrm{5}$ $s_{\mathrm{1}}=\mathrm{1.48}e-\mathrm{5},s_{\mathrm{2}}=-\mathrm{6.87}e-\mathrm{6}$

This method has a good estimation about parameters $k_{\mathrm{1}},p_{\mathrm{1}},p_{\mathrm{2}}$ , which are the main factors of distortion. The correction images are very similar to the origin image without distortion, shown in Fig. 8.

Synthetic data are used for further analysis for the performance of the proposed method. The synthetic data are setup as Gao and Yin did inRef [ 8], shown in Fig. 9(a), defining the optical axis as z-axis, the horizon direction on the image plane as x-axis, and vertical direction on the image plane as y-axis.

Fig.9 Performance of the proposed method with the synthetic data, dots are distorted points and circles are non-distorted points. (a) The distorted image; (b) the corrected image by the proposed method

Full size|PPT slide

Fig. 10 shows the three methods’ performance versus noises. In this experiment, Gaussian noise with zero mean and σ standard deviation is added to distorted points. We test the performance with noise standard deviation σ from 0.1 to 1.5 pixels. Define ERROR as root of mean square error in Eq. (13), which is to describe the correction accuracy. For each noise level, we used 100 images to evaluate the method’s performance. Brown’s method does not correct prism distortion. It indicates that the method without considering the entire distortion model could not make a good image correction.

(13) $$ERROR={\displaystyle \frac{{\displaystyle \sum _i\sqrt{{(x_d-x_r)}^{\mathrm{2}}+{(y_d-y_r)}^{\mathrm{2}}}}}{n}},$$

where n is the number of samples, (x_d,y_d) is test data and (x_r,y_r) is the real data.

Fig.10 Compare ERROR of the proposed method with Gao’s and Brown’s method versus different noises levels

Full size|PPT slide

To test the robust of the proposed method, the image data are rotated around z-axis or x-axis. The results are shown in Figs. 11 and 12. Figure 11 shows the ERROR results when the chessboard is rotated around z-axis. It shows that the correction accuracy is similar between different rotation angles. In Fig. 12, though the rotation around x-axis is even up to 7.5°, the ERROR is similar with Gao’s method.

Fig.11 ERROR of different rotations around the z-axis of image versus noises levels

Full size|PPT slide

Fig.12 ERROR versus different rotation around the x-axis of image

Full size|PPT slide

The proposed method is very fast to converge. In these experiments, the iteration times are nearly 2 to 4, and the total time consuming to solve distortion model is almost 10 ms (image size is 752 × 480, using PC with CPU of i3 2.5 GHz).

Real data experiment

Real distortion images are taken by the Pentax CCTV lens (F1.2) with the focal lengths of 4 and 8 mm, which are captured by smart camera based on TMS320DM6437 DSP. Figure 7 is the object for taking photos. Figure 13 shows the pictures with different lens.

Fig.13 (a) Pictures taken with 4 mm lens, image size is 752 ×480; (b) pictures is taken with 8 mm lens, image size is 640 ×480

Full size|PPT slide

Use the proposed method to correct the distortion images in Fig. 13, and results are shown in Fig. 14. It shows that proposed method makes the collinear point in straight pattern in image. Table 2 presents the solutions in detail. We evaluate distortion residual of the distance between point and line. Since the points on each columns and rows are actually collinear, the distortion residual is the average distance between points and line. First, least square method is used to estimate the line function. Second, the absolute distance from points to the line is computed. The average absolute distance of all the point indicates the error of distortion, named as root mean square error (RMSE). Maximum distortion rate is calculated by the maximum distance from point to line dividing the image diagonal length.

Fig.14 (a) Corrected image taken with 4 mm lens; (b) corrected image taken with 8 mm lens

Full size|PPT slide

Tab.2 Experiment results with proposed method

	len 1 focal f = 4 mm	len 2 focal f = 8 mm
correct model	$k_{\mathrm{1}}=\mathrm{8.86}e-\mathrm{7},k_{\mathrm{2}}=\mathrm{3.97}e-\mathrm{12}$ $p_{\mathrm{1}}=\mathrm{1.23}e-\mathrm{5},p_{\mathrm{2}}=\mathrm{5.64}e-\mathrm{5}$ $s_{\mathrm{1}}=\mathrm{1.43}e-\mathrm{5},s_{\mathrm{2}}=-\mathrm{2.23}e-\mathrm{5}$	$k_{\mathrm{1}}=\mathrm{7.72}e-\mathrm{7},k_{\mathrm{2}}=\mathrm{1.98}e-\mathrm{14}$ $p_{\mathrm{1}}=\mathrm{1.05}e-\mathrm{8},p_{\mathrm{2}}=\mathrm{2.89}e-\mathrm{8}$ $s_{\mathrm{1}}=\mathrm{1.37}e-\mathrm{8},s_{\mathrm{2}}=\mathrm{2.15}e-\mathrm{8}$
RMSE	0.34 pixel	0.25 pixel
maximum distortion rate	0.13%	0.10%

Table 2 indicates that proposed method is effective, which makes the image with little distortion and even the maximum distortion rather small.

Conclusions

This paper proposed a new method for lens distortion correction. Only one chessboard pattern image was used to rectify the three distortion models, which makes the implement convenient. To minimize the influence by the tiny distortion of the center four points, iterative procedure was introduced to optimize the distortion model and reevaluate the non-distorted points. Both synthetic data and real data have been adopted to test by the proposed method. Synthetic data experiment showed that the proposed method had good performance versus noises, and it was more accurate for lens distortion correction by considering the whole distortion models. Besides, chessboard rotation experiment was conducted to show the robustness of this method. The result of real data experiment also indicated that the proposed method was effective.

References

[1] Liang, J., Liu, J., Jin, Z.: All-inorganic halide perovskites for optoelectronics: progress and prospects. Solar RRL 1(10), 1700086(2017)

[2] Xiang, W., Tress, W.: Review on recent progress of all-inorganic metal halide perovskites and solar cells. Adv. Mater. 31(44), e1902851 (2019)

[3] Chen, W., Li, X., Li, Y., Li, Y.: A review: Crystal growth for highperformance all-inorganic perovskite solar cells. Energy Environ. Sci. 13(7), 1971–1996 (2020)

[4] Yuan, J., Hazarika, A., Zhao, Q., Ling, X., Moot, T., Ma, W., Luther, J.M.: Metal halide perovskites in quantum dot solar cells: progress and prospects. Joule 4(6), 1160–1185 (2020)

[5] Liu, P., Han, N., Wang, W., Ran, R., Zhou, W., Shao, Z.: High-quality ruddlesden-popper perovskite film formation for high-performance perovskite solar cells. Adv. Mater. 33(10), e2002582 (2021)

[6] Park, N.G.: Perovskite solar cells: an emerging photovoltaic technology. Mater. Today 18(2), 65–72 (2015)

[7] Leijtens, T., Bush, K.A., Prasanna, R., Mcgehee, M.D.: Opportunities and challenges for tandem solar cells using metal halide perovskite semiconductors. Nat. Energy 3(10), 828–838 (2018)

[8] Wu, T., Qin, Z., Wang, Y., Wu, Y., Chen, W., Zhang, S., Cai, M., Dai, S., Zhang, J., Liu, J., Zhou, Z., Liu, X., Segawa, H., Tan, H., Tang, Q., Fang, J., Li, Y., Ding, L., Ning, Z., Qi, Y., Zhang, Y., Han, L.: The main progress of perovskite solar cells in 2020–2021. Nano-Micro Lett. 13(1), 152(2021)

[9] Li, B., Li, Z., Wu, X., Zhu, Z.: Interface functionalization in inverted perovskite solar cells: from material perspective. Nano Res Energy 1, e9120011 (2022)

[10] Tan, Z.K., Moghaddam, R.S., Lai, M.L., Docampo, P., Higler, R., Deschler, F., Price, M., Sadhanala, A., Pazos, L.M., Credgington, D., Hanusch, F., Bein, T., Snaith, H.J., Friend, R.H.: Bright light-emitting diodes based on organometal halide perovskite. Nat. Nanotechnol. 9(9), 687–692 (2014)

[11] Liu, M., Wan, Q., Wang, H., Carulli, F., Sun, X., Zheng, W., Kong, L., Zhang, Q., Zhang, C., Zhang, Q., Brovelli, S., Li, L.: Suppression of temperature quenching in perovskite nanocrystals for efficient and thermally stable light-emitting diodes. Nat. Photonics 15(5), 379–385 (2021)

[12] Ji, K., Anaya, M., Abfalterer, A., Stranks, S.D.: Halide perovskite light-emitting diode technologies. Adv. Opt. Mater. 9(18), 2002128(2021)

[13] Dou, L., Yang, Y.M., You, J., Hong, Z., Chang, W.H., Li, G., Yang, Y.: Solution-processed hybrid perovskite photodetectors with high detectivity. Nat. Commun. 5(1), 5404(2014)

[14] Ramasamy, P., Lim, D.H., Kim, B., Lee, S.H., Lee, M.S., Lee, J.S.: All-inorganic cesium lead halide perovskite nanocrystals for photodetector applications. Chem. Commun. (Camb.) 52(10), 2067–2070 (2016)

[15] Wang, H.P., Li, S., Liu, X., Shi, Z., Fang, X., He, J.H.: Low-dimensional metal halide perovskite photodetectors. Adv. Mater. 33(7), e2003309 (2021)

[16] Li, Z., Peng, G., Chen, H., Shi, C., Li, Z., Jin, Z.: Metal-free PAZE-NH₄X₃·H₂O perovskite for flexible transparent X-ray detection and imaging. Angew. Chem. Int. Ed. 61(36), 202207198(2022)

[17] Chen, Q., Wu, J., Ou, X., Huang, B., Almutlaq, J., Zhumekenov, A.A., Guan, X., Han, S., Liang, L., Yi, Z., Li, J., Xie, X., Wang, Y., Li, Y., Fan, D., Teh, D.B.L., All, A.H., Mohammed, O.F., Bakr, O.M., Wu, T., Bettinelli, M., Yang, H., Huang, W., Liu, X.: All-inorganic perovskite nanocrystal scintillators. Nature 561(7721), 88–93 (2018)

[18] Cao, F., Yu, D., Ma, W., Xu, X., Cai, B., Yang, Y.M., Liu, S., He, L., Ke, Y., Lan, S., Choy, K.L., Zeng, H.: Shining emitter in a stable host: design of halide perovskite scintillators for X-ray imaging from commercial concept. ACS Nano 14(5), 5183–5193 (2020)

[19] Zhu, W., Ma, W., Su, Y., Chen, Z., Chen, X., Ma, Y., Bai, L., Xiao, W., Liu, T., Zhu, H., Liu, X., Liu, H., Liu, X., Yang, Y.M.: Low-dose real-time X-ray imaging with nontoxic double perovskite scintillators. Light Sci. Appl. 9(1), 112(2020)

[20] Chen, H., Wang, Q., Peng, G., Wang, S., Lei, Y., Wang, H., Yang, Z., Sun, J., Li, N., Zhao, L., Lan, W., Jin, Z.: Cesium lead halide nanocrystals based flexible X-ray imaging screen and visible dose rate indication on paper substrate. Adv. Opt. Mater. 10(8), 2102790(2022)

[21] Dong, H., Zhang, C., Liu, X., Yao, J., Zhao, Y.S.: Materials chemistry and engineering in metal halide perovskite lasers. Chem. Soc. Rev. 49(3), 951–982 (2020)

[22] Zhang, Q., Shang, Q., Su, R., Do, T.T.H., Xiong, Q.: Halide perovskite semiconductor lasers: materials, cavity design, and low threshold. Nano Lett. 21(5), 1903–1914 (2021)

[23] Li, Z., Zhou, F., Yao, H., Ci, Z., Yang, Z., Jin, Z.: Halide perovskites for high-performance X-ray detector. Mater. Today 48, 155–175 (2021)

[24] Zhou, F., Li, Z., Lan, W., Wang, Q., Ding, L., Jin, Z.: Halide perovskite, a potential scintillator for X-ray detection. Small Methods 4(10), 2000506(2020)

[25] Krishnamoorthy, T., Ding, H., Yan, C., Leong, W.L., Baikie, T., Zhang, Z., Sherburne, M., Li, S., Asta, M., Mathews, N., Mhaisalkar, S.G.: Lead-free germanium iodide perovskite materials for photovoltaic applications. J. Mater. Chem. A Mater. Energy Sustain. 3(47), 23829–23832 (2015)

[26] Yu, B.B., Chen, Z., Zhu, Y., Wang, Y., Han, B., Chen, G., Zhang, X., Du, Z., He, Z.: Heterogeneous 2D/3D tin-halides perovskite solar cells with certified conversion efficiency breaking 14%. Adv. Mater. 33(36), e2102055 (2021)

[27] Jiang, F., Yang, D., Jiang, Y., Liu, T., Zhao, X., Ming, Y., Luo, B., Qin, F., Fan, J., Han, H., Zhang, L., Zhou, Y.: Chlorine-incor-poration- induced formation of the layered phase for antimony-based lead-free perovskite solar cells. J. Am. Chem. Soc. 140(3), 1019–1027 (2018)

[28] Leng, M., Yang, Y., Zeng, K., Chen, Z., Tan, Z., Li, S., Li, J., Xu, B., Li, D., Hautzinger, M.P., Fu, Y., Zhai, T., Xu, L., Niu, G., Jin, S., Tang, J.: All-inorganic bismuth-based perovskite quantum dots with bright blue photoluminescence and excellent stability. Adv. Funct. Mater. 28(1), 1704446(2018)

[29] Li, M., Li, F., Gong, J., Zhang, T., Gao, F., Zhang, W.H., Liu, M.: Advances in TiN(II)-based perovskite solar cells: from material physics to device performance. Small Struct. 3(1), 2100102(2022)

[30] Tang, Y., Tang, S., Luo, M., Guo, Y., Zheng, Y., Lou, Y., Zhao, Y.: All-inorganic lead-free metal halide perovskite quantum dots: progress and prospects. Chem. Commun. (Camb.) 57(61), 7465–7479 (2021)

[31] Jun, T., Sim, K., Iimura, S., Sasase, M., Kamioka, H., Kim, J., Hosono, H.: Lead-free highly efficient blue-emitting Cs₃Cu₂I₅ with 0D electronic structure. Adv. Mater. 30(43), e1804547 (2018)

[32] Cao, L., Liu, X., Li, Y., Li, X., Du, L., Chen, S., Zhao, S., Wang, C.: Recent progress in all-inorganic metal halide nanostructured perovskites: materials design, optical properties, and application. Front. Phys. 16(3), 33201(2021)

[33] Hull, S., Berastegui, P.: Crystal structures and ionic conductivities of ternary derivatives of the silver and copper monohalides—II: ordered phases within the (AgX)_x(MX)_1–x and (CuX)_x(MX)_1–x (M=K, Rb and Cs; X=Cl, Br and I) systems. J. Solid State Chem. 177(9), 3156–3173 (2004)

[34] Li, Y., Zhou, Z., Tewari, N., Ng, M., Geng, P., Chen, D., Ko, P.K., Qammar, M., Guo, L., Halpert, J.E.: Progress in copper metal halides for optoelectronic applications. Mater. Chem. Front. 5(13), 4796–4820 (2021)

[35] Grandhi, G.K., Viswanath, N.S.M., Cho, H.B., Han, J.H., Kim, S.M., Choi, S., Im, W.B.: Mechanochemistry as a green route: Synthesis, thermal stability, and postsynthetic reversible phase transformation of highly-luminescent cesium copper halides. J. Phys. Chem. Lett. 11(18), 7723–7729 (2020)

[36] Lin, R., Guo, Q., Zhu, Q., Zhu, Y., Zheng, W., Huang, F.: All-inorganic CsCu2I3 single crystal with high-PLQY (approximately 15.7%) intrinsic white-light emission via strongly localized 1D excitonic recombination. Adv. Mater. 31(46), e1905079 (2019)

[37] Yang, B., Yin, L., Niu, G., Yuan, J.H., Xue, K.H., Tan, Z., Miao, X.S., Niu, M., Du, X., Song, H., Lifshitz, E., Tang, J.: Lead-free halide Rb2CuBr 3 as sensitive X-ray scintillator. Adv. Mater. 31(44), e1904711 (2019)

[38] Sun, X.J., Xia, M.L., Xu, Y.S., Tang, J., Niu, G.D.: Research progress of perovskite direct X-ray imaging. Chinese J. Luminescence 43(7), 1014–1026 (2022)

[39] Xu, Y., Li, Y., Wang, Q., Chen, H., Lei, Y., Feng, X., Ci, Z., Jin, Z.: Two-dimensional BA₂PbBr 4-based wafer for X-rays imaging application. Mater. Chem. Front. 6(10), 1310–1316 (2022)

[40] Zeng, J., Bi, L., Cheng, Y., Xu, B., Jen, A.K.Y.: Self-assembled monolayer enabling improved buried interfaces in blade-coated perovskite solar cells for high efficiency and stability. Nano Res Energy 1, e9120004 (2022)

[41] Zhang, F., Zhao, Z., Chen, B., Zheng, H., Huang, L., Liu, Y., Wang, Y., Rogach, A.L.: Strongly emissive lead-free 0D Cs₃Cu₂I₅ perovskites synthesized by a room temperature solvent evaporation crystallization for down-conversion light-emitting devices and fluorescent inks. Adv. Opt. Mater. 8(8), 1901723(2020)

[42] Lin, R., Zhu, Q., Guo, Q., Zhu, Y., Zheng, W., Huang, F.: Dual self-trapped exciton emission with ultrahigh photoluminescence quantum yield in CsCu₂I₃ and Cs₃Cu₂I₅ perovskite single crystals. J. Phys. Chem. C 124(37), 20469–20476 (2020)

[43] Zhou, Z., Li, Y., Xing, Z., Sung, H.H.Y., Williams, I.D., Li, Z., Wong, K.S., Halpert, J.E.: Rapid synthesis of bright, shape-controlled, large single crystals of Cs₃Cu₂X₅ for phase pure single (X=Br, Cl) and mixed halides (X=Br, Cl) as the blue and green components for printable white light-emitting devices. Adv. Mater. Interfaces 8(20), 2101471(2021)

[44] Mo, X., Li, T., Huang, F., Li, Z., Zhou, Y., Lin, T., Ouyang, Y., Tao, X., Pan, C.: Highly-efficient all-inorganic lead-free 1D CsCu₂I₃ single crystal for white-light emitting diodes and UV photodetection. Nano Energy 81, 105570(2021)

[45] Zhao, X., Niu, G., Zhu, J., Yang, B., Yuan, J.H., Li, S., Gao, W., Hu, Q., Yin, L., Xue, K.H., Lifshitz, E., Miao, X., Tang, J.: Allinorganic copper halide as a stable and self-absorption-free X-ray scintillator. J. Phys. Chem. Lett. 11(5), 1873–1880 (2020)

[46] Yang, J., Kang, W., Liu, Z., Pi, M., Luo, L.B., Li, C., Lin, H., Luo, Z., Du, J., Zhou, M., Tang, X.: High-performance deep ultraviolet photodetector based on a one-dimensional lead-free halide perovskite CsCu₂I₃ film with high stability. J. Phys. Chem. Lett. 11(16), 6880–6886 (2020)

[47] Ma, Z., Shi, Z., Qin, C., Cui, M., Yang, D., Wang, X., Wang, L., Ji, X., Chen, X., Sun, J., Wu, D., Zhang, Y., Li, X.J., Zhang, L., Shan, C.: Stable yellow light-emitting devices based on ternary copper halides with broadband emissive self-trapped excitons. ACS Nano 14(4), 4475–4486 (2020)

[48] Roccanova, R., Yangui, A., Nhalil, H., Shi, H., Du, M.H., Saparov, B.: Near-unity photoluminescence quantum yield in blue-emitting Cs₃Cu₂Br_5?xI_x (0 ≤ x ≤ 5). ACS Appl. Electron. Mater. 1(3), 269–274 (2019)

[49] Xie, L., Chen, B., Zhang, F., Zhao, Z., Wang, X., Shi, L., Liu, Y., Huang, L., Liu, R., Zou, B., Wang, Y.: Highly luminescent and stable lead-free cesium copper halide perovskite powders for UVpumped phosphor-converted light-emitting diodes. Photon. Res. 8(6), 768–775 (2020)

[50] Cheng, P., Sun, L., Feng, L., Yang, S., Yang, Y., Zheng, D., Zhao, Y., Sang, Y., Zhang, R., Wei, D., Deng, W., Han, K.: Colloidal synthesis and optical properties of all-inorganic low-dimensional cesium copper halide nanocrystals. Angew. Chem. Int. Ed. 58(45), 16087–16091 (2019)

[51] Li, Y., Vashishtha, P., Zhou, Z., Li, Z., Shivarudraiah, S.B., Ma, C., Liu, J., Wong, K.S., Su, H., Halpert, J.E.: Room temperature synthesis of stable, printable Cs₃Cu₂X₅ (X = I, Br/I, Br, Br/Cl, Cl) colloidal nanocrystals with near-unity quantum yield green emitters (X = Cl). Chem. Mater. 32(13), 5515–5524 (2020)

[52] Luo, Z., Li, Q., Zhang, L., Wu, X., Tan, L., Zou, C., Liu, Y., Quan, Z.: 0D Cs₃Cu₂X₅ (X = I, Br, and Cl) nanocrystals: colloidal syntheses and optical properties. Small 16(3), e1905226 (2020)

[53] Zhao, S., Chen, C., Cai, W., Li, R., Li, H., Jiang, S., Liu, M., Zang, Z.: Efficiently luminescent and stable lead-free Cs₃Cu₂Cl₅@ silica nanocrystals for white light-emitting diodes and communication. Adv. Opt. Mater. 9(13), 2100307(2021)

[54] Zhang, R., Mao, X., Zheng, D., Yang, Y., Yang, S., Han, K.: A lead-free all-inorganic metal halide with near-unity green luminescence. Laser Photonics Rev. 14(5), 2000027(2020)

[55] Han, L., Sun, B., Guo, C., Peng, G., Chen, H., Yang, Z., Li, N., Ci, Z., Jin, Z.: Photophysics in zero-dimensional potassium-doped cesium copper chloride Cs3Cu2Cl5 nanosheets and its application for high-performance flexible X-ray detection. Adv. Opt. Mater. 10(6), 2102453(2022)

[56] Zhang, B., Wu, X., Zhou, S., Liang, G., Hu, Q.: Self-trapped exciton emission in inorganic copper(I) metal halides. Front Optoelectron. 14(4), 459–472 (2021)

[57] Du, M.H.: Emission trend of multiple self-trapped excitons in luminescent 1D copper halides. ACS Energy Lett. 5(2), 464–469 (2020)

[58] Zhang, Z.X., Li, C., Lu, Y., Tong, X.W., Liang, F.X., Zhao, X.Y., Wu, D., Xie, C., Luo, L.B.: Sensitive deep ultraviolet photodetector and image sensor composed of inorganic lead-free Cs₃Cu₂I₅ perovskite with wide bandgap. J. Phys. Chem. Lett. 10(18), 5343–5350 (2019)

[59] Li, Y., Shi, Z., Liang, W., Wang, L., Li, S., Zhang, F., Ma, Z., Wang, Y., Tian, Y., Wu, D., Li, X., Zhang, Y., Shan, C., Fang, X.: Highly stable and spectrum-selective ultraviolet photodetectors based on lead-free copper-based perovskites. Mater. Horiz. 7(2), 530–540 (2020)

[60] Li, Y., Shi, Z., Wang, L., Chen, Y., Liang, W., Wu, D., Li, X., Zhang, Y., Shan, C., Fang, X.: Solution-processed one-dimensional CsCu₂I₃ nanowires for polarization-sensitive and flexible ultraviolet photodetectors. Mater. Horiz. 7(6), 1613–1622 (2020)

[61] Ma, J., Xia, X., Yan, S., Li, Y., Liang, W., Yan, J., Chen, X., Wu, D., Li, X., Shi, Z.: Stable and self-powered solar-blind ultraviolet photodetectors based on a Cs₃Cu₂I₅/β-Ga₂O₃ heterojunction prepared by dual-source vapor codeposition. ACS Appl. Mater. Interfaces 13(13), 15409–15419 (2021)

[62] Zhang, M., Zhu, J., Yang, B., Niu, G., Wu, H., Zhao, X., Yin, L., Jin, T., Liang, X., Tang, J.: Oriented-structured CsCu₂I₃ film by close-space sublimation and nanoscale seed screening for high-resolution X-ray imaging. Nano Lett. 21(3), 1392–1399 (2021)

[63] Zhao, X., Jin, T., Gao, W., Niu, G., Zhu, J., Song, B., Luo, J., Pan, W., Wu, H., Zhang, M., He, X., Fu, L., Li, Z., Zhao, H., Tang, J.: Embedding Cs₃Cu₂I₅ scintillators into anodic aluminum oxide matrix for high-resolution X-ray imaging. Adv. Opt. Mater. 9(24), 2101194(2021)

[64] Zhou, Q., Ren, J., Xiao, J., Lei, L., Liao, F., Di, H., Wang, C., Yang, L., Chen, Q., Yang, X., Zhao, Y., Han, X.: Highly efficient copper halide scintillators for high-performance and dynamic X-ray imaging. Nanoscale 13(47), 19894–19902 (2021)

[65] Wang, Q., Zhou, Q., Nikl, M., Xiao, J., Kucerkova, R., Beitlerova, A., Babin, V., Prusa, P., Linhart, V., Wang, J., Wen, X., Niu, G., Tang, J., Ren, G., Wu, Y.: Highly resolved X-ray imaging enabled by In(I) doped perovskite-like Cs₃Cu₂Cl₅ single crystal scintillator. Adv. Opt. Mater. 10(11), 2200304(2022)

[66] Li, X., Chen, J., Yang, D., Chen, X., Geng, D., Jiang, L., Wu, Y., Meng, C., Zeng, H.: Mn²⁺ induced significant improvement and robust stability of radioluminescence in Cs₃Cu₂I₅ for high-performance nuclear battery. Nat. Commun. 12(1), 3879(2021)

[67] Cheng, S., Nikl, M., Beitlerova, A., Kucerkova, R., Du, X., Niu, G., Jia, Y., Tang, J., Ren, G., Wu, Y.: Ultrabright and highly efficient all-inorganic zero-dimensional perovskite scintillators. Adv. Opt. Mater. 9(13), 2100460(2021)

RIGHTS & PERMISSIONS

2022 The Author(s) 2022