happytcj
1/21/2017 - 12:09 AM

quick rgbc spectra plot

quick rgbc spectra plot

spectra_path = '/Users/tjiang/Desktop/Processing/Mac/panel_validation/data/spectra/cal_ALSJ212143_011915.csv'
spectra = np.genfromtxt(spectra_path, delimiter=',')[:, 1:]
lambda_min = 300
lambda_max = 1100
num_peaks = 401
num_lights = 401
peaks = np.linspace(lambda_min, lambda_max, num_lights)
peaks.shape = (num_lights, 1)

# _ = plt.plot(peaks, rraw, color='red')
# _ = plt.plot(peaks, graw, color='green')
# _ = plt.plot(peaks, braw, color='blue')
# _ = plt.plot(peaks, craw, color='black')
plt.plot(peaks, spectra[:, -1]*1/np.max(spectra[:, -1]), color='black')
plt.plot(peaks, spectra[:, 0]*1/np.max(spectra[:, -1]), color='red')
plt.plot(peaks, spectra[:, 1]*1/np.max(spectra[:, -1]), color='green')
plt.plot(peaks, spectra[:, 2]*1/np.max(spectra[:, -1]), color='blue')

plt.ylim([-0.1,1.1])

plt.title('J212 Front Camera Aperture')
# _ = plt.ylabel('Spectral Irradiance (uW/cm^2)')
_ = plt.xlabel('Wavelength (nm)')
_ = plt.ylabel('Normalized Responsivity')