변종도
Variogram공간 통계학에서 도 2 ( s, ) 2_{ _는 공간 랜덤 필드 또는 확률 Z Z의 공간 의존도를 기술하는 함수이다.세미바리오그램 ( 1, ( \ \ _는 바리오그램의 절반입니다.
금광업계의 구체적인 사례의 경우, 바리오그램은 광구에서 채취한 2개의 시료 간의 거리에 따라 금의 비율이 얼마나 달라지는지를 나타낸다.멀리 떨어진 곳에서 채취한 샘플은 서로 가까이에서 채취한 샘플보다 더 다양합니다.
정의.
세미바리오그램 ( ) ( \ ( h )는 Matheron (1963)에 의해 최초로 정의되었으며, 에서 된 점 ( 1 \ [1][2])의 평균 제곱의 절반으로 정의되었다.
서 M M은 기하학 V(\ V의 점이고 는 해당 점의 값입니다.삼중 적분은 3차원 이상입니다.h는 관심 거리(미터 또는 km)입니다.예를 들어, 값 { f는 의 철분 함유량(위도, 경도, 표고 등의 지리적 좌표가 있는 위치 {M})을 D({의 요소로 나타낼 수 있습니다. () \ ( )를지정하면 정확한 거리에 있는 모든 점의 쌍이 샘플링됩니다.실제로 모든 곳에서 표본을 추출하는 것은 불가능하기 때문에 대신 경험적 변이도를 사용합니다.
바리안토그램은 세미바리오그램의 2배이며 두 위치( \및 2 \에서의 필드 값 차이의 분산으로 정의할 수 있으며, M{\ M에서 smathbf {}로 표기 변경되는 음표입니다. )에서 Z까지의 필드 구현(Cressie 1993):
공간랜덤 필드의이 일정할 경우 이는 1})과({2})(서s1 사이의 값의 제곱 증분에 대한 기대치에 해당합니다. 2는 공간상의 점이며 가능한 경우 시간:
정상 공정의 경우, 바리안그램과 세미바리오그램은 - 1\ h = s { s { } _ s ( 0 ,0 +) \ \ \ + h )} _ } { f ( , 0 , . 0 + h }로 나타낼 수 있다.1993):
프로세스가 더 등방성이면 바리안트 및 세미바리오그램은 함수 i ( ) : s ( 1){ _ -s { h2}-\ }의 거리(Cressie 1993):
는 으로되지 않습니다이 용어는 함수의 세 가지 형식 모두에 사용됩니다.또한 반변형도를 나타낼 때 "변형"이라는 용어를 사용하는 경우도 있고, 변형도에 기호 를 사용하는 경우도 있어 [3]혼동을 일으킨다.
특성.
(Cressie 1993, Chiles and Delfiner 1999, Wackernagel 2003)에 따르면 이론적 변이도에는 다음과 같은 특성이 있다.
- 반변형도는 정사각형의 기대치이므로 음이 아닌 ( s , 2 )0 \ ( \ } { { } \ } _ { 0 입니다.
- 세미바리오그램δ ( 1, 1) i ( ) ( ( ( ( )- ( )) ) ) \ display \ display \ disples ( \ }{ s}{ } , \ {s= \ { { ( \ ) \ ( 0 ) ( 0 ) Z(\ {s} _{
- 함수는 조건부로 음의 유한 함수인 경우에만 세미바리오그램입니다. 즉, 모든 w, N({,N})에 대해 0 {\ {의 을 받습니다유지:
- 이는 X i Z ( i 의 r () {X = \ _가 음수여야 한다는 사실에 해당한다.[disputed ]
- 정상 공정의 공분산 함수가 존재하는 경우 다음과 같이 변동도와 관련이 있습니다.
- 정지 랜덤 필드가 공간 의존성이 없는 경우(, C )이면 0(\C)= 세미바리오그램은 원점이 0인 경우를 제외하고 상수 r var가 됩니다.
- 2}}\right는 대칭 함수입니다.
- 따라서 { _)=\h)}는 짝수 함수이다.
- 랜덤 필드가 정상적이고 에르고딕한 경우 h s ( ) a ( )\_ {}\ _) →{s } } }는 필드의 분산에 대응합니다.세미바리오그램의 한계는 실이라고도 불린다.
- 따라서 세미바리오그램은 원점에서만 연속적이지 않을 수 있습니다.원점에서의 점프 높이를 너겟 또는 너겟 효과라고 부르기도 한다.
파라미터
요약하면, 다음 매개변수가 변동 그래프를 설명하는 데 자주 사용됩니다.
- n \ n:원점에서의 불연속점에서의 세미바리오그램 점프의 높이.
- s\ s: 무한대 지연 거리를 나타내는 바리안토그램의 한계.
- r \ r: 실과 바리오그램의 차이를 무시할 수 있는 거리.고정 실이 있는 모델에서는 처음 도달한 거리입니다. 점근 실이 있는 모델의 경우 일반적으로 반변위가 실의 95%에 처음 도달한 거리로 간주됩니다.
경험적 변이도
일반적으로 측정 데이터에는 샘플 Z(\ Z를 모든 위치에 사용할 수 있는 것은 아니기 때문에 경험적 변이도가 필요합니다.예를 들어 샘플 정보는 토양 샘플의 철 농도 또는 카메라의 픽셀 강도일 수 있습니다.각 샘플 정보에는 x x y(\가 지리 좌표인 2D 샘플 공간에 대한 s가 있습니다.토양 속 철의 경우 시료 공간은 3차원일 수 있습니다.시간적 변동(예: 호수 내 인 함유량)이 있는 경우 s는 4차원 벡터가 될 수 있습니다. 치수의 단위(예: 거리와 시간)가 다른 경우, 디스플레이 가 될 수 있습니다.수정된 유클리드 [4]거리를 얻기 위해 각각 d.
샘플 관측치는 Z( ) i { Z})=로 되며 샘플은k { k의 총 다른 에서 채취할 수 있다.은 샘플 z, \ 의 위치를 제공합니다.1 k \ \ 일반적으로 그림은 포인트 함수로 표시됩니다l 세미바리오그램, 정확한 거리보다는 분리 거리 h ± { h\displic이 사용되며, 일반적으로 등방성 조건을 가정한다(즉, 는 h{ h의 함수이며 중심 위치 등 다른 변수에 의존하지 않는다). 다음 각 빈에 대해 경험적 세미바리오그램 ( ± ) { {gamma }}()을 계산할 수 있습니다.
즉 로 구분된 각 포인트 쌍( 빈 폭 허용 범위 이 발견됩니다.은 N ) { ( , ) : , j ± ; , , , { N ( \ \ ) { \ { , \ } \ } \ s :이 빈에 있는 이러한 점의 는 ± {{ \delta입니다. 그러면 각 쌍 i {\ i에 대해 관측치의 차이 제곱(예: 토양 샘플 내용 또는 픽셀 강도이 발견됩니다( {이 제곱차이는 N± { N에 의해 합산되고 정규화됩니다.정의상 결과는 이 분리의 반변형도에 대해 2로 나뉩니다.
계산 속도를 위해 고유한 포인트 쌍만 필요합니다.예를 들어, 2개의 관찰 쌍 [( a , ), ( z , d) , ( z ( z _ { a , z _ { b } , ( _ { , z _ {} ) ]의 경우, 분리 ± \ h \ }만[ ( ,z , , b ) z_ display display display disteryle ) }, disteryle ) ( b, )] ( d , c )、 { ( z _ { _ { a } 、 ( _ { , z_ { c )]의 쌍은 추가 정보를 제공하지 않습니다.
변형도 모형
경험적 변동도는 모든 지연 hh에서 계산할 수 없으며, 추정의 변동으로 인해 위에서 정의한 바와 같이 유효한 변동 그래프인지 확실하지 않다.그러나 크리깅과 같은 일부 지질통계학적 방법에는 유효한 세미바리오그램이 필요합니다.따라서 적용된 지리 통계학에서 경험적 변동도는 유효성을 보장하는 모델 함수에 의해 근사된다(Chiles & Dellfiner 1999).몇 가지 중요한 모델은 다음과 같다(Chiles & Dellfiner 1999, Cressie 1993).
- 지수 변동도 모형
- 구형 변이도 모형
- 가우스 변이도 모형
aa는 범위 정의의 모호성으로 인해 참조마다 값이 다릅니다.예: a 3 a는 (Chiles&Delfiner 1999)에 사용된 값입니다. (h ){ _ { A h h 、 0 otherwise 。
논의
이 섹션은 어떠한 출처도 인용하지 않습니다.(2015년 5월 (이 및 ) |
관측치의 공간적 또는 시간적 상관관계를 설명하기 위해 지리 통계학에서 세 가지 함수가 사용된다. 즉, 상관도, 공분산 및 세미바리오그램이다.마지막은 좀 더 간단하게 바리오그램이라고도 불립니다.
변동도는 관측된 현상의 시간/공간 상관 관계 모델을 적합시키는 데 사용되므로 지리 통계학의 핵심 함수이다.따라서 하나는 가능한 공간/시간적 상관관계를 시각화하는 실험 변이도와 크리깅 함수의 가중치를 정의하는 데 추가로 사용되는 변이도 모델을 구별하는 것이다.실험 변동도는 가우스 공정의 공분산에 대한 경험적 추정치입니다.따라서 양의 확정값이 아니므로 제약이나 추가 처리 없이 크리깅에서 직접 사용할 수 없을 수 있다.이것은 왜 제한된 수의 변형도 모델만 사용되는지를 설명합니다. 가장 일반적으로 선형, 구형, 가우스 및 지수 모델입니다.
적용들
경험적 변이도는 크리깅에 의한 공간 보간에 필요한 변이도 모델의 첫 번째 추정치로 지리 통계학에서 사용된다.
- 기둥 평균 이산화탄소의 시공간 변동성에 대한 경험적 변이도를 사용하여 위성 [4]및 지상 기반 측정에 대한 일치 기준을 결정했다.
- 이질적인 물질(길소카본)[5]의 밀도에 대해 경험적 변이도를 계산했다.
- 경험적 변동도는 [6]지진의 강한 지반 운동을 관측하여 계산한다.이러한 모델은 공간적으로 분산된 [7]인프라의 지진 위험 및 손실 평가에 사용된다.
관련 개념
바리오그램의 제곱항(Z ( )- ( 2) { _는 다른 검정력으로 대체할 수 있습니다.마도그램은 절대차이 ( s)- ( 2) Z (\ _ - (\ _로 정의되며 로드그램은 절대차이 Z(의 제곱근으로 정의됩니다. Z _ 이러한 낮은 거듭제곱에 기초한 추정치는 특이치에 대한 저항성이 더 강하다고 한다.이들은 "차수 α의 변량"으로 일반화될 수 있다.
- [\ Z2}\right
이 경우 변동도는 2차, 매드그램은 1차 변동, 로드그램은 0.[8]5차 변동입니다.
여러 변수의 상관 관계를 설명하는 데 변동 도표가 사용되는 경우 교차 변동 도표라고 합니다.교차 변수 그래프는 공동 크리깅에 사용됩니다.변수가 2진수이거나 값의 클래스를 나타내는 경우에는 지시 변이도에 대해 설명합니다.인디케이터 변형도는 인디케이터 크리깅에 사용됩니다.
레퍼런스
- ^ Matheron, Georges (1963). "Principles of geostatistics". Economic Geology. 58 (8): 1246–1266. doi:10.2113/gsecongeo.58.8.1246. ISSN 1554-0774.
- ^ Ford, David. "The Empirical Variogram" (PDF). faculty.washington.edu/edford. Retrieved 31 October 2017.
- ^ Bachmaier, Martin; Backes, Matthias (2008-02-24). "Variogram or semivariogram? Understanding the variances in a variogram". Precision Agriculture. Springer Science and Business Media LLC. 9 (3): 173–175. doi:10.1007/s11119-008-9056-2. ISSN 1385-2256.
- ^ a b Nguyen, H.; Osterman, G.; Wunch, D.; O'Dell, C.; Mandrake, L.; Wennberg, P.; Fisher, B.; Castano, R. (2014). "A method for colocating satellite XCO2 data to ground-based data and its application to ACOS-GOSAT and TCCON". Atmospheric Measurement Techniques. 7 (8): 2631–2644. Bibcode:2014AMT.....7.2631N. doi:10.5194/amt-7-2631-2014. ISSN 1867-8548.
- ^ Arregui Mena, J.D.; et al. (2018). "Characterisation of the spatial variability of material properties of Gilsocarbon and NBG-18 using random fields". Journal of Nuclear Materials. 511: 91–108. Bibcode:2018JNuM..511...91A. doi:10.1016/j.jnucmat.2018.09.008.
- ^ Schiappapietra, Erika; Douglas, John (April 2020). "Modelling the spatial correlation of earthquake ground motion: Insights from the literature, data from the 2016–2017 Central Italy earthquake sequence and ground-motion simulations". Earth-Science Reviews. 203: 103139. Bibcode:2020ESRv..20303139S. doi:10.1016/j.earscirev.2020.103139.
- ^ Sokolov, Vladimir; Wenzel, Friedemann (2011-07-25). "Influence of spatial correlation of strong ground motion on uncertainty in earthquake loss estimation". Earthquake Engineering & Structural Dynamics. 40 (9): 993–1009. doi:10.1002/eqe.1074.
- ^ Olea, Ricardo A. (1991). Geostatistical Glossary and Multilingual Dictionary. Oxford University Press. pp. 47, 67, 81. ISBN 9780195066890.
추가 정보
- Cressie, N., 1993, 공간 데이터 통계, Wiley Interscience.
- Chiles, J. P., P. Delfiner, 1999, 지리통계학, 모델링 공간 불확실성, Wiley-Intercience.
- Wakernagel, H., 2003, 다변량지질통계학, 스프링거.
- P.A. Burrough 및 McDonnell, R.A., 1998, 지리정보시스템 원리.
- Isobel Clark, 1979년 응용과학출판사 실용지질통계학부
- Clark, I., 1979, 응용과학출판사 실용지질통계학.
- David, M., 1978, Elsevier Publishing, Geostatistic Ore Reserve Estimation.
- Hald, A., 1952년, Statistical Theory with Engineering Applications, John Wiley & Sons, New York.
- Journel, A. G. and Huijbregts, Ch. J., 1978 광산지질통계학, 학술신문.
외부 링크
