CN104330785A - Marine No. 2 microwave scatterometer surface element matching method - Google Patents
Marine No. 2 microwave scatterometer surface element matching method Download PDFInfo
- Publication number
- CN104330785A CN104330785A CN201410692674.2A CN201410692674A CN104330785A CN 104330785 A CN104330785 A CN 104330785A CN 201410692674 A CN201410692674 A CN 201410692674A CN 104330785 A CN104330785 A CN 104330785A
- Authority
- CN
- China
- Prior art keywords
- spherical distance
- wind vector
- substar
- observation station
- coordinate
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000008569 process Effects 0.000 claims description 31
- 238000005070 sampling Methods 0.000 claims description 4
- 230000002093 peripheral effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000010287 polarization Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000000049 pigment Substances 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/422—Simultaneous measurement of distance and other co-ordinates sequential lobing, e.g. conical scan
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
- G01S13/955—Radar or analogous systems specially adapted for specific applications for meteorological use mounted on satellite
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a marine No.2 microwave scatterometer surface element matching method and relates to the field of satellite observation data applications. The marine No.2 microwave scatterometer surface element matching method comprises the following steps of obtaining coordinates of a sub-satellite point and coordinates of an observation point; calculating the shortest spherical distance from the observation point to an appointed position in a sub-satellite point track and determining coordinates of an intersection point of a straight line which is corresponding to the shortest spherical distance and the sub-satellite point track; calculating the peripheral distance between the intersection coordinates to coordinates of a starting point of a wind vector unit; obtaining the distance along the track intersection direction and along the track following direction between the observation point and the starting point according to the width of the wind vector unit; further calculating the number of wind vector units being far from the starting point of the wind vector unit, namely determining the line number and the column number on which the observation point are arranged, and accordingly achieving a task for expressing a position of the observation point in a wind vector surface element position through the line number and the column number.
Description
Technical field
The present invention relates to Satellite Observations application, No. two microwave scatterometer bin matching process in particular to ocean.
Background technology
Using satellite to observe earth surface is the technology developed for many years, and earth observation satellite comprises earth resources satellite, military reconnaissance satellite, seasat and geodetic satellite (GEOS) etc.Seasat is wherein exactly the detection being mainly used in ocean color pigment, is the industry service such as halobiontic resources open utilization, marine pollution monitoring and control, coastal zone resources are developed, scientific research of seas, a kind of artificial earth satellite that design is launched.
Below HY-2 satellite scatterometer is simply introduced, in August, 2012 launch No. two, ocean, be equipped with first, China can businessization run microwave scatterometer HY2-SCAT.HY2-SCAT is mainly used in Global ocean wind field observation, and the survey wind wind speed range of scatterometer is 4 ~ 24m/s, and wind speed precision is 2m/s or 10%; Wind direction measurement range is 0 ~ 360 °, and wind direction precision is ± 20 °.HY2-SCAT frequency of operation is 13.256GHz, adopts pencil beam conical scanning mode, is rotated, form certain covered ground swath in the motion in satellite platform straight rail direction by pencil beam with fixed elevation around nadir direction.As shown in Figure 1, the satellite height that distance ground is same usually (if satellite is apart from ground H=963KM) is in operation, the flight track at satellite place projects to ground and also just defines sub-satellite track (in fact, satellite is when each scanning, all can produce these data of substar coordinate, substar coordinate i.e. satellite, when scanning, project to ground latitude and longitude coordinates, substar coordinate during Multiple-Scan are connected into sub-satellite track); Simultaneously, satellite when flight can revolving beyond the form of wave beam and interior wave beam carry out scanning that (outer wave beam and interior wave beam all scan in the mode of circumference, the radius of a circle formed as beam scanning interior in figure is 700KM, the radius of a circle that outer beam scanning is formed is 874KM), thus define outer wave beam footmark and interior wave beam footmark, wherein outer wave beam footmark and interior wave beam footmark are the sampling at interval respectively, and swath width just refers to that external wave bundle scans the diameter of the circumference formed.Scatterometer system comprises VV and HH two polarization modes, observe with different incidence angles respectively, different polarization mode can be obtained to same resolution element in the motion process of platform, the repetitive measurement result of different incidence angles degree is (as in Fig. 2, to same target respectively 1, 2, 3 and 4 these four different positions are observed, the observation carried out in position 1 is called as outer wave beam forward sight, the observation carried out in position 2 is called as interior wave beam forward sight, interior wave beam backsight is called as in the observation of position 3, outer wave beam backsight is called as) in the observation of position 4, to overcome the many-valued fuzziness problem of Ocean Wind-field direction inverting.Wherein interior wave beam adopts HH polarization mode, and incident angle is 41 °, and corresponding ground footmark size is about 23km × 31km, and swath width is 1400km.Outer wave beam adopts VV polarization mode, and incident angle is 48 °, and corresponding ground footmark size is about 25km × 38km, and swath width is 1700km.As Fig. 3 shows the schematic diagram of satellite in the enterprising line scanning of space orbit, the crucial point schematic diagram of sub-satellite track and outer wave beam footmark and interior wave beam footmark in Fig. 4.
Simply introduce HY-2 scatterometer data product below, No. two, the ocean current available data product of satellite scatterometer is divided into L1B level product data product, L2A DBMS product, L2B DBMS product and L3 DBMS product.
Wherein L1B data are with the scatterometer observation data of the time of telemetry frame for sequentially carrying out storing.Each telemetry frame comprises 96 scatterometer ranging pulses, each ranging pulse comprises backscattering coefficient, the geographic position of each pulse footprint (outer wave beam footprint and interior wave beam footprint) and being used for describes the parameter of the information such as the quality of measurement data and uncertainty, also comprise the latitude and longitude information of the sub-satellite track (substar is the point that satellite latitude and longitude coordinates place when not surveying transponder pulse, the track that sub-satellite track is made up of a large amount of substar) obtained by positioning system in this data file simultaneously.
L2A product documentation comprises each radar raster-displaying sigma0 (backscattering coefficient) measured value that satellite platform obtains in a space orbit.In addition, L2A product also comprises some assistance data element corresponding with each sigma0 measured value.These assistance data element list position corresponding to each sigma0 measured value, the relevant information such as quality and uncertainty.
In order to be described data more accurately and use, when L1B DBMS is converted into L2A DBMS, need the observation station (point at outer wave beam footprint place and the point at interior wave beam footprint place) in L1B DBMS to be matched respectively in corresponding wind vector unit.As shown in Figure 5, wind vector unit is for benchmark with sub-satellite track (straight rail), with respectively with parallel (straight rail) with vertical (cross rail) in the straight line of sub-satellite track, earth surface is divided into multiple cell, usual each cell is all square, and size is identical.
Sigma0 in L2A product divides into groups with wind vector unit.Each wind vector cell row corresponding ground measures a cross rail cutting of swath.Usually, each L2A wind vector unit is the square of a 25km.Therefore, needs 1624 wind vector cell row complete the once complete covering to the earth.
When L1B data transformations is L2A data, a very important step observation station is matched in the corresponding appropriate unit of wind, to complete the problem using wind vector unit to be described observation station position and to sort out.In correlation technique, traditional bin matching algorithm (the jet power laboratory JPL of the such as U.S.) adopts the mode of inclination Mercator projection (SOM), longitude and latitude first will be observed from terrestrial coordinate trajectory coordinates conversion earthward, carry out bin coupling again, to determine to specify observation station on which wind vector unit (usually using the line number in wind vector unit face and row number to represent the position of the wind vector unit of specifying), but find when using, adopt the method to calculate comparatively complicated.
Summary of the invention
In view of this, the object of the embodiment of the present invention is to provide No. two, ocean microwave scatterometer bin matching process, to improve the wind vector units match speed of satellite scatterometer data.
First aspect, embodiments provides No. two, ocean microwave scatterometer bin matching process, comprising:
According to the substar coordinate obtained in advance and observation station coordinate, calculating observation point arrives the shortest spherical distance of assigned address on sub-satellite track and determines the intersecting point coordinate of straight line corresponding to the shortest spherical distance and sub-satellite track, and sub-satellite track is connected to form by multiple substar coordinate;
According to starting point coordinate and the default wind vector cell width of intersecting point coordinate, wind vector unit, calculating observation point is along the wind vector line number in straight rail direction;
According to the shortest spherical distance and the width with wind vector unit, calculating observation point is along the wind vector row number in cross rail direction;
Determine wind vector line number and wind vector row number corresponding to wind vector cell position be observation station wind vector unit face on position.
In conjunction with first aspect, embodiments provide the first possible embodiment of first aspect, wherein,
According to the substar coordinate obtained in advance and observation station coordinate, the shortest spherical distance that calculating observation point arrives assigned address on sub-satellite track comprises:
Obtain multiple substar coordinate, and observation station coordinate;
The spherical distance of calculating observation point and each substar respectively;
Select spherical distance minimum in multiple spherical distance as the shortest spherical distance.
In conjunction with first aspect, embodiments provide the embodiment that the second of first aspect is possible, wherein, determine that the intersecting point coordinate of straight line corresponding to the shortest spherical distance and sub-satellite track comprises:
Determine that the coordinate of the substar corresponding to spherical distance minimum in the spherical distance of observation station and each substar is intersecting point coordinate.
In conjunction with first aspect, embodiments provide the third possible embodiment of first aspect, wherein, according to starting point coordinate and the default wind vector cell width of intersecting point coordinate, wind vector unit, calculating observation point comprises along the wind vector line number in straight rail direction:
The direct spherical distance of Two coordinate point is calculated according to intersecting point coordinate and starting point coordinate;
Calculate direct spherical distance and round business divided by first of the wind vector cell width preset;
Determine that observation station is first round business+1 along the wind vector line number in straight rail direction.
In conjunction with first aspect, embodiments provide the 4th kind of possible embodiment of first aspect, wherein, according to the shortest spherical distance and the width with wind vector unit, calculating observation point comprises along the wind vector row number in cross rail direction:
Calculate the shortest spherical distance divided by wind vector unit width second round business;
Determine that observation station number is second round business+1 along the wind vector row in cross rail direction.
In conjunction with first aspect, embodiments provide the 5th kind of possible embodiment of first aspect, wherein,
The spherical distance of calculating observation point and each substar comprises respectively:
Multiple continuous substar is divided into many groups;
According to the sampling rule preset, select at least one substar as a reference point in each group;
Calculate the spherical distance of each reference point and observation station respectively, to determine that the group at the reference point place nearest with observation station spherical distance is for reference group;
The reference spherical distance of each substar and observation station in difference computing reference group, to determine with reference to the minimum value in spherical distance for the shortest spherical distance.
In conjunction with first aspect, embodiments provide the 6th kind of possible embodiment of first aspect, wherein, the reference spherical distance of each substar and observation station in difference computing reference group, to determine to comprise for the shortest spherical distance with reference to the minimum value in spherical distance:
Judge whether be less than the first default reference threshold with reference to the minimum value in spherical distance;
If so, then determine that with reference to the minimum value in spherical distance be the shortest spherical distance.
In conjunction with first aspect, embodiments provide the 7th kind of possible embodiment of first aspect, wherein, according to the substar coordinate obtained in advance and observation station coordinate, the shortest spherical distance that calculating observation point arrives assigned address on sub-satellite track comprises:
Calculate the spherical distance of each substar apart from initial substar respectively, initial substar is the starting point of wind vector unit;
The coordinate of the substar that chosen distance observation station is nearest is intersecting point coordinate;
Spherical distance between calculating observation point and the substar nearest apart from observation station is as the shortest spherical distance.
In conjunction with first aspect, embodiments provide the 8th kind of possible embodiment of first aspect, wherein, calculate the spherical distance of each reference point and observation station respectively, to determine that the group at the reference point place nearest with observation station spherical distance also comprises for reference group:
Judge whether spherical distance minimum in the spherical distance of each reference point and observation station is less than the second default reference threshold;
If so, then determine that the group at the reference point place nearest with observation station spherical distance is reference group.
In conjunction with first aspect, embodiments provide the 9th kind of possible embodiment of first aspect, wherein, the first reference threshold is 900-1100KM; Second reference threshold is 900-1100KM.
No. two, the ocean that the embodiment of the present invention provides microwave scatterometer bin matching process, adopt the position determining place wind vector unit of observation station institute face based on the coordinate of substar and observation station, with the mode of employing inclination Mercator projection (SOM) of the prior art, longitude and latitude first will be observed from terrestrial coordinate trajectory coordinates conversion earthward, carry out bin coupling again, make the process of coupling comparatively complexity compare, it is by first getting substar coordinate and observation station coordinate, and calculating observation point arrives the shortest spherical distance of assigned address on sub-satellite track and determines the intersecting point coordinate of straight line corresponding to the shortest spherical distance and sub-satellite track, calculate the spherical distance of the starting point coordinate of intersecting point coordinate distance wind vector unit again, according to the width of wind vector unit, (wind vector unit is square again, wide namely its length of side of wind vector unit), just can draw observation station distance starting point along the Distance geometry in cross rail direction along the distance in straight rail direction, and distance wind vector unit starting point can be calculated further at a distance of how many wind vector unit, also line number and the row number at observation station place are just determined, thus complete and use line number and row number to represent the task of observation station in wind vector unit position, face.
For making above-mentioned purpose of the present invention, feature and advantage become apparent, preferred embodiment cited below particularly, and coordinate appended accompanying drawing, be described in detail below.
Accompanying drawing explanation
In order to be illustrated more clearly in the technical scheme of the embodiment of the present invention, be briefly described to the accompanying drawing used required in embodiment below, be to be understood that, the following drawings illustrate only some embodiment of the present invention, therefore the restriction to scope should be counted as, for those of ordinary skill in the art, under the prerequisite not paying creative work, other relevant accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 shows in correlation technique, and satellite carries out the first schematic diagram of earth scanning;
Fig. 2 shows in correlation technique, and satellite carries out the second schematic diagram of earth scanning;
Fig. 3 shows in correlation technique, and satellite carries out the third schematic diagram of earth scanning;
Fig. 4 shows in correlation technique, and satellite carries out the 4th kind of schematic diagram of earth scanning;
Fig. 5 shows in correlation technique, the straight rail of wind vector unit and cross rail schematic diagram;
Fig. 6 shows the basic flow sheet of No. two microwave scatterometer bin matching process in ocean that the embodiment of the present invention provides;
Fig. 7 shows the preferred obtaining step process flow diagram of the shortest spherical distance of No. two microwave scatterometer bin matching process in ocean that the embodiment of the present invention provides;
Fig. 8 shows the wind vector line number determining step process flow diagram of No. two microwave scatterometer bin matching process in ocean that the embodiment of the present invention provides;
Fig. 9 shows the shortest spherical distance obtaining step details process flow diagram of No. two microwave scatterometer bin matching process in ocean that the embodiment of the present invention provides;
Figure 10 shows the another kind the shortest preferred spherical distance obtaining step process flow diagram of No. two microwave scatterometer bin matching process in ocean that the embodiment of the present invention provides.
Embodiment
Below in conjunction with accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.The assembly of the embodiment of the present invention describing and illustrate in usual accompanying drawing herein can be arranged with various different configuration and design.Therefore, below to the detailed description of the embodiments of the invention provided in the accompanying drawings and the claimed scope of the present invention of not intended to be limiting, but selected embodiment of the present invention is only represented.Based on embodiments of the invention, the every other embodiment that those skilled in the art obtain under the prerequisite not making creative work, all belongs to the scope of protection of the invention.
The embodiment of the present invention 1 provides No. two, ocean microwave scatterometer bin matching process, as shown in Figure 6, comprises the steps:
S101, according to the substar coordinate obtained in advance and observation station coordinate, calculating observation point arrives the shortest spherical distance of assigned address on sub-satellite track and determines the intersecting point coordinate of straight line corresponding to the shortest spherical distance and sub-satellite track, and sub-satellite track is connected to form by multiple substar coordinate;
S102, according to starting point coordinate and the default wind vector cell width of intersecting point coordinate, wind vector unit, calculating observation point is along the wind vector line number in straight rail direction;
S103, according to the shortest spherical distance and the width with wind vector unit, calculating observation point is along the wind vector row number in cross rail direction;
S104, determine wind vector line number and wind vector row number corresponding to wind vector cell position be observation station wind vector unit face on position.
Wind vector unit is on wind vector unit face, for representing the unit of positional information.In L2A data, certain position (as substar position, observation station position, i.e. observation station footmark) particular location on wind vector unit face for convenience of explanation, and refer to this expression way of wind vector unit.
Below scatterometer wind vector retrieval is simply introduced,
Backscattering coefficient is wind speed, wind direction, the electromagnetic wavelength of radar emission, polarization mode, the function of incident angle.The general type of physical geography module function is:
σ
0=F(u,χ,...,f,p,θ) (1)
Satellite scatterometer is a kind of radar through calibration, and it to sea active emitting electromagnetic wave, and receives the echoed signal through sea modulation.Radar echo signal will by transmit and sea feature determines jointly.When wave wavelength and radar emission electromagnetic wavelength meet Bragg diffraction condition, the back scattering electromagnetic wave phase place that each corrugated produces is identical, thus produces resonance, and backward energy determines primarily of the electromagnetic wave producing resonance.Under the frequency of operation of microwave scatterometer, meet the surface wave of Bragg resonance condition for sea table capillary wave, the spectral density of ocean surface capillary wave is directly related with the wind speed on ocean surface.Therefore, the echoed signal recorded by radar can obtain the information of Ocean Wind-field.By the process to radar echo signal, only relevant with sea condition normalization backscattering coefficient (NRCS, or σ can be drawn
0), from the σ that scatterometer records
0can extract Ocean Wind-field further, the information extraction process of Ocean Wind-field is called wind vector retrieval.
The backscattering from ocean surfaces coefficient inverting sea surface wind vector recorded from scatterometer needs the problem of solution three aspects: set up geophysical model, wind vector derivation algorithm, and fuzzy solution removes algorithm.
Physical geography module function describes the relation between sea surface wind vector and Radar backscattering coefficients.The general type of physical geography module function is:
σ
0=F(V,χ,...,f,p,θ) (2)
Wherein σ
0represent the backscattering coefficient that sea is corresponding; V is wind speed; χ is the relative bearing of wind direction; F is the frequency of operation of scatterometer; P is polarization mode; θ is the incident angle of antenna.
Wind vector derivation algorithm mainly passes through the wind vector solution on the NRCS observed reading acquisition sea of the different orientations of physical geography module function and sea surface wind vector bin.The geophysical model represented from (2) formula, solves wind speed and direction (i.e. wind vector) by backscattering coefficient, needs the backscattering coefficient observed result of multiple (being no less than 4) different angles.
Concrete, in L1B data, directly can get each scatterometer that uses and carry out the whole parameters scanned, namely just uniquely can determine one group of mutual corresponding information after having scanned, as the information of the latitude and longitude information of substar, other measurement results such as latitude and longitude information, backscattering coefficient of observation station at every turn.Further, in once complete satellite scanning activity, need to carry out Multiple-Scan, namely can get the information that many groups are mutually corresponding.In step S101, obtaining substar coordinate is the substar coordinate that Multiple-Scan obtains, and observation station coordinate finally needs to determine the coordinate of its coordinate in wind vector unit position, face.The multiple substar coordinates obtained according to Multiple-Scan can determine the coordinate of sub-satellite track, and (as the equation be used under two-dimensional coordinate system expresses this track, or the equation under use three-dimensional system of coordinate expresses this track), and then can determine that observation station arrives the shortest spherical distance of assigned address on sub-satellite track, in order to determine that observation station arrives the shortest spherical distance of sub-satellite track accurately, assigned address can be through observation station, vertical with sub-satellite track, and make observation station arrive the shortest spherical curve (also can be line segment) of sub-satellite track distance, with the intersection point (this intersection point refers to the intersection point in approximate two-dimensional coordinate system) of sub-satellite track.Wherein, through observation station, and the straight line vertical with sub-satellite track can be that straight line in two-dimensional coordinate system is (when swath width is less, straight line sub-satellite track is interpreted as in the plane that can be similar to, nature, through the straight line in the line of this straight line i.e. two-dimensional coordinate system, in two-dimensional coordinate system, the method of calculation level and line bee-line is more, do not repeat them here), also can be in three-dimensional system of coordinate along the curve on earth's surface (when the earth is carried out understanding according to the spheroid of solid time, sub-satellite track is a curve namely, the bee-line that any point on basic mathematical algorithm calculating sphere can be used in three-dimensional system of coordinate equally to arrive a spherical curve does not repeat them here), also the straight line on earth's surface can be through (in time being carried out understanding according to the spheroid of solid by the earth, sub-satellite track is a curve namely, the bee-line of certain some arrival substar on sphere also just can be calculated by dotted line bee-line formula, its computing method are more, do not repeat them here).After the intersecting point coordinate of the straight line determined corresponding to the shortest spherical distance and the shortest spherical distance and sub-satellite track, just step S102 can be performed.
In step S102, needs carry out the wind vector line number of calculating observation point along straight rail direction according to the starting point coordinate of intersecting point coordinate, wind vector unit and default wind vector cell width, it should be noted that, the starting point of wind vector unit is the point that wind vector unit face starts to add up, and its line number and row number are 0.Needing precalculated is starting point coordinate (starting point coordinate of wind vector unit is default) calculating distance therebetween according to intersecting point coordinate and wind vector unit, the account form of sphere distance between two points can simply be calculated after known 2 spherical co-ordinates, earth radius, at this, concrete computation process repeats no more.After the distance calculating point-to-point transmission, come to accompany how many wind vector row between calculating observation point and wind vector unit at the width (normally 25KM) by wind vector unit.As the distance by calculating between intersecting point coordinate and the starting point coordinate of wind vector unit, obtain distance for 110KM, so just accompany 4 wind vector row between intersection point and wind vector unit starting point, so just can determine that the line number of observation station (alternatively intersection point) is 5 (the wind vector line number+1=4+1 accompanied between line number=intersection point and wind vector unit starting point).
Similar, in step s 103, after the width of the shortest known spherical distance and wind vector unit, the wind vector row number along cross rail direction just can be calculated.Wind vector cell columns is that benchmark arranges with sub-satellite track, and usual sub-satellite track is the initial of wind vector cell columns.So known, wind vector row number namely the shortest spherical distance/wind vector cell width round business's (if the shortest spherical distance is 123, wind vector cell width is 25, and so the business that rounds of 123/25 is 4, and namely wind vector row number are 4).
In step S104, determine observation station wind vector row number and wind vector line number after, just can identify the particular location of observation station in wind vector unit face according to this line number and row number, namely use the numbering of wind vector unit (line number and row number) to identify the position at observation station place.It should be noted that the order of step S102 and step S103 can be put upside down, the execution sequence of step S102 and step S103 can not affect net result of the present invention.
Concrete, in step S101, according to the substar coordinate obtained in advance and observation station coordinate, the shortest spherical distance that calculating observation point arrives assigned address on sub-satellite track can comprise the steps, as shown in Figure 7:
S1011, obtains multiple substar coordinate, and observation station coordinate;
S1012, the respectively spherical distance of calculating observation point and each substar;
S1013, selects spherical distance minimum in multiple spherical distance as the shortest spherical distance.
In step S1011, the substar coordinate got is multiple, can be all substar coordinates that satellite scans in once complete scanning activity action.In order to simplify calculating, also can be in all substar coordinates interval get multiple substar coordinate (as 1, interval substar carries out getting a little).
In step S1012, observation station refers to be needed to use wind vector unit to express the observation station of its position.Arrive substar spherical distance in order to observation station can be used and represent the shortest spherical distance, need to calculate the distance between each substar and observation station respectively, and by spherical distance the shortest in the distance that calculates in step S1013 selecting step S1012 as the shortest spherical distance.Satellite is when scanning, its sweep frequency is quite high, and the distance namely between two substars can not be separated by far, therefore, even if the straight line at the shortest spherical distance place is not perpendicular with sub-satellite track, these two lines are also approximately perpendicular.Therefore by finding the nearest substar of a distance observation station in multiple substar, and make the length along path of nearest substar and observation station line be rational as the shortest spherical distance, more because the width of wind vector unit is 25KM, therefore, the length along path of nearest substar and observation station line can be neglected as the error of the shortest spherical distance.
On the other hand, in step S101, determine that the intersecting point coordinate of straight line corresponding to the shortest spherical distance and sub-satellite track can comprise the steps:
Determine that the coordinate of the substar corresponding to spherical distance minimum in the spherical distance of observation station and each substar is intersecting point coordinate.
In step S1012, after calculating the distance of each substar and observation station, and further, in step S1013, after determining the shortest spherical distance, can be intersecting point coordinate by the coordinate of the substar corresponding to spherical distance (the shortest spherical distance) minimum in the spherical distance of observation station and each substar.The line number of the substar corresponding to the shortest spherical distance is identical with the line number of observation station, therefore, can use the line number of the substar corresponding to the shortest spherical distance, as the line number of observation station.
Concrete, step S102, according to starting point coordinate and the default wind vector cell width of intersecting point coordinate, wind vector unit, calculating observation point can comprise the steps along the wind vector line number in straight rail direction, as shown in Figure 8:
S1021, calculates the direct spherical distance of Two coordinate point according to intersecting point coordinate and starting point coordinate;
S1022, calculates direct spherical distance and rounds business divided by first of the wind vector cell width preset;
S1023, determines that observation station is first round business+1 along the wind vector line number in straight rail direction.
In step S1021, after intersecting point coordinate on known sphere and starting point coordinate (starting point coordinate of wind vector unit), just sphere distance between two points formula can be used to calculate the direct spherical distance of the two, concrete algorithm is more, exceeds description at this to concrete computation process.
After obtaining direct spherical distance by step S1021, step S1022 just calculates direct spherical distance according to the wind vector cell width obtained in advance and rounds business divided by first of wind vector cell width.Concrete, if direct spherical distance is 158, wind vector cell width is 25, so 156/25=6, remainder 8, and first rounds business namely 6.Round business+1 by step S1023 by get first again, just can obtain the wind vector line number along straight rail direction.
Similar, when calculating wind vector row number, can make in a like fashion, namely step S103, according to the shortest spherical distance and the width with wind vector unit, calculating observation point comprises along the wind vector row number in cross rail direction:
S1031, calculate the shortest spherical distance divided by wind vector unit width second round business;
S1032, determines that observation station number is second round business+1 along the wind vector row in cross rail direction.
In step S1031, calculate to be the shortest spherical distance obtain divided by wind vector cell width second rounds business, the account form that the second account form and first rounding business rounds business is identical.Again by step S1032, rounding business+1 by second, is the wind vector row number of observation station along cross rail direction.
Wherein, step S1013,
The spherical distance of calculating observation point and each substar respectively, specifically can be subdivided into following steps, as Fig. 9:
S201, is divided into many groups by multiple continuous substar;
S202, according to the sampling rule preset, selects at least one substar as a reference point in each group;
S203, calculates the spherical distance of each reference point and observation station respectively, to determine that the group at the reference point place nearest with observation station spherical distance is for reference group;
S204, the reference spherical distance of each substar and observation station in difference computing reference group, to determine with reference to the minimum value in spherical distance for the shortest spherical distance.
In order to can when making each substar of system-computed and observation station distance, reduce calculated amount, in step S201, calculating substar and observation station apart from front, first will be all, or multiple substar is divided into groups, the substar wherein in same group is all on the time, or continuous print point on position, if any 1-10, the substar that this 10 continuous sweeps obtain, one group, 4,5,6 and 7 so can be divided into be divided into one group by 1,2 and 3,8,9 and 10 are divided into one group.Certainly, when grouping, in order to ensure the balance of data, the number of the substar that can often comprise in group is identical, or approximate identical.Also can be selectively, improve the quantity that certain or decibel improve substar in multiple groups.In actual use, the quantity of the substar in same grouping is preferably about 100.
In step S202, need to sample the substar in each grouping.When the substar number in each grouping is identical time, all can get a substar in each grouping as a reference point, the mode of getting in each grouping is a little preferably identical, as a reference point in all got first point in each grouping, or it is as a reference point all to get last point; When the substar number in each grouping is identical time, can according to the difference of substar number in each grouping, the substar choosing different quantity is as a reference point.The quantity being substar in 20, B grouping as the quantity of substar in, A grouping is 10, and two substars just can got in A grouping are as a reference point, get B divide into groups in 1 substar as a reference point.Same, when the substar number in each grouping is identical time, when getting point, in order to ensure the reference point of taking out be the identical amplitude in interval (like this, be conducive to determining reference group accurately), when selecting reference point in each grouping, before needing basis not divide into groups, after spacing between each substar determines grouping, the position of reference point.
In step S203, calculate the reference point of each taking-up and the spherical distance of observation station, and the group selecting the reference point place nearest with observation station spherical distance is reference group.
Can be determined by step S203, the grouping of the nearest substar of distance observation station may be included.In order to ensure that the grouping gone out selected by step S203 is correct, therefore, when step S202 selects reference point, select position in same grouping, or the time is upper as a reference point near middle substar.Subsequently, in step S204, after determining reference group, the spherical distance of each substar and observation station in computing reference group can be distinguished, thus can select in reference group, the distance between the substar that the spherical distance of substar and observation station is nearest and observation station is as the shortest spherical distance.
By above-mentioned computation process, only need to use the data in L1B directly to calculate the shortest spherical distance, and do not need other data to calculate, reduce the acquisition difficulty of basic data and simplify the calculated amount of data.
Concrete, step S204, the reference spherical distance of each substar and observation station in difference computing reference group, to determine can comprise the steps: for the shortest spherical distance with reference to the minimum value in spherical distance
Judge whether be less than the first default reference threshold with reference to the minimum value in spherical distance;
If so, then determine that with reference to the minimum value in spherical distance be the shortest spherical distance.
When being less than the first default reference threshold with reference to the minimum value in spherical distance, then can thinking and determine that the minimum value with reference in spherical distance is rational as the shortest spherical distance.When system calculates, should arrange reference value, in case error appears in locking system, or there is deviation in the net result caused due to the error of raw data.
Concrete, in step S101, according to the substar coordinate obtained in advance and observation station coordinate, the shortest spherical distance that calculating observation point arrives assigned address on sub-satellite track can also comprise the steps, as shown in Figure 10:
S301, calculates the spherical distance of each substar apart from initial substar respectively, and initial substar is the starting point of wind vector unit;
S302, the coordinate of the substar that chosen distance observation station is nearest is intersecting point coordinate;
S303, the spherical distance between calculating observation point and the substar nearest apart from observation station is as the shortest spherical distance.
Namely directly can calculate the spherical distance of each substar and initial substar, like this after the wind vector unit line number and row number that complete an observation station calculate, just can directly use the substar coordinate nearest with next observation station as intersecting point coordinate, this intersecting point coordinate is calculated in step S301, does not need double counting.
Except determining whether be less than the first reference with it with reference to the minimum value in spherical distance, can also ensure that the wind vector unit line number that draws is rational with row number by other judgment mode, as in step S203, calculate the spherical distance of each reference point and observation station respectively, to determine that the group at the reference point place nearest with observation station spherical distance can comprise for reference group:
Judge whether spherical distance minimum in the spherical distance of each reference point and observation station is less than the second default reference threshold;
If so, then determine that the group at the reference point place nearest with observation station spherical distance is reference group.
Concrete, the first reference threshold is 900-1100KM; Second reference threshold is 900-1100KM.Preferably, the first reference threshold is 1000KM, and the second reference threshold also can be 1000KM.
The above; be only the specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; change can be expected easily or replace, all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should described be as the criterion with the protection domain of claim.
Claims (10)
1. No. two, ocean microwave scatterometer bin matching process, is characterized in that, comprising:
According to the substar coordinate obtained in advance and observation station coordinate, calculating observation point arrives the intersecting point coordinate of the shortest spherical distance of assigned address on sub-satellite track and the straight line described in determining corresponding to the shortest spherical distance and described sub-satellite track, and described sub-satellite track is connected to form by multiple substar coordinate;
According to starting point coordinate and the default wind vector cell width of described intersecting point coordinate, wind vector unit, calculate the wind vector line number of described observation station along straight rail direction;
According to the shortest described spherical distance and the width with wind vector unit, calculate the wind vector row number of described observation station along cross rail direction;
Determine wind vector line number and wind vector row number corresponding to wind vector cell position be described observation station wind vector unit face on position.
2. No. two, ocean according to claim 1 microwave scatterometer bin matching process, is characterized in that,
The substar coordinate that described basis obtains in advance and observation station coordinate, the shortest spherical distance that calculating observation point arrives assigned address on sub-satellite track comprises:
Obtain multiple substar coordinate, and observation station coordinate;
Calculate the spherical distance of described observation station and each substar respectively;
Select spherical distance minimum in multiple described spherical distance as the shortest spherical distance.
3. No. two, ocean according to claim 2 microwave scatterometer bin matching process, is characterized in that, described determine described in the intersecting point coordinate of straight line corresponding to the shortest spherical distance and described sub-satellite track comprise:
Determine that the coordinate of the substar corresponding to spherical distance minimum in the spherical distance of described observation station and each substar is intersecting point coordinate.
4. No. two, ocean according to claim 1 microwave scatterometer bin matching process, it is characterized in that, the described starting point coordinate according to described intersecting point coordinate, wind vector unit and default wind vector cell width, calculate described observation station and comprise along the wind vector line number in straight rail direction:
The direct spherical distance of Two coordinate point is calculated according to described intersecting point coordinate and described starting point coordinate;
Calculate direct spherical distance and round business divided by first of the wind vector cell width preset;
Determine that described observation station is first round business+1 along the wind vector line number in straight rail direction.
5. No. two, ocean according to claim 1 microwave scatterometer bin matching process, is characterized in that, the shortest spherical distance described in described basis and and the width of wind vector unit, calculate described observation station and comprise along the wind vector row number in cross rail direction:
The shortest spherical distance described in calculating divided by described wind vector unit width second round business;
Determine that described observation station number is second round business+1 along the wind vector row in cross rail direction.
6. No. two, ocean according to claim 2 microwave scatterometer bin matching process, is characterized in that,
The described spherical distance calculating described observation station and each substar respectively comprises:
Multiple continuous substar is divided into many groups;
According to the sampling rule preset, select at least one substar as a reference point in each group;
Calculate the spherical distance of each reference point and described observation station respectively, to determine that the group at the reference point place nearest with described observation station spherical distance is for reference group;
Calculate the reference spherical distance of each substar and observation station in described reference group respectively, with determine described with reference to the minimum value in spherical distance for the shortest spherical distance.
7. No. two, ocean according to claim 6 microwave scatterometer bin matching process, it is characterized in that, the described reference spherical distance calculating each substar and observation station in described reference group respectively, to determine that the minimum value in described reference spherical distance comprises for the shortest spherical distance:
Judge whether the minimum value in described reference spherical distance is less than the first default reference threshold;
If so, then determine that described is the shortest spherical distance with reference to the minimum value in spherical distance.
8. No. two, ocean according to claim 1 microwave scatterometer bin matching process, is characterized in that, the substar coordinate that described basis obtains in advance and observation station coordinate, and the shortest spherical distance that calculating observation point arrives assigned address on sub-satellite track comprises:
Calculate the spherical distance of each substar apart from initial substar respectively, described initial substar is the starting point of wind vector unit;
The coordinate of the substar that chosen distance observation station is nearest is intersecting point coordinate;
Spherical distance between calculating observation point and the substar nearest apart from observation station is as the shortest spherical distance.
9. No. two, ocean according to claim 6 microwave scatterometer bin matching process, it is characterized in that, the described spherical distance calculating each reference point and described observation station respectively, to determine that the group at the reference point place nearest with described observation station spherical distance comprises for reference group:
Judge whether spherical distance minimum in the spherical distance of each reference point and described observation station is less than the second default reference threshold;
If so, then determine that the group at the reference point place nearest with described observation station spherical distance is reference group.
10. No. two, the ocean according to claim 7 or 9 microwave scatterometer bin matching process, it is characterized in that, described first reference threshold is 900-1100KM; Described second reference threshold is 900-1100KM.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410692674.2A CN104330785B (en) | 2014-11-26 | 2014-11-26 | No. two, ocean microwave scatterometer bin matching process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410692674.2A CN104330785B (en) | 2014-11-26 | 2014-11-26 | No. two, ocean microwave scatterometer bin matching process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104330785A true CN104330785A (en) | 2015-02-04 |
CN104330785B CN104330785B (en) | 2015-11-18 |
Family
ID=52405546
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410692674.2A Expired - Fee Related CN104330785B (en) | 2014-11-26 | 2014-11-26 | No. two, ocean microwave scatterometer bin matching process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104330785B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700456A (en) * | 2015-03-18 | 2015-06-10 | 广州中海达卫星导航技术股份有限公司 | Method and system for pixel screening and filtering of SeaWinds scatterometer wind vector inversion |
CN104700457A (en) * | 2015-03-20 | 2015-06-10 | 广州中海达卫星导航技术股份有限公司 | Pixel screening and filtering method and system for wind vector retrieval under precipitation situations |
CN105676191A (en) * | 2016-01-08 | 2016-06-15 | 国家卫星海洋应用中心 | Marine microwave remote sensing technology-based scatterometer data converting method and device |
CN112257269A (en) * | 2020-10-23 | 2021-01-22 | 中国船舶科学研究中心 | Method for determining relative slippage form between layers of marine flexible pipeline |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020093450A1 (en) * | 2000-07-24 | 2002-07-18 | Agence Spatiale Europeenne | Aperture synthesis radiometer and method of controlling same |
CN101697009A (en) * | 2009-10-27 | 2010-04-21 | 武汉理工大学 | Sea wave surface reduction method |
CN101853335A (en) * | 2010-06-01 | 2010-10-06 | 国家卫星海洋应用中心 | Point inversion method for ocean surface wind field inversion |
CN103698750A (en) * | 2014-01-07 | 2014-04-02 | 国家卫星海洋应用中心 | HY-2 satellite scatterometer sea surface wind field retrieval method and device |
-
2014
- 2014-11-26 CN CN201410692674.2A patent/CN104330785B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020093450A1 (en) * | 2000-07-24 | 2002-07-18 | Agence Spatiale Europeenne | Aperture synthesis radiometer and method of controlling same |
CN101697009A (en) * | 2009-10-27 | 2010-04-21 | 武汉理工大学 | Sea wave surface reduction method |
CN101853335A (en) * | 2010-06-01 | 2010-10-06 | 国家卫星海洋应用中心 | Point inversion method for ocean surface wind field inversion |
CN103698750A (en) * | 2014-01-07 | 2014-04-02 | 国家卫星海洋应用中心 | HY-2 satellite scatterometer sea surface wind field retrieval method and device |
Non-Patent Citations (2)
Title |
---|
张毅,林明森,宋清涛,解学通,邹巨洪: "海洋二号卫星微波散射计数据预处理技术研究", 《中国工程科学》 * |
邹巨洪 , 林明森, 潘德炉,陈正华, 杨乐: "QuikSCAT风矢量快速反演的后向散射系数预处理算法", 《热带海洋学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700456A (en) * | 2015-03-18 | 2015-06-10 | 广州中海达卫星导航技术股份有限公司 | Method and system for pixel screening and filtering of SeaWinds scatterometer wind vector inversion |
CN104700457A (en) * | 2015-03-20 | 2015-06-10 | 广州中海达卫星导航技术股份有限公司 | Pixel screening and filtering method and system for wind vector retrieval under precipitation situations |
CN104700457B (en) * | 2015-03-20 | 2018-09-14 | 广州中海达卫星导航技术股份有限公司 | Picture dot screening, filtering method and the system of wind vector retrieval under rainfall situation |
CN105676191A (en) * | 2016-01-08 | 2016-06-15 | 国家卫星海洋应用中心 | Marine microwave remote sensing technology-based scatterometer data converting method and device |
CN105676191B (en) * | 2016-01-08 | 2016-11-30 | 国家卫星海洋应用中心 | Scatterometer data method for transformation based on ocean microwave remote sensing technology and device |
CN112257269A (en) * | 2020-10-23 | 2021-01-22 | 中国船舶科学研究中心 | Method for determining relative slippage form between layers of marine flexible pipeline |
CN112257269B (en) * | 2020-10-23 | 2023-05-12 | 中国船舶科学研究中心 | Method for determining relative slippage morphology between ocean flexible pipeline layers |
Also Published As
Publication number | Publication date |
---|---|
CN104330785B (en) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Balestrieri et al. | Sensors and measurements for unmanned systems: An overview | |
US10962650B2 (en) | Polyhedral geofences | |
CN105975763B (en) | Fusion method and device of multisource sea surface wind field | |
CN106526551B (en) | A kind of radar antenna dynamic performance testing system and method | |
US10915673B2 (en) | Device, method, apparatus, and computer-readable medium for solar site assessment | |
CN102495420B (en) | Underwater object precision positioning system and method | |
US5894450A (en) | Mobile underwater arrays | |
CN104330785B (en) | No. two, ocean microwave scatterometer bin matching process | |
Chen et al. | Review of AUV underwater terrain matching navigation | |
CN101246216A (en) | High-speed laser ranging system including a fiber laser | |
CN105115492B (en) | Underwater terrain matching navigation system based on acoustic Doppler log | |
CN107479045B (en) | Method and system for eliminating short vegetation based on full-waveform laser radar point cloud data | |
CN104407338A (en) | Chinese HY-2 satellite microwave scatterometer-based polar sea ice identification method | |
CN112558102A (en) | Airborne oblique laser three-dimensional measurement and composite imaging system and use method thereof | |
CN109319074B (en) | Multi-orthogonal signal emission unmanned submersible vehicle sound guiding and recycling system | |
CN104297743B (en) | Method and device for eliminating distance measuring ambiguity of high repetition frequency airborne laser radar system | |
CN103969639A (en) | Signal processing system and method of five-wave-beam fish finder | |
KR20210046783A (en) | System and method for determining the geographic-location of a signal source | |
CN115755071A (en) | Deep sea in-situ fine detection frame design method based on acousto-optic remote sensing and VR technology | |
Diaz et al. | The Bathy-Drone: An autonomous uncrewed drone-tethered sonar system | |
CN104897059A (en) | Portable irregular accumulation body volume measurement method | |
CN103777196B (en) | Based on terrain object distance single station measuring method and the measuring system thereof of geography information | |
CN206209109U (en) | Three-dimensional laser radar system | |
Grządziel | Application of remote sensing techniques to identification of underwater airplane wreck in shallow water environment: Case study of the Baltic Sea, Poland | |
CN116522548B (en) | Multi-target association method for air-ground unmanned system based on triangular topological structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CB03 | Change of inventor or designer information |
Inventor after: Zou Juhong Inventor after: Lin Mingsen Inventor after: Zou Bin Inventor after: Zeng Tao Inventor after: Cui Songxue Inventor before: Zou Juhong |
|
COR | Change of bibliographic data | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151118 Termination date: 20161126 |
|
CF01 | Termination of patent right due to non-payment of annual fee |