【vb精品源码栏目提醒】:网学会员鉴于大家对vb精品源码十分关注,论文会员在此为大家搜集整理了“【精品】GDAL学习笔记 - 其它资料”一文,供大家参考学习
可能你不玩 GIS,不懂这个库到底有什么用,或者和 python 有什么关系。
但是你要玩GIS,RS,你就应当知道这个库的价值。
就算你不玩 GIS,我想这个库对你也应该有致命的吸引力。
为什么?看下面的介绍吧! 先看看这段 GDAL 主页上的英文介绍吧! is a translator library for rastergeospatial data formats that is released under anX/MITstyleOpen Sourcelicense bytheOpen Source Geospatial Foundation. As a library it presents asingle abstract datamodelto the calling application for all supported formats. It also comes with a variety ofusefulcommandline utilitiesfor data translation and processing. 简单地说,GDAL 是一个操作各种栅格地理数据格式的库。
包括读取、写入、转换、处理各种栅格数据格式(有些特定 。
的格式对一些操作如写入等不支持) 它使用了一个单一的抽象数据模型就支持了大多数的栅格数据(GIS 对栅格,矢量,3D 数据模型的抽象能力实在令人叹服)。
当然除了栅格操作,这个库还同时包括了操作矢量数据的另一个有名的库 ogr(ogr 这个库另外介绍),这样这个库就同时具备了操作栅格和矢量数据的能力,买一送一,这么合算的买卖为什么不做。
最最最重要的是这个库是跨平台的,开源的!如今这个库对各种数据格式的支持强大到令人啧啧的地步了。
如果你对他的强大有什么怀疑的话,看看这里一大串的 GDAL 所支持格式清单,吓到了吧!再看看它的主页最后那些使用了它作为底层数据处理的软件列表吧!其中你可以不知道 GRASS,你也可以不知道 Quantum GIS QGIS,但是你总该知道 Google Earth吧!不知道?赶快下一个去玩玩--会当临绝顶,一览众山小! 有人说我又不玩 GIS。
不错,但是,你即使不玩 GIS,这个库也是满有用的。
首先,哪个库支持这么多栅格(图片)格式,哪个库在 C/C/python/ruby/
VB/java/C这个暂时不完全支持下都能用,而且都一样用?退一步讲,3S 软件又不一定要用在 3S 下(很多医学影像就是用 PCI 软件来处理的)。
再退一步,你的生活即使和 3S 一点关系都没有,栅格数据又不单单只有 GIS 下才用到。
你大可用这个库来读取 jpg,gif,tif,xpm 等格式。
而且对各种格式支持得不是一般的好,很大一部分非标准格式照样支持得非常好。
我曾经在 java下玩过 jai,以及一系列 jai 的扩展库,一些图像格式在很多图片浏览器中都可以正确读取(有的甚至不是非标准格式),用 jai 死活就读不出来! 这个库的 python 版和其他的 python 库结合的很好。
最直接、明显的支持是使用Numeric 库来进行数据读取和操作。
各种矩阵魔术可以发挥得淋漓尽致(图像其实就是矩阵)。
而且按我的观点,python 对矩阵的操作比其他的语言有明显的优势。
写出来的东西比其他语言写出来的短小的多,而且好看得多。
并且 python 的弱类型在处理栅格数据格式类型的时候代码量比强类型的语言少了数倍(不用 double,byte,short 等等分开处理,这简直就是先天上的优势)。
所以我就喜欢用 python 做图像的处理。
所以就连 GIS 界的微软ESRI 也直接在 ARCGIS9 中用 python 来作栅格数据的导入导出。
一句话,真是太方便啦! 2. 安装2.1. windows 下的安装 官方安装文档在这里。
下面是我自己的实践步骤: 先去 http://www.gdal.org/dl/下一个版本,解压。
打开控制台,输入: “D:ProgramFilesMicrosoft Visual Studio .NET 2003Vc7binvcvars32.batquot 注册 vc 的编译环境。
打 开 gdal 文件夹下的 nmake.opt 修改 GDAL_HOME quotC:warmerdabldquot把路径改到需要把 gdal 安装的地方。
不改也可以。
这里需要添加 python 支持,所以修改 PY_INST_DIR GDAL_HOMEpymod 把路径改成 python 下的 Libsite-packages 文件夹下。
PYDIR quotC:SoftwarePython24quot 改成 python 的安装路径。
下面的参数爱改什么就把前面的删除(要看您有没有那些库的
源码),注意一下路径就可以了。
我是都没改。
后面就依次运行Toggle line numbers Toggle line numbers 1 nmake /f makefile.vc 2 nmake /f makefile.vc install 3 nmake /f makefile.vc devinstall 最后最后,还要去 GDAL_HOME 目录下的 bin 文件夹下把 gdal13.dll (也有可能是gdal12.dll)copy 到 PY_INST_DIR 路径下 到此处就完成安装 gdal(python)的工作。
最后需要注意一下, 1.3 gdal 在.net2005 下只能顺利编译 1.2, 以上版本不能顺利编译,有一个地方指针转换出错。
可能是 2005 的编译器比以往的严厉一点吧。
另外安装了 QGIS对编译也有一些影响主要是 proj 库的冲突导致一个找不到quotd:/program.objquot文件的错误如果你有静态编译过 proj那么你可以打开 nmake.opt 修改有关 proj 的设置如果搞不定就卸载 QGIS然后编译编译后再安装 QGIS.呵呵还好 QGIS 的体积没有 ArcGIS 那么可怕.2.2. linux 下的安装 linux 下的安装就更简单了。
直接Toggle line numbers Toggle line numbers 1 ./configure 2 make 3 su 4 make install 5 ldconfig 就 ok默认就已经支持 python。
当然在第一步的时候需要看看是否依赖的库都安装了。
如果缺少,就去安装一个。
如果对 configure 的条件不理解,就用./configure --help 看看具体情况。
2.3. 安装其他驱动 这里讲一个安装 hdf4 的驱动的例子(默认情况下 gdal 是不安装 hdf4 的),其他驱动应该和这个也差不了多少吧,可以作为其他的参考。
完整步骤如下: 在 windows 下的安装: 从 ftp://ftp.ncsa.uiuc.edu/HDF/HDF/HDF5FCurrent/bin/windows/下载42r1-win.ZIP,解压。
编辑 gdal 根目录下的 nmake.opt,找到“ Uncomment the following and update toenable NCSA HDF Release 4 support.”这一行 把下面两行前面的去掉,然后改成: HDF4_DIR D:warmerda42r1-winrelease HDF4_LIB /LIBPATH:HDF4_DIRlib hd421m.lib HDF4_LIB HDF4_DIRdllhd421m.lib HDF4_DIRdllhm421m.lib HDF4_DIRlibhd421.lib HDF4_DIRlibhm421.lib 用 HDF4_LIB/LIBPATH:这种形式似乎可以建立 gdal 的库,但是往下编译会出错。
而且要把HDF4_DIRdll 和HDF4_DIRlib 拷贝到同一个目录下,不然会提示找不到库 你也可以试一试在 D:Program FilesMicrosoft Visual Studio .NET2003Common7Toolsvsvars32.bat 文件中添加 HDF4_LIB 路径到“set LIB”这行的末尾(不要忘记的分割符)。
然后找一下quotINCquot这行把 -IHDF4_DIRinclude 加到下一行的末尾(应该也可以在vsvars32.bat 中添加路径,不过要重启命令行)。
然后编译吧!祝你好运。
注意:上面的 HDF4_DIR 是我本机的路径,你要根据你自己的路径进行设置(想起我的一个老师说过的话:“抄人家的作业可以,不要连名字也一起抄走啊 ” ),下面的HDF4_DIR可以不用改,那个是变量,会自动替代 HDF4_DIR 路径。
编译成功后,要 HDF4 能运行,还需要两个库,一个是 zlib一个是 szip,可以到下面两个链接去下载一个 ftp://hdf.ncsa.uiuc.edu/lib-external/zlib/1.2/bin ftp://hdf.ncsa.uiuc.edu/lib-external/szip/2.0/bin/windows 把这两个库下载后解压,然后设置 PATH 系统变量,使得它们在默认状态下也可以被动态链接成功 。
在 Linux 下比在 Windows 下简单: 只要用./configure --help 察看一下打开 HDF4 的编译开关(包括库路径,头文件路径等,看清楚),然后在./configure 后面加上那个开关以及 hdf4 的安装路径后就可以了。
在configure 后 gdal 会提示是否支持 HDF4。
编译后也要把 zlib 和 szip 的动态链接库设置好 。
到此你已经可以用 C/C来操作 gdal 读写 hdf4 的格式了! 最后,为了让 Python 能够支持 hdf 的读写,别忘了把重新生成安装后的 pymod 目录下的内容,还有 gdal13,还有那两个 hdf 的库,还有 zlib,szip 的库拷贝到 Python 的Libsite-packages 目录下 。
2.4. 下载 如果你实在玩不转,可以在这里下载已经编译好的 gdal1.3.2 程序库 以及其依赖的其他库,其中包括 hdf4,hdf5 支持,以及 proj,geos 插件。
注意,这里的 geos 是静态链接的,注意版权(geos 是 LGPL 的 license)。
hdf4 和 hdf5 用的是 release 版本。
这里是我的nmake 配置文件,你可以对照你的实际情况参考一下。
3. 快速开始 其实在主站的教程里已经有 python 的示例了。
但是我们还是按照自己的思路来开始吧。
第一步就是打开一个数据集。
对于“数据集”这个名词大家可能不会太习惯,但是对于一般的格式来说,一个“数据集”就是一个文件,比如一个 gif 文件就是一个以 gif 为扩展名的文件。
但是对于众多 RS 数据来说,一个数据集包含的绝对不仅仅是一个文件。
对于很多 RS数据,他们把一张图像分成数个图像文件,然后放在一个文件夹中,用一些额外的文件来组织它们之间的关系,形成一个“数据集”。
如果你不理解,那么就算了,当成 jpg 或者 gif 文件好了。
下面我们打开一个 tiff 文件GeoTIFF。
这个文件是我从 GRASS 的示例数据 spearfish中导出的一个同名影像数据。
Toggle line numbers Toggle line numbers 1 gtgtgt import gdal 2 gtgtgt dataset gdal.Openquotj:/gisdata/gtif/spot.tifquot 3 gtgtgt dirdataset 4 AddBand AdviseRead BuildOverviews FlushCacheGetDescription Get 5 Driver GetGCPCount GetGCPProjection GetGCPsGetGeoTransform GetMe 6 tadata GetProjection GetProjectionRef GetRasterBandGetSubDatasets 7 RasterCount RasterXSize RasterYSize ReadAsArrayReadRaster Refr 8 eshBandInfo SetDescription SetGCPs SetGeoTransformSetMetadata Se 9 tProjection WriteRaster __del__ __doc__ __init____module__ _b 10 and _o 11 gtgtgt 这样我们就打开了这个文件。
并且我们可以看到可以供我们调用的函数们更具体的API 列表可以看这里。
现在我们不做修改,不做添加,所以只要带有 Set 开头的函数以及有 Write 开头的函数我们暂时都不管。
因为 RS 影像必然要和地理上的位置挂上钩,才能把图像正确铺展到一个坐标系中。
其中的信息和对应关系有点复杂,不适合在快速开始中介绍,我们暂时也先不管。
这里需要注意的就是几个函数。
GetDescription 获得栅格的描述信息。
Toggle line numbers Toggle line numbers 1 gtgtgt dataset.GetDescription 2 j:/gisdata/gtif/spot.tif 3 gtgtgt 看来这里的图像描述是图像的路径名,但是这是和各种不同数据集相关的,不同数据集可能有不同的描述。
这要看读取驱动的实现作者的高兴了。
RasterCount 获得栅格数据集的波段数。
GetRasterBand 获得栅格数据集的波段。
Toggle line numbers Toggle line numbers 1 gtgtgt dataset.RasterCount 2 1 3 gtgtgt band dataset.GetRasterBand1 4 gtgtgt 这里需要解释的是 Band 这个词。
这个词可以翻译成“波段”,“通道”等等。
我这里把它统一称为“波段”。
因为遥感卫星的传感器有很多个。
一个传感器只负责接收一个频率范围的地物反射光波,一个频率范围的光波记录称为一个波段。
是不是晕了?其实说得简单一点。
其实你可以把波段看成红绿蓝几种颜色。
图像不是分 RGB 三色吗?把 R,G,B 值都提取出来成为三个表。
R 值表就是波段一,G 值表就是波段二,B 值表就是波段三。
这里我们看到这张图只有一个波段一种颜色。
就可以把它看成是一个灰度图类似黑白照片。
如果 RasterCount 是 3,就有可能是彩色图。
如果 RasterCount 是比 3 大的数,恭喜你,你看到一张遥感影像。
有很多卫星的传感器大于 3 个,比如 TM 就有 7 个波段,不仅有可见光,还有红外等其他非可见光。
,所以,波段一般比 RGB 能表达的丰富地多。
不过这样一来就需要我们从中挑出 3 个波段然后组合成 RGB,当然这样就有可能使图像显示出来的东西不像平常我们看到的那样。
这样安排是因为对科学有帮助一些波段在科学家眼里比真实的彩色照片有价值。
不理解就跳过,很正常,我第一次听这种东西也觉得很玄: 这里我们获取了第一个波段红色值组成的表。
注意!这里的波段获取和通常的 C 数组获取不一样,开始是 1 不是 0。
获取了波段,我们就可以在下面的操作中读取这个波段的所有数值。
RasterXSize 图像的宽度X 方向上的像素个数 RasterYSize 图像的高度Y 方向上的像素个数Toggle line numbers Toggle line numbers 1 gtgtgt dataset.RasterXSize 2 950 3 gtgtgt dataset.RasterYSize 4 700 5 gtgtgt 可以看出我们的图像大小是 950700。
还是很小的一张图。
ReadRaster 读取图像数据以二进制的形式 ReadAsArray 读取图像数据以数组的形式Toggle line numbers Toggle line numbers 1 gtgtgt helpdataset.ReadRaster 2 Help on method ReadRaster in module gdal: 3 ReadRasterself xoff yoff xsize ysize buf_xsizeNonebuf_ysizeNone buf_t 4 ypeNone band_listNone method of gdal.Dataset instance 5 gtgtgt helpdataset.ReadAsArray 6 Help on method ReadAsArray in module gdal: 7 ReadAsArrayself xoff0 yoff0 xsizeNone ysizeNone method ofgdal.Dataset 8 instance 9 gtgtgt 这两个函数很重要,它们直接读取图像的数据,可以看到两个函数的帮助中有一大溜的参数。
解释一下: xoffyoffxsizeysize 你可能不想读取整张图像。
只想读取其中的一部分。
那么就用 xoff,yoff 指定想要读取的部分原点位置在整张图像中距离全图原点的位置。
用 xsize 和 ysize 指定要读取部分图像的矩形大小。
buf_xsize buf_ysize 你可以在读取出一部分图像后进行缩放。
那么就用这两个参数来定义缩放后图像最终的宽和高,gdal 将帮你缩放到这个大小。
buf_type 如果你要读取的图像的数据类型不是你想要的比如原图数据类型是 short,你要把它们缩小成 byte,就可以设置它。
band_list 这就适应上面多波段的情况。
你可以指定读取的波段序列。
要哪几个波段,不要哪几个波段,你说了算。
举个例子吧:Toggle line numbers Toggle line numbers 1 gtgtgt dataset.ReadAsArray2302701010 2 array255 255 255 232 232 255 255 255 255 222 3 255 255 255 255 255 255 210 110 11 122 4 255 255 255 255 255 255 210 255 11 243 5 201 255 255 255 255 200 200 110 122 243 6 111 211 255 201 255 255 100 11 132 243 7 255 100 100 100 110 100 110 111 122 243 8 255 255 255 255 255 255 122 222 255 255 9 255 255 255 255 255 255 243 243 255 255 10 255 255 255 255 255 255 255 255 255 255 11 255 255 255 255 255 255 255 255 255 255b 12 gtgtgt dataset.ReadRaster2302701010 13xffxffxffxe8xe8xffxffxffxffxdexffxffxffxffxffxffxd2nx0bzxff 14xffxffxffxffxffxd2xffx0bxf3xc9xffxffxffxffxc8xc8nzxf3oxd3xffx 15c9xffxffdx0bx84xf3xffdddndnozxf3xffxffxffxffxffxffzxdexffxffxff 16xffxffxffxffxffxf3xf3xffxffxffxffxffxffxffxffxffxffxffxffxff 17 xffxffxffxffxffxffxffxffxff 18 gtgtgt 我们就把图像中位于 230270,宽度 10 高度 10 的数据读取出来了。
我们看完了数据集的主要函数。
似乎已经够用了。
的确,如果只是为了显示图像,这些的确已经够了。
但是如果需要更多信息,我们就不得不进入波段操作数据实际上我们大多数时候都需要进入 band 获取信息。
下面我们现在来看看刚才读取出来的那个 band 有些什么东西可以供我们操作的具体的 API 列表看这里。
Toggle line numbers Toggle line numbers 1 gtgtgt dirband 2 AdviseRead Checksum ComputeBandStatsComputeRasterMinMax DataType 3 Fill FlushCache GetDefaultHistogram GetDescriptionGetHistogram 4 GetMaximum GetMetadata GetMinimum GetNoDataValueGetOffset GetO 5 verview GetOverviewCount GetRasterColorInterpretationGetRasterColorTab 6 le GetScale GetStatistics ReadAsArray ReadRasterSetDefaultHistog 7 ram SetDescription SetMetadata SetNoDataValueSetRasterColorInterpre 8 tation SetRasterColorTable WriteArray WriteRasterXSize YSize 9 __doc__ __init__ __module__ _o 10 gtgtgt 挑几个有用的吧。
Toggle line numbers Toggle line numbers 1 gtgtgt band.XSize 2 950 3 gtgtgt band.YSize 4 700 5 gtgtgt band.DataType 6 1 7 gtgtgt 不用解释了吧,波段图像的宽和高象元为单位。
DataType,图像中实际数值的数据类型。
具体数据类型定义在 gdalconst 模块里。
使用的时候用 import gdalconst 引入。
Toggle line numbers Toggle line numbers 1 gtgtgt import gdalconst 2 gtgtgt dirgdalconst 3 CE_Debug CE_Failure CE_Fatal CE_None CE_WarningCPLES_Backslash 4 Quotable CPLES_CSV CPLES_SQL CPLES_URL CPLES_XMLCPLE_AppDefined 5 CPLE_AssertionFailed CPLE_FileIO CPLE_IllegalArgCPLE_NoWriteAccess 6 CPLE_None CPLE_NotSupported CPLE_OpenFailedCPLE_OutOfMemory CPLE 7 _UserInterrupt CXT_Attribute CXT_Comment CXT_ElementCXT_Literal 8 CXT_Text DCAP_CREATE DCAP_CREATECOPYDMD_CREATIONDATATYPES DMD_CREAT 9 IONOPTIONLIST DMD_EXTENSION DMD_HELPTOPIC DMD_LONGNAMEDMD_MIMETYPE 10 GA_ReadOnly GA_Update GCI_AlphaBand GCI_BlackBandGCI_BlueBand 11 GCI_CyanBand GCI_GrayIndex GCI_GreenBand GCI_HueBandGCI_LightnessB 12 and GCI_MagentaBand GCI_PaletteIndex GCI_RedBandGCI_SaturationBand 13 GCI_Undefined GCI_YellowBand GDT_Byte GDT_CFloat32GDT_CFloat64 14 GDT_CInt16 GDT_CInt32 GDT_Float32 GDT_Float64GDT_Int16 GDT_Int 15 32 GDT_TypeCount GDT_UInt16 GDT_UInt32 GDT_UnknownGF_Read GF_ 16 Write GPI_CMYK GPI_Gray GPI_HLS GPI_RGBGRA_Bilinear GRA_Cubic 17 GRA_CubicSpline GRA_NearestNeighbour __builtins____doc__ __file 18 __ __name__ 19 gtgtgt 那些 GDT 开头的就是数值数据类型。
Toggle line numbers Toggle line numbers 1 gtgtgt band.GetNoDataValue 2 65535.0 3 gtgtgt band.GetMaximum 4 gtgtgt band.GetMinimum 5 gtgtgt band.ComputeRasterMinMax 6 1.0 255.0 7 gtgtgtMaximum 是表示在本波段数值中最大的值,Minimum 当然就是表示本波段中最小的值啦。
我们可以看到在一开始这两个都没有值。
因为对于文件格式不会有固有的最大最小值。
所以我们通过函数ComputeRasterMinMax 计算得到了。
注意!这里的最大最小值不包括“无意义值”!也就是上面显示的NoDataValue。
需要解释一下“无意义值”。
不要以为 0 或者 255 在任何情况下都无意义。
在很多情况下 0,255 需要和其他值一样表示一个实际意义。
虽然可能它最终会被显示得和黑色一样。
而一些位置上的点要表示的意思是“什么也不是”,它在那个位置上只是为了占一个位置,.
上一篇:
【精品】图书管理系统《软件测试技术》课程设计报告——测试计划
下一篇:
对普通机床的认识