Conoscopic interference pattern & Matlab analysis

ソースコード：

clear;
clc;
clf;
%Import image
myimage=imread('C:\Users\Administrator\Desktop\20210627140339.jpg');
gray=rgb2gray(myimage);
size(gray)
figure(1);
imshow(gray);%gray img which is converted from your picture
[x,y]=meshgrid(1:809,1:809);
z=double(gray(1:809,1:809));
[xq,yq]=meshgrid(1:809,1:809);
zq=griddata(x,y,z,xq,yq);
zq=rot90(zq,2);
zq=fliplr(zq);
figure(2);
mesh(x,y,zq);%distribution of light strength 
axis equal

CMOSで撮った写真:

結果を分析：