print('Mean error: {0} degrees'.format(df.az_error.mean()))
print('Standard deviation: {0} degrees'.format(df.az_error.std()))
print('Minimum error: {0} degrees'.format(df.az_error.min()))
print('Maximum error: {0} degrees'.format(df.az_error.max()))

"""
Explanation: Azimuth error
End of explanation
"""


print('Mean error: {0} degrees'.format(df.alt_error.mean()))
print('Standard deviation: {0} degrees'.format(df.alt_error.std()))
print('Minimum error: {0} degrees'.format(df.alt_error.min()))
print('Maximum error: {0} degrees'.format(df.alt_error.max()))

alt_error_plot = df.plot(x='alt1', y='alt_error', marker='.', linestyle='', color='#2980B9', figsize=(15,10))
alt_error_plot.set_xlabel('Altitude [degrees]', fontsize=16)
alt_error_plot.set_ylabel('Error [degrees]', fontsize=16)
alt_error_plot.set_title('Error in altitude between USNO and Pysolar', fontsize=16)

az_error_plot = df.plot(x='az1', y='az_error', marker='.', linestyle='', color='#2980B9', figsize=(15,10))
az_error_plot.set_xlabel('Azimuth [degrees]', fontsize=16)
az_error_plot.set_ylabel('Error [degrees]', fontsize=16)
az_error_plot.set_title('Error in azimuth between USNO and Pysolar', fontsize=16)
az_error_plot.set_xlim(0, 360)

longitude_error_plot = df.plot(kind='scatter', x='longitude', y='alt_error', color='#2980B9', figsize=(15,10))
df.plot(kind='scatter', x='longitude', y='az_error', color='#2980B9', ax=longitude_error_plot)
longitude_error_plot.set_xlabel('Longitude [degrees]', fontsize=16)
longitude_error_plot.set_ylabel('Error [degrees]', fontsize=16)
longitude_error_plot.set_title('Error between USNO and Pysolar as a function of longitude', fontsize=16)
longitude_error_plot.set_xlim(0, 360)

latitude_error_plot = df.plot(kind='scatter', x='latitude', y='alt_error', color='#2980B9', figsize=(15,10))
df.plot(kind='scatter', x='latitude', y='az_error', color='#2980B9', ax=latitude_error_plot)
latitude_error_plot.set_xlabel('Latitude [degrees]', fontsize=16)
latitude_error_plot.set_ylabel('Error [degrees]', fontsize=16)
latitude_error_plot.set_title('Error between USNO and Pysolar as a function of latitude', fontsize=16)
latitude_error_plot.set_xlim(-90, 90)

"""
Explanation: Altitude error
End of explanation
"""


