C# GDAL 学习一

最近一直琢磨如何用C#+GDAL读取栅格数据（.tif或.img），运气不错的在GDAL 的官网上找到一部分源码。经过本人测试，效果还不错。学习还将继续深入下去。

参考网址：http://trac.osgeo.org/gdal/browser/trunk/gdal/swig/csharp/apps/GDALRead.cs

开发环境：VS2008+GDAL1.5

所需dll ：gdal15.dll、gdal_csharp.dll、gdal_wrap.dll、gdalconst_csharp.dll、gdalconst_warp.dll

一、将以上.dll添加到工程bin\debug目录下。

二、建立控制台程序，添加gdal_csharp引用。如图

三、工程GDALRead处右键打开属性对话框，调试一栏添加命令行参数，如图：

这里的命令行参数在程序中直接被数组args[]调用。

四、任务与目标

（1）、读取栅格数据的一般参数，如坐标投影(Projection)、波段数(Rsatercount)、数据驱动、栅格大小(RasterSize)

(2)、每个波段的数据类型(DataType)、大小(Size)、PaletteInterp

五、完整代码

View Code

1 using System;
2 using System.Collections.Generic;
3 using System.Linq;
4 using System.Text;
5 using OSGeo.GDAL;
6 namespace GDALRead
7 {
8     class Program
9     {
10         public static void usage()
11         {
12             Console.WriteLine("usage");
13             System.Environment.Exit(-1);
14         }
15      public  static void Main(string[] args)
16         {
17             //Console.WriteLine(args[0]);
18             //Console.ReadLine();
19          try
20          {
21              Gdal.AllRegister();
22              Dataset ds = Gdal.Open(args[0],Access.GA_ReadOnly);
23              if (ds==null)
24              {
25                 Console.WriteLine("Can't open " + args[0]);
26                  System.Environment.Exit(-1);
27              }
28              Console.WriteLine("raster dataset parameters:");
29              Console.WriteLine("  Projection:" + ds.GetProjectionRef());
30              Console.WriteLine("  Rastercount:" + ds.RasterCount);//RasterCount是波段数
31              Console.WriteLine("  RasterSize (" + ds.RasterXSize + "," + ds.RasterYSize + ")");
32
33              /************************************************************************/
34              /* Get Driver                                                                     */
35              /************************************************************************/
36              Driver drv = ds.GetDriver();
37              if (drv ==null)
38              {
39                  Console.WriteLine("Can't get driver");
40                  System.Environment.Exit(-1);
41              }
42              Console.WriteLine("using driver" + drv.LongName);
43              /************************************************************************/
44              /* Get raster band                                                                    */
45              /************************************************************************/
46              for (int iBand=1;iBand<=ds.RasterCount;iBand++)
47              {
48                  Band band = ds.GetRasterBand(iBand);
49                  Console.WriteLine("Band" + iBand + ":");
50                  Console.WriteLine("   DataType:" + band.DataType);
51                  Console.WriteLine("   Size (" + band.XSize + "," + band.YSize + ")");
52                 Console.WriteLine("   PaletteInterp: " + band.GetRasterColorInterpretation().ToString()); //调色说明？
53
54                 for (int iOver = 0; iOver < band.GetOverviewCount(); iOver++)
55                 {
56                     Band over = band.GetOverview(iOver);
57                     Console.WriteLine("      OverView " + iOver + " :");
58                     Console.WriteLine("         DataType: " + over.DataType);
59                     Console.WriteLine("         Size (" + over.XSize + "," + over.YSize + ")");
60                     Console.WriteLine("         PaletteInterp: " + over.GetRasterColorInterpretation().ToString());
61                 }
62
63              }
64              /************************************************************************/
65              /* Processing the raster
66               * To be continued
67              /************************************************************************/
68
69          }
70
71          catch (System.Exception ex)
72          {
73              Console.WriteLine("Application error: " + ex.Message);
74          }
75          Console.ReadLine();
76         }
77     }
78 }

六、运行结果

我找了一幅.img的遥感影像，不含投影坐标，运行结果如下：

转载于:https://www.cnblogs.com/ming5536/archive/2011/05/30/ming5536.html

C# GDAL 学习一相关推荐

GDAL学习笔记02：GDAL基础知识
你的习惯决定了你会成为什么样的人. GDAL学习笔记02:GDAL基础知识前言 1. 版本 2. 摘要 3. 说明 4. 微信公众号GISRSGeography 一.GDAL简介二.导入GDAL ...
（C++）GDAL学习笔记—— 5 全色影像与多光谱影像的融合
任务将全色影像与多光谱影像融合. 原理使用的方法是HSI变换法,即先将多光谱影像选择RGB三个波段进行HSI变化得到H(色调).S(饱和度).I(亮度)三个波段.之后按照全色波段的尺寸将H和S两个 ...
Python GDAL学习笔记(一)
目录一.读取tif格式影像导入GDAL模块并查看版本/位置打开影像查看行列号和波段数查看影像文件描述和其他属性打开数据集波段并获取统计信息利用numpy将波段数据写入数组将波段数据写入 ...
【GDAL学习】用GDAL读取栅格数据
1.根据坐标读取遥感影像的单个像素值 # week 4: get pixel values at a set of coordinates by reading in one pixel at a t ...
GDAL学习笔记——OGR投影
解析坐标系当OGRSpatialReference对象被创建时,就可以解析该对象所包含的信息.可以使用OGRSpatialReference提供的IsProjected()和IsGeographic ...
Python+GDAL几何校正任意自带经纬度数据的遥感影像
目录前言构建虚拟数据集几何校正完整代码前言李民录老师在他的博客中使用C++ GDAL的gdaltranslate.exe.gdalwarp.exe等工具对FY3A_MERSI数据进行了几何 ...
windows下gdal的java开发环境搭建
本文接前文gis利器之Gdal(一),前面简单介绍了gdal工具是什么,具体有什么功能.对于gdal的具体功能,感兴趣的去官网看看.日常对于gis常见的矢量和栅格数据进行处理,gdal都是拿手的.今天 ...
iOS和Android上的本机心理健康监控：现在不是时候了吗？
Vinith Johnson and Girish Subramanyan Vinith Johnson和 Girish Subramanyan Imagine for a minute that i ...
高通芯片为什么输苹果_一位软件工程师解释了为什么苹果要转向自己的芯片
高通芯片为什么输苹果 Apple's Worldwide Developer Conference, which kicked off at the end of June, brought a sl ...

C# GDAL 学习一

C# GDAL 学习一相关推荐

最新文章

热门文章