CN106501790B - A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method - Google Patents

A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method Download PDF

Info

Publication number
CN106501790B
CN106501790B CN201610971539.0A CN201610971539A CN106501790B CN 106501790 B CN106501790 B CN 106501790B CN 201610971539 A CN201610971539 A CN 201610971539A CN 106501790 B CN106501790 B CN 106501790B
Authority
CN
China
Prior art keywords
water
target
sea
radiance
submarine target
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
Application number
CN201610971539.0A
Other languages
Chinese (zh)
Other versions
CN106501790A (en
Inventor
谷延锋
夏峙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Heilongjiang Industrial Technology Research Institute Asset Management Co ltd
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201610971539.0A priority Critical patent/CN106501790B/en
Publication of CN106501790A publication Critical patent/CN106501790A/en
Application granted granted Critical
Publication of CN106501790B publication Critical patent/CN106501790B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/4802Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method, the present invention relates to shallow sea submarine target EO-1 hyperion parameter characteristic extracting methods.The deficient in stability when present invention is to solve the problems, such as the prior art directly using full spectrum or single wavelength reflectivity as feature.Step of the present invention is:Step 1:Submarine target is established from water radiance simulation model;Step 2:Determining makes target from water radiance simulation value and the submarine target of actual measured value error minimum from water radiance simulation model input parameter;Step 3:Target state estimator depth in the input parameter that step 2 is determined is with reflectivity as the high-spectrum remote-sensing feature of target.The present invention proposes that receiving signal by spectrometer estimates that depth related with submarine target attribute and reflectivity are used as feature, improves the stability of feature, improves submarine target EO-1 hyperion detectivity while so that target signature is had physical significance.The present invention is applied to optical remote sensing field of signal processing.

Description

A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method
Technical field
The present invention relates to shallow sea submarine target EO-1 hyperion parameter characteristic extracting methods.
Background technology
Shallow Water Target is detected has important meaning with identification for marine monitoring, marine resources development and maritime search and rescue etc. Justice, but neritic environment Underwater Target Detection is currently a problem.Traditional detecting devices, such as radar and sonar, in Detecting Shallow Sea It is much limited when submarine target:Although radar can form target very high spatial resolution, electromagnetic wave can not be worn Pervious bed;Sonar can utilize the submarine target of acoustic resonance identification unlike material, but by the shadow of the strong interface reverberation in shallow sea It rings, detection performance is had a greatly reduced quality.In addition to these reasons, with the development of material, propulsion and information technology, modern subaqueous target All have certain electromagnetism stealth and Sound stealth ability, this further improves the difficulty of Underwater Target Detection.Therefore, develop Optical detection means, especially hyperspectral technique become the new solution route of current neritic environment Underwater Target Detection.
Hyperspectral technique has due to the main problem that Underwater Target Detection faces at present:1) marine background radiation spectrum easily by Marine Environment Factors are influenced and are changed;2) light round trip in water layer is propagated, and target scattering spectrum will produce uncertain The radiance absolute value difference of distortion, target and marine background is smaller.These problems cause using full spectrum or single wave band Deficient in stability when reflectivity is radiated from water radiation with marine background as feature differentiation submarine target.
Invention content
The present invention is to lack when solving the prior art directly using full spectrum or single wavelength reflectivity as feature The problem of stability, and a kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method proposed.
A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method is realized according to the following steps:
Step 1:Submarine target is established from water radiance simulation model;
Step 2:Determining makes target from water radiance simulation value and the submarine target of actual measured value error minimum from water spoke Penetrate rate simulation model input parameter;
Step 3:Target state estimator depth in the input parameter that step 2 is determined is with reflectivity as the EO-1 hyperion of target Remote sensing features.
Invention effect:
In view of the problems of the existing technology, the present invention proposes to receive signal estimation and submarine target attribute by spectrometer Related depth and reflectivity improve the stability of feature, carry as feature while so that target signature is had physical significance High submarine target EO-1 hyperion detectivity lacks steady when solving directly using full spectrum or single wavelength reflectivity as feature Qualitative problem.
Description of the drawings
Fig. 1 is to radiate simulation model technical solution figure from water;
Fig. 2 is submarine target from water radiation parameter feature extraction conceptual scheme;
Fig. 3 is the albedo image for emulating the detection of submarine target EO-1 hyperion in 530nm;
Fig. 4 is target sample and distribution map of the marine background sample in two-dimensional feature space under artificial water.
Specific implementation mode
Specific implementation mode one:A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method includes the following steps:
Step 1:Submarine target is established from water radiance simulation model;
Step 2:Determining makes target from water radiance simulation value and the submarine target of actual measured value error minimum from water spoke Penetrate rate simulation model input parameter;
Step 3:Target state estimator depth in the input parameter that step 2 is determined is with reflectivity as the EO-1 hyperion of target Remote sensing features.
Submarine target from water radiate simulation model technical solution as shown in Figure 1, sea from water radiance by Absorption by Sea Water Rate, reflectivity of sea water and sub-bottom reflection rate determine that Absorption by Sea Water rate in water in phytoplankton chlorophyll concentration, water by being dissolved with Machine object absorptivity and full sea water absorptivity determine that reflectivity of sea water is determined by suspended particulate reflectivity in water and full sea water reflectivity Fixed, sub-bottom reflection rate is determined by submarine sedimentary strata material;Error of the sea from water radiance simulation value and actual measured value is calculated, Full spectrum is divided into two wavelength bands according to the spectral window of seawater, is 400-600nm and 600-715nm respectively, is introduced flat The radiance absolute value difference of two wavelength bands of equalization parameter that weigh;Search makes submarine target from water radiance simulation value and reality The input parameter of measurement error minimum, extraction reflectivity and high-spectrum remote-sensing feature of the target depth as submarine target.
Specific implementation mode two:The present embodiment is different from the first embodiment in that:Water is established in the step 1 Detailed process of the lower target from water radiance simulation model be:
Sea determines from water radiance by Absorption by Sea Water rate, reflectivity of sea water and sub-bottom reflection rate, and light is transmitted through sea Can be affected when water-air interface, interface upper surface from water radiance RrsWith the uplink radiation rate r of lower surfacersHave following Relationship:
Wherein rrsSemi-analytical model be:
Wherein a is total absorptivity, bbFor total back reflection rate, θwFor the zenith angle of the sun, H is depth, and ρ is sub-bottom reflection Rate,B is the sub-bottom reflection rate at 550nm,It is using B as the seabed being normalized with reference to value Reflectance curve;
The total absorptivity of seawater dissolves organic in the absorption of chlorophyll and water contained by the absorption of full sea water, phytoplankton The absorption three of object codetermines, i.e.,:
A=aw+aφ+ag (3)
Wherein awFor the absorptivity of full sea water, obtained by tabling look-up;aφFor the absorptivity of phytoplankton chlorophyll, agFor water The absorptivity of middle dissolved organic matter, computational methods are respectively:
aφ(λ)=[a0(λ)+a1(λ)ln(P)]P (4)
P=aφ(440)=0.06 [chl-a]0.65 (5)
ag(λ)=Gexp [- 0.015 (λ -440)] (6)
Wherein λ represents wavelength, G=ag(440), i.e., dissolved organic matter absorptivity in the water at wavelength 440nm;
The total reflectivity b of seawaterbBy the backscattering ratio of suspended particulate in the backscattering ratio and water of two full sea waters, I.e.:
bb=bbw+bbp (7)
Wherein bbwFor the reflectivity of full sea water, obtained by tabling look-up;bbpIt is the backscattering ratio of suspended particulates in water, meter Calculating formula is:
bbp(λ)=[chl-a]0.62550/λ (8)。
Other steps and parameter are same as the specific embodiment one.
Specific implementation mode three:The present embodiment is different from the first and the second embodiment in that:In the step 2 really So that target is inputted from water radiance simulation model from water radiance simulation value and the submarine target of actual measured value error minimum surely The detailed process of parameter is:
Fig. 2 is submarine target from water radiation parameter feature extraction scheme schematic diagram.According to from water radiance simulation model, The absorptivity a of full sea waterwWith reflectivity bbwIt is to determine, marine phytoplankton absorptivity aφ, suspended particulate reflectivity b in waterbp It is related with water Determination of Chlorophyll density chl-a, dissolved organic matter absorptivity a in watergHave with its reference absorptivity G at 440nm It closes, sub-bottom reflection rate ρ is related with its reference reflectivity B at 550nm.From water radiance simulation model, there are four input ginsengs altogether Number, i.e.,:
Rrs(λ)=F (aw(λ),bbw(λ),chl-a,G,B,H) (9)
It brings different parameters into sea and obtains simulation result of the sea from water radiance from the semi-analytical model that water radiates, By the error of computer sim- ulation result and actual measured results, corresponding parameter value when error minimum is obtained;Sea is radiated from water Shown in the simulation result of rate and the error of actual measured results such as formula (10):
Wherein, Rrs(λ) is actual measured results of the sea from water radiance,It is emulation of the sea from water radiance As a result;It represents to wave band λ12Radiance in range is averaged, and parameter ω is to balance from water radiance in λ12 With λ34Strength difference in range;
Search for formula (11) it is minimum from water radiance simulation model input parameter, the seabed at 550nm therein is anti- Rate B and depth H are penetrated as target signature.
Other steps and parameter are the same as one or two specific embodiments.
Specific implementation mode four:Unlike one of present embodiment and specific implementation mode one to three:Described is 400- 600nm, λ34For 600-710nm.
Other steps and parameter are identical as one of specific implementation mode one to three.
Embodiment one:
Fig. 3 is the albedo image for emulating the detection of submarine target EO-1 hyperion in 530nm, from marine background radiation areas and water Respectively 50 pixels of random selection pass through above-mentioned flow as marine background sample and submarine target sample for lower target region The estimation reflectivity of each pixel and estimating depth are extracted as the corresponding feature of the pixel, Fig. 4 is submarine target sample and sea Foreign background sample is in estimation reflectivity and the distribution in the two-dimensional feature space of estimating depth composition, the estimation reflection of submarine target Rate is all higher than marine background sample with estimating depth, and the two has different and respective metastable distributed area in feature space Domain has reached the purpose of the present invention to realize the high-spectrum remote-sensing feature extraction of submarine target stabilization.

Claims (3)

1. a kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method, which is characterized in that the shallow sea submarine target is high Spectrum parameter feature extracting method includes the following steps:
Step 1:Submarine target is established from water radiance simulation model;
Step 2:Determining makes target from water radiance simulation value and the submarine target of actual measured value error minimum from water radiance Simulation model input parameter;
Step 3:Target state estimator depth in the input parameter that step 2 is determined is with reflectivity as the high-spectrum remote-sensing of target Feature;
Detailed process of the submarine target from water radiance simulation model is established in the step 1 is:
Sea determines from water radiance by Absorption by Sea Water rate, reflectivity of sea water and sub-bottom reflection rate, and light is transmitted through seawater-sky It can be affected when vapor interface, the upper surface of seawater-Air Interface is from water radiance RrsWith the uplink radiation rate r of lower surfacersHave Following relationship:
Wherein rrsSemi-analytical model be:
Wherein a is total absorptivity, bbFor total back reflection rate, θwFor the zenith angle of the sun, H is depth, and ρ is sub-bottom reflection rate,B is the sub-bottom reflection rate at 550nm,It is anti-using B as the seabed being normalized with reference to value Penetrate rate curve;
Dissolved organic matter in the total absorptivity of seawater absorption of chlorophyll and water contained by the absorption of full sea water, phytoplankton It absorbs three to codetermine, i.e.,:
A=aw+aφ+ag (3)
Wherein awFor the absorptivity of full sea water, aφFor the absorptivity of phytoplankton chlorophyll, agFor the absorption of dissolved organic matter in water Rate, computational methods are respectively:
aφ(λ)=[a0(λ)+a1(λ)ln(P)]P (4)
P=aφ(440)=0.06 [chl-a]0.65 (5)
ag(λ)=Gexp [- 0.015 (λ -440)] (6)
Wherein λ represents wavelength, G=ag(440), i.e., dissolved organic matter absorptivity in the water at wavelength 440nm;
The total reflectivity b of seawaterbBy the backscattering ratio of suspended particulate in the backscattering ratio and water of two full sea waters, i.e.,:
bb=bbw+bbp (7)
Wherein bbwFor the reflectivity of full sea water, bbpIt is the backscattering ratio of suspended particulates in water, calculation formula is:
bbp(λ)=[chl-a]0.62550/λ (8)。
2. a kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method according to claim 1, which is characterized in that Being determined in the step 2 makes target be radiated from water from water radiance simulation value and the submarine target of actual measured value error minimum The detailed process of rate simulation model input parameter is:
According to from water radiance simulation model, the absorptivity a of full sea waterwWith reflectivity bbwIt is to determine, is emulated from water radiance There are four input parameters altogether for model, i.e.,:
Rrs(λ)=F (aw(λ), bbw(λ), chl-a, G, B, H) (9)
It brings different parameters into sea and obtains simulation result of the sea from water radiance from the semi-analytical model that water radiates, pass through The error of computer sim- ulation result and actual measured results obtains corresponding parameter value when error minimum;Sea is from water radiance Shown in the error of simulation result and actual measured results such as formula (10):
Wherein, Rrs(λ) is actual measured results of the sea from water radiance,It is emulation knot of the sea from water radiance Fruit;It represents to wave band λ12Radiance in range is averaged, and parameter ω is to balance from water radiance in λ12With λ34Strength difference in range;
Search for formula (11) it is minimum from water radiance simulation model input parameter, by the sub-bottom reflection rate B at 550nm therein And depth H is as target signature.
3. a kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method according to claim 2, which is characterized in that The λ12For 400-600nm, λ34For 600-710nm.
CN201610971539.0A 2016-11-04 2016-11-04 A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method Active CN106501790B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610971539.0A CN106501790B (en) 2016-11-04 2016-11-04 A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610971539.0A CN106501790B (en) 2016-11-04 2016-11-04 A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method

Publications (2)

Publication Number Publication Date
CN106501790A CN106501790A (en) 2017-03-15
CN106501790B true CN106501790B (en) 2018-10-02

Family

ID=58323118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610971539.0A Active CN106501790B (en) 2016-11-04 2016-11-04 A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method

Country Status (1)

Country Link
CN (1) CN106501790B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110749861B (en) * 2019-10-24 2022-06-17 哈尔滨工程大学 Three-dimensional positioning method for underwater fixed target based on multiple assumed depths
CN111561916B (en) * 2020-01-19 2021-09-28 自然资源部第二海洋研究所 Shallow sea water depth uncontrolled extraction method based on four-waveband multispectral remote sensing image

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103714341B (en) * 2014-01-21 2016-09-28 北京航空航天大学 Hyper spectral reflectance data light spectrum signature extracting method based on overall situation sensitivity analysis
CN105865624B (en) * 2016-03-22 2018-02-09 清华大学 The spectrum extracting method and device of hyper-spectral data gathering system

Also Published As

Publication number Publication date
CN106501790A (en) 2017-03-15

Similar Documents

Publication Publication Date Title
US9639756B2 (en) Visibility of targets on and below the ocean surface
CN105954732B (en) A kind of the sediment reflectivity extracting method and system of laser radar depth measurement data
Liu et al. Wave height estimation from shipborne X-band nautical radar images
JP5501176B2 (en) Ship detection apparatus, ship detection program, and ship detection method for ship detection apparatus
CN106501790B (en) A kind of shallow sea submarine target EO-1 hyperion parameter characteristic extracting method
Romeiser et al. Advanced remote sensing of internal waves by spaceborne along-track InSAR—A demonstration with TerraSAR-X
Xue et al. Airborne optical polarization imaging for observation of submarine Kelvin wakes on the sea surface: Imaging chain and simulation
Negaharipour On 3-D scene interpretation from FS sonar imagery
Wang et al. Using multibeam backscatter strength to analyze the distribution of manganese nodules: A case study of seamounts in the Western Pacific Ocean
CN110554377A (en) Single-channel SAR two-dimensional flow field inversion method and system based on Doppler center offset
Greidanus Sub-aperture behavior of SAR signatures of ships
CN112836638B (en) High sea condition SAR ship detection method and application
Liu et al. Shadowing-analysis-based wave height measurement from ship-borne X-band nautical radar images
CN106908782A (en) The extracting method in the wave propagation direction based on buoyance continuous imaging system
Kim et al. Development of real-time HABs detection technique using unmanned aerial vehicle (UAV)
Moszyński et al. Reconstruction of 3D shape from sidescan sonar images using shape from shading technique
Shen et al. An alternative method for surface current extraction from X-band marine radar images
Woock et al. Approaches to acoustic and visual underwater sensing
Xu et al. Angular response classification of multibeam sonar based on multi-angle interval division
Rostovtseva et al. Application of shipborne three-channel hyperspectrometer for remote sensing of coastal marine areas and internal seas
Han et al. Acoustic multipath structure in direct zone of deep water and bearing estimation of tow ship noise of towed line array
Love et al. Target strength distributions of Pacific sardine schools: Model results at 500 Hz to 10 kHz
Wu 2.1 Principle of Multi-beam Sounding Technology
Gong Remote sensing of marine life and submerged target motions with ocean waveguide acoustics
Gruber et al. Scattering measurements of rocky seafloors using a split-beam echosounder

Legal Events

Date Code Title Description
C06 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
TR01 Transfer of patent right

Effective date of registration: 20210119

Address after: Building 9, accelerator, 14955 Zhongyuan Avenue, Songbei District, Harbin City, Heilongjiang Province

Patentee after: INDUSTRIAL TECHNOLOGY Research Institute OF HEILONGJIANG PROVINCE

Address before: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin

Patentee before: HARBIN INSTITUTE OF TECHNOLOGY

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230410

Address after: 150027 Room 412, Unit 1, No. 14955, Zhongyuan Avenue, Building 9, Innovation and Entrepreneurship Plaza, Science and Technology Innovation City, Harbin Hi tech Industrial Development Zone, Heilongjiang Province

Patentee after: Heilongjiang Industrial Technology Research Institute Asset Management Co.,Ltd.

Address before: Building 9, accelerator, 14955 Zhongyuan Avenue, Songbei District, Harbin City, Heilongjiang Province

Patentee before: INDUSTRIAL TECHNOLOGY Research Institute OF HEILONGJIANG PROVINCE

TR01 Transfer of patent right