1  /// <summary>
 2 /// 加载静止卫星数据
 3 /// </summary>
 4 /// <param name="sender"></param>
 5 /// <param name="e"></param>
 6 private void toolStripButton_AddStatic_Click(object sender, EventArgs e)
 7 {
 8 OpenFileDialog openFile = new OpenFileDialog();
 9 openFile.Filter = "HDF数据|*.hdf";
10 if (openFile.ShowDialog() != DialogResult.OK) return
11
12 string channelName = "NOMChannel13";//波段名称
13 string tempTif = System.IO.Path.GetDirectoryName(openFile.FileName) + "\\NOMChannel13.tiff";//输出tiff路径
14 ISpatialReference spatialReference = new ProjectedCoordinateSystem();//目标空间参考
15 spatialReference.ImportFromUserInput("+proj=geos +h=35785863 +a=6378137.0 +b=6356752.3 +lon_0=104.7 +no_defs");
16
17 IRasterLayer rasterLayer = OpenStaticData(openFile.FileName, channelName, tempTif, spatialReference);
18 if (rasterLayer == null) return;
19 mapControlMain.FocusMap.AddLayer(rasterLayer as ILayer);
20 mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);
21 }
22
23 /// <summary>
24 /// 打开风云4A、风云2G等静止卫星数据，读取指定波段数据为tiff
25 /// </summary>
26 /// <param name="filePath">hdf路径</param>
27 /// <param name="channelName">波段通道名称</param>
28 /// <param name="tiffPath">生成tiff路径</param>
29 /// <param name="spatialReference">空间参考</param>
30 /// <returns></returns>
31 private IRasterLayer OpenStaticData(string filePath, string channelName, string tiffPath, ISpatialReference spatialReference)
32 {
33 IRasterLayer rasteLayer = null;
34 //打开MultiDataset
35 IMultiDataset hdfDataset = PIE.DataSource.DatasetFactory.OpenDataset(filePath, OpenMode.ReadOnly) as IMultiDataset;
36 if (hdfDataset == null) return null;
37 //遍历，查找指定通道的Dataset，进行数据格式转换
38 for (int i = 0; i < hdfDataset.GetDatasetCount(); i++)
39 {
40     //1、获取操作数据
41     IDataset pTempDataset = hdfDataset.GetDataset(i);
42     if (pTempDataset.Name != channelName) continue;
43     IRasterDataset hdfRasterDatasetBand = pTempDataset as IRasterDataset;
44
45     //2、读写栅格数据形成新的栅格数据集
46     int nWidth = hdfRasterDatasetBand.GetRasterXSize();
47     int nHeight = hdfRasterDatasetBand.GetRasterYSize();
48     PixelDataType pixDataType = hdfRasterDatasetBand.GetRasterBand(0).GetRasterDataType();
49     int bandCount = hdfRasterDatasetBand.GetBandCount();
50     int[] bandMap = new int[bandCount];
51     for (int j = 0; j < bandCount; j++)
52     {
53         bandMap[j] = j + 1;
54     }
55     UInt16[] arr = new UInt16[nWidth * nHeight * bandCount];
56     bool IsOk = hdfRasterDatasetBand.Read(0, 0, nWidth, nHeight, arr, nWidth, nHeight, pixDataType, bandCount, bandMap);
57
58     IRasterDataset newRasterDataset = DatasetFactory.CreateRasterDataset(tiffPath, nWidth, nHeight, bandCount, pixDataType, "GTiff", null);
59     bool flag = newRasterDataset.Write(0, 0, nWidth, nHeight, arr, nWidth, nHeight, pixDataType, bandCount, bandMap);
60     newRasterDataset.SpatialReference = spatialReference;
61     newRasterDataset.GetRasterBand(0).SetNoDataValue(65535);
62
63     //六参数，根据输入坐标的不同需要进行动态设置，本示例代码以风云4A-4000m的数据作为实验数据。
64     int beginLineNum = 0;
65     int nReslution = 4000;
66     double[] geoTransform = new double[6];
67     geoTransform[0] = -5496000;
68     geoTransform[1] = nReslution;
69     geoTransform[2] = 0;
70     geoTransform[3] = 5496000 - beginLineNum * nReslution;
71     geoTransform[4] = 0;
72     geoTransform[5] = -nReslution;
73     newRasterDataset.SetGeoTransform(geoTransform);
74
75     (newRasterDataset as IDisposable).Dispose();
76     (hdfRasterDatasetBand as IDisposable).Dispose();
77     rasteLayer = PIE.Carto.LayerFactory.CreateDefaultLayer(tiffPath) as IRasterLayer;
78     break;
79 }
80 return rasteLayer;
81 }