• Journal of Infrared and Millimeter Waves
  • Vol. 40, Issue 1, 64 (2021)
Sai LI1、2、3, Qiu YIN4, Yong HU1、2、*, and Cai-Lan GONG1、2
Author Affiliations
  • 1Shanghai Institute of Technical Physics,Chinese Academy of Sciences,Shanghai 200083,China
  • 2CAS Key Laboratory of Infrared System Detection and Imaging Technology,Shanghai Institute of Technical Physics,Shanghai 200083,China
  • 3University of Chinese Academy of Sciences,Beijing 100049,China
  • 4Shanghai Meteorological Bureau,Shanghai 200030,China
  • show less
    DOI: 10.11972/j.issn.1001-9014.2021.01.011 Cite this Article
    Sai LI, Qiu YIN, Yong HU, Cai-Lan GONG. A push-sweep hyperspectral aerial image Mosaic method based on SPHP[J]. Journal of Infrared and Millimeter Waves, 2021, 40(1): 64 Copy Citation Text show less

    Abstract

    In order to solve the problem that the optical center of hyperspectral images does not coincide with each other, and the alignment of traditional global monotonic matrix may have error in the registration process, a splicing-method (H-SPHP) of hyperspectral images based on SPHP(shape-preserving half-projective) and considering spectral information is proposed. The main steps include as follows: 1) image correction using the visual vector method. 2) selection of reference band based on prior knowledge and PCA. 3) SPHP-based mesh optimization splicing method. 4) Weighted average fusion algorithm for fusion.5) splicing parameters applied to all bands to obtain the splicing hyperspectral data. By obtaining experimental images from Sanming, Fujian province and Nanchang,Jiangxi province, the research data splicing experimental results show that the proposed algorithm has strong robustness, eliminating parallax and geographic coordinates accuracy better than SIFT + single should transform algorithm. After splicing, the registration accuracy of bands is within one pixel, and spectral similarity of the overlapping region is above 90%.
    Li=[0 -c×j-N2+1 -f]T .

    View in Article

    xe=N+hcos latcos lon 

    View in Article

    ye=N+hcos latsin lon 

    View in Article

    ze=N(1-e)2+hsin lon 

    View in Article

    N=Ra/1-e2sin2(lat) 

    View in Article

    RexReyRez=xeyeze+λTgeTdgTcdTicAi

    View in Article

    RexReyRez=xeyeze+λwxwywz

    View in Article

    Rex2ae2+Rey2ae2+Rez2be2=1

    View in Article

    (xe+λwx)2ae2+(ye+λwy)2ae2+(ze+λwz)2be2=1 .

    View in Article

    wx2ae2+wy2ae2+wz2be2λ2+2xewxae2+yewyae2+zewzbe2λ+xe2ae2+ye2ae2+ze2be2-1=0 .

    View in Article

    A=wx2ae2+wy2ae2+wz2be2

    View in Article

    B=2xewxae2+yewyae2+zewzbe2

    View in Article

    C=xe2ae2+ye2ae2+ze2be2-1

    View in Article

    Aλ2+Bλ+C=0 .

    View in Article

    E=M×SD .

    View in Article

    H0=h1h2h3h4h7h5h8h61 .

    View in Article

    x0y0=cosθ-sinθsinθcosθμ0v0

    View in Article

    θ=atan2(-h8̑-h7̑) .

    View in Article

    x0'y0'1~h1h2h3h4-kh50h61μ0v01

    View in Article

    h1h2h4h5=h1̑h2̑h4̑h5̑cosθ-sinθsinθcosθ

    View in Article

    k=h8̑2+h7̑2

    View in Article

    h3,h6=(h3̑,h6̑) .

    View in Article

    x0'=Hxμ0,v0=h1μ0+h2v0+h31-ku0

    View in Article

    y0'=Hyμ0,v0=h4μ0+h5v0+h61-ku0 .

    View in Article

    x0'=Hxμ0,v0=h2v01-ku0+h1μ0+h31-ku0

    View in Article

    y0'=Hyμ0,v0=h5v01-ku0+h4μ0+h61-ku0 .

    View in Article

    S0μ0,v0=11-kμ01(h5h2h2h5μ0v0+h1-h5μ01+h3h4+h2μ01+h6) .

    View in Article

    W0μ0,v0=H0μ0,v0T0μ0,v0S0μ0,v0    if(μ0,v0)ϵRHif(μ0,v0)ϵRTif(μ0,v0)ϵRS

    View in Article

    S0μ0,v0=Sxμ0,v0Syμ0,v0=α-ββαμ0v0+txty

    View in Article

    T0μ0,v0=Txμ0,v0Tyμ0,v0=fx(μ0)fy(μ0)v0+gx(μ0)gy(μ0)

    View in Article

    W0μ0,v0=Wxμ0,v0Wyμ0,v0=Fxμ0v0+Gx(μ0)Fyμ0v0+Gy(μ0)

    View in Article

    Fxμ0=11-kμ0Fxμ0-β    if μ0μ01if μ01μ0<μ02if μ02μ0

    View in Article

    Gxμ0=h1μ0+h31-kμ0Gxμ0αμ0+tx    if μ0μ01if μ01μ0<μ02if μ02μ0 .

    View in Article

    fxμ01=h21-kμ01

    View in Article

    fx'μ01=kh2(1-kμ01)2

    View in Article

    fxμ02=-β

    View in Article

    fx'μ02=0 .

    View in Article

    Ix,y=I0x,yρI0x,y+(1-ρ)I1x,yS0μ0,v0    (x,y)ϵI0(x,y)ϵI0(x,y)ϵI1I1

    View in Article

    T= 0.9556-0.0637 -290.5153-0.0251 0.9944 27.6141001 .

    View in Article

    DSAMx=cos-1mTx(mTm)12(xTx)12

    View in Article

    cosθ=i=1nxiyii=1nxi2i=1nyi2

    View in Article

    SFF=Cov(P,ρ)D(P)D(ρ)

    View in Article

    SC=i=1n(xik-xi¯)(xjk-xj¯)i=1n(xik-xi¯)2i=1n(xjk-xj¯)2

    View in Article

    Sai LI, Qiu YIN, Yong HU, Cai-Lan GONG. A push-sweep hyperspectral aerial image Mosaic method based on SPHP[J]. Journal of Infrared and Millimeter Waves, 2021, 40(1): 64
    Download Citation