clear all
close all
clc

%Behoeft MCTools

filename={'cirkel_bestek.asc'};

for i=1:length(filename(:,1))
    clear noc nor xorg yorg dumgrid grid x y z
    
    fid=fopen([filename{i}]);

    fscanf(fid,'%s',1);
    noc=fscanf(fid,'%f',1);
    %
    fscanf(fid,'%s',1);
    nor=fscanf(fid,'%f',1);
    %
    fscanf(fid,'%s',1);
    xorg=fscanf(fid,'%f',1);
    %
    fscanf(fid,'%s',1);
    yorg=fscanf(fid,'%f',1);
    %
    fscanf(fid,'%s',1);
    grid_res = str2num(fscanf(fid,'%s',1));

    x = xorg:grid_res:xorg+grid_res*(noc-1);
    y = yorg+grid_res*(nor-1):-grid_res:yorg;
end
fgetl(fid);
fgetl(fid);

A = fscanf(fid,'%f',[noc nor]);
fclose(fid);
A = A';
[row,col] = find(A == 1);

cb_data.x = x(col);
cb_data.y = y(row);

A = get(h,'ContourMatrix');
cb_data.x = A(1,2:end);
cb_data.y = A(2,2:end);

figure
plot(cb_data.x,cb_data.y,'k.')
hold on

  

save cb_data cb_data