GDAL的初步介绍和使用

已于 2024-07-04 09:38:04 修改
阅读量696 收藏 17
点赞数 22
分类专栏： GDAL 文章标签： 开源软件
于 2024-07-04 09:26:05 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/kappan/article/details/140170201
版权
GDAL基本介绍
  1. GDAL全称是Geospatial Data Abstraction Library（地理空间数据抽象库），是使用C/C++语言编写的用于读写空间数据的一套跨平台开源库，它利用抽象数据模型来表达所支持的各种文件格式，还使用一系列命令行来进行数据转换和处理。现有的大部分GIS或者遥感平台，不论是商业软件ArcGIS，ENVI还是开源软件GRASS，QGIS，都使用了GDAL作为底层构建库。
  2. GDAL最初是由Frank Warmerdam于1998年开始开发的，在GDAL1.3.2版本之后，正式有开源空间信息基金会(Open Source Geospatial Foundation，OSGeo)下的GDAL/OGR项目委员会对其进行维护。
  3. GDAL库由OGR和GDAL项目合并而来，OGR主要用于空间要素矢量矢量数据的解析，GDAL主要用于空间栅格数据的读写。此外，空间参考及其投影转换使用开源库 PROJ.4进行。
  4. 目前，GDAL主要提供了三大类数据的支持：栅格数据，矢量数据以及空间网络数据（Geographic Network Model）。栅格数据格式包括Arc/Info ASCII Grid（ASC），GeoTiff（Tiff），Erdas Imagine Images（Img），ASCII DEM（DEM）等。矢量数据格式主要包括ESRI Shapefiles、S-57、SDTS、PostGIS、Oracle Spatial、Mapinfo mid/mif和Mapinfo TAB等。
  5. GDAL提供了C/C++接口，并且通过SWIG提供了Python，Java，C#等的调用接口。当我们在Java中调用GDAL的API函数时，其实底层执行的是C/C++编译的二进制文件。
GDAL 前置GIS基础
1. 栅格数据

图片格式的 GIS 数据我们统称为栅格数据，在 GIS 领域里，一切的图片，图像，影片，影像格式的数据都称作是栅格数据（图像数据）。

2. 矢量数据

矢量数据是指一些人工手绘的、测绘而得来的一些图形数据或者是由栅格数据中得出。这种数据往往只有一些简单的轮廓构成，比如一个行政区可能只用轮廓边界表示，一条河流可能只用一条线段表示。这就是我们的矢量数据（图形数据）。

GIS 开发过程中，矢量数据的主要格式有 shape 格式和 geojson 格式，wkt格式。一种是文件形式，一种是网络传输中的json形式，wkt是一种精简化数据格式，只记录地理数据，不记录属性数据；

3. 坐标系和投影

  1. 在地理学中，坐标系分为地理坐标系和投影坐标系。

  2. 地理坐标系又称为经纬度坐标系。地理坐标系的核心就是用经纬度来表示位置。不管是在球面上还是把球面投影到平面上，如果你想用地理坐标系来表示一幅地图的任意位置，那这个位置的表示方式一定是经纬度。投影坐标系是指基于地理坐标系进行的一些投影变换，投影坐标系的单位通常是米（m）或者千米，即距离长度为单位。每个地理坐标系都可以通过数据计算和变换产生出很多的投影坐标系。

  3. 在国内，常用的地理坐标系有两个，分别是 wgs84 坐标系和 CGCS2000 坐标系。常用的投影坐标系有很多，最主要的是 web 墨卡托投影坐标系，这个坐标系是经过 wgs84 坐标系投影而来的。

  4. 在实际开发过程中坐标系经常用EPSG 代码来表示，EPSG 是欧洲一个石油组织的简称，这个组织提出了一种规范，用代号的形式来表示坐标系。 wgs84 坐标系我们会表示成 EPSG:4326。CGCS2000 坐标系我们会表示成 EPSG:4490。web 墨卡托坐标系我们会表示成 EPSG:900913 或者 EPSG:3857 等等。所以当在开发的过程中遇到了这些代号，需要做一个转换时，清楚其要表示的坐标系即可。

GDAL体系结构
GDAL体系
  1. GDAL使用抽象数据模型(abstract data model)来解析它所支持的数据格式，抽象数据模型包括数据集(dataset)，坐标系统，仿射地理坐标转换(Affine Geo Transform)， 大地控制点(GCPs)， 元数据(Metadata)，栅格波段(Raster Band)，颜色表(Color Table)，子数据集域(Subdatasets Domain)，图像结构域(Image_Structure Domain)，XML域(XML:Domains)。
  2. GDALMajorObject类：带有元数据的对象。
  3. GDALDdataset类：通常是从一个栅格文件中提取的相关联的栅格波段集合和这些波段的元数据;GDALDdataset也负责所有栅格波段的地理坐标转换(georeferencing transform)和坐标系定义。
  4. GDALDriver类：GDAL会为每一个所支持的文件格式创建一个该类的实体，来管理该文件格式。
  5. GDALDriverManager类：用来管理GDALDriver类。
OGR体系
  1. Geometry类：Geometry (包括OGRGeometry等类)封装了OpenGIS的矢量数据模型，并提供了一些几何操作，WKB(Well Knows Binary)和WKT(Well Known Text)格式之间的相互转换，以及空间参考系统(投影)。
  2. Spatial Reference类：OGRSpatialReference封装了投影和基准面的定义。
  3. Feature类：OGRFeature封装了一个完整feature的定义，一个完整的feature包括一个geometry和geometry的一系列属性。
  4. Feature Definition类：OGRFeatureDefn里面封装了feature的属性，类型、名称及其默认的空间参考系统等。一个OGRFeatureDefn对象通常与一个层(layer)对应。
  5. Layer类：OGRLayer是一个抽象基类，表示数据源类OGRDataSource里面的一层要素(feature)。
  6. Data Source类：OGRDataSource是一个抽象基类，表示含有OGRLayer对象的一个文件或一个数据库。
  7. Drivers类：OGRSFDriver对应于每一个所支持的矢量文件格式。类OGRSFDriver由类 OGRSFDriverRegistrar来注册和管理。
GDAL的快速上手
本地开发环境的配置

​ 参考gdal的配置文档即可

demo1:对tiff的读取
  1. 使用GDAL第一步，注册gdal驱动
import org.gdal.gdal.gdal;
gdal.AllRegister();
  1. 读取数据源信息，以只读的方式打开文件
//载入栅格，读取相关信息
Dataset rasterDataset = gdal.Open("D:\\gdal\\class.tif", gdalconstConstants.GA_ReadOnly);

此时，Tiff文件中的数据已经读到了rasterDataset中，后续的数据操作将围绕rasterDataset进行。
在这里插入图片描述

  1. 最后需要释放资源
 rasterDataset.delete();
  1. 完整代码
    public static void main(String[] args) {
        String inRaster = "D:\\gdal\\class.tif";
        String outShp = "D:\\gdal\\class.shp";
        //注册驱动，否则下边执行报错
        gdal.AllRegister();
        //载入栅格，读取相关信息
        Dataset rasterDataset = gdal.Open(inRaster, gdalconstConstants.GA_ReadOnly);
        //栅格的波段信息
        Band band = rasterDataset.GetRasterBand(1);
        //设置坐标系
        SpatialReference prj = new SpatialReference();
         prj.ImportFromEPSG(4326);
        //释放资源
        rasterDataset.delete();
    }
demo2:对shp数据的读取

  1. ogr全局注册和全局配置

    // 指定文件的名字和路径
String strVectorFile = "F:\\vector_data\\道路中心线.shp";
// 注册所有的驱动
ogr.RegisterAll();
// 为了支持中文路径，请添加下面这句代码
gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");
// 为了使属性表字段支持中文，请添加下面这句
gdal.SetConfigOption("SHAPE_ENCODING", "CP936");

  2. 获取shp驱动，打开shp文件

    // 读取数据，这里以ESRI的shp文件为例
String strDriverName = "ESRI Shapefile";
// 创建一个文件，根据strDriverName扩展名自动判断驱动类型
 
org.gdal.ogr.Driver oDriver = ogr.GetDriverByName(strDriverName);
 
if (oDriver == null) {
	System.out.println(strDriverName + " 驱动不可用！\n");
	return;
}
DataSource dataSource = oDriver.Open(strVectorFile);

  3. 获取图层、空间参考、图层范围

    //Layer layer = dataSource.GetLayer("test");
Layer layer = dataSource.GetLayer(0);
		
for(int i = 0;i<dataSource.GetLayerCount();i++) {
	Layer layerIdx = dataSource.GetLayer(i);
	System.out.println("图层名称：<==>" + layerIdx.GetName());
}
		
String layerName = layer.GetName();
System.out.println("图层名称：" + layerName);
SpatialReference spatialReference = layer.GetSpatialRef();
//System.out.println(spatialReference);
System.out.println("空间参考坐标系：" + spatialReference.GetAttrValue("AUTHORITY", 0)
				+ spatialReference.GetAttrValue("AUTHORITY", 1));
 
double[] layerExtent = layer.GetExtent();
 
System.out.println("图层范围：minx:" + layerExtent[0] + ",maxx:" + layerExtent[1] + ",miny:" + layerExtent[2] + ",maxy:" + layerExtent[3]);

  4. 获取图层属性信息

    FeatureDefn featureDefn = layer.GetLayerDefn();
 
int fieldCount = featureDefn.GetFieldCount();
 
Map<String,String> fieldMap = new HashMap<String,String>();
for (int i = 0; i < fieldCount; i++) {
	FieldDefn fieldDefn = featureDefn.GetFieldDefn(i);
	// 得到属性字段类型
	int fieldType = fieldDefn.GetFieldType();
	String fieldTypeName = fieldDefn.GetFieldTypeName(fieldType);
	// 得到属性字段名称
	String fieldName = fieldDefn.GetName();
	fieldMap.put(fieldTypeName, fieldName);
}
		
System.out.println(fieldMap);

  5. 读取空间空间信息及属性列表

    System.out.println(layer.GetFeature(1).GetGeometryRef().ExportToJson());
System.out.println(layer.GetFeature(2).GetGeometryRef().ExportToJson());
System.out.println(layer.GetFeature(3).GetGeometryRef().ExportToJson());
		
for (int i = 0; i < featureCount; i++) {
	Feature feature = layer.GetFeature(i);
	Object[] arr = fieldMap.values().toArray();
	for (int k = 0; k < arr.length; k++) {
		String fvalue = feature.GetFieldAsString(arr[k].toString());
		System.out.println(" 属性名称:" + arr[k].toString() + ",属性值:" + fvalue);
	}
}

    结果如下

    img

  6. 关闭资源

    dataSource.delete();
demo3 栅格矢量化

说明：GDAL中内嵌了矢量化的算法，直接调用即可

    public static void main(String[] args) {
        String inRaster = "D:\\gdal\\class.tif";
        String outShp = "D:\\gdal\\class.shp";
        //注册驱动，否则下边执行报错
        gdal.AllRegister();
        //载入栅格，读取相关信息
        Dataset rasterDataset = gdal.Open(inRaster, gdalconstConstants.GA_ReadOnly);
        //栅格转矢量需要的波段信息
        Band band = rasterDataset.GetRasterBand(1);
        //设置坐标系
        SpatialReference prj = new SpatialReference();
        if (rasterDataset.GetProjectionRef().isEmpty()) {
            prj.ImportFromEPSG(4326);
        } else {
            //栅格数据的坐标系作为矢量化后的坐标系
            prj.ImportFromWkt(rasterDataset.GetProjectionRef());
        }
        //创建输出矢量
        org.gdal.gdal.Driver driver = gdal.GetDriverByName("ESRI Shapefile");
        Dataset shpDataset = driver.Create(outShp, rasterDataset.getRasterXSize(), rasterDataset.getRasterYSize(), 1, gdalconst.GDT_Byte);
        Layer layer = shpDataset.CreateLayer("layer", prj);
        FieldDefn field = new FieldDefn("value", ogr.OFTInteger);//创建一个字段用来存储栅格的像素值
        layer.CreateField(field, 1);
        //矢量化
        gdal.Polygonize(band, null, layer, 0);
        log.info("矢量化过程完成");
        layer.SyncToDisk();
        //释放资源
        shpDataset.delete();
        rasterDataset.delete();
    }
demo4 叠加之相交分析

说明：GDAL具有很多空间分析功能，感兴趣可以调用研究

public class VectorIntersection {
    public static void main(String[] args) {
        // 初始化GDAL库
        ogr.RegisterAll();

        // 加载输入的矢量数据
        DataSource src1 = ogr.Open("D:\\gdal\\geo_xzqh_cun_anhui\\GEO_XZQH_CUN.shp", 0);
        DataSource src2 = ogr.Open("D:\\gdal\\HS_H09_20240406_0900_anh_new.shp", 0);

        if (src1 == null || src2 == null) {
            System.out.println("无法打开矢量数据源");
            return;
        }

        // 获取操作的图层
        Layer layer1 = src1.GetLayerByIndex(0);
        Layer layer2 = src2.GetLayerByIndex(0);

        // 设置相交后结果的空间参考
        SpatialReference spatialReference = new SpatialReference();
        spatialReference.SetWellKnownGeogCS("WGS84");

        // 创建输出矢量数据源
        Driver driver = ogr.GetDriverByName("ESRI Shapefile");
        DataSource dst = driver.CreateDataSource("D:\\gdal\\outlay.shp", null);

        // 执行相交操作
        layer1.Intersection(layer2, dst.GetLayerByIndex(0));

        // 关闭数据源
        dst.SyncToDisk();
        dst.delete();
        src1.delete();
        src2.delete();

        System.out.println("相交操作完成，结果保存在 'outlay.shp'");
。。。后续待更新优化
开啊开阿凯
关注
专栏目录
新手使用GDAL详细教程
幸福的小骁仔的博客
10-02 4万+
文章是po主参考了很好的GDAL初学者资料，VS2010中项目配置引入GDAL，经过自己的实际操作检验归纳总结的，仅属于半原创。         首先介绍一下实验的环境，win7下32位系统，预先配置了VS2013。         步骤一：下载GDAL源码。 下载地址：http://trac.osgeo.org/gdal/wiki/DownloadSource，进入下面这个页面：
Python与GDAL初步（一）：栅格读写
sievert的博客
08-12 430
一、读 from osgeo import gdal dataset = gdal.Open(dir) #dir--文件全路径 将影像打开为osgeo.gdal.Dataset，可以通过一些方法或函数查看其属性 dataset.RasterCount #获取波段数目 dataset.RasterXSize #获取影像宽度 dataset.RasterYSize #获取影像高度 dataset.GetProjection() #获取投影信息 dataset.GetGeoTransform() #获取仿射变
GDAL库简介及函数说明
学习历程记录和分享
03-01 4890
GDAL（Geospatial Data Abstraction Library）是一个开源的地理空间数据处理库，提供了一系列用于读取、写入和处理各种地理空间数据格式的功能。它是一个跨平台的库，可以在多种操作系统上使用，包括 Windows、Linux 和 macOS。
【小沐学GIS】GDAL库安装和使用（C++、Python）
爱看书的小沐
05-25 6580
GDAL 是用于栅格和矢量地理空间数据格式的转换器库，由开源地理空间基金会在 MIT 风格的开源许可下发布。作为一个库，它向调用应用程序提供单个栅格抽象数据模型和单个矢量抽象数据模型，用于所有支持的格式。它还带有各种有用的命令行实用程序，用于数据转换和处理。╮(￣▽￣)╭如果您感觉方法或代码不咋地//(ㄒoㄒ)//，就在评论处留言，作者继续改进；o_O???如果您需要相关功能的代码定制化开发，可以留言私信作者；(✿◡‿◡)！！
QT/C++中的GDAL多线程应用（读取）：发生的问题以及解决方案
最新发布
NCHUkiss的博客
08-22 1849
详细记录了在使用GDAL库对TIF文件进行多线程切割和创建瓦片金字塔时遇到的问题及解决方案。从最初的常规多线程实现开始，经历了使用互斥锁、为每个线程创建独立数据集等多个优化阶段，最终发现通过设置GDAL库的环境变量来启用多线程读取功能才是最有效的解决方案。深入分析了每个阶段的问题和优化效果。
GDAL读写矢量数据-Python
热门推荐
忧郁王子的专栏
08-30 1万+
1、GDAL（Geospatial Data Abstraction Library），源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理 OGR提供对矢量数据格式的读写支持，它所支持的文件格式包括：ESRI Shapefiles， S-57， SDTS， PostGIS，Oracle Spatial， Mapinfo mid/mif
GDAL库使用方法.docx
11-01
GDAL库的使用文档，word版，包含了gdal库的简单使用（python入门学习）
GDAL的具体使用说明情况
05-14
GDAL使用情况，来自GDAL官方网站，详细流程见代码
python GDAL库的简单使用
HengXiaa
03-19 1362
python中GDAL库的3处基础用法
2.GDAL简介
XiaoCai
06-15 4877
1. 什么使GDAL GDAL全称是Geospatial Data Abstraction Library(地理空间数据抽象库),在一个X/MIT许可协议下读写空间数据(包括栅格数据和矢量数据)的开源库，它利用抽象数据模型来表达所支持的各种文件格式，还使用一系列命令行来进行数据转换和处理。 OGR(OGR Simple Features Library)是GDAL项目的一个分支，功能与GDAL类似，只是它提供对矢量数据的读写支持。同时它实现了一个对空间参考信息进行处理的类，用来对空间数据的空间信息进行处理。
gdal用法总结
蓝蓝的天
10-21 4210
USAGE OF GDAL RASTER API Import gdal from osgeo import gdal Open the file Dataset=gdal.Open(filename) Getting dataset information Xsize=dataset.RasterXsize Ysize=dataset.RasterYSize Count=dataset.Ra...
gdal库使用范例
12-25
本程序组要包含了GDAL库文件的使用，适合初学或者偶一定基础的遥感图像开发人员使用
gdal开发文档,gdal库的使用,C,C++
09-10
gdal中文版开发文档，帮助学习gdal开发库
GDAL入门教程（偏向使用文档）
03-18
入门教程，偏向使用文档。建议阅读者边看边写。找到自己所需要的那部分即可。
调用GDAL库实现图像的简单处理
01-03
调用GDAL库实现图像的读取、显示、均衡化、线性拉伸等功能
gdal工具exe
08-28
gdal exe工具包，可实现常规的影像处理操作，方便快捷，高效！
python gdal + skimage实现基于遥感影像的传统图像分割及合并外加矢量化
12-21
总的来说，这篇博客介绍了如何结合GDAL和skimage进行遥感影像的处理，包括读取、分割、合并以及初步的矢量化步骤。然而，实际的矢量化过程可能更为复杂，需要处理几何转换、拓扑关系以及数据优化等问题。对于初学者...
gdal 读取点要素_node-gdal实践-读写shp数据
weixin_35775778的博客
02-23 497
gdal用于读写矢量和栅格数据集，gdal常常使用c，c++做开发，也有部分其他语言扩展如py，node。但由于其他语言的官方例子不多，且存在很多问题，所以入门非常难。本文主要面向使用node语言做工具且对GIS软件有二次开发基础的同学，介绍使用gdal读写经典shp的入门教程。一 node-gdal安装Github地址：https://github.com/naturalatlas/node-g...
gdal学习程序
11-06
gdal初步学习程序，包含基本常用的代码，可供参考
GDAL maven依赖和使用示例
07-27
要在Maven项目中使用GDAL，你需要添加以下依赖项到你的`pom.xml`文件中： ```xml <dependency> <groupId>org.gdal</groupId> <artifactId>gdal</artifactId> <version>3.0.1</version> </dependency> ``` 这将添加GDAL的Java绑定到你的项目中。请确保将`version`属性设置为你想要使用的GDAL版本。 下面是一个简单的示例代码，演示如何使用GDAL进行坐标转换： ```java import org.gdal.gdal.Dataset; import org.gdal.gdal.gdal; import org.gdal.osr.CoordinateTransformation; import org.gdal.osr.SpatialReference; public class GDALExample { public static void main(String[] args) { // 初始化GDAL gdal.AllRegister(); // 打开输入数据集 Dataset inputDataset = gdal.Open("input.tif"); // 获取输入数据集的空间参考 SpatialReference inputSRS = new SpatialReference(inputDataset.GetProjection()); // 定义输出的EPSG:4479空间参考 SpatialReference outputSRS = new SpatialReference(); outputSRS.ImportFromEPSG(4479); // 创建坐标转换对象 CoordinateTransformation transform = new CoordinateTransformation(inputSRS, outputSRS); // 定义待转换的坐标 double[] srcCoord = new double[]{1234567.0, 9876543.0}; double[] targetCoord = new double[2]; // 进行坐标转换 transform.TransformPoint(targetCoord, srcCoord); // 输出转换后的坐标 System.out.println("Transformed coordinates: " + targetCoord[0] + ", " + targetCoord[1]); // 关闭数据集 inputDataset.delete(); } } ``` 在上述示例代码中，我们首先调用`gdal.AllRegister()`来初始化GDAL。然后，我们打开输入数据集（例如一个GeoTIFF文件），并获取其空间参考信息。接下来，我们创建一个输出的EPSG:4479空间参考，并使用`SpatialReference.ImportFromEPSG()`方法导入EPSG代码。然后，我们创建一个`CoordinateTransformation`对象，将输入和输出的空间参考传递给它。最后，我们定义待转换的坐标，并使用`CoordinateTransformation.TransformPoint()`方法进行坐标转换。 请确保将示例中的`input.tif`替换为你自己的输入文件路径，并根据需要调整EPSG代码和坐标转换的参数。 希望这可以帮助到你。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值