print path
# Get SKY REGION MASK
sky_region_mask = ImageProcess.getSkyRegion(path)
# Get Sun Orbit among three-day image stream using SUN REGION MASK
sun_orbit = []
sun_orbit = ImageProcess.getSunOrbit(path_1, path_2, path_3, sky_region_mask)
# Get Centroid of all the Detected Sun in sun_orbit.
centroid_list = []
centroid_list = ImageProcess.getCentroidList(sun_orbit)
# Fit the centroid list sample to a quadratic equation(general parabola).
theta, coeffs = ImageProcess.generalParabola(centroid_list)
# Detect sun of test image
img_centroid_radius = ImageProcess.sunDetect(filename, sky_region_mask)
# percentage = 0.0
if img_centroid_radius is not None and ImageProcess.sunDetectUsingOrbit(img_centroid_radius, theta, coeffs):
# Sun Detected!
percentage = ImageProcess.calculateCloudCoverage(filename, sky_region_mask, img_centroid_radius)
else:
percentage = ImageProcess.calculateCloudCoverage(filename, sky_region_mask)
print 'Cloud Coverage = ', percentage
cv2.imshow('final result', sky_region_mask)
cv2.waitKey(0)
