OR/14/014 Radiometric processing

Energy calibration

Energy recalibration, dead time correction and region of interest selection were performed in field during the field tape extraction process. A peak position statistics file is generated. This is used in the extraction process to select the appropriate channel (fractional) range for calculation of the standard channels, using Specdrift (Proprietary CGG Airborne Survey software).

1.1 ROI	Channel	KeV	Label
1	34 – 224	404 – 2805	Tc
2	115 – 131	1369 – 1574	K
3	139 – 155	1652 – 1863	U
4	202 – 233	2415 – 2812	Th
5	-	-	-
6	-	-	-
7	-	-	-
8	139 - 155	1652 – 1863	UP

Cs-137 was windowed from the 256 channel data using the limits below

	Energy Limits		Channel Limits
	Lower	Upper	Lower	Upper
Cs-137	648	732	54.08	61.09

Dead time correction

The calculated standard windows were time normalized to counts per second. This is done by dividing the count rate by the live time. In the case of the cosmic channel which has different counting circuitry with minimal dead time, the channel was normalised by dividing by the sample time. This is also done during the extraction process.

The extracted file was loaded into the database where subsequent radiometric processing was performed.

Cosmic and aircraft background removal

Before the radiometric processing below was performed, the data was filtered to reduce statistical noise. The filtering applied was typically as follows:

• Potassium, Uranium, Thorium, Total count channels were filtered with a Gaussian low pass filter with filter length of 5 samples to reduce statistical noise.
• Radar altitude channel was filtered with a 5 point Gaussian cut off filter so that the altimeter response time matches that of the radiometric data.
• Cosmic channel was filtered with a running average of 20 terms. Due to the increased volume of upward crystal the number of terms for uranium up may be decreased. This was decided after examination of data statistics.

The time normalized and filtered channels were corrected as follows:

Kbg	=	KF	-	KACBG	-	KCOS_F *	Cosmic
Ubg	=	UF	-	UACBG	-	UCOS_F *	Cosmic
Thbg	=	ThF	-	ThACBG	-	ThCOS F *	Cosmic
TCbg	=	TCF	-	TCACBG	-	TCCOS F *	Cosmic
UCAVbg	=	Uav	-	UACBG	-	UCOS_F *	Cosmic

where KF, UF, ThF, TCF, UUPF are time normalised filtered potassium, uranium, thorium, uranium up and cosmic channels; KACBG, UACBG, ThACBG, TCACBG are the aircraft backgrounds for potassium, uranium, thorium, total count; Uav and Thav are the channels filtered by means of a running average.

Radon background removal

The radon correction was then computed using Minty’s method.

Minty’s method (1992) is a technique for estimating the background atmospheric (radon) radiation in Airborne Gamma-ray Spectrometry.

The method uses the observations that:

in modern airborne spectrometric systems the resolution of the detectors at the 214Bi photopeaks at 0.609 MeV (from atmospheric radiation) and 1.76 MeV (from uranium in the ground) is well resolved above the Compton continuum, and;

due to the differences in spectral shapes between the two 214Bi photopeaks it is possible to differentiate low energy airborne 214Bi signals (at 0.609 MeV) from the higher energy terrestrial ones (from uranium, at 1.76 MeV): the low energy 214Bi photopeak at 0.609 MeV is less attenuated than the 214Bi peak at 1.76 MeV;

those photopeaks can then be used to estimate the contributions of radon and uranium to the observed spectrum because thorium and potassium sources do not contribute appreciably to these peak count rates.

The technique of determination of background by full spectrum analysis is applied after the aircraft and cosmic backgrounds are removed from the observed spectrum. It involves computing the radon contribution and the uranium contribution to the low energy (0.609 MeV) and high energy (1.76 MeV) peak count rates using the observed count rates in those energy peaks and the computation of three constants determined from the radon and uranium spectra. In turn, the radon contribution to the background in the standard uranium window is computed and added to the aircraft and cosmic backgrounds to get the total uranium channel background. And finally the computed uranium background is used to estimate the total count and potassium backgrounds. Thorium background is computed independently from the aircraft and cosmic backgrounds.

Reference

Minty, B R S. 1992, Airborne gamma-ray spectrometric background estimation using full spectrum analysis. Geophysics, Vol 57, No. 2, PP. 279–287.

The radon correction values were then removed from the filtered background and cosmic corrected values Kbg, Ubg, Thbg and Tcbg

Calculation of effective height

The height from the filtered height channel was then converted to effective height at standard temperature and pressure as per the following formula:

Where:

h = the observed radar altitude in metres.

T = the measured air temperature in degrees Celcius.

P = the barometric pressure in millibars.

Special stripping

The background corrected count rates in the three windows must be stripped to give the counts in the potassium, uranium and thorium windows that originate solely from potassium, uranium and thorium. The stripping ratios a, b, g, a and g must be determined from measurements over calibration pads. The three principal stripping ratios (a, b, and g) increase with altitude above the ground as shown in Table 2.The background corrected data is then stripped. Before stripping the coefficients are corrected for variation from height as follows:

Where:

The background corrected data was then stripped.

Before stripping the coefficients were corrected for variation from height as follows:

ae	=	a	+	00049 *	he
βe	=	b	+	00065 *	he
ge	=	g	+	00069 *	he

Table 6 – Stripping Ratios

Where:

he is equivalent height above ground at STP.

a, b and g are the stripping ratios calculated at ground level.

ae, be and ge are the corresponding co-efficients at height he and,

a, b, g are the reverse stripping co-efficients derived from the PAD calibrations. Stripping is then applied as follows:

N'K' ,K

= [N'TH' (a' 'e'g' 'e - b 'e) + N'u' (a'b' 'e - g 'e) + N'K' (1 - a'a 'e)] / A

N'U' ,U' = [N'TH' (g'b' 'e - a 'e) + N'u' (1 - b'b 'e) + N'K' (b'a 'e - g)] / A

N

TH ,TH

= [N'TH' (1 - G'g 'e) + N'u (b'g 'e - a) + N'K' (ag' - b)] / A

Where:

A = 1 - g ge - a (ae - g be) - b (be - ae ge).

NK = Observed potassium counts corrected for background. NU = Observed uranium counts corrected for background. NTH = Observed thorium counts corrected for background. NK,K = Stripped counts in K.

NU,U = Stripped counts in U. NTH,TH = Stripped counts in Th, and,

ae, be, ge are stripping coefficients derived from the pad calibration, and a, b, g are coefficients as described above.

name="_TOC_250029"

The background corrected and stripped count rates are corrected for variations in the altitude of the detector using the following equation:

N'corr

= N'obs' 'e - m (h'o' - h)

Where:

Ncorr = the count rate normalised to the nominal survey altitude, h0. Nobs = the background corrected, stripped count rate at STP height h and

m = the attenuation coefficient for that window.

File:C:\Users\dbk\AppData\Local\Temp\msohtmlclip1\01\clip image002.gif

name="_TOC_250028"

The corrected window count rate data will be converted to ground concentrations of potassium, uranium and thorium using the following expression:

File:C:\Users\dbk\AppData\Local\Temp\msohtmlclip1\01\clip image003.gifC = N

S

Where:

C = concentration of the radioelement (K%, U ppm or Th ppm). S = broad source sensitivity for the window; and,

N = count rate for each window, after dead-time, background, stripping and height correction.

name="_TOC_250027"

Some leveling was required after radiometric correction. This was done using proprietary tie-line leveling and micro-leveling software routines where required.

name="_TOC_250026"

Noise adjusted singular value decomposition (NASVD) was used to considerably enhance low signal to noise regions of the survey.

The technique is an enhancement of the principal component (PC) technique commonly used in Remote Sensing for processing of Landsat and Spot images. The PC technique is a linear transformation of multiband data that generates uncorrelated components.

Each successive component contributes successively less variance to the total response. The first PC can be considered the spectral shape that contributes most to the overall response.

This technique is enhanced by adjusting the variances to fit that expected for a Poisson distribution. This yields components that have physical meaning.

name="_TOC_250025"

The following products have been delivered. It is noted that all products are digital and have been presented in the project coordinates reference system (See Section 2.2).

name="_TOC_250024"

Magnetic Grids

§ tmi.grd

§ tmi_hg_enhanced.grd

§ tmi_igrf.grd

§ tmi_igrf_1VD.grd

§ tmi_igrf_analytical_signal.grd

§ tmi_igrf redp.grd

Radiometric Grids

§ Potassium masked.grd

§ Potassium.grd

§ Thorium masked.grd

§ Thorium.grd

§ Total Count masked.grd

§ Total Count.grd

§ Uranium masked.grd

§ Uranium.grd

Total Magnetic Intensity (TMI) - nT

TMI with Horizontal Gradient enhancement (nT) IGRF Correcxted TMI (nT)

First Vertical Derivative of TMI (nT/m) Analytical Signal of TMI-( nT/m) Reduction to the Pole of TMI (nT)

Processed Potassium with offshore data masked out (%) Processed Potassium (%)

Processed Thorium with offshore data masked out (ppm) Processed Thorium (ppm)

Processed Total Count with offshore data masked out (cps) Processed Total Count (cps)

Processed Uranium with offshore data masked out (ppm) Processed Uranium (ppm)

Other Grid

§ Terrain.grd Digital Terrain Model (m)

name="_TOC_250023"

'FORMAT' OF GEOSOFT MAGNETIC DATABASE

Channel	Description	Units
X	Easting (X) in WGS84 UTM Zone 30N	m
Y X_BNG Y_BNG Lat Lon fid flight date humidity Fluxgate X Fluxgate Y Fluxgate Z Igrf diurnal diurnal_igrf gps_time gps_height	Northing (Y) in WGS84 UTM Zone 30N Easting (X) OSGB 1936 Northing (Y) OSGB 1936 Latitude in WGS84 Longitude in WGS84 Fiducial number Survey flight number Survey flight date Humidity Fluxgate X Fluxgate Y Fluxgate Z) International Geomagnetic Reference Field Magnetic Ground Base Station IGRF Corrected Magnetic Ground Base Station Time in seconds after midnight Altitude above WGS84 Datum	m m m ddd.mm.ss.ss ddd.mm.ss.ss yyyy/mm/dd Degrees C millisecs millisecs millisecs nT nT nT seconds metres
altimeter pressure temperature compensated_mag_left_sensor compensated_mag _right_sensor compensated_mag _tail_sensor raw_horizontal_gradient processed_horizontal_gradient levelled_mag levelled_mag_igrf_corrected inclination declination terrain srtm_data	Radar altimeter height from surface Outside air pressure Outside air temperature Total Magnetic Intensity (compensated) Total Magnetic Intensity (compensated) Total Magnetic Intensity (compensated) Horizontal Gradient (compensated) Horizontal Gradient (processed) Total Magnetic Intensity (Levelled) Residual Magnetic Intensity (Levelled) Magnetic Inclination Magnetic Declination Calculated digital terrain model SRTM	metres mb Degrees C nT nT nT nT nT/m nT nT degrees degrees m m

'FORMAT' OF GEOSOFT RADIOMETRIC DATABASE

Channel	Description	Units
X	Easting (X) in WGS84 UTM Zone 30N	m
Y	Northing (Y) in WGS84 UTM Zone 30N	m
X_BNG	Easting (X) OSGB 1936	m
Y_BNG	Northing (Y) OSGB 1936	m
fid	Fiducial number
flight	Survey flight number
date	Survey flight date	yyyy/mm/dd
STime	Sample time	millisecs
LTime	Live time	millisecs
gps_time	Time in seconds after midnight	seconds
gps_height	Aircraft height above geoid	metres

File:C:\Users\dbk\AppData\Local\Temp\msohtmlclip1\01\clip image004.gif

pressure	Outside air pressure	mb
temperature	Outside air temperature	Degrees C
cosmic	Cosmic radiation	counts
uranium_up	Upward looking Uranium	counts
normalised_total_count	Total count (normalised)	counts
normalised_total_count_nasvd_processed	NASVD corrected total count (normalised)	counts
processed_total_count_nasvd	Final corrected total count	cps
normalised_potassium_nasvd	Potassium (normalised)	counts
normalised_potassium_nasvd_processed	NASVD corrected potassium (normalised)	counts
processed_potassium_nasvd	Final corrected potassium	%

normalised_uranium	Uranium (normalised)	counts
normalised_uranium_nasvd_processed	NASVD corrected uranium (normalised)	counts
processed_uranium_nasvd	Final corrected uranium	ppm
normalised_thorium	Thorium (normalised)	counts
normalised_thorium_nasvd_processed	NASVD corrected thorium (normalised)	counts
processed_thorium_nasvd	Final corrected thorium	ppm
raw_256_down	Downward looking radiometric spectra
raw_256_up	Upward looking radiometric spectra

'FORMAT' OF GEOSOFT TESTLINE DATABASES

Channel	Description	Units
X	Easting (X) in WGS84 UTM Zone 30N	m
Y	Northing (Y) in WGS84 UTM Zone 30N	m
Lat	Latitude in WGS84	ddd.mm.ss.ss
Lon	Longitude in WGS84	ddd.mm.ss.ss
fid	Fiducial number
flight	Survey flight number
date	Survey flight date	yyyy/mm/dd
gps_time	Time in seconds after midnight	seconds
gps_height	Aircraft height above geoid	metres
altimeter pressure	Radar altimeter height from surface Outside air pressure	metres mb
temperature	Outside air temperature	Degrees C
STime	Sample time	millisecs
LTime	Live time	millisecs
TC_nasvd	Total count (normalised)	counts
K_nasvd	Potassium (normalised)	counts
U_nasvd	Uranium (normalised)	counts
Th_nasvd	Thorium (normalised)	counts
TC_processed	Final corrected total count	cps
K_processed	Final corrected potassium	%
U_processed	Final corrected uranium	ppm
Th_processed	Final corrected thorium	ppm
raw_256_down	Downward looking radiometric spectra
raw_256_up	Upward looking radiometric spectra