PIE SDK打开静止卫星数据

项目路径

百度云盘地址下/PIE示例程序/05.打开静止卫星数据

数据路径

百度云盘地址下/PIE示例数据/科学数据集\HDF文件夹下数据

视频路径

百度云盘地址下/PIE视频教程/05.打开静止卫星数据.avi

示例代码

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);// 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);