Design of Matching Imaging on Agile Satellite with Wide-Swath Whiskbroom Payloads along the Coastal Zone

: The non-linearly curved coastal zone is very long and wide. The traditional satellite can’t cover the whole coastal zone in a single scan. So, the method of matching imaging combining the adjustment of satellite attitude and the scanning of the swing mirror are proposed. Firstly, based on the position of feature points of the coastal zone, the attitude when the optical axis of the satellite always points to the feature points is calculated. According to the width of the sea and land on both sides of the coastline, the parameters of the swing mirror of wide-swath whiskbroom payloads are analyzed. Secondly, the velocity vector model considering time-varying satellite attitude and the dynamic scanning of the swing mirror is constructed. The schemes of matching imaging such as adjustments of yaw angle and detector are developed. Finally, the precise experiment is designed to verify the correctness of the matching imaging. The experimental results show that the resolution of the matching imaging is less than 1 pixel and its modulation transfer function (MTF) is greater than the human eye’s minimum MTF of 0.026. The method of matching imaging using the adjustment of satellite attitude and the scanning of the swing mirror can realize wide imaging along the coastline and improve the temporal resolution at the same time.


Introduction
The coastal zone is not only the area with the most active interaction among land, ocean and atmosphere, but also the area with the most frequent human economic activities in the world.2/3 of the large and medium-sized cities in the world are concentrated along the coastal zone, with 60% of the population and a developed economy.In the coastal zone, 8% of the earth's surface provides 28% of the world's biological production and owns 90% of the global fishery industry.The total length of the global coastal zone is 504,000 km, with nonlinear bending and extension, including a coverage width of ±100 km extending from coastline to continental and oceanic.The coastal zone generally presents the characteristics of length, width and curvature.Its biological species are complex and diverse [1,2].
Satellites can quickly obtain comprehensive data on the coastal environment, ecosystems, landforms, coastline conditions, water quality, biological resources, natural disasters, and human activities [3][4][5][6].It has become an effective tool in regularly monitoring coastal resources and environmental changes.At present, there are two main types of satellites used in observation of coastal zone: wide-coverage push-broom satellite and agile highresolution satellite [7].The wide-coverage push-broom satellite has low resolution.For example, the resolution of SDGSAT-1 [8] is 30 m.It uses a one-dimensional swing mirror to expand coverage, which can obtain thermal infrared images of 300 km.But it can only expand the coverage of the nadir area.The coverage of agile high-resolution satellites is smaller.However, it can search for non-nadir targets by attitude adjustment, such as WorldView [9][10][11], STOP5/6 [12,13], Pleiades [14,15], etc.In this way, the fast revisiting of fine objects can be achieved by agile satellites.What's more, image stitching [16] can be used in agile high-resolution satellites to expand coverage.
Although the wide-coverage push-broom satellite has a large imaging range, it can only image the symmetrical areas on both sides of the satellite's nadir point.The coastal zone is the boundary between the ocean and the land, which shows an irregular curved distribution.It is difficult to match with the image trajectory of wide-coverage push-broom satellite [17].Only part of the coastal zone will be covered by the wide-coverage pushbroom satellite at single scanning.However, if agile high-resolution satellite and image stitching methods are used in imaging the coastal zone, the temporal resolution will be low.It is impossible to quickly cover both water and land on the sides of the coastline.
In addition, both wide-coverage push-broom satellites and agile high-resolution satellites have high-dynamic properties [18], which leads to time-varying and aeolotropic properties of velocity vectors on the image plane.The velocity on the image plane does not match the detector's imaging direction, causing image degradation [19].The velocity vector matching is difficult in those complex imaging modes [20].In 2020, Xu [21] et al. studied velocity vector on the image surface of a space camera for dynamic circular scanning, and carried out numerical and image simulation.In the same year, based on this dynamic imaging model, Wang [22,23] et al. explored the degradation of the image caused by rotation of the image plane.In 2021, Han [24] et al. analyzed the distribution characteristics of velocity vector when satellite scanning Antarctica by combining the characteristics of polar target distribution.In addition, the velocity vector is also applied to approve image quality, one of which is the image compensation method.For example, Sun [25] et al. proposed a method of segmental PSF compensation by calculating the velocity vector on the image plane, which effectively corrected the image motion blur of the KZ-1 agile satellite.The other is on-orbit compensation.Li [26] et al. predicted the yaw angle of satellites through the velocity vector, proposed a method to compensate for the average drift angle by adjusting the satellite yaw angle, and then used the image plane adjustment mechanism in the camera to eliminate the residual error.However, these models and compensation methods are all applicable to traditional push-broom It lacks velocity vectors considering satellite attitude adjustment and swing mirror.
To solve the difficulty of the full coverage of the coastal zone at single scanning by the traditional satellite, a new mode is proposed.The mode expands the coverage through the scanning of the swing mirror [27,28] when the agile satellite adjusts the attitude to image along the coastline [29].However, there will be an enormous mismatch between the velocity and the direction of detector because this new kind of imaging mode has large dynamic maneuvering.The image compensation method is not suitable.So, it is necessary to analyze the distribution of velocity vector on the image plane and develop an on-orbit matching strategy through some hardware adjustment.
The remainder of this article is organized as follows.In Section 2, the attitude of the agile satellites along the coastline is planned, and the maneuver scheme of the swing mirror is introduced.In Section 3, based on the proposed mode of imaging the coastal zone, a velocity vector model is constructed and the corresponding matching strategy is established.In Section 4, experiments are designed to verify the feasibility of the proposed matching imaging.In Section 5, conclusions are summarized.

Attitude Adjustment and Swing Mirror Scanning of Agile Satellites Imaging the Coastal Zone
As shown in Figure 1a, in order to cover the width of the coastal zone of 200 km and follow the natural contours of the coastline, it is necessary to adjust the satellite attitude angles and attitude angular velocities during the flight of the satellite, so that the optical axis always points to the extension direction of the coastline.At the same time, based on satellite attitude adjustment, the parameters of the swing mirror are designed in consideration of time, width, and overlap ratio of vertical-orbit coverage, which is shown in Figure 1b.This section introduces the full coverage technology of the coastal zone with both the water and land from two parts: attitude adjustment of agile satellite and seamless imaging of the swing mirror scanning.
As shown in Figure 1a, in order to cover the width of the coastal zone of 200 km and follow the natural contours of the coastline, it is necessary to adjust the satellite attitude angles and attitude angular velocities during the flight of the satellite, so that the optical axis always points to the extension direction of the coastline.At the same time, based on satellite attitude adjustment, the parameters of the swing mirror are designed in consideration of time, width, and overlap ratio of vertical-orbit coverage, which is shown in Figure 1b.This section introduces the full coverage technology of the coastal zone with both the water and land from two parts: attitude adjustment of agile satellite and seamless imaging of the swing mirror scanning.

Attitude Adjustment of Agile Satellite
There are two processes for the agile satellite scanning along the coast.One is the push-broom along the orbit.The other is side-swing scanning which is perpendicular to the satellite's flight direction.In order to ensure the whole coastal zone is within the field of view and attitude performance is optimal, only roll angle ϕ is adjusted under the condition that the pitch angle remains 0. It is designed to assures that the push-broom trajectory of the agile satellite is consistent with the shape of the coastline during flight.
As shown in Figure 1b, the right-hand Cartesian coordinate system ( ) , ,

Attitude Adjustment of Agile Satellite
There are two processes for the agile satellite scanning along the coast.One is the push-broom along the orbit.The other is side-swing scanning which is perpendicular to the satellite's flight direction.In order to ensure the whole coastal zone is within the field of view and attitude performance is optimal, only roll angle ϕ is adjusted under the condition that the pitch angle remains 0. It is designed to assures that the push-broom trajectory of the agile satellite is consistent with the shape of the coastline during flight.
As shown in Figure 1b, the right-hand Cartesian coordinate system G G x , G y , G z , called Earth centered fixed coordinate system, is established.Its origin is at the center of the earth.G x axis points to the intersection of the Prime Meridian and Equator.G z points to the North Pole.The positions of feature points in the coastal zone are (x T , y T , z T ).The equations of roll angle and roll angle velocity which versus time in G G x , G y , G z are solved.
The position of the satellite in G G x , G y , G z is calculated.The rotation matrices B 1 (•) and B 2 (•) are defined.The parameters of the agile satellite are (r, e, i, A N , ω, M 0 ), and r means the semi-major axis of the satellite.e represents the eccentricity of the satellite orbit.i means the orbital inclination.A N is short for the right ascension of ascending node.ω is the argument of perigee.M 0 is the true perigee angle.The position of the satellite in G G x , G y , G z can be expressed as (x s (t), y s (t), z s (t)), the expression is: (x s (t), y s (t), z s (t)) varies with time t.The distance d(t) between the feature points of the coastline and satellite at time t is iteratively calculated.
In order to keep the pitch angle as 0, the minimum value set of {d m (t m )} is calculated.The time set {t m } m=1,2,...,T can be figured out based on {d m (t m )}.Set {ϕ m {t m }} of roll agile is calculated by the geometry relationship of the triangle in Figure 1a. )) A polynomial formula of ϕ m (t m ) is fitted in [t 1 , t T ] using the least squares method.The planning of the roll angle along the coastline is completed.A n is the fitting coefficient and n is the fitting order.The angular velocity .ϕ(t) of the roll agile is also obtained. .

Seamless Imaging by the Swing Mirror
Adjusting the roll angle of the agile satellite during flight only ensures that the satellite scans along the coastline.As shown in Figure 2, the width S 1 of the coastal zone is large.However, in order to ensure high resolution, the camera of the satellite has a narrow field of view, and its coverage width S 2 is much smaller than S 1 .Therefore, we need to scan in the direction perpendicular to the coastline by the swing mirror when the optical axis always is along the coastline.It works when the swing mirror is added in front of the optical system.There will be a lot of scan strips on the earth and the imaging swath is enlarged by accumulating these strips.The scan strips are shown in Figure 2b.When the swing mirror scans the na the output image is similar to the classic push-broom image.When the tilt an swing mirror is large, due to the influence of the curvature of the earth, the gro of the single frame is approximately trapezoidal.In one period of the swing mi ning, the overall image is a concave quadrilateral image with wider sides.D flight of the agile satellite, the swing mirror is constantly scanning.All strips are lated to cover the entire coastal zone.
In order to avoid non-overlapping between adjacent scanning strips, the speed of the swing mirror needs to be reasonably calculated according t erage of the coastal zone and the exposure time.The coverage width of the sate The optical system with the swing mirror for observing the coastal zone is shown in Figure 2a.It uses both off-axis and on-axis reflective systems.The field of view is FOV.The focal length is f .A reflective mirror is set in front of the optical system to bend the light, which will not change the parameters of the detectors and optical system.The swing mirror scans alternately in the positive and negative directions at a constant speed, where the direction of scanning is perpendicular to the flight of the satellite.As shown in Figure 2b, the black arrows indicate the scanning directions.
The scan strips are shown in Figure 2b.When the swing mirror scans the nadir point, the output image is similar to the classic push-broom image.When the tilt angle of the swing mirror is large, due to the influence of the curvature of the earth, the ground strip of the single frame is approximately trapezoidal.In one period of the swing mirror scanning, the overall image is a concave quadrilateral image with wider sides.During the flight of the agile satellite, the swing mirror is constantly scanning.All strips are accumulated to cover the entire coastal zone.
In order to avoid non-overlapping between adjacent scanning strips, the scanning speed • /s of the swing mirror needs to be reasonably calculated according to the coverage of the coastal zone and the exposure time.The coverage width of the satellite along the orbit during scanning by the swing mirror can be obtained from ∆S A BO in Figure 2b.
The full length of the coastal zone which is imaged is S m.The time of imaging is T = t T − t 1 .The initial overlap rate of two adjacent strips is η.Then, n strips are required to cover the entire coastal zone.The equation of n is: The time for the swing mirror to scan to acquire one strip is ∆t.
As shown in Figure 2b.The extension width of the lands and water on both sides of the coastline is W y .The imaging angle FOV along the scanning direction is solved in the triangle ∆S A AO.
The angular velocity of the swing mirror is .
The directions of angular velocities are opposite when the swing mirror scans the adjacent strips.The time for changing scanning direction of the swing mirror is t .The satellite is still moving forward in time t .There is a forward distance S 1 = µ/(R + H) • t along the flight of the satellite.Then, the actual overlap rate between two adjacent strips is P.
The overlap rate needs to meet P ≥ 10%.It can achieve seamless splicing of scanning strips while expanding the coverage, which is shown in Figure 1b that the whole coastal zone is covered.

Matching Strategy Based on Image Velocity Vector Field
The difference between agile satellites with wide-swath whiskbroom payloads along the coastal zone and traditional push-broom satellites is illustrated as follows.The traditional push-broom satellite has no attitude changes.Its images correspond to the trajectory of the nadir point and its velocity vector are only affected by the satellite precession.However, for agile satellite with wide-swath whiskbroom payloads imaging the coastal zone, satellite attitude varies with time, and the imaging process is complicated after coupling with the swing mirror.The value of the velocity vector on the image plane is constantly changing with time.There is an angle between the velocity vector and the imaging direction of the detector.If measures are not taken to match the velocity vector, the image quality will be seriously affected.It is necessary to formulate on-orbit matching strategies such as yaw angle adjustment and changes of parameters of the detector based on velocity vector model on the image plane.

Velocity Vector Model on the Image Plane
An accurate velocity vector model on the image plane is the key to formulating the velocity vector matching strategy.Firstly, the precise positions of objects are determined through the object-image mapping relationship.Secondly, the velocity vector on the object space is described and decomposed.Finally, the velocity vector on the image plane is obtained by projecting the velocity vector from object to image.As shown in Figure 3 the process involves several right-handed coordinate systems: image coordinate system P P x , P y , P z and object coordinate systems P jk P jkx , P jky .The origin of P P x , P y , P z is at the center of the image plane.P x axis, on the image plane, points to the direction of the satellite flight.P z axis is parallel to the optical axis and points to the center of the earth.(j, k) stands for the serial number of pixels.The row number is j ∈ (−M, M), and the column number is k ∈ (−N, N).There are D pixels in total, corresponding to D object coordinate systems.Their origins are at objects.T jkx axis is tangential to the earth's surface and points to the direction of satellite flight.T jkx axis is perpendicular to the T jky axis, tangential to the earth's surface, and pointing in the direction of the earth's rotation.There is symbolic function:  3.

The Position of Objects
It is very important to solve the position of object which is conjugated to each pixel because that the linear velocities of the earth's rotation, satellite attitude and swing mirror change with the positions of objects.Starting from arrangement rule of pixels and combining pixel angle, optical axis angle and geocentric angle, the latitudes of the objects are acquired through coordinate transformation.
Pixel angle is β jk (t) shown in Figure 3.The focal length of the satellite is f , and the size of pixel is a.The expression for β jk (t) is as follows: There is symbolic function: Optical axis angle is composed of roll angle ϕ(t), pixel angle β jk (t) and scanning angle of the swing mirror η(t).There are η(t) = • t and −FOV /2 ≤ η(t) ≤ FOV /2.The equation of φ jk (t) is: Geocentric angle ρ jk (t) will be acquired by solving ∆OS A T in Figure 3.
In the Formulas ( 14) and ( 15), R is the radius of the earth.The expression of latitude of nadir point δ s (t) can be solved based on Formula (1).
Latitude is the symbol of the position of object.The latitudes of objects δ jk (t), which are affected by satellite attitude and the swing mirror, can be acquired by the spherical triangle ∆NTS A in Figure 3.
There is u (t) = arccos(cos i/ cos δ s ).Under the condition that the satellite is in the prograde orbit (i < 90 • ), when the latitude of the optical axis pointing point is higher than the latitude of the nadir point, '+' will be taken, otherwise, '−' will be taken.

Velocity Vector
Based on the latitudes of objects, it is more convenient to first solve the velocity vector on the object space.Four types of velocity vectors can be obtained in each object coordinate system.

1.
The projected linear velocity of the satellite's flight velocity at the project point: The linear velocity of the Earth's rotation at the object point: v ejk (t) = ω e R cos(δ jk );

3.
The linear velocity of roll angle at the object point:

4.
The line velocity of the swing mirror at the object point: L jk (t) = r 2 + R 2 − 2rR cos(ρ jk (t)) represents the distance from the satellite to object.
ω s is the angular velocity of satellite flight.ω e is the angular velocity of the earth's rotation.
Although the velocity vector in the object coordinate system is obtained, it is difficult to directly analyze the image quality change in the object space.Therefore, these velocity vectors are decomposed into the two axes of the object coordinate system, and then they are projected to the image plane for further analysis.v xjk and v yjk represent the velocity components in the direction of axis T x and axis T x , respectively.
u jk (t) = arccos cos(i)/cos(δ jk (t)) stands for the angle between v sjk and v ejk .Finally, v sjk and v ejk are projected to the image coordinate system by transformation matrix L TPjk (t) to solve the expressions of velocity vectors, v sjk and v ejk , on the image plane.
The value of resultant velocity on the image plane is:

Velocity Vector Matching
In order to take advantage of the widened width of the swing mirror as much as possible, the initial imaging direction of the detector is along the scanning direction of the swing mirror.However, through the above calculation of the velocity vector, it can be found that the direction of the resultant velocity on the image plane is not consistent with the direction of the initial imaging direction of the detector.It will seriously affect the imaging quality, especially when the image plane is stitched by several detectors.Firstly, the mismatch should be analyzed, and then methods such as adjustment of yaw angle, the tilt of detectors, and adjustment of imaging time of detectors will be taken to match velocity vector.
In the Equation (19), in addition to the velocity v pyjk (t) in the initial imaging direction of the detector, there is another velocity v pxjk (t) that is perpendicular to the initial imaging direction of the detector.There is an angle between the resultant velocity on the image plane and the initial imaging direction of the detector, which is called the drift angle.
ξ jk (t) is the mismatching angle of the velocity vector.In addition, the imaging time Ti of detectors should be fit with the value of resultant velocity V pjk (t) on the image plane.The relationship between V pjk (t) and Ti is: For the image plane stitched with Q detectors, the position of the center pixel of each detector in the image plane coordinate system is (q 1 , q 2 ).The following matching strategies are proposed.

1.
Adjustment of yaw angle: first, during imaging along the coast, the satellite continuously adjusts the yaw angle to ensure that the imaging direction of the detector coincides with the direction of resultant velocity in the center of the image plane.The value of yaw angle is Ψ(t) = ξ j=0,k=0 (t); 2. Tilt of detectors: adjustment of yaw angle only matches the velocity vector in the center pixel of the image plane.However, the residual mismatched angles of other pixels still exist.These residual mismatched angles can be compensated by tilting the detectors.The tilt angles of detectors are ∆ξ q 1 q 2 (t) = ξ q 1 q 2 (t) − Ψ(t);

3.
Adjustment of the imaging time of detectors: step 1 and step 2 only consider the direction of the velocity vector.The value of the velocity vector is compensated by adjusting the imaging time of detectors separately.

Experiment
The experiment is performed to verify the image quality of matching imaging on agile satellite with wide-swath whisk-broom payloads along the coastal zone.The experimental facilities are one-dimensional turntable, targets, and lens.The image plane of the lens is stitched by three detectors.The turntable drives the camera to rotate to simulate the scanning of the swing mirror.The lens and detectors simulate the camera of satellite.The curvature of target is proportional to Earth's curvature.The design and results of the experiment are introduced as follows.

Design of Experiment
In order to obtain the parameters of the experiment, we planned the curve of attitude versus time when the agile satellite with the swing mirror scanning the coastal zone of China.The angular velocity of the swing mirror is calculated.The distribution figure of the drift angle on the image plane is drawn.
Firstly, the agile satellite with the sun-synchronous orbit is selected to image the southeastern coastal zone of China.The satellite parameters are shown in Table 1.The focal length of satellite is 3.5 m.The size of single pixel of satellite is 7 µm.Its orbit altitude is 500 km.The ground resolution of the satellite can be calculated as GSD = (H • a)/ f = 1 m.The field of view of the satellite of is 2 • .The coverage width of satellite can be calculated as W x = 17.45 km.According to the satellite orbit parameters and target distribution, the curve of roll angle is demonstrated in Figure 4a.The imaging time region is [2242 s, 2485 s].The focal length of satellite 3.5 m a s The size of single pixel of satellite 7 µm FOV s The field of view of the satellite 2 The coverage width of satellite 17.45 km Secondly, the angular velocity of the swing mirror is calculated based on the curve of roll angle.The total length of China's coastal zone is S = 18, 000 km, with the coverage width of W y = 200 km.In the imaging time of 243 s, the angular velocity of the swing mirror can be obtained by Formula (10) as = 16 • /s.It ensures that the overlap rate of two adjacent strips is 10.33%.
Finally, the projected linear velocity of the satellite's flight, the linear velocity of the Earth's rotation, the linear velocity of roll angle, and the line velocity of the swing mirror are comprehensively analyzed and the velocity vector is calculated by Equation (19).The image plane is spliced by three detectors with size 1024 pixels × 1024 pixels.The number of overlapping pixels is 70. Figure 4c shows the distribution of velocity vectors on the image plane at 2367 s.There are drift angles between those velocity vectors and the imaging direction of the detector.Figure 4b   The experiment simulates the imaging process of the swing mirror at 2367 s 5 shows the imaging principle and facilities of the experiment.The turntable is eq with a horizontally placed camera and rotates at angular velocity to s scanning process of the swing mirror.The lens and detectors are horizontal and th table is tilted at an angle  The experiment simulates the imaging process of the swing mirror at 2367 s. Figure 5 shows the imaging principle and facilities of the experiment.The turntable is equipped with a horizontally placed camera and rotates at angular velocity = 16 • /s to simulate scanning process of the swing mirror.The lens and detectors are horizontal and the turntable is tilted at an angle σ = 1.40 • .The curvature of the spherical surface of target is proportional to the curvature of the earth.The scaling factor is k = H/R.The curvature of the spherical surface of target is r −1 = 28 −1 m −1 when the distance between the target and the lens is h = 2.2 m.All parameters in the experiment are shown in Table 2.
When the horizontally placed lens and detectors rotate with the tilt turntable, the velocity vector on the image plane is v .There is a drift angle σ = 1.40 • between v and the imaging direction of the center detector.It simulates the drift angle at the center pixel (j = 0, k = 0) at time 2367 s.In addition, affected by the curvature of target and the difference in object distances, the drift angle of the side detector is σ = 1.45 • , showing the asymmetrical distribution of the velocity vector on the image plane at time 2367 s.
with a horizontally placed camera and rotates at angular velocity to simulate scanning process of the swing mirror.The lens and detectors are horizontal and the turntable is tilted at an angle 1.40 The curvature of the spherical surface of target is proportional to the curvature of the earth.The scaling factor is k H R = . The curvature of the spherical surface of target is   The following three corresponding matching strategies are taken: 1.
The lens and detectors are rotated by an angle σ, which is equivalent to adjustment of yaw angle;

2.
The side detectors are tilted by an angle ε = −0.5 • to compensate for residual mismatched angles;

3.
Adjustment of the imaging time of detectors to compensate for the value of the velocity vector.The imaging time of the center detector is T i1 = 0.5 ms.The imaging time of the side detector is T i1 = 0.505 ms.

Experiment Procedure and Results
Resolution and modulation transfer function (MTF) are two important indicators for evaluating image quality.The resolution fringes, edge targets, and color targets are imaged through the above experimental methods to verify the feasibility of the matching imaging scheme of agile satellite with the swing mirror.The images acquired from the experiment are shown in Figure 6.

Resolution
As shown in (i) of Figure 5, light-transmitting and opaque strips with a width of A = 0.3 mm are arranged at intervals in resolution target.It is placed in a light-transmitting hole of the camera obscura, which forms bright and dark strip images at intervals.The lens and detectors scans and images the resolution target.The focal length f , the width of strip A, the size of a single pixel a and distance from the lens to target satisfy the resolution imaging relationship A/H =a/ f .Figure 6a,b are the images of resolution targets.Each black or white stripe has one pixel width.

Experiment Procedure and Results
Resolution and modulation transfer function (MTF) are two important indicator evaluating image quality.The resolution fringes, edge targets, and color targets imaged through the above experimental methods to verify the feasibility of the matc imaging scheme of agile satellite with the swing mirror.The images acquired from experiment are shown in Figure 6.

Modulation Transfer Function
As shown in (ii) of Figure 5b, the knife-edge target is attached to the spherical su which is illuminated by a light.The lens and detectors scan and image the knife target.MTF is calculated by one chosen knife-edge.Figure 6c,d

Color image
As shown in (iii) of Figure 5b, the color image is attached to the spherical su which is also illuminated by a light.The lens and detectors scans and images the image.Color images respectively taken by detector 1 and detector 2 are shown in F 6e,f, which correspond to the upper and lower parts of the target.

Discussion
Based on the above three sets of experiments and the experimental image discuss whether the image quality meets the imaging needs from three aspects resolution, MTF, and images of color targets respectively.Firstly, the resolut analyzed.If our matching scheme is valid, there will be clear strips image in the det and each strip occupies one pixel.It can be clearly found that each black or white str the images of detectors has one pixel width.It meets the requirement that resolu less than 1 pixel; Secondly, MTF is expounded.The image obtained by the optical s

Modulation Transfer Function
As shown in (ii) of Figure 5b, the knife-edge target is attached to the spherical surface, which is illuminated by a light.The lens and detectors scan and image the knife-edge target.MTF is calculated by one chosen knife-edge.Figure 6c,d are the images of the knife-edge targe.The chosen knife-edge of detector 1 is MTF 1 = 0.0271.The chosen knife-edge of detector 2 is MTF 2 = 0.0261.

Color image
As shown in (iii) of Figure 5b, the color image is attached to the spherical surface, which is also illuminated by a light.The lens and detectors scans and images the color image.Color images respectively taken by detector 1 and detector 2 are shown in Figure 6e,f, which correspond to the upper and lower parts of the target.

Discussion
Based on the above three sets of experiments and the experimental images.We discuss whether the image quality meets the imaging needs from three aspects from resolution, MTF, and images of color targets respectively.Firstly, the resolution is analyzed.If our matching scheme is valid, there will be clear strips image in the detectors and each strip occupies one pixel.It can be clearly found that each black or white stripe in the images of detectors has one pixel width.It meets the requirement that resolution is less than 1 pixel; Secondly, MTF is expounded.The image obtained by the optical system is finally observed by the human eye.The minimum MTF of the human eye is MTF eye = 0.026.If the MTF of the image is larger than 0.026, it can be considered that the image can be distinguished by the human eye.The MTF of the chosen knife-edge of detector 1 and detector 2 are obtained as MTF 1 = 0.0271 and MTF 2 = 0.0261, respectively.Both of them are larger than 0.026.It meets the requirement of MTF; Finally, images of color targets are discussed.The two images are taken by detector 1 and detector 2 at the same time.In the image, aircraft mounts and boundaries of house are clear.These indicators show that matching imaging has met the requirements of high quality.

Conclusions
The matching imaging method combining attitude adjustment of the agile satellites and the scanning of the swing mirror is proposed to realize full coverage of the coastal zone.Based on the comprehensive consideration of the length, width, and curvature of the coastal zone, an adjustment curve of the roll angle of the agile satellites is planned to ensure scanning along the coastline during the satellite flight.Then, a scanning scheme perpendicular to the coastal zone through the swing mirror is proposed to enlarge the coverage of the wide coastal zone by satellite.The paraments of the swing mirror are calculated to ensure the overlap rate between the adjacent strips is greater than 10%.Finally, the velocity vector model influenced by roll angle and the scanning of the swing mirror is established to calculate the mismatch of velocity vector on the image plane.The adjustment of yaw angle, the tilt of detectors, and the adjustment of the imaging time of detectors are used to eliminate the mismatch and ensure the imaging quality.Experimental results show that the matching imaging method meets the requirement that the resolution is less than 1 pixel.It also meets the requirement the MTF of our image is larger than the minimum MTF of 0.026 of the human eyes.The matching imaging mode of the satellite attitude adjustment and the scanning of the swing mirror in this paper is suitable for all remote sensing observations of curved and extended targets.The velocity vector matching method can better eliminate the image quality degradation which is caused by the inconsistency between the velocity vector and the imaging direction of the detectors in all dynamic imaging.It has practical value for improving the imaging quality of dynamic remote sensing.

Figure 1 .
Figure 1.The full coverage technology along the coastal zone: (a) Attitude adjustment of the agile satellite along the coastline; (b) The swing mirror is used to realize full coverage of the coastal zone based on the agile satellite.

.B
fixed coordinate system, is established.Its origin is at the center of the earth.x G axis points to the intersection of the Prime Meridian and Equator.z G points to the North Pole.The positions of feature points in the coastal zone are ( , , ) The equations of roll angle and roll angle velocity which versus time in ( ) ⋅ are defined.The parameters of the agile satellite are ( ) the semi-major axis of the satellite.e represents the eccentricity of the sat- ellite orbit.i means the orbital inclination.N A is short for the right ascension of

Figure 1 .
Figure 1.The full coverage technology along the coastal zone: (a) Attitude adjustment of the agile satellite along the coastline; (b) The swing mirror is used to realize full coverage of the coastal zone based on the agile satellite.

Figure 2 .
Figure 2. Wide-swath whiskbroom imaging with the swing mirror: (a) Optical system o imaging along the coastal zone; (b) The imaging trajectory of the swing mirror.

Figure 2 .
Figure 2. Wide-swath whiskbroom imaging with the swing mirror: (a) Optical system of scanning imaging along the coastal zone; (b) The imaging trajectory of the swing mirror.

Figure 4 .
Figure 4. Roll angel and velocity vector: (a) The curve of roll angle; (b) The value of drift a 2367 s; (c) the distribution of velocity vectors on the image plane at 2367 s.

.Figure 5 .
Figure 5. Experimental principle and facilities: (a) The imaging principle of experiment; facilities of experiment.

Figure 4 .
Figure 4. Roll angel and velocity vector: (a) The curve of roll angle; (b) The value of drift angles at 2367 s; (c) the distribution of velocity vectors on the image plane at 2367 s.

.Figure 5 .
Figure 5. Experimental principle and facilities: (a) The imaging principle of experiment; (b) The facilities of experiment.Figure 5. Experimental principle and facilities: (a) The imaging principle of experiment; (b) The facilities of experiment.

Figure 5 .
Figure 5. Experimental principle and facilities: (a) The imaging principle of experiment; (b) The facilities of experiment.Figure 5. Experimental principle and facilities: (a) The imaging principle of experiment; (b) The facilities of experiment.

1 0
are the images of the edge targe.The chosen knife-edge of detector 1 is

Figure 6 .
Figure 6.Experimental results: (a) Resolution image of detector 1; (b) Resolution image of detector 2; (c) Edge image of detector 1; (d) Edge image of detector 2; (e) Color image of detector 1; (f) Color image of detector 2.

Table 2 .
The parameters in the experiment.