In [85]:

    

# plot_DRT_profiles - quick look at beach profiles using UTM coords
import numpy as np
from matplotlib import pyplot as plt
%matplotlib inline

# This .csv file was edited to remove the headers, first point, and a weird point about halfway down
fn = r"C:\\crs\\proj\\chappy\\drt_profiles\\DRT-10312015.RW5_edited.CSV"
csv = np.genfromtxt(fn, delimiter=",")

# indices of points from first profile
p0 = csv[1:13,7:10]
# cumulative along-profile distance. There are other ways to do this.
dr = np.sqrt( np.diff(p0[:,0])**2+np.diff(p0[:,1])**2 )
d0 = np.zeros(len(dr)+1)
d0[1:] = np.cumsum(dr)

p1 = csv[14:34,7:10]
dr = np.sqrt( np.diff(p1[:,0])**2+np.diff(p1[:,1])**2 )
d1 = np.zeros(len(dr)+1)
d1[1:] = np.cumsum(dr)

p2 = csv[35:50,7:10]
dr = np.sqrt( np.diff(p2[:,0])**2+np.diff(p2[:,1])**2 )
d2 = np.zeros(len(dr)+1)
d2[1:] = np.cumsum(dr)

p3 = csv[1:13,7:10]
dr = np.sqrt( np.diff(p3[:,0])**2+np.diff(p3[:,1])**2 )
d3 = np.zeros(len(dr)+1)
d3[1:] = np.cumsum(dr)

p4 = csv[51:66,7:10]
dr = np.sqrt( np.diff(p4[:,0])**2+np.diff(p4[:,1])**2 )
d4 = np.zeros(len(dr)+1)
d4[1:] = np.cumsum(dr)


fig,ax=plt.subplots(5,1,sharex=True,figsize=(12,10),sharey=True)
ax[0].plot(d0,p0[:,2],'-k')
plt.ylim([-.5,3])
ax[1].plot(d1,p1[:,2],'-k')
ax[2].plot(d2,p2[:,2],'-k')
ax[3].plot(d3,p3[:,2],'-k')
ax[4].plot(d4,p4[:,2],'-k')
plt.ylabel(r'Elevation (m)',fontsize=14)
plt.xlabel(r'Cross-shore distance (m)',fontsize=14)
fig_name = r'C:\\crs\\proj\\chappy\\drt_profiles\\drt_profiles.png'
plt.savefig(fig_name)