local_bimodality — 2D texture op

데이터 종류: imageimage

호출: fullseye.apply(img, "local_bimodality", a=0.5, b=0.5)(2-D 는 이미지 1 장 + 스칼라 노브 2 개 a,b∈[0,1] 모델)

local_bimodality: input → output

*그림은 128×128 합성 입력에서 실제로 실행한 출력. 왼쪽이 입력, 오른쪽이 출력. 점군은 위에서 본 산점도(밝기 = z), 1-D 열은 꺾은선, 볼륨은 z 방향 최대값 투영, 동영상은 가운데 프레임, 복소 영상은 진폭이며, 그림이 되지 않는 반환값은 값 자체를 표시한다.*

노브 a 를 훑기(0.1 / 0.5 / 0.9, 다른 노브는 기본값):

local_bimodality: knob a sweep (docs site)

노브 b 를 훑기(0.1 / 0.5 / 0.9, 다른 노브는 기본값):

local_bimodality: knob b sweep (docs site)

다른 영상에서도(합성 장면 / 사진 / 동전. 윗줄이 입력, 아랫줄이 그 출력. 노브는 기본값):

local_bimodality: other inputs (docs site)

*4열은 컬러 (H,W,3) 입력. 이 op는 채널을 넘나들지 않고 컬러를 다룰 수 있다(색 채널을 세 번째 공간축으로 컨볼루션하지 않음).*

사용법

> 이 연산자의 설명은 아직 번역이 없습니다. 원문을 그대로 싣습니다.

局所の二峰性(bimodality coefficient)。閾値を 1 つも選ばずに「ここは 2 つの

値に分かれているか」だけを 0〜1 の地図で返す。

`a が窓の一辺を 9, 15, 25, 41` に振る(大きいほど推定は安定するが、細い構造は

周りに溶ける)。`b が**コントラストの床**を max(1e-5, 0.25b x 画像の値域)` に

振る ―― 窓の標準偏差がこれ未満なら 0 を返す(b=0 でも 1e-5 の絶対床は残る)。

`BC = (歪度^2 + 1) / 尖度` で、窓の 1〜4 次モーメント(箱型平均 4 回)だけから

閉形式に出る。平行移動にもスケールにも不変なので、薄い模様でも濃い模様でも

同じ値になる。参照値は解析解で決まっていて、実測もそこへ収束する(窓 41):

2 点分布(半々の 2 値)= 1.0(尖度 1 が下限)―― 段差の真上で実測 1.0000

一様分布 = 5/9 = 0.5556 ―― 一様乱数で実測 0.5598

正規分布 = 1/3 = 0.3333 ―― 正規乱数で実測 0.3306

`> 5/9` が「一様より尖った 2 山」の古典的な目安。★一峰でも 0 にはならない

―― 正規分布の床は 1/3 であって 0 ではないので、使えるのは 0.33〜1.0 の帯である。

出典と、引き継いでいる限界(Pfister, Schwarz, Janczyk, Dale & Freeman 2013,

*Frontiers in Psychology* 4:700)。★歪んだ単峰分布が偽陽性を出す: 論文の図では

明らかに単峰の分布が BC = 0.73本当に二峰の分布が BC = 0.67 で、順位が

逆転している。歪度が絶対値で大きいほど、山の数と無関係に BC が上がるため。

この op も同じ性質を持つ —— 片側に尾を引く窓(暗い地に明るい点が少しだけ、など)は

高く出る。二峰性の判定を 1 本で決めず、`local_stdotsu` の分離度と併せること。

なお論文の式は標本補正つき `(m3^2+1)/(m4 + 3(n-1)^2/((n-2)(n-3)))` で、この op は

母集団版を使う(窓は 81〜1,681 標本あり補正は無視できる。上の参照値 1/3・5/9・1.0 は

母集団版で厳密に成り立つ)。

閾値を選ぶ側の一般化との関係: Barron (2020) の Generalized Histogram Thresholding

(arXiv:2007.07350)は Otsu・Minimum Error Thresholding・重み付きパーセンタイル閾値を

特殊ケースとして包含し、どの閾値にするかを連続に補間する。この op は対になる側で、

そもそも閾値が在るかを測る。

**なぜ閾値を選ぶ op(`otsu / threshold`)の前に要るか**: 二値化はどんな画像

でも答えを返すが、山が 1 つしか無い窓では意味の無い位置で切る(`otsu` の

docstring が自ら認めている限界)。この op は切る前に「そもそも分かれるのか」を

返すので、「分かれない場所は切らない」と決められる。文字のマスが 1 色の地に 1 色の

字か、板が均一に照っているか、傷が本当に地と別の階調かの判定に使う。

床は必須。一様な面では分散が丸め屑になり、屑どうしの 3 次・4 次の比が構造に

化ける ―― 床を画像の値域に対する相対量だけで置いた版は、0.5 一色に振幅 1e-12〜1e-8

の屑を乗せただけの画像で BC の最大が 1.0000 になった(値域そのものが屑なので

相対床が効かない)。絶対床 1e-5 を併せて初めて 0 に落ちる(実測)。

正規化しない(値は BC そのもの)ので、画像間で比較できる。縁は反射

(`mode="reflect"`)。空フレームは全 0。★窓が広いのでタイル分割してはいけない

(`scale.py` に理由つきで登録済み)。

자세한 사용 가이드

gallery2d_texture_freq 패밀리 가이드

참고(샘플 데이터·문헌)

• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.

• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.

• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.

Studio에서 시도하기

아래 프로그램은 실제로 실행됨을 확인했습니다(그림과 같은 입력). Studio 도움말에서는 이 블록이 버튼이 되어 즉시 불러와 실행할 수 있습니다.

local_bimodality 0.50 0.50

▸ Load this pipeline  ·  Load & run

실행 가능한 예제(이 연산자를 실제로 호출하는 검증된 샘플)

gallery2d_texture_freqpy -3.11 examples/gallery2d_texture_freq.py

타입이 이어지는 다음 연산자(image 를 입력으로 받는 것)

identity · gaussian · mean_box · bilateral · unsharp · median · min_filter · max_filter

같은 카테고리(texture)

std_filter · local_std · scale_select_std · bootstrap_std_error · structure_tensor_orientation · structure_tensor_coherence · gabor · sk_frangi


*Provenance: ops.py — 2D 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.