A Comparative Study of Four Total Variational Regularization Reconstruction Algorithms for Sparse-View Photoacoustic Imaging

Table 1

Numerical results of the NCAT phantom with different sampling positions.

Experimental	CPU time (seconds)				SNR (dB)				NMAE
Positions	TVAL3	L1magic	NESTA	TwIST	TVAL3	L1magic	NESTA	TwIST	TVAL3	L1magic	NESTA	TwIST
5	2.69	202.84	1.64	1.06	10.01	2.28	6.5	0.76	26.78	65.22	40.13	77.69
10	4.29	298.61	2.23	1.98	15.03	6.3	8.42	0.72	15.02	41.05	32.16	78.04
15	6.16	442.24	2.81	1.8	20.35	8.44	9.85	1.4	8.14	32.09	27.29	72.14
20	8.11	470.25	3.47	3.01	25.28	10.56	10.95	1.93	4.62	25.14	24.04	67.93
25	9.02	425.4	3.89	4.05	29.79	12.27	11.83	2.74	2.75	20.65	21.71	61.83
30	9.93	474.47	4.33	3.13	33.3	13.3	12.31	4.19	1.83	18.34	20.56	52.33
35	10.59	446.89	4.78	3.71	37.09	14.55	13.11	6.15	1.18	15.87	18.74	41.77
40	10.93	444.98	5.06	4.78	40.83	15.62	13.86	9.3	0.77	14.04	17.18	29.06
45	10.94	532.55	5.44	2.2	42.53	17.29	14.33	7.16	0.63	11.58	16.28	37.17
50	10.99	509.73	5.89	3.92	45.71	18.27	15.06	14.36	0.44	10.35	14.98	16.23
55	12.16	474.01	6.27	3.61	48.34	19.4	15.7	18.33	0.32	9.09	13.91	10.28
60	11.8	559.14	6.81	4.33	51.05	20.87	16.3	27.72	0.24	7.67	12.98	3.49
Average	8.97	440.09	4.39	3.13	33.28	13.26	12.35	7.9	5.23	22.59	21.66	45.66