PIE SDK打开静止卫星数据

/// <summary>
/// 加载静止卫星数据
/// </summary>
/// <param name="sender"></param>
/// <param name="e"></param>
private void toolStripButton_AddStatic_Click(object sender, EventArgs e)
{
OpenFileDialog openFile = new OpenFileDialog();
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);
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;
//打开MultiDataset
IMultiDataset hdfDataset = PIE.DataSource.DatasetFactory.OpenDataset(filePath, OpenMode.ReadOnly) as IMultiDataset;
if (hdfDataset == null) return null;
//遍历，查找指定通道的Dataset，进行数据格式转换
for (int i = ; i < hdfDataset.GetDatasetCount(); i++)
{
//1、获取操作数据
IDataset pTempDataset = hdfDataset.GetDataset(i);
if (pTempDataset.Name != channelName) continue;
IRasterDataset hdfRasterDatasetBand = pTempDataset as IRasterDataset; //2、读写栅格数据形成新的栅格数据集
int nWidth = hdfRasterDatasetBand.GetRasterXSize();
int nHeight = hdfRasterDatasetBand.GetRasterYSize();
PixelDataType pixDataType = hdfRasterDatasetBand.GetRasterBand().GetRasterDataType();
int bandCount = hdfRasterDatasetBand.GetBandCount();
int[] bandMap = new int[bandCount];
for (int j = ; j < bandCount; j++)
{
bandMap[j] = j + ;
}
UInt16[] arr = new UInt16[nWidth * nHeight * bandCount];
bool IsOk = hdfRasterDatasetBand.Read(, , nWidth, nHeight, arr, nWidth, nHeight, pixDataType, bandCount, bandMap); IRasterDataset newRasterDataset = DatasetFactory.CreateRasterDataset(tiffPath, nWidth, nHeight, bandCount, pixDataType, "GTiff", null);
bool flag = newRasterDataset.Write(, , nWidth, nHeight, arr, nWidth, nHeight, pixDataType, bandCount, bandMap);
newRasterDataset.SpatialReference = spatialReference;
newRasterDataset.GetRasterBand().SetNoDataValue(); //六参数，根据输入坐标的不同需要进行动态设置，本示例代码以风云4A-4000m的数据作为实验数据。
int beginLineNum = ;
int nReslution = ;
double[] geoTransform = new double[];
geoTransform[] = -;
geoTransform[] = nReslution;
geoTransform[] = ;
geoTransform[] = - beginLineNum * nReslution;
geoTransform[] = ;
geoTransform[] = -nReslution;
newRasterDataset.SetGeoTransform(geoTransform); (newRasterDataset as IDisposable).Dispose();
(hdfRasterDatasetBand as IDisposable).Dispose();
rasteLayer = PIE.Carto.LayerFactory.CreateDefaultLayer(tiffPath) as IRasterLayer;
break;
}
return rasteLayer;
}