datainterpolate

clear all;
close all;
clc;
files = dir(fullfile('*600nm.csv'));
for i = 1:numel(files)
    filename = files(i).name;
    M = csvread(filename);
    x = M(:,1); y = M(:,2); field = M(:,3);
    % The next line is without interpolation
    fstr = filename(1:2);
    scatter = figure;
    plot3(x,y,field);
    xlabel('x');
    ylabel('y');
    saveas(scatter, sprintf('scatter%s.png',fstr));
    % Here is the interpolation
    interpolate = figure;
    % Lumerical grid was 5 nm apart
    [xi,yi] = meshgrid(-0.4:0.005:0.4, -0.4:0.005:0.4);
    Exi = griddata(x,y,field,xi,yi);
    plot3(xi,yi,Exi);
    xlabel('x');
    ylabel('y');
    zlabel('E_x');
    saveas(interpolate, sprintf('interpolate%s.png',fstr));
%     save(sprintf('%s_data.mat',fstr),sprintf('%si','xi','yi', fstr));
    close all;
    clearvars x y field xi yi
end