Geometric simulation and registration method for infrared hyperspectral atmospheric detector

Abstract

The invention relates to the technical field of satellite remote sensing, in particular to a geometric simulation and registration method for an infrared hyperspectral atmospheric detector, which comprises the following steps of: constructing a visual vector of a pixel in a reference ellipsoid coordinate system by acquiring detection time, ground longitude, ground latitude and satellite position of the pixel FOV of the HIRAS; simulating the imaging range of an HIRAS pixel on the ground based on the size of a pixel field angle by adopting a Rodrigues rotation formula, searching an MERSI pixel in the same platform by using a KD-tree, and respectively calculating the equivalent brightness temperature of the HIRAS pixel and the equivalent brightness temperature of an MERSI simulation image according to a spectrum convolution function and a blackbody radiation formula; and counting the root mean square error of the brightness temperature difference between the HIRAS pixel and the MERSI pixel through multi-row accumulation, and arranging the root mean square error according to an OR row time sequence to obtain the geometric positioning accuracy of the HIRAS in the along-rail and cross-rail directions. According to the invention, the rapid search of the pixel closest to the visual vector range is realized, and the accurate position determination of the registration error is realized.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of satellite remote sensing, and in particular to a geometric simulation and registration method for an infrared hyperspectral atmospheric sounder. BACKGROUND

[0002] The L1 data of various meteorological satellite infrared loads, such as a hyperspectral infrared atmospheric sounder (HIRAS) and a medium resolution spectral imager (MERSI), differ in geometric positioning, which makes it difficult to accurately match data from different sources under a unified spatial reference when performing fusion application or time series analysis. The L1 data is basic data that has been geolocated and radiometrically calibrated. This restricts the deep application efficiency of high-precision and consistent infrared data products in numerical weather prediction models and climate analysis and other meteorological core businesses.

[0003] The medium resolution image MERSI is a wide-spectrum imaging with a resolution of 250m and 1km. Generally, MERSI can be registered using the 250m resolution visible light band and higher resolution images, so that the geometric accuracy of the image can reach the sub-pixel level. The geometric relationship between the visible light and infrared channels of MERSI is stable, so the geometric accuracy of the infrared band is also high. The resolution of the FOV (Field of View, detection range) of one pixel of HIRAS is 14km, which is much larger than that of MERSI pixels, and the pixels of HIRAS are not continuously distributed, and there are gaps between the FOVs of each pixel. Therefore, the infrared channels of MERSI can be used to verify the geometric positioning accuracy of the infrared pixels of HIRAS.

[0004] The patent application No. 201910461383.5 discloses a positioning consistency verification method and device for infrared hyperspectral multi-pixel load, which discusses obtaining the deviation of low spatial resolution infrared hyperspectral multi-pixel load based on high spatial resolution infrared hyperspectral of the same platform. Some steps are similar to the patent, but the core process is different. The patent broadly gives the determination condition that the angle between the view vector of the second pixel in the ECEF coordinate system and the view vector of the first pixel is less than 1 / 2. In actual operation process, if each pixel of the high resolution image is determined, the calculation time cost is huge due to the large number of image pixels. The present project gives a fast search scheme and proposes several strategies to reduce the registration time. The patent application does not mention using the difference between the radiation brightness temperatures of the two matching samples as the basis for RMSE calculation, and many details are missing.

[0005] For example, the existing paper "Evaluation of HIRAS-II Positioning and Calibration Accuracy of FY-3E Same Platform Imager" discusses the method of positioning accuracy evaluation. However, this paper does not provide a method to accelerate pixel determination, nor does it mention how to use the view vector of the hyperspectral infrared atmospheric detector to simulate the process of pixel projection range on the ground.

[0006] Therefore, proposing a geometric simulation and registration method for infrared hyperspectral atmospheric detectors based on KD-tree search is of great practical significance. This method enables rapid search for the nearest pixel within the range of the view vector, thereby improving the registration speed and error calculation accuracy of the entire image. Summary of the Invention

[0007] To address the aforementioned issues, this invention provides a geometric simulation and registration method for an infrared hyperspectral atmospheric sounder. By constructing a KD-tree search method, it achieves rapid searching for the nearest pixel within the view vector range, thereby avoiding brute-force comparison of all pixels and reducing the time complexity from O(N) to an average of O(log N), thus improving the registration speed of the entire image. Furthermore, it sets sampling windows and step sizes. When the positioning accuracy of the original image is poor, a larger matching window and step size can be set to reduce registration time. When the accuracy after image correction is high, a smaller window and step size can be set to improve registration accuracy. The invention also clarifies that the simulated image and the original image are compared using temperature differences, and employs a multi-row overlay method to improve the accuracy of error calculation, achieving precise location determination of registration errors without the need for filtering the original data.

[0008] To achieve the above objectives, the embodiments of the present invention provide the following technical solutions:

[0009] This invention provides a geometric simulation and registration method for an infrared hyperspectral atmospheric sounder, comprising the following steps:

[0010] The detection time, ground longitude, ground latitude, and satellite position of a single pixel (FOV) from the Infrared Hyperspectral Atmospheric Sounder (HIRAS) are obtained, and the view vector of the pixel in the reference ellipsoidal coordinate system is constructed.

[0011] The imaging range of HIRAS pixels on the ground is simulated based on the pixel field of view using the Rodriguez rotation formula. The imaging range is formed by calculating the intersection points of multiple conical vectors with the ground to form a polygonal region.

[0012] Based on the imaging range, KD-tree is used to search for pixels of the same platform medium resolution imager (MERSI), including constructing an index of MERSI latitude and longitude data, reading block images, calculating MERSI pixel view vectors, and filtering pixels located within the HIRAS imaging range.

[0013] Based on the spectral convolution function and the blackbody radiation formula, the equivalent brightness temperature of the HIRAS pixel and the equivalent brightness temperature of the MERSI simulation image are calculated respectively.

[0014] The MERSI image is translated and registered by accumulating the root mean square error of the brightness temperature difference between HIRAS and MERSI pixels across multiple rows.

[0015] The root mean square error is arranged in a time sequence according to FOR (Field Of Regard) to obtain the geometric positioning accuracy of HIRAS in the along-track and across-track directions.

[0016] As a further aspect of the present invention, constructing the view vector of the pixel in the reference ellipsoidal coordinate system includes the following steps:

[0017] Based on the detection time of the obtained pixel , ground longitude Ground latitude and satellite position Convert latitude and longitude to ground point coordinates in the Earth-fixed system. , The three-dimensional coordinates in the Earth-fixed coordinate system are:

[0018] ;

[0019] ;

[0020] ;

[0021] in, ;

[0022] In the formula, Elevation represents the height of the observation point relative to the reference ellipsoid. The radius of curvature of the zonal circle is related to latitude. The semi-major axis of the reference ellipsoid is the WGS84 ellipsoid. The square of the first eccentricity of the ellipsoid; Indicates matrix transpose; Represents ground coordinates;

[0023] Calculate and normalize the view vector of the HIRAS FOV (Hyperspectral Infrared Atmospheric Sounder Field of View). In the formula, This represents the unit direction vector from the satellite to the ground point, used to represent the view vector of the HIRAS FOV; The satellite position is represented by the three-dimensional coordinate vector of the satellite in the Earth-fixed coordinate system. ; Coordinates of ground points under the Earth-Fixed System .

[0024] As a further aspect of the present invention, the Rodriguez rotation formula is used to simulate the imaging range of HIRAS pixels on the ground based on the pixel field of view size, including the following steps:

[0025] Based on the size of the pixel's field of view angle Construct any initial vector satisfy ;

[0026] Using Rodriguez's rotation formula Construct a conical vector that wraps around the view vector, where, The rotation angle is set to 10° to obtain 36 vectors; It is a conical vector;

[0027] All conical vectors, satellite positions, and ground intersection points together form a collinear equation, and the ground points satisfy the ellipsoidal equation of the reference ellipsoid.

[0028] The ground latitude and longitude corresponding to multiple conical vectors are solved sequentially, and the polygon formed by connecting multiple sets of latitude and longitude is represented as the imaging range of the FOV on the ground.

[0029] As a further aspect of the present invention, the collinearity equation is:

[0030] ;

[0031] The equation of the ellipsoid is:

[0032] In the formula, The coordinates of the ground intersection point represent the three-dimensional coordinates of the intersection point of the conical vector and the Earth's ellipsoid. This is the three-dimensional coordinate vector of the satellite's position; It is a scaling factor; With reference to the semi-major axis of the ellipsoid, The semi-minor axis of the ellipsoid Elevation.

[0033] As a further aspect of the present invention, using a KD-tree to search for pixels in a medium-resolution imager on the same platform includes the following steps:

[0034] Read MERSI imagery and latitude / longitude data captured within the same time period, build an index for all latitude / longitude data using a KD-tree, and find the nearest MERSI latitude / longitude coordinates to the target FOV view vector using a nearest neighbor search to obtain the image-side coordinates. ;

[0035] Based on the image coordinates, segmented images are read according to the reading range. The view vector of each MERSI pixel is calculated based on the satellite orbital position at the pixel imaging time. ;

[0036] Calculate view vectors and The angle between them, and filter those with an angle less than that of the HIRAS view vector. If the pixel is smaller than the HIRAS pixel FOV size If the MERSI pixel is within the HIRAS imaging range, calculate all pixels within the block image range and retain the pixels that meet the conditions.

[0037] As a further aspect of the present invention, the formula for calculating the reading range is:

[0038] ;

[0039] in, For the range to be read, This indicates the satellite's orbital altitude; for example, FY3F is 836km. This refers to the resolution of the corresponding MERSI band, such as the 1km resolution of CH22; This is the scaling factor, usually around 2, used to ensure a balance between memory overflow and speed. The size of the pixel's field of view angle.

[0040] As a further aspect of the present invention, the MERSI pixels of the MERSI image are: .

[0041] As a further aspect of the present invention, calculating the equivalent brightness temperature of HIRAS pixels and the equivalent brightness temperature of MERSI simulation images includes the following steps:

[0042] For HIRAS pixels, the spectral radiance of the HIRAS pixel is calculated using the continuous spectrum measured by HIRAS and combined with the MERSI spectral curve, employing a spectral convolution function: In the formula, The spectral radiance of a pixel, For the target at wavelength Spectral radiance at that location, wavelength Spectral response value at; and The band range represents the starting and ending wavelengths of the spectral band for integration.

[0043] The equivalent brightness temperature of HIRAS pixels is inverted using the blackbody radiation formula: In the formula, The equivalent brightness temperature of a HIRAS pixel. and The radiation constant is related to Planck's constant. The wave number corresponding to this band. The measured spectral radiance;

[0044] The temperature of each MERSI pixel is calculated using the blackbody radiation formula. Temperature inversion is then performed on all MERSI pixels within the HIRAS imaging range to calculate the equivalent brightness temperature of the MERSI simulation image. .

[0045] As a further aspect of the present invention, translational registration of MERSI images includes the following steps:

[0046] Comparison of HIRAS pixel equivalent brightness temperature Equivalent brightness temperature of MERSI pixels Differences as measured values ;

[0047] The window range set when MERSI pixels are translated. Internal translation, where, To register the window size, the step size is k, and the brightness temperature difference is calculated once for each translation.

[0048] The root mean square error (RMSE) of the brightness temperature difference at FOV is accumulated using a multi-row overlay method to obtain... The root mean square error matrix of size is calculated, and the position with the minimum value of the root mean square error matrix is ​​the optimal registration position.

[0049] Record the latitude and longitude information of the best registration position on MERSI, which is the actual latitude and longitude position of FOV in HIRAS at this time series. Compare it with the initial latitude and longitude to calculate the positioning error of the pixel.

[0050] Compared with existing technologies, the geometric simulation and registration method for infrared hyperspectral atmospheric detectors provided by this invention has the following advantages:

[0051] This invention improves the efficiency and accuracy of simulation verification. By employing a KD-tree search method for nearest pixel search, it achieves rapid searching for the nearest pixel within the view vector range. Based on the search results, it finds effective registration pixels within the field of view, improving the registration speed of the entire scene image. By introducing configurable sampling window size and step size, a dynamic balance between registration accuracy and efficiency is achieved. When the positioning accuracy of the original image is poor, a larger window and step size can be set for rapid registration; when the positioning accuracy is high, a smaller window and step size are used to improve detail accuracy. This invention also clarifies that the comparison between the simulated image and the original image is based on temperature difference. Brightness temperature inversion makes the registration results more reliable. By using multi-row superposition of the root mean square error of the brightness temperature difference, the precise location of the registration error is determined. Precise location of the registration error can be achieved without filtering the original data, realizing high-precision simulation of the HIRAS pixel imaging range.

[0052] These or other aspects of the invention will become more apparent from the following description of embodiments. It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only, and are not intended to limit the invention. Attached Figure Description

[0053] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. In the drawings:

[0054] Fig. 1 This is a flowchart of a geometric simulation and registration method for an infrared hyperspectral atmospheric detector according to the present invention.

[0055] Fig. 2 This is a schematic diagram showing the overlap between the FOV conical surface of the infrared hyperspectral atmospheric detector and the mid-resolution imager in the geometric simulation and registration method of the infrared hyperspectral atmospheric detector of the present invention.

[0056] Fig. 3 This is a schematic diagram showing the correspondence between the simulated HIRAS imaging range and the MERSI image obtained in the geometric simulation and registration method of an infrared hyperspectral atmospheric detector according to the present invention.

[0057] Fig. 4 This is a schematic diagram of the calculated position and registration position in the geometric simulation and registration method of an infrared hyperspectral atmospheric detector according to the present invention. Detailed Implementation

[0058] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0059] The technical solutions in the exemplary embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described exemplary embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0060] See Figs. 1 to 4 As shown in the figure, embodiments of this application also provide a geometric simulation and registration method for an infrared hyperspectral atmospheric sounder, including the following steps:

[0061] Step S10: Obtain the detection time, ground longitude, ground latitude and satellite position of a pixel (FOV) of the Infrared Hyperspectral Atmospheric Sounder (HIRAS), and construct the view vector of the pixel in the reference ellipsoidal coordinate system.

[0062] This step involves constructing the view vector of the pixel in the reference ellipsoidal coordinate system, including the following steps:

[0063] Step S101: Based on the obtained detection time of the pixel , ground longitude Ground latitude and satellite position Convert latitude and longitude to ground point coordinates in the Earth-fixed system. , The three-dimensional coordinates in the Earth-fixed coordinate system are:

[0064] ;

[0065] ;

[0066] ;

[0067] in, ;

[0068] In the formula, Elevation represents the height of the observation point relative to the reference ellipsoid. The radius of curvature of the zonal circle is related to latitude. The semi-major axis of the reference ellipsoid is the WGS84 ellipsoid. The square of the first eccentricity of the ellipsoid; Indicates matrix transpose; Represents ground coordinates;

[0069] Step S102: Calculate the view vector of HIRAS FOV (Hyperspectral Infrared Atmospheric Sounder Field of View) and normalize it.

[0070] Among them, view vector In the formula, This represents the unit direction vector from the satellite to the ground point, used to represent the view vector of the HIRAS FOV; The satellite position is represented by the three-dimensional coordinate vector of the satellite in the Earth-fixed coordinate system. ; Coordinates of ground points under the Earth-Fixed System .

[0071] Step S20: Using the Rodriguez rotation formula, the imaging range of HIRAS pixels on the ground is simulated based on the pixel field of view size. The imaging range is formed by calculating the intersection points of multiple conical vectors with the ground to form a polygonal region.

[0072] In this step, the Rodriguez rotation formula is used to simulate the imaging range of HIRAS pixels on the ground based on the pixel field of view size, including the following steps:

[0073] Step S201: Based on the pixel field of view angle size Construct any initial vector satisfy .

[0074] Step S202: The root is rotated using the Rodrigues formula. Construct a conical vector that wraps around the view vector, where, The rotation angle is set to 10° to obtain 36 vectors; It is a conical vector.

[0075] Among them, the initial vector about the axis of rotation Rotation angle To obtain the next vector, when the interval is 10°, 36 conical vectors can be obtained.

[0076] Step S203: All conical vectors, satellite positions, and ground intersection points together form a collinear equation, and the ground points satisfy the ellipsoid equation of the reference ellipsoid.

[0077] The collinear equation is:

[0078] ;

[0079] The equation of the ellipsoid is:

[0080] In the formula, The coordinates of the ground intersection point represent the three-dimensional coordinates of the intersection point of the conical vector and the Earth's ellipsoid. This is the three-dimensional coordinate vector of the satellite's position; It is a scaling factor; With reference to the semi-major axis of the ellipsoid, The semi-minor axis of the ellipsoid Elevation.

[0081] Step S204: Solve the ground latitude and longitude corresponding to multiple conical vectors in sequence, and connect the multiple sets of latitude and longitude in sequence to form a polygon. The polygon represents the imaging range of the FOV on the ground.

[0082] By combining the above collinearity equation and ellipsoid equation formula, the ground latitude and longitude corresponding to multiple conical vectors can be solved sequentially. The polygon formed by connecting multiple sets of latitude and longitude can be represented as the imaging range of the FOV on the ground, thereby realizing the geometric positioning simulation of HIRAS pixels.

[0083] Step S30: Based on the imaging range, use KD-tree to search for pixels of the same platform Medium Resolution Imager (MERSI), including constructing an index of MERSI latitude and longitude data, reading block images, calculating MERSI pixel view vectors, and filtering pixels located within the HIRAS imaging range.

[0084] This step involves using a KD-tree to search for pixels in a medium-resolution imager on the same platform, including the following steps:

[0085] Step S301: Read MERSI images and latitude / longitude data captured in the same time period, build an index for all latitude / longitude data using a KD-tree, and find the MERSI latitude / longitude closest to the target FOV view vector using nearest neighbor search to obtain the image coordinates. .

[0086] Step S302: Read the segmented image according to the image coordinates within the reading range, wherein the calculation formula for the reading range is:

[0087] ;

[0088] in, For the range to be read, This indicates the satellite's orbital altitude; for example, FY3F is 836km. This refers to the resolution of the corresponding MERSI band, such as the 1km resolution of CH22; This is the scaling factor, usually around 2, used to ensure a balance between memory overflow and speed. The size of the pixel's field of view angle.

[0089] Step S303: Calculate the view vector of the MERSI pixel based on the satellite orbit position at the pixel imaging time. .

[0090] Step S304: Calculate the view vector and The angle between them, and filter those with an angle less than that of the HIRAS view vector. If the pixel is smaller than the HIRAS pixel FOV size If the MERSI pixel is within the HIRAS imaging range, calculate all pixels within the block image range and retain the pixels that meet the conditions.

[0091] Wherein, the MERSI pixels of the MERSI image are .

[0092] Step S40: Calculate the equivalent brightness temperature of the HIRAS pixel and the equivalent brightness temperature of the MERSI simulation image based on the spectral convolution function and the blackbody radiation formula.

[0093] This step involves calculating the equivalent brightness temperature of HIRAS pixels and the equivalent brightness temperature of the MERSI simulation image, including the following steps:

[0094] Step S401: For HIRAS pixels, using the continuous spectrum measured by HIRAS and combined with the MERSI spectral curve, calculate the spectral radiance of the HIRAS pixel using the spectral convolution function: In the formula, The spectral radiance of a pixel, For the target at wavelength Spectral radiance at that location, wavelength Spectral response value at; and The band range represents the starting and ending wavelengths of the spectral band for integration.

[0095] Step S402: Invert the equivalent brightness temperature of the HIRAS pixel according to the blackbody radiation formula: In the formula, The equivalent brightness temperature of a HIRAS pixel. and The radiation constant is related to Planck's constant. The wave number corresponding to this band. The measured spectral radiance;

[0096] Step S403: Calculate the temperature of each MERSI pixel using the blackbody radiation formula, perform temperature inversion on all MERSI pixels within the HIRAS imaging range, and calculate the equivalent brightness temperature of the MERSI simulation image. .

[0097] Step S50: Perform translation registration on the MERSI image by accumulating the root mean square error of the brightness temperature difference between HIRAS pixels and MERSI pixels across multiple rows.

[0098] This step involves translational registration of the MERSI image, including the following steps:

[0099] Step S501: Compare the equivalent brightness temperature of HIRAS pixels. Equivalent brightness temperature of MERSI pixels Differences as measured values ;

[0100] Step S502: Set the window range for MERSI pixels during translation. Internal translation, where, To register the window size, the step size is k, and the brightness temperature difference is calculated once for each translation.

[0101] Step S503: Accumulate and statistically analyze the root mean square error (RMSE) of the FOV brightness temperature difference using a multi-row overlay method to obtain... The root mean square error matrix of size is calculated, and the position with the minimum value of the root mean square error matrix is ​​the optimal registration position.

[0102] Step S504: Record the latitude and longitude information of the best registration position on MERSI, which is the actual latitude and longitude position of FOV in HIRAS at this time series. Compare it with the initial latitude and longitude to calculate the positioning error of the pixel.

[0103] Step S60: Arrange the root mean square error in a time sequence according to FOR (Field Of Regard) to obtain the geometric positioning accuracy of HIRAS in the along-track and across-track directions.

[0104] In this step, arranging the error according to the FOR line sequence allows us to obtain the geometric positioning accuracy of HIRAS in the along-rail and across-rail directions, which is beneficial for subsequent positioning accuracy verification and correction.

[0105] This invention improves the efficiency and accuracy of simulation verification. By employing a KD-tree search method for nearest pixel search, it achieves rapid searching for the nearest pixel within the view vector range. Based on the search results, it finds effective registration pixels within the field of view, improving the registration speed of the entire scene image. By introducing configurable sampling window size and step size, a dynamic balance between registration accuracy and efficiency is achieved. When the positioning accuracy of the original image is poor, a larger window and step size can be set for rapid registration; when the positioning accuracy is high, a smaller window and step size are used to improve detail accuracy. This invention also clarifies that the comparison between the simulated image and the original image is based on temperature difference. Brightness temperature inversion makes the registration results more reliable. By using multi-row superposition of the root mean square error of the brightness temperature difference, the precise location of the registration error is determined. Precise location of the registration error can be achieved without filtering the original data, realizing high-precision simulation of the HIRAS pixel imaging range.

[0106] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A geometric simulation and registration method for an infrared hyperspectral atmospheric sounder, characterized in that, Includes the following steps: Obtain the detection time, ground longitude, ground latitude, and satellite position of a HIRAS pixel's FOV, and construct the view vector of the pixel in the reference ellipsoidal coordinate system; The imaging range of HIRAS pixels on the ground is simulated based on the pixel field of view using the Rodriguez rotation formula. The imaging range is formed by calculating the intersection points of multiple conical vectors with the ground to form a polygonal region. Based on the imaging range, KD-tree is used to search for MERSI pixels on the same platform, including constructing an index of MERSI latitude and longitude data, reading block images, calculating MERSI pixel view vectors, and filtering pixels located within the HIRAS imaging range. Based on the spectral convolution function and the blackbody radiation formula, the equivalent brightness temperature of the HIRAS pixel and the equivalent brightness temperature of the MERSI simulation image are calculated respectively. The MERSI image is translated and registered by accumulating the root mean square error of the brightness temperature difference between HIRAS and MERSI pixels across multiple rows. The root mean square error is arranged in a FOR row time sequence to obtain the geometric positioning accuracy of HIRAS in the along-rail and across-rail directions.

2. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 1, characterized in that, Constructing the view vector of the pixel in the reference ellipsoidal coordinate system includes the following steps: Based on the detection time of the obtained pixel , ground longitude Ground latitude and satellite position Convert latitude and longitude to ground point coordinates in the Earth-fixed system. , The three-dimensional coordinates in the Earth-fixed coordinate system are: ； ； ； in, ; In the formula, Elevation represents the height of the observation point relative to the reference ellipsoid. The radius of curvature of the circumpolar region is related to latitude. The semi-major axis of the reference ellipsoid is the WGS84 ellipsoid. The square of the first eccentricity of the ellipsoid; Calculate the view vector of the HIRAS FOV and normalize it, where the view vector... In the formula, This represents the unit direction vector from the satellite to the ground point, used to represent the view vector of the HIRAS FOV; The satellite position is represented by the three-dimensional coordinate vector of the satellite in the Earth-fixed coordinate system. ; Coordinates of ground points in the Earth-Fixed System .

3. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 2, characterized in that, Using the Rodriguez rotation formula, the imaging range of HIRAS pixels on the ground is simulated based on the pixel field of view size, including the following steps: Based on the size of the pixel's field of view angle Construct any initial vector satisfy ; Using Rodriguez's rotation formula Construct a conical vector that wraps around the view vector, where, The rotation angle is set to 10° to obtain 36 vectors; It is a conical vector; All conical vectors, satellite positions, and ground intersection points together form a collinear equation, and the ground points satisfy the ellipsoidal equation of the reference ellipsoid. The ground latitude and longitude corresponding to multiple conical vectors are solved sequentially, and the polygon formed by connecting multiple sets of latitude and longitude is represented as the imaging range of the FOV on the ground.

4. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 3, characterized in that, The collinear equation is: ； The equation of the ellipsoid is: In the formula, The coordinates of the ground intersection point represent the three-dimensional coordinates of the intersection point of the conical vector and the Earth's ellipsoid. This is the three-dimensional coordinate vector of the satellite's position; It is a scaling factor; With reference to the semi-major axis of the ellipsoid, The semi-minor axis of the ellipsoid Elevation.

5. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 4, characterized in that, Using a KD-tree to search for pixels in a medium-resolution imager on the same platform includes the following steps: Read MERSI imagery and latitude / longitude data captured within the same time period, build an index for all latitude / longitude data using a KD-tree, and find the nearest MERSI latitude / longitude coordinates to the target FOV view vector using a nearest neighbor search to obtain the image-side coordinates. ; Based on the image coordinates, segmented images are read according to the reading range. The view vector of each MERSI pixel is calculated based on the satellite orbital position at the pixel imaging time. ; Calculate view vectors and The angle between them, and filter those with an angle less than that of the HIRAS view vector. If the pixel is smaller than the HIRAS pixel FOV size If the MERSI pixel is within the HIRAS imaging range, calculate all pixels within the block image range and retain the pixels that meet the conditions.

6. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 5, characterized in that, The formula for calculating the reading range is: ； in, For the range to be read, Indicates the satellite's orbital altitude; This represents the resolution of the corresponding MERSI band; This refers to the scaling ratio; The size of the pixel's field of view angle.

7. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 6, characterized in that, The MERSI pixels of the MERSI image are .

8. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 7, characterized in that, Calculating the equivalent brightness temperature of HIRAS pixels and the equivalent brightness temperature of the MERSI simulation image involves the following steps: For HIRAS pixels, the spectral radiance of the HIRAS pixel is calculated using the continuous spectrum measured by HIRAS and combined with the MERSI spectral curve, employing a spectral convolution function: In the formula, The spectral radiance of a pixel, For the target at wavelength Spectral radiance at that location, wavelength Spectral response value at; and The band range represents the starting and ending wavelengths of the spectral band for integration. The equivalent brightness temperature of HIRAS pixels is inverted using the blackbody radiation formula: In the formula, The equivalent brightness temperature of a HIRAS pixel. and The radiation constant is related to Planck's constant. The wave number corresponding to this band. The measured spectral radiance; The temperature of each MERSI pixel is calculated using the blackbody radiation formula. Temperature inversion is then performed on all MERSI pixels within the HIRAS imaging range to calculate the equivalent brightness temperature of the MERSI simulation image. .

9. The geometric simulation and registration method for an infrared hyperspectral atmospheric sounder as described in claim 8, characterized in that, Translational registration of MERSI images includes the following steps: Comparison of HIRAS pixel equivalent brightness temperature Equivalent brightness temperature of MERSI pixels Differences as measured values ; The window range set when MERSI pixels are translated. Internal translation, where, To register the window size, the step size is k, and the brightness temperature difference is calculated once for each translation. The root mean square error of the FOV brightness temperature difference is accumulated and statistically analyzed using a multi-row overlay method. The root mean square error matrix of size is calculated, and the position with the minimum value of the root mean square error matrix is ​​the optimal registration position. Record the latitude and longitude information of the best registration position on MERSI, which is the actual latitude and longitude position of FOV in HIRAS at this time series. Compare it with the initial latitude and longitude to calculate the positioning error of the pixel.

Images