专利名称：遥感图像数据重采样方法
技术领域：
本发明涉及图像处理技术领域，特别涉及一种遥感图像数据重采样方法， 适用于对极轨气象卫星平台上光学遥感仪器获得的遥感图像的重采样。
背景技术：
数据重采样一般包括两个步骤1、像元坐标变换也即参与像元内插所对 应的原始像元坐标的确定，2、像元灰度值重新计算，有关像元灰度值的计算 方法有临近像元填充法、双线性内插法、双三次巻积法、正态采样法、反距 离加权法等；在重采样过程中，消耗时间最大的为步骤l像元坐标变换也即内
插点的确定，有关极轨卫星上仪器的数据重采样方法，《遥感学报》期刊2002 年3月第6巻第2期中'遥感图像数据重采样的一种快速算法' 一文中公开了 对半査找法逐点找到输出像元对应的原始点，该方法对每一个输出像元，都需 要在原始图像中查找多次，计算量大，对于一幅大小为NXM个像元的原始图 像，若采用对半查找方法，它的时间复杂度约为0(A^JkPbg2iVHog2M);对于 极轨气象卫星上的光学遥感仪器，其成像往往具有以下几个规律1、空间分 辨率不高，2、.扫描线边缘几何畸变大，3、仪器穿轨扫描，对于同一扫描线， 采样点序号与经纬度值之间的关系有一定的规律，如对于升轨图像，同一扫描 行内，随着采样点序号的增大，其经度值逐渐减小，处于西半球跨越到东半球 的像元除外；若从此规律出发，改进坐标变换算法将会提高速度，减少时间复 杂度。 发明内容如上所述，本发明的主要目的是为了建立一种快速确定内插点坐标的遥感 数据重采样方法以降低重采样时间复杂度提高重采样速度。
本发明的技术方案是提供一种遥感图像数据重采样方法，适用于对极轨 气象卫星平台上光学遥感仪器获得的遥感图像的重采样，包含下列步骤
1、 输入原始图象；
2、 根据输入图像的经纬度信息和空间分辨率大小计算目标图像的像元大 小和行列数；
3、 逐行根据相邻两像元经度值的大小关系确定参与输出像元内插的原始 图像像元坐标，然后进行线性内插输出等经度间隔的像素值数据集和纬度值数 据集；
4、 将步骤3输出的等经度间隔图像作为输入图像，对输入图像逐列根据 相邻两像元纬度值的大小关系确定参与输出像元内插的输入图像的像元坐标 然后进行线性内插；
5、 输出等经度间隔X等纬度间隔的图像。
步骤1中所述的输入图像为星上数据下传后经过定位处理含有定位信息的 HDF图像；
步骤2中，所述的输出等经度间隔X等纬度间隔的图像的像元大小的计算 公式为
pix = l度x IFOV(米)/100000(米) 公式中IFOV为像元空间分辨率大小，Pix为等经度间隔X等纬度间隔图 像的像元大小，计量单位为度；输出等经度间隔X等纬度间隔的图像的行列数 计算公式为-
Co/ = ce// ((丄max —丄min) / p/x )
6Z/"e = ce// ((_S max — 5 min) / / /x ) 其中，表示朝正方向取整，。/为列数，"恥为行数，图像经度最 大、最小值，纬度最大、最小值分别为Zmax，Z min,5max，5min 。
步骤3中，所述的逐行根据相邻两像元经度值的大小关系确定参与中间输
出的等经度间隔图像像元内插的原始图像像元坐标的方法及插值公式如下 令原始图像任一像元7在输出图像中的列坐标为;c(y)，对应的输出图像像
元列坐标为y(y)， x(y)， y(y)的计算公式分别为-
x(y') = (/o" (y) - Z min) / p/jc ， = ce// O(y))，其中，/o"为原始图像像
元经度值数组；
关系l、同一扫描行内，原始图像左右相邻两像元/， 7 + l的经度值大小 满足
/o" (_/) < /o" (_/ + 1); 在此种关系下，存在两种情况
A、 /卵(y)在西半球-180。附近，/。"(y + l)在东半球180°附近，经度值存
在跳跃；此时，原始像元与输出像元的空间坐标关系的确定及插值方法如下 1)、输出像元y())， y(j' + l)有直接映射点;c(y) , :c(y + l)，
a、 O) < y(y)，则参与y(y)插值的邻近点为j-1， j;输出像元像素值
.和纬度值的插值公式为
g"yC/)) = g(y -1) + -gC/' -1)) x (/o"丄CK力)-""C/ -1)) /(/。"(■/) - —1)) 其中^为输出图像像元的像素值或纬度值；g为输入图像像元的像素值或纬 度值；/o"为输出图像像元的经度值数组，对每行中像元列坐标为"的输出像 元的经度值计算公式为
/o"i(w)=丄min+(" —0.5)x p!'xb、 则参与少(y')插值的邻近点为y，+ 插值公式为-
g丄CK力)=g(力+ (gC/' +1) — x (/o"丄(y(/)) - /。"(_/ 》/(/。"(/ +1) — 360 — /o"(力)
2)、对于任一扫描线，越到扫描线边缘，几何畸变越明显，导致输出像元
y(y),j^+i)之间可能存在着没有原始像元对应的空白点，需进行空间坐标的确
定及插值，无空白点则不需要，对空白点的空间坐标确定及插值公式如下
a、 序号为l到y(/)之间对应的插值点为"+ l,插值公式为
= + (gC/ +1) - x (/一") - /—_/)) /(/ +1) - 3 60 - /—_/))
("-i,…X/)-i)
b、 序号为KZ + i)到。/之间对应的插值点为y， y + i;插值公式为
= g(y+1)+-g(/+1)) x - /。"(/+1)) /(3 60+/。"(/) - +1))
(少"-X7 + l) + l,…,Co/)
B、经度值不存在跳跃
1)、对于任一直接映射点X/),对应的原始像元的确定方法为
a、 O力')，则参与新像元力)插值的邻近点为j'J + l;插值公式为
g勒'))=g(Z)+(gC/+1) - x (,o""力))-/o"(力)/(/。"(_/ +1) - /。"(; ))
b、 力')"C/')，则参与新像元力)插值的邻近点为y-l,力插值公式
=g(y -1)+(gC/)-gC/-1)) x (/o""X;)) -/o"C/-禱o"(;') - -1))
2 )、相邻映射点jK几X/ +1)之间的空白点
a、 X/)与力'+ l)之间的存在着的无直接对应的原始像元点的空白像元，像
元序号从x/)+i到x;+i)-i的点对应参与插值的邻近点为，插值公式为
g丄O) = +(g(7+1) - x - /。</))/(/。</+1) - /。"(/))
("=></)+i,'",X/+i)-0 关系2、相邻两像元的经度值大小满足+在此种关系下，存在两种情况
A、 在东半球180°附近，/。"(j' + l)在西半球-180°附近，经度值存在跳
跃；
1)、对于直接映射点xy), ;^ + l),对应的原始像元的确定方法为
a、 力')〈X/)，则参与新像元力')插值的邻近点为W + 1;插值公式为
福力)=g(力+ +1) - x (/o"ZO(力)-/。</))/(360+/o"C/ +1) - /《))
b、 x(力〉乂力，则参与新像元X/)插值的邻近点为卜W;
=g(y -1)+-g(y -1)) x (/o""力.))-/《-1)) /(/—乂) - -1))
2)、空白点
a、 X/)与。/之间的存在着的无直接对应的原始像元点的空白像元，像元
序号从X力+ l到Co/的点对应参与插值的邻近点为乂y + l，插值公式为
由")=g(力+ +1) - x (/o"Z(y") - /—力)/(/on(y +1) + 360- /—;))
_y" = X/) + l，."，Co/ ;
b、 输出图像像元序号1至bcz+i)-l空白像元参与插值的原始像元点为
+ 插值公式为
gZ(一 = g(_/ +1) + -gC/ +1)) x (,o""") - +1)) /(/—力-3 60 - +1))
(戶=1,-.,力.+ 1)-1)
B、 经度值不存在跳跃
1)、对于任一直接映射点X力，对应的原始像元的确定方法为
a、 <力<7(/)，则参与新像元X/)插值的邻近点为7-W;插值公式
扭(力〕)=g(y -1)+-gc/ -1)) x (/o""x/》--1)) /(/—/) - /。"(/ -1))
b、 xC/)>K/)，则参与新像元x力插值的邻近点为力y+i;插值公式
=g(力+(g(/+1)—x (/。"z(Xy'))—/—力)/(/0</+1)—/。"(/))2)、相邻映射点力')，力+ l)之间的空白点
a、力+l)与力)之间的存在着的无直接对应的原始像元点的空白像元，像
元序号从x/)+i到xy+i)-l的点对应参与插值的邻近点为力y+i，插值公式为-gZ(X少"))=g(力+(gC/+1) - gC;.)) x —/o</))/(/。"(/+1)—/。《(_/))
(，X/'+i)+i,…,X/)-1)
以上所有情形中原始像元刚好落在输出像元中心上，即x(y')-y(y)则直接 用/填充y(/)，插值公式为
=， g丄(;^'+1)) = g(y+1)
逐行插值完后没有插值点的输出图像像元的纬度值都用所在行的纬度最大
值赋值；开辟一个数组将所有找到插值点的像元标记为1，记录输出图像中每
一列第一个和最后一个像元标记为1的点行号。
步骤(4)中，所述逐列利用上下邻近点内插过程中上下邻近点的确定及相
应的内插公式为
令等经度向隔图像像元/在输出的等经度间隔X等纬度间隔目标图像中的 坐标为z(/)，对应的输出图像像元坐标为^)， z(O、 g(/)的计算公式分别为
z(/) = (/a/X(/) — 5max)/pfx ; g。 = cC/(z(0) 其中，to￡为输入图像像元的纬度值；对每一列，循环从第一个标记为l 的点开始到最后一个像元标记为1的点为止；
同一列内，输入图像上下相邻两像元!'， / + 1的纬度值大小满足关系 A、 /fl^(0</a^(; + l)
1)、若</)<《(/)，则参与像元"/)插值的点为/-U;插值公式为
g"ew(柳)=g丄(/ _ 1) + (gZ:(z') - gZ(/ —1)) x to"e咖(/)) — /a即—1))/(/a柳一 to"!' -1))
g"w为输出图像的像素值数组，gZ:为输入图像的像素值数组，to"ew为输出图像的纬度值数组，对每列中行坐标为m的输出像元讳度值的计算公式为
/"加ew(附)二 5 max—(附—0,5) x / /x
2)、力)>9(/)，则参与像元《(/)插值的点为/,/ + 1;插值公式为. g"ew(《(/)) = gZ(/) + +1) - gZ(/)) x /"加ew(《(/)) 一 /a/Z(/)) /(/a仏(/ +1) — /a仏(/))
幼=柳，直接采用! 像元的值填充；
1)、 ？(/)<《(/)，则参与像元《(/)插值的点为/,/ + 1;插值公式为
g,( W》=柳)++1) — g柳x to"一(/))-/(/",Z(/ +1) - /a柳)
2) z(0>g(0，则参与像元g(/)插值的点为"U;插值公式为-.
g"ew(《(Z)) = g丄(/ — 1) + (g丄(/) 一 gZ(/ —1)) x /a加ew(《(/)) — /WZ(/ —1)) — to,Z(/ -1))
z(!〕 = g(/)，直接采用/像元的值填充；
步骤5中，输出经纬度间隔都为;^的等经度间隔X等纬度间隔图像。
本发明的优点坐标变换过程中不需要对每一个点进行多次查找，对于一 幅大小为N*M的图像处理的时间复杂度约为0(2*^*^)，时间复杂度明显小 于对半査找方法；同时输出的等经度间隔X等纬度间隔图像可供选择多种投影 方式满足不同用户的需要；还能在重采样过程中解决如45°镜加多元并扫引起 的图象像旋现象。
图1、流程图。
图2、实施例中输入的原始图。 图3、实施例经重采样后的输出图。 具体实施方案
下面根据图1-图3给出本发明一个较好实施例，详述如下1、选择实例所用卫星数据为风云三号气象卫星上搭载的中分辨率光谱成
像仪于UTC时间2008年7月6日01时55分通道18所获得的一幅5分钟HDF 图2作为处理流程1中步骤1的输入图像；由图2中的太湖明显可见，图 像几何畸变以及像旋现象非常明显；
2、 剔除缺失行，.图象大小为1980行2048列;
3、 该图空间分辨率大小为IOOO米，步骤2中计算得目标图像的像元大小 Pix值为0.01度；计算图像的经度最大最小值分别为149°， 111°，纬度最大最 小值分别为46。、 23°，利用公式
Co/ = ce// ((丄max — Z min) / j /;c) 丄/we = ce// (( 5 max 一 5 min) / / /x )
计算输出图像的行列数为2300行3800列；
4、 逐行内插过程开始(步骤3):
判断相邻两像元的经纬度值大小关系如果/朋U)^/。"(y + l)，但/。0)在 东半球180。附近，/。"(/+1)在西半球_180°附近，即扫描从东半球跳至西半球， 经度值存在跳跃，首先判断;c(刀与y(7)的关系，确定y(乂)的插值点和插值公
式
1) 、如果0) < ;K7)，则用+ i (pw， w为原始图象总采样点数，若_/ = ^
直接用y填充)的纬度值和像素值内插得到少(/)的像素值和纬度值，插值公式
为<formula>formula see original document page 12</formula>然后对空白点插值，如果Co/ -少(/)^l表示序号为^C/) + l到Co/的所有点 都为空白点，采用原始像元力j + l内插，内插公式为
=gc/+1)++1) - * (/o"ick") - /《+1)) /(/《+1)+3 60 - +1))
(，X/) + l,…，Co/)
如果力+1)>1，参与输出图像像元序号i至iJxz+i)-i的空白像元插值的原
始像元点为/J + l;插值公式为
处")=+1)+, - g(/+1)) x (/o"ZO) - /0</+1)) /崎)-360- /—/+1))
(yw = V",X/+l)-l)
对于经度值关系满足/卵(y) >g +1)的其他情况也即经度值不存在跳跃
的情况，则首先判断x(y)与y(/)的关系，确定^(/)的插值点和插值公式
1) 、如果x(y) > ;Ky')，则用力"i (_/W， iv为原始图象总采样点数，若j:iv 直接用y填充)的纬度值和像素值内插得到少(y)的像素值和纬度值，插值公式
为
g丄(力.))=g(力+ (gG +1) — x (〖0"1(力》-/(/o心+1) - /o"C0)
2) 、如果;c(7)q(/)，则用卜W (7."，若^1直接用卜1填充)的纬度 值和像素值内插得到少G)的像素值和纬度值，插值公式为
=g(y -1)+- g(y -1)) x (/o""力))--1)) /(/—_/) - /。《0' -1))
然后判断:Ky)与yG + i)之间是否存在空白点，如果;K7 + i) - :K0 > i表示 输出图像像元序号为y(y) +1到y(7 +1) -1的所有点都为空白点，采用像元
y,y + l内插，内插公式为
=g(力+(g(/+1) - x (/o""少")-/o o〕) /(/o"(y+1) -
(_y" = X/) + l,H.，X/ + l)-l)
如果相邻两像元的经度值大小关系满足/o" (_/) < /o" (7 + 1)，但/。</ + 1)在东半球180°附近，/。n(力在西半球-180°附近，经度值存在跳跃；则首先判断;c(乂)
与少C/)的关系，确定少(y)的插值点和插值公式
1)、如果xC/)^C0，则参与输出图像像元X/)插值的邻近点为+ 1 ( _/W， W为原始图象总采样点数，若y-W直接用y填充)；插值公式为
g"X力)=+ (gC/ +1) - gC/)) x (/o""X力)-W)) /(/。《0' +1) - 3 60 - /。</)) 2、 </ )<></)，则参与输出像元X力插值的邻近点为/-W (Wl，若/ = 1直
接用卜l填充)；插值公式为
= i)+(gC/)—gC/—1))x(/o"":K/))—1))
然后计算空白点，如果c。/ - w + i)〉i表示序号为yO + i) + i到c。/的所
有点都为空白点，采用原始像元y,y+i内插，内插公式为
g丄O) = g(y+1)+- g(;+1)) x- /《+1)) /(360+/o"co - /—乂+1))
(，-XZ+i)+i,…,Co/)
如果X/)w，参与输出图像像元序号1至Ux/)-i的空白像元插值的原始像 元点为丄j+l;插值公式为
WCy") = gC0+(gC/+i)-gC/))x(/o""")-/o</)Wo</+i)-360-/。"(/)).
(少"-l,…,x/)-l)
对于经度值关系满足/0"(/)</0"(_/ + 1)的其他情况，即经度值不存在跳
跃
1) 、若x(y)〉力)，则用y-w 若_/ = 1直接用卜i填充)插值，插值
公式为
g"x/))=gC/'-i)+(gc/)-ga-i))x—io</'))-,o o o')-1))
2) 、力')<力')，则参与输出图像像元X/)插值的邻近点为"'+l;插值公式
gz(x/)) =+(gC/+1) - x —z(X/)) - W)) /(W+1) - /。</))然后判断y(y)与力'+1)之间是否存在空白点，如果力.)—w +1) > i表示 输出图像像元序号为xz+i)+i到力vi的所有点都为空白点，采用像元az+i内 插，内插公式为
' (，K/)+i,…,xy+i)-1)
以上所有情形中原始像元刚好落在输出像元中心上，即;cC/) = yC/)则直接
用/填充y(y)，插值公式为:woc/))=g(y)， giw+i))=gc/+i);逐行插值完后
没有插值点的输出图像像元的纬度值都用所在行的讳度最大值赋值，开辟一个 数组将所有找到插值点的像元标记为1，记录输出图像中每一列第一个和最后
一个像元标记为1的点行号；
5、逐列内插开始(步骤4):
对每一列，循环从第一个标记为1的点开始到最后一个像元标记为1的点
为止，逐点计算像元/在输出图像中的列坐标z(/)和对应输出图象的像元列坐标 W)，计算公式分别为
z(/) = (5max— toZ(/))/p/jc ， 《(/) = ce//(z(/))
判断它与相邻点/ + 1的纬度值关系，如果/加印)</^1^+1)即9(/)>^+1)，则
首先判断z(/)与《(/)的关系，确定^)的插值点和插值公式
1)、如果z(/) 则用/-U ("1，若/ = 1直接用7点填充)的纬度值和
像素值内插得到《(i)的像素值，插值公式为
g"—《(/))=洲-1) + (洲)—gZ(/ — 1》x /a/w—柳)—-1)) /(/a柳—/fl邵—1))
2)、如果2(/)>^)，则用W + 1 ("""e，若 丄^直接用/ + 1点填充)的
纬度值和像素值内插得到W)的像素值和炜度值，插值公式为 g"ew(《(/)) = gZ(/) + (gZ(/ +1) — g丄(/)) x /a加ew(《(/)) 一 /aA(/))/(/W丄(/ +1) -然后判断《(/ + 1)与^)之间是否存在空白点，如果《(刀-《(/ + 1)>1表示 ^ + l)到《(/)之间所有的点都为空白点，空白点处于原始像元/,"l之间，因此 采用像元/,/ + 1内插，内插公式为
g"evi(，) = g丄(/) + (g(/ +1) — g(/)) x (/a/"eM(，) — toZ(/)) /(/" W +1) — /。 W))
(拜"(/ + 1) + 1，-..,柳-1)
如果是/"仏(/)>/"^(; + 1)即^)<^ + 1)，则首先判断z(/)与《(/)的关系，确定 W)的插值点和插值公式
1) 、如果z(/) 则用/,/ + 1 ("Z/""若/ =丄/^直接用/ + 1点填充)的 纬度值和像素值内插得到^)的像素值，插值公式为
gwew(g0.)) = g厶(z.) + (g丄(/ +1) — g丄(/)) x /flmew(《(6) - /a/丄0')) /(/a"(i +1) - /"/丄(i))
2) 如果"0>^/)，则用卜1,/ ("1，若 l直接用/点填充)的纬度值和
像素值内插得到《(/)的像素值和讳度值，插值公式为
g"ew(《(/)) = g丄(/ 一 1) + (g丄(/) — g"/ —1)) x /W"evv(《(/)) — /a/Z(/ —1)) /(/a,丄(/) — /a,丄(/ -1))
然后判断^)与《(；+ l)之间是否存在空白点，如果《(/ + 1)-9(/)> 1表示《(/ + 1) 到q(/)之间所有的点都为空白点，空白点处于原始像元U + 1之间，因此采用像 元/,/ + 1内插，内插公式为^
gwen(戸)=柳)+ (g(/ +1) — g(/)) x (/a/"ew(，) — /。,"/)) +1) —
(拜=柳+1, .，《(/ +1) -1)
对于z(/) =的像元直接采用像元/填充像元g(i);
6、步骤5输出消除了像旋的等经度间隔X等纬度间隔3，由图3可看 出像旋现象已消除，图像严重压縮的情况己不存在
权利要求
1. 一种遥感图像数据重采样方法，其特征在于它包括以下步骤A. 输入星上数据下传后经过定位处理含有定位信息的HDF图像；B. 根据输入图像的经纬度信息和空间分辨率大小计算输出目标图像的像元大小和行列数，像元大小的计算公式为pix＝1度×IFOV(米)/100000(米)其中，IFOV为像元空间分辨率大小，pix为目标图像的像元大小，计量单位为度；目标图像的行列数计算公式为Col＝ceil((L max-L min)/pix)Line＝ceil((B max-B min)/pix)其中，ceil表示朝正方向取整，Col为列数，Line为行数，图像经度最大值、最小值，纬度最大值、最小值分别为L max，L min，B max，B min；C. 逐行根据相邻两像元经度值的大小关系确定参与中间输出的等经度间隔图像像元内插的原始图像像元坐标，然后进行线性内插输出等经度间隔的像素值数据集和纬度值数据集；D. 将步骤C输出的等经度间隔图像作为输入图像，对输入图像逐列根据相邻两像元纬度值的大小关系确定参与输出像元内插的输入图像的像元坐标然后进行线性内插；E. 输出等经度间隔×等纬度间隔的图像。
2. 根据权利要求1所述的一种遥感图像数据重采样方法，其特征在于所 述的步骤C中逐行根据相邻两像元经度值的大小关系确定参与中间输出的等 经度间隔图像像元内插的原始图像像元坐标的方法如下令原始图像任一像元y在输出图像中的列坐标为x(y)，对应的输出图像像元列坐标为^(y), x(y)，少(y)的计算公式分别为<formula>formula see original document page 3</formula> 其中，为原始图像像元经度值数组；(1) 同一扫描行内，原始图像相邻两像元h j + l的经度值大小满足<formula>formula see original document page 3</formula>在此种关系下，存在两种情况情况A:扫描线从西半球跨越到东半球，/。《(;')在西半球-180°附近， /o"(y' + l)在东半球180°附近，经度值存在跳跃，此时像元y， y' + l对应的映 射点y(y)和+ 1)以及处在两映射点之间空白像元的内插点坐标确定如下若<，则参与输出图像像元yC/)插值的原始图像像元坐标为y'-i，_/;若1(_/)>:^_/)，—则参与少(7')插值的原始像元坐标为y， + 1 ;若 x(y) = K/)，则直接用y填充"y);输出图像像元序号l至iJ力)-i和序号 力+ l) + l到Co/的所有像元为映射点之间的空白点，对应的插值点为j, y + i;情况B:相邻两点经度值不存在跳跃，若</)<力)，则参与X/)插值的邻近点为_/,_/+1;若xC0"a)，则参X/)插值的邻近点为卜W;若4y) = 则 直接用7填充y(/);像元序号从xy)+i到XZ+i)-i的点对应参与插值的点为 A/.+i;(2) 同一扫描行内，原始图像相邻两像元y， _/ + 1的经度值大小满足<formula>formula see original document page 3</formula>在此种关系下，存在两种情况情况A:扫描从东半球跨越到西半球，/。"(y)在东半球180。附近，/0"c/+i) 在西半球-180°附近，经度值存在跳跃；若x(yxx刀，则参与输出图像像元X/)插值的邻近点为力y+i;若力')>^/)，则参与输出图像像元:K刀插值的邻近点为 若</) = ></)，则直接用7填充少C/');像元序号:K/)+i到。/和1到 i的像元参与插值的原始像元点为;情况B:经度值不存在跳跃，若</)<少(力，则参与输出图像像元X/)插值 的邻近点为卜W;若x(/oX力，则参与输出图像像元X/)插值的邻近点为 yj + l;若Oy(/)，则直接用/填充yC/);对于相邻映射点X/)， XZ + D之间的空白点，像元序号从力'+i)+i到X/)-i的点对应参与插值的邻近点为
3.根据权利要求1所述的一种遥感图像数据重采样方法，其特征在于所 述的步骤D中逐列根据相邻两像元纬度值的大小关系确定参与输出图像像元 内插的输入图像的像元坐标然后进行线性内插的方法如下令等经度间隔图像像元/在输出的等经度间隔X等纬度间隔目标图像中的 坐标为W)，对应的输出图像像元坐标为《(/)， z(/)、 g(/)的计算公式分别为z(/) = (to丄0') — 5 max) / / /x=其中，to丄为输入图像像元的绵度值；同一列内，输入图像上下相邻两像元/， "l的纬度值大小满足关系-若力)>^)，则参与像元《(/)插值的点为^ + l;若力')<^)，则参与像元《。插值的点为/-U;若2 =《(/)，直接采用/像元的值填充；同一列内，输入图像上下相邻两像元/, / + 1的纬度值大小满足关系若2(0<9(0，则参与像元^)插值的点为M+1;若4!')>9(!')，则参与像元W)插 值的点为/-U;对于^/)=《(/)，直接采用/像元的值填充。
全文摘要
本发明公开一种遥感图像数据重采样方法，适用于对极轨气象卫星平台上光学遥感仪器获得的遥感图像的重采样。本发明先对输入的原始图像逐行根据同一扫描行内相邻两点经度值的大小关系确定参与输出像元内插的原始图像点坐标，然后进行等经度间隔线性内插，接下来对输出的等经度间隔数据逐列根据上下邻近点的纬度值大小关系确定参与输出像元内插的输入像元的坐标进行线性内插，输出等经度间隔×等纬度间隔的图象数据。本发明所提供的方法减少了确定内插点坐标的计算量，对于有像旋现象的图像可同时消除像旋。本发明方法输出的图像为像元大小规则的图像，为接下来图像进行适当的投影打下了基础，进而为卫星遥感数据进一步的数据处理与应用奠定了基础。
文档编号G06T3/40GK101533510SQ20091004911
公开日2009年9月16日 申请日期2009年4月10日 优先权日2009年4月10日
发明者荣 叶, 球 尹, 巩彩兰, 勇 胡, 陈海燕 申请人:中国科学院上海技术物理研究所