검색 전체 메뉴
PDF
맨 위로
OA 학술지
Imaging Performance Analysis of an EO/IR Dual Band Airborne Camera
  • 비영리 CC BY-NC
  • 비영리 CC BY-NC
ABSTRACT
Imaging Performance Analysis of an EO/IR Dual Band Airborne Camera
KEYWORD
MTF , Airborne sensor , Optomechanics , EO , IR , (120.4880) Optomechanics , (120.4820) Optical systems , (280.4788) Optical sensing and sensors
  • I. INTRODUCTION

    Satellites or aircrafts have been used for remote sensing for many years. The most common payloads for both spaceborne and air-borne remote sensing are optical cameras or telescopes, which form images (of Earth’s surface) on image detectors [1-4]. The most common detectors for optical imaging are charge coupled devices (CCDs) and infrared detectors [5-8]. The distinction between the two is working wavelength or spectral quantum efficiency: A CCD can sense light in the near ultraviolet (UV) down to near IR, which is often called the electro-optic (EO) spectral band, and IR detectors can sense the IR spectrum. Fig. 1 shows a typical spectral response of a CCD (e2V AT71201M [6]) and a mid-wave HgCdTe IR detector [8].

    An airborne sensor is developed for remote sensing on an aerial vehicle (AV). The sensor is an optical payload for an electro-optical/infrared (EO/IR) dual band camera that combines visible and IR imaging capabilities in a compact and lightweight manner, as similarly reported by the previous studies [2, 7, 9-11]. It adopts a Ritchey-Chretien telescope for the common front end optics with several relay optics that divide and deliver the EO and IR bands to a CCD and an IR detector, respectively. The EO/IR camera for dual bands is mounted on a two-axis gimbal that provides stabilized imaging and precision pointing in both the along and crosstrack directions. Fig. 2 shows a schematic of the expanded view of the EO/IR dual band camera.

    The UV payload should be designed to sustain its imaging requirements over hostile environments, such as temperature change between on-ground and in-air at its flying altitude, gravity acceleration due to the UAV’s maneu-

    vering, and jitter movement caused by the AV propulsion[12]. In addition, stringent constraints are imposed on the weight and volume of the payload design.We first investigate the mechanical deformations, displacements and stress of the EO/IR camera through finite element analysis (FEA). The optical performance, represented by the modulation transfer function (MTF), is then predicted by integrating the FEA results into optics design/analysis software(ZEMAX) [13]. The analysis is done for gravity and thermoelastic deformations.

    We first investigate the mechanical deformations, displacements and stress of the EO/IR camera through finite element analysis (FEA). The optical performance, represented by the modulation transfer function (MTF), is then predicted by integrating the FEA results into optics design/analysis software(ZEMAX) [13]. The analysis is done for gravity and thermoelastic deformations.

    II. FEA MODEL & ANALYSIS

       2.1. Optical Model

    The EO/IR dual band camera is basically a Ritchey-Chretien telescope consisting of a primary mirror and a secondary mirror. The entrance pupil of the camera is about 30 cm in diameter. The telescope is then followed by a beam splitter which splits the incoming beam to EO and IR channels.Each beam is then forwarded by relay optics to its corresponding detector. Fig. 3 shows the schematic optical layout of the camera. Table 1 lists the primary parameters of the EO and IR channels.

       2.2. FEA Model and Boundary Conditions

    In order to calculate the surface deformation and displacements precisely, a FEA model of the EO/IR camera is expressed in the commercial software I-DEAS. Fig. 4 shows the isometric and bottom views of the FEA model. The boundary conditions of the model are set by fixing the locations of the gimbal mounts, which are expressed as dotted circles in Fig. 4. The directions of the three axes that are referred to later in this paper are also drawn in Fig. 4.

    [TABLE 1.] Principal parameters of the EO and IR channels

    label

    Principal parameters of the EO and IR channels

       2.3. Analysis

    The analysis is performed for five cases: three gravitational effects and two thermal conditions. For investigating gravitational effects, one gravitational acceleration (1 g) is given along each of the +X, +Y and +Z directions. The two thermal conditions are the overall temperature change to 30℃ from 20℃ and the temperature gradient across the y axis from -5℃ to +5℃.

    The FEA model is solved by a commercial FEA solver called ABAQUS. The optical surface deformation is then analyzed using the Zernike polynomial decomposition method [14, 15]. A wavefront distortion φ (r, θ) over a circular aperture of an arbitrary radius R is expanded as a linear sum of Zernike polynomials Zi(ρ, θ) given by

    image

    with ρ=r/R and the coefficients α i? being given by

    image

    where the weighting function W(ρ) defining the unit circular aperture is given by

    image

    [TABLE 2.] Deformations of the primary and secondary mirrorsdue to 1 g gravitational force.

    label

    Deformations of the primary and secondary mirrorsdue to 1 g gravitational force.

    The first three Zernike polynomials are called the piston,the tip and the tilt. Since those terms represent the rigid-body motion (RBM) of a lens or mirror, those terms are often removed prior to optical calculations such as peak-to-valley(PV) and root-mean-squared (RMS) wavefront errors of a single optical surface.

    III. FEA RESULTS

       3.1. Gravitational Effects

    The first investigation checks the structural stiffness. Fig. 5 shows the stress due to 1g gravity along the +X direction.The maximum stress, which occurred at the flexure mount,is about 3.4 MPa, which is quite far away from the yield stress of the flexure material (invar), i.e., 145 Gpa. The FEA analysis also confirms the mechanical stiffness of the model over the +Y and +Z gravity.

    Table 2 summarizes the deformations of the primary and

    secondary mirrors. The primary mirror is found to be more vulnerable to gravity. Fig. 6 shows the surface deformations of the primary and secondary mirrors with the rigid body motion removed.

       3.2. Thermal Effects

    The thermal effects are investigated for two conditions:a) the overall temperature change to 30℃ from 20℃ and b) the temperature gradient across the Y axis from -5℃ to+5℃. The analysis shows that the maximum stress is smaller than the yield strength of the material for all thermal conditions. Fig. 8 shows the stress due to the overall temperature change. Fig. 7 shows the contour plots of the primary and secondary mirrors for the two thermal conditions. The RMS displacements of the primary and secondary mirrors are

    9.1 and 1.11 nm for the overall temperature change, and 1.12 and 0.23 nm for the temperature gradient condition.

    IV. IMAGING PERFORMANCE PREDICTION

       4.1. Method

    Conventionally FEA results are used to check whether the designed tolerances are satisfied [16]. This traditional approach is a top-down process which considers how simultaneous changes in multiple parameters interact in a statistical way using the Monte Carlo approach [17]. The Monte Carlo approach determines an optimum set of tolerance for a given optical product by taking into account the manufacturing and assembly variations that occur in mass production. We investigate optical performance by integrating the FEA results into optics design/analysis software. This way predicts the imaging performance of the camera more precisely than statistical prediction.

    Most lens design and analysis software can import optical surface deformations, which are from FEA analysis or metrological measurement, into a lens design. In this paper we use ZEMAX [13]. In ZEMAX, optical surface deformation can be defined by a surface type called Zernike Standard Sag [18]. The Zernike Standard Sag surface is defined as an even aspheric surface plus aspheric terms defined by the Zernike Standard coefficients.

    image

    where c is the curvature, r is the radial coordinate, k is the conic constant, αi is the even aspheric coefficient, and N is the number of Zernike coefficients in the series. The rigid body motions of mirrors or lens are expressed by the change of thickness along the optical axis, decenter and tilt of the elements. Table. 3 and 4 list sample data fo

    [TABLE 3.] Rigid body motions of the optical elements of the EO channel due to 1g gravity along the X axis

    label

    Rigid body motions of the optical elements of the EO channel due to 1g gravity along the X axis

    [TABLE 4.] Zernike polynomial coefficients or optical deformations of the IR channel optics due to 1 g gravity along the Y axis

    label

    Zernike polynomial coefficients or optical deformations of the IR channel optics due to 1 g gravity along the Y axis

    integrating FEA results with ZEMAX.

       4.2 Predicted Imaging Performance

    The modulation transfer function (MTF) is used as the performance metric in this paper, and refocus or the longitudinal shift of the image plane is applied to allow compensation. Fig. 9 shows the MTF variations of the EO channel for 1g gravity along the Z axis. The imaging performance of the EO channel drops by 9% and then back to the design value with refocus. Both the EO and IR channels are simulated to keep the design performance with refocus capability regardless of the gravity direction.Fig. 10 shows the MTF variations of the EO channel for the temperature gradient across the y axis (ΔT=10℃). The EO channel shows about a 20% MTF drop after the surface deformations are considered, which means that some optical surfaces are distorted. However, the distortion is shown to mostly be compensated for by refocusing the focal plane.Unlike the EO channel, the IR channel does not show any significant change.

    Figs. 11 and 12 show the MTF variation of the EO and IR channels for the overall temperature change (ΔT=20℃).The EO and IR channels experience a significant performance drop due to the overall temperature change of the camera before applying refocus. The performance drop is worse in the EO channel than in the IR channel. Those results require that the EO/IR camera be refocused on air, and an application of active thermal control is desirable.

    Table 5 summarizes the predicted MTF at the Nyquist frequency of each detector. The Nyquist frequencies for the EO and IR channels are 27.3 and 13.0 line-pairs-per-mm (lpm),respectively. First, the IR channel can sustain imaging perfor-

    [TABLE 5.] Predicted imaging performance

    label

    Predicted imaging performance

    mance as good as designed with refocus capability for all gravitation and thermal conditions. The EO channel experiences slight performance degradation (16% MTF drop) for a 20℃overall temperature change. This result might require the EO/IR camera to be actively thermal controlled.

    V. CONCLUSION

    This paper reports an integrated performance prediction of an EO/IR dual band camera for an AV. The camera is basically a Ritchey-Chretien telescope with ~ 30 cm entrance pupil diameter with relay optics delivering visible and IR images to the detectors (a CCD and an IR detector).

    The reported EO/IR camera is shown by finite element analysis to be mechanically stiff over any directional gravitational force and two thermal conditions. The two thermal conditions are overall temperature change to 30℃ from 20℃and temperature gradient across the y axis from -5℃ to +5℃.

    An investigation over the MTF variation of the camera is performed by integrating the FEA results into lens design/analysis software ZEMAX. The analysis shows both the EO and IR channels mostly keep imaging performance as good as designed with refocus capability. However, the EO channel is shown to be slightly susceptible to overall temperature change and it is expected to experience slight performance degradation from 74% to 59% in terms of MTF for 20℃ overall temperate change. Further investigation will be performed over active thermal control on the EO/IR camera.

참고문헌
  • 1. Schowengerdt R. A 2007 Remote Sensing google
  • 2. Cox C, Kishner S, Whittesy R 2005 “Reconnaissance payloadsfor responsive space” google
  • 3. Zhang D, Ma G, Liu H “Aerial reconnaissance camerausing optical drum and mirror rotation mechanicsm” P.415-418 google
  • 4. Lee J. H, Jang T. S, Yang H.-S, Rhee S.-W 2008 “Optical design of a compact imaging spectrometer for STSAT3” [J. Opt.Soc. Korea] Vol.12 P.262-268 google
  • 5. Stanton M, Phillips W. C, Li Y, Kalata K 1992 The detective quantum efficiency of CCD and vidicon-based detectors for X-ray crystallographic applications [J. Appl. Cryst.] Vol.25 P.638-645 google cross ref
  • 6.
  • 7.
  • 8. Henini M, Razeghi M 2002 Handbook of Infrared Detection Technologies google
  • 9. Beran S. R, Partynski A. J, Baker A. M, Jones J 2002 “Cassegrain optical system for framing aerial reconnaissance camera” google
  • 10. Petrushevsky V, Karklinsky Y, Chernobrov A 2003 “ELOpEO/IR LOROP camera: image stabilization for dual-band whiskbroom scanning photograpy” [Proc. SPIE] Vol.4820 P.607-617 google
  • 11. Seong S, Yu J, Ryu D, Hong J, Yoon J.-Y, Kim S.-W, Lee J.-H, Shin M.-J 2009 “Imaging and radiometric performance simulation for a new high performance dual band airborne reconnaissance camera” [Proc. SPIE] Vol.7307 P.730705-1-730705-11 google
  • 12. An Y, Song C.-P, Kuang R.-J, Jin G 2010 “Research onthe device of non-angular vibration for opto-eltronic platform” [Proc. SPIE] Vol.7659 P.76590H-1-76590H-6 google
  • 13.
  • 14. Noll R. J 1976 Zernike polynomials and atmospheric turbulence [J. Opt. Soc. Am.] Vol.66 P.207-211 google cross ref
  • 15. Roddier N 1990 “Atmospheric wavefront simulation and Zernike polynomials” [Proc. SPIE] Vol.1237 P.668-679 google
  • 16. Smith W. J 2005 Modern Lens Design google
  • 17. Sasak T, Shinkai M, Higashiyama K, Tanaka F, Kishinami T 1998 “Development of statistical tolerancing system for optical product - virtual PT system and mass product simulation” [Proc. SPIE] Vol.3482 P.528-537 google
  • 18. 2005 ZEMAX Optical Design Program User’'s Guide google
OAK XML 통계
이미지 / 테이블
  • [ FIG. 1. ]  Spectral response of the two most common imagingdetectors for remote sensing; (a) CCD and(b) IR detector.
    Spectral response of the two most common imagingdetectors for remote sensing; (a) CCD and(b) IR detector.
  • [ FIG. 2. ]  Expanded view of the EO/IR dual band sensor.
    Expanded view of the EO/IR dual band sensor.
  • [ FIG. 3. ]  Schematic optical layout of the EO/IR dual band IRcamera.
    Schematic optical layout of the EO/IR dual band IRcamera.
  • [ TABLE 1. ]  Principal parameters of the EO and IR channels
    Principal parameters of the EO and IR channels
  • [ FIG. 4. ]  FEA model of the EO/IR camera.
    FEA model of the EO/IR camera.
  • [ FIG. 5. ]  Stress due to +x gravity.
    Stress due to +x gravity.
  • [ TABLE 2. ]  Deformations of the primary and secondary mirrorsdue to 1 g gravitational force.
    Deformations of the primary and secondary mirrorsdue to 1 g gravitational force.
  • [ Fig. 6 ]  Contour plots of the primary and secondary mirror deformations due to gravitational effects. Rigid body motion removed.Contour maps are in meter scale.
    Contour plots of the primary and secondary mirror deformations due to gravitational effects. Rigid body motion removed.Contour maps are in meter scale.
  • [ FIG. 7. ]  Contour plots of the primary and secondary mirror deformations due to overall temperature rise and gradient across the Y axis.Rigid body motion removed. Contour maps are in meter scale.
    Contour plots of the primary and secondary mirror deformations due to overall temperature rise and gradient across the Y axis.Rigid body motion removed. Contour maps are in meter scale.
  • [ FIG. 8. ]  Stress due to overall temperature change to 30℃ from 20℃.
    Stress due to overall temperature change to 30℃ from 20℃.
  • [ TABLE 3. ]  Rigid body motions of the optical elements of the EO channel due to 1g gravity along the X axis
    Rigid body motions of the optical elements of the EO channel due to 1g gravity along the X axis
  • [ TABLE 4. ]  Zernike polynomial coefficients or optical deformations of the IR channel optics due to 1 g gravity along the Y axis
    Zernike polynomial coefficients or optical deformations of the IR channel optics due to 1 g gravity along the Y axis
  • [ FIG. 9. ]  MTF variation of the EO channel as the FEA results are considered and then refocused due to 1 g gravity along the Z axis.
    MTF variation of the EO channel as the FEA results are considered and then refocused due to 1 g gravity along the Z axis.
  • [ FIG. 10 ]  MTF variation of the EO channel as the FEA results are considered and refocused due to the temperature gradient across y axis (ΔT=10℃).
    MTF variation of the EO channel as the FEA results are considered and refocused due to the temperature gradient across y axis (ΔT=10℃).
  • [ FIG. 11. ]  MTF variation of the EO channel as the FEA results are considered and refocused due to the overall temperature change(ΔT=20℃).
    MTF variation of the EO channel as the FEA results are considered and refocused due to the overall temperature change(ΔT=20℃).
  • [ FIG. 12. ]  MTF variation of the IR channel as the FEA results are considered and refocused due to the overall temperature change (ΔT=20℃).
    MTF variation of the IR channel as the FEA results are considered and refocused due to the overall temperature change (ΔT=20℃).
  • [ TABLE 5. ]  Predicted imaging performance
    Predicted imaging performance
(우)06579 서울시 서초구 반포대로 201(반포동)
Tel. 02-537-6389 | Fax. 02-590-0571 | 문의 : oak2014@korea.kr
Copyright(c) National Library of Korea. All rights reserved.