1. GDAL附帶工具

1.1. GDAL創(chuàng)建了下面的工具程序

• gdalinfo - 統(tǒng)計(jì)文件的信息。

• gdal_translate - 拷貝一個(gè)可以控制控制輸出格式的柵格文件。

• gdaladdo - 為文件添加略縮圖。

• gdalwarp - 把一個(gè)圖像轉(zhuǎn)換到一個(gè)新的坐標(biāo)系統(tǒng)。

• gdalindex - 為MapServer建立一個(gè)柵格索引。

• gdal_contour - 給DEM創(chuàng)建等高線。

• rgb2pct.py - 把一個(gè)24位的圖像轉(zhuǎn)化成8位調(diào)色板圖像。

• pct2rgb.py - 把一個(gè)8位調(diào)色板圖像轉(zhuǎn)化成24位位圖。

• gdal_merge.py - 把一系列的圖像創(chuàng)建一個(gè)快速的地圖集。

• gdal_rasterize - 把一個(gè)矢量圖層?xùn)鸥窕?/P>

• nearblack - 轉(zhuǎn)換臨近的黑/白邊框?yàn)榫珳?zhǔn)數(shù)值。

• gdal-config - 獲取建立使用gdal的軟件需要的參數(shù)。

2. 創(chuàng)建新的文件

存取一個(gè)已存在的文件來(lái)讀取是一件很容易的事情，只要在命令行中指定文件或者數(shù)據(jù)集的名字。但是，創(chuàng)建一個(gè)文件是一件非常復(fù)雜的事情。你可能需要指定創(chuàng)建格式，各種創(chuàng)建參數(shù)，以及指定一個(gè)坐標(biāo)系統(tǒng)。在不同的GDAL工具中有許多參數(shù)都是差不多的，這里列舉出來(lái)。

-of format

• 選擇要?jiǎng)?chuàng)建新的文件的格式。這個(gè)格式被指定為類似GTiff（GeoTIFF格式）或者HFA（ERDAS格式）。所有的格式列表可以用--formats參數(shù)列出來(lái)。只有格式列表“(rw)”可以被寫入。許多工具如果沒有指定，默認(rèn)是創(chuàng)建GeoTIFF文件。文件擴(kuò)展名不能用來(lái)猜測(cè)輸出格式，如果沒有指定文件名，gdal一般不會(huì)添加任何擴(kuò)展名。

-co NAME=VALUE

• 許多格式會(huì)有一個(gè)或者更多的創(chuàng)建參數(shù)來(lái)控制文件創(chuàng)建的細(xì)節(jié)。比如GeoTIFF可以用創(chuàng)建參數(shù)控制壓縮，或者控制是否用分片還是分帶來(lái)進(jìn)行存儲(chǔ)。

  可以使用的創(chuàng)建參數(shù)根據(jù)格式驅(qū)動(dòng)不同而不同。而一些簡(jiǎn)單的格式根本就沒有創(chuàng)建參數(shù)。雖然某個(gè)格式可以用"--format <format>"參數(shù)列出所有可用的參數(shù)列表，但是更詳細(xì)的信息可以在格式介紹網(wǎng)頁(yè)中查到。

-a_srs SRS

• 有幾個(gè)工具(gdal_translate,gdalwarp)可以在命令行中通過類似-a_srs(分配輸出SRS)、-s_srs(源SRS)、-t_srs(目標(biāo)SRS)來(lái)指定坐標(biāo)系統(tǒng)這些工具允許以一系列格式定義坐標(biāo)系統(tǒng)（SRS就是空間參考系統(tǒng)）。

• NAD27/NAD83/WGS84/WGS72:這些常見的地理坐標(biāo)系統(tǒng)可以通過名字來(lái)直接使用。

• EPSG:n坐標(biāo)系統(tǒng)（投影或者地理坐標(biāo)）可以通過EPSG碼來(lái)選擇。例如EPSG 27700是英國(guó)國(guó)家網(wǎng)格。一系列的EPSG坐標(biāo)系統(tǒng)可以在GDAL數(shù)據(jù)文件gcs.csv和pcs.csv中找到。

• PROJ.4定義： 一個(gè)PROJ4定義字符串可以用作坐標(biāo)系統(tǒng)定義。例如“+proj=utm +zone=11 +datum=WGS84”。注意在命令行中要保持Proj4字符串在一起作為一個(gè)單獨(dú)的參數(shù)（一般用雙引號(hào)引起來(lái)）。

• OpenGIS 知名文本: OpenGIS標(biāo)準(zhǔn)定義了一個(gè)文本格式來(lái)描述坐標(biāo)系統(tǒng)作為簡(jiǎn)單要素規(guī)范的一個(gè)部分。這個(gè)格式是gdal中使用的坐標(biāo)系統(tǒng)的內(nèi)部工作格式。包含wkt坐標(biāo)系統(tǒng)描述的文件的文件名可以被用來(lái)作為坐標(biāo)系統(tǒng)參數(shù)，或者坐標(biāo)系統(tǒng)元素本身也可以被用來(lái)作為命令行參數(shù)（雖然wkt中釋放所有引號(hào)是很有爭(zhēng)議的）。

• ESRI知名文本：ESRI在他們的ArcGIS產(chǎn)品（ArcGIS .prj文件）中使用了一種經(jīng)過精簡(jiǎn)OGC WKT的格式，而且這個(gè)格式被用在一個(gè)和wkt相似的風(fēng)格的文件中。但是文件名要被加以ESRI::前綴。比如"ESRI::NAD 1927 StatePlane Wyoming West FIPS 4904.prj"。

3. 總的命令行參數(shù)

所有的GDAL命令行工具程序都支持下面的“總的”命令行參數(shù)。

--version

• 登記gdal版本并退出。

--formats

• 列出所有g(shù)dal支持的柵格格式（只讀和讀寫的）并退出。格式支持描述如下："ro"是只讀驅(qū)動(dòng)，'rw'是讀寫驅(qū)動(dòng)（比如支持CreateCopy）；"rw+"是讀寫和更新驅(qū)動(dòng)（比如支持Create）。

--format format

• 列出單個(gè)格式驅(qū)動(dòng)的細(xì)節(jié)信息。格式名需要是在--formats 中列出的格式名。比如GTiff。

--optfile file

• 讀取指定名稱的文件并把其中的內(nèi)容當(dāng)成參數(shù)傳入命令行列表。行首以#開頭的行將被忽略。多字組成的參數(shù)需要用雙引號(hào)來(lái)保正其為單一的參數(shù)。

--config key value

• 設(shè)置配置把指定鍵設(shè)置為某個(gè)值，而不必把他們?cè)O(shè)置為環(huán)境變量。一些命令參數(shù)鍵是GDAL_CACHEMAX(用于緩存的內(nèi)存有多少M(fèi))以及GDAL_DATA（gdal的數(shù)據(jù)路徑）。單一的驅(qū)動(dòng)會(huì)被其他配置參數(shù)影響。

--debug value

• 控制調(diào)試信息的打印。ON值表示允許調(diào)試信息輸出，OFF值表示不要輸出調(diào)試信息。只有在調(diào)試信息包含了調(diào)試前綴編碼中的字符串后其他值才會(huì)有選擇。

--help-general

• 給出一個(gè)簡(jiǎn)短的普通的GDAL命令行參數(shù)用法信息然后退出。

4. 各工具程序詳解

4.1. gdalinfo

列出柵格數(shù)據(jù)集的信息

用法：

gdalinfo [--help-general] [-mm] [-stats] [-nogcp] [-nomd] [-mdd domain] * datasetname

gdalinfo程序列出gdal支持的柵格格式的一系列信息。

-mm

• 強(qiáng)制計(jì)算柵格每個(gè)波段的最大最小值

-stats

• 牽制計(jì)算柵格的統(tǒng)計(jì)值

-nogcp

• 抑制地面控制點(diǎn)列表打印。這可能對(duì)大量的GCP的數(shù)據(jù)集們來(lái)說(shuō)是十分有用的。比如L1B AVHRR或者h(yuǎn)df4 MODIS數(shù)據(jù)，這些數(shù)據(jù)包含了成千上萬(wàn)的地面控制點(diǎn)。

-nomd

• 抑制原數(shù)據(jù)打印，一些數(shù)據(jù)集可能包含極多的原數(shù)據(jù)字符串。

-mdd domain

• 給出指定域的原數(shù)據(jù)信息。

gdalinfo會(huì)給出如下的信息（如果有的話）：

• 用于存取的文件格式驅(qū)動(dòng)
• 柵格大小(行列數(shù))
• 文件的坐標(biāo)系統(tǒng)(OGC WKT形式)
• 圖像關(guān)聯(lián)到地理的轉(zhuǎn)換參數(shù)(當(dāng)前不包含旋轉(zhuǎn)系數(shù))
• 地理上的邊界坐標(biāo)，如果可能的話還有基于經(jīng)緯度的完整的地理轉(zhuǎn)換參數(shù)（如果是GCPs就沒有）
• 地面控制點(diǎn)(GCPs)
• 廣義的（包括子?xùn)鸥竦脑瓟?shù)據(jù)）文件原數(shù)據(jù)
• 波段數(shù)據(jù)類型
• 波段光度解析
• 波段瓦片大小
• 波段描述
• 波段最大最小值（已經(jīng)經(jīng)過計(jì)算的）
• 波段無(wú)意義值
• 波段可獲得的略縮圖分辨率
• 波段單位類型(如：波段的高程是米制還是英制)
• 波段的假顏色列表

例子:

gdalinfo ~/openev/utm.tif Driver: GTiff/GeoTIFF Size is 512, 512 Coordinate System is: PROJCS["NAD27 / UTM zone 11N", GEOGCS["NAD27", DATUM["North_American_Datum_1927", SPHEROID["Clarke 1866",6378206.4,294.978698213901]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-117], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["metre",1]] Origin = (440720.000000,3751320.000000) Pixel Size = (60.000000,-60.000000) Corner Coordinates: Upper Left ( 440720.000, 3751320.000) (117d38'28.21"W, 33d54'8.47"N) Lower Left ( 440720.000, 3720600.000) (117d38'20.79"W, 33d37'31.04"N) Upper Right ( 471440.000, 3751320.000) (117d18'32.07"W, 33d54'13.08"N) Lower Right ( 471440.000, 3720600.000) (117d18'28.50"W, 33d37'35.61"N) Center ( 456080.000, 3735960.000) (117d28'27.39"W, 33d45'52.46"N) Band 1 Block=512x16 Type=Byte, ColorInterp=Gray

4.2. gdal_translate

不同格式間的柵格數(shù)據(jù)轉(zhuǎn)換

用法：

gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-not_strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing]* [-mo "META-TAG=VALUE"]* [-quiet] [-sds] [-co "NAME=VALUE"]* src_dataset dst_dataset

gdal_translate 工具可以用來(lái)在不同格式間轉(zhuǎn)換柵格數(shù)據(jù)。并且在處理過程中暗地里做一些諸如子?xùn)鸥裨O(shè)置，重采樣和象元值集體變化等勾當(dāng)。

-ot type

• 指定輸出波段的數(shù)據(jù)類型

-not_strict

• 當(dāng)轉(zhuǎn)換到輸出格式時(shí)遇到不匹配或者丟失數(shù)據(jù)時(shí)，寬大處理。

-of format

• 選擇輸出格式。默認(rèn)是GeoTiff（GTiff）。注意，指定的時(shí)候用格式簡(jiǎn)稱。

-b band

• 選擇一個(gè)波段來(lái)輸出。波段號(hào)從1開始編號(hào)。多個(gè) -b 參數(shù)可以用于選擇輸出某幾個(gè)波段或者重新進(jìn)行波段的排序。

-outsize xsize[%] ysize[%]

• 設(shè)置輸出文件的大小。大小以象元為單位，除非用“%”來(lái)標(biāo)記，這時(shí)，表示的是輸出為輸入圖層大小的百分比。

-scale [src_min src_max [dst_min dst_max] ]

• 重新組織輸入象元的值。將它們從src_min~src_max范圍內(nèi)縮放到dst_min ~ dst_max范圍內(nèi)。如果省略，輸出范圍將為0~255。輸入范圍將由源數(shù)據(jù)自動(dòng)計(jì)算。

-srcwin xoff yoff xsize ysize

• 選擇一個(gè)取值窗口，通過該窗口在原圖像中的行列位置來(lái)拷貝數(shù)值。

-projwin ulx uly lrx lry

• 選擇一個(gè)取值窗口，通過該窗口在原圖像中地理坐標(biāo)范圍來(lái)拷貝數(shù)據(jù)（類似srcwin）。

-a_srs srs_def

• 把輸出文件投影強(qiáng)制指定坐標(biāo)系。srs_def可以是任何常用的GDAL/OGR格式，WKT，Proj4，EPSG：n 或者一個(gè)包含著wkt的文件的文件名。

-a_ullr ulx uly lrx lry

• 強(qiáng)制指定輸出文件的空間轉(zhuǎn)換邊界范圍。而不需要和原圖像有關(guān)系。

-a_nodata value

• 分配指定一個(gè)無(wú)意義值到輸出波段。

-mo "META-TAG=VALUE"

• 如果可能的話通過設(shè)置一個(gè)原數(shù)據(jù)的鍵來(lái)指定輸出數(shù)據(jù)集原數(shù)據(jù)對(duì)應(yīng)的值。

-co "NAME=VALUE"

• 通過一個(gè)創(chuàng)建參數(shù)來(lái)指定輸出格式特殊創(chuàng)建要求。多個(gè)-co 參數(shù)可以組合起來(lái)使用。創(chuàng)建參數(shù)可以參考個(gè)個(gè)數(shù)據(jù)格式本身說(shuō)明。

-gcp pixel line easting northing

• 添加指定地面控制點(diǎn)到輸出數(shù)據(jù)集。這個(gè)參數(shù)可能被提供多次以提供一系列的GCPs

-quiet

• 忽略處理進(jìn)程監(jiān)視和其他不是錯(cuò)誤的輸出。

-sds

• 拷貝文件中所有子數(shù)據(jù)集到各自的輸出文件中。通常這個(gè)參數(shù)用在HDF或者OGDI這樣有子數(shù)據(jù)集的格式中。

src_dataset

• 源柵格數(shù)據(jù)文件名

dst_dataset

• 目標(biāo)柵格數(shù)據(jù)文件民

例子：

gdal_translate -of GTiff -co "TILED=YES" utm.tif utm_tiled.tif

4.3. gdaladdo

建立或者重建圖像金字塔。

用法：

gdaladdo [-r {nearest,average,average_mp,average_magphase,mode}] filename levels

gdaladdo工具可以用于為大多數(shù)支持的格式建立或者重建金字塔?？梢允褂孟旅鎺追N重采樣算法中的一種來(lái)進(jìn)行縮小重采樣操作。

-r {nearest, average,average_mp, average_magphase,mode}:

• 選擇重采樣算法。

filename

• 要建立金字塔的文件名。

levels:

• 要建立略縮圖的層號(hào)的列表。

Mode不是一個(gè)確切的實(shí)現(xiàn)。而且average_mp不適合使用。Average_magphase在一個(gè)mag/phase空間中平均聯(lián)合數(shù)據(jù)（譯者：這句涉及算法，我不懂這個(gè)是什么平均方法）。最臨近和平均值法可以用于普通的影像數(shù)據(jù)。最臨近法用一個(gè)最臨近的鄰居(簡(jiǎn)單象元)重采樣，而平均值法計(jì)算所有非無(wú)意義值的象元的平均值來(lái)進(jìn)行重采樣。

選擇一個(gè)縮放水平如2表示略縮圖縮放程度是源圖像每個(gè)維上分辨率的1/2。如果文件在所選縮放水平上已經(jīng)存在略縮圖，那么那個(gè)縮放程度將被重新計(jì)算并覆蓋寫入。

一些格式根本不支持金字塔。許多格式在文件以外以擴(kuò)展名.ovr存儲(chǔ)金字塔，TIFF就是如此。GeoTIFF格式直接把金字塔存儲(chǔ)到原有的文件中。

在TIFF中創(chuàng)建金字塔可以通過用COMPRESS_OVERVIEW配置參數(shù)進(jìn)行壓縮。所有GeoTIFF支持的壓縮方法，可以在這里獲得（如： --config COMPRESS_OVERVIEW DEFLATE）。

大多數(shù)驅(qū)動(dòng)也支持一個(gè)備用的略縮圖格式（使用的是Erdas圖像格式）。引發(fā)這個(gè)備用格式使用 USE_RRD=YES 來(lái)設(shè)置參數(shù)。這樣做會(huì)把GDAL程序創(chuàng)建的金字塔放到一個(gè)輔助的.aux文件中使得可以該金字塔可以直接在Erdas中使用或者也可以在ArcGIS中使用。

例子：

在所提供的TIFF文件內(nèi)部創(chuàng)建金字塔

gdaladdo -r average abc.tif 2 4 8 16

從一個(gè)ERDAS.IMG文件中創(chuàng)建一個(gè)外部的壓縮的金字塔文件。

gdaladdo --config COMPRESS_OVERVIEW DEFLATE erdas.img 2 4 8 16

為給定JPEG文件創(chuàng)建一個(gè)Erdas Imagine 格式金字塔

gdaladdo --config USE_RRD YES airphoto.jpg 3 9 27 81

4.4. gdalwarp

簡(jiǎn)單圖像重投影和綁定工具。

用法

gdalwarp [-s_srs srs_def] [-t_srs srs_def] [-order n] [-et err_threshold] [-te xmin ymin xmax ymax] [-tr xres yres] [-ts width height] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16] [-srcnodata "value [value...]"] [-dstnodata "value [value...]"] [-rn] [-rb] [-rc] [-rcs] [-wm memory_in_mb] [-multi] [-q] [-of format] [-co "NAME=VALUE"]* srcfile* dstfile

gdalwarp工具是一個(gè)圖像鑲嵌、重投影、和綁定的工具。程序可以重投影到任何支持的投影，而且如果圖像("raw" with)控制信息也可以把GCPs和圖像存儲(chǔ)在一起。

-s_srs srs_def

• 源空間參考集。所有可以使用OGRSpatialReference.SetFromUserInput()來(lái)調(diào)用的坐標(biāo)系統(tǒng)都可以使用。包括EPSG PCS和GCSes（例如EPSG:4296），PROJ.4描述（如上）?；蛘甙谋镜囊詐rf為擴(kuò)展名的文件。

-t_srs srs_def

• 目標(biāo)坐標(biāo)系統(tǒng)集。（可以參考上面的解釋）

-order n

• 用于綁定的多項(xiàng)式規(guī)則。默認(rèn)的是選擇一個(gè)基于GCPs數(shù)量的多項(xiàng)式。

-tps

• 允許根據(jù)已有的GCPs使用薄板內(nèi)插轉(zhuǎn)換方法?？梢杂眠@個(gè)來(lái)替代-order參數(shù)。

-et err_threshold

• 轉(zhuǎn)換的錯(cuò)誤臨界的近似值。（以象元單位- 默認(rèn)到0.125）

-te xmin ymin xmax ymax

• 設(shè)置被創(chuàng)建的輸出文件的地理邊界范圍

-tr xres yres

• 設(shè)置輸出文件分辨率（單位以目標(biāo)地理參考為準(zhǔn)）

-ts width height

• 設(shè)置輸出文件大?。ㄒ孕辛卸嗌傧笤?jì)量）

-wo "NAME=VALUE"

• 設(shè)置綁定參數(shù)。GDALWarpOptions::papszWarpOptions 不會(huì)顯示所有的參數(shù)，多個(gè)-wo參數(shù)可以并列。

-ot type

• 設(shè)置輸出波段的數(shù)據(jù)類型。

-wt type

• 工作的象元數(shù)據(jù)類型。包括在源圖像和目標(biāo)圖像緩沖中的象元數(shù)據(jù)類型。

-rn

• 用最臨近發(fā)進(jìn)行重采樣（默認(rèn)，最快，但是內(nèi)插質(zhì)量最差）

-rb

• 用雙線性法進(jìn)行重采樣

-rc

• 用立方體發(fā)進(jìn)行重采樣

-rcs

• 用立方曲線發(fā)進(jìn)行重采樣（最慢的方法）

-srcnodata value [value...]

• 設(shè)置輸入波段的無(wú)意義值掩膜（如果設(shè)置多個(gè)值，就可以支持每個(gè)波段設(shè)一個(gè)）。如果有多個(gè)值，就需要把他們用雙引號(hào)括起來(lái)，以保持在命令參數(shù)中作為單一參數(shù)輸入。掩膜值不會(huì)在內(nèi)插中處理。

-dstnodata value [value...]

• 設(shè)置輸出波段的無(wú)意義值掩膜（可以參考上面）。注意，如果可能的化，新的文件將初始化這個(gè)值，并記錄在輸出文件中。

-wm memory_in_mb

• 設(shè)置內(nèi)存總量（以兆字節(jié)為單位），如果綁定API允許使用高速緩存。

-multi

• 使用多線程綁定實(shí)現(xiàn)。多線程將被用來(lái)處理打快的圖像并同時(shí)處理輸入輸出操作。

-q

• 不打印輸出。

-of format

• 選擇輸出格式，默認(rèn)是GeoTiff

-co "NAME=VALUE"

• 通過一個(gè)創(chuàng)建參數(shù)來(lái)控制創(chuàng)建新的文件。多個(gè)'-co'可以排列起來(lái)使用。參數(shù)可用的值可以察看每個(gè)驅(qū)動(dòng)。

srcfile

• 源文件名

dstfile

• 目標(biāo)文件名

如果輸出文件已經(jīng)存在，則在已存在的文件中創(chuàng)建地圖集。

例子：

一個(gè)帶有用經(jīng)緯度標(biāo)記的邊界控制點(diǎn)的8位GeoTiff格式的Spot影像可以通過下面的命令綁定到一個(gè)UTM投影上。

gdalwarp -t_srs '+proj=utm +zone=11 +datum=WGS84' raw_spot.tif utm11.tif

一個(gè)帶有用經(jīng)緯度標(biāo)記的邊界控制點(diǎn)的兩個(gè)波段的HDF格式的ASTER影像可以通過下面的命令綁定到UTM投影上。

gdalwarp HDF4_SDS:ASTER_L1B:"pg-PR1B0000-2002031402_100_001":2 pg-PR1B0000-2002031402_100_001_2.tif