[{"id":13742,"web_url":"https://patchwork.libcamera.org/comment/13742/","msgid":"<CAEmqJPq=kWdPrK4yJhuycUT2j9=YCAY4ENXpcm0-3SgeKm7sXA@mail.gmail.com>","date":"2020-11-17T11:12:35","subject":"Re: [libcamera-devel] [PATCH 09/10] libcamera: src: ipa:\n\traspberrypi: agc: Improve gain update calculation for partly\n\tsaturated images","submitter":{"id":34,"url":"https://patchwork.libcamera.org/api/people/34/","name":"Naushir Patuck","email":"naush@raspberrypi.com"},"content":"Hi David,\n\n\nOn Mon, 16 Nov 2020 at 16:49, David Plowman <david.plowman@raspberrypi.com>\nwrote:\n\n> When parts of an image saturate then the image brightness no longer\n> increases linearly with increased exposure/gain. Having calculated a\n> linear gain value it's better then to try it, allowing for saturating\n> regions, and if necessary increase the gain some more. We repeat this\n> several times.\n>\n> Signed-off-by: David Plowman <david.plowman@raspberrypi.com>\n>\n\nReviewed-by: Naushir Patuck <naush@raspberrypi.com>\n\n\n> ---\n>  src/ipa/raspberrypi/controller/rpi/agc.cpp | 34 +++++++++++++++-------\n>  1 file changed, 24 insertions(+), 10 deletions(-)\n>\n> diff --git a/src/ipa/raspberrypi/controller/rpi/agc.cpp\n> b/src/ipa/raspberrypi/controller/rpi/agc.cpp\n> index 1037f968..93b46a28 100644\n> --- a/src/ipa/raspberrypi/controller/rpi/agc.cpp\n> +++ b/src/ipa/raspberrypi/controller/rpi/agc.cpp\n> @@ -422,17 +422,21 @@ void Agc::fetchAwbStatus(Metadata *image_metadata)\n>  }\n>\n>  static double compute_initial_Y(bcm2835_isp_stats *stats, AwbStatus const\n> &awb,\n> -                               double weights[])\n> +                               double weights[], double gain)\n>  {\n>         bcm2835_isp_stats_region *regions = stats->agc_stats;\n>         // Note how the calculation below means that equal weights give you\n>         // \"average\" metering (i.e. all pixels equally important).\n>         double R_sum = 0, G_sum = 0, B_sum = 0, pixel_sum = 0;\n>         for (int i = 0; i < AGC_STATS_SIZE; i++) {\n> -               R_sum += regions[i].r_sum * weights[i];\n> -               G_sum += regions[i].g_sum * weights[i];\n> -               B_sum += regions[i].b_sum * weights[i];\n> -               pixel_sum += regions[i].counted * weights[i];\n> +               double counted = regions[i].counted;\n> +               double r_sum = std::min(regions[i].r_sum * gain, ((1 <<\n> PIPELINE_BITS) - 1) * counted);\n> +               double g_sum = std::min(regions[i].g_sum * gain, ((1 <<\n> PIPELINE_BITS) - 1) * counted);\n> +               double b_sum = std::min(regions[i].b_sum * gain, ((1 <<\n> PIPELINE_BITS) - 1) * counted);\n> +               R_sum += r_sum * weights[i];\n> +               G_sum += g_sum * weights[i];\n> +               B_sum += b_sum * weights[i];\n> +               pixel_sum += counted * weights[i];\n>         }\n>         if (pixel_sum == 0.0) {\n>                 LOG(RPiAgc, Warning) << \"compute_initial_Y: pixel_sum is\n> zero\";\n> @@ -476,11 +480,21 @@ void Agc::computeGain(bcm2835_isp_stats *statistics,\n> Metadata *image_metadata,\n>         target_Y =\n>\n> config_.Y_target.Eval(config_.Y_target.Domain().Clip(lux.lux));\n>         target_Y = std::min(EV_GAIN_Y_TARGET_LIMIT, target_Y * ev_gain);\n> -       double initial_Y = compute_initial_Y(statistics, awb_,\n> -                                            metering_mode_->weights);\n> -       gain = std::min(10.0, target_Y / (initial_Y + .001));\n> -       LOG(RPiAgc, Debug) << \"Initially Y \" << initial_Y << \" target \" <<\n> target_Y\n> -                          << \" gives gain \" << gain;\n> +\n> +       // Do this calculation a few times as brightness increase can be\n> +       // non-linear when there are saturated regions.\n> +       gain = 1.0;\n> +       for (int i = 0; i < 8; i++) {\n> +               double initial_Y = compute_initial_Y(statistics, awb_,\n> +\n> metering_mode_->weights, gain);\n> +               double extra_gain = std::min(10.0, target_Y / (initial_Y +\n> .001));\n> +               gain *= extra_gain;\n> +               LOG(RPiAgc, Debug) << \"Initial Y \" << initial_Y << \"\n> target \" << target_Y\n> +                                  << \" gives gain \" << gain;\n> +               if (extra_gain < 1.01) // close enough\n> +                       break;\n> +       }\n> +\n>         for (auto &c : *constraint_mode_) {\n>                 double new_target_Y;\n>                 double new_gain =\n> --\n> 2.20.1\n>\n> _______________________________________________\n> libcamera-devel mailing list\n> libcamera-devel@lists.libcamera.org\n> https://lists.libcamera.org/listinfo/libcamera-devel\n>","headers":{"Return-Path":"<libcamera-devel-bounces@lists.libcamera.org>","X-Original-To":"parsemail@patchwork.libcamera.org","Delivered-To":"parsemail@patchwork.libcamera.org","Received":["from lancelot.ideasonboard.com (lancelot.ideasonboard.com\n\t[92.243.16.209])\n\tby patchwork.libcamera.org (Postfix) with ESMTPS id B0E85BE081\n\tfor <parsemail@patchwork.libcamera.org>;\n\tTue, 17 Nov 2020 11:12:55 +0000 (UTC)","from lancelot.ideasonboard.com (localhost [IPv6:::1])\n\tby lancelot.ideasonboard.com (Postfix) with ESMTP id 437C363324;\n\tTue, 17 Nov 2020 12:12:55 +0100 (CET)","from mail-lj1-x229.google.com (mail-lj1-x229.google.com\n\t[IPv6:2a00:1450:4864:20::229])\n\tby lancelot.ideasonboard.com (Postfix) with ESMTPS id 2DC7C63282\n\tfor <libcamera-devel@lists.libcamera.org>;\n\tTue, 17 Nov 2020 12:12:54 +0100 (CET)","by mail-lj1-x229.google.com with SMTP id 142so10162330ljj.10\n\tfor <libcamera-devel@lists.libcamera.org>;\n\tTue, 17 Nov 2020 03:12:54 -0800 (PST)"],"Authentication-Results":"lancelot.ideasonboard.com;\n\tdkim=fail reason=\"signature verification failed\" (2048-bit key;\n\tunprotected) header.d=raspberrypi.com header.i=@raspberrypi.com\n\theader.b=\"BJaTqqsE\"; dkim-atps=neutral","DKIM-Signature":"v=1; a=rsa-sha256; c=relaxed/relaxed;\n\td=raspberrypi.com; s=google;\n\th=mime-version:references:in-reply-to:from:date:message-id:subject:to\n\t:cc; bh=nZko4+CzbAXGvspG9820ih8c9qaRpRy3TBWJlQpw1qU=;\n\tb=BJaTqqsEENPnnqqdCwC6JTZKsEzxFw3+q88RWxsyhudaxvRmJA0zncrEnQQoQjmjUf\n\t7rXukGw2WWIx5x0HJ47oyqPTb9czE1URaHPjQnUAWvx9oAPWa9Gxin6ow3zSBp5d1KCy\n\tEHzgwIuMA62Hj0zRQqVuos8BRNxanruP/qOorv0JLhwt8YUumfTHTYKvLqnNesGs9HPD\n\tKF277BXZ9q2WaxQLimhD+1G7yqbL2TkgmUZEaIT8/F9vEYd8g1vHuW+c1ApB0bU5Y6Oe\n\t/ZrDDBXFQ4jgXFK400SKpqACF9TCoKo9GkhT6qZ1ILcU4vJpqGdO/tGOsiPw/YxcWUG/\n\tAHgw==","X-Google-DKIM-Signature":"v=1; a=rsa-sha256; c=relaxed/relaxed;\n\td=1e100.net; s=20161025;\n\th=x-gm-message-state:mime-version:references:in-reply-to:from:date\n\t:message-id:subject:to:cc;\n\tbh=nZko4+CzbAXGvspG9820ih8c9qaRpRy3TBWJlQpw1qU=;\n\tb=K8EGmjbaORTBENCl9Vxt3wBGLbaS/w3q28vahdg6E5NaKtDLprZ/0sQd9nEBsaQTg2\n\t44zsoOD2919Ww2zzSFSKFiorzOOCjBO7YfdlkESMEv0J01eEtAJ9iqEtJXcltq1QjCQM\n\t4sFPzJPCfzq2ZeNX192TF3RGL13BRLXBpc8idtCwM8tQ9KKQFCTtg1WGIqW66sdK8l/P\n\tK9kWIFoQZw2hRp9kxK7YUBIDm2yOZKKhaBAcCgwExcoY+icXZ/pz0rdrBTKqkfsQQcsz\n\tZZJcaW989uI5Cg/V9RQk1o8KilFgDz1cmSePp8c1NA3JChtsbBZ/9bjw28ey9cj1C6Dn\n\tOLYw==","X-Gm-Message-State":"AOAM530YZ0HShhTi0Vs5yEvigKRbE3fsl2/NzfbETuxYFTBB1E4JUTQt\n\tZJsw3aY+PivAaySzeHdR+vo/qVED8CKAr/XwZ+feKw==","X-Google-Smtp-Source":"ABdhPJxjp8EmUxVOVJVILFrAD2bG+eu95PO/r6BF3Sd6QIYudVmShoBuPbb+yD/BhoR//539qnS4EiYxOKiS9jZwJS4=","X-Received":"by 2002:a2e:9a02:: with SMTP id o2mr1606495lji.1.1605611572172; \n\tTue, 17 Nov 2020 03:12:52 -0800 (PST)","MIME-Version":"1.0","References":"<20201116164918.2055-1-david.plowman@raspberrypi.com>\n\t<20201116164918.2055-10-david.plowman@raspberrypi.com>","In-Reply-To":"<20201116164918.2055-10-david.plowman@raspberrypi.com>","From":"Naushir Patuck <naush@raspberrypi.com>","Date":"Tue, 17 Nov 2020 11:12:35 +0000","Message-ID":"<CAEmqJPq=kWdPrK4yJhuycUT2j9=YCAY4ENXpcm0-3SgeKm7sXA@mail.gmail.com>","To":"David Plowman <david.plowman@raspberrypi.com>","Subject":"Re: [libcamera-devel] [PATCH 09/10] libcamera: src: ipa:\n\traspberrypi: agc: Improve gain update calculation for partly\n\tsaturated images","X-BeenThere":"libcamera-devel@lists.libcamera.org","X-Mailman-Version":"2.1.29","Precedence":"list","List-Id":"<libcamera-devel.lists.libcamera.org>","List-Unsubscribe":"<https://lists.libcamera.org/options/libcamera-devel>,\n\t<mailto:libcamera-devel-request@lists.libcamera.org?subject=unsubscribe>","List-Archive":"<https://lists.libcamera.org/pipermail/libcamera-devel/>","List-Post":"<mailto:libcamera-devel@lists.libcamera.org>","List-Help":"<mailto:libcamera-devel-request@lists.libcamera.org?subject=help>","List-Subscribe":"<https://lists.libcamera.org/listinfo/libcamera-devel>,\n\t<mailto:libcamera-devel-request@lists.libcamera.org?subject=subscribe>","Cc":"libcamera-devel@lists.libcamera.org","Content-Type":"multipart/mixed;\n\tboundary=\"===============8101475379340968087==\"","Errors-To":"libcamera-devel-bounces@lists.libcamera.org","Sender":"\"libcamera-devel\" <libcamera-devel-bounces@lists.libcamera.org>"}}]