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

五、完整代码

 

 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的遥感影像，不含投影坐标，运行结果如下：