Knowledge

Horn–Schunck method

Source 📝

713: 579: 1886: 1689: 585: 451: 1027: 1487:
and may be solved for each pixel in the image. However, since the solution depends on the neighboring values of the flow field, it must be repeated once the neighbors have been updated. The following iterative scheme is derived using
1442: 1315: 1695: 1498: 940: 839: 1124: 227: 708:{\displaystyle {\frac {\partial L}{\partial v}}-{\frac {\partial }{\partial x}}{\frac {\partial L}{\partial v_{x}}}-{\frac {\partial }{\partial y}}{\frac {\partial L}{\partial v_{y}}}=0} 574:{\displaystyle {\frac {\partial L}{\partial u}}-{\frac {\partial }{\partial x}}{\frac {\partial L}{\partial u_{x}}}-{\frac {\partial }{\partial y}}{\frac {\partial L}{\partial u_{y}}}=0} 1166: 395: 46:
The Horn-Schunck algorithm assumes smoothness in the flow over the whole image. Thus, it tries to minimize distortions in flow and prefers solutions which show more smoothness.
948: 1915:
Advantages of the Horn–Schunck algorithm include that it yields a high density of flow vectors, i.e. the flow information missing in inner parts of homogeneous objects is
439: 419: 442: 311: 284: 257: 1485: 1465: 1321: 1194: 1186: 736: 1881:{\displaystyle v^{k+1}={\overline {v}}^{k}-{\frac {I_{y}(I_{x}{\overline {u}}^{k}+I_{y}{\overline {v}}^{k}+I_{t})}{\alpha ^{2}+I_{x}^{2}+I_{y}^{2}}}} 1684:{\displaystyle u^{k+1}={\overline {u}}^{k}-{\frac {I_{x}(I_{x}{\overline {u}}^{k}+I_{y}{\overline {v}}^{k}+I_{t})}{\alpha ^{2}+I_{x}^{2}+I_{y}^{2}}}} 845: 744: 1036: 1188:
calculated in a neighborhood around the pixel at location (x,y). Using this notation the above equation system may be written
59: 1973: 1958: 1129: 316: 1022:{\displaystyle \Delta ={\frac {\partial ^{2}}{\partial x^{2}}}+{\frac {\partial ^{2}}{\partial y^{2}}}} 1033:. In practice the Laplacian is approximated numerically using finite differences, and may be written 1919:
from the motion boundaries. On the negative side, it is more sensitive to noise than local methods.
313:
are the derivatives of the image intensity values along the x, y and time dimensions respectively,
53:
which is then sought to be minimized. This function is given for two-dimensional image streams as:
50: 1928: 424: 404: 289: 262: 235: 1907:, applied to the large, sparse system arising when solving for all pixels simultaneously. 1437:{\displaystyle I_{x}I_{y}u+(I_{y}^{2}+\alpha ^{2})v=\alpha ^{2}{\overline {v}}-I_{y}I_{t}} 1310:{\displaystyle (I_{x}^{2}+\alpha ^{2})u+I_{x}I_{y}v=\alpha ^{2}{\overline {u}}-I_{x}I_{t}} 8: 1470: 1450: 1171: 721: 441:
lead to a smoother flow. This functional can be minimized by solving the associated
1900: 1489: 1030: 30: 1967: 1904: 1945: 35: 21: 935:{\displaystyle I_{y}(I_{x}u+I_{y}v+I_{t})-\alpha ^{2}\Delta v=0} 834:{\displaystyle I_{x}(I_{x}u+I_{y}v+I_{t})-\alpha ^{2}\Delta u=0} 1119:{\displaystyle \Delta u(x,y)=({\overline {u}}(x,y)-u(x,y))} 24:
is a global method which introduces a global constraint of
1940:
B.K.P. Horn and B.G. Schunck, "Determining optical flow."
1895:
denotes the next iteration, which is to be calculated and
945:
where subscripts again denote partial differentiation and
1698: 1501: 1473: 1453: 1324: 1197: 1174: 1132: 1039: 951: 848: 747: 724: 588: 454: 427: 407: 319: 292: 265: 238: 62: 1899:
is the last calculated result. This is in essence a
222:{\displaystyle E=\iint \left{{\rm {d}}x{\rm {d}}y}} 1880: 1683: 1479: 1459: 1436: 1309: 1180: 1160: 1118: 1021: 934: 833: 738:is the integrand of the energy expression, giving 730: 707: 573: 433: 413: 397:is the optical flow vector (which is to be solved 389: 305: 278: 251: 221: 1965: 421:is a regularization constant. Larger values of 180: 170: 158: 148: 443:multi-dimensional Euler–Lagrange equations 49:The flow is formulated as a global energy 41: 1966: 1161:{\displaystyle {\overline {u}}(x,y)} 390:{\displaystyle {\vec {V}}=^{\top }} 13: 1040: 1003: 993: 971: 961: 952: 920: 819: 683: 675: 663: 659: 638: 630: 618: 614: 600: 592: 549: 541: 529: 525: 504: 496: 484: 480: 466: 458: 382: 210: 200: 173: 151: 14: 1985: 1952: 1824: 1751: 1627: 1554: 1382: 1351: 1229: 1198: 1155: 1143: 1113: 1110: 1098: 1089: 1077: 1064: 1058: 1046: 904: 859: 803: 758: 378: 374: 362: 353: 341: 335: 326: 189: 145: 123: 77: 1: 1944:, vol 17, pp 185–203, 1981. 1934: 1910: 1800: 1770: 1724: 1603: 1573: 1527: 1406: 1279: 1138: 1072: 7: 1922: 10: 1990: 38:for further description). 1974:Motion in computer vision 1168:is a weighted average of 1948:available on MIT server. 1942:Artificial Intelligence 1903:method, similar to the 434:{\displaystyle \alpha } 414:{\displaystyle \alpha } 1891:where the superscript 1882: 1685: 1481: 1461: 1438: 1311: 1182: 1162: 1120: 1023: 936: 835: 732: 709: 575: 435: 415: 391: 307: 280: 253: 223: 1959:OpenCV implementation 1883: 1686: 1482: 1462: 1439: 1312: 1183: 1163: 1121: 1024: 937: 836: 733: 710: 576: 436: 416: 401:), and the parameter 392: 308: 306:{\displaystyle I_{t}} 281: 279:{\displaystyle I_{y}} 254: 252:{\displaystyle I_{x}} 224: 1696: 1499: 1471: 1451: 1322: 1195: 1172: 1130: 1037: 949: 846: 745: 722: 586: 452: 425: 405: 317: 290: 263: 236: 60: 42:Mathematical details 1929:Lucas–Kanade method 1874: 1856: 1677: 1659: 1447:which is linear in 1368: 1215: 18:Horn–Schunck method 1878: 1860: 1842: 1681: 1663: 1645: 1477: 1457: 1434: 1354: 1307: 1201: 1178: 1158: 1116: 1019: 932: 831: 728: 705: 571: 431: 411: 387: 303: 276: 249: 219: 1876: 1803: 1773: 1727: 1679: 1606: 1576: 1530: 1480:{\displaystyle v} 1460:{\displaystyle u} 1409: 1282: 1181:{\displaystyle u} 1141: 1075: 1017: 985: 731:{\displaystyle L} 697: 670: 652: 625: 607: 563: 536: 518: 491: 473: 329: 1981: 1901:Matrix splitting 1887: 1885: 1884: 1879: 1877: 1875: 1873: 1868: 1855: 1850: 1838: 1837: 1827: 1823: 1822: 1810: 1809: 1804: 1796: 1793: 1792: 1780: 1779: 1774: 1766: 1763: 1762: 1750: 1749: 1739: 1734: 1733: 1728: 1720: 1714: 1713: 1690: 1688: 1687: 1682: 1680: 1678: 1676: 1671: 1658: 1653: 1641: 1640: 1630: 1626: 1625: 1613: 1612: 1607: 1599: 1596: 1595: 1583: 1582: 1577: 1569: 1566: 1565: 1553: 1552: 1542: 1537: 1536: 1531: 1523: 1517: 1516: 1486: 1484: 1483: 1478: 1466: 1464: 1463: 1458: 1443: 1441: 1440: 1435: 1433: 1432: 1423: 1422: 1410: 1402: 1400: 1399: 1381: 1380: 1367: 1362: 1344: 1343: 1334: 1333: 1316: 1314: 1313: 1308: 1306: 1305: 1296: 1295: 1283: 1275: 1273: 1272: 1257: 1256: 1247: 1246: 1228: 1227: 1214: 1209: 1187: 1185: 1184: 1179: 1167: 1165: 1164: 1159: 1142: 1134: 1125: 1123: 1122: 1117: 1076: 1068: 1031:Laplace operator 1028: 1026: 1025: 1020: 1018: 1016: 1015: 1014: 1001: 1000: 991: 986: 984: 983: 982: 969: 968: 959: 941: 939: 938: 933: 919: 918: 903: 902: 887: 886: 871: 870: 858: 857: 840: 838: 837: 832: 818: 817: 802: 801: 786: 785: 770: 769: 757: 756: 737: 735: 734: 729: 714: 712: 711: 706: 698: 696: 695: 694: 681: 673: 671: 669: 658: 653: 651: 650: 649: 636: 628: 626: 624: 613: 608: 606: 598: 590: 580: 578: 577: 572: 564: 562: 561: 560: 547: 539: 537: 535: 524: 519: 517: 516: 515: 502: 494: 492: 490: 479: 474: 472: 464: 456: 440: 438: 437: 432: 420: 418: 417: 412: 396: 394: 393: 388: 386: 385: 331: 330: 322: 312: 310: 309: 304: 302: 301: 285: 283: 282: 277: 275: 274: 258: 256: 255: 250: 248: 247: 228: 226: 225: 220: 218: 214: 213: 204: 203: 196: 192: 188: 187: 166: 165: 144: 143: 131: 130: 121: 120: 105: 104: 89: 88: 31:aperture problem 1989: 1988: 1984: 1983: 1982: 1980: 1979: 1978: 1964: 1963: 1955: 1937: 1925: 1913: 1869: 1864: 1851: 1846: 1833: 1829: 1828: 1818: 1814: 1805: 1795: 1794: 1788: 1784: 1775: 1765: 1764: 1758: 1754: 1745: 1741: 1740: 1738: 1729: 1719: 1718: 1703: 1699: 1697: 1694: 1693: 1672: 1667: 1654: 1649: 1636: 1632: 1631: 1621: 1617: 1608: 1598: 1597: 1591: 1587: 1578: 1568: 1567: 1561: 1557: 1548: 1544: 1543: 1541: 1532: 1522: 1521: 1506: 1502: 1500: 1497: 1496: 1472: 1469: 1468: 1452: 1449: 1448: 1428: 1424: 1418: 1414: 1401: 1395: 1391: 1376: 1372: 1363: 1358: 1339: 1335: 1329: 1325: 1323: 1320: 1319: 1301: 1297: 1291: 1287: 1274: 1268: 1264: 1252: 1248: 1242: 1238: 1223: 1219: 1210: 1205: 1196: 1193: 1192: 1173: 1170: 1169: 1133: 1131: 1128: 1127: 1067: 1038: 1035: 1034: 1010: 1006: 1002: 996: 992: 990: 978: 974: 970: 964: 960: 958: 950: 947: 946: 914: 910: 898: 894: 882: 878: 866: 862: 853: 849: 847: 844: 843: 813: 809: 797: 793: 781: 777: 765: 761: 752: 748: 746: 743: 742: 723: 720: 719: 690: 686: 682: 674: 672: 662: 657: 645: 641: 637: 629: 627: 617: 612: 599: 591: 589: 587: 584: 583: 556: 552: 548: 540: 538: 528: 523: 511: 507: 503: 495: 493: 483: 478: 465: 457: 455: 453: 450: 449: 426: 423: 422: 406: 403: 402: 381: 377: 321: 320: 318: 315: 314: 297: 293: 291: 288: 287: 270: 266: 264: 261: 260: 243: 239: 237: 234: 233: 209: 208: 199: 198: 197: 183: 179: 161: 157: 139: 135: 126: 122: 116: 112: 100: 96: 84: 80: 76: 72: 61: 58: 57: 44: 12: 11: 5: 1987: 1977: 1976: 1962: 1961: 1954: 1953:External links 1951: 1950: 1949: 1936: 1933: 1932: 1931: 1924: 1921: 1912: 1909: 1889: 1888: 1872: 1867: 1863: 1859: 1854: 1849: 1845: 1841: 1836: 1832: 1826: 1821: 1817: 1813: 1808: 1802: 1799: 1791: 1787: 1783: 1778: 1772: 1769: 1761: 1757: 1753: 1748: 1744: 1737: 1732: 1726: 1723: 1717: 1712: 1709: 1706: 1702: 1691: 1675: 1670: 1666: 1662: 1657: 1652: 1648: 1644: 1639: 1635: 1629: 1624: 1620: 1616: 1611: 1605: 1602: 1594: 1590: 1586: 1581: 1575: 1572: 1564: 1560: 1556: 1551: 1547: 1540: 1535: 1529: 1526: 1520: 1515: 1512: 1509: 1505: 1476: 1456: 1445: 1444: 1431: 1427: 1421: 1417: 1413: 1408: 1405: 1398: 1394: 1390: 1387: 1384: 1379: 1375: 1371: 1366: 1361: 1357: 1353: 1350: 1347: 1342: 1338: 1332: 1328: 1317: 1304: 1300: 1294: 1290: 1286: 1281: 1278: 1271: 1267: 1263: 1260: 1255: 1251: 1245: 1241: 1237: 1234: 1231: 1226: 1222: 1218: 1213: 1208: 1204: 1200: 1177: 1157: 1154: 1151: 1148: 1145: 1140: 1137: 1115: 1112: 1109: 1106: 1103: 1100: 1097: 1094: 1091: 1088: 1085: 1082: 1079: 1074: 1071: 1066: 1063: 1060: 1057: 1054: 1051: 1048: 1045: 1042: 1013: 1009: 1005: 999: 995: 989: 981: 977: 973: 967: 963: 957: 954: 943: 942: 931: 928: 925: 922: 917: 913: 909: 906: 901: 897: 893: 890: 885: 881: 877: 874: 869: 865: 861: 856: 852: 841: 830: 827: 824: 821: 816: 812: 808: 805: 800: 796: 792: 789: 784: 780: 776: 773: 768: 764: 760: 755: 751: 727: 716: 715: 704: 701: 693: 689: 685: 680: 677: 668: 665: 661: 656: 648: 644: 640: 635: 632: 623: 620: 616: 611: 605: 602: 597: 594: 581: 570: 567: 559: 555: 551: 546: 543: 534: 531: 527: 522: 514: 510: 506: 501: 498: 489: 486: 482: 477: 471: 468: 463: 460: 430: 410: 384: 380: 376: 373: 370: 367: 364: 361: 358: 355: 352: 349: 346: 343: 340: 337: 334: 328: 325: 300: 296: 273: 269: 246: 242: 230: 229: 217: 212: 207: 202: 195: 191: 186: 182: 178: 175: 172: 169: 164: 160: 156: 153: 150: 147: 142: 138: 134: 129: 125: 119: 115: 111: 108: 103: 99: 95: 92: 87: 83: 79: 75: 71: 68: 65: 43: 40: 20:of estimating 9: 6: 4: 3: 2: 1986: 1975: 1972: 1971: 1969: 1960: 1957: 1956: 1947: 1943: 1939: 1938: 1930: 1927: 1926: 1920: 1918: 1908: 1906: 1905:Jacobi method 1902: 1898: 1894: 1870: 1865: 1861: 1857: 1852: 1847: 1843: 1839: 1834: 1830: 1819: 1815: 1811: 1806: 1797: 1789: 1785: 1781: 1776: 1767: 1759: 1755: 1746: 1742: 1735: 1730: 1721: 1715: 1710: 1707: 1704: 1700: 1692: 1673: 1668: 1664: 1660: 1655: 1650: 1646: 1642: 1637: 1633: 1622: 1618: 1614: 1609: 1600: 1592: 1588: 1584: 1579: 1570: 1562: 1558: 1549: 1545: 1538: 1533: 1524: 1518: 1513: 1510: 1507: 1503: 1495: 1494: 1493: 1491: 1490:Cramer's rule 1474: 1454: 1429: 1425: 1419: 1415: 1411: 1403: 1396: 1392: 1388: 1385: 1377: 1373: 1369: 1364: 1359: 1355: 1348: 1345: 1340: 1336: 1330: 1326: 1318: 1302: 1298: 1292: 1288: 1284: 1276: 1269: 1265: 1261: 1258: 1253: 1249: 1243: 1239: 1235: 1232: 1224: 1220: 1216: 1211: 1206: 1202: 1191: 1190: 1189: 1175: 1152: 1149: 1146: 1135: 1107: 1104: 1101: 1095: 1092: 1086: 1083: 1080: 1069: 1061: 1055: 1052: 1049: 1043: 1032: 1011: 1007: 997: 987: 979: 975: 965: 955: 929: 926: 923: 915: 911: 907: 899: 895: 891: 888: 883: 879: 875: 872: 867: 863: 854: 850: 842: 828: 825: 822: 814: 810: 806: 798: 794: 790: 787: 782: 778: 774: 771: 766: 762: 753: 749: 741: 740: 739: 725: 702: 699: 691: 687: 678: 666: 654: 646: 642: 633: 621: 609: 603: 595: 582: 568: 565: 557: 553: 544: 532: 520: 512: 508: 499: 487: 475: 469: 461: 448: 447: 446: 444: 428: 408: 400: 371: 368: 365: 359: 356: 350: 347: 344: 338: 332: 323: 298: 294: 271: 267: 244: 240: 215: 205: 193: 184: 176: 167: 162: 154: 140: 136: 132: 127: 117: 113: 109: 106: 101: 97: 93: 90: 85: 81: 73: 69: 66: 63: 56: 55: 54: 52: 47: 39: 37: 33: 32: 28:to solve the 27: 23: 19: 1941: 1916: 1914: 1896: 1892: 1890: 1446: 1029:denotes the 944: 717: 445:. These are 398: 231: 48: 45: 36:Optical Flow 29: 25: 22:optical flow 17: 15: 1946:Manuscript 1935:References 1911:Properties 51:functional 26:smoothness 1917:filled in 1831:α 1801:¯ 1771:¯ 1736:− 1725:¯ 1634:α 1604:¯ 1574:¯ 1539:− 1528:¯ 1412:− 1407:¯ 1393:α 1374:α 1285:− 1280:¯ 1266:α 1221:α 1139:¯ 1093:− 1073:¯ 1041:Δ 1004:∂ 994:∂ 972:∂ 962:∂ 953:Δ 921:Δ 912:α 908:− 820:Δ 811:α 807:− 684:∂ 676:∂ 664:∂ 660:∂ 655:− 639:∂ 631:∂ 619:∂ 615:∂ 610:− 601:∂ 593:∂ 550:∂ 542:∂ 530:∂ 526:∂ 521:− 505:∂ 497:∂ 485:∂ 481:∂ 476:− 467:∂ 459:∂ 429:α 409:α 383:⊤ 327:→ 181:‖ 174:∇ 171:‖ 159:‖ 152:∇ 149:‖ 137:α 70:∬ 1968:Category 1923:See also 1126:where 718:where 232:where 34:(see 1467:and 286:and 16:The 1893:k+1 399:for 1970:: 1492:: 259:, 1897:k 1871:2 1866:y 1862:I 1858:+ 1853:2 1848:x 1844:I 1840:+ 1835:2 1825:) 1820:t 1816:I 1812:+ 1807:k 1798:v 1790:y 1786:I 1782:+ 1777:k 1768:u 1760:x 1756:I 1752:( 1747:y 1743:I 1731:k 1722:v 1716:= 1711:1 1708:+ 1705:k 1701:v 1674:2 1669:y 1665:I 1661:+ 1656:2 1651:x 1647:I 1643:+ 1638:2 1628:) 1623:t 1619:I 1615:+ 1610:k 1601:v 1593:y 1589:I 1585:+ 1580:k 1571:u 1563:x 1559:I 1555:( 1550:x 1546:I 1534:k 1525:u 1519:= 1514:1 1511:+ 1508:k 1504:u 1475:v 1455:u 1430:t 1426:I 1420:y 1416:I 1404:v 1397:2 1389:= 1386:v 1383:) 1378:2 1370:+ 1365:2 1360:y 1356:I 1352:( 1349:+ 1346:u 1341:y 1337:I 1331:x 1327:I 1303:t 1299:I 1293:x 1289:I 1277:u 1270:2 1262:= 1259:v 1254:y 1250:I 1244:x 1240:I 1236:+ 1233:u 1230:) 1225:2 1217:+ 1212:2 1207:x 1203:I 1199:( 1176:u 1156:) 1153:y 1150:, 1147:x 1144:( 1136:u 1114:) 1111:) 1108:y 1105:, 1102:x 1099:( 1096:u 1090:) 1087:y 1084:, 1081:x 1078:( 1070:u 1065:( 1062:= 1059:) 1056:y 1053:, 1050:x 1047:( 1044:u 1012:2 1008:y 998:2 988:+ 980:2 976:x 966:2 956:= 930:0 927:= 924:v 916:2 905:) 900:t 896:I 892:+ 889:v 884:y 880:I 876:+ 873:u 868:x 864:I 860:( 855:y 851:I 829:0 826:= 823:u 815:2 804:) 799:t 795:I 791:+ 788:v 783:y 779:I 775:+ 772:u 767:x 763:I 759:( 754:x 750:I 726:L 703:0 700:= 692:y 688:v 679:L 667:y 647:x 643:v 634:L 622:x 604:v 596:L 569:0 566:= 558:y 554:u 545:L 533:y 513:x 509:u 500:L 488:x 470:u 462:L 379:] 375:) 372:y 369:, 366:x 363:( 360:v 357:, 354:) 351:y 348:, 345:x 342:( 339:u 336:[ 333:= 324:V 299:t 295:I 272:y 268:I 245:x 241:I 216:y 211:d 206:x 201:d 194:] 190:) 185:2 177:v 168:+ 163:2 155:u 146:( 141:2 133:+ 128:2 124:) 118:t 114:I 110:+ 107:v 102:y 98:I 94:+ 91:u 86:x 82:I 78:( 74:[ 67:= 64:E

Index

optical flow
aperture problem
Optical Flow
functional
multi-dimensional Euler–Lagrange equations
Laplace operator
Cramer's rule
Matrix splitting
Jacobi method
Lucas–Kanade method
Manuscript
OpenCV implementation
Category
Motion in computer vision

Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.