モチベーション

画像のコントラストをいじるトーンカーブをワンパラメータで実装したかった
画像のコントラスト調整のためのヒストグラムに適切な関数を探しています。ではシグモイドが推されまくっているが、シグモイドは値域が$[0,1]$ではない
$0 \rightarrow 0$, $1 \rightarrow 1$であってほしい

PseudoSigmoidの歴史

『シグモイド曲線の定義域が有界でないので、有界なシグモイド曲線を勝手に作ってみた』を見つけた
『シグモイド関数でコントラスト強調』も見つけた
『ガンマ曲線とシグモイド曲線による画像補正』でも似たようなことをしている

1. シグモイドをベースにちょっと縦に伸ばす

上記『シグモイド関数でコントラスト強調』や『ガンマ曲線とシグモイド曲線による画像補正』での実装

s(x) = \frac{1}{1+e^{a(b-x)}}\\
f(x) = \frac{s(x)-s(0)}{s(1)-s(0)}

結果

https://www.geogebra.org/m/ydqszvqw

いいところ

滑らか
ナイーブ

悪いところ

コントラスト下げる方には対応できなさそう
と思ったら『ガンマ曲線とシグモイド曲線による画像補正』ではできてるっぽい。どうして。。。

2. PseudoSigmoid

『シグモイド曲線の定義域が有界でないので、有界なシグモイド曲線を勝手に作ってみた』の実装

f(x) = \frac{sinh^{−1}[(2x−1)sinh(t)]}{2t} + \frac{1}{2}

結果

https://www.geogebra.org/m/ubnfhnzn

いいところ

綺麗

悪いところ

難しい
コントラスト下げる方には対応できない

3. 指数アプローチ

指数で実装してみる。
$x^t$は$(0,0)$と$(1,1)$を通るので、これを二つひっくり返して繋げれば点対称なS字カーブができる。
以下$h$は変曲点。($[0,1]$にしたければ$h=0.5$)

f(x) = \left\{
\begin{array}{ll}
h \Big( 1 + \big( \frac{x-h}{h} \big)^{2^t} \Big) & (x \geq h) \\
h \Big( 1 - \big( \frac{h-x}{h} \big)^{2^t} \Big) & (x \lt h)
\end{array}
\right.

結果

https://www.geogebra.org/m/d4zzqtqy

なんかだいぶキモい
pHジャンプだこれ

いいところ

負の値が実装できる

悪いところ

変曲点付近に張り付く

ここで気づく

f(x) = \frac{sinh^{−1}[(2x−1)sinh(t)]}{2t} + \frac{1}{2}

これ、逆関数出せるのでは？

4. 拡張PseudoSigmoid

$t$が負の時には逆関数に移行することで、コントラストを落とす方にも対応する

f(x) = \left\{
\begin{array}{ll}
0 & (t = 0)\\
\frac{sinh^{−1}[(2x−1)sinh(t)]}{2t} + \frac{1}{2} & (t \gt 0) \\
\frac{sinh[(2x−1)t]}{2sinh(t)} + \frac{1}{2} & (t \lt 0)
\end{array}
\right.

結果

https://www.geogebra.org/m/fxfwepyu

実装

そのまま書いてみた ($[0,MAX] \rightarrow [0,MAX]$)

extended_pseudosigmoid.js

function extended_pseudosigmoid(x, t, MAX){
	if (t === 0) return x;
	else if (t > 0) return (Math.asinh((2 * x / MAX - 1) * Math.sinh(t)) / 2 / t + 0.5) * MAX;
	else return (Math.sinh(t * (2 * x / MAX - 1)) / 2 / Math.sinh(t) + 0.5) * MAX;
};

結論

$\star$Sigmoidは偉大––––––！