Patent Publication Number: US-11388339-B2

Title: Anti-shake method for panoramic video, and portable terminal

Description:
CROSS REFERENCE TO RELATED APPLICATIONS 
     The present application is a 35 U.S.C. § 371 National Phase conversion of International (PCT) Patent Application No. PCT/CN2019/078329, filed on Mar. 15, 2019, which claims priority of Chinese Patent Application No. 201810217087.6, filed on Mar. 16, 2018, the disclosure of which is incorporated by reference herein. The PCT International Patent Application was filed and published in Chinese. 
     FIELD OF THE INVENTION 
     The present invention relates generally to the field of video processing, and particularly relates to an anti-shake method for a panoramic video and a portable terminal. 
     BACKGROUND OF THE INVENTION 
     CMOS and CCD sensors are two types of image sensors that are currently commonly used. Both use photodiodes for photoelectric conversion to convert images into digital data. The main difference is the way in which digital data is transmitted. 
     The charge data of each pixel in each line of the CCD sensor will be transferred to the next pixel in turn, and be outputted from the bottom thereof, then be amplified and outputted by the amplifier at the edge of the sensor; while in the CMOS sensor, each pixel will be adjacent to one amplifier and an A/D conversion circuit, and data is outputted in a way similar to a memory circuit. 
     Generally, CMOS cameras use rolling shutter mode to obtain fisheye images through progressive exposure. 
     CMOS chips share the workload through many parallel A/Ds, but the entire array in the sensor must be converted one line at a time, which results in a small time delay between the readings of each line. Each individual line is usually able to start the exposure of the next frame when the readout of the previous frame is completed. Although the speed is very fast, the time delay between the readings of each line will be converted into the delay between the start of the exposure of each line, exposure of each line will not happen at the same time. The result is that each line in the frame will be exposed for the same time length, but the exposure will start at a different time point. Two frames are allowed to overlap exposure, and the final frame rate depends on the speed at which the scrolling readout process can be completed. 
     Such exposure mode causes a time difference between different lines of the same frame of image. If a handheld panoramic camera is used to shoot, when moving at high speeds, a jelly effect will occur due to the progressive exposure of the CMOS camera. The jelly effect is a problem that has not been well solved, especially for a VR panoramic video, and the anti-shake effect of video is poor. 
     SUMMARY OF THE INVENTION 
     The problem solved by the present invention is to provide an anti-shake method for a panoramic video, a computer-readable storage medium and a portable terminal, which aims to solve the jelly effect caused by the shaking of the picture due to the rolling shutter mode of the CMOS chip, and the poor anti-shake effect of the camera. 
     According to a first aspect, the present invention provides an anti-shake method for a panoramic video, comprises steps of: 
     obtaining, in real time, an output video frame, a fisheye image corresponding to the output video frame, a timestamp of the pixel in the video frame, and a corresponding camera gyroscope timestamp; 
     synchronizing the timestamp of the pixel in the video frame with the corresponding camera gyroscope timestamp, and calculating a rotation matrix of the camera movement in the camera gyroscope timestamp; 
     smoothing the camera movement and establishing a coordinate system of a smooth trajectory; 
     correcting the fisheye image distortion; and 
     rendering the fisheye image by means of forward rendering to generate a stable video. 
     According to a second aspect, the present invention provides an anti-shake method for a panoramic video, comprises steps of: 
     S 201 , obtaining, in real time, an output video frame, a fisheye image corresponding to the output video frame, a timestamp of the pixel in the video frame, and a corresponding camera gyroscope timestamp; 
     S 202 , synchronizing the timestamp in the video frame to the camera gyroscope acquisition system, and calculating a rotation matrix of the camera movement in the timestamp of the pixel in the current video frame as the approximate rotation matrix of the camera when collecting the current pixel; 
     S 203 , projecting the output video frame pixels to the corresponding grid points of the spherical model, and using the approximate rotation matrix to rotate the corresponding grid points of the spherical model to the grid points in the sensor coordinate system; using the fisheye image distortion correction method to establish the relationship between the grid points in the sensor coordinate system and the pixels in the fisheye image to obtain approximate reverse-mapping pixels; 
     S 204 , calculating the camera gyroscope timestamp of the approximate reverse-mapping pixels in the fisheye image, obtaining the accurate rotation matrix of the camera when collecting the current pixel, and using the accurate rotation matrix to rotate the corresponding grid points of spherical model again to obtain the second grid points in the sensor coordinate system; 
     S 205 , using the fisheye image distortion correction method to establish the relationship between the second grid points in the sensor coordinate system and the pixels in the fisheye image to obtain accurate pixels in the fisheye image, and using the mapping relationship between the output video frame pixels and the accurate pixels in the fisheye image, rendering the fisheye image in reverse rendering to generate a stable video. 
     According to a third aspect, the present invention provides a computer-readable medium that stores a computer program or computer programs including set of computer-executable instructions, which when being executed by a processor or processors, cause the processor or processors to perform steps of the above-mentioned anti-shake method for a panoramic video. 
     According to a fourth aspect, the present invention provides a portable terminal, comprising: 
     one or more processors; 
     a memory; and 
     one or more computer programs where the one or more computer programs including set of computer-executable instructions are stored in the memory and are configured to be executed by the one or more processors, and when executed by the one or more processors, cause the one or more processors to perform the above-mentioned anti-shake method for a panoramic video. 
     In the present invention, based on the timestamp of the pixel in the video frame, the fisheye image is pixel by pixel converted into the coordinate system of the smooth trajectory using the rotation matrix of the camera motion under the timestamp of the current video frame pixel in real time, and then rendered and output; thereby the distortion of the rolling shutter of the panoramic video sequence can be corrected. The invention can correct the image distortion caused by the CMOS rolling shutter and eliminate the jelly effect, and achieve better video image anti-shake effect; 
     Further, in the present invention, according to the timestamp of the video frame pixel, calculating the rotation matrix of the camera movement under the timestamp of the current video frame pixel as the approximate rotation matrix of the camera when collecting the current pixel, projecting the output video frame pixels to the corresponding grid points of the spherical model, and using the approximate rotation matrix to rotate the corresponding grid points of the spherical model to the grid points in the sensor coordinate system; using the fisheye image distortion correction method to establish the relationship between the grid points in the sensor coordinate system and the pixels in the fisheye image to obtain approximate reverse mapping pixels, calculating the camera gyroscope timestamp of the approximate reverse-mapping pixels in the fisheye image, obtaining the accurate rotation matrix of the camera when collecting the current pixel, and using the accurate rotation matrix to rotate the corresponding grid points of spherical model again to obtain the second grid points in the sensor coordinate system; using the fisheye image distortion correction method to establish the relationship between the second grid points in the sensor coordinate system and the pixels in the fisheye image to obtain accurate pixels in the fisheye image, and using the mapping relationship between the output video frame pixels and the accurate pixels in the fisheye image, rendering the fisheye image in reverse rendering to generate a stable video. Thereby, the distortion of the rolling shutter of the panoramic video sequence can be corrected. The invention can correct the image distortion caused by the CMOS rolling shutter and eliminate the jelly effect, thereby achieve a better video image anti-shake effect. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIG. 1  is a flowchart of an anti-shake method for a panoramic video in accordance with a first embodiment of the present invention; 
         FIG. 2  is a schematic diagram of the timestamp of the scan line sequence of the CMOS camera in accordance with the first embodiment of the present invention; 
         FIG. 3  is a mapping process diagram of a forward mapping rendering method in the anti-shake method for a panoramic video in accordance with the first embodiment of the present invention; 
         FIG. 4  is a flowchart of the anti-shake method for a panoramic video in accordance with a second embodiment of the present invention; 
         FIG. 5  is a mapping process diagram of a backward mapping rendering method in the anti-shake method for a panoramic video in accordance with the second embodiment of the present invention; and 
         FIG. 6  is a schematic structural diagram of a portable terminal in accordance with a fourth embodiment of the present invention. 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     The foregoing objects, technical solutions and advantages of the invention will be much clearer from the following detail description taken with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. 
     In order to explain the technical solutions of the present invention, the following will be described by specific embodiments. 
     First Embodiment 
     Referring to  FIG. 1 , an anti-shake method for a panoramic video provided in the first embodiment of the present invention comprising the following steps: 
     S 101 , obtaining, in real time, an output video frame, a fisheye image corresponding to the output video frame, a timestamp of the pixel in the video frame, and a corresponding camera gyroscope timestamp; 
     in the first embodiment of the present invention, S 101  may specifically comprise the following steps: 
     S 1011 , obtaining the output video frame and the fisheye image corresponding to the output video frame; 
     S 1012 , obtaining the timestamp t frame   k     p(x,y)    of the pixel p(x, y) in the k video frame in the time axis of a video frame acquisition system, where k is a natural number; and 
     S 1013 , obtaining the gyroscope timestamp t gyro   k     p(x,y)    in the time axis of the corresponding camera gyroscope acquisition system. 
     S 102 , synchronizing the timestamp of the pixel in the video frame with the corresponding camera gyroscope timestamp, and calculating a rotation matrix of the camera movement in the camera gyroscope timestamp; 
     referring to  FIG. 2 , in the first embodiment of the present invention, S 102  may specifically comprise the following steps: 
     S 1021 , calculating the timestamp t frame   k     p(x,y)    of the pixel p(x, y) in the k video frame by formula (1): 
     
       
         
           
             
               
                 
                   
                     
                       t 
                       frame 
                       
                         k 
                         
                           p 
                           ⁡ 
                           
                             ( 
                             
                               x 
                               , 
                               y 
                             
                             ) 
                           
                         
                       
                     
                     = 
                     
                       
                         t 
                         frame 
                         k 
                       
                       + 
                       
                         
                           
                             Δ 
                             ⁢ 
                             t 
                           
                           H 
                         
                         ⁢ 
                         y 
                       
                       - 
                       
                         
                           1 
                           2 
                         
                         ⁢ 
                         Δ 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                     
                   
                   ; 
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
     in formula (1), t frame   k     p(x,y)    is the timestamp of the pixel at the beginning of the k video frame, Δt is a sampling time interval of the progressive scan of the video frame, and H is the number of lines of the image; 
     S 1022 , synchronizing the timestamp of the pixel in the k video frame with the camera gyroscope timestamp, a conversion relationship between the two as formula (2):
 
 t   gyro   k     p(x,y)     =ΔT   0   +t   frame   k     p(x,y)     (2)
         in formula (2), ΔT 0 =t gyro   k −t frame   k , t gyro   k  is the camera gyroscope timestamp at the beginning of the k video frame;       

     S 1023 , calculating a rotation matrix of the camera movement in the camera gyroscope timestamp at time t gyro   k     p(x,y)    by formula (3):
 
 R   w2s ( t   gyro   k     p(x,y)   )= R   w2g ( t   gyro   k     p(x,y)   )· R   g2s ( t   gyro   k     p(x,y)   )  (3)
 
     in formula (3), R w2g (t gyro   k     p(x,y)   ) is a rotation matrix measured by a gyroscope sensor at time t gyro   k     p(x,y)   , R g2s (t gyro   k     p(x,y)   ) is a rotation matrix from the calibrated gyroscope coordinate system to the camera, and R w2s (t gyro   k     p(x,y)   ) is a 3*3 matrix. 
     S 103 , smoothing the camera movement and establishing a coordinate system of a smooth trajectory; 
     in the first embodiment of the present invention, S 103  may specifically comprise the following steps: 
     S 1031 , where P s  representing the 3D coordinates in the camera coordinate system, let t=t gyro   k     p(x,y)   , converting the 3D coordinates P t   w  in the world coordinate system to the 3D coordinates P t   s  in the camera coordinate system at time t by a formula:
 
 P   t   s   =R   w2s ( t )· P   t   w +{right arrow over ( t   0 )}  (4)
 
     in the formula (4), R w2s (t) is a matrix of the camera coordinate system relative to the world coordinate system at time t, and {right arrow over (t 0 )} is a translation amount of the two frames before and after of the video frame of the camera; 
     S 1032 , smoothing the camera movement, where setting {right arrow over (t 0 )}=0 and then obtaining 
     
       
         
           
             
               
                 P 
                 t 
                 w 
               
               = 
               
                 
                   P 
                   t 
                   s 
                 
                 
                   
                     R 
                     
                       w 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       2 
                       ⁢ 
                       s 
                     
                   
                   ⁡ 
                   
                     ( 
                     t 
                     ) 
                   
                 
               
             
             ; 
           
         
       
     
     S 1033 , establishing the coordinate system P t ′ of the smooth trajectory, performing a 3D grid rendering processing on the coordinate system P t ′ of the camera&#39;s smooth trajectory by OpenG1 Rendering formula:  P t   g1 =K·R mvp ·P t ′ ; where K is a perspective matrix, R mvp  is a rotation matrix of the view direction manually controlled, and P t   g1  is the 3D coordinates rendered by OpenG1 in the coordinate system of the smooth trajectory; and 
     S 1034 , in the coordinate system of the smooth trajectory, obtaining the minimum square difference of any pixel in P t ′ in the two frames before and after by a formula: min(P t+1 ′−P t ′) 2 +(P t ′−P t   w ) 2 , where P t ′ and P t+1 ′ are the coordinate systems of the smooth trajectory of the two frames before and after, respectively; setting P t ′=P t+1 ′=P t   w , and obtaining the rendering formula as formula (5):
 
 P   t   g1   =K·R   mvp   ·R   x2z(t)   −1   ·P   t   s   (5)
 
calculating the 3D coordinates of the camera rendered to the coordinate system of the smooth trajectory by formula (5).
 
     S 104 , correcting the fisheye image distortion. 
     In the first embodiment of the present invention, S 104  may specifically comprise: 
     where let 
                 P   s     =     P     t   gyro     k     p   ⁡     (     x   ,   y     )           s       ,         
obtaining a panoramic expanded image of P s  in the camera coordinate system:
 
     
       
         
           
             
               
                 
                   
                     P 
                     s 
                   
                   = 
                   
                     { 
                     
                       
                         
                           
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 u 
                                 ) 
                               
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 v 
                                 ) 
                               
                             
                           
                         
                       
                       
                         
                           
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 u 
                                 ) 
                               
                             
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 v 
                                 ) 
                               
                             
                           
                         
                       
                       
                         
                           
                             sin 
                             ⁡ 
                             
                               ( 
                               u 
                               ) 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
     in formula (6), u and v are the abscissa and ordinate of the panoramic expanded image in the camera coordinate system, respectively; 
     a fisheye distortion model ρ=f(θ)=θ(1+k 0 θ 2 +k 1 θ 4 +k 2 θ 5 +k 3 θ 8 ), where θ is an incident angle of the light, k 0 , k 1 , k 2 , and k 3  are the camera calibration coefficients, obtaining the corrected coordinates (x, y) of the fisheye image: 
     
       
         
           
             
               
                 
                   { 
                   
                     
                       
                         
                           x 
                           = 
                           
                             
                               
                                 f 
                                 x 
                               
                               · 
                               ρ 
                               · 
                               
                                 
                                   P 
                                   x 
                                   s 
                                 
                                 r 
                               
                             
                             + 
                             
                               x 
                               c 
                             
                           
                         
                       
                     
                     
                       
                         
                           y 
                           = 
                           
                             
                               
                                 f 
                                 y 
                               
                               · 
                               ρ 
                               · 
                               
                                 
                                   P 
                                   y 
                                   s 
                                 
                                 r 
                               
                             
                             + 
                             
                               y 
                               c 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
           
         
       
     
     in formula (7), r=√{square root over (P x   s2 +P y   s2 )}, and θ=atan (r,P z   s ); the point (x c , y c ) is a projection center of the sensor, x and y are the 2D coordinates of the projection position of the camera sensor respectively. 
     S 105 , rendering the fisheye image by means of forward rendering to generate a stable video. 
     In the first embodiment of the present invention, S 105  may specifically comprise the following steps of: 
     establishing a correspondence between the fisheye image and the grid points of the spherical model in the camera coordinate system by means of the fisheye image distortion correction method; 
     generating the grid points of the spherical model in the coordinate system of the smooth trajectory by rendering; and 
     performing panoramic expansion of the grid points of the spherical model in the coordinate system of the smooth trajectory to generate a stable video. Where, specifically, color information can be input by interpolation, and projection can generate perspective images or asteroid images. 
     Referring to  FIG. 3 , in the first embodiment of the present invention: 
     A 1  corresponds to A 2 , where A 1  is the coordinate in the fisheye image, A 2  corresponds to A 1 , is the coordinate of the grid point of the spherical model in the camera coordinate system; A 3  corresponds to A 2 , is the coordinate of the grid point of the spherical model in the coordinate system of the smooth trajectory generated by rendering; and A 4  is the coordinate point generated by the panoramic expansion projection of A 3 ; 
     Similarly, B 1  corresponds to B 2 , where B 1  is the coordinate in the fisheye image, B 2  corresponds to B 1 , is the coordinate of the grid point of the spherical model in the camera coordinate system; B 3  corresponds to B 2 , is the coordinate of the grid point of the spherical model in the coordinate system of the smooth trajectory generated by rendering; and B 4  is the coordinate point generated by the panoramic expansion projection of B 3 . 
     In the first embodiment of the present invention, based on the timestamp of the pixel in the video frame, the fisheye image is pixel by pixel converted into the coordinate system of the smooth trajectory using the rotation matrix of the camera motion under the timestamp of the current video frame pixel in real time, and then rendered and output; thereby the distortion of the rolling shutter of the panoramic video sequence can be corrected. The invention can correct the image distortion caused by the CMOS rolling shutter and eliminate the jelly effect, and achieve better video image anti-shake effect. 
     Second Embodiment 
     Referring to  FIG. 4 , an anti-shake method for a panoramic video provided in is the second embodiment of the present invention comprising the following steps: 
     S 201 , obtaining, in real time, an output video frame, a fisheye image corresponding to the output video frame, a timestamp of the pixel in the video frame, and a corresponding camera gyroscope timestamp; 
     in the second embodiment of the present invention, S 201  may specifically comprise steps of: 
     S 2011 , obtaining the output video frame and the fisheye image corresponding to the output video frame; 
     S 2012 , obtaining the timestamp t frame   k     p(x,y)    of the pixel p(x, y) in the k video frame in the time axis of a video frame acquisition system, where k is a natural number; and 
     S 2013 , obtaining the gyroscope timestamp t gyro   k     p(x,y)    in the time axis of the corresponding camera gyroscope acquisition system. 
     S 202 , synchronizing the timestamp of the pixel in the video frame to the camera gyroscope acquisition system, and calculating the rotation matrix of the camera movement in the timestamp of the pixel in the current video frame as the approximate rotation matrix of the camera when collecting the current pixel; 
     in the second embodiment of the present invention, S 202  may specifically include the following steps: 
     S 2021 , synchronizing the timestamp of the pixel in the video frame to the camera gyroscope acquisition system, the conversion relationship between the two as shown in formula (8):
 
 t   gyro   k   =ΔT   0   +t   frame   k   (8)
 
     in formula (8), ΔT 0 =t gyro   k −t frame   k , t gyro   k  is the camera gyroscope timestamp at the beginning of the k video frame, and t frame   k  is the timestamp of the pixel at the beginning of the k video frame; 
     S 2022 , calculating a rotation matrix of the camera movement under the timestamp of the pixel in the current video frame at time t gyro   k  by formula (9):
 
 R   w2s ( t   gyro   k )= R   w2g ( t   gyro   k )· R   g2s ( t   gyro   k )  (9)
 
     in formula (9), R w2g (t gyro   k ) is the rotation matrix measured by the gyroscope sensor at time t gyro   k , R w2g (t gyro   k ) is a 3*3 matrix; and R g2s (t gyro   k ) is the rotation matrix from the calibrated gyroscope coordinate system to the camera. 
     S 203 , projecting the output video frame pixels to the corresponding grid points of the spherical model, and using the approximate rotation matrix to rotate the corresponding grid points of the spherical model to the grid points in the sensor coordinate system; using the fisheye image distortion correction method to establish the relationship between the grid points in the sensor coordinate system and the pixels in the fisheye image to obtain approximate reverse-mapping pixels; 
     in the second embodiment of the present invention, S 203  may specifically be: 
     where let 
                 P   s     =     P     t   gyro     k     p   ⁡     (     x   ,   y     )           s       ,         
obtaining a panoramic expanded image of P s  in the camera coordinate system:
 
                 P   s     =     {             cos   ⁡     (   u   )       ⁢     cos   ⁡     (   v   )                     cos   ⁡     (   u   )       ⁢     sin   ⁡     (   v   )                   sin   ⁡     (   u   )                 ,         
where u and v are the abscissa and ordinate of the panoramic expanded image in the camera coordinate system, respectively;
 
     a fisheye distortion model  92 =f(θ)=θ(1+k 0 θ 2 +k 1 θ 4 +k 2 θ B +k 3 θ B ), where θ is an incident angle of the light, k 0 , k 1 , k 2 , and k 3  are the camera calibration coefficients, obtaining the corrected coordinates (x, y) of the fisheye image: 
     
       
         
           
               
             
               { 
               
                 
                   
                     
                       x 
                       = 
                       
                         
                           
                             f 
                             x 
                           
                           · 
                           ρ 
                           · 
                           
                             
                               P 
                               x 
                               s 
                             
                             r 
                           
                         
                         + 
                         
                           x 
                           c 
                         
                       
                     
                   
                 
                 
                   
                     
                       y 
                       = 
                       
                         
                           
                             f 
                             y 
                           
                           · 
                           ρ 
                           · 
                           
                             
                               P 
                               y 
                               s 
                             
                             r 
                           
                         
                         + 
                         
                           y 
                           c 
                         
                       
                     
                   
                 
               
             
           
         
       
     
     where r=√{square root over (P x   s2 +P y   s2 )} and θ=atan (r,P z   s ); the point (x c ,y c ) is a projection center of the sensor, x and y are the 2D coordinates of the projection position of the camera sensor respectively. 
     S 204 , calculating the camera gyroscope timestamp of the approximate reverse-mapping pixels in the fisheye image, obtaining the accurate rotation matrix of the camera when collecting the current pixel, and using the accurate is rotation matrix to rotate the corresponding grid points of spherical model again to obtain the second grid points in the sensor coordinate system; 
     in the second embodiment of the present invention, the step of calculating the camera gyroscope timestamp of the approximate reverse-mapping pixels in the fisheye image, obtaining the accurate rotation matrix of the camera motion when collecting the current pixel, specifically is: 
     using the approximate reverse-mapping pixels in the fisheye image to calculate the camera gyroscope timestamp of the pixel by the formula: 
     
       
         
           
             
               
                 
                   
                     t 
                     gyro 
                     
                       k 
                       
                         p 
                         ⁡ 
                         
                           ( 
                           
                             x 
                             , 
                             y 
                           
                           ) 
                         
                       
                     
                   
                   = 
                   
                     
                       Δ 
                       ⁢ 
                       
                         T 
                         0 
                       
                     
                     + 
                     
                       t 
                       frame 
                       
                         k 
                         
                           p 
                           ⁡ 
                           
                             ( 
                             
                               x 
                               , 
                               y 
                             
                             ) 
                           
                         
                       
                     
                     + 
                     
                       
                         
                           Δ 
                           ⁢ 
                           t 
                         
                         H 
                       
                       ⁢ 
                       y 
                     
                     - 
                     
                       
                         1 
                         2 
                       
                       ⁢ 
                       Δ 
                       ⁢ 
                       t 
                     
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
     in formula (10), Δt is a sampling time interval of the progressive scan of the video frame, and H is the number of lines of the image; 
     further calculating the accurate rotation matrix R w2s (t gyro   k     p(x,y)   ) of the camera motion under the camera gyroscope timestamp, where the calculation formula is:
 
 R   w2s ( t   gyro   k     p(x,y)   )= R   w2g ( t   gyro   k     p(x,y)   )· R   g2s ( t   gyro   k     p(x,y)   ).
 
     S 205 , using the fisheye image distortion correction method to establish the relationship between the second grid points in the sensor coordinate system and the pixels in the fisheye image to obtain accurate pixels in the fisheye image, and using the mapping relationship between the output video frame pixels and the accurate pixels in the fisheye image, rendering the fisheye image in reverse rendering to generate a stable video. 
     Referring to  FIG. 5 , in the second embodiment of the present invention: 
     projecting the output video frame pixel A 4  to the corresponding grid point A 3  of the spherical model, using the approximate rotation matrix to rotate A 3  to obtain A 2  in the sensor coordinate system, using the fisheye image distortion correction method to project A 2  onto the fisheye image to obtain the pixel A 1 , using the accurate rotation matrix to rotate A 1  to obtain A′ 2  in the sensor coordinate system, and then using the fisheye image distortion correction method to project A′ 2  to the fisheye image to obtain the accurate pixel A′ 1 ; using the mapping relationship between the output video frame pixel A 4  and the accurate pixel A′ 1  in the fisheye image to render the fisheye image in reverse rendering to generate a stable video; 
     Similarly, projecting the output video frame pixel B 4  to the corresponding grid point B 3  of the spherical model, using the approximate rotation matrix to rotate B 3  to obtain B 2  in the sensor coordinate system, using the fisheye image distortion correction method to project B 2  onto the fisheye image to obtain pixel B 1 , using the accurate rotation matrix to rotate B 1  to obtain B′ 2  in the sensor coordinate system, and then using the fisheye image distortion correction method to project B′ 2  to the fisheye image to obtain the accurate pixel B′ 1 ; using the mapping between the output video frame pixel B 4  and the accurate pixel B′ 1  in the fisheye image to render the fisheye image in reverse rendering to generate a stable video; 
     It should also be noted that the perspective view generated by the projection or the color information of the asteroid image is input by interpolation, and finally a stable video can be generated. 
     In the second embodiment of the present invention, according to the timestamp of the video frame pixel, calculating the rotation matrix of the camera movement under the timestamp of the current video frame pixel as the approximate rotation matrix of the camera when collecting the current pixel, projecting the output video frame pixels to the corresponding grid points of the spherical model, and using the approximate rotation matrix to rotate the corresponding grid points of the is spherical model to the grid points in the sensor coordinate system; using the fisheye image distortion correction method to establish the relationship between the grid points in the sensor coordinate system and the pixels in the fisheye image to obtain approximate reverse-mapping pixels, calculating the camera gyroscope timestamp of the approximate reverse-mapping pixels in the fisheye image, obtaining the accurate rotation matrix of the camera when collecting the current pixel, and using the accurate rotation matrix to rotate the corresponding grid points of spherical model again to obtain the second grid points in the sensor coordinate system; using the fisheye image distortion correction method to establish the relationship between the second grid points in the sensor coordinate system and the pixels in the fisheye image to obtain accurate pixels in the fisheye image, and using the mapping relationship between the output video frame pixels and the accurate pixels in the fisheye image, rendering the fisheye image in reverse rendering to generate a stable video. Thereby, the distortion of the rolling shutter of the panoramic video sequence can be corrected. 
     The invention can correct the image distortion caused by the CMOS rolling shutter and eliminate the jelly effect, thereby achieve a better video image anti-shake effect. 
     Third Embodiment 
     The third embodiment of the present invention provides a computer-readable storage medium that stores a computer program or computer programs including set of computer-executable instructions, which when being executed by a processor or processors, cause the processor or processors to perform steps of the anti-shake method for a panoramic video provided in the first or second embodiment of the present invention. 
     Fourth Embodiment 
       FIG. 6  shows a structural block diagram of a portable terminal provided in the fourth embodiment of the present invention. A portable terminal  100  comprises: 
     one or more processors  101 , a memory  102 , and one or more computer programs, where the one or more processors  101  and the memory  102  are connected by a bus, the one or more computer programs including set of computer-executable instructions are stored in the memory  102 , and are configured to be executed by the one or more processors  101 . The one or more processors  101  execute the one or more computer programs to perform steps of the anti-shake method for a panoramic video provided in the first or second embodiment of the present invention. 
     A person of ordinary skill in the art may understand that all or part of the steps in the method of the above-mentioned embodiments can be implemented by a program or programs instructing relevant hardware. The program or programs can be stored in a computer-readable storage medium, and the storage media can be, such as ROM/RAM, magnetic disk, optical disk, etc. 
     The above descriptions are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modification, equivalent replacements and improvement made within the spirit and principle of the present invention shall be included in the protection of the present invention.