spatial prediction of radioactivity using general regression neural ...

2 downloads 0 Views 477KB Size Report
not a very good instrument for predicting extreme and out1ier va1ues, but it a11ows the user to detect their presence. Resu1ts from the second (“joker“) data set ...
INTRODUCTION TO THE SPATIAL INTERTICLES THE DA

ARTICLES

SPATIAL PREDICTION OF RADIOACTIVITY USING GENERAL REGRESSION NEURAL NETWORK

_______________________________________________________________________________________

s~Çáã=qáãçåáå= kìÅäÉ~ê=p~ÑÉíó=fåëíáíìíÉ=Ef_o^bF=çÑ=o^pI=jçëÅçïI=oìëëá~= `çêêÉëéçåÇÉåÅÉ=íçW=î~Çáã]áÄê~ÉK~ÅKêì= bäÉå~=p~îÉäáÉî~= kìÅäÉ~ê=p~ÑÉíó=fåëíáíìíÉ=Ef_o^bF=çÑ=o^pI=jçëÅçïI=oìëëá~= `çêêÉëéçåÇÉåÅÉ=íçW=Éë~î]áÄê~ÉK~ÅKêì=

_______________________________________________________________________________________

This work describes an application of General Regression Neural Network (GRNN) to spatial predictions of radioactivity. GRNN belongs to a class of neural networks widely used for mapping continuous functions. It is based on a non-parametric (kernel) Parzen-Rosenblatt density estimator. The kernel size is the only tuning parameter, and it allows the user to implement a GRNN in an automatic mode. An important advantage of the GRNN is its very simple and fast training procedure. The most important drawbacks are high smoothing and dependence on the spatial density of the monitoring data set. The current case study is performed on the SIC2004 data sets, and the results obtained here can be compared with those obtained by the other participants using other approaches.

1. INTRODUCTION qÜÉ=pf`OMMQ=ÉñÉêÅáëÉ=ÇÉ~äë=ïáíÜ=íÉëíáåÖ=ÇáÑÑÉêÉåí=~ééêç~ÅÜÉë=íç=~ìíçã~íáÅ=ëé~íá~ä=éêÉÇáÅJ íáçåë= çÑ= ê~Çáç~ÅíáîÉ= Åçåí~ãáå~íáçåK= qÜÉëÉ= ~ééêç~ÅÜÉë= ~êÉ= ÅçåëáÇÉêÉÇ= ~ë= éçëëáÄäÉ= Å~åÇáÇ~íÉë=Ñçê=áãéäÉãÉåí~íáçå=áåíç=~=ëóëíÉã=Ñçê=ãçåáíçêáåÖ=êÉÖìä~êäó=ê~Çáç~Åíáîáíó=äÉîÉäë= çÑ= éçíÉåíá~ääó= Ç~åÖÉêçìë= çÄàÉÅíëK= pìÅÜ= ãÉíÜçÇë= êÉèìáêÉ= çåJäáåÉ= ÑìåÅíáçåáåÖI= áå= çíÜÉê= ïçêÇëW= ÜáÖÜ=êÉ~Åíáçå=ëéÉÉÇ=~åÇ=~ÄëçäìíÉ=~ìíçã~íáçåK=^äëçI=áí=áë=îÉêó=áãéçêí~åí=íÜ~í=íÜÉ= éêÉÇáÅíáçå= ÄÉ= ìåáèìÉW= ïÉ~â= Eçê= ÉîÉå= åçF= ÇÉéÉåÇÉåÅÉ= çÑ= íÜÉ= êÉëìäíë= çå= íÜÉ= áåáíá~ä= ÅçåÇáJ íáçåë=çÑ=íÜÉ=íê~áåáåÖ=éêçÅÉÇìêÉK=qÜáë=áë=åçí=íÜÉ=Å~ëÉ=Ñçê=íÜÉ=ÄÉëíJâåçïå=kÉìê~ä=kÉíïçêâ=Ó= ~= jìäíái~óÉê= mÉêÅÉéíêçå= EjimFI= ïÜáÅÜ= áë= ~= îÉêó= éçïÉêÑìä= ~åÇ= ÑäÉñáÄäÉ= áåëíêìãÉåí= Ñçê= ëé~íá~ä=éêÉÇáÅíáçåë=Eh~åÉîëâó=NVVRFK=pç=áå=íÜÉ=ÅìêêÉåí=ëíìÇó=~åçíÜÉê=âáåÇ=çÑ=åÉìê~ä=åÉíJ ïçêâë= ï~ë= ìëÉÇ= Ó= dÉåÉê~ä= oÉÖêÉëëáçå= kÉìê~ä= kÉíïçêâ= EdokkFK= ^= Å~ëÉ= ëíìÇó= ïáíÜ= dokk= Ñçê= ãçåáíçêáåÖ= ÉåîáêçåãÉåí~ä= î~êá~ÄäÉë= ï~ë= ~äêÉ~Çó= éêÉëÉåíÉÇ= Äó= h~åÉîëâó= ENVVUFK= dokk= áë= ~= åçåäáåÉ~ê= áåíÉêéçä~íçê= Ä~ëÉÇ= çå= ~= ëçìåÇ= ã~íÜÉã~íáÅ~ä= ÑçìåÇ~íáçåK= fíë= ã~áå=~Çî~åí~ÖÉ=áë=áíë=ëáãéäáÅáíó=~åÇ=Ñ~ëí=íê~áåáåÖI=ÄÉáåÖ=áåÇÉéÉåÇÉåí=çÑ=íÜÉ=áåáíá~ä=ÅçåÇáJ íáçåK= qÜÉ= ãçëí= áãéçêí~åí= Çê~ïÄ~Åâ= áë= ëãççíÜáåÖK= fí= ãÉ~åë= íÜ~í= éêÉÇáÅíáçå= ~äï~óë= çîÉêÉëíáã~íÉë= íÜÉ=ãáåáãìã= ~åÇ=~äï~óë=ìåÇÉêÉëíáã~íÉë=íÜÉ=ã~ñáãìã= î~äìÉëK=qÜìëI=áí=áë= åçí= ~= îÉêó= ÖççÇ= áåëíêìãÉåí= Ñçê= éêÉÇáÅíáåÖ= ÉñíêÉãÉ= ~åÇ= çìíäáÉê= î~äìÉëI= Äìí= áí= ~ääçïë= íÜÉ= ìëÉê= íç= ÇÉíÉÅí= íÜÉáê= éêÉëÉåÅÉK= oÉëìäíë= Ñêçã= íÜÉ= ëÉÅçåÇ= E“àçâÉêÒF= Ç~í~= ëÉí= ÅçåÑáêã= íÜáë= éçáåíK=^åÇ=áí=ëÉÉãë=íç=ÄÉ=ÉåçìÖÜ=Ñçê=~å=~ìíçã~íáÅ=ãçåáíçêáåÖ=ëóëíÉã=áÑ=áí=áë=~ÄäÉ=íç=ÇÉíÉÅí=

APPLIED GIS, VOLUME 1, NUMBER 2, 2005 MONASH UNIVERSITY EPRESS

19-1

íÜ~í=“ëçãÉíÜáåÖ=áë=ïêçåÖÒ=áå=íÜÉ=ëíìÇó=~êÉ~I=éçëáíáçå=íÜÉ=Ñäìñ=~åÇ=éêçîáÇÉ=~=éêÉäáãáå~êó= äÉîÉä=çÑ=î~äìÉëK=jçêÉ=ÇÉí~áäÉÇ=éêÉÇáÅíáçå=çÑ=ÉñíêÉãÉë=Å~å=ÄÉ=ã~ÇÉ=ä~íÉê=EáÑ=åÉÉÇÉÇF=ìëáåÖ= çíÜÉê=~ééêç~ÅÜÉë=~åÇLçê=Äó=áåÅäìÇáåÖ=~ÇÇáíáçå~ä=ãÉ~ëìêÉãÉåíëK= gìëí=~ë=Ñçê=~ää=áåíÉêéçä~íáçå=íÉÅÜåáèìÉëI=dokk=êÉëìäíë=ÇÉéÉåÇ=çå=íÜÉ=ëé~íá~ä=ÇÉåëáíó=çÑ= íÜÉ=ãçåáíçêáåÖ=Ç~í~=ëÉíK=qÜÉ=ãçêÉ=ÜçãçÖÉåÉçìëäó=íÜÉ=áåéìí=ë~ãéäÉë=~êÉ=ÇáëíêáÄìíÉÇI=íÜÉ= ÄÉííÉê=éÉêÑçêãë=íÜÉ=dokkK=aÉåëÉ=ë~ãéäÉë=ÅäìëíÉêë=Eçê=Éãéíó=ÜçäÉëF=ã~ó=Çáëíçêí=íÜÉ=ÉÑÑáJ ÅáÉåÅó= çÑ= dokkK= qÜÉ= pf`OMMQ= ÉñÉêÅáëÉ= ë~ãéäÉë= EOMM= íê~áåáåÖ= éçáåíëF= ~êÉ= ÇáëíêáÄìíÉÇ= ê~íÜÉê=ÜçãçÖÉåÉçìëäóI=ëç=íÜáë=ÑÉ~íìêÉ=çÑ=dokk=ÇçÉë=åçí=áåÑäìÉåÅÉ=íÜÉ=ÖççÇ=êÉëìäíK= qÜÉ= pf`OMMQ= ÉñÉêÅáëÉ= áåÅäìÇÉÇ= ëçãÉ= éêáçê= áåÑçêã~íáçåK= ^= èìáÅâ= äççâ= ~í= íÜÉëÉ= Ç~í~= ëÉíë= äÉÇ= íç= íÜÉ= ÅçåÅäìëáçå= íÜ~í= íÜÉ= dokk= ÅçìäÇ= ÄÉ= ~ééäáÉÇ= áå= íÜáë= ÉñÉêÅáëÉK= ^äëç= éêáçê= áåÑçêã~íáçå= ~ääçïÉÇ= ìë= íç= ÇÉÑáåÉ= áåéìí= é~ê~ãÉíÉêë= Ñçê= íÜÉ= dokk= çéíáãáë~íáçå= éêçÅÉJ ÇìêÉK=qÜÉ=ÇÉí~áäÉÇ=é~ê~ãÉíÉê=íìåáåÖ=ï~ë=éÉêÑçêãÉÇ=çå=êÉ~ä=ÉñÉêÅáëÉ=Ç~í~=áå=~å=~ìíçã~íáÅ= ãçÇÉK=

2. METHODOLOGY 2.1. KERNEL REGRESSION WITH NEURAL NETWORK: GENERAL REGRESSION NEURAL NETWORK

`çåëáÇÉê=~=åçåJäáåÉ~ê=êÉÖêÉëëáçå=éêçÄäÉãI=ÇÉëÅêáÄÉÇ=Äó=~=ãçÇÉä=ïÜçëÉ=çÄëÉêî~ÄäÉ=çìíéìí= òá=áå=êÉëéçåëÉ=íç=~å=áåéìí=îÉÅíçê=ñá=áë=ÇÉÑáåÉÇ=Äó=

z i = f ( xi ) + ε i

i = 1,2K, N

(1)

ïÜÉêÉ=ÑEñiF=áë=~=“ëãççíÜ=ÅìêîÉÒI=~åÇ= ε i áë=~=ê~åÇçã=êÉëáÇì~ä=Çê~ïå=Ñêçã=~=ïÜáíÉ=åçáëÉ= éêçÅÉëë=çÑ=òÉêç=ãÉ~å=~åÇ=î~êá~åÅÉ=σ K=qÜ~í=áë= 2

E[ε i ] = 0 ∀i

and

⎧σ 2 i = j E[ε i ε j ] = ⎨ ⎩0 otherwise

(2)

qÜÉ=éêçÄäÉã=áë=íç=êÉÅçåëíêìÅí=íÜÉ=ìåÇÉêäóáåÖ=ÑìåÅíáçåI=ÑEñáFI=ÖáîÉå=íÜÉ=íê~áåáåÖ=ë~ãéäÉ=

{( x i , z i )}iN=1 K=^ë=~=êÉ~ëçå~ÄäÉ=Éëíáã~íÉ=çÑ=~å=ìåâåçïå=êÉÖêÉëëáçå=ÑìåÅíáçå=ÑEñFI=çåÉ=ã~ó= í~âÉ=~=ãÉ~å=çÑ=çÄëÉêî~íáçåë=EáKÉK=î~äìÉë=çÑ=íÜÉ=ãçÇÉä=çìíéìí=òF=Ñêçã=~=åÉáÖÜÄçêÜççÇ=çÑ= íÜÉ=éçáåí=ñK=qÜáë=~ééêç~ÅÜ=áë=ëìÅÅÉëëÑìä=áÑ=íÜÉ=äçÅ~ä=~îÉê~ÖÉ=áë=ÅçåÑáåÉÇ=íç=çÄëÉêî~íáçåë=áå= ~=ëã~ää=åÉáÖÜÄçêÜççÇ=EáKÉK=êÉÅÉéíáîÉ=ÑáÉäÇF=çÑ=íÜÉ=éçáåí=ñ=~ë=çÄëÉêî~íáçåë=ÅçêêÉëéçåÇáåÖ=íç= éçáåíë= ~ï~ó= Ñêçã= ñ= ïáää= ÖÉåÉê~ääó= Ü~îÉ= ÇáÑÑÉêÉåí= ãÉ~å= î~äìÉëK= jçêÉ= éêÉÅáëÉäóI= ÑEñF= áë= Éèì~ä=íç=íÜÉ=ÅçåÇáíáçå~ä=ãÉ~å=çÑ=ò=ÖáîÉå=ñ=EáKÉKI=íÜÉ=êÉÖêÉëëáçå=çÑ=ò=çå=ñFK=rëáåÖ=íÜÉ=ÑçêJ ãìä~=Ñçê=íÜÉ=ÉñéÉÅí~íáçå=çÑ=~=ê~åÇçã=î~êá~ÄäÉI=áí=Å~å=ÄÉ=ïêáííÉå= ∞

f ( x ) = E[ z | x ] =

∫ zf

Z

( z | x)dz

(3)

−∞

ïÜÉêÉ= ÑwEòöñF= áë= íÜÉ= ÅçåÇáíáçå~ä= éêçÄ~Äáäáíó= ÇÉåëáíó= ÑìåÅíáçå= EéKÇKÑKF= çÑ= òI= ÖáîÉå= ñK= cêçã= éêçÄ~Äáäáíó=íÜÉçêóI=ïÉ=Ü~îÉ=

f Z ( z | x) =

19-2

f X ,Z ( x , z ) f X ( x)

(4)

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

ïÜÉêÉ=ÑuIIwEñIòF=áë=íÜÉ=àçáåí=éKÇKÑK=çÑ=ñ=~åÇ=ò=~åÇ=ÑuEñF=áë=íÜÉ=ã~êÖáå~ä=éKÇKÑK=çÑ=ñK=eÉåÅÉI= ìëáåÖ=EQF=áå=EPFI=ïÉ=çÄí~áå=íÜÉ=ÑçääçïáåÖ=Ñçêãìä~=Ñçê=íÜÉ=ÖÉåÉê~ä=êÉÖêÉëëáçå=ÑìåÅíáçå= ∞

f ( x) =

∫ zf

X ,Z

( x , z )dz

−∞

(5)

f X ( x)

kçíÉ=íÜ~í=ïÉ=Å~å=ÇÉêáîÉ=íÜÉ=ã~êÖáå~ä=éKÇKÑK=ÑuEñF=Äó=áåíÉÖê~íáåÖ=íÜÉ=àçáåí=éKÇKÑK=ÑuIwEñIòF= çîÉê=ÇòW= ∞

f X ( x) =

∫f

X ,Z

( x , z )dz

(6)

−∞

lìê=é~êíáÅìä~ê=áåíÉêÉëí=áë=áå=~=ëáíì~íáçå=ïÜÉêÉ=íÜÉ=àçáåí=éKÇKÑK=ÑuIwEñIòF=áë=ìåâåçïåK=^ää= N

íÜ~í= ïÉ= Ü~îÉ= ~î~áä~ÄäÉ= áë= ~= ëÉí= çÑ= íê~áåáåÖ= ë~ãéäÉëÄ= {( x i , z i )}i =1 K= qç= Éëíáã~íÉ= ÑuIwEñIòF= ~åÇ= ÑuEñF= íÜÉêÉÑçêÉI= ïÉ= ã~ó= ìëÉ= ~= åçåJé~ê~ãÉíêáÅ= Éëíáã~íçê= âåçïå= ~ë= íÜÉ= m~êòÉåJ oçëÉåÄä~íí= ÇÉåëáíó= Éëíáã~íçê= EoçëÉåÄä~íí= NVRSI= NVTMX= m~êòÉå= NVSOFK= qÜÉ= Ä~ëáë= Ñçê= íÜÉ= Ñçêãìä~íáçå=áë=~=âÉêåÉäI=ÇÉåçíÉÇ=Äó=hEñFI=ïÜáÅÜ=Ü~ë=éêçéÉêíáÉë=ëáãáä~ê=íç=íÜçëÉ=~ëëçÅá~íÉÇ= ïáíÜ=~=éKÇKÑK=E~=ÅçåíáåìçìëI=ÄçìåÇÉÇI=~åÇ=êÉ~ä=ÑìåÅíáçå=çÑ=ñ=ïáíÜ=íÜÉ=íçí~ä=îçäìãÉ=ìåÇÉê= íÜÉ=ëìêÑ~ÅÉ=Éèì~ä=íç=~=ìåáíFK=tÜÉêÉ~ë=íÜÉ=âÉêåÉä=áë=~=ëóããÉíêáÅ=ÑìåÅíáçå=~êçìåÇ=íÜÉ=çêáÖáå= ïÜÉêÉ=áí=~íí~áåë=áíë=ã~ñáãìã=î~äìÉI=íÜáë=áë=åçí=~=êÉèìáêÉãÉåí=Ñçê=~=éKÇKÑK= ^ëëìãáåÖ=íÜ~í=EñiIòiF=~êÉ=áåÇÉéÉåÇÉåíI=áÇÉåíáÅ~ääó=ÇáëíêáÄìíÉÇ=ê~åÇçã=îÉÅíçêëI=ïÉ=ã~ó= Ñçêã~ääó=ÇÉÑáåÉ=íÜÉ=m~êòÉåJoçëÉåÄä~íí=ÇÉåëáíó=Éëíáã~íçê=çÑ=àçáåí=éKÇKÑK=ÑuIwEñIòF=áå=íÜÉ=àçáåí= áåéìíJçìíéìí=ëé~ÅÉ=~ë=Em~êòÉå=NVSOX=oçëÉåÄä~íí=NVRSFW=

ˆf X ,Z ( x , z ) =

1

σ

p +1

N

⎛ x − xi ⎞ ⎛ z − zi ⎞ ⎟K ⎜ ⎟ σ ⎠ ⎝ σ ⎠

∑ K ⎜⎝ N i =1

for

x ∈ ℜ p and z ∈ ℜ1

(7)

ïÜÉêÉ=íÜÉ=ëãççíÜáåÖ=é~ê~ãÉíÉê=σ=áë=~=éçëáíáîÉ=åìãÄÉê=Å~ääÉÇ=Ä~åÇïáÇíÜ=çê=ëáãéäó=ïáÇíÜK= fí= Åçåíêçäë= íÜÉ= ëáòÉ= çÑ= íÜÉ= âÉêåÉä= hEñFK= ^ëëìãáåÖ= íÜ~í= íÜÉ= ëÅ~äáåÖ= Ñçê= áåéìí= ~åÇ= çìíéìí= î~êá~ÄäÉë=áë=íÜÉ=ë~ãÉ=ÇçÉë=åçí=~ÑÑÉÅí=íÜÉ=êÉëìäíI=~ë=ïáää=ÄÉ=ëÉÉå=ä~íÉêK=rëáåÖ=ESF=çê=áåíÉÖê~íJ áåÖ=

ˆf X ,Z ( x , z ) = áå= ETF= ïáíÜ= êÉëéÉÅí= íç= òI= ïÉ= Å~å= ÇÉÑáåÉ= íÜÉ= m~êòÉåJoçëÉåÄä~íí= ÇÉåëáíó=

Éëíáã~íçê=çÑ=éKÇKÑK=ÑuEñF=~ë=

fˆX ( x) =

1 σ pN

N

⎛ x − xi ⎞ ⎟ σ ⎠

∑ K ⎜⎝ i =1

for

x ∈ℜp

(8)

kçíÉ= íÜ~í= ~å= áãéçêí~åí= éêçéÉêíó= çÑ= íÜáë= Éëíáã~íçê= áë= íÜ~í= áí= áë= ÅçåëáëíÉåí= EáKÉK= ~ëóãéíçJ íáÅ~ääó=ìåÄá~ëÉÇF=áå=íÜÉ=ëÉåëÉ=íÜ~í=áÑ=σZσEkF=áë=ÅÜçëÉå=~ë=~=ÑìåÅíáçå=çÑ=k=ëìÅÜ=íÜ~í=

lim σ ( N ) = 0 N →∞

íÜÉå=

lim E[ fˆX ( x)] = f X ( x) N →∞

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

19-3

cçê=íÜáë=ä~ííÉê=Éèì~íáçå=íç=ÜçäÇI= x =ëÜçìäÇ=ÄÉ=~=éçáåí=çÑ=Åçåíáåìáíó=Ñçê=

fˆX ( x ) K=

qÜìëI=ìëáåÖ=ETF=ïÉ=Å~å=êÉïêáíÉ=íÜÉ=åìãÉê~íçê=çÑ=ERF=~ë= ∞

∫ zˆf X ,Z ( x , z )dz =

−∞



⎛ x − xi ⎞ ⎛ z − zi ⎞ K⎜ ⎟ ∫ zK ⎜ ⎟dz ∑ p +1 σ N i =1 ⎝ σ ⎠ −∞ ⎝ σ ⎠ N

1

(9)

`Ü~åÖáåÖ=íÜÉ=î~êá~ÄäÉ=çÑ=áåíÉÖê~íáçå=Äó=ëÉííáåÖ=~=Z=Eò=Ó=òi=F=L=σ=I=~åÇ=ìëáåÖ=íÜÉ=ëóããÉíJ êáÅ=éêçéÉêíó=çÑ=íÜÉ=âÉêåÉä=hE⋅FI=ïÉ=çÄí~áå=íÜÉ=êÉëìäí= ∞

∫ zˆf X ,Z ( x , z )dz =

−∞

1 σ pN

N

⎛ x − xi ⎞ ⎟ σ ⎠

∑ z K ⎜⎝ i =1

i

(10)

qÜÉ=ÇÉåçãáå~íçê=çÑ=ERF=áë=EUFK=`~åÅÉäáåÖ=Åçããçå=íÉêãëI=ïÉ=çÄí~áå=íÜÉ=ÑçääçïáåÖ=ÉëJ íáã~íÉ=çÑ=íÜÉ=êÉÖêÉëëáçå=ÑìåÅíáçå=ÑEñFW=

⎛ x − xi ⎞ ⎟ σ ⎠ ˆf ( x ) = i =1 N ⎛ x − xi ⎞ K⎜ ⎟ ∑ ⎝ σ ⎠ i =1 N

∑ z K ⎜⎝ i

(11)

tÉ=Å~å=ÇÉÑáåÉ=íÜÉ=åçêã~äáòÉÇ=ïÉáÖÜíáåÖ=ÑìåÅíáçå=~ë=~=ÑìåÅíáçå=çÑ=ñ=

⎛ x − xj ⎞ ⎟ K ⎜⎜ σ ⎟⎠ ⎝ W j ( x) = N ⎛ x − xi ⎞ K⎜ ⎟ ∑ ⎝ σ ⎠ i =1

j = 1,2K, N

(12)

qÜÉ=ÇÉåçãáå~íçê=çÑ=ENOF=ÖáîÉë=ìë=íÜÉ=åçêã~äáò~íáçå=éêçéÉêíó= N

∑W ( x ) = 1 j =1

j

∀x

(13)

kçï=ïÉ=Å~å=êÉïêáíÉ=ENPF=áå=íÜÉ=ëáãéäáÑáÉÇ=Ñçêã= N

ˆf ( x ) = ∑ W j ( x )z j

(14)

j =1

få=íÜáë=ÑçêãI=Éèì~íáçå=ENNF=~åÇ=éêçéÉêíó=ENPF=ÇÉëÅêáÄÉ=

fˆ ( x) =~ë=~=ïÉáÖÜíÉÇ=~îÉê~ÖÉ=çÑ=

íÜÉ= ò= çÄëÉêî~íáçåëK= qÜÉ= é~êíáÅìä~ê= Ñçêã= çÑ= ïÉáÖÜíáåÖ= ÑìåÅíáçå= tjEñF= ÖáîÉå= áå= ENOF= ï~ë= çêáÖáå~ääó= éêçéçëÉÇ= Äó= k~Ç~ê~ó~= ENVSQF= ~åÇ= t~íëçå= ENVSQFK= ^ÅÅçêÇáåÖäóI= íÜÉ= ~ééêçñáJ ã~íáçå=ÑìåÅíáçå=áå=ENQF=áë=çÑíÉå=Å~ääÉÇ=íÜÉ=k~Ç~ê~ó~Jt~íëçå=hÉêåÉä=oÉÖêÉëëáçå=bëíáã~íçê= EkthobFI=çê=êÉÑÉêêáåÖ=íç=íÜÉ=ëí~êíáåÖ=éçáåíë=EPF=~åÇ=ERFI=íÜÉ=ÖÉåÉê~ä=êÉÖêÉëëáçå=Éëíáã~íçêK= qÜÉêÉ=áë=~=î~êáÉíó=çÑ=éçëëáÄäÉ=âÉêåÉä=ÑìåÅíáçåë=hE⋅FK=qÜÉ=ãçëí=ïáÇÉäó=ìëÉÇ=~ãçåÖ=íÜÉã= áë=~=ãìäíáî~êá~íÉ=d~ìëëá~å=âÉêåÉäW=

19-4

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

⎛ x 2⎞ 1 ⎟ x = (x , K , x ) K ( x) = ∏ K i ( xi ) = exp⎜ − 1 p p/2 ⎜ ⎟ ( 2 ) 2 π i =1 ⎝ ⎠ p

(15)

ïáíÜ= ööñöö= ÜÉêÉ= ~åÇ= ÑìêíÜÉê= ãÉ~åáåÖ= íÜÉ= åçêã= ÇÉÑáåÉÇ= áå= íÜÉ= éJëáòÉ= ëé~ÅÉK= `ÉåíÉêáåÖ= íÜÉ= âÉêåÉä= çå= ~= Ç~í~= éçáåí= ñi= ~åÇ= ëÅ~äáåÖ= íÜÉ= ïáÇíÜ= ïáíÜ= íÜÉ= ëãççíÜáåÖ= é~ê~ãÉíÉê= σI= ïÉ== çÄí~áå=íÜÉ=ÑçääçïáåÖ=Ñçêã=

⎛ x − xi 1 ⎛ x − xi ⎞ ⎜− K⎜ exp ⎟= 2 p/2 ⎜ 2σ 2 ⎝ σ ⎠ ( 2πσ ) ⎝

2

⎞ ⎟ i = 1,2K, N ⎟ ⎠

(16)

táíÜ=íÜÉ=âÉêåÉä=ENSFI=íÜÉ=k~Ç~ê~ó~Jt~íëçå=hÉêåÉä=oÉÖêÉëëáçå=bëíáã~íçê=äççâë=äáâÉW=

⎛ x − xi 2 ⎞ ⎟ y i exp⎜ − ∑ 2 ⎟ ⎜ 2 σ i =1 ⎠ ⎝ fˆ ( x) = 2 N ⎛ x − xi ⎞ ⎟ exp⎜ − ∑ ⎜ 2σ 2 ⎟ i =1 ⎝ ⎠ N

(17)

ïÜÉêÉ=íÜÉ=ÇÉåçãáå~íçê=íÉêã=êÉéêÉëÉåíë=íÜÉ=m~êòÉåJoçëÉåÄä~íí=ÇÉåëáíó=Éëíáã~íçê=EïáíÜçìí= åçêã~äáò~íáçåF= ~ë= ~= ëìã= çÑ= k= ãìäíáî~êá~íÉ= d~ìëëá~å= ÇáëíêáÄìíáçåë= ÅÉåíÉêÉÇ= çå= íÜÉ= Ç~í~= éçáåíë=ñi=EpéÉÅÜí=NVVNFK=qç=ÄÉ=ìåÄá~ëÉÇ=íÜáë=íÉêã=ëÜçìäÇ=ÄÉ=ÇáîáÇÉÇ=Äó=íÜÉ=åçêã~äáòáåÖ= Åçåëí~åí=

(2πσ ) 2

p 2

N K= qÜÉ= Ñçêã= ENTF= éêÉëÉåíë= íÜÉ= k~Ç~ê~ó~Jt~íëçå= d~ìëëá~å= hÉêåÉä=

oÉÖêÉëëáçå=bëíáã~íçêK= qÜÉ=Ñçêã=ENTF=Å~å=ÄÉ=îáÉïÉÇ=~ë=~=åçêã~äáòÉÇ=Ñçêã=çÑ=íÜÉ=d~ìëëá~å=ê~Çá~ä=Ä~ëáë=ÑìåÅJ íáçåë=E_áëÜçé=NVVRI=éK=NTTF=ÇÉÑáåÉÇ=áå=íÜÉ=áåéìí=ëé~ÅÉK=b~ÅÜ=Ä~ëáë=ÑìåÅíáçå=áë=ÅÉåíÉêÉÇ=çå= ~=Ç~í~=éçáåíI=~åÇ=íÜÉ=ÅçÉÑÑáÅáÉåíë=áå=íÜÉ=Éñé~åëáçå=~êÉ=ÖáîÉå=Äó=íÜÉ=í~êÖÉí=î~äìÉë=òiK= 2.2. NWKRE AS A NEURAL NETWORK: GRNN

qÜÉ=kthob=áå=íÜÉ=Ñçêã=ENTF=Å~å=ÄÉ=áãéäÉãÉåíÉÇ=áå=íÉêãë=çÑ=~=åÉìê~ä=åÉíïçêâ=EpéÉÅÜí= NVVNFI=~åÇ=íÜÉ=ÅçêêÉëéçåÇáåÖ=åÉìê~ä=åÉíïçêâ=áë=Å~ääÉÇ=~=dÉåÉê~ä=oÉÖêÉëëáçå=kÉìê~ä=kÉíJ ïçêâ= EdokkFK= cáÖìêÉ= N= áë= ~= Öê~éÜáÅ~ä= êÉéêÉëÉåí~íáçå= çÑ= dokk= áå= íÉêãë= çÑ= ~= åÉìê~ä= åÉíïçêâK= dokk=Åçåëáëíë=çÑ=Ñçìê=ä~óÉêëW=áåéìíI=é~ííÉêåI=ëìãã~íáçå=~åÇ=çìíéìíK=qÜÉ=áåéìí=ä~óÉê= íê~åëÑÉêë= ~å= áåéìí= ëáÖå~ä= Ó= ~= éJÇáãÉåëáçå~ä= îÉÅíçê= ñI= áåíç= íÜÉ= åÉñí= é~ííÉêå= ä~óÉêK= få= íÜÉ= ëé~íá~ä=ã~ééáåÖ=Å~ëÉ=éZOK=qÜÉ=åìãÄÉê=çÑ=åÉìêçåë=EâÉêåÉäëF=áå=íÜÉ=é~ííÉêå=ä~óÉê=áë=Éèì~ä= íç=íÜÉ=åìãÄÉê=çÑ=íê~áåáåÖ=ë~ãéäÉë=EkFI=É~ÅÜ=åÉìêçå=ÅçêêÉëéçåÇáåÖ=íç=~=íê~áåáåÖ=ë~ãéäÉK= cçê=~å=áåéìí=îÉÅíçê=~ää=é~ííÉêå=ä~óÉê=åÉìêçåë=ÅçãéìíÉ=íÜÉ=bìÅäáÇá~å=Çáëí~åÅÉ=ÄÉíïÉÉå=íÜÉ= áåéìí= îÉÅíçê= ~åÇ= íÜÉ= ÅçêêÉëéçåÇáåÖ= åÉìêçå= äçÅ~íáçåK= qÜÉëÉ= Çáëí~åÅÉë= ~êÉ= ~ëëìãÉÇ= ~ë= ~å= ~Åíáî~íáçå=ÑìåÅíáçå=E~=d~ìëëá~å=âÉêåÉä=áå=çìê=Å~ëÉF=~åÇ=íê~åëÑÉêêÉÇ=íç=íÜÉ=ÑçääçïáåÖ=ä~óÉêK= qÜÉ=ëìãã~íáçå=ä~óÉê=Åçåëáëíë=çÑ=íïç=åÉìêçåë=íÜ~í=Å~äÅìä~íÉ=íÜÉ=åçãáå~íçê=~åÇ=ÇÉåçãáJ å~íçê=êÉëéÉÅíáîÉäó=áå=Éèì~íáçå=NTK=b~ÅÜ=çÑ=íÜÉëÉ=åÉìêçåë=ÅçãéìíÉë=~=ïÉáÖÜíÉÇ=ëìã=çÑ=íÜÉ= çìíéìí=Ñêçã=~=éêÉîáçìë=ä~óÉêK=qÜÉ=ïÉáÖÜíë=ÅçêêÉëéçåÇ=íç=íÜÉ=äáåâë=ÄÉíïÉÉå=íÜÉ=åÉìêçåëK= b~ÅÜ= äáåâ= Ñêçã= íÜÉ= é~ííÉêå= ä~óÉê= åÉìêçå= íç= íÜÉ= åçãáå~íçê= åÉìêçå= áë= ~= í~êÖÉí= î~äìÉ==

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

19-5

~ëëçÅá~íÉÇ=ïáíÜ=íÜÉ=ÅçêêÉëéçåÇáåÖ=åÉìêçå= äçÅ~íáçåK=^ää=äáåâë=ÄÉíïÉÉå=é~ííÉêå=ä~óÉê=åÉìJ êçåë= ~åÇ= ~= ÇÉåçãáå~íçê= åÉìêçå= ~êÉ= Éèì~ä= íç= ~= ìåáíK= qÜÉ= çìíéìí= ä~óÉê= åÉìêçå= éÉêÑçêãë= Çáîáëáçå=çéÉê~íáçå=~åÇ=íê~åëÑÉêë=íÜÉ=çìíéìí=î~äìÉK=

Figure 1 NWKRE as a general regression neural network (GRNN)

qÜÉ= çåäó= åçåJåÉìê~ä= çéÉê~íáçå= EÇáîáëáçåF= áå= íÜáë= ~êÅÜáíÉÅíìêÉ= áë= íÜ~í= çÑ= íÜÉ= çìíéìí= ä~óÉêK= qÜÉêÉ= áë= ~= ÖÉåÉê~ä= EíÜçìÖÜ= åçí= ìåáîÉêë~äF= ~ÖêÉÉãÉåí= íÜ~í= ÄáçäçÖáÅ~ä= åÉìêçåë= ~êÉ= áåÅ~é~ÄäÉ= çÑ= ÇáîáëáçåK= qÜáë= áë= ~= éçáåí= çÑ= ÅêáíáÅáëã= ~Ö~áåëí= áåíÉêéêÉíáåÖ= íÜÉ= kthob= ~ë= ~= åÉìê~ä=åÉíïçêâK=píáääI=ëìÅÜ=~å=áåíÉêéêÉí~íáçå=áë=îÉêó=~ííê~ÅíáîÉ=ÄÉÅ~ìëÉ=çÑ=íÜÉ=ëáãéäÉ=ëÅÜÉJ ã~íáÅ= áääìëíê~íáçå= çÑ= íÜÉ= áåíÉêéçä~íáçå= éêçÅÉëë= EcáÖìêÉ= NF= ~åÇ= íÜÉ= éçëëáÄáäáíó= íç= ìëÉ= ã~ÅÜáåÉ=äÉ~êåáåÖ=íÉêãáåçäçÖóK= 2.3. TRAINING OF THE GRNN

qÜÉ=çåäó=~Ç~éíáîÉ=EÑêÉÉF=é~ê~ãÉíÉê=áå=íÜÉ=dokk=ãçÇÉä=ïáíÜ=d~ìëëá~å=âÉêåÉä=EÉèì~íáçå= NTF=áë= σI=íÜÉ=ïáÇíÜ=çÑ=íÜÉ=âÉêåÉäK=qç=ÇÉãçåëíê~íÉ=Üçï=~=dokk=ïçêâëI=äÉí=ìë=ÅçåëáÇÉê=~= ëáãéäÉ=çåÉJÇáãÉåëáçå~ä=éêçÄäÉã=EcáÖìêÉ=OFK=^=ëáãéäÉ=ëáåÉ=ÑìåÅíáçå=êÉéêÉëÉåíë=~=íêìÉ=ÑìåÅJ íáçå= çÑ= ~å= ìåÇÉêäóáåÖ= ëíêìÅíìêÉ= çÑ= ~= ÅçääÉÅíáçå= çÑ= ë~ãéäÉ= Ç~í~K= få= çêÇÉê= íç= ÖÉåÉê~íÉ= ~= íê~áåáåÖ=ëÉíI=íÜáë=ëáåÉ=ÑìåÅíáçå=áë=ë~ãéäÉÇ=~í=~=ïáÇÉ=î~êáÉíó=çÑ=éçáåíëI=~åÇ=ê~åÇçã=åçáëÉ=áë= ~ÇÇÉÇK=cìêíÜÉêãçêÉI=çåÉ=ïáäÇ=éçáåí=EçìíäáÉêF=áë=áåÅäìÇÉÇK=qÜáë=íêìÉ=ÑìåÅíáçå=~åÇ=íê~áåáåÖ= ë~ãéäÉë=~êÉ=Öê~éÜÉÇ=áå=cáÖìêÉ=O~K= kçï= äÉí= ìë= Éñ~ãáåÉ= íÜÉ= ÉÑÑÉÅí= éêçÇìÅÉÇ= Äó= ~= î~êáÉíó= çÑ= σ= î~äìÉëI= íÜÉ= ëãççíÜáåÖ== é~ê~ãÉíÉêK= få= cáÖìêÉ= OÄ= ïÉ= ëÉÉ= ïÜ~í= Ü~ééÉåë= ïÜÉå= ~= îÉêó= ëã~ää= î~äìÉ= áë= ìëÉÇK= qÜÉ= dokk=Ñçääçïë=íÜÉ=íê~áåáåÖ=Ç~í~=ÅäçëÉäóI=~äãçëí=ãçîáåÖ=Ñêçã=éçáåí=íç=éçáåíK=fÑ=íÜÉ=Ç~í~=áë= âåçïå=íç=ÄÉ=ÅäÉ~êI=íÜÉ=dokk=ã~âÉë=~å=ÉñÅÉääÉåí=áåíÉêéçä~íáçå=~äÖçêáíÜãI=~å~äçÖçìë=íç= íÜÉ=åÉ~êÉëí=åÉáÖÜÄçê=ãÉíÜçÇK=kçíÉI=ÜçïÉîÉêI=íÜ~í=íÜáë=êÉëìäí=ïáää=ÄÉ=~ÅÅÉéí~ÄäÉ=çåäó=áÑ=íÜÉ= ÇÉåëáíó=çÑ=íÜÉ=íê~áåáåÖ=Ç~í~=áë=ÜáÖÜ=ÉåçìÖÜK=få=çíÜÉê=Å~ëÉëI=~å=“çîÉêÑáííáåÖÒ=ÉÑÑÉÅíI=ïÜáÅÜ= áë= ïÉääJâåçïå= áå= åÉìê~ä= åÉíïçêâëI= ã~ó= ~ééÉ~ê= ~åÇ= ëìÅÜ= ëçäìíáçåë= ïáää= åçí= ÄÉ= ~= Ñ~áäìêÉK= pçI=ëáåÅÉ=ïÉ=âåçï=íÜ~í=íÜÉ=Ç~í~=~êÉ=ÇáëíçêíÉÇ=Äó=åçáëÉ=áå=ãçëí=Å~ëÉëI=ëíê~áÖÜíÑçêï~êÇ=áåJ íÉêéçä~íáçå=áë=åçí=ïÜ~í=ïÉ=ï~åíK= cáÖìêÉ=OÅ=áë=çÄí~áåÉÇ=Äó=ìëáåÖ=~=ãçÇÉê~íÉäó=ëã~ää=ëãççíÜáåÖ=é~ê~ãÉíÉêK=lÄëÉêîÉ=íÜ~í= íÜÉ=ïáäÇ=éçáåí=E~å=çìíäáÉêF=ëìÅÅÉÉÇë=áå=éìääáåÖ=íÜÉ=ÑìåÅíáçå=~=Äáí=çìí=çÑ=íÜÉ=äáåÉK=fÑ=ïÉ=~êÉ= ÅçåÑáÇÉåí=íÜ~í=~åó=ëìééçëÉÇ=ïáäÇ=éçáåí=áë=~Åíì~ääó=~=î~äáÇ=Ç~íìãI=íÜÉå=íÜáë=áë=àìëí=ïÜ~í=ïÉ= ï~åíK=_ìí=ìëì~ääó=áí=áë=ÇìÉ=íç=ëçãÉ=âáåÇ=çÑ=Ç~í~=ÇáëíçêíáçåK=^=ëçãÉïÜ~í=ä~êÖÉê=ëãççíÜáåÖ= é~ê~ãÉíÉê=ÖáîÉë=ìë=íÜÉ=êÉëìäí=ëÜçïå=áå=cáÖìêÉ=OÇI=ïÜáÅÜ=äççâë=áÇÉ~äK= cáå~ääóI=ïÉ=ãìëí=ÄÉ=~ï~êÉ=çÑ=í~âáåÖ=~=ÖççÇ=íÜáåÖ=íçç=Ñ~êK=cáÖìêÉ=OÉ=áääìëíê~íÉë=íÜÉ=ÉÑJ ÑÉÅí= çÑ= íçç= ä~êÖÉ= ~= ëãççíÜáåÖ= é~ê~ãÉíÉêK= qÜÉ= ÖäçÄ~ä= ëÜ~éÉ= çÑ= íÜÉ= íê~áåáåÖ= ëÉí= Ü~ë= ëï~ãéÉÇ=çìí=~ää=íÜÉ=ÇÉí~áäëK=dê~åíÉÇI=ïÉ=ìëì~ääó=ï~åí=íÜÉ=ÑáåÉ=ÇÉí~áäëI=ïÜáÅÜ=~êÉ=ÇìÉ=íç= åçáëÉI=íç=ÄÉ=Éäáãáå~íÉÇK=

19-6

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

qÜìë=ïÉ=Å~å=ã~âÉ=íÜÉ=ÅçåÅäìëáçå=íÜ~í=ÅÜççëáåÖ=íÜÉ=î~äìÉ=çÑ=íÜÉ=ïáÇíÜ=é~ê~ãÉíÉê= σ=áë= îáí~ä=Ñçê=íÜÉ=dokk=ãçÇÉä=~åÇ=éêçÄäÉãJÇÉéÉåÇÉåíK=

Figure 2 A simple problem to illustrate the influence of the parameter on the GRNN result: a) true function and noised training samples with an outlier; b) too small σ; c) perhaps a bit small σ; d) a perfect σ; e) too large σ

qÜÉ=çêÇáå~êó= dokk=íê~áåáåÖ=éêçÅÉÇìêÉ=áë=~=ãÉ~åJëèì~êÉJÉêêçê=EjpbF=ãáåáãáò~íáçå= éêçÅÉÇìêÉI= ~ÅÅçãéäáëÜÉÇ= ìëáåÖ= ~= ÅêçëëJî~äáÇ~íáçå= EâJÑçäÇI= äÉ~îÉJçåÉJçìíF= ~ééêç~ÅÜ= Ee~óâáå= NVVVX= _áëÜçé= NVVRFK= ^= í~êÖÉí= ÑìåÅíáçå= EjpbF= áë= ÅçãéìíÉÇ= Ñêçã= êÉëìäíë= Ñçê= ~ää= ë~ãéäÉë=í~âÉå=çìíI=çåÉ=Äó=çåÉI=~åÇ=Éëíáã~íÉÇ=Ä~ëÉÇ=çå=íÜÉ=çíÜÉê=ë~ãéäÉë=~åÇ=~=ÖáîÉå= σ= î~äìÉK=qÜÉ=î~äìÉë=çÑ=σ=ÅÜ~åÖÉI=í~âáåÖ=î~äìÉë=Ñêçã=~å=áåíÉêî~ä=xσlowI=σhighzW=

σ i = σ low + (i − 1)

σ high − σ low M

(18)

ïÜÉêÉ=j=áë=íÜÉ=åìãÄÉê=çÑ=ëíÉéë=íç=éÉêÑçêãK=qÜÉ=áåáíá~ä=áåíÉêî~ä=çÑ=éçëëáÄäÉ= σ=î~äìÉë=~åÇ= íÜÉ= åìãÄÉê= çÑ= ëíÉéë= ãìëí= ÄÉ= ÇÉÑáåÉÇK= qÜÉ= áåíÉêî~ä= ~åÇ= íÜÉ= åìãÄÉê= çÑ= ëíÉéë= Ü~îÉ= íç= ÄÉ= ÅçåëáëíÉåí= áå= çêÇÉê= íç= Å~íÅÜ= íÜÉ= ÉñéÉÅíÉÇ= çéíáã~ä= î~äìÉK= qç= ã~âÉ= íÜÉ= áåíÉêî~ä= ëÉäÉÅíáçå= ãçêÉ=ÖÉåÉê~äI=áí=áë=ìëÉÑìä=íç=íê~åëÑçêã=íÜÉ=áåéìí=ëé~ÅÉ=î~äìÉë=íç=íÜÉ=xJNI=Nz=áåíÉêî~ä=ìëáåÖI= Ñçê= Éñ~ãéäÉI= ~= ëáãéäÉ= äáåÉ~ê= íê~åëÑçêãK= ^ÑíÉê= íÜ~í= ~ää= éçëëáÄäÉ= ïáÇíÜë= ~êÉ= ïáíÜáå= ~= ìåáíK= qÜìëI=íÜÉ=ëí~êíáåÖ=áåíÉêî~ä=~åÇ=íÜÉ=åìãÄÉê=çÑ=ëíÉéë=~êÉ=åçí=îÉêó=áãéçêí~åíK=`ÜççëáåÖ=íçç= ïáÇÉ=~å=áåíÉêî~ä=~åÇLçê=~=ä~êÖÉ=åìãÄÉê=çÑ=ëíÉéë=çåäó=äÉ~Ç=íç=ëéÉåÇáåÖ=ãçêÉ=íáãÉ=çå=Å~äÅìJ ä~íáçåëK=_ìí=ïáíÜ=ãçÇÉêå=ÅçãéìíÉêë= ~åÇ=~=Ç~í~=ëÉí=ïáíÜ=~Äçìí=OMM=éçáåíë=íÜáë=Ñ~Åíçê=áë= åçí=îÉêó=ëáÖåáÑáÅ~åíK= ^= î~äìÉ= çÑ= σ= ïÜáÅÜ= ãáåáãáòÉë= íÜÉ= jpb= Å~å= ÄÉ= ~ÇçéíÉÇ= ~ë= ~= ëçäìíáçåK= cìêíÜÉêãçêÉI= ~åó=Öê~ÇáÉåí=ëÉ~êÅÜ=ãÉíÜçÇ=Å~å=ÄÉ=ìëÉÇ=Ñçê=íìåáåÖ=íÜáë=î~äìÉK=táíÜ=~å=ÉñÅÉääÉåí=ëí~êíáåÖ= éçáåí= Ó= ~= ïÉääJÇÉÑáåÉÇ= ãáåáãìãI= ~åÇ= ~= ëã~ää= ëÉ~êÅÜ= ~êÉ~= Ó= íÜÉ= Öê~ÇáÉåí= ëÉ~êÅÜ= ãÉíÜçÇ= ÑáåÇë=íÜÉ=ëçäìíáçå=èìáÅâäó=~åÇ=É~ëáäóK=bîÉå=áÑ=~å=çéíáã~ä=î~äìÉ=çÑ= σ=áë=çìíëáÇÉ=íÜÉ=áåáíá~ä= áåíÉêî~äI=áí=áë=åçí=~=éêçÄäÉã=Ñçê=íÜÉ=Öê~ÇáÉåí=ëÉ~êÅÜ=ãÉíÜçÇ=íç=Çç=ëçãÉ=ëíÉéë=çìíI=~ë=áí=áë= åçí=äáãáíÉÇ=Äó=~åó=ÄçìåÇ~êáÉëK= qÜÉ=pf`OMMQ=ÉñÉêÅáëÉ=ÇÉ~äë=ïáíÜ=íïçJÇáãÉåëáçå~ä=Ç~í~=áå=íÜÉ=EñIóF=ÅççêÇáå~íÉ=ëé~ÅÉK= få= íÜáë= Å~ëÉI= ~å= ~Çî~åÅÉÇ= dokk= êÉ~äáò~íáçå= ã~ó= ÄÉ= ~ééäáÉÇI= ìëáåÖ= ~å= ~åáëçíêçéáÅ= Oa= é~ê~ãÉíÉê= σZEσxIσyFK= `çåëÉèìÉåíäóI= ~= Oa= EêÉÅí~åÖìä~êF= áåíÉêî~ä= áë= ìëÉÇ= Ñçê= σ= çéíáãáòáåÖ= íÜÉ=~äÖçêáíÜã=~å~äçÖçìë=íç=çåÉ=ÇÉëÅêáÄÉÇ=~ÄçîÉ=Ñçê=~å=áëçíêçéáÅ=Å~ëÉK=låäó=íÜÉêÉ=ïçìäÇ=

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

19-7

ÄÉ=jxGjó=î~äìÉë=Ñçê= σ=íç=ÅÜÉÅâK=rëì~ääó=íÜáë=Oa=Éêêçê=ÑìåÅíáçå=Ü~ë=çåÉ=ëãççíÜÉÇI=ïÉääJ ÇÉÑáåÉÇ=ãáåáãìãK=qÜÉ=î~äìÉ=çÑ=σ=ÖáîáåÖ=íÜÉ=ãáåáãìã=Éêêçê=ã~ó=ÄÉ=~ÇçéíÉÇ=~ë=~=ëçäìíáçå= çê=íìåÉÇ=ìëáåÖ=~=Öê~ÇáÉåí=ëÉ~êÅÜ=ãÉíÜçÇK=få=íÜÉ=ÅìêêÉåí=ëíìÇó=~=ÅçåàìÖ~íÉ=Öê~ÇáÉåí=ëÉ~êÅÜ= ãÉíÜçÇ=ï~ë=ìëÉÇ=Ñçê=Ñáå~ä=íìåáåÖ=çÑ=σK=qÜáë=~äÖçêáíÜã=áë=~ÄëçäìíÉäó=~ìíçã~íáÅI=åçí=êÉèìáêJ áåÖ= ~åó= çéíáãáò~íáçå= é~ê~ãÉíÉêë= ~åÇ= ~äëç= îÉêó= Ñ~ëí= ~åÇ= ÉÑÑÉÅíáîÉK= fíë= ëçìêÅÉ= ÅçÇÉ= ~åÇ= ~ééäáÅ~íáçå=ÇÉí~áäë=Å~å=ÄÉ=ÑçìåÇ=áå=j~ëíÉêë=ENVVRFK=^ééäóáåÖ=~=Öê~ÇáÉåí=ëÉ~êÅÜ=Ñçê=Ñáå~ä= é~ê~ãÉíÉê=íìåáåÖ=~ääçïÉÇ=ìë=íç=Å~íÅÜ=~=î~äìÉ=Ñçê= σ=ÖççÇ=ÉåçìÖÜ=Ñçê=íÜÉ=ëÉÅçåÇ=E“àçâÉêÒF= Ç~í~=ëÉíI=ïÜáÅÜ=~ééÉ~êÉÇ=íç=Ñ~ää=çìíëáÇÉ=íÜÉ=áåáíá~ä=áåíÉêî~äK= cáå~ääóI= íÜÉ= éêçÅÉÇìêÉ= Ñçê= ~ééäóáåÖ= dokk= Ñçê= íÜÉ= pf`OMMQ= ÉñÉêÅáëÉ= ï~ë= íÜÉ== ÑçääçïáåÖW= NK kçêã~äáòÉ=íÜÉ=áåéìí=ëé~ÅÉ=íç=íÜÉ=áåíÉêî~ä=xJNI=Nz=EÄó=äáåÉ~ê=íê~åëÑçêã~íáçåFK= OK aÉÑáåÉ=σlow=~ë=MKMP=~åÇ=σhigh=~ë=MKQ=Ñçê=ÄçíÜ=ÇáêÉÅíáçåë=Eáå=íê~åëÑçêãÉÇ=ÅççêÇáå~íÉëFK= qÜÉ= ÄçìåÇë= ïÉêÉ= ëÉäÉÅíÉÇ= Ä~ëÉÇ= çå= ~å~äóëáë= çÑ= íÜÉ= éêáçê= Ç~í~K= qÜÉ= åìãÄÉê= çÑ== Å~äÅìä~íáçåë=áë=NM=Ñçê=ÄçíÜ=ÇáêÉÅíáçåëK=pç=ïÉ=Ü~îÉ=NMGNMZNMM=ÇáÑÑÉêÉåí=î~äìÉë=çÑ=σK= PK açïåäç~Ç=Ç~í~=ëÉíë=Ñêçã=íÜÉ=ïÉÄ=ëáíÉK=pí~êí=íÜÉ=ÅçãéÉíáíáçåK= QK pí~êí=íÜÉ=äÉ~îÉJçåÉJçìí=ÅêçëëJî~äáÇ~íáçå=êçìíáåÉ=Ñçê=É~ÅÜ=NMM=î~äìÉë=çÑ= σK=`~äÅìä~íÉ= íÜÉ=ÅêçëëJî~äáÇ~íáçå=Éêêçê=Ñçê=É~ÅÜ=î~äìÉK=cáÖìêÉ=P=éêÉëÉåíë=ÅêçëëJî~äáÇ~íáçå=Éêêçêë=~ë= ~=ÑìåÅíáçå=çÑ=σK=qÜÉ=äçïÉê=áë=íÜÉ=jpbI=íÜÉ=äáÖÜíÉê=áë=íÜÉ=ÅçäçêK=qÜÉ=ëÅ~äÉ=Ñçê=Éêêçêë=áë= åçí=áãéçêí~åí=ÄÉÅ~ìëÉ=ïÉ=ëÉ~êÅÜ=Ñçê=ãáåáãìã=çåäó=~åÇ=Çç=åçí=é~ó=~åó=~ííÉåíáçå=íç= íÜÉ=î~äìÉë=íÜÉãëÉäîÉëK= RK pÉäÉÅí=íÜÉ= σ=î~äìÉ=íÜ~í=ãáåáãáòÉë=íÜÉ=ÅêçëëJî~äáÇ~íáçå=ÉêêçêK=rëÉ=íÜáë=î~äìÉ=~ë=~=ëí~êíJ áåÖ=éçáåí=Ñçê=íÜÉ=ÅçåàìÖ~íÉÇ=Öê~ÇáÉåí=íìåáåÖ=éêçÅÉÇìêÉK= SK rëáåÖ= íÜÉ= íê~åëÑçêã~íáçå= é~ê~ãÉíÉêë= Ñêçã= ëíÉé= NI= Çç= Ä~Åâ= íê~åëÑçêã~íáçå= Ñçê= íÜÉ= çéíáã~ä= σ=î~äìÉK=q~ÄäÉ=N=éêÉëÉåí=íÜÉ=ÅçêêÉëéçåÇáåÖ=êÉëìäíëK=`çãé~êÉ=ïáíÜ=q~ÄäÉ=P= Ñçê=êÉëìäíë=çÑ=σ=íìåáåÖ=Ñçê=éêáçê=Ç~í~K= TK mÉêÑçêã=éêÉÇáÅíáçå=Ñçê=UMU=éçáåíë=ìëáåÖ=íÜÉ=çéíáã~ä=σ=î~äìÉK= UK pÉåÇ=êÉëìäíë=íç=íÜÉ=ÉñÉêÅáëÉ=çêÖ~åáòÉêëK=båÇ=çÑ=íÜÉ=ÅçãéÉíáíáçåK=

Figure 3 Contour lines of the cross-validation errors (MSE) for the 1st (left) and the 2nd data set (right). Crosses indicate locations of the values from the initial interval

σ

σx

Data sets First data set Second data set

Original 27 287 21 861

σy Transformed 0.16 0.10

Original 34 285 6 959

Transformed 0.13 0.02

Table 1 Optimal anisotropic σ=(σx,σy) values after step 6

19-8

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

2.4. USE OF PRIOR INFORMATION

qÜÉ= dokk= íê~áåáåÖ= éêçÅÉÇìêÉ= ïáíÜ= ~å= áåáíá~ä= áåíÉêî~ä= çÑ= xMKMNI= MKSz= ï~ë= ~ééäáÉÇ= íç= ~ää= éêáçê=Ç~í~=ëÉíëK=q~ÄäÉ=O=éêÉëÉåíë=íÜÉ=êççíJãÉ~åJëèì~êÉ=Éêêçê=EojpbF=~åÇ=mÉ~êëçåÛë=ÅçêêÉJ ä~íáçå= ÅçÉÑÑáÅáÉåí= ê= Ñçê= ~ää= NM= Ç~í~= ëÉíëK= léíáã~ä= ~åáëçíêçéáÅ= σ= î~äìÉë= ~êÉ= éêÉëÉåíÉÇ= áå= q~ÄäÉ=PK= Prior day

Std. Dev.

1

2

3

4

5

6

7

8

9

10

Mean

RMSE

8.85

8.87

8.41

8.52

8.37

7.99

8.34

8.39

8.72

8.70

8.52

0.26

r

0.89

0.89

0.89

0.87

0.87

0.87

0.86

0.87

0.88

0.87

0.88

0.01

Table 2 Training errors for all prior data sets

Prior

Std. Dev.

1

2

3

4

5

6

7

8

9

10

Mean

σx

0.16

0.16

0.16

0.15

0.17

0.17

0.20

0.16

0.16

0.14

0.16

0.01

σy

0.09

0.09

0.09

0.10

0.10

0.09

0.08

0.10

0.10

0.11

0.10

0.01

Table 3 Optimal σ values for all prior data sets (in transformed coordinates)

qÜÉ=ã~áå=éìêéçëÉ=çÑ=íÜÉ=éêáçê=Ç~í~=~å~äóëáë=ï~ë=íç=ÇÉÑáåÉ=íÜÉ=áåáíá~ä=áåíÉêî~ä=Ñçê=íÜÉ= çéíáãáë~íáçå= éêçÅÉÇìêÉK= qÜÉ= áåíÉêî~ä= ï~ë= ëÉäÉÅíÉÇ= ~ë= xMKMPI= MKQz= Eáå= íÜÉ= íê~åëÑçêãÉÇ= ÅçJ çêÇáå~íÉëFI= íÜÉ= ë~ãÉ= Ñçê= ÄçíÜ= ÇáêÉÅíáçåëK= qÜÉ= ÉñéÉÅíÉÇ= çéíáã~ä= âÉêåÉä= ïáÇíÜ= ï~ë= ~Äçìí= MKNS=Ñçê=σñ=~åÇ=~Äçìí=MKN=Ñçê=σóK=fí=ï~ë=ÅçãéäÉíÉäó=ÅçåÑáêãÉÇ=Ñçê=íÜÉ=Ñáêëí=Ç~í~=ëÉíI=ëáãáä~ê= áå=ÑÉ~íìêÉë=íç=íÜÉ=éêáçê=Ç~í~=ëÉíëK= cçê=íÜÉ=ëÉÅçåÇ=Ç~í~=ëÉí=íÜÉ=çéíáã~ä= σ=î~äìÉë=ÇáÑÑÉê=ëáÖåáÑáÅ~åíäóW=MKNP=~åÇ=MKMOK=fí=áë=~= ÖççÇ= áåÇáÅ~íçê= íÜ~í= íÜÉêÉ= áë= “ëçãÉíÜáåÖ= ïêçåÖÒ= ïáíÜ= íÜáë= Ç~í~= ëÉíK= _ìí= ~ìíçã~íáÅ= ÅêçëëJ î~äáÇ~íáçå= ÅçãÄáåÉÇ= ïáíÜ= ~= Öê~ÇáÉåí= íìåáåÖ= éêçÅÉÇìêÉ= ~ääçïÉÇ= ìë= íç= ÅÜççëÉ= ÖççÇ== é~ê~ãÉíÉêë=ÉîÉå=áÑ=íÜÉó=~ééÉ~êÉÇ=íç=ÄÉ=çìíëáÇÉ=íÜÉ=áåáíá~ä=áåíÉêî~äK=lìê=ÖìÉëë=áë=íÜ~í=é~êíáÅáJ é~åíë=ïÜç=ìëÉÇ=çåäó=íÜÉ=éêáçê=áåÑçêã~íáçå=Ñçê=íê~áåáåÖ=~åÇ=ÇáÇ=åçí=éÉêÑçêã=~åó=~ìíçã~íáÅ= íìåáåÖ=~äÖçêáíÜãë=ïçìäÇ=åçí=çÄí~áå=îÉêó=ÖççÇ=êÉëìäíë=Ñêçã=íÜÉ=“àçâÉêÒ=Ç~í~=ëÉíK= 2.5. MEASUREMENT OF THE UNCERTAINTY

råÇÉê= íÜÉ= ~ëëìãéíáçå= çÑ= áåÇÉéÉåÇÉåíI= áÇÉåíáÅ~ääó= ÇáëíêáÄìíÉÇ= EáKáKÇKF= ê~åÇçã= î~êá~ÄäÉë= çîÉê=íÜÉ=ïÜçäÉ=áåéìí=ëé~ÅÉ=Eáå=çìê=Å~ëÉ=~=åçêã~ä=Ed~ìëëá~åF=ÇáëíêáÄìíáçåFI=~=dokk=ÉëíáJ ã~íÉ=EwGF=Å~å=ÄÉ=íêÉ~íÉÇ=~ë=~=ãÉ~å=çÑ=~=ê~åÇçã=î~êá~ÄäÉK=s~êá~íáçå=Å~å=ÄÉ=ÅçåëáÇÉêÉÇ=~ë= O íÜÉ=ëèì~êÉÇ=çÑ=íÜÉ=êÉëáÇì~äë=Ó=EwGJwF K=qÜìëI=ìåÅÉêí~áåíó=Eî~êá~íáçåF=áë=ÅçãéìíÉÇ=Ñêçã=íÜÉ= êÉëáÇì~äë=~í=íÜÉ=íê~áåáåÖ=äçÅ~íáçåë=EïÜÉêÉ=íÜÉ=êÉ~ä=î~äìÉë=~êÉ=âåçïåF=Äó=áåíÉêéçä~íáåÖ=çîÉê= íÜÉ=ïÜçäÉ=ëé~ÅÉK= fåíÉêéçä~íáçå= Å~å= ÄÉ= éÉêÑçêãÉÇ= ïáíÜ= íÜÉ= ÜÉäé= çÑ= ~å= ~äêÉ~Çó= íê~áåÉÇ= dokk= Eσ= = é~J ê~ãÉíÉê=íìåÉÇ=Ñçê=Ç~í~=éêÉÇáÅíáçåFK=fí=áë=åÉÅÉëë~êó=íç=âÉÉé=áå=ãáåÇ=íÜ~í=ëìÅÜ=~å=Éëíáã~íáçå= çÑ= î~êá~íáçåëI= ÑáêëíI= ÇÉéÉåÇë= çå= íÜÉ= Ç~í~= éçáåíëÛ= ÇáëíêáÄìíáçå= EÇÉåëáíóFI= ~åÇI= ëÉÅçåÇI= áë== Äá~ëÉÇ=E~ë=Ä~ëÉÇ=çå=íÜÉ=Äá~ëÉÇ=ãÉ~å=Éëíáã~íáçåFK=qÜÉ=áKáKÇK=~ëëìãéíáçå=áë=íçç=ëíêçåÖI=~åÇ= áí=~ääçïë=ìë=íç=çÄí~áå=çåäó=êìÇÉI=~ééêçñáã~íÉ=Éëíáã~íÉë=çÑ=íÜÉ=êÉ~ä=ìåÅÉêí~áåíóK=få=êÉ~äáíó= áí=çåäó=~ääçïë=ìë=íç=ÇÉíÉÅí=~êÉ~ë=ïÜÉêÉ=ìåÉñéÉÅíÉÇ=EÑ~ê=Ñêçã=éêÉÇáÅíÉÇF=î~äìÉë=ã~ó=çÅÅìê= ïáíÜ=éçëáíáîÉ=éêçÄ~ÄáäáíóK=

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

19-9

qÜÉ= çíÜÉê= éçëëáÄäÉ= ~ééêç~ÅÜ= áë= íç= ÄìáäÇ= ~= åÉï= ëéÉÅá~ä= dokk= íç= éêÉÇáÅí= íÜÉ= ìåÅÉêJ í~áåíóI=~ë=ï~ë=ÇçåÉ=áå=íÜÉ=ÅìêêÉåí=Å~ëÉK=qÜÉ=éêçÅÉëë=áë=Ñìääó=~ìíçã~íáÅ=E~ë=ïÉää=~ë=íÜÉ=ã~áå= êçìíáåÉFI= ~åÇ= áí= ï~ë= áãéäÉãÉåíÉÇ= ~ë= é~êí= çÑ= íÜÉ= ~ìíçã~íáÅ= ÇÉÅáëáçå= ã~âáåÖ= ëóëíÉã= íçJ ÖÉíÜÉê=ïáíÜ=~=dokk=áíëÉäÑK=fí=çåäó=~ÇÇë=~å=~ÇÇáíáçå~ä=íê~áåáåÖ=éêçÅÉÇìêÉ=ïÜÉêÉ=ëèì~êÉÇ= êÉëáÇì~äë=~í=íÜÉ=íê~áåáåÖ=éçáåíë=~êÉ=ìëÉÇ=~ë=~=í~êÖÉí=ÑìåÅíáçåK=qÜÉ=çéíáã~ä=î~äìÉ=çÑ= σ=Ñçê= íÜáë=ìåÅÉêí~áåíó=dokk=áë=ÑçìåÇ=Äó=íÜÉ=ë~ãÉ=éêçÅÉÇìêÉ=~ë=ÇÉëÅêáÄÉÇ=~ÄçîÉK=qÜÉ=íê~áåÉÇ= ãçÇÉä=Å~å=ÄÉ=ìëÉÇ=íç=Éëíáã~íÉ=íÜÉ=ìåÅÉêí~áåíó=~í=íÜÉ=UMU=ìåâåçïå=äçÅ~íáçåë=çê=~åóïÜÉêÉ= ÉäëÉK=

3. RESULTS 3.1. OVERALL RESULTS

q~ÄäÉ= Q= éêÉëÉåíë= íÜÉ= Ä~ëáÅ= ëí~íáëíáÅë= EãáåáãìãI= ã~ñáãìãI= ãÉ~åI= ãÉÇá~å= ~åÇ= ëí~åÇ~êÇ= ÇÉîá~íáçåF=çÑ=íÜÉ=UMU=Éëíáã~íÉÇ=~åÇ=çÄëÉêîÉÇ=î~äìÉë=Ñçê=íÜÉ=íïç=Ç~í~=ëÉíëK= N = 808 Observed (first data set) Estimates (first data set) Observed (second data set) Estimates (second data set)

Min. 57.0 69.9 57.0 58.7

Max. 180.0 125.9 1528.2 1330.5

Mean 98.0 96.8 105.4 104.9

Median 98.8 98.8 99.0 100.3

Std. Dev. 20.0 13.9 83.7 71.5

Table 4 Comparison of the estimated and measured values (nSv/h)

qÜÉ=ãÉ~å=~ÄëçäìíÉ=Éêêçê=Ej^bFI=íÜÉ=Äá~ë=Eçê=ãÉ~å=ÉêêçêI=jbFI=íÜÉ=êççíJãÉ~åJëèì~êÉÇ= Éêêçê= EojpbFI= mÉ~êëçåÛë= ÅçêêÉä~íáçå= ÅçÉÑÑáÅáÉåí= ê= ÄÉíïÉÉå= íÜÉ= Éëíáã~íÉÇ= ~åÇ= íêìÉ= î~äìÉë= ~åÇ= åçêã~äáòÉÇ= ojpb= EkojpbF= Ñçê= íÜÉ= éêÉÇáÅíáçåë= ~í= íÜÉ= å= Z= UMU= äçÅ~íáçåë= ~êÉ= éêÉJ ëÉåíÉÇ=áå=q~ÄäÉ=RK=kojpb=áë=ÅçãéìíÉÇ=Ñêçã=

NRMSE=RMSE / Std.Dev. ïÜÉêÉ=píÇKaÉîK=áë=íÜÉ=ëí~åÇ~êÇ=ÇÉîá~íáçå=çÑ=íÜÉ=çÄëÉêîÉÇ=î~äìÉëK=qÜÉ=kojpb=ëÜçïë=íÜÉ= êÉä~íáçå=ÄÉíïÉÉå=íÜÉ=ojpb=çÑ=íÜÉ=Éëíáã~íÉÇ=î~äìÉë=~åÇ=íÜÉ=ëí~åÇ~êÇ=ÇÉîá~íáçå=çÑ=íÜÉ=çÄJ ëÉêîÉÇ=î~äìÉëK=qÜìë=çåäó=î~äìÉë=çÑ=kojpb=äÉëë=íÜ~å=N=ã~ó=ÄÉ=~ÅÅÉéíÉÇ=~ë=ë~íáëÑ~ÅíçêóK= líÜÉêïáëÉI= áÑ= íÜÉ= éêÉÇáÅíçê= éêÉëÉåíë= íÜÉ= ãÉ~å= çÑ= íÜÉ= çÄëÉêîÉÇ= î~äìÉë= ~ë= íÜÉ= êÉëìäíI= íÜÉå= kojpbZNKM=áë=çÄí~áåÉÇK= Data sets

MAE

First data set Second data set

ME

RMSE

Pearson’s r

NRMSE

9.40

– 1.25

12.59

0.78

0.63

14.85

– 0.51

45.46

0.84

0.54

Table 5 Comparison of errors

qÜÉ=Ñáêëí=íÜêÉÉ=Éêêçê=ÅÜ~ê~ÅíÉêáëíáÅë=Ej^bI=jbI=ojpbF=ÇÉéÉåÇ=çå=î~äìÉë=çÑ=íÜÉ=Ç~í~= ëÉíë=éêÉëÉåíÉÇK=pç=áí=áë=ÇáÑÑáÅìäí=íç=Åçãé~êÉ=íÜÉ=èì~äáíó=çÑ=éêÉÇáÅíáçå=Ñçê=ÇáÑÑÉêÉåí=Ç~í~=ëÉíë= ïáíÜ=íÜÉëÉ=ÉêêçêëK=^åçíÜÉê=íïç=Éêêçêë=EêI=kojpbF=~êÉ=êÉä~íáîÉ=~åÇ=Çç=åçí=ÇÉéÉåÇ=çå=íÜÉ= ê~åÖÉ=çÑ=íÜÉ=Ç~í~K=s~äìÉë=çÑ=j^b=~åÇ=ojpb=Ñçê=íÜÉ=ëÉÅçåÇ=Ç~í~=ëÉí=~êÉ=ïçêëÉI=Äìí=íÜÉ= ÅçêêÉä~íáçå= ÅçÉÑÑáÅáÉåí= ~åÇ= kojpb= ~êÉ= ãìÅÜ= ÄÉííÉêK= qÜìëI= íÜÉ= ëÉÅçåÇ= Ç~í~= ëÉí= áë= éêÉJ ÇáÅíÉÇ=ÄÉííÉê=íÜ~å=íÜÉ=Ñáêëí=çåÉK=kÉÖ~íáîÉ=î~äìÉë=çÑ=íÜÉ=jb=áåÇáÅ~íÉ=íÜ~í=dokkI=çå=íÜÉ=

19-10

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

~îÉê~ÖÉI= ãçêÉ= íÉåÇë= íç= êÉÇìÅÉ= EìåÇÉêÉëíáã~íÉF= ÜáÖÜ= î~äìÉë= íÜ~å= áåÅêÉ~ëÉ= EçîÉêÉëíáã~íÉF= äçï=çåÉëK=fí=áë=íÜÉ=Å~ëÉ=Ñçê=ÄçíÜ=Ç~í~=ëÉíëK= qÜêÉÉ=ëÉíë=çÑ=ã~éë=ïáíÜ=áëçäáåÉë=~êÉ=éêÉëÉåíÉÇW=íÜÉó=ëÜçï=íÜÉ=Éëíáã~íÉÇ=î~äìÉë=EcáÖìêÉ= QF=~åÇ=íÜÉ=~ëëçÅá~íÉÇ=ìåÅÉêí~áåíáÉë=EcáÖìêÉë=R=~åÇ=SFK=cáÖìêÉ=S=éêÉëÉåíë=íÜÉ=ìåÅÉêí~áåíáÉë= Ñêçã=íÜÉ=Ää~Åâ=ëéçí=çÑ=cáÖìêÉ=R=EêáÖÜíFI=ïÜáÅÜ=~êÉ=åçí=ÇáëíáåÖìáëÜ~ÄäÉ=áå=íÜÉ=çêáÖáå~ä=ëÅ~äÉK= få=~ää=íÜÉëÉ=ÑáÖìêÉë=ÅêçëëÉë=áåÇáÅ~íÉ=éçáåí=äçÅ~íáçåë=çÑ=íÜÉ=Éëíáã~íÉÇ=î~äìÉë=~åÇ=ÄçñÉë=ÇÉJ éáÅí=äçÅ~íáçåë=çÑ=íÜÉ=áåéìí=î~äìÉëK=

Figure 4 Isoline levels (nSv/h) for the 1st set (left) and the 2nd set (right)

Figure 5 Isoline levels showing the uncertainty associated to the estimations obtained st nd for the 1 set (left) and the 2 set (right) ) in the same colour scale for comparison

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

19-11

Figure 6 Isoline levels showing the uncertainty for the black spot of Figure 5 (right)

3.2. DETECTING ANOMALIES AND OUTLIERS

Figure 7 3D map showing extreme values found in the 2nd set by GRNN åÇ

cáÖìêÉ=T=éêÉëÉåíë=íÜÉ=éÉ~â=ÑçìåÇ=áå=íÜÉ=O =Ç~í~=ëÉí=Äó=íÜÉ=dokk=~ééêç~ÅÜK=qÜÉ=ëãççíÜJ áåÖ=ÅÜ~ê~ÅíÉêáëíáÅ=Ñçê=dokk=äçïÉêÉÇ=íÜÉ=ÜÉáÖÜí=çÑ=íÜÉ=éÉ~âK=_ìí=îáëì~ääó=áí=ÇçÉë=åçí=äççâ= ëãççíÜÉÇ= Ó= ~= ÅçääÉÅíáçå= çÑ= ëÉîÉê~ä= ê~íÜÉê= ëÜ~êé= íçéëK= ^ìíçã~íáÅ= dokk= ~ääçïë= ìë= íç= ~Ç~éí=èìáÅâäó=íç=ÅçãéäÉíÉäó=åÉï=Ç~í~=~åÇ=íç=ã~âÉ=îÉêó=ÖççÇ=éêÉÇáÅíáçåëK=

4. DISCUSSION ^= dÉåÉê~ä= oÉÖêÉëëáçå= kÉìê~ä= kÉíïçêâ= EdokkF= áë= ~å= áãéäÉãÉåí~íáçå= çÑ= ~= åçåJ é~ê~ãÉíêáÅ=âÉêåÉä=Éëíáã~íçêK=få=íÜÉ=ÅìêêÉåí=ïçêâ=áí=áë=íêÉ~íÉÇ=~ë=~å=~ìíçã~íáÅ=ëé~íá~ä=éêÉJ ÇáÅíçê= çÑ= ê~Çáç~ÅíáîáíóK= fí= ï~ë= ëìÅÅÉëëÑìääó= ~ééäáÉÇ= íç= íïç= ÇáÑÑÉêÉåí= Ç~í~= ëÉíë= Ó= ìëì~ä= êçìíáåÉ=Ç~í~=~åÇ=ìåìëì~ä=ÉãÉêÖÉåÅó=Ç~í~K=dokk=áë=îÉêó=Ñ~ëí=Ó=ê~ï=íê~áåáåÖ=~åÇ=éêÉÇáÅJ íáçå=Å~äÅìä~íáçåë=íççâ=~Äçìí=PM=ëÉÅçåÇë=Ñçê=É~ÅÜ=Ç~í~=ëÉí=çå=~=mÉåíáìã=fff=RMM=jeò=m`K=

19-12

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

qÜÉ=ãçëí=áãéçêí~åí=éêçÄäÉã=ïáíÜ=dokk=éêÉÇáÅíçê=áë=ëãççíÜáåÖK=qÜÉ=ëãççíÜáåÖ=ÇÉJ ÖêÉÉ= ÇÉéÉåÇë= çå= dokk= é~ê~ãÉíÉêë= íìåÉÇ= Ñçê= íÜÉ= Ç~í~K= qÜÉ= éêÉëÉåÅÉ= çÑ= ëãççíÜáåÖ= Ñçê= íÜÉ=ÅìêêÉåí=í~ëâë=Å~å=ÄÉ=ëÉÉå=áå=q~ÄäÉ=QK=jáåáã~=~êÉ=çîÉêÉëíáã~íÉÇI=ã~ñáã~=~êÉ=ìåÇÉêÉëJ íáã~íÉÇI=~åÇ=ëí~åÇ~êÇ=ÇÉîá~íáçåë=Ñçê=íÜÉ=éêÉÇáÅíÉÇ=ÇáëíêáÄìíáçåë=~êÉ=ãìÅÜ=äçïÉê=íÜ~å=íÜÉ= çåÉë=çÑ=íÜÉ=êÉ~ä=ÑáÉäÇëK=_ìí=íÜÉ=ãÉ~åë=~åÇ=ãÉÇá~åë=~êÉ=êÉéêçÇìÅÉÇ=îÉêó=ïÉääK=fí=ãÉ~åë=íÜ~í= dokk=éêçîáÇÉÇ=~=ÖäçÄ~ääó=ìåÄá~ëÉÇI=ëãççíÜÉÇ=Éëíáã~íÉK= qÜÉ=ã~áå=éêçÄäÉãë=ïáíÜ=~ééäóáåÖ=dokk=~êÉ=ïÜÉíÜÉê=ëìÅÜ=~=ÇÉÖêÉÉ=çÑ=ëãççíÜáåÖ=áë= ë~íáëÑ~Åíçêó=~åÇ=ïÜÉíÜÉê=íÜÉ=êÉëìäíë=Ñêçã=íÜÉ=ãçÇÉäáåÖ=~êÉ=~ÅÅÉéí~ÄäÉ=Ñçê=ÇÉÅáëáçå=ã~âJ áåÖK=fí=ÉîáÇÉåíäó=ÇÉéÉåÇë=çå=íÜÉ=ã~áå=Öç~ä=çÑ=íÜÉ=ëíìÇóW=áÑ=ïÉ=~êÉ=áåíÉêÉëíÉÇ=áå=éêÉÇáÅíáåÖ= ÉñíêÉãÉ=î~äìÉë=éêÉÅáëÉäóI=íÜÉå=dokk=áë=åçí=ïÜ~í=ïÉ=åÉÉÇK=_ìí=áÑ=ïÉ=åÉÉÇ=~=Ñ~ëíI=É~ëáäó= íê~áåÉÇ=íççä=Ñçê=êÉÅçêÇáåÖ=~=Ñäìñ=~ë=~å=ÉîÉåíI=íÜÉå=dokk=áë=ÖççÇ=ÉåçìÖÜK=dokk=ÇÉÑáåÉÇ= íÜÉ=éêÉëÉåÅÉ=çÑ=~=Üçí=ëéçí=~åÇ=éçëáíáçåÉÇ=áí=ïÉääK=qç=ÄÉ=ãçêÉ=ÅÉêí~áå=~Äçìí=áíë=~ééäáÅ~ÄáäJ áíóI=ïÉ=ëÜçìäÇ=íÉëí=íÜÉ=ãçÇÉä=çå=çíÜÉê=í~ëâë=áåîçäîáåÖ=ÉãÉêÖÉåÅó=ëáíì~íáçåëK= iÉí= ìë= êÉã~êâ= ÜÉêÉ= íÜ~í= áí= áë= åçí= ~= éêçÄäÉã= íç= ÄìáäÇ= ~= dokk= EÇÉÑáåÉ= ~= σF= ïáíÜ= ~ää= íê~áåáåÖ=Éêêçêë=ÄÉáåÖ= òÉêç=~åÇ=~=ÅçêêÉä~íáçå=ÅçÉÑÑáÅáÉåí=çÑ=NK=fí=áë= ÉåçìÖÜ=íç=ìëÉσ=î~äìÉë= ÅäçëÉ=íç=òÉêçX=~=ëáãáä~ê=Éñ~ãéäÉ=ï~ë=ÇáëÅìëëÉÇ=Ñçê=~=ëáãéäÉ=éêçÄäÉã=EcáÖìêÉ=OÄFK=_ìí=ãçëí= éêçÄ~Ääó=áí=ïçìäÇ=áãéäó=~å=çîÉêÑáííáåÖI=ÄÉÅ~ìëÉ=ëìÅÜ=~=dokk=ÅçìäÇ=åçí=ÄÉ=ÖÉåÉê~äáòÉÇI= áKÉK=ìëÉÇ=íç=éêÉÇáÅí=î~äìÉë=Ñçê=åÉïI=ìåâåçïå=éçáåíëK=qÜÉ=ÖÉåÉê~äáò~íáçå=~Äáäáíó=áë=íÜÉ=ãçëí= áãéçêí~åí=ÑÉ~íìêÉ=çÑ=~åó=éêÉÇáÅíáåÖ=ãçÇÉäK= _ìí=~=ëã~ää=σ=ÇçÉë=åçí=~äï~óë=áãéäó=çîÉêÑáííáåÖK=cçê=íÜÉ=ëÉÅçåÇ=Ç~í~=ëÉí=íÜÉ=çéíáã~ä= σ= áë=ãìÅÜ=ëã~ääÉê=íÜ~å=Ñçê=íÜÉ=Ñáêëí=Ç~í~=ëÉí=E~åÇ=Ñçê=~ää=éêáçê=Ç~í~=ëÉíëFK=_ìí=ëíáää=íÜÉ=dokk= éêÉëÉåíë=îÉêó=ÖççÇ=ÖÉåÉê~äáò~íáçå=~ÄáäáíáÉëI=ïÜáÅÜ=Å~å=ÄÉ=çÄëÉêîÉÇ=íÜêçìÖÜ=~=ÜáÖÜ=ÅçêêÉä~J íáçå= ÅçÉÑÑáÅáÉåí= ÄÉíïÉÉå= êÉ~ä= ~åÇ= Éëíáã~íÉÇ= î~äìÉë= Eq~ÄäÉ= RFK= ^äëç= íÜÉ= Ä~ëáÅ= ëí~íáëíáÅë= Eq~ÄäÉ=QF=ëÜçïë=ÉîÉå=ÄÉííÉê=êÉéêçÇìÅíáçå=çÑ=ÖäçÄ~ä=ÇáëíêáÄìíáçå=íÜ~å=Ñçê=íÜÉ=Ñáêëí=Ç~í~=ëÉíK= qÜÉêÉ=áë=ëçãÉ=íêáÅâ=ïáíÜ=ìåÅÉêí~áåíó=Éëíáã~íáçå=ïÜáäÉ=ïçêâáåÖ=ïáíÜ=~=dokk=Ó=ìåÅÉêJ í~áåíó= áë= Éëíáã~íÉÇ= ~ë= ~= êÉÖêÉëëáçå= çÑ= âåçïå= ëèì~êÉÇ= êÉëáÇì~äëK= mÉêÜ~éë= ~= ãçêÉ= íÜÉçêÉíáÅ~ääó=éêçîÉå=~ééêç~ÅÜ=áë=íç=é~ëë=íç=~=éêçÄ~ÄáäáëíáÅ=íêÉ~íãÉåí=çÑ=íÜÉ=ìåÅÉêí~áåíó=Ó= íç= Éëíáã~íÉ= íÜÉ= éêçÄ~Äáäáíó= çÑ= ÉñÅÉÉÇáåÖ= ~= ëÉí= çÑ= ÅêáíáÅ~ä= äÉîÉäëK= fí= ãÉ~åë= íç= ÅÜ~åÖÉ= ~= êÉJ ÖêÉëëáçå=éêçÄäÉã=íç=~=Åä~ëëáÑáÅ~íáçå=çåÉK=^=êÉèìáêÉÇ=Åä~ëëáÑáÉê=Å~å=ÄÉ=ÅçåëíêìÅíÉÇ=Ä~ëÉÇ=çå= íÜÉ=ë~ãÉ=ã~íÜÉã~íáÅ~ä=éêáåÅáéäÉë=Ó=m~êòÉåÛë=åçåJé~ê~ãÉíêáÅ=ÇÉåëáíó=Éëíáã~íáçåK=^=éêçÄJ ~ÄáäáëíáÅ= ÇÉëÅêáéíáçå= çÑ= ìåÅÉêí~áåíó= áë= ãçêÉ= áåíÉêéêÉí~ÄäÉ= ~åÇI= íÜìëI= Å~å= ÄÉ= ÉîÉå= ãçêÉ= ~ííê~ÅíáîÉ=Ñçê=ÇÉÅáëáçå=ã~âÉêë=íÜ~å=ÇáêÉÅí=î~äìÉë=çÑ=ìåÅÉêí~áåíóK=

5. CONCLUSIONS qÜÉ= ÅìêêÉåí= ëíìÇó= äÉÇ= íç= íÜÉ= ã~áå= ÅçåÅäìëáçå= íÜ~í= dokk= áë= ~= éêçãáëáåÖ= íççä= Ñçê= ~ìíçJ ã~íáÅ= ëé~íá~ä= éêÉÇáÅíáçå= çÑ= ê~Çáç~Åíáîáíó= áå= ÄçíÜ= êçìíáåÉ= ~åÇ= ÉãÉêÖÉåÅó= ëáíì~íáçåëK= fí= áë= îÉêó= Ñ~ëíI= É~ëó= íç= ~ééäó= ~åÇ= áåíÉêéêÉíK= fí= éêçîáÇÉë= êÉ~ëçå~ÄäÉ= êÉëìäíëI= ÉîÉå= Ñçê= ÅçãéäÉñI= ìåÉñéÉÅíÉÇ=Ç~í~K=qÜÉ=ëãççíÜáåÖ=ÉÑÑÉÅí=áë=~=îáëáÄäÉ=Çê~ïÄ~ÅâK=_ìí=áí=Å~å=ÄÉ=áÖåçêÉÇ=áÑ=ïÉ= ~êÉ=åçí=íçç=ãìÅÜ=áåíÉêÉëíÉÇ=éêÉÇáÅíáåÖ=~ÅÅìê~íÉäó=ÜáÖÜ=~åÇ=äçï=î~äìÉëI=Äìí=ï~åí=íç=ÇÉíÉÅí= íÜÉ=éêÉëÉåÅÉ=çÑ=~=Üçí=ëéçí=~åÇ=áåÇáÅ~íÉ=íÜ~í=~=ÅêáíáÅ~ä=äÉîÉä=áë=ÉñÅÉÉÇÉÇK=tÉ=çåäó=åÉÉÇ=íç= âÉÉé=áå=ãáåÇ=íÜ~í=ã~ñáãìãë=ïçìäÇ=~äï~óë=ÄÉ=ìåÇÉêÉëíáã~íÉÇ=Ó=ëãççíÜÉÇK= qÜÉ=ÇáêÉÅí=Éëíáã~íáçå=çÑ=íÜÉ=éêÉÇáÅíáçå=ìåÅÉêí~áåíó=Ü~ë=ëçãÉ=íÜÉçêÉíáÅ~ä=éêçÄäÉãëK=fåáJ íá~ääó=íÜÉ=ãÉíÜçÇ=ÇçÉë=åçí=éêçîáÇÉ=~=ãÉ~ëìêÉ=çÑ=ìåÅÉêí~áåíóX=áí=áë=Å~äÅìä~íÉÇ=ìåÇÉê=ê~íÜÉê= ëíêçåÖ= ~ëëìãéíáçåëK= _ìí= íÜÉêÉ= ~êÉ= ëçãÉ= ï~óëI= ~äëç= ÇáëÅìëëÉÇ= áå= íÜáë= ïçêâI= íç= áåíêçÇìÅÉ= ìåÅÉêí~áåíó= Éëíáã~íáçå= áå= íÜÉ= éêçÅÉëëK= qÜÉ= çíÜÉê= éçëëáÄáäáíó= áë= íç= ãçîÉ= íç= ~= éêçÄ~Äáäáíó= íêÉ~íãÉåí=çÑ=ìåÅÉêí~áåíóK=

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES

19-13

CODES qÜÉ= ëçìêÅÉ= ÅçÇÉë= çÑ= ~ää= dokkJêÉä~íÉÇ= ~äÖçêáíÜãë= ~åÇ= éêçÅÉÇìêÉë= ïÉêÉ= í~âÉå= Ñêçã== j~ëíÉêë= ENVVRFK= qÜÉó= ïÉêÉ= ~ëëáãáä~íÉÇ= áå= ~= ëéÉÅá~äI= ~ìíçã~íáÅ= dokk= íê~áåáåÖ= ~åÇ== éêÉÇáÅíáçå=ëóëíÉãK=

REFERENCES _áëÜçéI=`KjKX=kÉìê~ä=kÉíïçêâë=Ñçê=m~ííÉêå=oÉÅçÖåáíáçåK=lñÑçêÇW=`ä~êÉåÇçå=mêÉëëX=NVVRK= e~óâáåI=pKX=kÉìê~ä=kÉíïçêâW=^=`çãéêÉÜÉåëáîÉ=cçìåÇ~íáçåK=pÉÅçåÇ=bÇáíáçåK=kÉï=vçêâW== j~Åãáääá~å=mìÄäáëÜáåÖ=`çãé~åóX=NVVVK= h~åÉîëâóI=jKX=Dpé~íá~ä=mêÉÇáÅíáçåë=çÑ=pçáä=`çåí~ãáå~íáçå=rëáåÖ=dÉåÉê~ä=oÉÖêÉëëáçå=kÉìê~ä== kÉíïçêâëDK=póëíÉãë=oÉëÉ~êÅÜ=~åÇ=fåÑçK=póëíÉãëI=NVVU=j~êX=UW=éé=OQNJORSK= h~åÉîëâóI=jKX=D^êíáÑáÅá~ä=kÉìê~ä=kÉíïçêâë=~åÇ=pé~íá~ä=fåíÉêéçä~íáçåëK=`~ëÉ=ëíìÇóW=`ÜÉêåçÄóä== Ñ~ääçìíDK=mêÉéêáåí=VRJM=jçëÅçïW=kìÅäÉ~ê=p~ÑÉíó=fåëíáíìíÉX=NVVRK= j~ëíÉêëI=qKX=^Çî~åÅÉÇ=^äÖçêáíÜãë=Ñçê=kÉìê~ä=kÉíïçêâëK=^=`HH=pçìêÅÉÄççâK=qçêçåíçW=gçÜå=táäÉó= C=pçåëI=fåÅX=NVVRK= k~Ç~ê~ó~I=bK^KX=Dlå=Éëíáã~íáåÖ=êÉÖêÉëëáçåDK=qÜÉçêó=çÑ=mêçÄ~Äáäáíó=~åÇ=áíë=^ééäáÅ~íáçåëI=NVSQX=VW=ééK= NQNJNQOK= m~êòÉåI=bKX=Dlå=Éëíáã~íáçå=çÑ=~=éêçÄ~Äáäáíó=ÇÉåëáíó=ÑìåÅíáçå=~åÇ=ãçÇÉDK=^åå~äë=çÑ=j~íÜÉã~íáÅ~ä= pí~íáëíáÅëX=NVSOX=PPW=ééK=NMSRJNMTSK= oçëÉåÄä~ííI=jKX=DaÉåëáíó=Éëíáã~íÉë=~åÇ=j~êâçî=ëÉèìÉåÅÉëDK=fåW=mìêáI=jKI=ÉÇáíçêK=kçåé~ê~ãÉíêáÅ= qÉÅÜåáèìÉë=áå=pí~íáëíáÅ~ä=fåÑÉêÉåÅÉK=içåÇçåW=`~ãÄêáÇÖÉ=råáîÉêëáíó=mêÉëëX=NVTMK=ééK=NVVJONP= oçëÉåÄä~ííI=jKX=DoÉã~êâë=çå=ëçãÉ=åçåé~ê~ãÉíêáÅ=Éëíáã~íÉë=çÑ=~=ÇÉåëáíó=ÑìåÅíáçåDK=^åå~äë=çÑ= j~íÜÉã~íáÅ~ä=pí~íáëíáÅëI=NVRSX=OTW=ééK=UPOJUPTK= péÉÅÜíI=aKbKX=D^=dÉåÉê~ä=oÉÖêÉëëáçå=kÉìê~ä=kÉíïçêâDK=fbbb=qê~åë~Åíáçåë=çå=kÉìê~ä=kÉíïçêâëI= NVVNX=OW=ééK=RSUJRTSK= t~íëçåI=dKpKX=DpãççíÜ=êÉÖêÉëëáçå=~å~äóëáëDK=p~åâÜó~W=qÜÉ=fåÇá~å=gçìêå~ä=çÑ=pí~íáëíáÅëI=pÉêáÉë=^I= NVSQX=OSW=ééK=PRVJPTOK=

Cite this article as: Timonin, Vadim; Savelieva, Elena. ‘Spatial prediction of radioactivity using general regression neural network’. Applied GIS, Vol 1, No 2, 2005. pp. 19-01 to 19-14. DOI: 10.2104/ag050019 Copyright © 2005 Vadim Timonin and Elena Savelieva

19-14

SPATIAL PREDICTION OF RADIOACTIVITY ARTICLES