Method GDALMDArrayFromDataset

gcore/gdaldataset.cpp:11863–12156

  protected:
    GDALMDArrayFromDataset(GDALDataset *poDS, CSLConstList papszOptions)
        : GDALAbstractMDArray(std::string(),
                              std::string(poDS->GetDescription())),
          GDALMDArray(std::string(), std::string(poDS->GetDescription())),
          m_poDS(poDS), m_dt(GDALExtendedDataType::Create(
                            poDS->GetRasterBand(1)->GetRasterDataType())),
          m_osFilename(poDS->GetDescription()), m_aosOptions(papszOptions)
    {
        m_poDS->Reference();

        const int nBandCount = poDS->GetRasterCount();
        for (int i = 1; i <= nBandCount; ++i)
        {
            const auto poBand = poDS->GetRasterBand(i);
            if (i == 1)
                m_osUnit = poBand->GetUnitType();
            else if (m_osUnit != poBand->GetUnitType())
                m_osUnit.clear();

            std::vector<GByte> abyNoData;
            int bHasNoData = false;
            switch (poBand->GetRasterDataType())
            {
                case GDT_Int64:
                {
                    const auto nNoData =
                        poBand->GetNoDataValueAsInt64(&bHasNoData);
                    if (bHasNoData)
                    {
                        abyNoData.resize(m_dt.GetSize());
                        GDALCopyWords64(&nNoData, GDT_Int64, 0, &abyNoData[0],
                                        m_dt.GetNumericDataType(), 0, 1);
                    }
                    break;
                }

                case GDT_UInt64:
                {
                    const auto nNoData =
                        poBand->GetNoDataValueAsUInt64(&bHasNoData);
                    if (bHasNoData)
                    {
                        abyNoData.resize(m_dt.GetSize());
                        GDALCopyWords64(&nNoData, GDT_UInt64, 0, &abyNoData[0],
                                        m_dt.GetNumericDataType(), 0, 1);
                    }
                    break;
                }

                default:
                {
                    const auto dfNoData = poBand->GetNoDataValue(&bHasNoData);
                    if (bHasNoData)
                    {
                        abyNoData.resize(m_dt.GetSize());
                        GDALCopyWords64(&dfNoData, GDT_Float64, 0,
                                        &abyNoData[0],
                                        m_dt.GetNumericDataType(), 0, 1);