GF-1数据前期处理步骤:以16m分辨率,锡盟地图图像为例一、数据读取在ENVI中,选择file-OpenAs-CRESDA-GF-1,直接选择.xml文件打开影像数据。补丁更新了ENVI可以识别的传感器列表,加入了国产卫星。如下图所示:图1新增对国产卫星传感器的支持在Datamanager中可以看到ENVI自动识别了相应的RPC文件。如下图所示。图2DataManager中显示的RPC信息在LayerManager中选中影像数据,单击右键,选择ViewMetadata,可以查看元数据信息,包括传感器类型、RPC信息、坐标系统、光谱信息、获取时间等,都可以自动读取。如下图所示,是多光谱图像各个波段的详细信息。图3Metadata查看二、辐射定标1.在Toolbox中,选择RadiometricCorrection/RadiometricCalibration,打开辐射定标工作流,选择全色影像作为InputFile,点击Ok。2.进入RadiometricCalibration对话框。CalibrationType(定标类型)有两个下拉选项:Radiance(辐射亮度)和Reflectance(反射率),这里选择Radiance;单击ApplyFLAASHSetting,因为下面我们将进行Flaash大气校正,所以这里应用为Flaash大气校正对图像文件的要求:OutputInterleave(输出数据存储类型)为BIL格式,OutputDataType(输出数据类型)为Float浮点型,ScaleFactor(缩放尺度)为0.10;最后选择输出路径和文件名,点击OK。完成定标。图9辐射定标参数设置3.如下图所示,得到辐射亮度数数据。图10辐射亮度结果三、Flaash大气校正1.在Toolbox中,选择RadiometricCorrection/AtmosphericCorrectionModule/FLAASHAtmosphericCorrection,打开FLAASH大气校正的工作流。2.在FLAASH大气校正参数面板,有三大部分参数,首先设置第一部分:选择输入文件为经过辐射定标后的多光谱文件,设置如下图:选择输出路径和文件名。3.FLAASH大气校正参数面板第二部分,传感器基本信息设置:SceneCenterLocation(成像中心点经纬度)自动从影像中获取;SensorAltitude(传感器高度):在SensorType(传感器类型)中选择Multispectral/GF-1,可以自动获取,PixelSize(像元大小):16m;GroundElevation(成像区域平均高度):0.9(锡盟为0.9)FlightDate(成像时间):可以在数据上点击右键,选择ViewMetadata,从Time中读取为2013-06-13T10:35:07Z,实际上是2013年6月13日2点35分07秒,减去8换算为GMT时间。4.FLAASH大气校正参数面板第三部分:大气模型根据影像时间和地区选择,气溶胶类型选择Rural,AerosolRetrieval选择None。整个FLAASH大气校正参数设置如下图所示。图12FLAASH大气校正参数5.单击MultispectralSetting按钮,可以看到FilterFunctionFile已经自动读取入光谱响应曲线“gf1.sli”。6.单击AdvancedSettings,在高级设置中:TileSize默认的是Cashsize的大小,根据自己机器内存设定,可以设定为内存的75%,8G内存设定6000M。讨论:一般中低分辨率的影像近似天顶角(ZenithAngle)为180,方位角(AzimuthAngle)为0,即垂直观测。ENVIFLAASH的天顶角:90-180度,方位角:-180and180,没有查到高分一号的角度说明,这里选择默认。图13大气校正高级参数设置7.设置好后,在大气校正模块面板中,单击Apply。8.大气校正完成后,检查大气校正的结果,分别加载校正前后的图像(选择CIR假彩色方式加载,可以更好的识别植被),查看典型地物的大气校正前后的光谱曲线,大致可以看出大气校正之后消除了大气散射的影响。如下图所示。图14高分一号FAASH大气校正结果图15大气校正前后的植被波谱曲线(上-校正前,下-校正后)四、正射校正ENVI自动识别了RPC文件,直接选择ENVI中的正射校正的workflow进行正射校正。由于缺少控制点信息这里直接使用无控制点正射校正功能。1.在Toolbox中,选择GeometricCorrection/Orthorectification/RPCOrthorectificationWorkflow,打开正射校正工作流。选择待校正影像文件。点击Next。2.进入Orthorectification参数设置。Outputpixelsize像元大小为16m。如下图所示。图5正射校正参数3.选择输出目录