辐射定标

辐射定标（像元亮度值，辐射亮度/亮温）、表观反射率、地表反射率、反照率、比辐射率(转)
(2012-11-28 13:58:29)
转载▼
分类：科研
标签：
杂谈
(2012-01-26 01:18:44)
标签：
校园分类：工作篇
定标系数为：增益53.473，单位：DN/（W⋅m-2⋅sr-1⋅μm-1）；截距26.965，单位：DN。

利用绝对定标系数将DN值图像转换为辐亮度图像的公式为L=（DN-b）/coe，式中coe为绝对定标系数的增益，b为截距，转换后辐亮度单位为
W⋅m-2⋅sr-1⋅μm-1。

HJ1B红外相机中红外波段则条带较为严重，不利于定量化应用。

遥感数字图像
遥感数字图像是以数字形式记录的二维遥感信息，即其内容是通过遥感手段获得的，通常是地物不同波段的电磁波谱信息。

其中的像素值称为亮度值（或称为灰度值、DN值）。

遥感概念DN值（Digital Number ）是遥感影像像元亮度值，记录的地物的灰度值。

无单位，是一个整数值，值大小与传感器的辐射分辨率、地物发射率、大气透过率和散射率等有关。

遥感图像量化image quantification。

释文：按一定的函数关系将图像所代表的物理量分割成有限的离散等级,以使观测数据可用一定字长的二进制码表示,因此又称为数据编码。

量化后的级别称为图像的像元值、灰度或亮度,记为
DN(digital number)。

DN值没有单位，数量级与像素深度有关，如果是无符号整型的就是0-255，符点型，无符号16位均根据其类型确定。

在遥感领域，定标一般分为几何定标和辐射定标两种。

几何定标即指对遥感图像几何特性进行校正，以还原为真实情况。

辐射定标指对遥感图像的辐射度进行校准，以实现定量遥感。

辐射定标一般也可称为校准，其主要目的是保证传感器获取遥感数据的准确性。

通常，采用系统自身内部监视环路和外部标准目标方法对系统链路中的各个环节进行误差修正，来实现辐射定标过程。

一般在主动式遥感系统中，辐射定标可以作得很好，可以认为在一定误差范围内实现了定量遥感。

而被动式遥感系统相对困难些。

几何定标相对简单，就不多说了。

辐射定标是对传感器引起的误差校正，将影像校正为星上反射率
辐射定标和辐射校正——遥感数据定量化的最基本环节
由于遥感图像成像过程的复杂性，传感器接收到的电磁波能量与目标本身辐射的能量是不一致的。

传感器输出的能量包含了由于太阳位置、大气条件、地形影响和传感器本身的性能等所引起的各种失真，这些失真不是地面目标的辐射，因此对图像的使用和理解造成影响，必须加以校正和消除，而校正和消除的基本方法就是辐射定标和辐射校正。

辐射定标是指传感器探测值的标定过程方法，用以确定传感器入口处的准确辐射值。

辐射校正是指消除或改正遥感图像成像过程中附加在传感器输出的辐射能量中的各种噪声的过程。

一般情况下，用户得到的遥感图像在地面接收站处理中心已经作了辐射定标和辐射校正。

大气对电磁波的作用基本可以归纳为两种物理过程，即吸收和散射。

对环境遥感来说，大气的吸收与散射作用均可使电磁波受到削弱，此谓大气效应。

大气纠正，即大气效应纠正，其基本目标是＂如何将星载传感器所测得的辐射亮度值（ＤＮ）中的大气效应消除掉，从而获得地表真实的辐射亮度值＂。

辐射校正是指消除或改正遥感图像成像过程中附加在传感器输出的辐射能量中的各种噪声的过程。

辐射校正包括三个方面：
１、影像的辐射校正
２、太阳高度角和地形影响引起的辐射误差校正
３、大气校正
大气校正只是辐射校正的一个方面。

光学遥感图像例如TM的一些辐射定标工作，总的来说可以有广义和狭义之分：狭义的是指将图像的DN值转化为有物理量纲的观测目标大气上界的亮度值或者反射率
广义的除了狭义的步骤外还包括去除大气影响，计算地表目标真实的反射亮度或者是反射率
假如用６Ｓ模型校正，第一步是辐射定标，根据ＤＮ值、影像的ｏｆｆｓｅｔ和ｇａｉｎ计算辐射量度Ｌ；第二步用６Ｓ模型得到大气校正参数ｘａ、ｘｂ、ｘｃ；第三步根据三个大气校正参数和Ｌ就可计算校正后的反射率。

从计算方程来看，如果具有相同ＤＮ值的像元，校正后其反射率也相同。

地表反射电磁波，被卫星的传感器记录下来就得到DN值，DN值进过定标和相关的公式能够转变为地表的反射率值。

DN值就是遥感传感器的数字量化输出值，相当于图像的灰度值。

可以根据记录的原始DN值进行辐射定标，转换为大气外层表面反射率。

DN值就是像素值。

DN值的大小代表地物反射电磁波的能力，通俗点就是分辨地物。

因为不同的地物应该有不同的DN值。

表观反射率
用多辐射校正水平遥感数据提取植被叶面积指数的精度分析
顾祝军刘咏梅陆俊英南京晓庄学院生物化工与环境工程学院南京211171 土壤与农业可持续发展国家重点实验室（中国科学院南京土壤研究所）南京210008 选用南京市SPOT5 HRG图像的地物反射率（PAC）、表观反射率（TOA）、星上辐射率（SR）和灰度值（DN）影像,提取了2种植被指数（VI）,即归一化植被指数
（NDVI）和比值植被指数（RVI）,与地面实测的植被叶面积指数（LAI）进行了相关分析,并建立了157个LAI-VI关系模型。

结果显示：LAI与VI呈显著的正相关关系（r=0.303~0.927,p〈0.01）,对应不同植被的优选模型自变量包括了3个辐射校正水平的2种植被指数,可见基于不同辐射校正水平的植被指数在LAI 遥感反演中具有一定的应用潜力。

这些优选模型为：阔叶林：LAI=-3.345＋
5.378RVISR＋7.329NDVISR（R2=0.818,RMSE=0.527）,针阔混交林：LAI=1.696＋17.076NDVIDN＋137.684（NDVIDN）2-288.240（NDVIDN）3
（R2=0.919,RMSE=0.440）,灌木：LAI=-0.065＋19.112NDVISR-113.820（NDVISR）2＋184.207（NDVISR）3（R2=0.900,RMSE=0.448）,草地：LAI=-5.905＋6.446RVISR ＋9.477NDVISR（R2=0.944,RMSE=0.378）,植被总体：LAI=-1.615＋7.199NDVIDN ＋2.640NDVISR＋2.105RVIPAC（R2=0.801,RMSE=0.668）。

研究表明,基于不同植被类型、不同辐射校正水平影像的LAI遥感估算有利于充分挖掘遥感影像信息,进而提高LAI估算的精度。

由于植被遥感应用定量化和监测等的需求，光学遥感数据的辐射校正更加受到重视该文论述了辐射校正，辐射定标和大气校正的慨念以及它们之间的区别及关系特别对辐射定标的结果之一，大气层顶表观反射率，简称表观反射率(Apparent reflectance)的定义、慨念、计算和它在植被遥感中的应用等方面，进行了详细的论述。

暗像元算法基于表观反射率的大气贡献项,即利用卫星观测的路径辐射反演气溶胶光学厚度。

它是目前陆地上空气溶胶遥感应用最为广泛的算法。

遥感反射率的定义：地物表面反射能量与到达地物表面的入射能量的比值。

遥感表观反射率的定义：地物表面反射能量与近地表太阳入射能量的比值。

大气校正就是将辐射亮度或者表观反射率转换为地表实际反射率，目的是消除大气散射、吸收、反射引起的误差。

1、反射率：是指任何物体表面反射阳光的能力。

这种反射能力通常用百分数来表示。

比如说某物体的反射率是45%，这意思是说，此物体表面所接受到的太阳辐射中，有45%被反射了出去.英文表示：Reflectance
2、地表反射率：地面反射辐射量与入射辐射量之比，表征地面对太阳辐射的吸收和反射能力。

反射率越大，地面吸收太阳辐射越少；反射率越小，地面吸收太阳辐射越多，表示：surface albedo
3、表观反射率：表观反射率就是指大气层顶的反射率,辐射定标的结果之一，大气层顶表观反射率，简称表观反射率，又称视反射率。

英文表示为：apparent reflectance（=地表反射率+大气反射率。

所以需要大气校正为地表反射率）。

‚5S‛和‚6S‛模型输入的是表观反射率而MODTRAN模型要求输入的是辐射亮度。

4、行星反射率：从文献‚一种实用大气校正方法及其在ＴＭ影像中的应用‛中看到‚卫星所观测的行星反射率（未经大气校正的反射率）‛；在‚基于地面耦合的TM影像的大气校正-以珠江口为例‛一文有‚该文应用1998年的LANDSAT5
TM影像,对原始数据进行定标、辐射校正,求得地物的行星反射率‛。

因此行星反射率就是表观反射率。

英文表示：planetary albedo
5、反照率：反照率是指地表在太阳辐射的影响下，反射辐射通量与入射辐射通量的比值。

它是反演很多地表参数的重要变量，反映了地表对太阳辐射的吸收能力。

英文表示：albedo
它与反射率的概念是有区别的：反射率（reflectance)是指某一波段向一定方向的反射，因而反照率是反射率在所有方向上的积分；反射率是波长的函数，不同波长反射率不一样，反照率是对全波长而言的。

反照率的定义是地物全波段的反射比,反射率为各个波段的反射系数。

因此,反照率为地物波长从0 到∞的反射比。

6.地表比辐射率(Surface Emissivity)，又称发射率，指在同一温度下地表发射的辐射量与一黑体发射的辐射量的比值，与地表组成成分，地表粗糙度，波长等因素有关。

比辐射率的直接测量。

理论上,比辐射率的测定有两种
途径,一种是比色法,这种方法目前只能使用在被测物的温度大于50 ℃的场合。

因为信噪比太小,不适合常温地球表面的测量。

然而,随着传感器技术的发展,如果能测量零度以下物体的话,这种比色法似可取得突破性的发展; 另一种是亮度法。

也是目前人们所采用的办法。

在实验室里,利用封闭式黑体筒可以成功地测量地物的比辐射率。

也可以利用主动和被动相结合的方法测量比辐射率,这种方法已在实验室里取得成功。

利用二氧化碳激光,可以远距离测量地物的比辐射率,目前,已经开始把这一技术向航空和航天遥感扩展,它的可行性已经得到证实,其目标是对区域范围的地物比辐射率进行直接测定。

我们深信这种高技术的实现已为期不远了。

这种比辐射率的直接测定,不仅可以直接获得比辐射率区域分布,而且可以获得比辐射率的多角度以及地物性质的有关信息。

这种研究思路的实现,对定量热红外遥感的推动作用是巨大的。

辐射亮度表示面辐射源上某点在一定方向上的辐射强弱的物理量。

辐射亮度的SI单位为瓦/( 球面度.米2 ) 。

面辐射源的辐射亮度L：辐射源在某一方向的单位投影面积在单位立体角内的辐射通量，称为辐射亮度（Radiance）L。

Le=dLe/dScosθθ为给定方向和辐射源面元法线间的夹角。

单位是W/(sr*m2)。

辐射
基本内容
表示面辐射源上某点在一定方向上的辐射强弱的物理量。

辐射亮度的SI单位为瓦／( 球面度.米2 ) 。

面辐射源的辐射亮度L：辐射源在某一方向的单位投影面积在单位立体角内的辐射通量，称为辐射亮度
（Radiance）L。

单位是W/(sr*m2)。

自然界中的一切物体,只要温度在绝对温度零度以上,都以电磁波的
形式时刻不停地向外传送热量,这种传送能量的方式称为辐射。

物体通过辐射所放出的能量，称为辐射能。

辐射按伦琴/小时(R)计算辐射有一个重要的特点，就是它是‚对等的‛。

不论物体(气体)温度高低都向外辐射，甲物体可以向乙物体辐射，同时乙也可向甲辐射。

这一点不同于传导，传导是单向进行的。

任何已经遭遇辐射的人都应用肥皂和大量清水彻底冲洗整个身体,并立即寻求医生或专家的帮助 !（图为"放射性物质危险，小心辐射‛的警示标志）
辐射能被物体吸收时发生热的效应，物体吸收的辐射能不同，所产生的温度也不同。

因此，辐射是能量转换为热量的重要方式。

辐射传热（radiant heat transfer）依靠电磁波辐射实现热冷物体间热量传递的过程，是一种非接触式传热，在真空中也能进行。

物体发出的电磁波，理论上是在整个波谱范围内分布，但在工业上所遇到的温度范围内，有实际意义的是波长位于0.38～1000μm之间的热辐射,而且大部分位于红外线(又称热射线)区段中0.76～20μm的范围内。

所谓红外线加热,就是利用这一区段的热辐射。

研究热辐射规律，对于炉内传热的合理设计十分重要，对于高温炉操作工的劳动保护也有积极意义。

当某系统需要保温时，即使此系统的温度不高，辐射传热的影响也不能忽视。

如保温瓶胆镀银，就是为了减少由辐射传热造成的热损失。

热辐射的基本概念任何物体在发出辐射能的同时，也不断吸收周围物体发来的辐射能。

一物体辐射出的能量与吸收的能量之差，就是它传递出去的净能量。

物体的辐射能力（即单位时间内单位表面向外辐射的能量），随温度的升高增加很快。

一般说来，当一物体受到其他物体投来的辐射（能量为Q）时,其中被吸收转为热能的部分为QA,被反射的部分为QR，透过物体的部分为QD,显然这些部分与总能量之间有下式所示的关系： QA+QR+QD=Q如果把A=QA/Q称为吸收
率,R=QR/Q称为反射率,D=QD/Q称为穿透率，则有： A+R+D=1 若物体的A=1，R=D=0，即到达该物体表面的热辐射的能量完全被吸收，此物体称为绝对黑体,简称黑体。

若R=1，A=D=0,即到达该物体表面的热辐射的能量全部被反射；当这种反射是规则的，此物体称为镜体；如果是乱反射，则称为绝对白体。

若D=1，A=R=0，即到达物体表面的热辐射的能量全部透过物体，此物体称为透热体。

实际上没有绝对黑体和绝对白体，仅有些物体接近绝对黑体或绝对白体。

例如：没有光泽的黑漆表面接近于黑体，其吸收率为0.97～0.98；磨光的铜表面接近于白体，其反射率可达0.97。

影响固体表面的吸收和反射性质的，主要是表面状况和颜色，表面状况的影响往往比颜色更大。

固体和液体一般是不透热的。

热辐射的能量穿过固体或液体的表面后只经过很短的距离（一般小于1mm,穿过金属表面后只经过1μm）,就被完全吸收。

气体对热辐射能几乎没有反射能力，在一般温度下的单原子和对称双原子气体（如 Ar、He、H2、N2、O2等），可视为透热体，多原子气体（如CO2、H2O、SO2、NH3、CH4等）在特定波长范围内具有相当大的吸收能力。

辐射以电磁波和粒子（如阿尔法粒子、贝塔粒子等）的形式向外放散。

无线电波和光波都是电磁波。

它们的传播速度很快，在真空中的传播速度与光波(3×1010厘米/秒)相同，在空气中稍慢一些。

电磁波是由不同波长的波组成的合成波。

它的波长范围从10E-10微米(1微米=10E-4厘米)的宇宙线到波长达几公里的无线电波。

Υ射线、X
射线、紫外线、可见光、红外线，超短波和长波无线电波都属于电磁波的范围。

肉眼看得见的是电磁波中很短的一段，从0.4-0.76微米这部分称为可见光。

可见光经三棱镜分光后，成为一条由红、橙、黄、绿、青、蓝、紫七种颜色组成的光带，这光带称为光谱。

其中红光波长最长，紫光波长最短，其它各色光的波长则依次介于其间。

波长长于红光的(>0.76微米)有红外线有无线电波；波长短于紫色光的(<0.4微米)有紫外线，Υ射线、X射线等。

这些辐射虽然肉眼看不见，但可用仪器测出。

太阳辐射波长主要为0.15-4微米，其中最大辐射波长平均为0.5微米；地面和大气辐射波长主要为3-120微米，其中最大辐射波长平均为
10微米。

习惯上称前者为短波辐射，后者为长波辐射。

！！！！表观反射
率/view/e332f6d53186bceb19e8bbfe.html
图像的DN转化为图像的幅亮度转化为表观反射率！
！！！！！！光度与辐射度基础
/view/d4c2001959eef8c75fbfb3a3.html
如何将TM影像DN值转换为反射率？
L（幅亮度）=gain*DN +bias
gain/bias在头文件中可以读取。

R=（pi*L*d^2）/（Esun*cosA）R为反射率
pi=3.14,?d表示日地距离，A表示太阳天顶角
Esun值(landsat 5)
B1:1957
B2:1826
B3:1554
B4:1036
B5:215
B7:80.67
不知道对不对？另外太阳天顶角A和方位角有何区别？TM头文件中怎么只能获取方位角？日地距离从哪里获取？
日地距离一般如果不是很需要的话，可以选取1。

还有一个公式可以计算。

TM中应该有一个高度角，天顶角＝90－高度角，
你的公式可以用
大气校正有很多种，一般如果简单点的话，就是用你上面的这个公式中的
L*=L - 瑞利散射，再将L*代入你上面的公式就可以了
如果要求严格点的话国际上用6s比较多，不过里面6s中很多参数没法得
我现在也在尝试使用6S模型做,但是参数确实比较多,很头疼.我主要做内陆湖泊水质遥感
上面的公式是建立在晴朗天气下可以用，如果你是专门搞大气校正的话，估计还是要精确的用6s做，我的研究方向不在预处理这块，我主要是用它来做初步的大气校正，利于后期的影像处理
TM数据预处理（DN转表观反射率）
一般我们拿到的TM数据都是灰度值（DN值），必须转换为反射率才能进行运算（比如NDVI运算），否则是不严密的。

由灰度值转换为反射率的过程为：
具体过程：
1、DN转辐射能量值
公式为： =DN*gain +bias
其中L为地物在大气顶部的辐射能量值，单位为；DN为样本的灰度值，gain 和bias分别为图像的增益与偏置，可从图像的头文件中读取(需要经过转换)，头文件一般与原始数据一起提供。

在ENVI中可以这么做：
打开原始影像，用basic tools->preprocessing ->general purpose utilities->apply gain and offset ,
并选中要进行转换的波段，弹出如下对话框：
从头文件中读取该波段的gain和biases值：
也可查看固定值，如下表:
Landsat5和Landsat7各波段光谱通道的增益和偏置单位：W/m2. ster.μm
Table4 the Gains and Biases of L5 and L7 unit: W/m2. ster.μm
《Landsat5图像的增益、偏置取值及其对行星反射率计算分析》）
!PI*(1.01250756^2)*b7/(74.52*(COS(0.436332306)))
!PI*(1.01250756^2)*b5/(219.3*(COS(0.436332306)))）
当然也可以将上述步骤写成一个整体：
（2004年5月21日公式：
!PI*(1.01250756^2)*( b7*0.0568-0.1481)/(74.52*(COS(0.436332306))) !PI*(1.01250756^2)*( b5*0.1168-0.4000)/(219.3*(COS(0.436332306))) !PI*(1.01250756^2)*( b4*0.7418-1.3786)/(1047*(COS(0.436332306))) !PI*(1.01250756^2)*( b3*0.8962-1.3)/(1557*(COS(0.436332306)))
!PI*(1.01250756^2)*( b2*1.1998-2.9)/(1829*(COS(0.436332306)))
!PI*(1.01250756^2)*( b1*0.5658-1.4286)/(1957*(COS(0.436332306)))）
对于ETM数据可以使用ERDAS，提供了很方便的处理接口。

（在不知道增益等参数的情况下，ENVI也有默认值，不过默认值存在较大偏差使用的是1986年的参数，最好用自己数据的参数做，envi中 basic tools->preprocess-> calibration utilities ->landsat TM 3、大气纠正（略）
可用ENVI的FLAASH，或参考<</FONT>基于TM/ETM+遥感数据的地面相对反射率反演
据预处理是遥感应用的第一步，也是非常重要的一步。

目前的技术也非常成熟，大多数的商业化软件都具备这方面的功能。

预处理的大致流程在各个行业中有点差异，而且注重点也各有不同。

本小节包括以下内容：
l数据预处理一般流程介绍
l预处理常见名词解释
lENVI中的数据预处理
1、数据预处理一般流程
数据预处理的过程包括几何精校正、配准、图像镶嵌与裁剪、去云及阴影处
理和光谱归一化几个环节，具体流程图如图所示。

图1数据预处理一般流程
各个行业应用会有所不同，比如在精细农业方面，在大气校正方面要求会高点，因为它需要反演；在测绘方面，对几何校正的精度要求会很高。

2、数据预处理的各个流程介绍
（一）几何精校正与影像配准
引起影像几何变形一般分为两大类：系统性和非系统性。

系统性一般有传感器本身引起的，有规律可循和可预测性，可以用传感器模型来校正；非系统性几何变形是不规律的，它可以是传感器平台本身的高度、姿态等不稳定，也可以是地球曲率及空气折射的变化以及地形的变化等。

在做几何校正前，先要知道几个概念：
地理编码：把图像矫正到一种统一标准的坐标系。

地理参照：借助一组控制点，对一幅图像进行地理坐标的校正。

图像配准：同一区域里一幅图像（基准图像）对另一幅图像校准
影像几何精校正，一般步骤如下，
（1）GCP（地面控制点）的选取
这是几何校正中最重要的一步。

可以从地形图（DRG）为参考进行控制选点，也可以野外GPS测量获得，或者从校正好的影像中获取。

选取得控制点有以下特征：
1、GCP在图像上有明显的、清晰的点位标志，如道路交叉点、河流交叉点等；
2、地面控制点上的地物不随时间而变化。

GCP均匀分布在整幅影像内，且要有一定的数量保证，不同纠正模型对控制点个数的需求不相同。

卫星提供的辅助数据可建立严密的物理模型，该模型只需9
个控制点即可；对于有理多项式模型，一般每景要求不少于30个控制点，困难地区适当增加点位；几何多项式模型将根据地形情况确定，它要求控制点个数多于上述几种模型，通常每景要求在30-50个左右，尤其对于山区应适当增加控制点。

（2）建立几何校正模型
地面点确定之后，要在图像与图像或地图上分别读出各个控制点在图像上的像元坐标（x，y）及其参考图像或地图上的坐标（X，Y），这叫需要选择一个合理的坐标变换函数式（即数据校正模型），然后用公式计算每个地面控制点的均方根误差（RMS）
根据公式计算出每个控制点几何校正的精度，计算出累积的总体均方差误差，也叫残余误差，一般控制在一个像元之内，即RMS<1。

（3）图像重采样
重新定位后的像元在原图像中分布是不均匀的,即输出图像像元点在输入图像
中的行列号不是或不全是正数关系。

因此需要根据输出图像上的各像元在输入图像中的位置，对原始图像按一定规则重新采样，进行亮度值的插值计算，建立新的图像矩阵。

常用的内插方法包括：
1、最邻近法是将最邻近的像元值赋予新像元。

该方法的优点是输出图像仍然保持原来的像元值，简单，处理速度快。

但这种方法最大可产生半个像元的位置偏移，可能造成输出图像中某些地物的不连贯。

2、双线性内插法是使用邻近4个点的像元值，按照其距内插点的距离赋予不同的权重，进行线性内插。

该方法具有平均化的滤波效果，边缘受到平滑作用，而产生一个比较连贯的输出图像，其缺点是破坏了原来的像元值。

3、三次卷积内插法较为复杂，它使用内插点周围的16个像元值，用三次卷积函数进行内插。

这种方法对边缘有所增强，并具有均衡化和清晰化的效果，当它仍然破坏了原来的像元值，且计算量大。

一般认为最邻近法有利于保持原始图像中的灰级，但对图像中的几何结构损坏较大。

后两种方法虽然对像元值有所近似，但也在很大程度上保留图像原有的几何结构，如道路网、水系、地物边界等。

（二）数字图像镶嵌与裁剪
l镶嵌
当研究区超出单幅遥感图像所覆盖的范围时，通常需要将两幅或多幅图像拼接起来形成一幅或一系列覆盖全区的较大的图像。

在进行图像的镶嵌时，需要确定一幅参考影像，参考图像将作为输出镶嵌图像的基准，决定镶嵌图像的对比度匹配、以及输出图像的像元大小和数据类型等。

镶嵌得两幅或多幅图像选择相同或相近的成像时间，使得图像的色调保持一致。

但接边色调相差太大时，可以利用直方图均衡、色彩平滑等使得接边尽量一致，但用于变化信息提取时，相邻影像的色调不允许平滑，避免信息变异。

l裁剪
图像裁剪的目的是将研究之外的区域去除，常用的是按照行政区划边界或自然区划边界进行图像的分幅裁剪。

它的过程可分为两步：矢量栅格化和掩膜计算（Mask）。

矢量栅格化是将面状矢量数据转化成二值栅格图像文件，文件像元大小与被裁剪图像一致；把二值图像中的裁剪区域的值设为1，区域外取0值，与被裁剪图像做交集运算，计算所得图像就是图像裁剪结果。

（三）大气校正
遥感图像在获取过程中，受到如大气吸收与散射、传感器定标、地形等因素的影响，且它们会随时间的不同而有所差异。

因此，在多时相遥感图像中，除了地物的变化会引起图像中辐射值的变化外，不变的地物在不同时相图像中的辐射值也会有差异。

利用多时相遥感图像的光谱信息来检测地物变化状况的动态监测，其重要前提是要消除不变地物的辐射值差异。

辐射校正是消除非地物变化所造成的图像辐射值改变的有效方法，按照校正后的结果可以分为2种，绝对辐射校正方法和相对辐射校正方法。

绝对辐射校正方法是将遥感图像的DN(Digital Number)值转换为真实地表反射率的方法，它需要获取影像过境时的地表测量数据，并考虑地形起伏等因素来校正大气和传感器的影响，因此这类方法一般都很复杂，目前大多数遥感图像都无法满足上述条件。

相对辐射校正是将一图像作为参考（或基准）图像，调整另一图像的DN值，使。