其他

在部落格本文開始之前，先祝大家在新的一年裡工作順利、萬事如意。新年新氣象，希望在新的一年裡能夠和大家多多交流，共同學習進步。GDAL虛擬檔案大概從1.6.0版本開始，最開始主要是用來訪問zip或gz等壓縮包歸檔檔案中的資料。距離詳細可以參考網址：http://trac.osgeo.org/gdal/wiki/UserDocs/ReadInZip。最開始只支援兩種虛擬檔案系統，/vsigzip和/viszip。對於.zip檔案只支援讀取，而/visgzip支援讀取和順序寫入。但是寫入的速度有些慢。

原文地址：http://blog.thematicmapping.org/2012/07/terrain-mapping-with-mapnik.html作者：作者：Bj?rnSandvik（Google

在使用GDAL編寫演算法的時候，經常會將計算的中間結果存在一個臨時的影像檔中，然後使用完再將其刪除，如果臨時檔案就一個的話，建立一個也無所謂，但是當一個複雜的演算法中可能會出現很多個臨時檔案的時候（我在編寫Hariss角點自動匹配演算法的時候有4個臨時檔案），這種情況下總覺得臨時檔案很不爽，此外第一個不爽的地方；第二個映像太大的時候，臨時檔案也會佔用很大的空間，假如空間不足或者給定的臨時檔案路徑不可寫等問題會讓人頭疼；第三，在建立臨時檔案的讀寫上會耗用比較多的時間，尤其在磁碟的IO時，耗時比較多

在上一篇部落格中，大概說了下怎麼使用gdal提供的gdalwarp工具來進行校正處理。其實質與envi的glt校正應該是一樣的。我把gdalwarp的代碼封裝了一下，寫了一個類來進行geoloc處理。希望對大家有用。先是標頭檔，函數的注釋很詳細，就不多說了。後面的源檔案就是從gdalwarp.cpp中摘錄出來的，有興趣的可以看gdalwarp.cpp ，下面的代碼只是把這個檔案中沒有用到的代碼刪除了，同時有些參數直接寫死了。不太清楚的可以直接看原來的代碼。/*******************

在去年的時候，在GDAL的郵件清單中發現有人提出要在GDAL中添加自動提取匹配點（GCP）的演算法，在經過了將近半年的時間，GDAL中的匹配演算法終於有了實質性的進展。在目前GDAL的開發版本（GDAL1.10beta）中已經提交，該匹配演算法是基於SURF演算法來實現的。實現該演算法的人叫Andrew

    在工作中，會經常使用一個行政區的向量圖去裁剪一個遙感影像圖，得到該行政區的影像圖，一般的商業遙感軟體都具有這個功能。今天就是用GDAL來實現這一個很實用的功能。首先用到的是GDAL中的gdalwarp，又是warp，呵呵，上一篇就是使用warp進行重採樣的。   首先需要用到gdal源碼目錄裡面的app檔案夾下的gdalwarp.cpp檔案中的幾個函數，大概行數是1651行，直到檔案結尾，代碼如下：/********************************************

前言：一直在使用和研究GDAL的相關東西，發現網上對GDAL的內容倒是不少，但是很少有系統的介紹說明，以及內部的一些結構說明，基於這些原因，將本人的一些粗淺的理解放在此處，形成一個系列，暫時名為《GDAL源碼剖析》（名稱有點大言不慚，歡迎大家口水吐之，板磚拍之），供大家交流參考，有什麼錯誤之處，望大家不吝指正，本系列對於GDAL的使用均是在Windows平台下，對於Linux平台下的不在此系列討論範圍之內。此外，轉載本部落格內容，請註明出處，強烈鄙視轉載後不註明出處的“類剽竊”行為。一、GDAL

    GDAL1.8.0發布很久了，一直沒有將自己的工程中的版本更新到1.80。今天將其更新到1.80發現含有中文路徑的檔案都不能開啟了，影像和向量檔案都是。仔細對比了GDAL1.72和GDAL1.80的代碼，終於發現了問題的所在之處，詳細代碼在GDAL_HOME\port\cpl_vsil_win32.cpp檔案中的類VSIWin32FilesystemHandler中，以Stat()函數為例(435行)，其他函數類似。代碼如下：       GDAL1.8.0代碼（部分）：/*******

使用GDAL開啟裸資料首先需要指定裸資料的大小，即波段數，長度，寬度，以及裸資料之間的相片順序等。我們可以使用GDAL提供的虛擬檔案格式來進行對裸資料的讀取。首先，需要建立一個叫vrt的xml格式的檔案，內容如下：<VRTDataset rasterXSize="12000" rasterYSize="12000"> <VRTRasterBand dataType="Byte" band="1" subClass="VRTRawRasterBand">

GDAL中提供了一個訪問網路資料的介面，具體的協助文檔可以訪問這裡（http://gdal.org/cpl__http_8h.html）。不過前提是GDAL需要依賴一個第三方的開源庫libcurl，libcurl主要功能就是用不同的協議串連和溝通不同的伺服器~也就是相當封裝了的sockPHP 支援libcurl(允許你用不同的協議串連和溝通不同的伺服器)。, libcurl當前支援http, https, ftp, gopher, telnet, dict, file, 和ldap 協議。lib

       在處理遙感映像中，發現往往比較耗時的是在資料的IO中，尤其是在O（寫入）的時候更加耗時。GDAL可以支援映像的多線程寫入，下面結合執行個體進行簡單的測試，看看實際效果會不會提高。       在這裡我使用的是boost庫的thread庫來進行多線程建立。下面先使用計算PI來對boost的thread有個簡單的說明。這裡的計時使用的是boost的progress_timer。下面是使用多線程計算PI的一段小代碼，對於多線程計算的10部分沒有加起來。#include

         這篇算是RasterIO的一個示範程式吧，原始碼可以去資源中心去下載。在之前有一篇部落格中（地址為：http://blog.csdn.net/liminlu0314/article/details/7072224），GDALRasterIO函數的使用中，有很多人留言，之後又在文章後面更新了一些例子代碼。但是還是有很多人不是很清楚，就寫個顯示的小程式進行說明下吧。        首先簡單的介紹下，我的介面用的QT，VS2008

之前寫了一篇使用GDAL下載SRTM的資料，只是大概說明了下怎麼下載，文章後面說要寫代碼實現一下，最近太慢，一直拖著，這事肯定不會忘記的，呵呵。今天就把上次剩下的尾巴處理一下。之前的部落格地址：http://blog.csdn.net/liminlu0314/article/details/8068715。在上篇部落格中大概分析了一下，首先要通過代碼構造一個VRT檔案。函數代碼如下：int CreateSRTMVRTFile(const char* pszSrcFile, const char*

一、簡介英文網址：http://www.gdal.org/grid_tutorial.html。網格插值的意思就是從離散的資料點建立一個柵格映像的過程。通常情況下，你有一系列研究地區的離散點，如果你想將這些點轉換為規則的網格資料來進行進一步的處理，或者和其他網格資料進行合并等處理。是網格插值的一個：網格插值使用資料插值和逼近演算法可以用來解決這個問題，但是插值的作用並不僅僅用來處理這個問題。有時候，你並不需要對資料進行插值處理，你需要做的就是計算一下資料覆蓋地區的統計值和其他指標。統計值本身是十

寫在開頭：最近工作太忙，沒時間更新，今天抽空更新一篇，上篇部落格是四，這篇是六，中間還有一篇第五篇，主要將GDAL中的Python工具集，儘快給大家更新。一、簡單的調用關於GDAL的使用，網上的資料都很多，主要還是要熟悉GDAL的組織圖，類以及類的函數等，熟悉了，使用GDAL就不在話下了。最常用的就是動態庫的GDAL，當然你也可以使用靜態庫，這裡只是簡單的介紹使用動態GDAL庫來做開發。首先開啟VS，建立一個工程，控制台的就成。然後在工程的屬性對話方塊中，找到【配置屬性】-【C/C++】-【常規

        使用OGR庫建立DXF格式的資料和建立Shp格式的資料基本一樣，代碼除了註冊驅動不一樣之外，其他的完全一樣。        需要注意的是，OGR建立DXF格式需要data檔案夾下的head.dxf，所以首先要設定GDAL_DATA的目錄。使用下面代碼進行設定，也可以直接設定環境變數，但是比較麻煩，還是用代碼設定比較方便：CPLSetConfigOption("GDAL_DATA","E:\\gdal\\data");       

一、簡介本文原文地址：http://www.gdal.org/warptut.htmlGDAL Warp

    很多遙感衛星資料使用的量化層級都要比8bit高，比如常用的WorldView用的是12bit的量化，對於一些影像處理軟體，是不能直接處理12bit量化的映像，所以需要將12bit的資料轉為16bit資料或者8bit資料來進行處理。   

        GDAL/OGR

本文主要講怎麼使用qmake工具產生的VS的工程檔案，對於QT的庫不做說明。qt的安裝很簡單，一路下一步就OK，安裝完看看qt的bin目錄是否添加中環境變數path中，如果沒有，請手動添加。首先，使用的代碼例子是QT安裝目錄下的qmake的樣本檔案，路徑為：C:\QT\4.8.1\examples\qmake\tutorial。將裡面的代碼拷貝到一個新的檔案夾，以免破壞原來的代碼。我的目錄是F:\Work\Test\QtDemo。開啟cmd，輸入下面的命令“qmake