| 主权项 |
一种环境一号卫星CCD影像正射校正方法，其特征在于步骤如下：(1)导入环境一号卫星CCD的一级影像；导入格网间距为M1米的SRTM高程数据与分辨率为M2米的LandSat全色参考数据；所述M1的范围为[30,150]，M2的范围为[7.5,30]；(2)在步骤(1)导入的一级影像中均匀采集N个地面控制点，获得该控制点在参考影像的空间坐标(P，L，H)和控制点在一级影像中的像素坐标(X，Y)；所述N大于等于19；(3)将步骤(2)采集的地面控制点在参考影像的空间坐标(P，L，H)和控制点在一级影像中的像素坐标(X，Y)导入有理函数模型，建立4N个等式，计算出有理多项式系数，即RPC，确定有理函数模型，具体计算公式如下：<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>P</mi><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>,</mo><mi>H</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>,</mo><mi>H</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mi>L</mi><mo>=</mo><mfrac><mrow><msub><mi>R</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>,</mo><mi>H</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>R</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>,</mo><mi>H</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000623308860000011.GIF" wi="421" he="357" /></maths>其中，P、L和H分别为采集的控制点在参考影像的横坐标、纵坐标和高程值；X和Y分别为采集的控制点在一级影像中的像素横坐标和像素纵坐标；P<sub>1</sub>、P<sub>2</sub>、P<sub>3</sub>和P<sub>4</sub>均为形式相同的多项式，展开后分别为：P<sub>1</sub>＝a<sub>1</sub>+a<sub>2</sub>X+a<sub>3</sub>Y+a<sub>4</sub>H+a<sub>5</sub>XY+a<sub>6</sub>YH+a<sub>7</sub>XH+a<sub>8</sub>Y<sup>2</sup>+a<sub>9</sub>X<sup>2</sup>+a<sub>10</sub>H<sup>2</sup>+a<sub>11</sub>XYH+a<sub>12</sub>Y<sup>3</sup>+a<sub>13</sub>YX<sup>2</sup>+a<sub>14</sub>YH<sup>2</sup>+a<sub>15</sub>Y<sup>2</sup>X+a<sub>16</sub>X<sup>3</sup>+a<sub>17</sub>XH<sup>2</sup>+a<sub>18</sub>Y<sup>2</sup>H+a<sub>19</sub>X<sup>2</sup>H+a<sub>20</sub>H<sup>3</sup>P<sub>2</sub>＝b<sub>1</sub>+b<sub>2</sub>X+b<sub>3</sub>Y+b<sub>4</sub>H+b<sub>5</sub>XY+b<sub>6</sub>YH+b<sub>7</sub>XH+b<sub>8</sub>Y<sup>2</sup>+b<sub>9</sub>X<sup>2</sup>+b<sub>10</sub>H<sup>2</sup>+b<sub>11</sub>XYH+b<sub>12</sub>Y<sup>3</sup>+b<sub>13</sub>YX<sup>2</sup>+b<sub>14</sub>YH<sup>2</sup>+b<sub>15</sub>Y<sup>2</sup>X+b<sub>16</sub>X<sup>3</sup>+b<sub>17</sub>XH<sup>2</sup>+b<sub>18</sub>Y<sup>2</sup>H+b<sub>19</sub>X<sup>2</sup>H+b<sub>20</sub>H<sup>3</sup>P<sub>3</sub>＝c<sub>1</sub>+c<sub>2</sub>X+c<sub>3</sub>Y+c<sub>4</sub>H+c<sub>5</sub>XY+c<sub>6</sub>YH+c<sub>7</sub>XH+c<sub>8</sub>Y<sup>2</sup>+c<sub>9</sub>X<sup>2</sup>+c<sub>10</sub>H<sup>2</sup>+c<sub>11</sub>XYH+c<sub>12</sub>Y<sup>3</sup>+c<sub>13</sub>YX<sup>2</sup>+c<sub>14</sub>YH<sup>2</sup>+c<sub>15</sub>Y<sup>2</sup>X+c<sub>16</sub>X<sup>3</sup>+c<sub>17</sub>XH<sup>2</sup>+c<sub>18</sub>Y<sup>2</sup>H+c<sub>19</sub>X<sup>2</sup>H+c<sub>20</sub>H<sup>3</sup>P<sub>4</sub>＝d<sub>1</sub>+d<sub>2</sub>X+d<sub>3</sub>Y+d<sub>4</sub>H+d<sub>5</sub>XY+d<sub>6</sub>YH+d<sub>7</sub>XH+d<sub>8</sub>Y<sup>2</sup>+d<sub>9</sub>X<sup>2</sup>+d<sub>10</sub>H<sup>2</sup>+d<sub>11</sub>XYH+d<sub>12</sub>Y<sup>3</sup>+d<sub>13</sub>YX<sup>2</sup>+d<sub>14</sub>YH<sup>2</sup>+d<sub>15</sub>Y<sup>2</sup>X+d<sub>16</sub>X<sup>3</sup>+d<sub>17</sub>XH<sup>2</sup>+d<sub>18</sub>Y<sup>2</sup>H+d<sub>19</sub>X<sup>2</sup>H+d<sub>20</sub>H<sup>3</sup>式中，a<sub>j</sub>,b<sub>j</sub>,c<sub>j</sub>,d<sub>j</sub>，j＝1,2,...,20为有理多项式系数；(4)将环境一号卫星CCD的一级影像所有像素坐标(P，L)以及对应高程值H代入步骤(3)中确定的有理函数模型，得到校正后图像的坐标值(X，Y)；(5)采用内插值方法对步骤(4)中得到的校正后图像的坐标值(X，Y)进行处理，得到整数坐标值像素的灰度值，从而得到校正后的环境一号卫星CCD影像。 |
