PIE-SDK For C++打开静止卫星数据

第一步

打开静止卫星数据为多数据集

第二步

获取指定通道的栅格数据集

第三步

读取第二步中的数据集的数据至内存中

第四步

创建与静止卫星同数据类型、同宽高、同波段数的目标栅格文件

第五步

将数据写入目标栅格数据文件

第六步

对目标栅格数据赋值空间参考和六参数

接口/类

方法

说明

SysDataSource::DatasetFactory

OpenDataset

打开数据集

CreateRasterDataset

创建栅格数据集

SysDataSource::RasterDatasetPtr

Read

将栅格数据读取至内存中

Write

将内存数据写入至栅格数据中

项目路径

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

数据路径

百度云盘地址下/PIE示例数据/栅格数据/00.FY/FY4A/**.tif

视频路径

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

示例代码

//加载静止卫星数据

void PIEMainWindow::On_OpenStaticData_Triggered(bool checked)

{

    QString filter = "HDF Files (*.hdf *.h5)";

    QString lstFile = QFileDialog::getOpenFileName(nullptr, "添加数据", "", filter);

    if (lstFile.isEmpty()) return;

    QFileInfo fileInfo(lstFile);

    QString desDir = fileInfo.absoluteDir().absolutePath();

    QString fileName = fileInfo.baseName();

    QString channelName = "NOMChannel13";

    QString desTif = desDir + "/" + fileName+"_" +channelName+ ".tiff";

    QFileInfo desTifInfo(desTif);

    if (desTifInfo.exists()) return;

    SysGeometry::SpatialReferencePtr spatialReference = new SysGeometry::ProjectedCoordinateSystem();

    QString prj4 = "+proj=geos +h=35785863 +a=6378137.0 +b=6356752.3 +lon_0=104.7 +no_defs";

    spatialReference->ImportFromProj4(prj4);

    SysDataSource::MultiDatasetPtr mulDataSetPtr = SysDataSource::DatasetFactory::Instance()->OpenDataset(lstFile, SysDataSource::GA_ReadOnly);

    if (mulDataSetPtr == nullptr) return;

    SysDataSource::RasterDatasetPtr channelDataset= mulDataSetPtr->GetDataset(channelName);

    if (channelDataset!=nullptr)

    {

        //1、获取栅格数据集的相关参数

        int nWidth = channelDataset->GetRasterXSize();

        int nHeight = channelDataset->GetRasterYSize();

        SysDataSource::PixelDataType pixDataType = channelDataset->GetRasterBand(0)->GetRasterDataType();

        int bandCount = channelDataset->GetBandCount();

        QVector<int> bandMap;

        int* bandMapTarget = new int[bandCount];

        for (int i=0;i<bandCount;i++)

        {

            bandMap.insert(i, i + 1);

            bandMapTarget[i] = i + 1;

        }

        //2、读取数据至内存中

        SysDataSource::PixelBufferPtr bufferPtr = channelDataset->Read(0, 0, nWidth, nHeight, nWidth, nHeight, bandMap);

        //3、创建目标栅格数据并写入

        SysDataSource::RasterDatasetPtr targetDataSetPtr = SysDataSource::DatasetFactory::Instance()->CreateRasterDataset(desTif, nWidth, nHeight, bandCount,

            pixDataType, "GTiff", nullptr);

        bool flag= targetDataSetPtr->Write(0, 0, nWidth, nHeight, bufferPtr->GetData(), nWidth, nHeight, pixDataType, bandCount, bandMapTarget);

        if (flag)

        {

            //设置空间投影并设置无效值

            targetDataSetPtr->SetSpatialReference(spatialReference);

            targetDataSetPtr->GetRasterBand(0)->SetNoDataValue(65535);

            //设置六参数

            double* geoTransform = new double[6];

            int beginLineNum = 0;

            int nReslution = 4000;

            geoTransform[0] = -5496000;

            geoTransform[1] = nReslution;

            geoTransform[2] = 0;

            geoTransform[3] = 5496000 - beginLineNum * nReslution;

            geoTransform[4] = 0;

            geoTransform[5] = -nReslution;

            targetDataSetPtr->SetGeoTransform(geoTransform);

            targetDataSetPtr->Flush();

            SysCarto::LayerPtr channelLayer = SysCarto::LayerFactory::Instance()->CreateDefaultRasterLayer(targetDataSetPtr);

            m_pCurrentControl->GetMap()->AddLayer(channelLayer);

            m_pCurrentControl->GetActiveView()->PartialRefresh(SysCarto::ViewAll);

        }

    }

}