遥感原理与应用算法论文.pdf

基于GDAL 库的遥感图像处理相关算法实现 1201 班，2012302590190，包文韬 （武汉大学遥感信息工程学院） 摘 要：本文介绍了在VC++6.0 集成开发环境下，利用GDAL 开发库实现几个遥感影像处理 算法的MFC 程序设计。具体内容包括GDAL 开发库的使用、图像的RGB 与HIS 相互转换、图 像的灰度化处理、多光谱影像K 均值聚类算法以及基于HIS 变换的遥感影像融合算法。 关键词：GDAL；MFC；HIS 变换；K 均值聚类；图像灰度化；遥感影像融合 1 原理介绍 1.1 GDAL 库简介 GDAL(Geospatial Data Abstraction Library)是一个在X/MIT 许可协议下 的开源栅格空间数据转换库。提供对多种栅格数据的支持，包括Arc/Info ASCII Grid(asc)，GeoTiff (tiff)，Erdas Imagine Images(img)，ASCII DEM(dem) 等 格式。GDAL 使用抽象数据模型(abstract data model)来解析它所支持的数据格 式，其抽象数据模型包括数据集(dataset)，坐标系统，仿射地理坐标转换(Affine Geo Transform)， 大地控制点(GCPs)， 元数据(Metadata)，栅格波段(Raster Band)，颜色表(Color Table)，子数据集域(Subdatasets Domain)，图像结构域 (Image_Structure Domain)，XML 域(XML:Domains)。常用的 GDAL 类包括 GDALDataset 类、GDALRasterBand 类、GDALDriver 类、GDALDriverManager 类、 GDALMajorObject 类等。核心类结构如下： GDALMaj orObj ect GDALDataset GDALDriver GDALDriverManager GDALRasterBand 图1.1 GDAL 核心类结构 1.2 RGB 与HIS 正逆变换 在计算机上定量处理色彩时通常采用RGB 表色系统，但在视觉上定性描述色 彩时，采用HIS 系统则更为直观，IHS （Intensity，Hue，Saturation）分别表 示亮度或强度（I）、色调（H）和饱和度（S），它们是从人眼中认识颜色的三个 特征。HIS 定义的假彩色空间如图1.2.1。明度I 沿着轴线从底部的黑变到顶部 的白。色度H 由圆柱底面圆的角度表示。假定0°为红色，120°为绿色，240° 为蓝色，则色度0°到360°覆盖了所有可见光谱的彩色。饱和度S 是色度环的 原点（圆心）到彩色点的半径的长度。圆心的饱和度为零，圆周上的饱和度为I。 图1.2.1 RGB 颜色空间 图1.2.2 IHS 颜色空间 为了实现RGB 颜色空间域HIS 颜色空间的相互转换，需要借助一定的转换模 型。常用的转换模型有球体彩色变换、圆柱体彩色变换、三角形彩色变换以及单 六角锥彩色变换4 种。根据前人的大量试验表明，球体变换和三角形变换是比较 好的转换模型，因此，本文主要实现这两种转换模型。 图1.2.2 球体变换模型与三角形变换模型 1.2.1 IHS 正变换 M  R M  G M  B r ,g ,b 令M = max(R,G,B),m = min(R,G,B)， M  m M  m M  m ，其中r,g,b 1 1/ 2[(R  G)  (R  B )]  cos 2 1/ 2