CN112017217B - Registration method and device between remote sensing data channels of multi-frequency terahertz detector - Google Patents
Registration method and device between remote sensing data channels of multi-frequency terahertz detector Download PDFInfo
- Publication number
- CN112017217B CN112017217B CN202010907597.3A CN202010907597A CN112017217B CN 112017217 B CN112017217 B CN 112017217B CN 202010907597 A CN202010907597 A CN 202010907597A CN 112017217 B CN112017217 B CN 112017217B
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- sensing data
- radio
- source
- data channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20132—Image cropping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30192—Weather; Meteorology
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The disclosure discloses a registration method and a registration device among remote sensing data channels of a multi-frequency terahertz detector. The method comprises the following steps: acquiring a basic radio frequency source observation image of each remote sensing data channel, wherein a radio frequency source disc in the basic radio frequency source observation image is an excessive bulge or a deficient bulge; respectively carrying out preliminary registration on the observation images of the basic radio source, then respectively cutting out the images of the radio source areas and amplifying the images by preset times; detecting the edge position of the radio source in the radio source region image amplified by each remote sensing data channel; and acquiring registration parameters of each remote sensing data channel according to the edge position of the remote sensing data channel radio source. The registration parameters of the remote sensing data channels acquired by the registration method among the remote sensing data channels are used for registering the remote sensing data channels, so that the radio sources in the radio source observation images of the remote sensing data channels are the same, and the positions of the radio sources are the same.
Description
Technical Field
The disclosure relates to a registration method between remote sensing data channels of a multi-frequency terahertz detector, and also relates to a corresponding registration device, belonging to the technical field of satellite remote sensing.
Background
The numerical weather forecast depends on geophysical parameters, such as temperature and humidity profiles, wind profiles and the like, which are mostly obtained by means of joint inversion of observation data of a plurality of channels (also called frequency bands). When joint inversion is carried out on the observation data of a plurality of channels, high-precision registration processing of the observation data is the basis. If the spatial registration (the geographic positioning of the same target by different channels) among a plurality of channels has large errors, the measurement result cannot be used for the joint inversion of the same position parameter. When the atmospheric parameters are inverted by using the observation data of a plurality of channels, the absolute geographic positioning of the remote sensing data is required to be very accurate, and the registration among the remote sensing data of different channels must also reach very high accuracy. Therefore, determining the registration parameters among the channels of the microwave remote sensing instrument is important.
The factors causing registration errors among channels of the microwave remote sensing instrument are many, the registration among the channels of the existing microwave remote sensing instrument is mostly carried out on the basis of laboratory measurement data of instrument design in a satellite development stage, and the influence caused by the change of thermal and mechanical environments of the microwave remote sensing instrument during satellite emission or in-orbit is not considered. Therefore, it is very necessary to design a method suitable for determining the registration parameters between the remote sensing data channels of the satellite-borne multi-frequency terahertz detector.
Although there are a lot of documents discussing image registration techniques in the fields of medical images and remote sensing, there are few articles related to registration errors between remote sensing data channels of a multi-frequency terahertz detector. In the existing inter-channel registration technology, landmark optical images are used as research objects, each channel image is resampled to the same spatial resolution, and an inter-channel registration parameter is calculated from the image matching angle by using a certain channel image as a reference channel. The calculation results of the inter-channel registration parameters depend on the resampling algorithm to a large extent, so that there is no persuasion to systematically evaluate the inter-channel registration errors. In addition, when the inter-channel registration of the GOES satellite imager is researched, it is found that the reliability of the inter-channel registration result obtained by utilizing landmark calculation cannot be ensured under the condition of cloud occlusion.
Disclosure of Invention
The primary technical problem to be solved by the present disclosure is to provide a registration method between remote sensing data channels of a multi-frequency terahertz detector.
Another technical problem to be solved by the present disclosure is to provide a registration device between remote sensing data channels of a multi-frequency terahertz detector.
In order to achieve the purpose, the following technical scheme is adopted in the disclosure:
according to a first aspect of the embodiments of the present disclosure, a registration method between remote sensing data channels of a multi-frequency terahertz detector is provided, which includes the following steps:
acquiring a basic radio frequency source observation image of each remote sensing data channel, wherein a radio frequency source disc in the basic radio frequency source observation image is an excessive bulge or a deficient bulge;
respectively cutting out the images of the radio source areas and amplifying the images by preset times after the basic radio source observation images are respectively subjected to preliminary registration;
detecting the edge position of the radio source in the radio source region image amplified by each remote sensing data channel;
and acquiring registration parameters of each remote sensing data channel according to the edge position of the remote sensing data channel radio frequency source.
Preferably, the process of obtaining the basic radio source observation image of each remote sensing data channel is as follows:
determining the radio frequency power source observation data of the excess or deficiency convexity in the radio frequency power source original observation data of each remote sensing data channel, and storing the task number of the radio frequency power source observation data in a record;
analyzing the radio source observation data corresponding to the task number screened out from the record to obtain a multi-row radio source observation image of each remote sensing data channel;
and correcting the sawtooth boundaries of the radio sources in the multi-row radio source observation images of the remote sensing data channels to obtain basic radio source observation images of the remote sensing data channels.
Preferably, the process of correcting the sawtooth boundary of the radio source in the multi-row radio source observation image of each remote sensing data channel comprises:
determining the edge position of the radio source in each row of radio source observation images of each remote sensing data channel by using an edge detection algorithm;
respectively calculating the radio source edge pixel deviation of the rest rows of radio source observation images by taking the radio source edge pixel of the first row of radio source observation images of each remote sensing data channel as a reference;
and aligning the radio source edge pixels of the rest row radio source observation images of each remote sensing data channel with the radio source edge pixels of the first row radio source observation images respectively to obtain the basic radio source observation images of each remote sensing data channel.
Preferably, the preliminary registration of the basic radio frequency source observation image is completed according to registration parameters between each remote sensing data channel measured in a laboratory, and the registration parameters between each remote sensing data channel measured in the laboratory are offset of the basic radio frequency source observation image of each remote sensing data channel relative to a central view field.
Preferably, the process of cutting the image of the radio source region from the preliminarily registered radio source image of each remote sensing data channel comprises the following steps:
on the basis of the radio source image after the preliminary registration of the remote sensing data channel with low resolution, cutting a radio source area containing a small number of cold air, storing the radio source area image into a radio source area image, and recording a scanning line number and a packet sequence number corresponding to the cut radio source area;
and selecting the radio frequency source area of the rest remote sensing data channel according to the scanning line number and the packet sequence number corresponding to the cut radio frequency source area image, and storing the radio frequency source area in the radio frequency source area image.
Preferably, the edge position of the radio source in the radio source region image amplified by each remote sensing data channel comprises coordinates of a plurality of fitting points enclosing a synthetic radio source edge profile.
Preferably, the disc center coordinate of the radio source in the basic radio source observation image of a certain remote sensing data channel closest to the central view field is taken as a reference, and the offset between the disc center coordinate of the radio source in the basic radio source observation image of the rest remote sensing data channels and the disc center coordinate of the radio source of the reference remote sensing data channel is respectively calculated to obtain the registration parameter of each remote sensing data channel.
Preferably, the center of the central view field is used as a reference, and the offset between the disk center coordinate of the radio source in the basic radio source observation image of each remote sensing data channel and the center of the central view field is respectively calculated to obtain the registration parameter of each remote sensing data channel.
Preferably, the formula is adopted according to the coordinates (x, y) of a plurality of fitting points enclosing and synthesizing the edge contour of the radio source of each remote sensing data channel2+y2And + Dx + Ey + F is 0, and the coordinates of the center of the radio source disc in the radio source area image amplified by each remote sensing data channel are calculated as follows:
wherein D, E, F is a fitting parameter.
According to a second aspect of the embodiments of the present disclosure, there is provided an apparatus for registration between remote sensing data channels of a multi-frequency terahertz detector, including a processor and a memory, where the processor reads a computer program or instructions in the memory to perform the following operations:
acquiring a basic radio frequency source observation image of each remote sensing data channel, wherein a radio frequency source disc in the basic radio frequency source observation image is an excessive bulge or a deficient bulge;
respectively cutting out the images of the radio source areas and amplifying the images by preset times after the basic radio source observation images are respectively subjected to preliminary registration;
detecting the edge position of the radio source in the radio source region image amplified by each remote sensing data channel;
and acquiring registration parameters of each remote sensing data channel according to the edge position of the remote sensing data channel radio frequency source.
And acquiring registration parameters of each remote sensing data channel according to the edge position information of the radio frequency source of each remote sensing data channel.
The registration method and the registration device for the remote sensing data channels of the multi-frequency terahertz detector are characterized in that after primary registration is carried out on a basic radio-frequency source observation image of each remote sensing data channel based on whether a radio-frequency source disc is convex or deficient, a radio-frequency source region image is cut out and amplified by preset times, and then registration parameters of each remote sensing data channel are obtained according to edge position information of a radio-frequency source in the radio-frequency source region image amplified by each remote sensing data channel obtained through detection, so that registration is carried out on each remote sensing data channel according to the registration parameters of each remote sensing data channel, and the radio-frequency sources in the radio-frequency source observation images of each remote sensing data channel are the same and have the same position.
Drawings
Fig. 1 is a flowchart of a registration method between remote sensing data channels of a multi-frequency terahertz detector provided in an embodiment of the present disclosure;
fig. 2 is a schematic diagram of a radio source center fitting in the registration method between remote sensing data channels of the multi-frequency terahertz detector provided by the embodiment of the present disclosure;
fig. 3 is a schematic structural diagram of a registration device between remote sensing data channels of a multi-frequency terahertz detector provided by an embodiment of the present invention.
Detailed Description
The technical contents of the present disclosure are further described in detail below with reference to the accompanying drawings and specific embodiments.
The moon, the planet, the radio astrology and the like can be used as a radio source for the multi-frequency terahertz detector, and the cold air with the universe background is arranged around the multi-frequency terahertz detector, so that the influence of natural environment factors such as cloud and atmosphere does not exist in the radio source observation image, and the registration parameters among remote sensing data channels of the multi-frequency terahertz detector can be determined based on the radio source observation image without being influenced by the factors. Therefore, as shown in fig. 1, an embodiment of the present disclosure provides a method for registering between remote sensing data channels of a multi-frequency terahertz detector, where the method determines registration parameters between remote sensing data channels of the multi-frequency terahertz detector based on a radio source observation image, and includes the following steps:
and step S1, obtaining a basic radio source observation image of each remote sensing data channel.
Because the field of view of the stationary orbit satellite-borne multi-frequency terahertz detector is small, the multi-frequency terahertz detector can transversely scan the radio source in east and west directions in each remote sensing data channel through joint adjustment of attitude maneuver of a satellite platform and antenna pointing direction, and can move line by line in north and south directions, so that the full coverage scanning of the radio source is realized, a plurality of radio source observation images corresponding to each remote sensing data channel are obtained, and the complete radio source observation image of each remote sensing data channel can be obtained only by splicing the plurality of radio source observation images of each remote sensing data channel together. Therefore, one radio source observation task is composed of N scanning lines, the sampling data of each scanning line is sequentially composed into p data packets according to the actual sampling number p of the low spatial resolution channel, and correspondingly, one data packet of the high spatial resolution channel comprises a plurality of sampling data. Thus, each data packet may be determined by both the scan line number and the packet sequence number.
In the embodiment of the disclosure, in the obtained basic radio source observation image of each remote sensing data channel, the area of the radio source occupies more than half of the whole size of the radio source on the image surface, that is, the radio source disc in the basic radio source observation image of each remote sensing data channel is an excessive projection or a deficient projection, so that the registration parameters between the remote sensing data channels of the multi-frequency terahertz detector can be calculated. The process of obtaining the basic radio source observation image of each remote sensing data channel is as follows:
and step S11, determining the radio source observation data of the excess or deficiency convexity of the radio source disc in the radio source original observation data of each remote sensing data channel, and storing the task number of the radio source observation data in a record.
The implementation process of the step is as follows:
and S110, acquiring original observation images of the radio sources of all remote sensing data channels.
Information such as the sampling data, observation time, observation mode and the like of the radio frequency source of each remote sensing data channel acquired by the multi-frequency terahertz detector can be packaged into a corresponding remote sensing data packet to be sent to the ground data processing device. And the ground data processing system sequentially unpacks the remote sensing data packets according to the remote sensing data packet format to form the original observation data of the radio source of each remote sensing data channel.
And step S111, screening observation data meeting the requirements based on the original observation image, and storing the task number of the observation data in a record.
The radio source is regarded as a sphere, and the radio source observed by the geostationary satellite is a disc surface. If the field angle of the radio power supply to the satellite is theta and the angular resolution of the satellite-borne remote sensor is Res, the number of pixels lightened on the image surface by the whole disc of the radio power supply is about
N＝π(θ/2Res)2 (1)
And if the number of the lighted pixels of the radio source in the image exceeds N/2 in the original observation data of the radio source of each remote sensing data channel, judging that the radio source of the observation data is excessive or deficient, and storing the task number of the observation data in a record.
And step S12, analyzing the radio source observation data corresponding to the task number screened from the record to obtain a multi-row radio source observation image of each remote sensing data channel.
And the task number of the observation data of the radioactive source with the excessive or deficient convexity is determined according to the step S11. And analyzing the original observation data of the radio source corresponding to the task number to obtain a multi-row radio source observation image of each remote sensing data channel.
And S13, correcting the sawtooth boundary of the radio sources in the multi-row radio source observation image of each remote sensing data channel to obtain a basic radio source observation image of each remote sensing data channel.
In the integration time, the relative motion of the radio source and the satellite causes saw-toothed shapes to appear between adjacent rows of radio source observation images of each remote sensing data channel acquired by the multi-frequency terahertz detector, so that the radio source observation images of each remote sensing data channel correspondingly composed of the multiple rows of radio source observation images of each remote sensing data channel acquired by the multi-frequency terahertz detector are saw-toothed, for example, the moon moves towards the east at a speed of about 15.04 degrees/hour. Due to the fact that saw-tooth shapes exist in the radio frequency source observation images of the remote sensing data channels and seriously affect the calculation of registration parameters among the remote sensing data channels of the multi-frequency terahertz detector, the saw-tooth boundaries of the radio frequency sources in the multi-row radio frequency source observation images of the remote sensing data channels need to be corrected respectively, and basic radio frequency source observation images of the remote sensing data channels are obtained. The process of correcting the sawtooth boundary of the radio source in the multi-row radio source observation image of each remote sensing data channel comprises the following steps:
and S130, determining the position of the radio source edge in each line of the radio source observation image of each remote sensing data channel by using an edge detection algorithm.
Determining the edge position of the radio source in each row of radio source observation image of each remote sensing data channel by using an edge detection algorithm is the prior mature technology and is not described herein again.
And S131, respectively calculating the radio source edge pixel deviation of the rest rows of radio source observation images by taking the radio source edge pixel of the first row of radio source observation images of each remote sensing data channel as a reference.
And respectively calculating the deviation between the radio source edge pixel of the rest row of radio source observation images and the radio source edge pixel of the first row of radio source observation images by taking the radio source edge pixel of the first row of radio source observation images of each remote sensing data channel as a reference. For example, the west edge pixel of the radio source of the first row of the radio source observation image is the 5 th pixel, and the west edge pixels of the radio sources of the second and third row of the radio source observation images are the 3 rd pixel and the 7 th pixel, respectively, then the west edge pixel of the radio source of the second row of the radio source observation image is shifted by 2 pixels to the right compared with the west edge pixel of the radio source of the first row of the radio source observation image, and the west edge pixel of the radio source of the third row of the radio source observation image is shifted by 2 pixels to the left compared with the west edge pixel of the radio source of the first row of the radio source observation image.
And S132, aligning the radio frequency source edge pixels of the rest row radio frequency source observation images of each remote sensing data channel with the radio frequency source edge pixels of the first row radio frequency source observation images respectively to obtain basic radio frequency source observation images of each remote sensing data channel.
And translating the rest of the remote sensing data channel radio source observation images line by line according to the radio source edge pixel deviation of the rest of the remote sensing data channel radio source observation images calculated in the step S131, so that the radio source edge pixels of the rest of the remote sensing data channel radio source observation images are respectively aligned with the radio source edge pixels of the first row of the radio source observation images, and thus obtaining the basic radio source observation images of the remote sensing data channels.
And step S2, respectively cutting out the images of the radio source areas and amplifying the images by preset times after respectively carrying out primary registration on the basic radio source observation images of each remote sensing data channel.
When the preliminary registration is respectively carried out on the basic radio source observation images of each remote sensing data channel, the preliminary registration of the basic radio source observation images needs to be completed according to registration parameters measured in a laboratory among each remote sensing data channel. The registration parameter measured by the laboratory among each remote sensing data channel is the offset of a basic radio source observation image of each remote sensing data channel relative to a central view field. According to the offset of the basic radio source observation image of each remote sensing data channel relative to the central view field, the basic radio source observation image of each remote sensing data channel is translated integrally, so that the basic radio source observation image of each remote sensing data channel is aligned to the central view field, and the preliminarily registered radio source image { I ] of each remote sensing data channel is obtainedN×MAnd the radio source image consists of N rows and M columns of pixels, and N, M is a positive integer.
Since the registration influence of the change of the thermal and mechanical environments of the multi-frequency terahertz detector during satellite transmission or in-orbit on the remote sensing data channels is small, the relative deviation between the radio source images of the channels subjected to preliminary registration is small.
In order to ensure that all remote sensing data channels can obtain a complete radio frequency source image, when a radio frequency source observation task is arranged, enough time and space are reserved by considering registration factors of all the remote sensing data channels, so that a preliminarily registered radio frequency source image { I ] of each remote sensing data channelN×MOften, a large number of cold air areas are contained in the heat exchanger. To reduce the amount of computation of registration parameters for each remote sensing data channel and improve computational efficiency, the preliminary registered radio source image { I ] of each remote sensing data channel can be obtainedN×MCutting out interested radio source area image { J }n×m}，The image of the radio frequency source area consists of n rows of pixels and m columns of pixels, wherein n and m are positive integers. Primary registered RF source images I from remote sensing data channelsN×MThe process of cutting the image of the radio source area in the previous step is as follows:
and step S21, on the basis of the radio source image after the preliminary registration of the remote sensing data channel with low resolution, cutting a radio source area containing a small amount of cold air, storing the radio source area in the radio source area image, and recording the scanning line number and the packet sequence number corresponding to the cut radio source area.
For faster and more accurate preliminary registration of the RF source image I from each telemetry data channelN×MCutting the radio source area image, and obtaining the preliminarily registered radio source image (I) of all remote sensing data channelsN×MSelecting a radio image with low spatial resolution { I }N×MAnd at a selected low spatial resolution radio source image { I }N×MBased on the image, cutting a radio source area containing a small amount of cold air, and storing the radio source area image into a radio source area image { J }n×mAnd recording the scanning line number and the packet sequence number corresponding to the cut image. As the primary registered radio source images { I ] of all remote sensing data channelsN×MThe remote sensing data channel comprises a plurality of channels, wherein the channels are divided into a plurality of low-spatial resolution channels, the low-spatial resolution channels are divided into a plurality of low-spatial resolution channels, and the low-spatial resolution channels are divided into a plurality of low-spatial resolution channels.
Wherein, the size of the cut radio source region containing a small amount of cold air is determined according to experience, and simultaneously, the radio source region images { J ] of all remote sensing data channels are ensuredn×mAll contain a complete radio image and a small amount of cold air.
And S22, selecting the radio frequency source area of the rest remote sensing data channel according to the scanning line row number and the packet sequence number corresponding to the cut radio frequency source area image, and storing the radio frequency source area image.
According to the scanning line number and the packet number corresponding to the radio frequency source region cut in the step S21, determining the radio frequency source region size of the residual remote sensing data channel corresponding to the radio frequency source region cut in the step S21 from the original observation data of the radio frequency source corresponding to the task number, and according to the determined radio frequency source region size of the residual remote sensing data channel, performing primary registration on the residual remote sensing data channel to obtain a radio frequency source image { I }N×MBased on the image, cutting a radio source area containing a small amount of cold air, and storing the radio source area image into a radio source area image { J }n×m}。
Because the multi-frequency terahertz detector has low spatial resolution, the area occupied by the radio source in the observation image is small, and the radio source region image { J } based on each remote sensing data channel is not favorablen×mCalculating registration parameters among remote sensing data channels of the multi-frequency terahertz detector, and therefore, acquiring radio source area images { J } of each remote sensing data channeln×mAnd (6) carrying out interpolation and amplification by preset times. Wherein, the radio source region image { J ] of each remote sensing data channeln×mAnd the preset times of interpolation amplification are specifically determined according to the image of the radio power region.
And step S3, detecting the edge position of the radio source in the radio source area image amplified by each remote sensing data channel.
And detecting the edge of the radio source in the radio source region image amplified by each remote sensing data channel by adopting common edge detection operators such as Canny, Sobel and the like to obtain the edge position information of the radio source in the radio source region image amplified by each remote sensing data channel. The edge position information of the radio source in the radio source region image amplified by each remote sensing data channel comprises coordinates of a plurality of fitting points enclosing and synthesizing a radio source edge outline.
And step S4, acquiring registration parameters of each remote sensing data channel according to the edge position of the radio frequency source of each remote sensing data channel.
As shown in fig. 2, according to the coordinates (xj, yj) j of the plurality of fitting points which enclose the edge contour of each remote sensing data channel radio source, which are 1,2, … l, j is the number of fitting points which enclose the edge contour of each remote sensing data channel radio source, l is a positive integer, the radius R of the radio source disc and the coordinates (xci, yci) of the center of the radio source disc in the radio source area image amplified by each remote sensing data channel are calculated by using the fitting solving formula (2), i is 1,2, …, p, i is the number of remote sensing data channels, and p is a positive integer.
x2+y2+Dx+Ey+F＝0 (2)
In the above formula, (x, y) is the coordinates of a certain fitting point around the edge contour of the radio source of each telemetric data channel, and D, E, F is a fitting parameter.
According to the formula 2, the radius R of the radio source disc and the coordinates (xci, yci) of the center of the radio source disc in the radio source area image amplified by each remote sensing data channel are respectively:
and respectively calculating the offset between the disk center coordinates of the radio sources in the basic radio source observation images of the rest remote sensing data channels and the disk center coordinates of the radio sources in the reference remote sensing data channels by taking the disk center coordinates of the radio sources in the basic radio source observation images of one remote sensing data channel closest to the central view field as a reference, namely the registration parameters of each remote sensing data channel. And correspondingly translating the disc center of the radio source in the basic radio source observation image of the rest remote sensing data channel according to the calculated registration parameter of each remote sensing data channel, so that the disc center of the radio source in the basic radio source observation image of the rest remote sensing data channel is aligned to the disc center of the radio source in the reference remote sensing data channel, thereby ensuring that the positions of the radio sources in the radio source observation images of each remote sensing data channel are the same.
Or respectively calculating the offset between the disk center coordinate of the radio source in the basic radio source observation image of each remote sensing data channel and the center of the central view field by taking the center of the central view field as a reference, namely the registration parameter of each remote sensing data channel. And correspondingly translating the disc center of the radio source in the basic radio source observation image of each remote sensing data channel according to the calculated registration parameter of each remote sensing data channel, so that the disc center of the radio source in the basic radio source observation image of each remote sensing data channel is aligned to the center of the central view field, thereby ensuring that the positions of the radio sources in the radio source observation images of each remote sensing data channel are the same.
It should be noted that the inter-channel registration parameters obtained by calculation are often in a decimal form, and resampling of observed data is required during translation, and an existing resampling algorithm is mature and is not described herein again.
Further, as shown in fig. 3, the present disclosure also provides a device for registering between remote sensing data channels of a multi-frequency terahertz detector, which includes a processor 32 and a memory 31, and may further include a communication component, a sensor component, a power supply component, a multimedia component, and an input/output interface according to actual needs. The memory, communication components, sensor components, power components, multimedia components, and input/output interfaces are all connected to the processor 32. As mentioned above, the memory 31 may be a Static Random Access Memory (SRAM), an Electrically Erasable Programmable Read Only Memory (EEPROM), an Erasable Programmable Read Only Memory (EPROM), a Programmable Read Only Memory (PROM), a Read Only Memory (ROM), a magnetic memory, a flash memory, etc.; the processor may be a Central Processing Unit (CPU), Graphics Processing Unit (GPU), Field Programmable Gate Array (FPGA), Application Specific Integrated Circuit (ASIC), Digital Signal Processing (DSP) chip, etc. Other communication components, sensor components, power components, multimedia components, etc. may be implemented using common components found in existing smartphones and are not specifically described herein.
On the other hand, in the registration device between the remote sensing data channels of the multi-frequency terahertz detector, the processor 32 reads a computer program or instructions in the memory 31 for executing the following operations:
and acquiring a basic radio source observation image of each remote sensing data channel.
And respectively cutting out the images of the radio source areas and amplifying the images by preset times after respectively carrying out primary registration on the basic radio source observation images of each remote sensing data channel.
And detecting the edge position of the radio source in the radio source region image amplified by each remote sensing data channel.
And acquiring registration parameters of each remote sensing data channel according to the edge position information of the radio frequency source of each remote sensing data channel.
The registration method and the registration device for the remote sensing data channels of the multi-frequency terahertz detector provided by the embodiment of the disclosure are based on the fact that the radio frequency power disc is used for carrying out primary registration on basic radio frequency power source observation images of all the remote sensing data channels which are full of protrusions or lack of protrusions, then the radio frequency power source region images are respectively cut out and amplified by preset times, and then registration parameters of all the remote sensing data channels are obtained according to the edge position information of the radio frequency power source in the radio frequency power source region images obtained by detection and amplification of all the remote sensing data channels, so that the registration of all the remote sensing data channels is conveniently carried out according to the registration parameters of all the remote sensing data channels, and the radio frequency power sources in the radio frequency power source observation images of all the remote sensing data channels are the same and have the same position. And moreover, based on the basic radio-frequency source observation image of each remote sensing data channel with the convex or deficient radio-frequency source disc, the calculated registration parameters among the remote sensing data channels are more accurate, and the on-orbit state of the multi-frequency terahertz detector can be more accurately reflected.
The registration method and device between the remote sensing data channels of the multi-frequency terahertz detector provided by the present disclosure are explained in detail above. It will be apparent to those skilled in the art that various modifications can be made without departing from the spirit of the disclosure, and the scope of the disclosure is to be accorded the full scope of the claims appended hereto.
Claims (10)
1. A registration method between remote sensing data channels of a multi-frequency terahertz detector is characterized by comprising the following steps:
acquiring a basic radio frequency source observation image of each remote sensing data channel, wherein a radio frequency source disc in the basic radio frequency source observation image is an excessive bulge or a deficient bulge;
respectively cutting out the images of the radio source areas and amplifying the images by preset times after the basic radio source observation images are respectively subjected to preliminary registration;
detecting the edge position of the radio source in the radio source region image amplified by each remote sensing data channel;
and acquiring registration parameters of each remote sensing data channel according to the edge position of the remote sensing data channel radio frequency source.
2. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 1, wherein:
the process of obtaining the basic radio source observation image of each remote sensing data channel comprises the following steps:
determining the radio frequency power source observation data of the excess or deficiency convexity in the radio frequency power source original observation data of each remote sensing data channel, and storing the task number of the radio frequency power source observation data in a record;
analyzing the radio source observation data corresponding to the task number screened out from the record to obtain a multi-row radio source observation image of each remote sensing data channel;
and correcting the sawtooth boundaries of the radio sources in the multi-row radio source observation images of the remote sensing data channels to obtain basic radio source observation images of the remote sensing data channels.
3. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 2, wherein:
the process of correcting the sawtooth boundary of the radio source in the multi-row radio source observation image of each remote sensing data channel comprises the following steps:
determining the edge position of the radio source in each row of radio source observation images of each remote sensing data channel by using an edge detection algorithm;
respectively calculating the radio source edge pixel deviation of the rest rows of radio source observation images by taking the radio source edge pixel of the first row of radio source observation images of each remote sensing data channel as a reference;
and aligning the radio source edge pixels of the rest row radio source observation images of each remote sensing data channel with the radio source edge pixels of the first row radio source observation images respectively to obtain the basic radio source observation images of each remote sensing data channel.
4. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 1, wherein:
and completing preliminary registration of the basic radio frequency source observation image according to registration parameters between each remote sensing data channel measured in a laboratory, wherein the registration parameters between each remote sensing data channel measured in the laboratory are the offset of the basic radio frequency source observation image of each remote sensing data channel relative to a central view field.
5. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 1, wherein:
the process of cutting the radio source region image from the radio source image after the preliminary registration of each remote sensing data channel comprises the following steps:
on the basis of the radio source image after the preliminary registration of the remote sensing data channel with low resolution, cutting a radio source area containing a small number of cold air, storing the radio source area image into a radio source area image, and recording a scanning line number and a packet sequence number corresponding to the cut radio source area;
and selecting the radio frequency source area of the rest remote sensing data channel according to the scanning line number and the packet sequence number corresponding to the cut radio frequency source area image, and storing the radio frequency source area in the radio frequency source area image.
6. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 1, wherein:
the edge position of the radio source in the radio source region image amplified by each remote sensing data channel comprises the coordinates of a plurality of fitting points which surround the edge outline of the synthetic radio source.
7. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 6, wherein:
and respectively calculating the offset between the disc center coordinate of the radio source in the basic radio source observation image of the rest remote sensing data channels and the disc center coordinate of the radio source in the reference remote sensing data channel by taking the disc center coordinate of the radio source in the basic radio source observation image of one remote sensing data channel closest to the central view field as a reference, and obtaining the registration parameters of each remote sensing data channel.
8. The method for registering the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 6, wherein:
and respectively calculating the offset between the disk center coordinates of the radio sources in the basic radio source observation images of all the remote sensing data channels and the center of the central view field by taking the center of the central view field as a reference to obtain the registration parameters of all the remote sensing data channels.
9. The method for registering between the remote sensing data channels of the multi-frequency terahertz detector as claimed in claim 7 or 8, wherein:
according to the coordinates (x, y) of a plurality of fitting points which surround and synthesize the edge contour of the radio source of each remote sensing data channel, adopting a formula, x2+y2And + Dx + Ey + F is 0, and the coordinates of the center of the radio source disc in the radio source area image amplified by each remote sensing data channel are calculated as follows:
wherein D, E, F is a fitting parameter.
10. A registration device between remote sensing data channels of a multi-frequency terahertz detector comprises a processor and a memory, wherein the processor reads a computer program or instructions in the memory and is used for executing the following operations:
acquiring a basic radio frequency source observation image of each remote sensing data channel, wherein a radio frequency source disc in the basic radio frequency source observation image is an excessive bulge or a deficient bulge;
respectively cutting out the images of the radio source areas and amplifying the images by preset times after the basic radio source observation images are respectively subjected to preliminary registration;
detecting the edge position of the radio source in the radio source region image amplified by each remote sensing data channel;
and acquiring registration parameters of each remote sensing data channel according to the edge position of the remote sensing data channel radio frequency source.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010907597.3A CN112017217B (en) | 2020-09-02 | 2020-09-02 | Registration method and device between remote sensing data channels of multi-frequency terahertz detector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010907597.3A CN112017217B (en) | 2020-09-02 | 2020-09-02 | Registration method and device between remote sensing data channels of multi-frequency terahertz detector |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112017217A CN112017217A (en) | 2020-12-01 |
| CN112017217B true CN112017217B (en) | 2022-05-03 |
Family
ID=73516683
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010907597.3A Active CN112017217B (en) | 2020-09-02 | 2020-09-02 | Registration method and device between remote sensing data channels of multi-frequency terahertz detector |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112017217B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107301661A (en) * | 2017-07-10 | 2017-10-27 | 中国科学院遥感与数字地球研究所 | High-resolution remote sensing image method for registering based on edge point feature |
| CN107392951A (en) * | 2017-06-06 | 2017-11-24 | 上海卫星工程研究所 | Remote sensing images high accuracy rapid registering method |
| CN111145231A (en) * | 2019-12-30 | 2020-05-12 | 国家卫星海洋应用中心 | Method and device for determining waveband offset of remote sensing image and electronic equipment |
| CN111161186A (en) * | 2019-12-30 | 2020-05-15 | 国家卫星气象中心(国家空间天气监测预警中心) | Push-broom remote sensor channel registration method and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11024044B2 (en) * | 2016-09-26 | 2021-06-01 | Digitalglobe, Inc. | Techniques for image co-registration |
-
2020
- 2020-09-02 CN CN202010907597.3A patent/CN112017217B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107392951A (en) * | 2017-06-06 | 2017-11-24 | 上海卫星工程研究所 | Remote sensing images high accuracy rapid registering method |
| CN107301661A (en) * | 2017-07-10 | 2017-10-27 | 中国科学院遥感与数字地球研究所 | High-resolution remote sensing image method for registering based on edge point feature |
| CN111145231A (en) * | 2019-12-30 | 2020-05-12 | 国家卫星海洋应用中心 | Method and device for determining waveband offset of remote sensing image and electronic equipment |
| CN111161186A (en) * | 2019-12-30 | 2020-05-15 | 国家卫星气象中心(国家空间天气监测预警中心) | Push-broom remote sensor channel registration method and device |
Non-Patent Citations (1)
| Title |
|---|
| TDICCD相机大视场多通道遥感图像自动拼接方法;禄金波 等;《空间科学学报》;20121231;第32卷(第1期);154-160 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112017217A (en) | 2020-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schläpfer et al. | Geo-atmospheric processing of airborne imaging spectrometry data. Part 1: Parametric orthorectification | |
| Tang et al. | Triple linear-array image geometry model of ZiYuan-3 surveying satellite and its validation | |
| Gugan | Practical aspects of topographic mapping from SPOT imagery | |
| CN102968788B (en) | Wave band registering method based on regular grid surface element | |
| US20160343152A1 (en) | Methods and systems for object based geometric fitting | |
| Tang et al. | Verification of ZY-3 satellite imagery geometric accuracy without ground control points | |
| Poli | A rigorous model for spaceborne linear array sensors | |
| CN106887016B (en) | Automatic relative registration method for GF-4 satellite sequence images | |
| US20120063668A1 (en) | Spatial accuracy assessment of digital mapping imagery | |
| CN114972545B (en) | On-orbit data rapid preprocessing method for hyperspectral satellite | |
| JPH1183478A (en) | Global information supply system | |
| Poli et al. | Orientation of satellite and airborne imagery from multi-line pushbroom sensors with a rigorous sensor model | |
| CN113469896A (en) | Method for improving geometric correction precision of geosynchronous orbit satellite earth observation image | |
| CN116754076B (en) | Inversion method for high-heterogeneity surface temperature of urban complex three-dimensional scene | |
| CN113313769B (en) | Seamless geometric calibration method between optical satellite multi-area array sensor chips | |
| Chmiel et al. | Orthorectification and geometric quality assessment of very high spatial resolution satellite imagery for Common Agricultural Policy purposes | |
| CN112017217B (en) | Registration method and device between remote sensing data channels of multi-frequency terahertz detector | |
| CN111121729B (en) | Optical remote sensing satellite internal checking method and system based on flat terrain | |
| Erdogan et al. | The georeferencing of RASAT satellite imagery and some practical approaches to increase the georeferencing accuracy | |
| JPS62500055A (en) | Method and system for mapping actual image data from a geostationary satellite to an idealized image | |
| CN112213750B (en) | Hyperspectral guard sheet full-spectrum pixel-by-pixel imaging angle parameter processing method and medium | |
| EP1144946A1 (en) | Improved method of processing data and data processing apparatus | |
| Boukerch et al. | Geometry based co-registration of ALSAT-2A panchromatic and multispectral images | |
| Long et al. | Block adjustment with relaxed constraints from reference images of coarse resolution | |
| Richter et al. | Geo-atmospheric processing of wide-FOV airborne imaging spectrometry data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |