1 # Copyright 2013 The Android Open Source Project 2 # 3 # Licensed under the Apache License, Version 2.0 (the "License"); 4 # you may not use this file except in compliance with the License. 5 # You may obtain a copy of the License at 6 # 7 # http://www.apache.org/licenses/LICENSE-2.0 8 # 9 # Unless required by applicable law or agreed to in writing, software 10 # distributed under the License is distributed on an "AS IS" BASIS, 11 # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 12 # See the License for the specific language governing permissions and 13 # limitations under the License. 14 15 import its.image 16 import its.caps 17 import its.device 18 import its.objects 19 import its.target 20 import pylab 21 import numpy 22 import os.path 23 import matplotlib 24 import matplotlib.pyplot 25 26 def main(): 27 """Test that a constant exposure is seen as ISO and exposure time vary. 28 29 Take a series of shots that have ISO and exposure time chosen to balance 30 each other; result should be the same brightness, but over the sequence 31 the images should get noisier. 32 """ 33 NAME = os.path.basename(__file__).split(".")[0] 34 35 THRESHOLD_MAX_OUTLIER_DIFF = 0.1 36 THRESHOLD_MIN_LEVEL = 0.1 37 THRESHOLD_MAX_LEVEL = 0.9 38 THRESHOLD_MAX_ABS_GRAD = 0.001 39 40 mults = [] 41 r_means = [] 42 g_means = [] 43 b_means = [] 44 45 with its.device.ItsSession() as cam: 46 props = cam.get_camera_properties() 47 its.caps.skip_unless(its.caps.compute_target_exposure(props) and 48 its.caps.per_frame_control(props)) 49 50 e,s = its.target.get_target_exposure_combos(cam)["minSensitivity"] 51 expt_range = props['android.sensor.info.exposureTimeRange'] 52 sens_range = props['android.sensor.info.sensitivityRange'] 53 54 m = 1 55 while s*m < sens_range[1] and e/m > expt_range[0]: 56 mults.append(m) 57 req = its.objects.manual_capture_request(s*m, e/m) 58 cap = cam.do_capture(req) 59 img = its.image.convert_capture_to_rgb_image(cap) 60 its.image.write_image(img, "%s_mult=%02d.jpg" % (NAME, m)) 61 tile = its.image.get_image_patch(img, 0.45, 0.45, 0.1, 0.1) 62 rgb_means = its.image.compute_image_means(tile) 63 r_means.append(rgb_means[0]) 64 g_means.append(rgb_means[1]) 65 b_means.append(rgb_means[2]) 66 m = m + 4 67 68 # Draw a plot. 69 pylab.plot(mults, r_means, 'r') 70 pylab.plot(mults, g_means, 'g') 71 pylab.plot(mults, b_means, 'b') 72 pylab.ylim([0,1]) 73 matplotlib.pyplot.savefig("%s_plot_means.png" % (NAME)) 74 75 # Check for linearity. For each R,G,B channel, fit a line y=mx+b, and 76 # assert that the gradient is close to 0 (flat) and that there are no 77 # crazy outliers. Also ensure that the images aren't clamped to 0 or 1 78 # (which would make them look like flat lines). 79 for chan in xrange(3): 80 values = [r_means, g_means, b_means][chan] 81 m, b = numpy.polyfit(mults, values, 1).tolist() 82 print "Channel %d line fit (y = mx+b): m = %f, b = %f" % (chan, m, b) 83 assert(abs(m) < THRESHOLD_MAX_ABS_GRAD) 84 assert(b > THRESHOLD_MIN_LEVEL and b < THRESHOLD_MAX_LEVEL) 85 for v in values: 86 assert(v > THRESHOLD_MIN_LEVEL and v < THRESHOLD_MAX_LEVEL) 87 assert(abs(v - b) < THRESHOLD_MAX_OUTLIER_DIFF) 88 89 if __name__ == '__main__': 90 main() 91 92