%%% Script needs:
% oetsettings
% Vlugtenburg scripts folder added to path!


% Survey 3
walks = Read_XYZ('walking\raw\2009_08_05\Products\XYZ_RDNAP.txt');

if 1 %% rotate yes(1) or no(0)
    x0=69141.4936;      %% coordinates northern Jarkus pole (exracted from the 5 Aug 2009 survey)
    y0=447448.7293;
    theta0 =-131.65534968173412+90;  %% angle ref line (fence) in add_surveygrid!

    Xrect = -walks.data(:,1).*sind(theta0) + walks.data(:,2).*cosd(theta0); 
    Yrect =  -walks.data(:,1)*cosd(theta0) - walks.data(:,2)*sind(theta0);    
end



figure

scatter(Xrect,Yrect,3,walks.data(:,3))
hold on
%add_surveygrid
axis image
% xlim([6.69e4 6.946e4])
% ylim([4.459e5 4.484e5])
xlabel('RD coordinates [10^3 m]');
ylabel('RD coordinates [10^3 m]');
colorbar
title('Survey Grid');


