Vol49_1_2006def 45 ANNALS OF GEOPHYSICS, VOL. 49, N. 1, February 2006 Key words radiometric calibration – MIVIS – at- mospheric correction – MODTRAN 1. Introduction A quality index of radiometrically calibrat- ed and atmospherically corrected data can be obtained by comparing spectral ground truths acquired at the same time of flight (Clark et al., 2002). This ground truths area must be characterized by a good spectral and spatial homogeneity and a high reflectance to assure a high signal perceived at the sensor. Such comparison between spectra measured by the sensor and ground truths has been applied to data collected during last MIVIS (Multispectral Infrared and Visibile Imaging Spectrometer, AA500, Daedalus Inc.) campaigns performed by CNR. MIVIS is a whiskbroom scanner in which the light collimated by the scan head is distributed by means of dichroic filters in- to 4 spectrometers that simultaneously record the radiation coming from the Earth (table I). The scanner geometry is characterized by a FoV of 71.06°, an IFoV 2 mrad and a scan rotation fre- quency varying from 6.25 Hz to 25 Hz; data are digitized at 12 bit (Daedalus, 1993). MIVIS standard radiometric calibration se- quence is carried out on a Test Bench (TB) made up of a black painted aluminium box with two calibration halogen lamps and a reference panel of barium sulphate (Spectraflect, reflectance cer- Laboratory activity for a new procedure of MIVIS calibration and relative validation with test data Cristiana Bassani (1), Rosa Maria Cavalli (2), Angelo Palombo (1), Stefano Pignatti (1)(2) and Fabio Madonna (1) (1) Istituto di Metodologie per l’Analisi Ambientale (IMAA), CNR, Tito Scalo (PZ), Italy (2) Laboratorio Aereo Ricerche Ambientali (LARA), IIA-CNR, Tor Vergata (RM), Italy Abstract Remotely sensed data, recorded by means of the MIVIS hyperspectral scanner in the framework of the research activity of the CNR Institutes IIA-LARA and IMAA, have been calibrated to reflectance values and then quan- titatively compared with ground data. A new procedure for radiometric calibration has been defined by utilizing the MIVIS test-bench and applying a wider radiance range with respect to the one provided by the manufactur- ing company. New calibration curves have been determined and applied in the pre-processing chain. For valida- tion purpose ground spectra were measured during the campaign by means of a portable spectroradiometer. The atmospheric correction has been carried out by implementing an IDL procedure to manage MODTRAN4 input and output cards. MIVIS test data acquired over Passo Corese (Roma) have shown how the new calibration co- efficients significantly improve the radiometric accuracy. In particular, in the VIS spectral region the percentage error, with respect to a ground truth spectrum, is about half of that occurring if the standard calibration coeffi- cients are used. Mailing address: Dr. Cristiana Bassani, Istituto di Metodologie per l’Analisi Ambientale (IMAA), CNR, C.da S. Loja, 85050 Tito Scalo (PZ), Italy; e-mail: bassa- ni@imaa.cnr.it 46 Cristiana Bassani, Rosa Maria Cavalli, Angelo Palombo, Stefano Pignatti and Fabio Madonna tified in 1994) of known reflectance. The MIVIS sensor is mounted on the bench in the lab for the automatic Acceptance Test Procedure (ATP), which is run for a fast and periodic scanner cali- bration and a performance check (signal to noise ration – SNR in fig. 1). Such standard calibration associates the DN, measured during the test, to the panel radi- ance values tabulated for each channel, so pro- ducing the calibration curves, straight lines in- tersecting the axis at the origin, relating DN with tabulated radiance values. The straight line slope determines the Radiance Factor (RF) ap- plied to the radiometric calibration to convert the DN to radiance while the standard deviation of the DN values has been considered an esti- mation of the SNR. From 1994 to date, besides the common checks performed at the Deadalus labs, a trial (Lechi, 2000) of reliability of the RF values de- fined by ATP has been carried out as regards possible variations of the TB parameters (lamps and panels). Lechi’s work is not to be consid- ered resolutive, but it certainly represents an outstanding experimental basis to define a more precise procedure. The present work describes a non standard procedure, implemented in the lab using the TB and a reference spectroradiometer to obtain new RF curves on the basis of a wider range of lighting conditions to better describe the natural irradiance values. The new calibration coeffi- cients have been applied to the test image recorded on an cultivated area close to Passo Corese (Rome, Italy) on March 7, 2002 at an al- titude of 1500 m a.s.l (pixel ground resolution 3 m/pixel) where spectral ground truth campaign has been performed for image validation pur- pose. Furthermore, the paper includes an Appen- dix to demonstrate analytically and experimen- tally how the optical defocusing occurring in the TB during ATP does not invalidate the accu- racy of the MIVIS calibration. Thus, the pres- ence of a collimator facing the scan head is not necessary for a correct calibration. Table I. MIVIS spectral characteristics. Spectrometer Bands Spectral range Bandwidth (nm) (nm) 1 1-20 430-830 20 2 21-28 1150-1550 50 3 29-92 1983-2478 9 4 93-102 8180-12700 340÷440 Fig. 1. MIVIS SNR @ 25 Hz, plus the radiance values from which it has been calculated. 47 Laboratory activity for a new procedure of MIVIS calibration and relative validation with test data 2. MIVIS radiometric calibration In the ATP the method used to calculate the panel radiance data for the TB is based on the following equation: (2.1) where Nλ is the panel radiance at wavelength λ, Hλ the panel irradiance at wavelength λ; ϕ the lamp to panel illumination angle; eλ the panel reflectance at wavelength λ, and π the factor for Lambert surface. The physical quantities in- volved in eq. (2.1) are all monochromatic and the λ index corresponds to the central wave- length of each channel of the sensor. In this procedure the Radiance Factors for the first 92 channels are derived from this sim- ple relation (2.2) where DNλ REF1 is the offset of dark current (DN of the first MIVIS internal blackbodies). During the analysis it was noted that the ra- diance received by the instrument during the ATP is much lower than the one measured throughout instrument’s normal operation con- DN DN K N REF1= - m m m m sin N H e2 $ $= r { m m mc m; E dition (summer clear sky) (fig. 2). Therefore, to verify the behavior of the scanner operating un- der radiance conditions similar to those occur- ring during surveys, a new calibration proce- dure was tuned. The developed calibration procedure con- sists in building calibration curves with the use of different radiance values obtained by either varying the power of the lighting source (intro- ducing additional lamps besides those normally used for ATP) or employing different re- flectance panels. In the experimental phase the combination of different irradiating powers and different panels produced 18 experimental con- ditions measured by both MIVIS and the portable spectrometer. The panels employed during non standard calibration are: the Spectraflect (BaSO4) pro- vided with the TB; the Spectralon certified by LabSphere Inc. up to December 2002 (∆ρ ≤ ≤± 0.005 between 300-2200 nm and ∆ρ ≤ ≤± 0.020 between 2250-2500 nm); a 18% Kodak grey panel and the TB black panel. They have been measured by means of the portable ASD FieldSpec Pro FR spectroradiometer (ASD). ASD is a portable spectrometer covering the range between 350 nm and 2500 nm with a spec- tral resolution of 3 nm for the range 350-1000 Fig. 2. Radiances inherent to the radiometric calibration process: Rad_WP: radiance of Spectralon panel meas- ured in open air at our latitude in summer time; Rad_max (maximum radiance) and Rad_min (lowest radiance) measured in TB in the new procedure; Rad_ATP: radiance in the TB during ATP. 48 Cristiana Bassani, Rosa Maria Cavalli, Angelo Palombo, Stefano Pignatti and Fabio Madonna nm, and of 10 nm for the 1000-2500 nm (Milton et al., 1997). In accordance with the last calibra- tion certificate the error associated with the ra- diance measures is: 1.4 × 10−9 W/cm2 nm sr @ @700 nm; 2.4 × 10−9 W/cm2 nm sr @ 1400 nm; and of 10.4 × 10−9 W/cm2 nm sr @ 2100 nm. The DNλ used to express the function rela- tionship with the radiance are not those utilized by the TB software performing the ATP, but they were obtained off-line by reading the DN directly from the raw files recorded on the MIVIS data tape. Each raw file recorded for each panel differently illuminated has a length of about 400 scan lines at a scan rotation fre- quency of 25 Hz. The DNλ used in the RF calculation are the mean computed on a 20 × 20 pixel square, se- lected in the middle of the scanner scan line, once the dark currents are subtracted. This win- dow size was chosen to avoid image blurring due to the reduced distance between the panel and scan head, which does not permit panel images to be well focused, and to allow a significant sta- tistical assessment. The square standard devia- tion takes into consideration both measurement variations typical of the scanner, which is a least Fig. 3. Best fit of the 18 radiance values for channel 21 (second polynomial fit) and channel 31 (linear fit) to derive the calibration curves. Fig. 4. Comparison of Spectralon radiance measured in the TB. Rad_ATP curve is the MIVIS radiance cali- brated with RF from ATP; Rad_new is the MIVIS radiance calibrated with the new procedure; Rad_ASD is the ASD measured inside the TB. 49 Laboratory activity for a new procedure of MIVIS calibration and relative validation with test data contribution, and the possible error due to the non homogeneities in the panel lighting. The 18 couples of values (Radiance/DN) so defined for each channel, positioned on the graph of fig. 3, show that the linear trend, pre- sumed as valid by the TB software during the ATP, is not always perfectly verified. For exam- ple, some points of the channel 21 (fig. 3 on the left) wander off from the hypothesized linear trend at higher radiance values (close to those occurring under MIVIS normal working condi- tions). In the first 28 channels (first and second spectrometer, VIS-NIR regions) the best fitting curve is a second polynomial, while for the MIR the trend is a perfect straight line (fig. 4 on the right). To determine the RF of the non lin- ear channels, instead of a linear relationship, the following equation is used (2.3) where a, b and c are the fit coefficients (Bassani et al., 2002). To ascertain the accuracy of the MIVIS ra- diance values obtained by applying the new cal- ibration they were compared with the radiance values obtained by applying the RF from the ATP, and the measured panel radiance by ASD. Figure 4 depicts how the panel radiance curves, obtained with the new calibration coefficients, better fits the reference ASD curve to those cal- culated by using the ATP. The difference be- tween the ATP and ASD curves could be prima- rily due to variations of the TB lamp irradiance and the weathering conditions of the Spec- traflect panel (Rollin, 2000). As a consequence, MIVIS radiometric correction was performed using the new calibration coefficients. 3. MIVIS atmospheric correction and validation To convert the instrument perceived radiance values to reflectance ρg(λ), an intrinsic property of the surface independent of the illumination and atmospheric conditions at the time of meas- urements, an atmospheric correction procedure was applied to the MIVIS data. In order to cor- DN RAD RADa b c2$ $= + + rect the atmospheric effects, the relationship be- tween the upward radiance, L0(λ), measured by MIVIS and the surface reflectance, ρg(λ), has to be established by using the MODTRAN4 radia- tive transfer model (Slater, 1980; Tanré, 1992; Vermote et al., 1997). Under the assumption of a cloudless sky, a plane-parallel atmosphere and a lambertian sur- face, the general model applied is given by (3.1) where L0(λ) is the radiance measured at the sen- sor; Ld↑(λ) represents the upwelling atmosphere spectral radiance, scattered from the direct sun- beam into the sensor by the atmosphere (path radiance); Ed ↓(λ) is the downwelling spectral irradiance at the surface due to the scattered so- lar flux in the atmosphere; τv↑(λ) is the atmos- pheric transmittance along the path from the ground surface to the sensor; τz↓(λ) is the at- mospheric transmittance along the path from the sun to the ground surface; θz is the solar zenith angle, E0↓(λ) is the solar spectral irradi- ance at the top of the Earth’s atmosphere. For applying the atmospheric correction to MIVIS data, an IDL procedure has been devel- oped to manage input and output cards, to run MODTRAN4 (Berk et al., 1999). The IDL code is a stepwise procedure in which the different atmospheric parameters are computed separate- ly. The first step takes into consideration the contribution of path radiance, Ld↑(λ) along the scan line, while in the following step the quan- tities Ld↑(λ), Ed↓(λ), τv↑(λ), τz↓(λ), E0↓(λ) are computed. The procedure automatically re- samples each MODTRAN output quantity Fi(λ) to MIVIS spectral channels by using the following formula: (3.2) where the index i is the wavelength step select- ed in the MODTRAN simulations (1 nm in our case), j is the MIVIS channels index (i.e. 1-92), ( ) ( ) ( ) ( ) F r F r j i j i n i i n i j 1 1 $ =m m m m = = 6 6 @ @ / / ( ) ( ) ( ) ( ) ( ) ( ) ( ) cosE E L L g v z z d d0 0 = + - t m x m m i x m m r m m - . - .6 6 @ @ 50 Cristiana Bassani, Rosa Maria Cavalli, Angelo Palombo, Stefano Pignatti and Fabio Madonna and [ri (λ)]j is the spectral response function for the j MIVIS channel in the MODTRAN wave- length step. The determined Fj(λ) are applied in eq. (3.1) to obtain the MIVIS reflectance value. The IDL code was applied to correct the MIVIS survey on Passo Corese configuring MODTRAN cards with a mid-latitude summer profile and with a rural aerosol model (visibili- ty equal to 23 km). The validation process of the atmospheric correction is based on the comparison of the MIVIS reflectance with the spectral ground re- flectance spectra acquired with the ASD. To fulfil this task the selection of a wide Areas Of Interest (AOI), as spectrally and spatially ho- mogeneous as possible, is required, so to mini- mize the adjacency effects; moreover high re- flectance values are required in order to have high SNR at the sensor. As regards the image of Passo Corese, the cultivated fields were not used as target reference (mixture of soil and bean plants) because of the Fig. 6. Comparison of APT and new calibration spectrum of the AOI in concrete with respect to ASD re- flectance. ASD measure is plotted with one sigma errors bars. Fig. 5. Comparison of APT and new calibration spectrum of the AOI in concrete with respect to ASD re- flectance. ASD measure is graphed with one sigma errors bars. 51 Laboratory activity for a new procedure of MIVIS calibration and relative validation with test data growing season and, therefore, the AOI were se- lected on a concrete platform covering a total area of about 200 m2. The following graphs in figs. 5, 6 and 7 show the MIVIS reflectance spectra ob- tained respectively by applying the radiometric calibration with the RF derived from the ATP and the new calibration coefficients from the above described new procedure. On the graphs MIVIS spectra are compared with the ASD measure- ments of the AOI. ASD measurements have been averaged and the standard deviation has been considered like the spectral variability of the con- crete spectrum (error bars on the graphs). The comparison of MIVIS and ground spec- tra outlined how the ATP tends to overestimate the reflectance values in the VIS region (fig. 5) where the reflectance value determined by means of the new procedures is to be considered equal within the error to the mean of the ground truths. In terms of percentage with the new cali- bration the reflectance errors are of about 15% with respect to the ASD ground truth, while the ATP carries an error of about 30%. The NIR region is manifestly subject to a contribution of water vapour (fig. 6) as regards the ground truths as well. The high variability of the water vapour content in the ground spec- tra leads to a great uncertainty in determining the reference reflectance value so much that both MIVIS spectra fall within the ASD meas- urements confidence interval. In the SWIR re- gion (fig. 7), the two calibration procedures in- stead show comparable reflectance values, even though both highlight an underestimation of the image reflectance with respect to the reference ground truth spectra. 4. Results and discussion In order to test the accuracy of the RF for the entire instrument dynamic range, a new proce- dure has been devised. Such non-standard proce- dure is able to produce radiance values with a ra- diometric precision higher than that reachable by the ATP procedure. The non-standard procedure, being developed not taking into account the ATP tabulated values, decouples the scanner perform- ance variation from that related to the TB. By comparing the new procedure with the ATP in the data pre-processing chain to obtain re- flectance values, it has been observed that the ATP in the VIS region overestimates the spec- trum of the selected AOI. This result is congruent with the experimental results obtained from the calibration of the TB DN 18 images, showing the overestimation of the panels’ radiance in the ATP. In the NIR and SWIR spectral regions, the ATP and the new calibration procedure attain two sim- ilar curves both occurring nevertheless in the sec- Fig. 7. Comparison of APT and new calibration spectrum of the AOI in concrete with respect to ASD re- flectance. ASD measure is plotted with one sigma errors bars. 52 Cristiana Bassani, Rosa Maria Cavalli, Angelo Palombo, Stefano Pignatti and Fabio Madonna ond spectrometer within the intrinsic AOI vari- ability, while in the third MIVIS optical port they both appear below the ASD reflectance curve. The result is satisfactory if compared to the purpose of this work: the new procedure brings about an indeterminacy in the reflectance meas- urement which is about the half of that occur- ring if the ATP RF are used. Once the experimental phase is over, this procedure can be applied at regular intervals throughout the year to verify, and possibly up- date, the ATP routine calibration. The new radiometric calibration is going to be applied in the next MIVIS airborne campaign to be deployed in Sweden in June 2003, where atmospheric correction will be performed using in situ atmospheric measures (i.e. radiosounding, solar direct measurements for optical thickness), thus permitting the uncertainty of MODTRAN simulations to be reduced. Acknowledgements The authors are grateful to Maurizio Poscol- ieri and Roberto Ponzianelli for valuable dis- cussions and to Angela Mirabelli for helping to improve the paper. Appendix. ATP accuracy as function of the optical defocusing within the TB. To demonstrate how the distance between the TB panel and the MIVIS scan head is not affecting the calibration procedure, a twofold procedure has been followed: the former, analytic, according to the rules of geometric optics; the latter, experimental, to prove the independence of the radiance measures on the distance from a uniform and wide source. A telescope could be schematised with a thin lens (fig. A.1), in which, for distances smaller than the infinite, the image will be displayed at a distance I . The difference between I and the focal dis- tance (F) is named defocalization (D). The distance I, for a point C located at Z distance from the lens, is defined by the following equation (formula valid for thin lens) . (A.1) Therefore, the defocalization is defined as . (A.2) From the above equations, it derives that the defocalization goes asymptotically to zero by increas- ing the distance. The defocusing phenomenon causes a lost of energy passing the field stop, coming from the single point of the viewed surface, that is due to the non coincidence of the image disc and the field stop located on the focal plane (fig. A.1-β). Nevertheless, it is possible to demonstrate that: from a surface, sufficiently large and uniformly illuminated, the energy reaching the sensor after the field stop is not correlated to the distance. A generic element of the surface C at distance Z from the lens, forming an angle θ with the opti- cal axis, is focalised at position Cl and it is spread on the focal plane F, as a circle (disc image) of di- ameter EH. The diameter of the above mentioned circle (2Rc) is yielded by . (A.3) By means of a geometric consideration . (A.4) D I R OCGH t= - l D EG I R OCt= + l EG GHEH = + D Z F F = - 2 I Z F ZF = - 53 Laboratory activity for a new procedure of MIVIS calibration and relative validation with test data Where Rt is the lens radius and OCl is the parallax (P) which is determined by the following relation: . (A.5) Including eq. (A.4) in eq. (A.3), the new relationship is obtained . (A.6) From which it derives that (A.7) Assuming E0 as the energy falling on the telescope from the surface viewed under the infinitesimal solid angle δω, it is spread on the focal plane within a circle of radius Rc. The fraction of E0 passing . EH D I R EH Z FR 2 2 t t= = EH D I R OC D I R OCt t= + + -l l ( )tanP I= j Fig. A.1. Telescope scheme with one convergent lens. On the right (α) is the front view of the field stop (grey disc) and the image of C on the focal plane (light disc patterned with squares); (β) show the relative position and dimension of the discs for condition τ1; (γ) show the relative position and dimension of the discs for τ3; (δ) show the relative position and dimension of the discs for τ2. On β, γ, δ white area point out the two discs intersection. 54 Cristiana Bassani, Rosa Maria Cavalli, Angelo Palombo, Stefano Pignatti and Fabio Madonna the field stop is described by the τ function. The τ function is defined according to three different for- mulas (fig. A.1). In the condition of partial overlap between field stop and disc image τ1 (fig. A.1-β) (A.8) where A represents the overlap area between disc image and field stop and is yielded by the relation (A.9) where P is the parallax and corresponds to the distance between the centres of field stop and disc image. In the condition of total overlap between field stop and disc image and in the case of Rc ≥ Rd (fig. A.1-γ) the τ function is expressed by . (A.10) In the situation of total overlap between field stop and disc image with Rc ≤ Rd (fig. A.1-δ) the τ func- tion is expressed by . (A.11) The condition of Rc = Rd implies that and, therefore, for higher Z values the disc image has dimensions lower than the field stop. The energy fraction ρ| , arriving to the sensor from the whole surface located at distance Z, is ob- tained by integrating the τ function for the entire solid angle, lying within the sensor viewing angle. For energy fraction ρ| is obtained by (A.12) while for by (A.13) The integration boundaries θlmax and θ mmax are the limit angles beyond which the disc image on the focal plane: 1) does not totally overlap the field stop (θlmax); 2) is totally out of the field stop aperture (θ mmax). The geometric condition that identifies θlmax and θ mmax are respectively (A.14) Utilizing eqs. ((A.1), (A.2)) and the first expression of eq. (A.4) and considering that OCl= P = Itan(θ), the next relation is derived (A.15)arctan= . Z R F R max t d+im b larctan= Z R max til b l . EG OC EG OC R 0 d - = - = - l l max max .sen send d2 23 0 1 max = +t r x i i r x i i i i il l m ! # # Z FR Rt d2 maxmax sensen dd 22 12 0 max = +t r x i ir x i i i iim l m ! ## Z FR Rt d# Z FR Rt d= 13 =x R R 2 c d 2 2 =x ( ) ( ) ( ) ( ) arccos arccosA R dR P R R R dR P R R P R R P R R P R R P R R 2 2 2 1 c c c d d d c d c d c d c d c d 2 2 2 2 2 2 2 2 = + - + - + - - + + + + + + + - - c cm m R A c 1 2=x r 55 Laboratory activity for a new procedure of MIVIS calibration and relative validation with test data Numerically resolving the integral of eq. (A.12) and eq. (A.13), the trend of ρ| energy fraction as function of the distance Z is obtained. This ρ| trend, only for distances very close to F, shows a dif- ference from the asymptotic value. The difference is related to the methods applied to resolve the in- tegral; in fact it decreases with the increase of the steps number used in the integration increase. At the distance in the TB used for the MIVIS calibration ≅ 70 cm (distance panel – scan head ≅ 40 cm plus the optical path within the MIVIS from the scan head to primary mirror) the error is lower of 0.2% (fig. A.2). To obtain an experimental verification, an apparatus composed by a Spectralon panel illuminated by two halogen lamps, and the ASD has been arranged. The ASD has been used with two foreoptics (focal distance respectively of 98 mm and 32 mm with radius Rt of 25 mm and 9 mm) mounted on the fiber optic which operate as field stop with a diameter of 1.7 mm. Radiance measurements have been performed varying the foreoptics and panel distance. In table A.I the Rc values, calculated for both foreoptics at the various distances, are shown. The Rc values are ten times the field stop (Rd) dimension. The graph of radiance @ 750 nm (fig. A.2) depicts the non- Fig. A.2. Graph of the radiance trend with respect to the distance from the panel to ASD foreoptics. The curves show the ASD measured trend for FoV 1° (F = 98 mm) and FoV 3° (F = 32 mm), while energy fraction curve rep- resents the numerical solution of eq. (A.12) and eq. (A.13). Table A.I. Rc as a function of the distance between field stop source for the two ASD optics. Z (cm) Rc 98 (mm) Rc 32 (mm) 17.5 14.00 1.65 45.0 5.44 0.64 93.0 2.63 0.31 222.0 1.10 0.13 437.0 0.56 0.07 56 Cristiana Bassani, Rosa Maria Cavalli, Angelo Palombo, Stefano Pignatti and Fabio Madonna functional dependence on the radiance from the distance and the foreoptic field of view. Also in the limit condition (i.e. distance of 10.5 cm optics of 98 mm) in which Rc is equal to 14 mm (≅ 10 times Rd) the radiance ASD measure is basically the same as those measured in the other experimental con- figurations (Rd >>Rc). Therefore, both analytically and experimentally demonstrations show how the optical defocusing, occurring in the TB during ATP, does not invalidate the accuracy of the MIVIS calibration. distribuiti sulle bande dal visibile all’infrarosso termi- co dello scanner aerotrasportato MIVIS del CNR, Rap- porto Interno CNR IIA. MILTON, E.J., D.R. EMERY and C.H. KERR (1997): NERC equipment pool for field spectroscopy. The 21st centu- ry observations and interactions, in Proceedings of the 23rd Annual Conference of the Remote Sensing Soci- ety, Nottingham, U.K., 159-164. ROLLIN, E.M. (2000): Monitoring the average spectral radi- ance of the Daedalus AB532A calibration bench for Da- edalus 1268 ATM scanner, in Activities of the NERC E- quipment Pool for Field Spectroscopy (EPFS) in support of the NERC Airborne Remote Sensing Facility, British Geological Survey, University of Southampton, 33-43. SLATER, P.N. (1980): Remote Sensing – Optics and Optical Systems (Addison-Wesley Publishing Company, Read- ing, Massachusetts, U.S.A.), pp. 575. TANRÉ, D., B.N. HOLBEN and Y.J. KAUFMAN (1992): Atmos- pheric correction algorithm for NOAA-AVHRR Pro- ducts: theory and application, IEEE Trans. Geosci. Re- mote Sensing, 30, 231-248. VERMOTE, E.F., N. EL SALEOUS, C.O. JUSTICE, Y.J. KAUF- MAN, J.L. PRIVETTE, L.A. REMER, J.C. ROGER and D. TANRÉ (1997): Atmospheric correction of visible to middle-infrared EOS-MODIS data over land surfaces: Background, operational algorithm and validation, J. Geophys. Res., 102 (D14), 17131-17141. REFERENCES BASSANI, C., R.M. CAVALLI, A. PALOMBO and S. PIGNATTI (2002): La calibrazione dello scanner MIVIS (Multi- spectral Imaging Visibile and Infrared Spectrometer) nelle campagne di misura realizzate dal CNR nel 2001-2002, in Atti VI Conferenza Nazionale ASITA «Geomatica per l’Ambiente, il Territorio e il Patrimo- nio culturale», vol. I, p. 337. BERK, A., G.P. ANDERSON, L.S. BERNSTEIN, P.K. ACHARYA, H. DOTHE, M.W. MATTHEW, S.M. ADLER-GOLDEN, J.H. CHETWYND JR., S.C. RICHTSMEIER, B. PUKALL, C.L. ALLRED, L.S. JEONG and M.L. HOKE (1999): MOD- TRAN4 radiative transfer modeling for atmospheric correction, in Optical Spectroscopic Techniques and Instrumentation for Atmospheric and Space Research III, SPIE Proc., 3756, 348-353. CLARK, R.N., G.A. SWAYZE, K.E. LIVIO, R.F. KOKALY, T.VV. KING, J.B. DAL TON, J.S. VANCE, BW. ROCWELL, T. HOEFEN and R.R. MCDOUGAL (2002): Surface re- flectance calibration of terrestrial imaging spec- troscopy data: a tutorial using AVIRIS, in Proceedings of the AVIRIS Conference, Pasadena, JPL, NASA. DAEDALUS (1993): Daedalus Manual’s. Maintenance Man- ual AA5000 Multispectral Infrared and Visibile Imag- ing Spectrometer (MIVIS), vol. 2. LECHI, G., (2000): Calibrazione radiometrica su 102 canali