# Patent

## Title
Method of precisely correcting geometrically distorted satellite images and computer-readable storage medium for the method

## Abstract

Disclosed herein is a method of precisely correcting geometrically distorted satellite images. In order to correct a geometrically distorted satellite image, the position and attitude of a satellite must be accurately recognized. In this case, a correlation between the position and attitude values of the satellite exists, so accurate position and attitude values cannot be obtained. Therefore, the present invention accurately obtains position and attitude information of the satellite by separating variables having a high correlation to minimize the correlation between the position and attitude of the satellite and then calculating corresponding variables, and precisely corrects a geometrically distorted satellite image using the position and attitude information of the satellite.

## Background
PRIORITY CLAIM
This application claims priority from Korean patent application No. 2002-81171, filed 18 Dec. 2002, which is herein incorporated by reference.
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates generally to a method of precisely correcting geometrically distorted satellite images, and more particularly to a method of precisely correcting geometrically distorted satellite images, which measures position and attitude values of a satellite by separating variables having a high correlation therebetween, and which precisely corrects geometrically distorted satellite images using the position and attitude values of the satellite. In addition, the present invention relates to a computer-readable storage medium in which a program for executing the satellite image correcting method with a computer is stored.
2. Description of the Prior Art
Generally, remote sensing data obtained by an artificial satellite or the like may include many errors and distortions according to equipment allowances and atmospheric conditions at the time of observation, the trajectory and photographing position and attitude of a photographing satellite, the curvature and rotation of the earth, etc. Technologies for correcting geometrically distorted satellite images are very important in application fields, such as precise map production, multi-temporal image applications, and remote sensing applications.
Satellite image data including errors and distortions are designated as geometrically distorted satellite images. Such a geometrical distortion can be corrected by geometrically modeling the position of a satellite as shown by the dashed line FIG. 1, the photographing angle of a sensor, the curvature and rotating speed of the earth, etc. Such a geometric distortion correcting method is designated as geometric correction. However, even after the geometric correction is carried out, a remaining error (a difference between an actual ground control point and a corrected ground control point) still exists. The reason for existence of the remaining error is that the determination of the position and attitude of the satellite when the satellite image is photographed is inaccurate.
Therefore, ground control points are extracted from the satellite image, and the position and attitude of the satellite at the time of photographing are obtained using the extracted ground control points, such that the remaining error is eliminated. This process is called precision geometric correction. In this case, the ground control points are coordinates obtained by extracting points, having the same shape in both the satellite image and a reference map, from the satellite image so as to allow the satellite image to correspond to the reference map.
As described above, in order to precisely correct the satellite image, more accurate position and attitude values of the satellite must be obtained. The position and attitude values of the satellite have a high correlation therebetween. This means that the position and attitude values are related to each other, so the variation of one of them influences the other.
Due to the correlation between the position and attitude values of the satellite, it is impossible to accurately estimate position and attitude values using a typical least square method. In the prior art, the correlation therebetween is completely ignored, or a part of variables are set to constants and a plurality of ground control points are used, such that the satellite image is precisely corrected. Therefore, if plural accurate ground control points are used, the remaining error can be reduced, but accurate position and attitude values of the satellite cannot be obtained. In this case, the position and attitude values of the satellite, which are previously obtained, cannot be utilized for precision correction of another satellite image adjacent to the satellite image. Further, precision correction cannot be carried out for a satellite image obtained by photographing an area from which ground control points cannot be extracted.
FIG. 1 illustrates a geometric relationship between image coordinates (x, y) of a satellite image photographed by a camera and coordinates (X, Y, Z) of a ground control point “A”. The camera photographs a three-dimensional photographing area as a two-dimensional image, so the coordinates (x, y) of the photographed satellite image and the ground control point (X, Y, Z) have a certain relation. However, since the satellite image is geometrically distorted, a remaining error ε exists between actual coordinates (x, y) corresponding to an arbitrary ground control point and calculated coordinates (x′, y′) from the geometric modeling. Hereinafter, a method of correcting the remaining error of the satellite image is described in detail.
First, there is set a model that represents a geometric relationship between a reference coordinate system of a camera sensor of a satellite at the time of photographing a satellite image and a reference ground coordinate system of a photographed area (photographed ground surface), or between image coordinates (x, y) and a ground control point A (X, Y, Z). In this case, a geometric model of a linear pushbroom sensor is used as an example. The geometric model of the linear pushbroom sensor is expressed by the following collinearity Equations [1] and [2],




  
    
      
        x
        =
        
          
            
              -
              f
            
            ⁢
            
              
            
            ⁢
            
              
                
                  
                    r
                    11
                  
                  ⁢
                  
                    (
                    
                      X
                      -
                      
                        X
                        s
                      
                    
                    )
                  
                
                +
                
                  
                    r
                    21
                  
                  ⁡
                  
                    (
                    
                      Y
                      -
                      
                        Y
                        s
                      
                    
                    )
                  
                
                +
                
                  
                    r
                    31
                  
                  ⁡
                  
                    (
                    
                      Z
                      -
                      
                        Z
                        s
                      
                    
                    )
                  
                
              
              
                
                  
                    r
                    13
                  
                  ⁡
                  
                    (
                    
                      X
                      -
                      
                        X
                        s
                      
                    
                    )
                  
                
                +
                
                  
                    r
                    23
                  
                  ⁡
                  
                    (
                    
                      Y
                      -
                      
                        Y
                        s
                      
                    
                    )
                  
                
                +
                
                  
                    r
                    33
                  
                  ⁡
                  
                    (
                    
                      Z
                      -
                      
                        Z
                        s
                      
                    
                    )
                  
                
              
            
          
          =
          0
        
      
    
    
      
        [
        1
        ]
      
    
  
  
    
      
        y
        =
        
          
            -
            f
          
          ⁢
          
            
          
          ⁢
          
            
              
                
                  r
                  12
                
                ⁢
                
                  (
                  
                    X
                    -
                    
                      X
                      s
                    
                  
                  )
                
              
              +
              
                
                  r
                  22
                
                ⁡
                
                  (
                  
                    Y
                    -
                    
                      Y
                      s
                    
                  
                  )
                
              
              +
              
                
                  r
                  32
                
                ⁡
                
                  (
                  
                    Z
                    -
                    
                      Z
                      s
                    
                  
                  )
                
              
            
            
              
                
                  r
                  13
                
                ⁡
                
                  (
                  
                    X
                    -
                    
                      X
                      s
                    
                  
                  )
                
              
              +
              
                
                  r
                  23
                
                ⁡
                
                  (
                  
                    Y
                    -
                    
                      Y
                      s
                    
                  
                  )
                
              
              +
              
                
                  r
                  33
                
                ⁡
                
                  (
                  
                    Z
                    -
                    
                      Z
                      s
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        [
        2
        ]
      
    
  




where f is a focal length of the camera, and (Xs, Ys, Zs) are focal positions of the linear pushbroom sensor and (X, Y, Z) is coordinate of a ground control point. In this case, the focal positions of the linear pushbroom sensor correspond to desired position values of the satellite. Further, r11 to r33 are elements of a matrix R for rotating a coordinate system of the linear pushbroom sensor so as to allow the linear pushbroom sensor coordinate system to correspond to a ground coordinate system, and the matrix R is expressed by the following Equation [3],




  
    
      
        R
        =
        
          (
          
            
              
                
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                  ⁢
                  
                    
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  κ
                
              
              
                
                  
                    -
                    cos
                  
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                  ⁢
                  
                    
                  
                  ⁢
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  κ
                
              
              
                
                  sin
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                
              
            
            
              
                
                  
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    φcos
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                  +
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                
              
              
                
                  
                    
                      -
                      sin
                    
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    φsin
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                  +
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                
              
              
                
                  
                    -
                    sin
                  
                  ⁢
                  
                    
                  
                  ⁢
                  ω
                  ⁢
                  
                    
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                
              
            
            
              
                
                  
                    
                      -
                      cos
                    
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    φ
                    ⁢
                    
                      
                    
                    ⁢
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                  +
                  
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                
              
              
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    φ
                    ⁢
                    
                      
                    
                    ⁢
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                  +
                  
                    sin
                    ⁢
                    
                      
                    
                    ⁢
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    κ
                  
                
              
              
                
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  ω
                  ⁢
                  
                    
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  φ
                
              
            
          
          )
        
      
    
    
      
        [
        3
        ]
      
    
  




[3] where κ, φ and ω are rotating angles at which the linear pushbroom sensor coordinate system is rotated around Z, Y and X axes of the satellite coordinate system, respectively, so as to allow the linear pushbroom sensor coordinate system to correspond to the ground coordinate system. Unlike a perspective projection sensor, the linear pushbroom sensor is designed such that its focal positions are determined depending upon lines or parts of a satellite image, and its attitude can vary depending upon lines or parts thereof. Therefore, the focal position values (Xs, Ys, Zs) of the linear pushbroom sensor in Equations [1] and [2], and the rotating angles κ, φ and ω thereof can be represented by linear or non-linear polynomial expressions for an image coordinate value x according to a photographing manner or a scanning manner of the linear pushbroom sensor. In this case, κ, φ and ω are desired attitude values of the satellite to be obtained.
As described above, in the collinearity Equations expressed in Equations [1] and [2], the satellite photographs an image while moving along a direction of its movement, so the focal positions of the linear pushbroom sensor, that is, position values (Xs, Ys, Zs) of the satellite, and the elements r11 to r33 of the rotation matrix R are expressed as functions of time, not as constants. In the prior art, the position values (Xs, Ys, Zs) and the attitude values (κ, φ and ω) of the satellite are modeled as quadratic equations of time as indicated in the following Equation [4].




  
    
      
        
          
            
              
                X
                z
              
              =
              
                
                  X
                  0
                
                +
                
                  
                    a
                    1
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    1
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
          
            
              
                Y
                z
              
              =
              
                
                  Y
                  0
                
                +
                
                  
                    a
                    2
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    2
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
        
        
          
            
              
                Z
                z
              
              =
              
                
                  Z
                  0
                
                +
                
                  
                    a
                    3
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    3
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
          
            
              
                κ
                z
              
              =
              
                
                  κ
                  0
                
                +
                
                  
                    a
                    4
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    4
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
        
        
          
            
              
                φ
                z
              
              =
              
                
                  φ
                  0
                
                +
                
                  
                    a
                    5
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    5
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
          
            
              
                ω
                s
              
              =
              
                
                  ω
                  0
                
                +
                
                  
                    a
                    6
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    6
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
        
      
    
    
      
        [
        4
        ]
      
    
  




In order to obtain the position and attitude of the satellite in the prior art, φ and ω are assumed to be constants, and quadratic equations of the remaining four variables are applied to the collinearity equations, a linearizing operation is carried out using Taylor's theorem, and unknown values, such as X0, Y0, Z0, κ0, a1, a2, a3, a4, b1, b2, b3, and b4, are obtained. Further, using the obtained unknown values, a geometrically distorted satellite image is precisely corrected. In this case, the reason for assuming the functions of time φ and ω to be constants is that a high correlation exists between the functions φ and ω and the position variables of the satellite.
In this way, variables having high correlation are assumed to be constants and equations are solved in the prior art. In this case, if many ground control points are used, the remaining error of the satellite image can be reduced, but accurate position and attitude information of the satellite cannot be obtained.
Generally, if accurate position and attitude information of a satellite is obtained for a satellite image, more accurate position and attitude information of the satellite can be obtained for a satellite image adjacent to the processed satellite image by interpolation. However, in the prior art, accurate position and attitude information of the satellite cannot be obtained, so the obtained position and attitude information cannot be utilized for precision correction of the adjacent satellite image.
SUMMARY OF THE INVENTION
Accordingly, the present invention has been made keeping in mind the above problems occurring in the prior art, and an object of the present invention is to provide a method of precisely correcting geometrically distorted satellite images, which obtains accurate position and attitude information of a satellite by separating variables with a high correlation therebetween to obtain the solutions of non-linear equations, and which precisely corrects geometrically distorted satellite images using the position and attitude information of the satellite.
Another object of the present invention is to provide a computer-readable storage medium in which a program for executing the satellite image correcting method with a computer is stored.
In order to accomplish the above object, the present invention provides a method of precisely correcting geometrically distorted satellite images, using a satellite image and initial position and attitude values of a satellite at the time of photographing the satellite image received from the satellite, comprising the steps of a) establishing a geometric relationship between the satellite image and a ground control point as non-linear equations comprised of position variables and attitude variables of the satellite, and respectively modeling the position variables and the attitude variables to equations of time; b) initializing the attitude variables of the satellite to the initial attitude values provided by the satellite; c) applying the initial attitude values of the satellite and position variable equations of the satellite, obtained by the modeling, to the non-linear equations to obtain position values of the satellite; d) applying the obtained position values of the satellite and attitude variable equations of the satellite, obtained by the modeling, to the non-linear equations to obtain attitude values of the satellite; and e) precisely correcting geometric distortion of the satellite image using the position values and the attitude values of the satellite obtained at steps c) and d).



BRIEF DESCRIPTION OF THE DRAWINGS
The above and other objects, features and other advantages of the present invention will be more clearly understood from the following detailed description taken in conjunction with the accompanying drawings, in which:
FIG. 1 illustrates a typical geometric relationship between coordinates of a satellite image and a ground control point; and
FIG. 2 is a flowchart of a method of precisely correcting satellite images according to the present invention.



DESCRIPTION OF THE PREFERRED EMBODIMENTS
Hereinafter, embodiments of a method of precisely correcting geometrically distorted satellite images and computer-readable storage medium for the method according to the present invention will be described in detail with reference to the attached drawings.
The present invention is a method implemented through an algorithm executed by a typical computer system. The typical computer system comprises a typical input/output device, a microprocessor, a read only memory (ROM), a random access memory (RAM), and a database. The microprocessor entirely controls the computer system. The ROM stores an algorithm for precisely correcting satellite images according to the present invention and provides the algorithm to the microprocessor. The RAM stores data temporarily generated while the microprocessor executes the satellite image correcting method of the present invention. The database stores data required to execute the satellite image correcting method of the present invention.
The microprocessor executes the method of precisely correcting satellite images according to the present invention. In this case, the present invention separates variables having a correlation to minimize the correlation, and then calculates the values of the variables, thus accurately obtaining the position and attitude of a satellite.
Referring to FIG. 2, the method of precisely correcting satellite images according to the present invention is described in detail.
First, non-linear equations consisting of position variables and attitude variables of the satellite are established so as to correct a geometrically distorted satellite image at step S11, wherein the non-linear equations are geometric transforms of a distance between a camera of the satellite and a point on the ground. In the embodiment of the present invention, the geometric model of the linear pushbroom sensor, described as the example in the prior art, is used. However, the present invention is not limited to the geometric model, but can employ any equations as long as the equations are geometric transforms of a distance between camera sensor of a satellite and an arbitrary ground control point.
In the geometric model of the linear pushbroom sensor expressed in Equations [1] and [2], position values (Xs, Ys, Zs) of the satellite and attitude values (κs, φs, ωs) of the satellite are modeled as quadratic equations of time at step S12. The quadratic equations are arranged in the following Equations [5] and [6].




  
    
      
        
          
            
              
                X
                z
              
              =
              
                
                  X
                  0
                
                +
                
                  
                    a
                    1
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    1
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
          
            
              
                Y
                z
              
              =
              
                
                  Y
                  0
                
                +
                
                  
                    a
                    2
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    2
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
        
        
          
            
              
                Z
                z
              
              =
              
                
                  Z
                  0
                
                +
                
                  
                    a
                    3
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    3
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
          
            
              
            
          
        
      
    
    
      
        [
        5
        ]
      
    
  
  
    
      
        
          
            
              
                
                  
                    
                      κ
                      z
                    
                    =
                    
                      
                        κ
                        0
                      
                      +
                      
                        
                          a
                          4
                        
                        ⁢
                        t
                      
                      +
                      
                        
                          b
                          4
                        
                        ⁢
                        
                          t
                          2
                        
                      
                    
                  
                
              
              
                
                  
                    
                      φ
                      z
                    
                    =
                    
                      
                        φ
                        0
                      
                      +
                      
                        
                          a
                          5
                        
                        ⁢
                        t
                      
                      +
                      
                        
                          b
                          5
                        
                        ⁢
                        
                          t
                          2
                        
                      
                    
                  
                
              
            
          
        
        
          
            
              
                ω
                s
              
              =
              
                
                  ω
                  0
                
                +
                
                  
                    a
                    6
                  
                  ⁢
                  t
                
                +
                
                  
                    b
                    6
                  
                  ⁢
                  
                    t
                    2
                  
                
              
            
          
        
      
    
    
      
        [
        6
        ]
      
    
  




Next, the position values of the satellite of Equation [5] are initialized to (Xs0, Ys0, Zs0) using position values contained in auxiliary information provided by the satellite together with the satellite image, and the attitude values of the satellite of Equation [6] are initialized to (κs0, φs0, ωs0) using attitude values contained in the auxiliary information at step S13.
Thereafter, a variable “i” is set to “0” (the variable “i” is used to repeatedly perform the later steps S15 to S18) at step S14. The position variable equations of the satellite in Equation [5] and the attitude values (κsi, φsi, ωsi) of the satellite, converted into constants at step S13, are applied to the collinearity equations of Equations [1] and [2], where if “i” is “0”, the attitude values (κsi, φsi, ωsi) of the satellite become the initial attitude values (κs0, φs0, ωs0). Thereafter, the applied results are linearized using Taylor's theorem, and then unknown values, such as X0, Y0, Z0, a1, a2, a3, b1, b2, and b3 used to obtain geometrically corrected position values (Xsi+1, Ysi+1, Zsi+1) of the satellite, are obtained. Properly, position values (Xsi, Ysi, Zsi) of the satellite are given as initial values while the solutions of the collinearity equations are obtained. If the unknown values are obtained, the position values (Xsi+1, Ysi+1, Zsi+1) of the geometrically corrected satellite can be obtained at step S15.
Next, the position values (Xsi+1, Ysi+1, Zsi+1) of the geometrically corrected satellite, obtained at step S15, and the attitude variable equations of the satellite in Equation [6] are applied to the collinearity equations of Equations [1] and [2]. Thereafter, the applied results are linearized using Taylor's theorem, and then unknown values, such as t, κ0, φ0, ω0, a4, a5, a6, b4, b5, and b6 used to obtain the attitude values (κsi+1, φsi+1, ωs+1) of the geometrically corrected satellite, are obtained. Properly, even in this case, the attitude values (κsi, φsi, ωsi) of the satellite are given as initial values while the solutions of the collinearity equations are obtained. If the unknown values are obtained, the attitude values (κsi+1, φsi+1, ωsi+1) of the geometrically corrected satellite can be obtained at step step S16.
Next, a difference value between coordinates of the satellite image calculated using the geometrically corrected position values (Xsi+1, Ysi+1, Zsi+1) and geometrically corrected attitude values (κsi+1, φsi+1, ωsi+1) of the satellite respectively obtained at steps S15 and S16, and actual coordinates (X, Y) of the satellite image corresponding to an actual ground control point, is obtained at step S17. It is determined whether the difference value is less than a threshold value at step S18.
If it is determined that the difference value is less than the threshold value at step S18, the geometrically corrected position values (Xsi+1, Ysi+1, Zsi+1) and the geometrically corrected attitude values (κsi+1, φsi+1, ωsi+1) of the satellite, respectively obtained at steps S15 and S16, are applied to a typical precision correction method, thus enabling the geometrically distorted satellite image to be precisely corrected at step S19. Further, if it is determined that the difference value is not less than the threshold value at step S18, the variable “i” is increased by “1” at step S20, and the processing step returns to step S15.
As described above, the present invention provides a method of precisely correcting a geometrically distorted satellite image and computer-readable storage medium for the method, in which variables having a correlation are separated in non-linear equations represented by the position and attitude values of a satellite, and the non-linear equations are simplified to obtain solutions through respective linearizing operations for the separated variables, thus enabling the precise position and attitude of the satellite to be obtained.
Although the preferred embodiments of the present invention have been disclosed for illustrative purposes, those skilled in the art will appreciate that various modifications, additions and substitutions are possible, without departing from the scope and spirit of the invention as disclosed in the accompanying claims.

## Claims

1. A method of precisely correcting geometrically distorted satellite images, using a satellite image and position and attitude values of a satellite at the time of photographing the satellite image received from the satellite, comprising the steps of:
a) establishing a geometric relationship between the satellite image and a ground control point as non-linear equations comprised of position variables and attitude variables of the satellite, and respectively modeling the position variables and the attitude variables to equations of time;
b) initializing the attitude variables of the satellite to the initial attitude values provided by the satellite;
c) applying the initial attitude values of the satellite and position variable equations of the satellite, obtained by the modeling, to the non-linear equations to obtain position values of the satellite;
d) applying the obtained position values of the satellite and attitude variable equations of the satellite, obtained by the modeling, to the non-linear equations to obtain attitude values of the satellite; and
e) precisely correcting geometric distortion of the satellite image using the position values and the attitude values of the satellite obtained at steps c) and d), if a difference value between coordinates of the satellite image, calculated using the obtained position and attitude values of the satellite, and actual coordinates of the satellite image is less than a threshold value; and
f) repeatedly performing the steps c) and d) using the obtained attitude values of the satellite if the difference value therebetween is equal to or greater than the threshold value.

2. The method of precisely correcting geometrically distorted satellite images according to claim 1, wherein the step c) is performed such that the position variable equations of the satellite and the initial attitude values of the satellite are applied to the non-linear equations, the non-linear equations are linearized, and then the position values of the satellite are obtained.

3. The method of precisely correcting geometrically distorted satellite images according to claim 1, wherein the step d) is performed such that the attitude variable equations of the satellite and the position values of the satellite are applied to the non-linear equations, the non-linear equations are linearized, and then the attitude values of the satellite are obtained.

4. A computer-readable storage medium storing a program for executing the steps of:
a) establishing a geometric relationship between a satellite image and a ground control point as non-linear equations comprised of position variables and attitude variables of a satellite, and respectively modeling the position variables and the attitude variables to equations of time, if the satellite image and position and attitude values of the satellite at the time of photographing are received from the satellite;
b) initializing the attitude variables of the satellite to the initial attitude values provided by the satellite;
c) applying the initial attitude values of the satellite and position variable equations of the satellite, obtained by the modeling, to the non-linear equations to obtain position values of the satellite;
d) applying the obtained position values of the satellite and attitude variable equations of the satellite, obtained by the modeling, to the non-linear equations to obtain attitude values of the satellite; and
e) precisely correcting geometric distortion of the satellite image using the position values and the attitude values of the satellite obtained at steps c) and d), if a difference value between coordinates of the satellite image, calculated using the obtained position and attitude values of the satellite, and actual coordinates of the satellite image is less than a threshold value; and
f) repeatedly performing the steps c) and d) using the obtained attitude values of the satellite if the difference value therebetween is equal to or greater than the threshold value.