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)