img/brainphantom256x256_pdps_known_zerodual_20_10000_9000.txt

changeset 31
cbedcfcc0557
parent 30
dec72b794b08
child 32
88632284396f
equal deleted inserted replaced
30:dec72b794b08 31:cbedcfcc0557
1 # params = (ρ = 0, verbose_iter = 100, maxiter = 4000, save_results = true, save_images = true, save_images_iters = Set([500, 1800, 2300, 3000, 1300, 3800, 1000, 2800, 300, 4000, 3500, 2000, 2500, 3300, 800, 100, 1500]), pixelwise_displacement = false, dual_flow = true, prox_predict = true, handle_interrupt = true, init = :zero, plot_movement = false, stable_interval = Set([1144, 3634, 1175, 1953, 1546, 1703, 1956, 3601, 1028, 1299, 1438, 1074, 3672, 3705, 3631, 1812, 1199, 1090, 1437, 1662, 1370, 1991, 3930, 3608, 1704, 1663, 1957, 1266, 1646, 1823, 1015, 1373, 3981, 3719, 3646, 3982, 3655, 1168, 3855, 1039, 1085, 3914, 3794, 1050, 1276, 1285, 3540, 3524, 1356, 3762, 1466, 3711, 1036, 3549, 1171, 1739, 1625, 1669, 1058, 1814, 3963, 1102, 1388, 1905, 3765, 3943, 1570, 3651, 1745, 1574, 1864, 1117, 1859, 1338, 1636, 3915, 1861, 3754, 3570, 1896, 3510, 1075, 1342, 3557, 1725, 1873, 1526, 1769, 1185, 1928, 3937, 1922, 1939, 1504, 1052, 1374, 1807, 1724, 3587, 1245, 3791, 3652, 3628, 1138, 1095, 1604, 1848, 1155, 1111, 1916, 1906, 1251, 3938, 1341, 1753, 1278, 1221, 1621, 3731, 1497, 1819, 1845, 3728, 3809, 3677, 3617, 1521, 3713, 1539, 1670, 3803, 1672, 1808, 1661, 3729, 1056, 1610, 1236, 3504, 3571, 3743, 3722, 3507, 1070, 1994, 1548, 1357, 1903, 1580, 1678, 1004, 1091, 1147, 1045, 1255, 1635, 3526, 1001, 3786, 1972, 1274, 1336, 3940, 3792, 1157, 1836, 1073, 1612, 1282, 3835, 1448, 3856, 3757, 3740, 3648, 3815, 1578, 3691, 1304, 2000, 1372, 1849, 3877, 3979, 1475, 1585, 1783, 3859, 1378, 3602, 1719, 3614, 1259, 1528, 3956, 1325, 1287, 1935, 3574, 1019, 1913, 3770, 1354, 1238, 1559, 1249, 1818, 1021, 1519, 3844, 1201, 1313, 3851, 1362, 1326, 3696, 1449, 1054, 3825, 1120, 1908, 3657, 1923, 1302, 1637, 3690, 1232, 1395, 1344, 3965, 1698, 3530, 1208, 1811, 1333, 1941, 3817, 1458, 1865, 3886, 1116, 3931, 3926, 1268, 1732, 3568, 3828, 1714, 1494, 1846, 3702, 3923, 1642, 1398, 3616, 1883, 1079, 1506, 1186, 1582, 3773, 3806, 1043, 1535, 3858, 1872, 1103, 1319, 3972, 1246, 1152, 1409, 1560, 1996, 3903, 1145, 1767, 1680, 1792, 1141, 1009, 1566, 1137, 1360, 1565, 1854, 1588, 1538, 3519, 1315, 1681, 1723, 3681, 3917, 1283, 3509, 1825, 3868, 3870, 3819, 1770, 3593, 1063, 1235, 1614, 1166, 1064, 3641, 1355, 1151, 1169, 3700, 3973, 3986, 4000, 3784, 1452, 1479, 3962, 3508, 1576, 1830, 3644, 1183, 3599, 1842, 1377, 3566, 1046, 1564, 1493, 1180, 3550, 1007, 3837, 1256, 1993, 3850, 1759, 1998, 3922, 3697, 1909, 3609, 1040, 1937, 1728, 3525, 3990, 3759, 1384, 1562, 3658, 3925, 3971, 3649, 1675, 3928, 1600, 3693, 1543, 1172, 1135, 1607, 3737, 3977, 3852, 1966, 1471, 1618, 3680, 1927, 1309, 1529, 1797, 3764, 3909, 1664, 1999, 3801, 1476, 1403, 1041, 1839, 1813, 1595, 1900, 3620, 1442, 1296, 3751, 3505, 1123, 3832, 1743, 1942, 1228, 1381, 3501, 1899, 1901, 3703, 1727, 1177, 1850, 3613, 1393, 1886, 1434, 1469, 3880, 1591, 1673, 1027, 1581, 1002, 1113, 3882, 3686, 1391, 1260, 1109, 1383, 3622, 3576, 1429, 1435, 3556, 1634, 3861, 1676, 1012, 1987, 1730, 1425, 1347, 1620, 1831, 3821, 3905, 1198, 3952, 3625, 1140, 1404, 3955, 1720, 1414, 3901, 1875, 3684, 3933, 1544, 3692, 1394, 1761, 1006, 1748, 3541, 3788, 1207, 3862, 1523, 1067, 3579, 3727, 1284, 1044, 1066, 1488, 1520, 1961, 1016, 1617, 1763, 1706, 1159, 3590, 1337, 1781, 1771, 1766, 3607, 1912, 1483, 3747, 3920, 1252, 1167, 1324, 1793, 1294, 1127, 1192, 1352, 3753, 1944, 3775, 1988, 3744, 1156, 1895, 1921, 1917, 1170, 1631, 3582, 1068, 1701, 1542, 3793, 3569, 1705, 1345, 1597, 1321, 1290, 1666, 1189, 3538, 3633, 1623, 1505, 1744, 1622, 1047, 1834, 3769, 3534, 3707, 1779, 1218, 3584, 3807, 1272, 1958, 1639, 1982, 1132, 1816, 1667, 1887, 1379, 3567, 3585, 1173, 1051, 1729, 1303, 1507, 1650, 1205, 1862, 1871, 3954, 1463, 1656, 1777, 1918, 3547, 3627, 1240, 1334, 1653, 3708, 1129, 3823, 1416, 1482, 1624, 1413, 3799, 1933, 1665, 1086, 1496, 1555, 1243, 3748, 3522, 1866, 3863, 3883, 3983, 1525, 1204, 1380, 1225, 1919, 1550, 1035, 1470, 1184, 1649, 1216, 3725, 1194, 1693, 3517, 3621, 1963, 1048, 1514, 1660, 1440, 1231, 1545, 1286, 3798, 1758, 1125, 1584, 1700, 1852, 1332, 3826, 1220, 1415, 1734, 1615, 3812, 1951, 3527, 3985, 1776, 1267, 1472, 1013, 3756, 1323, 3957, 1754, 3726, 1485, 1930, 1212, 1645, 1997, 1134, 1847, 1322, 3787, 1311, 1551, 1078, 1081, 1101, 1490, 1110, 1764, 3660, 1003, 1697, 1230, 3736, 1197, 1643, 1911, 1298, 1527, 3811, 1242, 1420, 3839, 1696, 1474, 3561, 1806, 3797, 1568, 1269, 1348, 3624, 3539, 1209, 1153, 3642, 3988, 3704, 1626, 3896, 1737, 1178, 1486, 3712, 1423, 1353, 1890, 3518, 3885, 3996, 3721, 1014, 1694, 1860, 3818, 1952, 1387, 1879, 1124, 3829, 3664, 3595, 3610, 1695, 3667, 1461, 1804, 1567, 3827, 1389, 1587, 3881, 1328, 1182, 1349, 1386, 1426, 1445, 3654, 1306, 1112, 1782, 1084, 1513, 1467, 1736, 1133, 1262, 1931, 1340, 3942, 1794, 1060, 3558, 1264, 1271, 1596, 3783, 1962, 1162, 3560, 3929, 1833, 1508, 3647, 1316, 3813, 1291, 1206, 1613, 3531, 1837, 1281, 1740, 3840, 1799, 3932, 1417, 1428, 1258, 1254, 1358, 1707, 1577, 3750, 3968, 1072, 1686, 3678, 3573, 1108, 3589, 1644, 1750, 1838, 3565, 3854, 1973, 1853, 3715, 3594, 3842, 1751, 1950, 3669, 1011, 3998, 1008, 1396, 1093, 1983, 1148, 1450, 1314, 1176, 1583, 3867, 1757, 3846, 3626, 1843, 1010, 1422, 3830, 3710, 1503, 3778, 1457, 1722, 3910, 1965, 3911, 1468, 1616, 1191, 1293, 1351, 1094, 3913, 1515, 3665, 1844, 1401, 1193, 1037, 1092, 1932, 1877, 3761, 3875, 1477, 3884, 1366, 1335, 1297, 3777, 3969, 1702, 1361, 1088, 1473, 1594, 3895, 3572, 3802, 1510, 1163, 1301, 1142, 1265, 1684, 1975, 1022, 1224, 1139, 1436, 3772, 1752, 3833, 3675, 1233, 1214, 1237, 1465, 1652, 1367, 3987, 1690, 1893, 3597, 1431, 1130, 1055, 1870, 3745, 1000, 1305, 1970, 3831, 1257, 1629, 1460, 1196, 1363, 1876, 1327, 3536, 1190, 1885, 1762, 3716, 1924, 1029, 1627, 1891, 3662, 1076, 1034, 1061, 1462, 1599, 3637, 1682, 1119, 1765, 1534, 1026, 1312, 1948, 1943, 1598, 1146, 1300, 1484, 1042, 1784, 1331, 1097, 1310, 1824, 1203, 1499, 1685, 1810, 1082, 3876, 3580, 3653, 1721, 3537, 1914, 1747, 1657, 1986, 1715, 3709, 1981, 3992, 1934, 1980, 1679, 1571, 3927, 1089, 3623, 3533, 3650, 1279, 1018, 3553, 3741, 3984, 1775, 1925, 1495, 3735, 1023, 3600, 1211, 3611, 1822, 1524, 1829, 3892, 1926, 3699, 1418, 1031, 1080, 1840, 1263, 3908, 1270, 1756, 3659, 1248, 3562, 1350, 3544, 1537, 3874, 3894, 3897, 3532, 1478, 1215, 1411, 1553, 1638, 1038, 1364, 3959, 1406, 3635, 1541, 1976, 1441, 1407, 1710, 3845, 1516, 3780, 3666, 3824, 1057, 1857, 3720, 1408, 3888, 1330, 1894, 1114, 3763, 1531, 1969, 3934, 3906, 1785, 1558, 1726, 1960, 3790, 1532, 3695, 3995, 1795, 1716, 1668, 3596, 3732, 3618, 1907, 3563, 1502, 1511, 1633, 3918, 1791, 3630, 3766, 1371, 3948, 1802, 1946, 1187, 1421, 1552, 3723, 1889, 1755, 1397, 3640, 1005, 1995, 1165, 1659, 3834, 1400, 3890, 1971, 1164, 1104, 1234, 1881, 3991, 1244, 1647, 3871, 1222, 3964, 3506, 3774, 3916, 1121, 1498, 1878, 1547, 1025, 1640, 1424, 3742, 1863, 3724, 3836, 3843, 1540, 3632, 1405, 3670, 3980, 1343, 1059, 1671, 3873, 1049, 1735, 3767, 1630, 3779, 1821, 1464, 3588, 3978, 1786, 1455, 1851, 1979, 3545, 1856, 3760, 3717, 1832, 3639, 3535, 1805, 3528, 3872, 1892, 3781, 3730, 3543, 1365, 1053, 3796, 3899, 1788, 3752, 1920, 3866, 1122, 3936, 1945, 1586, 1984, 1898, 3701, 3604, 3838, 3921, 3583, 1512, 1107, 1247, 1787, 1033, 1855, 1651, 3638, 1087, 1882, 1188, 3941, 1077, 3513, 3999, 1692, 3820, 1619, 1359, 1742, 1253, 1609, 1593, 1590, 3970, 3516, 1131, 1569, 1501, 1017, 1772, 1874, 3912, 1654, 3864, 1223, 1241, 3661, 3668, 1150, 1768, 3961, 1250, 1632, 1280, 1884, 1712, 1575, 1105, 3521, 1433, 1202, 1481, 3656, 3689, 1938, 1954, 3860, 1213, 1641, 3902, 1557, 1158, 3629, 1796, 3552, 1320, 3945, 3575, 1688, 1803, 3958, 1790, 3810, 1154, 1292, 3603, 3785, 3967, 1992, 1800, 3694, 1382, 1308, 3949, 3994, 1275, 1674, 3869, 1940, 1369, 3673, 1691, 1518, 3889, 1071, 1897, 1985, 1603, 1967, 1024, 3887, 1867, 3739, 1446, 1841, 1427, 3768, 3577, 1317, 3808, 3619, 1974, 1798, 3891, 1815, 3546, 1959, 1746, 1443, 3822, 1915, 3733, 1128, 1454, 3947, 1329, 1227, 3879, 1554, 3679, 1828, 1376, 3795, 3520, 3953, 3612, 1717, 3848, 1563, 1536, 1447, 3800, 1868, 1160, 1989, 3771, 3898, 1990, 1556, 1683, 1880, 3893, 1273, 1605, 1749, 1492, 1936, 1738, 1573, 1606, 1219, 1709, 1589, 1713, 1339, 3636, 1239, 1677, 1459, 1731, 1968, 1977, 1118, 3924, 1020, 1910, 1390, 1100, 3904, 3758, 1106, 1419, 1307, 3682, 1161, 1385, 3849, 3976, 1820, 3688, 3746, 1126, 1530, 3605, 3551, 1500, 1444, 1689, 3559, 1195, 1229, 3663, 1522, 1432, 1453, 1099, 3789, 3542, 1602, 3975, 3841, 1655, 1083, 3512, 1947, 1592, 3878, 1733, 1835, 3564, 3685, 3997, 1489, 3529, 1136, 1062, 3907, 1289, 1069, 3615, 1611, 1904, 1741, 3503, 1773, 1929, 3586, 1978, 1412, 1699, 3865, 3714, 1658, 3554, 1098, 1517, 1760, 1096, 1949, 1708, 3950, 1491, 3951, 1277, 1858, 1509, 1456, 1030, 1346, 1179, 1572, 1817, 3993, 1780, 1210, 3674, 1648, 1801, 1964, 3782, 3814, 1392, 3804, 1809, 3944, 3581, 3749, 3776, 3598, 3671, 3935, 1826, 3755, 3989, 1955, 1410, 3939, 3853, 1439, 1628, 1288, 3974, 3738, 1149, 3645, 1200, 1181, 3523, 3511, 3718, 1778, 1295, 3606, 3578, 1174, 1789, 3946, 3687, 1601, 1430, 1718, 1902, 3919, 1487, 1261, 3500, 1888, 3816, 1399, 3706, 3548, 3676, 1533, 1869, 3900, 1368, 3643, 1032, 1318, 1579, 3592, 1226, 3805, 3502, 1143, 3960, 3857, 1065, 1451, 1774, 1375, 1549, 1561, 3698, 1480, 1687, 3515, 3966, 3734, 1608, 1217, 3683, 3555, 3591, 3847, 1711, 3514, 1827, 1402, 1115]), noise_level = 0.5, shake_noise_level = 0.1, shake = 1.0, rotation_factor = 0.075, rotation_noise_level = 0.0075, α = 0.2, ρ̃₀ = 1.0, σ̃₀ = 1.0, δ = 0.9, σ₀ = 1.0, τ₀ = 0.9, λ = 1, radondims = [128, 64], sz = (256, 256), scale = 1, c = 1.0, sino_sparsity = 0.5, L = 300.0, L_experiment = false, save_prefix = "img/brainphantom256x256_pdps_known_zerodual_20_10000_9000", dynrange = 1.0, Λ = 1.0)
2 iter time function_value psnr ssim
3 1 8.749979315325618e-6 563529.9047468598 -0.02636674332749045 0.0041734697881560985
4 2 0.027876249951077625 478856.08122394385 0.4727406757380256 0.007392494550781689
5 3 0.0412792079441715 446629.7517429079 0.8539140607374829 0.010069665864595012
6 4 0.06899879095726646 400269.21028804214 1.1814775051209048 0.01206703098275115
7 5 0.08268387493444607 362962.4296867711 1.4922690080185277 0.014042451163532112
8 6 0.11036995792528614 323591.6316846677 1.7754944575834577 0.01508482569419906
9 7 0.12366241592098959 284365.7473419522 2.0830435884394536 0.015941529087566995
10 8 0.15149541592109017 246354.00136495073 2.397592669599506 0.01769874438581332
11 9 0.1650198329007253 205571.0802674675 2.7246366424854465 0.019819222915664732
12 10 0.19340866690617986 166785.36886898818 3.046645922327711 0.021419885750702854
13 11 0.20701108290813863 129674.12709935356 3.383895663760555 0.023311111836831093
14 12 0.23496912492555566 92653.28690360574 3.7170275054469064 0.026093821873056477
15 13 0.24838299994007684 56452.42106017006 4.060330686248957 0.02889385552418672
16 14 0.27641983394278213 20857.101703946475 4.406071274468021 0.034948265482444334
17 15 0.28997666694340296 -10438.770138025126 4.7464343299667675 0.03570432755419898
18 16 0.31794804194942117 -42505.71385518951 5.09999542278511 0.04368141902332796
19 17 0.33136675096466206 -63126.3486904598 5.417650255125764 0.03836780958154917
20 18 0.35913070995593444 -105540.98876235927 5.815058755545951 0.06961002429483774
21 19 0.37249391796649434 -135974.1889506175 6.182079030583943 0.08845381424367746
22 20 0.4002704179729335 -165136.5459505766 6.549971469614685 0.11495759048821691
23 30 0.6057929589878768 -350415.0314217395 10.252427224594552 0.4471126464336167
24 40 0.8108901259838603 -421351.0510403995 13.745084438160491 0.5359227561921107
25 50 1.0156440009886865 -448479.0073656686 16.53690856952951 0.5606411527744863
26 60 1.2197823339665774 -457336.8096845315 18.27931489930829 0.5904973995238473
27 70 1.4234844589664135 -459398.0840489778 19.179624185637728 0.6093611914294952
28 80 1.627001124958042 -460559.5591591289 19.6561073243845 0.6226485701401986
29 90 1.8303712499618996 -460974.94356359495 19.952200148056296 0.6280188369990043
30 100 2.033651915960945 -460976.9330205219 20.158006520071332 0.6310120146415426
31 110 2.223528998962138 -460823.35601861565 20.28045260848213 0.6323698991507318
32 120 2.4310180399625096 -461399.4134079207 20.371377467926813 0.6309711207896302
33 130 2.6399422899703495 -461647.8006028427 20.30448183839954 0.6193954448838336
34 140 2.8482378729677293 -459669.0638747321 20.033771146525563 0.5999480920533595
35 150 3.056630247971043 -461382.3516274438 20.488673631379672 0.6275619497138916
36 160 3.261942831974011 -459584.76990198955 20.117378070081273 0.6064410771422082
37 170 3.4653351229499094 -461262.8949247521 20.632698623775173 0.6349720698087877
38 180 3.6690817899652757 -460776.93593953113 20.508040251631897 0.6220042768785026
39 190 3.873239622975234 -460380.741208267 20.506010644425544 0.6229426538602285
40 200 4.0777439149678685 -460321.5768809585 20.405097040206886 0.6158879739860437
41 300 6.1171862069750205 -461481.0754801085 20.74744573932619 0.6339098676634498
42 400 8.15905512298923 -461076.03712357854 20.80750624836734 0.6350963001174728
43 500 10.19770728997537 -460897.45434653264 20.785116329988476 0.6361090967152319
44 600 12.221155372972135 -461715.21774791036 20.8885812245061 0.6447264460010231
45 700 14.252476414956618 -461174.46579402883 20.889219459912386 0.6433037109515745
46 800 16.28088845696766 -461252.9004520723 20.843357967092256 0.6427375585920718
47 900 18.369778872962343 -461358.4785052773 20.93078477800351 0.6456148736715817
48 1000 20.434952913958114 -461372.8941923818 20.718958180955326 0.6327190681059426
49 1100 22.510358997969888 -461414.9701224553 21.766102959638015 0.673821653762958
50 1200 24.568871122959536 -461041.9166435804 22.003929647010448 0.6807280440248624
51 1300 26.620753414958017 -461372.63751805713 22.090299527481303 0.6834402398892566
52 1400 28.689906873943983 -461170.01896055706 22.128976681629965 0.68542545689322
53 1500 30.758403623942286 -461160.2126787338 22.14814323034254 0.686039807025556
54 1600 32.84527916595107 -461181.97184875305 22.158852066450976 0.6873869724378071
55 1700 34.89148724995903 -461081.1592314979 22.166042041787154 0.68789123111372
56 1800 36.95888554197154 -461185.26824670244 22.17030504773458 0.6875923962828542
57 1900 39.00933399997302 -461184.5751019353 22.173225250094085 0.6859667063517108
58 2000 41.08802108396776 -461258.5980084426 22.17592448055875 0.6878784679003866
59 2100 43.132536291959696 -460771.53124172956 20.927094120595704 0.6373595114078217
60 2200 45.182526332937414 -461620.79945889197 20.701856244574948 0.6236782213101687
61 2300 47.241347124945605 -460856.5242060274 20.79806687347351 0.6369615329273195
62 2400 49.273692499962635 -461346.24083914765 20.86233673773571 0.6422327739324779
63 2500 51.30414429094526 -460707.61172871926 20.662297654560877 0.6274252244047018
64 2600 53.32980708294781 -461549.6827994127 20.89159447838363 0.64551323986934
65 2700 55.377858207939425 -460944.5930327483 20.840278730844854 0.6412001573412833
66 2800 57.3990448739205 -461607.95922741835 20.89508058479119 0.6435383338540954
67 2900 59.44118787391926 -462069.6194101821 20.868371095578603 0.6423312207735403
68 3000 61.509484915906796 -461939.7728975037 20.863462578975806 0.6400346980387652
69 3100 63.57368387392489 -459662.8110055467 20.473002640122978 0.6119351687778724
70 3200 65.65636933292262 -459492.1410088092 19.92206478983645 0.5818023435679106
71 3300 67.69076954090269 -460777.87024429097 20.512673853422765 0.6128868277123647
72 3400 69.73519895688514 -460807.84031042387 20.59940201450136 0.6182434157116807
73 3500 71.76092108190642 -460866.9791645211 20.721798843383645 0.6319543021296905
74 3600 73.82680079090642 -461865.55813633965 21.721787076828377 0.6718406246758974
75 3700 75.96071062490228 -461811.69497260876 21.96519136008503 0.6772952205049021
76 3800 78.05097733391449 -461697.5668382382 22.052473958616435 0.6810839610534091
77 3900 80.13434595891158 -461790.21643061546 22.09139141098718 0.6825147958910878
78 4000 82.19546791689936 -461565.1222510449 22.11075340450595 0.6839747971613462

mercurial