Least-square　reverse　time　migration　using　plane-wave　coding　always　has　two　problems:　encoding　data　will　introduce　crosstalk　noiseand　the　excessive　number　of　plane　records　will　reduce　the　computational　efficiency.　In　this　paper,　the　Seislet　transform,　which　is　suitable　for　seismic　data,　is　combined　with　the　fractional　threshold　function　based　on　the　Riemann-Liouville　fractional　integration　theory,　which　is　then　introduced　into　the　least-square　reverse　time　migration　using　plane-wave　coding.　Numerical　tests　on　the　simple　layered　and　complex　model　show　that　seismic　data　have　good　sparsity　in　the　Seislet　domain　and　the　plane-wave　reverse　time　migration　based　on　Seislet　fractional　threshold　algorithm　can　effectively　suppress　the　crosstalk　noise　caused　by　multi-source　data.　Compared　with　the　traditional　method,　this　proposed　method　uses　less　number　of　plane-wave　gathers　to　obtain　the　same　imaging　effect,　and　can　improve　the　computational　efficiency.　©　2020,　Editorial　Office　of　Journal　of　China　University　of　Petroleum(Edition　of　Natural　Science).　All　right　reserved.