A new CFAR ship target detection method in SAR imagery

Many ship target detection methods have been developed since it was verified that ship could be imaged with the space-based SAR systems. Most developed detection methods mostly emphasized ship detection rate but not computation time. By making use of the advantages of the K-distribution CFAR method and two-parameter CFAR method, a new CFAR ship target detection algorithm was proposed. In that new method, we use the K-distribution CFAR method to calculate a global threshold with a certain false-alarm rate. Then the threshold is applied to the whole SAR imagery to determine the possible ship target pixels, and a binary image is given as the preliminary result. Mathematical morphological filter are used to filter the binary image. After that step, we use the two-parameter CFAR method to detect the ship targets. In the step, the local sliding window only works in the possible ship target pixels of the SAR imagery. That step avoids the statistical calculation of the background pixels, so the method proposed can much improve the processing speed. In order to test the new method, two SAR imagery with different background were used, and the detection result shows that that method can work well in different background circumstances with high detection rate. Moreover, a synchronous ship detection experiment was carried out in Qingdao port in October 28, 2005 to verify the new method and one ENVISAT ASAR imagery was acquired to detect ship targets. It can be concluded from the experiment that the new method not only has high detection rate, but also is time-consuming, and is suitable for the operational ship detection system.     

Numerical simulation of non-gaussian process of wind waves

In this paper further mathematical analysis on "correlation transfer technique" by Polge el al. is carried out, the tenable conditions and the extent of suitability for the said method are proved as well. In consideration of the influence of skewness of the sea surface elevation on spectral shape, a "quasi-correlation transfer techique" is developed by the modification of the simulated target spectrum. Meanwhile, the numerical simulation of the non-Gaussian process of wind waves is carried out in view of the two conditions of the surface elevation probability distribution and the spectrum. By using its simulated results, the influence of skewness of the sea surface elevation on two parameters in the distribution of wave heights (which had been fitted by using the Weibull distribution) is analysed. The "quasi- correlation tranfer technique" is verified and compared with the selection wave data observed in the Jiaozhou Bay in the period of 1980 to 1981. Results make clear (hat, as far as the statistical d     

Response Characteristics of Load on Vessels in Waves

Considering the requirement of direct design and fatigue test for ships and floating structures byuse of FEM technique,a computational procedure of spectral analysis for wave load on the hull surface is de-veloped in this paper.The response of hydrodynamic pressure on the body surface to a designated sea state forships and floating structures is calculated by use of the revised strip method with the hull bound perturbationflow concept.The spectral function of wave load for the defined point on the body surface can be determinedfrom the Wiener-Khinhin theorem and the characteristic load value can be also obtained from spectral mo-ment analysis.A container ship is taken as a computational example and the sample of wave load with a cer-tain probability and corresponding encountered frequency is provided.     

Homotopy Method for Inverse Design of the Bulbous Bow of A Container Ship

The homotopy method is utilized in the present inverse hull design problem to minimize the wave-making coefficient of a 1300 TEU container ship with a bulbous bow. Moreover, in order to improve the computational efficiency of the algorithm, a properly smooth function is employed to update the homotopy parameter during iteration. Numerical results show that the homotopy method has been successfully applied in the inverse design of the ship hull. This method has an advantage of high performance on convergence and it is credible and valuable for engineering practice.     

NON-LINEAR STOCHASTIC RESPONSE AND ITS EXTREME VALUE OF JACK UP PLATFORMS

The non-linear stochastic response of a jack-up platform subjected to wave load has been analyzed dynamically in this paper, and the analysis method in time domain is considered. Monte Carlo simulation is used to generate random sea. An emphasis is placed on the nonlinear hydrodynamic force. Several distributions for the statistical estimation of extreme responses are compared. For Gumble distribution, the parameters of its asymptotic distribution expression have been checked. The results show that the Gumbel distribution agrees well with the simulated values of the responses.     

Study on the Configuration of Towed Flexible Cables

Based on the fundamental equation of flexible cable dynamics for a towed system, an easily solved mathematical model is set up in this paper by means of appropriate simplification. Several regular patterns of spatial motion of towed flexible cables in water are obtained through numerical simulation with the finite difference method, and then modification and verification by trial results at sea. A technical support is provided for the towing ship to maneuver properly when a flexible cable is towed. Furthermore, the relations between two towed flexible cables, which are towed simultaneously by a ship, are investigated. The results show that the ship towing two flexible cables is safe under the suggested arrangement of two winches for the towing system, and the coiling/uncoiling sequences of the cables as well as the suggested way of maneuvering.     

A maximum-entropy compound distribution model for extreme wave heights of typhoon-affected sea areas

A new compound distribution model for extreme wave heights of typhoon-affected sea areas is proposed on the basis of the maximum-entropy principle.The new model is formed by nesting a discrete distribution in a continuous one,having eight parameters which can be determined in terms of observed data of typhoon occurrence-frequency and extreme wave heights by numerically solving two sets of equations derived in this paper.The model is examined by using it to predict the N-year return-period wave height at two hydrology stations in the Yellow Sea,and the predicted results are compared with those predicted by use of some other compound distribution models.Examinations and comparisons show that the model has some advantages for predicting the N-year return-period wave height in typhoon-affected sea areas.     

A coastal ocean model of semi-implicit finite volume unstructured grid

A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.     

Probability distribution of surface slope for nonlinear random sea waves and calculation of whitecap coverage

Based on the nonlinear model of two-dimensional random sea waves, a statistical distribution of wave surface slope exact to the third order is derived by using the expansion of the characteristic function and direct calculations of each order moment. Based on the distribution of wave surface slope derived in this paper, a whitecap coverage is proposed by using the limit surface slope as a criterion of wave breaking. The whitecap coverage expressed by the model depends on three parameters which can be determined in principle by the linear wave spectrum and three kinds of wave-wave interaction.     

Shallow sea topography detection using fully polarimetric Gaofen-3 SAR data based on swell patterns

Compared to single-polarization synthetic aperture radar(SAR) data, fully polarimetric SAR data can provide more detailed information of the sea surface, which is important for applications such as shallow sea topography detection. The Gaofen-3 satellite provides abundant polarimetric SAR data for ocean research. In this paper, a shallow sea topography detection method was proposed based on fully polarimetric Gaofen-3 SAR data. This method considers swell patterns and only requires SAR data and ...     

基于模型相似度拟合的海杂波统计方法

本文提出一种基于模型相似度拟合的海杂波统计方法。首先根据合成孔径雷达(SAR)图像计算瑞利分布、对数正态分布、韦布尔分布、K分布、G0分布5种经典的海杂波分布的概率密度函数,然后根据模型间的相似度准则拟合得到新的海杂波分布模型。文章利用四景不同类型的真实SAR数据对算法的拟合性能进行了评价,结果显示利用该算法得到的拟合模型与真实SAR数据的平均Kullback-Leibler距离仅为0.015 84,远优于其他分布模型。基于该拟合模型的恒虚警率舰船检测算法对四景SAR数据的平均检测精度达到95.75%,在控制虚警和漏检方面均优于采用其他模型的同类方法。     

一种新的用于极化SAR图像船只检测的散射相似性测度

本文提出了一种新的用于极化合成孔径雷达（SAR）海上船只目标检测的测度。首先利用散射相似性参数研究船只与海杂波散射机制的差异。基于这些差异，提出了一种新的船只检测测度，该测度能够有效区分船只目标和海杂波。然后提出了利用核密度估计（KDE）方法对该测度进行建模的方法。基于统计模型，实现了自适应恒虚警率（CFAR）的检测方案。本文利用多景C波段RADARSAT-2极化SAR（Pol-SAR）数据上进行实验，系统分析了所提出测度的海杂波拟合性能与船只检测性能，并与两个经典的极化测度进行了比较，实验和比较结果证明了所提的测度的可行性。     

基于AIS距离-多普勒投影的地波雷达CFAR检测验证方法

恒虚警(CFAR)检测是地波雷达船只目标探测的主要方法。目前基于船舶自动识别系统(AIS)信息的CFAR检测验证方法是一种间接验证方式,容易受地波雷达系统测向误差的影响,且不具备对错检与漏检目标的分析能力。本文提出了一种基于AIS信息的评价地波雷达CFAR检测结果的直接验证方法。该方法将有效的AIS信息转换到地波雷达的距离-多普勒谱中,通过在该谱中AIS信息与CFAR检测结果的关联分析,实现CFAR检测结果的直接评价。论文首先给出了方法的原理和处理流程,然后利用实测数据给出了该方法在CFAR检测结果评价中的实际应用,验证了方法有效性,而且该方法也为低可观测目标的CFAR检测提供了参数优化调整的依据。     

基于一维综合孔径微波辐射计的海面温度反演研究

Due to the low spatial resolution of sea surface temperature(T_S) retrieval by real aperture microwave radiometers,in this study, an iterative retrieval method that minimizes the differences between brightness temperature(T_B)measured and modeled was used to retrieve sea surface temperature with a one-dimensional synthetic aperture microwave radiometer, temporarily named 1 D-SAMR. Regarding the configuration of the radiometer, an angular resolution of 0.43° was reached by theoretical calculation. Experiments on sea surface temperature retrieval were carried out with ideal parameters; the results show that the main factors affecting the retrieval accuracy of sea surface temperature are the accuracy of radiometer calibration and the precision of auxiliary geophysical parameters. In the case of no auxiliary parameter errors, the greatest error in retrieved sea surface temperature is obtained at low T_S scene(i.e., 0.710 6 K for the incidence angle of 35° under the radiometer calibration accuracy of0.5 K). While errors on auxiliary parameters are assumed to follow a Gaussian distribution, the greatest error on retrieved sea surface temperature was 1.330 5 K at an incidence angle of 65° in poorly known sea surface wind speed(W)(the error on W of 1.0 m/s) over high W scene, for the radiometer calibration accuracy of 0.5 K.     

Estimation of standing stock of chlorophylla and primary production from remote-sensed ocean color in the Oyashio region, the Western subarctic Pacific, during the spring bloom in 1997

In order to examine the applicability of remotely-sensed ocean color for the estimation of phytoplankton biomass and primary production in the Oyashio region, the western subarctic Pacific, vertical distributions of chlorophylla concentration and primary production were observed in April and May 1997. Spring bloom was observed in both April and May, and the surface concentration of chlorophylla exceeded 40 mg m⁻³. The relationship between the standing stocks of chlorophylla within the layer from the sea surface to one optical depth (0–1/k layer) and the surface chlorophylla concentration is expressed as a Michaelis-Menten equation. The mean ratio of the standing stock of chlorophylla in the euphotic layer to that in the 0–1/k layer was 4.41, this ratio did not significantly differ from 4.61 which was obtained at homogeneous distribution of chlorophylla within the euphotic layer. These facts suggest that the distribution of chlorophylla could be assumed to be homogeneous in the euphotic layer during the spring bloom. Results of primary production measurements by simulatedin situ method were compared with those by an algorithm with two variables; chlorphylla and non-spectral PAR. Daily primary production in the euphotic layer estimated by the algorithm varied in a range of 38–274% of that estimated by incubation, although the primary productions by the algorithm agreed with those by the incubation at a half of stations. Primary production within the euphotic layer calculated using simply the surface data was the same as that estimated using vertical distribution of chlorophylla. These results show that the primary production in the euphotic layer may be estimated from the remote sensed measurements during the spring bloom in the Oyashio region.     

Habitat suitability index of Chub mackerel (<Emphasis Type="Italic">Scomber japonicus</Emphasis>) from July to September in the East China Sea

The habitat quality of Chub mackerel (Scomber japonicus) in the East China Sea has been a subject of concern in the last 10 years due to large fluctuations in annual catches of this stock. For example, the Chinese light-purse seine fishery recorded 84000 tons in 1999 compared to 17000 tons in 2006. The fluctuations have been attributed to variability in habitat quality. The habitat suitability Index (HSI) has been widely used to describe fish habitat quality and in fishing ground forecasting. In this paper we use catch data and satellite derived environmental variables to determine habitat suitability indices for Chub mackerel during July to September in the East China Sea. More than 90% of the total catch was found to come from the areas with sea surface temperature of 28.0°–29.4°C, sea surface salinity of 33.6–34.2 psu, chlorophyll-a concentration of 0.15–0.50 mg/m³ and sea surface height anomaly of −0.1–1.1 m. Of the four conventional models of HSI, the Arithmetic Mean Model (AMM) was found to be most suitable according to Akaike Information Criterion analysis. Based on the estimation of AMM in 2004, the monthly HSIs in the waters of 123°–125°E and 27°30′–28°00′ N were more than 0.6 during July to September, which coincides with the catch distribution in the same time period. This implies that AMM can yield a reliable prediction of the Chub mackerel’s habitat in the East China Sea.     

高精度新方法测量计算船舶磁场分布

将高精度质子磁探仪和位场变换理论用于测量计算船舶磁场，以研究其在海洋中的分布。由质子磁探仪实测海面上１０ｍ处磁场的模（船的磁场与地磁场矢量的和的绝对值），用位场变换理论计算出船的磁场三个分量在水面和水下的分布；并采用一系列有效措施，提高计算速度和结果的精度。将该方法和磁球理论严格解比较，在ｚ＝０面上误差小于０．１％，在ｚ＝－２０ｍ面上误差小于１％；还给出船的磁场强度在船的纵横向的分布曲线和平面立体     

Mixed layer depth front and subduction of low potential vorticity water in an idealized ocean GCM

It is known that there is a front-like structure at the mixed layer depth (MLD) distribution in the subtropical gyre, which is called the MLD front, and is associated with the formation region of mode water. In the present article, the generation mechanism of the MLD front is studied using an idealized ocean general circulation model with no seasonal forcing. First, it is shown that the MLD front occurs along a curve where u _g ·∇T _s = 0 is satisfied (u _g is the upper ocean geostrophic velocity vector, T _s is the sea surface temperature and ∇ is the horizontal gradient operator). In other words, the front is the boundary between the subduction region (u _g ·∇T _s > 0) and the region where subduction does not occur (u _g ·∇T _s < 0). Second, we have investigated subduction of low potential vorticity water at the MLD front, which has been pointed out by past studies. Since u _g ·∇T _s = 0 at the MLD front, the water particles do not cross the outcrop at the MLD front. The water that is subducted at the MLD front has come from the deep mixed layer region where the sea surface temperature is higher than that at the MLD front. The temperature of the water in the deep mixed layer region decreases as it is advected eastward, attains its minimum at the MLD front where u _g ·∇T _s = 0, and then subducts under the warmer surface layer. Since the deep mixed layer water subducts beneath a thin stratified surface layer, maintaining its thickness, the mixed layer depth changes abruptly at the subduction location.     

海洋二号微波散射计极地海冰识别算法研究

An Antarctic sea ice identification algorithm on the HY-2A scatterometer(HSCAT) employs backscattering coefficient(σ0) and active polarization ratio(APR) for a preliminary sea ice identification.Then standard deviation(STD) filtering and space filtering are carried out.Finally,it is used to identify sea ice.A process uses a σ0,STD threshold and an APR as sea ice indicators.The sea ice identification results are verified using the sea ice distribution data of the ASMR2 released by the National Snow and Ice Data Center as a reference.The results show very good consistence of sea ice development trends,seasonal changes,area distribution,and sea ice edge distribution of the sea ice identification results obtained by this algorithm relative to the ASMR2 sea ice results.The accuracy of a sea ice coverage is 90.8% versus the ASMR2 sea ice results.This indicates that this algorithm is reliable.