תקציר

Mathematica 11.0 code

ghost3=(*Import or create a square image to use*)
dim = Dimensions[ghost3][[1]]
zpoly[x_, y_, n_, m_] := -Which[m >= 0, ZernikeR[n, m, Sqrt[x^2 + y^2]]*Cos[m ArcTan[x, y]], m < 0, ZernikeR[n, -m, Sqrt[x^2 + y^2]]*Sin[-m ArcTan[x, y]] ];
Z = Table[
   N@Table[
     If[x == y == 0, 0, zpoly[x, y, n, m]], {x, -1, 1, 2/(dim - 1)}, {y, -1, 1, 2/(dim - 1)}], {n, 0, 3}, {m, -n, n, 2}];
Z = Flatten[Z, 1];
ClearAll[fpsf]
fpsf[t_] := 
  Table[RandomReal[{-1, 1}] Sin[RandomReal[{0, 1}] t + RandomReal[{0, 2 \[Pi]}] ], {10}];
p1 = Table[
   psfcoeff = fpsf[t];
   psf = psfcoeff.Z;
   psfcoefferr = RandomReal[{-0.1, 0.1}, 10];
   GraphicsRow[{ArrayPlot[Abs[Fourier[RotateRight[InverseFourier[ghost3], {128/2, 128/2}] *(Exp /@ (I 20*psf))] ], PlotRange -> All, PlotLabel -> "Distorted image"]
     ,
     MatrixPlot[psf, PlotRange -> All, Frame -> False, PlotLabel -> "Measured distortion"]
     ,
     ArrayPlot[Abs[Fourier[RotateRight[InverseFourier[ghost3], {128/2, 128/2}] *(Exp /@ (I 20*(psfcoefferr.Z) ))] ], PlotRange -> All, PlotLabel -> "Recovered image"]
     }, ImageSize -> Large]
   , {t, 0, 1, 0.01}];
ListAnimate[p1]

רישיון

אני, בעל זכויות היוצרים על עבודה זו, מפרסם בזאת את העבודה תחת הרישיון הבא:

קובץ זה זמין לפי תנאי הקדשה עולמית לנחלת הכלל CC0 1.0 של Creative Commons.

האדם ששייך יצירה להיתר הזה הקדיש את היצירה לנחלת הכלל על־ידי ויתור על כל הזכויות שלו או שלה על היצירה בכל העולם לפי חוק זכויות יוצרים, לרבות כל הזכויות הקשורות או הסמוכות כקבוע בחוק. באפשרותך להעתיק, לשנות, להפיץ, או להציג את היצירה, אפילו למטרות מסחריות, וכל זה אפילו מבלי לבקש רשות. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse