Trade-off of resolution and variance computed from ensembles of solutions, with application to Markov Chain Monte Carlo methods

SUMMARY We consider the case where the ‘solution’ to an inverse problem is an ensemble (e.g. drawn from the conditional probability density function $p( {{{\bf m}}|{{{\bf d}}^{obs}}} )$ of M model parameters ${{\bf m}}$ given observed data ${{{\bf d}}^{obs}}$). Here we presume that the ${{\bf m}}$s have a natural ordering, say in position x, so that ‘resolution’ means the ability of the inverse problem to distinguish physically adjacent model parameters. The trade-off curve for resolution and variance is constructed using the following steps: (1) the single solution ${{{\bf m}}^{est}}$ and its covariance ${{{\bf C}}_m}$ are estimated as the ensemble mean and covariance; (2) the eigenvalue decomposition ${{{\bf C}}_m} = {\rm{\ }}{{{\bf V} {\boldsymbol \Lambda} }}{{{\bf V}}^{\rm{T}}}$ is computed and the submatrix ${{{\boldsymbol \Lambda }}^{( N )}}$ of the N smallest eigenvalues, and submatrix ${{{\bf V}}^{( N )}}$of the N corresponding eigenvectors, are formed; (3) the equation ${{{\boldsymbol \mu }}^{( N )}} = {{{\boldsymbol \Phi }}^{( N )}}\ {{\bf m}}$ with ${{{\boldsymbol \mu }}^{( N )}} = [ {{{{\bf V}}^{( N )}}} ]{\boldsymbol{\ }}{{{\bf m}}^{est}}$and${{{\boldsymbol \Phi }}^{( N )}} = {[ {{{{\bf V}}^{( N )}}} ]^{\rm{T}}}$ is formed, as is its covariance ${{\bf C}}_\mu ^{( N )} = {{{\boldsymbol \Lambda }}^{( N )}}{\boldsymbol{\ }}$; (4) the equation is solved to yield a localized average ${\langle {{\bf m}} \rangle ^{( N )}} = \ {{{\boldsymbol \Phi }}^{ - g}}{{{\boldsymbol \mu }}^{( N )}}$, where ${{{\boldsymbol \Phi }}^{ - g}}$ is either the minimum length or Backus–Gilbert generalized inverse of ${{\boldsymbol \Phi }}$; (5) the resolution and covariance are computed as ${{{\bf R}}^{( N )}} = {{{\boldsymbol \Phi }}^{ - g}}{\boldsymbol{\ }}{{{\boldsymbol \Phi }}^{( N )}}$ and ${{\bf C}}_m^{( N )} = {{{\boldsymbol \Phi }}^{ - g}}{\boldsymbol{\ }}{{\bf C}}_\mu ^{( N )}{( {{{{\boldsymbol \Phi }}^{ - g}}} )^{\rm{T}}}$; (6) the spread ${K^{( N )}}$ of resolution and size ${J^{( N )}}\ $of covariance are computed using either the Dirichlet or Backus–Gilbert measures and (7) the process is repeated for $1 \le N \le M$ to build up the trade-off curve $K( J )$. We show that, in the Dirichlet case, ${K^{( N )}} = \ M - N$ and ${J^{( N )}} = \ {\rm{tr}}( {{{{\boldsymbol \Lambda }}^{( N )}}} )$. We also consider the case where the model parameters correspond to spline coefficients and a sequence ${y_i}( {{{\bf m}},{x_i}} )$ derived from these coefficients possesses natural ordering. Layered models are an example of such a parametrization. We construct the trade-off curve for ${{\bf y}}$ by converting each member of the ensemble from ${{\bf m}}$ to ${{\bf y}}$ and applying the above procedure to them. We demonstrate the method by applying it to several simple examples.

A Model of Extracting Building from High Resolution Remote Sensing Image Based on Bayesian Convolutional Neural Networks

When extract building from high resolution remote sensing image with meter/sub-meter accuracy, the shade of trees and interference of roads are the main factors of reducing the extraction accuracy. Proposed a Bayesian Convolutional Neural Networks(BCNET) model base on standard fully convolutional networks(FCN) to solve these problems. First take building with no shade or artificial removal of shade as Sample-A, woodland as Sample-B, road as Sample-C. Set up 3 sample libraries. Learn these sample libraries respectively, get their own set of feature vector; Mixture Gauss model these feature vector set, evaluate the conditional probability density function of mixture of noise object and roofs; Improve the standard FCN from the 2 aspect:(1) Introduce atrous convolution. (2) Take conditional probability density function as the activation function of the last convolution. Carry out experiment using unmanned aerial vehicle(UVA) image, the results show that BCNET model can effectively eliminate the influence of trees and roads, the building extraction accuracy can reach 97%.

Effects of Compound K-Distributed Sea Clutter on Angle Measurement of Wideband Monopulse Radar

The effects of compound K-distributed sea clutter on angle measurement of wideband monopulse radar are investigated in this paper. We apply the conditional probability density function (pdf) of monopulse ratio (MR) error to analyze these effects. Based on the angle measurement procedure of the wideband monopulse radar, this conditional pdf is first deduced in detail for the case of compound K-distributed sea clutter plus noise. Herein, the spatial correlation of the texture components for each channel clutter and the correlation of the texture components between the sum and difference channel clutters are considered, and two extreme situations for each of them are tackled. Referring to the measured sea clutter data, angle measurement performances in various K-distributed sea clutter plus noise circumstances are simulated, and the effects of compound K-distributed sea clutter on angle measurement are discussed.