A method and device for precise registration of sentinel satellite data

Abstract

The application provides a kind of method and device for the fine registration of sentinel satellite data.The method comprises the following steps: obtaining the upper and lower band data of the main image after coarse registration and the upper and lower band data of the auxiliary image;according to the upper and lower band data of the main image and the auxiliary image, the coherence and differential phase of each upper and lower band are calculated iteratively;the coherence and differential phase of the upper and lower band are processed by multi-view, and the differential phase of enhanced spectral diversity is determined based on the coherence and differential phase after multi-view processing;according to the differential phase of enhanced spectral diversity, the pulse repetition frequency and the change value of adjacent burst doppler center frequency, the residual offset in azimuth direction of the enhanced spectral diversity is calculated to complete the fine registration of the sentinel satellite data.The present scheme can improve the registration accuracy of the sentinel satellite data.

Description

Technical Field

[0001] This invention relates to the field of signal processing technology, and in particular to a method and apparatus for precise registration of Sentinel satellite data. Background Technology

[0002] The Sentinel satellite's IW mode employs a progressive scan SAR (TOPS) mode for ground target observation and imaging. During operation, the beam pointing rotates periodically in both the azimuth and range directions, enabling wide-swath imaging of three sub-mapped zones and 27 bursts (image slices). The azimuth beam rotation causes a frequency shift in the Doppler center frequency of individual burst data, reaching a maximum of 5200Hz. To avoid phase jumps in the interferograms during subsequent time-series processing, fine registration is required on top of coarse registration, with a registration accuracy of at least one-thousandth of a pixel to meet the registration requirements.

[0003] In related technologies, enhanced spectral diversity is generally used to achieve fine registration when processing Sentinel satellite data. Enhanced spectral diversity is based on the principle of spectral diversity, but it utilizes the overlapping region between adjacent bursts as upper and lower frequency bands, and differentially obtains the ESD phase to achieve fine registration and compensate for the residual azimuth offset of coarse registration. However, the fine registration method in related technologies performs poorly when processing a large number of adjacent bursts with incoherent images, resulting in poor registration accuracy.

[0004] Therefore, there is an urgent need for a new method and device for precise registration of Sentinel satellite data. Summary of the Invention

[0005] To address the problem that traditional fine registration methods perform poorly when processing a large number of adjacent burst-incoherent images, resulting in low registration accuracy, this invention provides a fine registration method and apparatus for Sentinel satellite data.

[0006] In a first aspect, embodiments of the present invention provide a method for fine registration of Sentinel satellite data, the method comprising:

[0007] Obtain the upper and lower frequency band data of the main image and the auxiliary image after coarse registration;

[0008] Based on the upper and lower frequency band data of the main image and the auxiliary image, the coherence and differential phase of each upper and lower frequency band are calculated traversally.

[0009] The coherence and differential phase of the upper and lower frequency bands are subjected to multi-view processing, and the differential phase of the enhanced spectrum diversity is determined based on the coherence and differential phase after multi-view processing.

[0010] Based on the differential phase, pulse repetition frequency, and changes in the center frequencies of adjacent burst Doppler signals of the enhanced spectral diversity, the azimuth residual offset of the enhanced spectral diversity is calculated to complete the fine registration of the Sentinel satellite data.

[0011] Secondly, embodiments of the present invention also provide a fine registration device for Sentinel satellite data, the device comprising:

[0012] The acquisition unit is used to acquire the upper and lower frequency band data of the main image and the upper and lower frequency band data of the auxiliary image after coarse registration.

[0013] The first calculation unit is used to calculate the coherence and differential phase of each of the upper and lower frequency bands based on the upper and lower frequency band data of the main image and the auxiliary image.

[0014] The determining unit is used to perform multi-view processing on the coherence and differential phase of the upper and lower frequency bands, and determine the differential phase of the enhanced spectrum diversity based on the coherence and differential phase of the upper and lower frequency bands after multi-view processing.

[0015] The second calculation unit is used to calculate the azimuth residual offset of the enhanced spectrum diversity based on the differential phase, pulse repetition frequency, and changes in the center frequencies of adjacent burst Doppler signals, so as to complete the fine registration of the Sentinel satellite data.

[0016] Thirdly, embodiments of the present invention also provide a computing device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, it implements the method described in any embodiment of this specification.

[0017] Fourthly, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to perform the methods described in any embodiment of this specification.

[0018] This invention provides a method and apparatus for fine registration of Sentinel satellite data. First, using the upper and lower frequency band data of the main image and the auxiliary image, the coherence and differential phase of the upper and lower frequency bands are determined. Then, multi-view processing is performed on the coherence and differential phase of each upper and lower frequency band to suppress noise. Based on the correspondence between the coherence and differential phase after multi-view processing, the differential phase of the enhanced spectrum diversity is determined. This can effectively improve the estimation accuracy of the azimuth residual offset of the enhanced spectrum diversity, thereby improving the registration accuracy of Sentinel satellite data. Attached Figure Description

[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0020] Figure 1 This is a flowchart of a method for fine registration of Sentinel satellite data provided in an embodiment of the present invention;

[0021] Figure 2 It is a differential phase distribution map of the upper and lower frequency band interferograms of adjacent burst overlap regions after fine registration of Sentinel satellite data using the traditional enhanced spectrum diversity method;

[0022] Figure 3 This is a differential phase distribution map of the upper and lower frequency band interferograms of adjacent burst overlapping regions after the Sentinel satellite data is finely registered using the method in the embodiments of the present invention.

[0023] Figure 4 This is a hardware architecture diagram of a computing device provided in an embodiment of the present invention;

[0024] Figure 5 This is a structural diagram of a fine registration device for Sentinel satellite data provided in an embodiment of the present invention. Detailed Implementation

[0025] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0026] The following describes the specific implementation of the above concept.

[0027] Please refer to Figure 1 This invention provides a method for fine registration of Sentinel satellite data, the method comprising:

[0028] Step 100: Obtain the upper and lower frequency band data of the main image and the upper and lower frequency band data of the auxiliary image after coarse registration;

[0029] Step 102: Based on the upper and lower frequency band data of the main image and the auxiliary image, calculate the coherence and differential phase of each upper and lower frequency band.

[0030] Step 104: Perform multi-view processing on the coherence and differential phase of the upper and lower frequency bands, and determine the differential phase of the enhanced spectrum diversity based on the coherence and differential phase after multi-view processing.

[0031] Step 106: Calculate the azimuth residual offset of the enhanced spectrum diversity based on the differential phase, pulse repetition frequency, and changes in the center frequencies of adjacent burst Doppler signals, in order to complete the fine registration of the Sentinel satellite data.

[0032] In this embodiment of the invention, the coherence and differential phase of the upper and lower frequency bands are first determined using the upper and lower frequency band data of the main image and the auxiliary image. Then, multi-view processing is performed on the coherence and differential phase of each upper and lower frequency band to suppress noise. Based on the correspondence between the coherence and differential phase after multi-view processing, the differential phase of the enhanced spectrum diversity is determined. This can effectively improve the estimation accuracy of the azimuth residual offset of the enhanced spectrum diversity, thereby improving the registration accuracy of the Sentinel satellite data.

[0033] For step 100:

[0034] In some implementations, the upper and lower frequency band data of the main image includes data of each adjacent burst overlapping region in the main image, and the upper and lower frequency band data of the auxiliary image includes data of each adjacent burst overlapping region in the auxiliary image.

[0035] It should be noted that in this embodiment, before acquiring the upper and lower frequency band data of the main image and the auxiliary image, coarse registration of the main image and the auxiliary image is required. The coarse registration method can be geometric registration.

[0036] Regarding step 102:

[0037] In some implementations, step 102 includes:

[0038] Based on the upper frequency band data of the main image and the upper frequency band data of the auxiliary image, the interference data of the upper and lower frequency bands are calculated respectively.

[0039] Based on the interference data of the upper and lower frequency bands, the coherence and differential phase of the upper and lower frequency bands are calculated respectively.

[0040] Unlike traditional enhanced spectral diversity methods that typically estimate the azimuth residual offset of data from each adjacent burst overlap region in the primary and secondary images, this invention treats adjacent burst overlap regions in the primary and secondary images as a whole. Specifically, it uses the data from each adjacent burst overlap region in the primary image and the data from each adjacent burst overlap region in the secondary image as upper and lower frequency band data. The differential phase and coherence of the upper and lower frequency band interferometric data are calculated iteratively. This effectively avoids the virtual coherence phenomenon between adjacent burst overlap regions, allowing the upper and lower frequency band interferometric data to be better suited for subsequent coherence and multi-view processing algorithms, thereby effectively improving the estimation accuracy of the azimuth offset in enhanced spectral diversity.

[0041] In some implementations, the interference data of the upper and lower frequency bands are calculated using the following formula:

[0042]

[0043]

[0044] In the formula, I up For the interferometric data of the upper frequency band, I down For the interferometric data of the lower frequency band, M overlapUp S represents the upper frequency band data in the main image. overlapUp M is the upper frequency band data in the auxiliary image. overlapDown S represents the lower frequency band data in the main image. overlapDown This refers to the lower frequency band data in the auxiliary image.

[0045] In some implementations, the coherence of the upper and lower frequency bands and the differential phase of the upper and lower frequency bands are calculated using the following formulas:

[0046]

[0047]

[0048] In the formula, γ represents the coherence of the upper and lower frequency bands, M represents the size of the coherence window in the azimuth direction, N represents the size of the coherence window range image, and I represents the size of the range image. up For the interferometric data of the upper frequency band, I down Let i be the coordinate number of the pixel in the azimuth coherence window, j be the coordinate number of the pixel in the range coherence window, φ be the differential phase of the upper and lower frequency bands, and arg{·} be the phase extraction operation.

[0049] Regarding step 104:

[0050] In the process of estimating the azimuth residual offset of an image directly through differential interferometric phase of interferometric data in the upper and lower frequency bands of adjacent burst overlap regions in enhanced spectral diversity, noise in the interferometric data can easily interfere with the estimation accuracy of the azimuth residual offset. In this embodiment of the invention, after obtaining the coherence and differential phase of all adjacent burst overlap regions in the main image and the auxiliary image in step 102, multi-view processing is then performed on the overall upper and lower frequency band coherence and the overall upper and lower frequency band differential phase. This not only improves the processing efficiency but also effectively suppresses noise interference in the interferometric data, thereby ensuring the accuracy of the azimuth residual offset calculation in enhanced spectral diversity.

[0051] In some implementations, the coherence after multi-visualization processing and the differential phase after multi-visualization processing are calculated using the following formulas:

[0052]

[0053]

[0054] In the formula, γ mul The coherence of the upper and lower frequency bands after the multi-viewing processing is denoted by A, where A is the multi-view factor in the azimuth direction of the coherence window, R is the multi-view factor in the range direction of the coherence window, γ is the coherence of the upper and lower frequency bands, i is the coordinate number of the pixel in the azimuth multi-view window, j is the coordinate number of the pixel in the range multi-view window, and φ is the coordinate number of the pixel in the range multi-view window. mul The differential phase of the upper and lower frequency bands after multi-visualization processing, φ is the differential phase of the upper and lower frequency bands.

[0055] In some implementations, the differential phase of the enhanced spectral diversity is determined as follows:

[0056] After multi-visualization processing, the differential phase corresponding to coherence less than a preset threshold is set to zero, and the differential phase after setting to zero and the differential phase corresponding to coherence greater than the preset threshold are taken as the differential phase set; wherein, the preset threshold is 0.5;

[0057] The mode of the differential phase set is selected as the differential phase of the enhanced spectral diversity.

[0058] In this embodiment of the invention, when calculating the azimuth residual offset for enhanced spectral diversity, a preset coherence threshold is set to filter the differential phase. Pixels in the primary and secondary images with coherence below 0.5 in adjacent burst overlap areas are set to zero to prevent them from affecting the differential phase decision of enhanced spectral diversity. Pixels in the primary and secondary images with coherence above 0.5 in adjacent burst overlap areas are used for enhanced spectral diversity differential phase estimation, thereby deriving the azimuth residual offset. This avoids introducing some low-coherence pixels into the differential phase estimation, thus simplifying the estimation process and effectively improving the registration accuracy of Sentinel satellite data.

[0059] Regarding step 106:

[0060] In some implementations, the azimuth residual offset of the enhanced spectral diversity is calculated using the following formula:

[0061]

[0062] In the formula, Δaz is the azimuth residual offset of the enhanced spectral diversity, and φ ESD Δf is the differential phase of the enhanced spectrum diversity, Δf is the change in the center frequency of the adjacent burst Doppler, and PRF is the pulse repetition frequency in the azimuth direction.

[0063] For example, in this embodiment of the invention, based on the completion of geometric coarse registration, the above method is used to perform fine registration on the sub-mapped strips of the real sentry data. The relevant parameters used in the fine registration process are shown in Table 1:

[0064] Table 1

[0065]

[0066] Figure 2 This describes the differential phase distribution of the upper and lower frequency band interferograms in the overlapping region of adjacent bursts after fine registration of Sentinel satellite data using traditional enhanced spectral diversity methods. Figure 3 Table 2 shows the differential phase distribution of the upper and lower frequency band interferograms in the overlapping region of adjacent bursts after fine registration of sentinel satellite data according to the method in this embodiment of the invention. Figure 2 and Figure 3 The average phase value.

[0067] Table 2

[0068]

[0069] Combination Figure 2 , Figure 3As shown in Table 2, compared to the traditional enhanced spectral diversity method, the mean differential phase of the upper and lower frequency band interferograms obtained by the fine registration method in this embodiment of the invention is reduced by approximately 12 times, meaning the phase error is closer to 0. In summary, the method in this embodiment of the invention can effectively improve registration accuracy when performing fine registration of Sentinel satellite data.

[0070] like Figure 4 , Figure 5 As shown, this embodiment of the invention provides a fine registration device for Sentinel satellite data. The device embodiment can be implemented through software, hardware, or a combination of both. From a hardware perspective, as... Figure 4 The diagram shown is a hardware architecture diagram of a computing device housing a fine registration device for Sentinel satellite data provided in an embodiment of the present invention. (Except for...) Figure 4 In addition to the processor, memory, network interface, and non-volatile memory shown, the computing device in the embodiment may also include other hardware, such as a forwarding chip responsible for processing packets. Taking software implementation as an example, such as... Figure 5 As shown, a device in a logical sense is formed by the CPU of its computing device reading the corresponding computer program from non-volatile memory into memory for execution. This embodiment provides a fine registration device for Sentinel satellite data, the device comprising:

[0071] The acquisition unit 501 is used to acquire the upper and lower frequency band data of the main image and the upper and lower frequency band data of the auxiliary image after coarse registration.

[0072] The first calculation unit 502 is used to calculate the coherence and differential phase of each upper and lower frequency band based on the upper and lower frequency band data of the main image and the auxiliary image.

[0073] The determining unit 503 is used to perform multi-view processing on the coherence and differential phase of the upper and lower frequency bands, and determine the differential phase of the enhanced spectrum diversity based on the coherence and differential phase of the upper and lower frequency bands after multi-view processing.

[0074] The second calculation unit 504 is used to calculate the azimuth residual offset of the enhanced spectrum diversity based on the differential phase, pulse repetition frequency and the change value of the adjacent burst Doppler center frequency of the enhanced spectrum diversity, so as to complete the fine registration of the Sentinel satellite data.

[0075] In one embodiment of the present invention, the upper and lower frequency band data of the main image in the acquisition unit 501 includes the data of each adjacent burst overlapping region in the main image, and the upper and lower frequency band data of the auxiliary image includes the data of each adjacent burst overlapping region in the auxiliary image.

[0076] In one embodiment of the present invention, the first computing unit 502 is configured to perform the following operations:

[0077] Based on the upper frequency band data of the main image and the upper frequency band data of the auxiliary image, the interference data of the upper and lower frequency bands are calculated respectively.

[0078] Based on the interference data of the upper and lower frequency bands, the coherence and differential phase of the upper and lower frequency bands are calculated respectively.

[0079] In one embodiment of the present invention, in the first calculation unit 502, the interference data of the upper and lower frequency bands are calculated by the following formula:

[0080]

[0081]

[0082] In the formula, I up For the interferometric data of the upper frequency band, I down For the interferometric data of the lower frequency band, M overlapUp S represents the upper frequency band data in the main image. overlapUp M is the upper frequency band data in the auxiliary image. overlapDown S represents the lower frequency band data in the main image. overlapDown This refers to the lower frequency band data in the auxiliary image.

[0083] In one embodiment of the present invention, in the first calculation unit 502, the coherence of the upper and lower frequency bands and the differential phase of the upper and lower frequency bands are calculated by the following formulas:

[0084]

[0085]

[0086] In the formula, γ represents the coherence of the upper and lower frequency bands, M represents the size of the coherence window in the azimuth direction, N represents the size of the coherence window range image, and I represents the size of the range image. up For the interferometric data of the upper frequency band, I down Let i be the coordinate number of the pixel in the azimuth coherence window, j be the coordinate number of the pixel in the range coherence window, φ be the differential phase of the upper and lower frequency bands, and arg{·} be the phase extraction operation.

[0087] In one embodiment of the present invention, in the determining unit 503, the coherence after multi-visualization processing and the differential phase after multi-visualization processing are calculated by the following formulas:

[0088]

[0089]

[0090] In the formula, γ mul The coherence of the upper and lower frequency bands after the multi-viewing processing is denoted by A, where A is the multi-view factor in the azimuth direction of the coherence window, R is the multi-view factor in the range direction of the coherence window, γ is the coherence of the upper and lower frequency bands, i is the coordinate number of the pixel in the azimuth multi-view window, j is the coordinate number of the pixel in the range multi-view window, and φ is the coordinate number of the pixel in the range multi-view window. mul The differential phase of the upper and lower frequency bands after multi-visualization processing, φ is the differential phase of the upper and lower frequency bands.

[0091] In one embodiment of the present invention, in the determining unit 503, the differential phase of the enhanced spectral diversity is determined in the following manner:

[0092] After multi-visualization processing, the differential phase corresponding to coherence less than a preset threshold is set to zero, and the differential phase after setting to zero and the differential phase corresponding to coherence greater than the preset threshold are taken as the differential phase set; wherein, the preset threshold is 0.5;

[0093] The mode of the differential phase set is selected as the differential phase of the enhanced spectral diversity.

[0094] In one embodiment of the present invention, in the second calculation unit 504, the azimuth residual offset of the enhanced spectral diversity is calculated by the following formula:

[0095]

[0096] In the formula, Δaz is the azimuth residual offset of the enhanced spectral diversity, and φ ESD Δf is the differential phase of the enhanced spectrum diversity, Δf is the change in the center frequency of the adjacent burst Doppler, and PRF is the pulse repetition frequency in the azimuth direction.

[0097] It is understood that the structures illustrated in the embodiments of the present invention do not constitute a specific limitation on a fine registration device for Sentinel satellite data. In other embodiments of the present invention, a fine registration device for Sentinel satellite data may include more or fewer components than illustrated, or combine some components, split some components, or have different component arrangements. The illustrated components may be implemented in hardware, software, or a combination of software and hardware.

[0098] The information interaction and execution process between the modules in the above-mentioned device are based on the same concept as the method embodiment of the present invention, and the specific details can be found in the description of the method embodiment of the present invention, and will not be repeated here.

[0099] This invention also provides a computing device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, it implements a method for fine registration of sentinel satellite data according to any embodiment of this invention.

[0100] This invention also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform a fine registration method for sentinel satellite data according to any embodiment of this invention.

[0101] Specifically, a system or apparatus equipped with a storage medium may be provided, on which software program code implementing the functions of any of the embodiments described above is stored, and the computer (or CPU or MPU) of the system or apparatus may read and execute the program code stored in the storage medium.

[0102] In this case, the program code read from the storage medium can itself implement the function of any of the above embodiments, and therefore the program code and the storage medium storing the program code constitute part of the present invention.

[0103] Examples of storage media used to provide program code include floppy disks, hard disks, magneto-optical disks, optical disks (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), magnetic tapes, non-volatile memory cards, and ROMs. Alternatively, program code can be downloaded from a server computer via a communication network.

[0104] Furthermore, it should be clear that not only can the program code read by the computer be executed, but also the operating system or other components operating on the computer can be instructed based on the program code to perform some or all of the actual operations, thereby realizing the function of any of the embodiments described above.

[0105] Furthermore, it is understood that the program code read from the storage medium is written to the memory set in the expansion board inserted into the computer or to the memory set in the expansion module connected to the computer. Then, based on the instructions of the program code, the CPU or other components installed on the expansion board or expansion module execute some and all of the actual operations, thereby realizing the function of any of the above embodiments.

[0106] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus.

[0107] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media that can store program code, such as ROM, RAM, magnetic disk, or optical disk.

[0108] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for precise registration of Sentinel satellite data, characterized in that, include: Obtain the upper and lower frequency band data of the main image and the auxiliary image after coarse registration; Based on the upper and lower frequency band data of the main image and the auxiliary image, the coherence and differential phase of each upper and lower frequency band are calculated traversally. The coherence and differential phase of the upper and lower frequency bands are subjected to multi-visualization processing. The differential phases corresponding to coherence less than a preset threshold after multi-visualization processing are set to zero. The differential phases after setting to zero and the differential phases corresponding to coherence greater than the preset threshold are used as a differential phase set. The preset threshold is 0.

5. The mode of the differential phase set is selected as the differential phase of the enhanced spectral diversity; Based on the differential phase, pulse repetition frequency, and changes in the center frequencies of adjacent burst Doppler signals of the enhanced spectral diversity, the azimuth residual offset of the enhanced spectral diversity is calculated to complete the fine registration of the Sentinel satellite data.

2. The method according to claim 1, characterized in that, The upper and lower frequency band data of the main image includes the data of each adjacent burst overlapping region in the main image, and the upper and lower frequency band data of the auxiliary image includes the data of each adjacent burst overlapping region in the auxiliary image.

3. The method according to claim 2, characterized in that, The step of calculating the coherence and differential phase of the upper and lower frequency bands based on the upper and lower frequency band data of the main image and the auxiliary image includes: Based on the upper frequency band data of the main image and the upper frequency band data of the auxiliary image, the interference data of the upper and lower frequency bands are calculated respectively. Based on the interference data of the upper and lower frequency bands, the coherence and differential phase of the upper and lower frequency bands are calculated respectively.

4. The method according to claim 3, characterized in that, The interference data of the upper and lower frequency bands are calculated using the following formula: ， In the formula, The interference data for the upper frequency band, The interferometric data for the lower frequency band, This refers to the upper frequency band data in the main image. is the conjugate complex number of the upper frequency band data in the auxiliary image. This refers to the lower frequency band data in the main image. The conjugate complex number of the lower frequency band data in the auxiliary image; and / or The coherence of the upper and lower frequency bands and the differential phase of the upper and lower frequency bands are calculated using the following formulas: In the formula, The coherence of the upper and lower frequency bands is given by M, where M is the size of the coherence window in the azimuth direction and N is the size of the coherence window in the range image. The interference data for the upper frequency band, The interferometric data for the lower frequency band, Let be the conjugate complex number of the interferometric data in the lower frequency band, i be the coordinate number of the pixel in the azimuth coherence window, and j be the coordinate number of the pixel in the range coherence window. The differential phase of the upper and lower frequency bands, This is a phase-taking operation.

5. The method according to claim 1, characterized in that, The coherence after multi-visualization processing and the differential phase after multi-visualization processing are calculated using the following formulas: In the formula, The coherence of the upper and lower frequency bands after the multi-viewing processing is given by A, where A is the multi-view factor in the azimuth direction and R is the multi-view factor in the range direction. The coherence of the upper and lower frequency bands is given by , where i is the coordinate number of the pixel in the azimuth multi-view window and j is the coordinate number of the pixel in the range multi-view window. The differential phase of the upper and lower frequency bands after multi-visualization processing The differential phase of the upper and lower frequency bands.

6. The method according to any one of claims 1 to 5, characterized in that, The azimuth residual offset of the enhanced spectral diversity is calculated using the following formula: In the formula, The azimuth residual offset of the enhanced spectral diversity is given. The differential phase of the enhanced spectral diversity, The value representing the frequency variation of the adjacent burst Doppler centers. The pulse repetition frequency is the direction in which the pulse is directed.

7. A fine registration device for Sentinel satellite data, characterized in that, For implementing the method as described in any one of claims 1 to 6, comprising: The acquisition unit is used to acquire the upper and lower frequency band data of the main image and the upper and lower frequency band data of the auxiliary image after coarse registration. The first calculation unit is used to calculate the coherence and differential phase of each of the upper and lower frequency bands based on the upper and lower frequency band data of the main image and the auxiliary image. The determining unit is used to perform multi-visualization processing on the coherence and differential phase of the upper and lower frequency bands, and determine the differential phase of the enhanced spectrum diversity based on the coherence and differential phase of the upper and lower frequency bands after multi-visualization processing. The second calculation unit is used to calculate the azimuth residual offset of the enhanced spectrum diversity based on the differential phase, pulse repetition frequency, and changes in the center frequencies of adjacent burst Doppler signals, so as to complete the fine registration of the Sentinel satellite data.

8. A computing device comprising a memory and a processor, wherein the memory stores a computer program, and the processor, when executing the computer program, implements the method as described in any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to perform the method of any one of claims 1-6.

Images