00001 /*************************************************************************** 00002 * Copyright (C) 1998-2009 by authors (see AUTHORS.txt ) * 00003 * * 00004 * This file is part of LuxRender. * 00005 * * 00006 * Lux Renderer is free software; you can redistribute it and/or modify * 00007 * it under the terms of the GNU General Public License as published by * 00008 * the Free Software Foundation; either version 3 of the License, or * 00009 * (at your option) any later version. * 00010 * * 00011 * Lux Renderer is distributed in the hope that it will be useful, * 00012 * but WITHOUT ANY WARRANTY; without even the implied warranty of * 00013 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * 00014 * GNU General Public License for more details. * 00015 * * 00016 * You should have received a copy of the GNU General Public License * 00017 * along with this program. If not, see <http://www.gnu.org/licenses/>. * 00018 * * 00019 * This project is based on PBRT ; see http://www.pbrt.org * 00020 * Lux Renderer website : http://www.luxrender.net * 00021 ***************************************************************************/ 00022 00023 // nonlinear.cpp* 00024 #include "tonemap.h" 00025 #include "paramset.h" 00026 00027 namespace lux 00028 { 00029 00030 // NonLinearOp Declarations 00031 class NonLinearOp : public ToneMap { 00032 public: 00033 // NonLinearOp Public Methods 00034 NonLinearOp(float my) { maxY = my; } 00035 virtual ~NonLinearOp() { } 00036 virtual void Map(vector<XYZColor> &xyz, int xRes, int yRes, float maxDisplayY) const { 00037 const int numPixels = xRes * yRes; 00038 float invY2; 00039 if (maxY <= 0.f) { 00040 // Compute world adaptation luminance, _Ywa_ 00041 float Ywa = 0.; 00042 for (int i = 0; i < numPixels; ++i) 00043 if (xyz[i].Y() > 0) Ywa += logf(xyz[i].Y()); 00044 Ywa = expf(Ywa / (xRes * yRes)); 00045 invY2 = 1.f / (Ywa * Ywa); 00046 } 00047 else 00048 invY2 = 1.f / (maxY * maxY); 00049 for (int i = 0; i < numPixels; ++i) { 00050 float ys = xyz[i].c[1]; 00051 xyz[i] *= (1.f + ys * invY2) / (1.f + ys); 00052 } 00053 } 00054 00055 static ToneMap *CreateToneMap(const ParamSet &ps); 00056 private: 00057 float maxY; 00058 }; 00059 00060 }//namespace lux 00061