이미지가 공간적으로 연결된 별도의 영역으로 구성되어 있는지에 대한 통계 측정


14

이 두 가지 회색조 이미지를 고려하십시오.

강 무작위

첫 번째 이미지는 구불 구불 한 강 패턴을 보여줍니다. 두 번째 이미지는 랜덤 노이즈를 보여줍니다.

이미지에 강 패턴이 나타나는지 여부를 확인하는 데 사용할 수있는 통계적 방법을 찾고 있습니다.

강 이미지에는 강 = 높은 가치와 다른 곳 = 낮은 가치의 두 영역이 있습니다.

결과는 히스토그램이 바이 모달이라는 것입니다.

여기에 이미지 설명을 입력하십시오

따라서 강 패턴의 이미지는 분산이 높아야합니다.

그러나 위의 임의 이미지도 마찬가지입니다.

River_var = 0.0269, Random_var = 0.0310

반면에, 랜덤 이미지는 공간 연속성이 낮고, 강 이미지는 공간 연속성이 높으며, 이는 실험용 바리오 그램에서 명확하게 보여진다 : 여기에 이미지 설명을 입력하십시오

분산이 히스토그램을 하나의 숫자로 "요약"하는 것과 같은 방식으로, 나는 실험적 바리오 그램을 "요약"하는 공간적 연속성의 척도를 찾고 있습니다.

이 측정 값이 큰 지연보다 작은 지연에서 더 높은 반 분산을 "처분"하기를 원하므로 다음을 생각해 냈습니다.

 svar=h=1nγ(h)/h2

lag = 1에서 15까지만 합하면 다음과 같습니다.

River_svar = 0.0228, Random_svar = 0.0488

강 이미지에는 분산이 높아야하지만 공간 분산이 낮아야 분산 비가 발생한다고 생각합니다.

 ratio=var/svar

결과는 다음과 같습니다.

River_ratio = 1.1816, Random_ratio = 0.6337

내 생각은 이미지가 강 이미지인지 아닌지에 대한 결정 기준으로이 비율을 사용하는 것입니다. 높은 비율 (예 :> 1) = 강.

내가 어떻게 개선 할 수 있는지에 대한 아이디어가 있습니까?

모든 답변에 미리 감사드립니다!

편집 : whuber와 Gschneider의 조언에 따라 Felix Hebeler의 Matlab 함수를 사용하여 15x15 역 거리 가중치 행렬로 계산 된 두 이미지의 Morans I이 있습니다 .

River_M 랜덤 _M

결과를 각 이미지 당 하나의 숫자로 요약해야합니다. Wikipedia에 따르면 "값은 -1 (완벽한 분산 표시)에서 +1 (완벽한 상관 관계)까지입니다. 0은 임의의 공간 패턴을 나타냅니다." Morans의 제곱을 요약하면 모든 픽셀에 대해 다음을 얻습니다.

River_sumSqM = 654.9283, Random_sumSqM = 50.0785 

여기에 큰 차이가 있으므로 Morans I는 공간 연속성의 매우 좋은 척도 인 것 같습니다 :-).

강 이미지의 20,000 개의 순열에 대한이 값의 히스토그램은 다음과 같습니다. 순열 히스토그램

분명히 River_sumSqM 값 (654.9283)은 거의 없으므로 River 이미지는 공간적으로 랜덤하지 않습니다.


4
흥미로운 질문입니다. 즉각적으로 떠오르는 조언 중 하나는 바리오 그램의 단거리 부분에 초점을 맞추는 것입니다. 그것은 중요한 부분이며 이와 같은 이미지 쌍을 구별하는 최선의 작업을 수행 할 것입니다. (가까운 관련 통계는 Moran 's IGeary 's C 입니다.) 처리 할 것으로 예상되는 이미지 종류를보다 명확하게 특성화 할 수 없으면 구체적인 조언을 제공하기가 어렵습니다.
whuber

2
+1 whuber, Moran 's이 경우 좋은 출발 인 것 같습니다. 그런 다음 이미지가 얼마나 "최고"인지 확인하기 위해 순열 테스트를 고려하십시오.
Gschneider

@Gschneider 순열 테스트는 좋은 생각입니다. 특히 순열을 수행 할 필요가 없기 때문입니다. 이미지의 값을 치환하여 가변도 (또는 Moran 's I 또는 Geary 's C 등)의 분포를 계산하는 것은 간단합니다. CLT에 적용 할 값이 너무 많습니다. (예를 들어, 바리오 그램은 녹색 점처럼 보일 것입니다. 일정한 높이의 두 배는 이미지 값의 분산입니다.) "강"패턴을 "호수"또는 "레이크"와 같은 다른 패턴과 구별해야 할 경우 문제가 더욱 어려워집니다. "강"은 폭이 매우 다양 할 수 있습니다.
whuber

죄송하지만 확실하지는 않습니다. 테스트 할 이미지를 임의의 방식으로 픽셀 당 픽셀로 순열 한 다음 순열 된 이미지의 Moran I 값과 테스트중인 이미지의 I 값을 비교하라고 지시하십니까?
Andy

내가 생각한 것은 먼저 일종의 (아마도 간단한) 이웃 구조를 정의하고 Moran 's I을 계산하는 것이 었습니다. 그런 다음 K, 20 만 개의 가능한 픽셀 순열을 계산하고 각 순열에 대해 Moran 's I를 계산할 수 있습니다. 이 200,000 Moran 's I를 얻으면 관찰 된 통계가 어디에 있는지보십시오. 그러나 whuber의 방법은 더 쉽습니다. :).
Gschneider 2016 년

답변:


1

가우시안 블러는 대규모 구조를 남겨두고 높은 파수의 구성 요소를 제거하는 저역 통과 필터 역할을한다고 생각했습니다.

이미지 생성에 필요한 웨이블릿의 스케일을 볼 수도 있습니다. 모든 정보가 소규모 웨이블릿에 있다면 강이 아닐 가능성이 높습니다.

강 한 줄의 자체와 어떤 종류의 자동 상관 관계를 고려할 수 있습니다. 따라서 잡음이있는 강의 픽셀 행을 가져 와서 다음 행과의 상호 상관 함수를 찾으면 피크의 위치와 값을 모두 찾을 수 있습니다. 이 값은 랜덤 노이즈로 얻을 수있는 것보다 훨씬 높습니다. 강이있는 지역에서 무언가를 선택하지 않으면 픽셀 열은 많은 신호를 생성하지 않습니다.

http://en.wikipedia.org/wiki/Gaussian_blur

http://en.wikipedia.org/wiki/Cross-correlation


1
몇 가지 흥미로운 아이디어가 있습니다! 샘플 이미지에 대한 하나 이상의 접근 방식을 (1) 분석법의 작동 방식을 보여주고 (2) 분석법의 성능을 평가하여이 답변을 구체화하도록 설득 할 수 있습니까?
whuber

1
그것은 내 논문에서 온 것입니다. (기계 공학) 나는 비선형 진자의 위치에 가우시안 블러를 넣은 다음 컨볼 루션 방법을 사용하여 이미지에서 위치를 제거하고 분석과 비교했습니다. 몇 가지 소음 수준으로 시도했습니다. 잡음이 가우스 크기와 관련된 임계 값 미만인 한, 재구성이 매우 양호했다. link 그림 11은 재건 관련 그래프입니다. 그림 6과 방정식 2는 잡음에 대한 반응을 나타냅니다.
EngrStudent-복직 모니카

감사! 따라서이 답변을 완성하고 그것이 실제로 얼마나 효과적인지를 보여줄 수있는 좋은 위치에있는 것 같습니다. :-)
whuber

이중 의견에 대해 죄송합니다. 객체의 일부를 의사 "어머니-웨이블릿"자체로 사용하는 것이 "영리한"것이라고 생각합니다. 그것은 자기 참조 적이지만 시적이기도합니다. 자기 인접.
EngrStudent-복직 모니카

1

이것은 조금 늦었지만 하나의 제안과 하나의 관찰에 저항 할 수는 없습니다.

먼저 히스토그램 / 바로 그램 분석보다 "이미지 처리"접근 방식이 더 적합 할 것으로 생각합니다. EngrStudent의 "부드러운"제안이 올바른 길에 있지만 "흐림"부분은 비생산적입니다. 무엇이 요구되는 것은 인 에지 보존 등으로 매끄러운 양측 필터메디안 필터 . 필연적으로 비선형 이므로 이동 평균 필터보다 더 정교 합니다.

여기에 내가 의미하는 바가 있습니다. 아래는 히스토그램과 함께 두 시나리오를 근사한 두 이미지입니다. (이미지는 표준화 된 강도로 각각 100 x 100입니다).

원시 이미지 원시 이미지

그런 다음 각 이미지에 대해 5 x 5 중앙값 필터를 15 번 * 적용 하여 가장자리유지하면서 패턴을 부드럽게합니다 . 결과는 아래와 같습니다.

스무딩 된 이미지 부드러운 이미지

(* 더 큰 필터를 사용하면 가장자리 전체의 선명한 대비가 유지되지만 위치가 부드럽게됩니다.)

"강"이미지에 여전히 바이 모달 히스토그램이 포함되어 있지만 이제 2 개의 구성 요소로 멋지게 분리되었습니다 *. 한편, "백색 노이즈"이미지에는 여전히 단일 구성 요소의 히스토그램이 있습니다. (* Otsu의 방법을 통해 쉽게 임계 값을 지정 하여 마스크를 만들고 분할을 완료합니다.)


xf[y]

(미안해 줘서 미안해 ... 내 훈련은 원래 지형 학자 였어)


참고로, 일반적으로 고정 되지 않은 자연 이미지에 "variogram method"를 적용 할 때는 매우 신중한 것이 좋습니다 . 이것은 내 대답에 터치 여기 .
GeoMatt22

0

빠른 승리 일 수도 있고 전혀 효과가 없을 수도 있지만 쉽게 제거 할 수있는 제안-이미지 강도 히스토그램의 평균 대 분산의 비율을 보셨습니까?

랜덤 노이즈 이미지를 촬영하십시오. 임의로 방출 된 광자 (또는 이와 유사한)가 카메라에 닿아 생성되고 각 픽셀이 똑같이 맞을 가능성이 있고 원시 판독 값이 있다고 가정합니다 (예 : 값의 크기가 조정되지 않거나 실행 취소 할 수있는 알려진 방법으로 크기가 조정 됨). 그러면 각 픽셀의 판독 수는 포아송 분포되어야합니다. 고정 된 기간 (노출 시간)에서 여러 번 (모든 픽셀에 걸쳐) 발생하는 이벤트 수 (픽셀에 도달하는 광자)를 계산합니다.

두 개의 다른 강도 값의 강이있는 경우 두 개의 포아송 분포가 혼합되어 있습니다.

이미지를 테스트하는 가장 빠른 방법은 강도 대 평균의 비율을 보는 것입니다. 포아송 분포의 경우 평균은 분산과 거의 같습니다. 두 개의 포아송 분포가 혼합 된 경우 분산이 평균보다 큽니다. 미리 설정된 임계 값에 대해 둘의 비율을 테스트해야합니다.

매우 조잡합니다. 그러나 작동하면 이미지의 각 픽셀을 한 번만 통과하여 필요한 충분한 통계를 계산할 수 있습니다. :)

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.