--- a/img/brainphantom256x256_pdps_known_proximal_15_10000_9000.txt Wed Apr 24 16:55:04 2024 +0300 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,78 +0,0 @@ -# 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.15, ρ̃₀ = 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, phantom_ρ = 100, save_prefix = "img/brainphantom256x256_pdps_known_proximal_15_10000_9000", dynrange = 1.0, Λ = 1.0) -iter time function_value psnr ssim -1 7.208000170066953e-6 561764.4291806474 -0.02636674332749045 0.0041734697881560985 -2 0.029386874986812472 477772.60352856905 0.47416210950828774 0.007415531622632549 -3 0.04420633299741894 445831.54201752564 0.8556395922189581 0.010119589144049099 -4 0.07331924999016337 399652.6986019082 1.1834220640538975 0.0121473521051133 -5 0.08836183298262767 362400.4389922261 1.494333619204808 0.014149753113820093 -6 0.11764654098078609 323024.03425067704 1.778254249533942 0.01522774635184404 -7 0.13247116599814035 283882.8640395682 2.0858838868164566 0.016114624787739226 -8 0.1616962069820147 245919.64670498838 2.4005183696836827 0.017897547454711907 -9 0.17658429098082706 205195.3352253661 2.7276762982458527 0.020041397138472895 -10 0.2055951239890419 166431.96416691126 3.0500610179333396 0.021669226677280964 -11 0.22041212397743948 129350.44266024753 3.387480496295 0.023541047843757906 -12 0.24945612397277728 92339.67907049721 3.7208727047749894 0.02638172254501852 -13 0.26427566597703844 56147.77029549168 4.064487757284027 0.02929538506547372 -14 0.29331962496507913 20572.118110416806 4.410528453062292 0.03540857231991382 -15 0.3081306659732945 -10717.21148112784 4.7510965050743 0.03626512680470327 -16 0.3384145409800112 -42786.43826475472 5.10495633423716 0.04426236117452481 -17 0.3533202069811523 -63384.87202863484 5.4230117332797025 0.03893011155319675 -18 0.38246474898187444 -105837.90759620132 5.8205728962821714 0.07071721320744478 -19 0.39733216498279944 -136277.43848340202 6.18796899455673 0.09007742956462768 -20 0.4263431649596896 -165454.1526176333 6.556268797608462 0.1167238847931111 -30 0.645438123960048 -350596.48582448537 10.261929469722748 0.4515319317673516 -40 0.8636804159614258 -421404.5116904354 13.762296794871032 0.5367584291902869 -50 1.0820691239496227 -448518.3978258499 16.564121685319297 0.5614030823582099 -60 1.3005964159383439 -457373.9085048995 18.30588057221346 0.5914100338476127 -70 1.5189499579137191 -459432.86409579293 19.19921858592214 0.609679583569413 -80 1.7374591249099467 -460585.59002685547 19.66333186496991 0.6220467099606226 -90 1.9560254999087192 -460996.87652548985 19.947649623108553 0.6266827743602128 -100 2.173608874902129 -460994.3826579701 20.140385260812828 0.6290103992314371 -110 2.378357624897035 -460837.70738044876 20.25010008894009 0.6300870357362052 -120 2.5969945419055875 -461412.0924345887 20.332815310297896 0.6285638546521161 -130 2.8144480418995954 -461658.89207825245 20.272447037818814 0.6182707705242505 -140 3.03174116689479 -459697.5887098162 20.010457281062187 0.5991754925400398 -150 3.248936666874215 -461388.0745708661 20.432047592212626 0.6249815949447837 -160 3.4659844998677727 -459598.46030530456 20.086587765554206 0.6060105063583754 -170 3.682812124839984 -461263.03743456514 20.55083790223796 0.6314303911125991 -180 3.8998824578593485 -460785.1983080981 20.428615008101904 0.6190426413310206 -190 4.116731790883932 -460386.24951323256 20.425940217084047 0.6199531030423062 -200 4.333990249899216 -460332.37921665737 20.328014325985503 0.6132154651987625 -300 6.510508208884858 -461470.7027605368 20.63198556687535 0.6297499279720288 -400 8.689072999870405 -461068.78618832177 20.676103651816092 0.6304578785624024 -500 10.872397541883402 -460889.52990054723 20.65336294609731 0.6313994099699308 -600 13.04934995790245 -461705.69892908004 20.7511004631057 0.6399725002028397 -700 15.216933290910674 -461165.2126459676 20.7517820270491 0.6380474790019816 -800 17.381680123886326 -461245.2000950096 20.714654637638752 0.638138517691965 -900 19.5476342489128 -461349.7385930514 20.79399066338013 0.6406764080286608 -1000 21.74007312388858 -461368.10407306807 20.605447117399546 0.6289342386097674 -1100 23.924216457875445 -461527.83026946115 21.616291773315723 0.6817806732254383 -1200 26.09317041587201 -461248.647813059 21.98721232568007 0.6993059160637243 -1300 28.260992207884556 -461636.5369824387 22.197349355010175 0.7079662332492418 -1400 30.429287915903842 -461471.7551362793 22.326098831792347 0.714994318112899 -1500 32.5962968749227 -461477.0226392162 22.405944486933382 0.7185020441894621 -1600 34.76356783290976 -461508.7001345785 22.46188308603088 0.7223598761565722 -1700 36.94523666691384 -461416.29635482014 22.503393805210333 0.7246931297489942 -1800 39.1128853749251 -461526.79680164054 22.535319903902668 0.7259730883269055 -1900 41.295438707922585 -461530.28649795457 22.555507666566527 0.7254438064507713 -2000 43.46189366694307 -461609.84374154016 22.57423108919398 0.7276674269361936 -2100 45.65475829193019 -460790.5980782592 20.882845551605985 0.6342910307319234 -2200 47.82176862593042 -461616.32627589116 20.604385686471993 0.6197987399008086 -2300 49.99434533395106 -460848.1997743013 20.67723909016766 0.6324957287172485 -2400 52.1647558749537 -461337.6764989341 20.728983824634383 0.6374979506212082 -2500 54.331151999969734 -460707.1907771798 20.55463800214148 0.6241803058581371 -2600 56.500439666968305 -461541.28551366477 20.754409053805468 0.6406513388218529 -2700 58.66917033394566 -460935.4430990912 20.708950402739337 0.6366894086388818 -2800 60.88020183396293 -461599.3130413398 20.755683434617033 0.638431925728777 -2900 63.04865087597864 -462060.66440661007 20.732091954982465 0.6376257718833535 -3000 65.21420025097905 -461931.5317004886 20.73139839092415 0.6351303235745192 -3100 67.39173433399992 -459670.6183546943 20.373206920737836 0.6084590027724093 -3200 69.56313725101063 -459507.5774207035 19.87401427579497 0.5815126734639297 -3300 71.73393629200291 -460772.1549017727 20.401709350018148 0.6095554652807589 -3400 73.92400345901842 -460797.7804659181 20.480097357608685 0.614274134915554 -3500 76.1067618750385 -460861.64256642124 20.597212522349935 0.6278432623094587 -3600 78.2815276660258 -461978.03524543456 21.585960742251714 0.6802457280868743 -3700 80.45292812405387 -462026.5804320709 21.973028053232863 0.6966235300338117 -3800 82.62317233305657 -461970.6727958898 22.188329138762313 0.7076371472070837 -3900 84.7941199990455 -462094.1608076924 22.318997047023824 0.7137776150322142 -4000 86.97895983204944 -461888.7879098893 22.40499723244197 0.7182709936334599