原创 Python | 圖像座標、地理座標和投影座標的相互轉化

from osgeo import gdal from osgeo import osr import numpy as np def getSRSPair(dataset): ''' 獲得給定數據的投影參考系和地理

原创 Python 30段代碼分享

學 Python 怎樣才最快,當然是實戰各種小項目,只有自己去想與寫,才記得住規則。本文是 30 個極簡任務,初學者可以嘗試着自己實現;本文同樣也是 30 段代碼,Python 開發者也可以看看是不是有沒想到的用法。   Python 是

原创 Python 庫(詳細整理)簡要說明用途

庫名稱簡介 Chardet 字符編碼探測器,可以自動檢測文本、網頁、xml的編碼。colorama 主要用來給文本添加各種顏色,並且非常簡單易用。Prettytable 主要用於在終端或瀏覽器端構建格式化的輸出。difflib,[Pyth

原创 GDAL+Python | 實現柵格影像處理之柵格矢量化及矢量柵格化

柵格矢量化:將柵格數據生成矢量數據。 矢量柵格化:將矢量數據生成柵格數據。 使用方法 柵格矢量化 使用gdal.Polygonize()方法 參數說明(未列完) 參數    說明 srcBand    輸入數據波段 maskBand 

原创 Python學習知識清單(基礎+進階)

Python學習知識清單 基礎語法。 Pyhton與C語言有一定的聯繫,但要注意一下兩種思維方式的不同。 列表的使用。 定義空列表。 List = [ ]   1 列表的長度 len(list) 1 列表的組合 list1 += list

原创 GDAL & Python | Chapter 3. Reading and writing vector data

Chapter 3. Reading and writing vector data This chapter covers Understanding vector data Introducing OGR Reading vector

原创 4. GDAL python教程(3)——過濾器,簡單的空間分析,函數和模塊

  4.1. 屬性過濾器Attribute filters Layer對象有一個方法叫SetAttributeFilter(<where_clause>)可以將Layer中符合某一條件的Feature過濾出來。設定了Filter之後就可以

原创 3. GDAL python教程(2)——幾何形狀geometry與投影projection

  3.1. 建立新的幾何形狀 建立空的geometry對象:ogr.Geometry 定義各種不同的geometry使用的方法是不一樣的(point, line, polygon, etc) 新建點point,使用方法AddPoint(

原创 5. GDAL python教程(4)——用GDAL讀取柵格數據

  GDAL原生支持超過100種柵格數據類型,涵蓋所有主流GIS與RS數據格式,包括 ArcInfo grids, ArcSDE raster, Imagine, Idrisi, ENVI, GRASS, GeoTIFF HDF4, HD

原创 1. 地理數據處理軟件包GDAL簡介

  GDAL(Geospatial Data Abstraction Library)是一個在X/MIT許可協議下的開源柵格空間數據轉換庫。該項目由Frank Warmerdam教授於1998年發起。 它利用抽象數據模型來表達所支持的各種

原创 6. GDAL python教程(5)——地圖代數與柵格數據的寫入

  6.1. 以計算NDVI爲例: NDVI=(NIR-RED)/(NIR+RED) 其中NIR爲波段3,RED爲波段2 編程要點如下: 將波段3讀入數組data3,將波段2讀入數組data2 計算公式爲: ndvi = (data3

原创 GDAL & Python Chapter 2. Python basics

  This chapter covers Using the Python interpreter vs. writing scripts Using the core Python data types Controlling the

原创 2. GDAL python教程(1)——用OGR讀寫矢量數據

  本教程的講義和源碼都是取自Utah State University的openGIS課程 相關資料,包括講義、源碼、數據樣例,請從此處下載 本人只是做點翻譯,寫寫學習體會而已,版權屬於原作者。 歡迎轉載,不過別忘了上面這段話。 2.1

原创 7. GDAL python教程(6)——更多柵格數據處理函數

  柵格數據的投影 首先要知道輸入投影和輸出投影的WKT(Well Known Text),可以通過GetProjection()讀到,可以用SpatialReference對象創建。 用下面的語句新建柵格數據集並重新投影,投影結果輸出到

原创 8. GDAL python教程(7)——雜七雜八的話

  介紹工具FWTools,命令行操作等等 FWTools包括如下命令: ogrinfo:打印矢量圖層的信息 ogr2ogr:矢量數據格式轉換 gdalinfo:打印柵格圖層的信息 gdal_translate:柵格數據格式轉換 gdal