""" Capacitance sensor Calibration """ # 2972 y_cal = np.array([100, 500, 800, 1200, 1800, 2400, 3000, 3600, 3700]) x_cal = np.array([1901, 2176, 2393, 2668, 3095, 3496, 3914, 4330, 4403]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] print coeff_cal """ read tank data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_015_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_015_002_22_08_2015.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_010_001.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_010_002.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_007_001_24_12_2015.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # remove error value when logger is removed for downloading # http://stackoverflow.com/a/16613835/2632856 water_level_4.drop(pd.Timestamp('2015-09-10 12:00:00'), inplace=True, axis=0) # drop last and first value for i in range(1, 6, 1): eval(
__author__ = 'Tejas' import numpy as np import pandas as pd import matplotlib.pyplot as plt import checkdam.checkdam as cd from datetime import timedelta """ Read Check dam data """ block_1 = 'E:\Tejas\JK\\a.csv' water_level_1 = cd.read_correct_ch_dam_data(block_1, 0.922, -1893) block_2 = 'E:\Tejas\JK\\b.csv' water_level_2 = cd.read_correct_ch_dam_data(block_2, 0.922, -1893) block_3 = 'E:\Tejas\JK\\c.csv' water_level_3 = cd.read_correct_ch_dam_data(block_2, 0.922, -1893) for i in range(1, 3, 1): eval("water_level_{0}.drop(water_level_{0}.tail(1).index, inplace=True, axis=0)".format(i)) eval("water_level_{0}.drop(water_level_{0}.head(1).index, inplace=True, axis=0)".format(i)) # for i in range(1, 11, 1): # print "water_level_{0}".format(i) # print eval("water_level_{0}.head()".format(i)) fig = plt.figure() for i in range(1, 3, 1): x = eval("water_level_{0}.index".format(i)) y = eval("water_level_{0}['stage(m)']".format(i)) plt.plot(x, y)
import matplotlib.colors as colors from datetime import timedelta import math # latex parameters rc('font', **{'family': 'sans-serif', 'sans-serif': ['Helvetica']}) rc('text', usetex=True) plt.rc('text', usetex=True) plt.rc('font', family='serif', size=18) """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_002.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_005.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_006.CSV' water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_001.CSV' water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # water_level_7['stage(m)'] += 0.05 block_8 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_002_12_12_2014.CSV' water_level_8 = cd.read_correct_ch_dam_data(block_8, slope, intercept)
# drop date time column rain_df = rain_df.drop('Date_Time', 1) """ Check dam calibration """ y_cal = np.array([100, 400, 1000, 1600, 2250, 2750]) x_cal = np.array([1987, 2454, 3344, 4192, 5104, 5804]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_002_25_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_005_001.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_005_002_03_12_2014.CSV' water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) for i in range(1, 7, 1): eval("water_level_{0}.drop(water_level_{0}.tail(1).index, inplace=True, axis=0)".format(i)) eval("water_level_{0}.drop(water_level_{0}.head(1).index, inplace=True, axis=0)".format(i))
# raise SystemExit(0) """ Check dam calibration """ y_cal = np.array([100, 500, 1000, 2000, 2500, 3000]) x_cal = np.array([2018, 2761, 3492, 4924, 5609, 6320]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept, stage_cutoff=stage_cutoff) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_002_25_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept, stage_cutoff=stage_cutoff) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept, stage_cutoff=stage_cutoff) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept,
__author__ = 'kiruba' import numpy as np import matplotlib.pyplot as plt import pandas as pd import itertools import checkdam.checkdam as cd # calibration y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1875, 2516, 3212, 3901, 4605, 5280]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] interecept = coeff_cal[1] # depth correction based on capacitance block_1 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_012_002_21_08_2015.CSV' block_2 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_015_001.CSV' block_3 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_010_001_08_01_2015.CSV' block_1_df = cd.read_correct_ch_dam_data(block_1, calibration_slope=slope, calibration_intercept=interecept) block_2_df = cd.read_correct_ch_dam_data(block_2, calibration_slope=slope, calibration_intercept=interecept) block_3_df = cd.read_correct_ch_dam_data(block_3, calibration_slope=slope, calibration_intercept=interecept) fig = plt.figure() plt.plot(block_1_df.index, block_1_df['stage(m)'], 'g-') plt.plot(block_2_df.index, block_2_df['stage(m)'], 'g-') plt.plot(block_3_df.index, block_3_df['stage(m)'], 'g-') plt.show()
rain_df = rain_df.drop('Date_Time', 1) """ Check dam calibration """ y_cal = np.array([100, 400, 1000, 1500, 2000, 3000]) x_cal = np.array([1971, 2336, 3083, 3720, 4335, 5604]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_005_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) # block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_003_008.CSV' # water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) # water_level_2 = water_level_2[ :"2014-08-23"] # block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_001.CSV' # water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) # block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_002_12_12_2014.CSV' # water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) # block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_003_24_12_2014.CSV' # water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_004_24_1_2014.CSV' # water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) # fig = plt.figure() # for i in range(1, 7, 1): # x = eval("water_level_{0}.index".format(i))
# raise SystemExit(0) """ Check dam calibration """ y_cal = np.array([100, 500, 1000, 2000, 2500, 3000]) x_cal = np.array([2018, 2761, 3492, 4924, 5609, 6320]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept, stage_cutoff=stage_cutoff) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_002_25_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept, stage_cutoff=stage_cutoff) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept, stage_cutoff=stage_cutoff) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept, stage_cutoff=stage_cutoff) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_001.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept, stage_cutoff=stage_cutoff) block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_002_03_12_2014.CSV' water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept, stage_cutoff=stage_cutoff) block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_003_11_12_2014.CSV' water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept, stage_cutoff=stage_cutoff) block_8 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_004_23_12_2014.CSV' water_level_8 = cd.read_correct_ch_dam_data(block_8, slope, intercept, stage_cutoff=stage_cutoff) block_9 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_005_3_1_2015.CSV'
y_cal_1 = [100, 400, 1000, 1500, 2000, 3000] x_cal_1 = [1971, 2336, 3083, 3720, 4335, 5604] coeff = np.polyfit(x_cal_1, y_cal_1, 1) print coeff # print coeff[0]*3020 + coeff[1] # print coeff[0]*2810 + coeff[1] # raise SystemExit(0) fig = plt.figure() plt.plot(x_cal_1, y_cal_1, 'ro-') plt.show() raise SystemExit(0) """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_002.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_005.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_006.CSV' water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_001.CSV' water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # water_level_7['stage(m)'] += 0.05 block_8 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_002_12_12_2014.CSV' water_level_8 = cd.read_correct_ch_dam_data(block_8, slope, intercept)
# drop date time column rain_df = rain_df.drop('Date_Time', 1) """ Check dam calibration """ y_cal = np.array([100, 400, 1000, 1500, 2000, 3000]) x_cal = np.array([1971, 2336, 3083, 3720, 4335, 5604]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_005_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) # block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_003_008.CSV' # water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) # water_level_2 = water_level_2[ :"2014-08-23"] # block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_001.CSV' # water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) # block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_002_12_12_2014.CSV' # water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) # block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_003_24_12_2014.CSV' # water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_004_24_1_2014.CSV' # water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) # fig = plt.figure() # for i in range(1, 7, 1): # x = eval("water_level_{0}.index".format(i))
x_cal_1 = np.array([1894, 2563, 3298, 4049, 4794, 5548]) a_stage_1 = cd.polyfit(x_cal_1, y_cal_1, 1) coeff_cal_1 = a_stage_1['polynomial'] slope_1 = coeff_cal_1[0] intercept_1 = coeff_cal_1[1] y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1864, 2540, 3313, 4078, 4835, 5582]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope_1, intercept_1) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_00325_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope_1, intercept_1) # block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_001.CSV' # water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) # block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_002_03_12_2014.CSV' # water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) # block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_003_11_12_2014.CSV' # water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_004_23_12_2014.CSV' # water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) # block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_005_03_01_2015.CSV' # water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # Take data upto july 17, reason : mouse bite water_level_2 = water_level_2[:'2014-07-17'] # print water_level_2.tail()
a_stage_1 = cd.polyfit(x_cal_1, y_cal_1, 1) coeff_cal_1 = a_stage_1['polynomial'] slope_1 = coeff_cal_1[0] intercept_1 = coeff_cal_1[1] y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1864, 2540, 3313, 4078, 4835, 5582]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope_1, intercept_1) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_00325_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope_1, intercept_1) # block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_001.CSV' # water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) # block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_002_03_12_2014.CSV' # water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) # block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_003_11_12_2014.CSV' # water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_004_23_12_2014.CSV' # water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) # block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_005_03_01_2015.CSV' # water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # Take data upto july 17, reason : mouse bite water_level_2 = water_level_2[:'2014-07-17'] # print water_level_2.tail()
import numpy as np import matplotlib.pyplot as plt import pandas as pd import itertools import checkdam.checkdam as cd # calibration y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1875, 2516, 3212, 3901, 4605, 5280]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] interecept = coeff_cal[1] # depth correction based on capacitance block_1 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_012_002_21_08_2015.CSV' block_2 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_015_001.CSV' block_3 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_010_001_08_01_2015.CSV' block_1_df = cd.read_correct_ch_dam_data(block_1, calibration_slope=slope, calibration_intercept=interecept) block_2_df = cd.read_correct_ch_dam_data(block_2, calibration_slope=slope, calibration_intercept=interecept) block_3_df = cd.read_correct_ch_dam_data(block_3, calibration_slope=slope, calibration_intercept=interecept) fig = plt.figure() plt.plot(block_1_df.index, block_1_df['stage(m)'], 'g-') plt.plot(block_2_df.index, block_2_df['stage(m)'], 'g-') plt.plot(block_3_df.index, block_3_df['stage(m)'], 'g-') plt.show()