# test_extract = rs.zonal_stats(vector, raster, stats='mean', all_touched=False, weights=False, geojson_out=True)
# test_data = [i['properties'] for i in test_extract]
# test_df = pd.DataFrame(test_data)
# test_df.rename(columns = {'mean':'ad_extract'}, inplace=True)
# test_df['ad_extract'].fillna('NA', inplace=True)
# test_df.to_csv("/home/userz/Desktop/test_pd_out.csv", sep=",", encoding="utf-8", index=False)
print '-----'
Ts1 = time.time()
nmean = rs.zonal_stats(vector, raster, weights=False, all_touched=False, geojson_out=True)
T_run1 = time.time() - Ts1
print '-----'
Ts2 = time.time()
wmean = rs.zonal_stats(vector, raster, weights=False, all_touched=True, geojson_out=True)
T_run2 = time.time() - Ts2
print '-----'
Ts3 = time.time()
wmean = rs.zonal_stats(vector, raster, weights=True, all_touched=True, geojson_out=True)
# test_df.to_csv("/home/userz/Desktop/test_pd_out.csv", sep=",", encoding="utf-8", index=False)
# print '-----'
# Ts1 = time.time()
# stats1 = rs.zonal_stats(vector, raster, stats='mean', weights=False, all_touched=False)#, geojson_out=True)
# print stats1
# T_run1 = time.time() - Ts1
print '-----'
Ts2 = time.time()
stats2 = rs.zonal_stats(vector, raster2, stats='mean', weights=True, all_touched=True)#, geojson_out=True)
print stats2
T_run2 = time.time() - Ts2
print '-----'
Ts3 = time.time()
stats3 = rs.zonal_stats(vector, raster3, stats='mean', weights=True, all_touched=True)#, geojson_out=True)
print stats3
T_run3 = time.time() - Ts3
# print 'Normal Center: ' + str(T_run1)
