PIE SDK打开静止卫星数据

openFile.Filter = "HDF数据|*.hdf";

if (openFile.ShowDialog() != DialogResult.OK)

return;

string channelName = "NOMChannel13";//波段名称

string tempTif = System.IO.Path.GetDirectoryName(openFile.FileName) + "\\NOMChannel13.tiff";//输出tiff路径

ISpatialReference spatialReference = new ProjectedCoordinateSystem();//目标空间参考

spatialReference.ImportFromUserInput("+proj=geos +h=35785863 +a=6378137.0 +b=6356752.3 +lon_0=104.7 +no_defs");

IRasterLayer rasterLayer = OpenStaticData(openFile.FileName, channelName, tempTif, spatialReference);// OpenStaticData方法的定义在下面

if (rasterLayer == null)

return;

//添加至地图并刷新视图

mapControlMain.FocusMap.AddLayer(rasterLayer as ILayer);

mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);

/// <summary>

/// 打开风云4A、风云2G等静止卫星数据，读取指定波段数据为tiff

/// </summary>

/// <param name="filePath">hdf路径</param>

/// <param name="channelName">波段通道名称</param>

/// <param name="tiffPath">生成tiff路径</param>

/// <param name="spatialReference">空间参考</param>

/// <returns></returns>

private IRasterLayer OpenStaticData(string filePath, string channelName, string tiffPath, ISpatialReference spatialReference)

{

IRasterLayer rasteLayer = null;

try

{

//打开MultiDataset

IMultiDataset hdfDataset = PIE.DataSource.DatasetFactory.OpenDataset(filePath, OpenMode.ReadOnly) as IMultiDataset;

if (hdfDataset != null)

{

//遍历，查找指定通道的Dataset，进行数据格式转换

for (int i = 0; i < hdfDataset.GetDatasetCount(); i++)

{

IDataset pTempDataset = hdfDataset.GetDataset(i);

if (pTempDataset.Name == channelName)

{

//将RasterDataset写入指定tiff

IRasterDataset hdfRasterDatasetBand = pTempDataset as IRasterDataset;

int nWidth = hdfRasterDatasetBand.GetRasterXSize();

int nHeight = hdfRasterDatasetBand.GetRasterYSize();

PixelDataType eDateType = hdfRasterDatasetBand.GetRasterBand(0).GetRasterDataType();

int count = hdfRasterDatasetBand.GetBandCount();

int[] bandMap = new int[count];

for (int j = 0; j < count; j++)

{

bandMap[j] = j + 1;

}

string[] tempList = null;

IPixelBuffer pixBuffer = hdfRasterDatasetBand.Read(0, 0, nWidth, nHeight, nWidth, nHeight, bandMap);

//创建输出栅格数据集

IRasterDataset tifRasterDataset = DatasetFactory.CreateRasterDataset(tiffPath, nWidth, nHeight, count, eDateType, "GTiff", tempList);

bool flag = tifRasterDataset.Write(0, 0, nWidth, nHeight, pixBuffer.GetData_Ref(), nWidth, nHeight, eDateType, count, bandMap);

tifRasterDataset.SpatialReference = spatialReference;

tifRasterDataset.GetRasterBand(0).SetNoDataValue(65535);

//六参数，根据输入坐标的不同需要进行动态设置，本示例代码以风云4A-4000m的数据作为实验数据。

int beginLineNum = 0;

int nReslution = 4000;

string beginlineNumStr = hdfDataset.GetMetadataItem("Begin_Line_Number", "");

if (string.IsNullOrEmpty(beginlineNumStr))

{

beginlineNumStr = hdfDataset.GetMetadataItem("geospatial_lat_lon_extent_begin_line_number", "");

if (string.IsNullOrEmpty(beginlineNumStr)) beginLineNum = 183;

}

if (!string.IsNullOrEmpty(beginlineNumStr))

{

StringBuilder sb = new StringBuilder();

foreach (char charStr in beginlineNumStr)

{

if ((charStr >= ‘0‘ && charStr <= ‘9‘) || charStr == ‘ ‘ || charStr == ‘-‘) sb.Append(charStr);

}

bool result = int.TryParse(sb.ToString(), out beginLineNum);

}

double[] geoTransform = new double[6];

geoTransform[0] = -5496000;

geoTransform[1] = nReslution;

geoTransform[2] = 0;

geoTransform[3] = 5496000 - beginLineNum * nReslution;

geoTransform[4] = 0;

geoTransform[5] = -nReslution;

tifRasterDataset.SetGeoTransform(geoTransform);

(tifRasterDataset as IDisposable).Dispose();

(hdfRasterDatasetBand as IDisposable).Dispose();

(pixBuffer as IDisposable).Dispose();

rasteLayer = PIE.Carto.LayerFactory.CreateDefaultLayer(tiffPath) as IRasterLayer;

break;

}

}

}

}

catch (Exception ex)

{

rasteLayer = null;

}

return rasteLayer;

}

mapControlMain.FocusMap.AddLayer(fLayer as ILayer);

mapControlMain.ActiveView.PartialRefresh(PIE.Carto.ViewDrawPhaseType.ViewAll);