16 - Geoestadística No Lineal
16 - Geoestadística No Lineal
16 - Geoestadística No Lineal
GEOESTADÍSTICA NO LINEAL
• Geoestadística no Lineal
• Krigeaje Lognormal
• Krigeaje Multigaussiano
• Krigeaje con Rango Estandarizado (Rank Kriging)
• Krigeaje de Indicadores
2 / 22
Geoestadística No Lineal
La Geoestadística no Lineal son algoritmos del Krigeaje
Lineal (KS o KO) aplicados a transformaciones no lineales
de los datos originales.
3 / 22
Krigeaje Lognormal (KL): Krigeaje aplicado a transformaciones logarítmicas de los
datos.
4 / 22
Transformación inversa
Estimación Puntual:
Krigeaje Simple
Z * ( x o ) exp L ( x o ) KS LKS
2
2
Krigeaje Ordinario Z* ( x ) exp L ( x )
o o
KO
LKO
2
2
Estimación de Bloques:
Krigeaje Simple
Z * ( v ) exp L ( v ) KS KS
2
2 (v, v ) 2
Krigeaje Ordinario Z ( v ) exp L ( v )
* KO
2
KO (v, v ) 2
5 / 22
Krigeaje Multigaussiano (KM): Krigeaje aplicado a transformaciones gaussianas de los
datos.
(Es una generalización del Krigeaje Lognormal)
6 / 22
7 / 22
Pasos para realizar una Transformada Normal
1.- Ordenar los datos originales z(xi) i = 1,2,…, n, es decir, obtener z(xi)
tal que z(x1) < z(x2) < z(xn)
Z(x) Pk G(y)
1 0,135 0,05 -1,6449
2 0,275 0,15 -1,0364
3 0,570 0,25 -0,6745
4 0,910 0,35 -0,3853
5 1,310 0,45 -0,1257
6 1,325 0,55 0,1257
7 1,485 0,65 0,3853
8 2,120 0,75 0,6745
9 2,805 0,85 1,0364
10 3,925 0,95 1,6449 9 / 22
2
1 Z
Z(x) N(μ,σ) g ( y)
1
2
e
2
z2
Z(x) N(0,1) g ( y)
1
2
e
2
z2
e
y
1
Pr(Z(x) < y) G( y) 2
dz
2
10 / 22
Graficamente
11 / 22
Datos.dat
NSCORE.EXE Nscore.out
Nscore.par
Normal
Geostatistical
Score
Parameters forEnvironmental
Transform:Geostatistical
NSCORE Assessment Environmental
Software Assessment
5****************************
6
Easting
START OF PARAMETERS:
Northing
Datos.dat - file with data
Arsenic
4 0 - columns for variable and weight
Cadmium
-1.0e21 1.0e21 - trimming limits
Lead
0 - 1=transform according to specified ref. dist.
288.0
Normal 311.0
Score.850
histsmth.out value11.5
- file18.25
with reference dist.
285.6
288.0
1 2 311.0
288.0 .850
.630 11.5
8.50 18.25
30.25 for -.54852
- columns variable and weight
273.6
285.6 269.0 .630
288.0
Nscore.out 1.02 8.50
7.00
- file 30.25
20.00
for output-.06271
280.8
273.6 249.0 1.02 7.00
269.0
Tabla.trn 10.7
- file 20.00
19.25
for output-.40792
transformation table
273.6 249.0
280.8 231.0 1.02
1.01 10.7
11.2 151.5
19.25 -.50058
276.0 231.0
273.6 206.0 1.01
1.47 11.2
11.6 151.5
37.50 2.39398
…
12 / 22
Nscore.out
BACKTR.EXE Backtr.out
Backtr.par
Krigeaje de Probabilidades.
14 / 22
Krigeaje de Indicadores (KI): Krigeaje aplicado a transformaciones
de indicadores de los datos.
1 si Z ( x ) Zc
Y ( x)
0 si no
1 si Z ( x ) Zc La posibilidad de que
Y ( x) la variable Z(x) sea
0 si no
mayor o menor que Zc
1 si Z ( x ) Zc es del 1Z(x)Zc %
Y ( x)
0 si no
E{1Z(x)Zc} = Pr{Z(x)Zc} 15 / 22
Ejemplo:
Z(x) Contenido en %
19 17 11 E{Z(x)12%} = Pr (Z(x) 12%}
40m x x x
= 4/8 =0.5
14 ? 10
x • x
12 5 8 Datos transformados
x x x
1 1 0
x x x
40m
1 ? 0
x • x
1 0 0
1Y ( x )1 2 0.25 Sph (100 m ) x x x
40m
21
C31 C32 C33 C38 3 C30
x4 x0 x5
C81 C82 C83 C88 n C80 x6 x7 x8
A = -0.011 B = 0.286
18 / 22
Finalmente:
1
1 i E 1Z ( x )12
8 8
i 1Z ( xi )12
KS
Z ( xo ) 12
i 1 i 1
1 Z ( xo ) 12 KS
- 0.011 x 1.000 + 0.286 x 1.000
- 0.011 x 0.000 + 0.286 x 1.000
+ 0.286 x 0.000 - 0.011 x 1.000
+ 0.286 x 0.000 - 0.011 x 0.000
+ (1 - 1.100) x 0.5000
1 Z ( xo ) 12 KS
0.489 48.9 %
19 / 22
19 17 11 Mínimo = 5.0 100 %
x x x
Cuartil 25% = 9.0 75 %
40m
14 ? 10
x • x Mediana 50% = 11.5 50 %
Cuartil 75% = 15.5 25 %
12 5 8
x x x Máximo = 19.0 0%
40m
1 1 1 1 1 0 1 1 0
x x x x x x x x x
1 ? 1 1 ? 0 0 ? 0
x • x x • x x • x
1 0 0 1 0 0 0 0 0
x x x x x x x x x
n
Z *
KD ( x ) f i Z ( xi )
i 1
KD puede ser usado, igual que KI, para derivar una función de distribución
específica de los datos para caracterizar la incertidumbre sobre valores
desconocidos Z(x).
El KD puede ser sustituido por el KI que resuelve el mismo problema de forma no
paramétrica.
21 / 22