/****************************************************************************** * * Project: OpenGIS Simple Features Reference Implementation * Purpose: C API to read a MiraMon layer * Author: Abel Pau, a.pau@creaf.uab.cat, based on the MiraMon codes, * mainly written by Xavier Pons, Joan Maso (correctly written * "Mas0xF3"), Abel Pau, Nuria Julia (N0xFAria Juli0xE0), * Xavier Calaf, Lluis (Llu0xEDs) Pesquer and Alaitz Zabala, from * CREAF and Universitat Autonoma (Aut0xF2noma) de Barcelona. * For a complete list of contributors: * https://www.miramon.cat/eng/QuiSom.htm ****************************************************************************** * Copyright (c) 2024, Xavier Pons * * SPDX-License-Identifier: MIT ****************************************************************************/ #ifdef GDAL_COMPILATION #include "ogr_api.h" // For CPL_C_START #include "mm_wrlayr.h" #include "mm_wrlayr.h" // For MMReadHeader() #include "mm_gdal_functions.h" #include "mm_gdal_constants.h" #else #include "CmptCmp.h" // Compatibility between compilers #include "mm_gdal\mm_wrlayr.h" // For MMReadHeader() #include "mm_gdal\mm_gdal_functions.h" // For int MM_GetArcHeights() #include "mm_constants.h" #endif #include "mm_rdlayr.h" #ifdef GDAL_COMPILATION CPL_C_START // Necessary for compiling in GDAL project #endif /* -------------------------------------------------------------------- */ /* Reading MiraMon format file functions */ /* -------------------------------------------------------------------- */ // Initializes a MiraMon vector layer for reading int MMInitLayerToRead(struct MiraMonVectLayerInfo *hMiraMonLayer, FILE_TYPE *m_fp, const char *pszFilename) { char szResult[MM_MAX_ID_SNY + 10]; char *pszSRS; memset(hMiraMonLayer, 0, sizeof(*hMiraMonLayer)); if (MMReadHeader(m_fp, &hMiraMonLayer->TopHeader)) { MMCPLError(CE_Failure, CPLE_NoWriteAccess, "Error reading header of file %s", pszFilename); return 1; } hMiraMonLayer->ReadOrWrite = MM_READING_MODE; strcpy(hMiraMonLayer->pszFlags, "rb"); hMiraMonLayer->pszSrcLayerName = strdup_function(pszFilename); hMiraMonLayer->LayerVersion = (char)MMGetVectorVersion(&hMiraMonLayer->TopHeader); if (hMiraMonLayer->LayerVersion == MM_UNKNOWN_VERSION) { MMCPLError(CE_Failure, CPLE_NotSupported, "MiraMon version file unknown."); return 1; } if (hMiraMonLayer->LayerVersion == MM_LAST_VERSION) hMiraMonLayer->nHeaderDiskSize = MM_HEADER_SIZE_64_BITS; else if (hMiraMonLayer->LayerVersion == MM_32BITS_VERSION) hMiraMonLayer->nHeaderDiskSize = MM_HEADER_SIZE_32_BITS; else hMiraMonLayer->nHeaderDiskSize = MM_HEADER_SIZE_64_BITS; if (hMiraMonLayer->TopHeader.aFileType[0] == 'P' && hMiraMonLayer->TopHeader.aFileType[1] == 'N' && hMiraMonLayer->TopHeader.aFileType[2] == 'T') { if (hMiraMonLayer->TopHeader.Flag & MM_LAYER_3D_INFO) { hMiraMonLayer->TopHeader.bIs3d = 1; hMiraMonLayer->eLT = MM_LayerType_Point3d; } else hMiraMonLayer->eLT = MM_LayerType_Point; hMiraMonLayer->bIsPoint = TRUE; } else if (hMiraMonLayer->TopHeader.aFileType[0] == 'A' && hMiraMonLayer->TopHeader.aFileType[1] == 'R' && hMiraMonLayer->TopHeader.aFileType[2] == 'C') { if (hMiraMonLayer->TopHeader.Flag & MM_LAYER_3D_INFO) { hMiraMonLayer->TopHeader.bIs3d = 1; hMiraMonLayer->eLT = MM_LayerType_Arc3d; } else hMiraMonLayer->eLT = MM_LayerType_Arc; hMiraMonLayer->bIsArc = TRUE; } else if (hMiraMonLayer->TopHeader.aFileType[0] == 'P' && hMiraMonLayer->TopHeader.aFileType[1] == 'O' && hMiraMonLayer->TopHeader.aFileType[2] == 'L') { // 3D if (hMiraMonLayer->TopHeader.Flag & MM_LAYER_3D_INFO) { hMiraMonLayer->TopHeader.bIs3d = 1; hMiraMonLayer->eLT = MM_LayerType_Pol3d; } else hMiraMonLayer->eLT = MM_LayerType_Pol; hMiraMonLayer->bIsPolygon = TRUE; if (hMiraMonLayer->TopHeader.Flag & MM_LAYER_MULTIPOLYGON) hMiraMonLayer->TopHeader.bIsMultipolygon = 1; } //hMiraMonLayer->Version = MM_VECTOR_LAYER_LAST_VERSION; if (MMInitLayerByType(hMiraMonLayer)) return 1; hMiraMonLayer->bIsBeenInit = 1; // Get the basic metadata pszSRS = MMReturnValueFromSectionINIFile( hMiraMonLayer->pszMainREL_LayerName, "SPATIAL_REFERENCE_SYSTEM:HORIZONTAL", "HorizontalSystemIdentifier"); if (pszSRS) hMiraMonLayer->pSRS = pszSRS; else hMiraMonLayer->pSRS = nullptr; if (!hMiraMonLayer->pSRS && hMiraMonLayer->bIsPolygon) { pszSRS = MMReturnValueFromSectionINIFile( hMiraMonLayer->MMPolygon.MMArc.pszREL_LayerName, "SPATIAL_REFERENCE_SYSTEM:HORIZONTAL", "HorizontalSystemIdentifier"); hMiraMonLayer->pSRS = pszSRS; } if (!ReturnEPSGCodeSRSFromMMIDSRS(hMiraMonLayer->pSRS, szResult)) { if (MMIsEmptyString(szResult)) hMiraMonLayer->nSRS_EPSG = 0; else hMiraMonLayer->nSRS_EPSG = atoi(szResult); } else hMiraMonLayer->nSRS_EPSG = 0; if (hMiraMonLayer->nSRS_EPSG == 0) { if (hMiraMonLayer->pSRS && strcmp(hMiraMonLayer->pSRS, "plane")) { MMCPLWarning(CE_Warning, CPLE_NotSupported, "The MiraMon layer SRS has no equivalent " "in EPSG code"); } } // If more nNumStringToOperate is needed, it'll be increased. hMiraMonLayer->nNumStringToOperate = 0; if (MMResizeStringToOperateIfNeeded(hMiraMonLayer, 5000)) return 1; return 0; } // Reads stringline coordinates and puts them in a buffer static int MMAddStringLineCoordinates(struct MiraMonVectLayerInfo *hMiraMonLayer, MM_INTERNAL_FID i_elem, uint32_t flag_z, MM_N_VERTICES_TYPE nStartVertice, MM_BOOLEAN bAvoidFirst, unsigned char VFG) { FILE_TYPE *pF; struct MM_AH *pArcHeader; struct MiraMonArcLayer *pMMArc; struct MM_ZD *pZDescription = nullptr; if (hMiraMonLayer->bIsPolygon) pMMArc = &hMiraMonLayer->MMPolygon.MMArc; else pMMArc = &hMiraMonLayer->MMArc; pF = pMMArc->pF; pArcHeader = pMMArc->pArcHeader; if (hMiraMonLayer->TopHeader.bIs3d) pZDescription = pMMArc->pZSection.pZDescription; fseek_function(pF, pArcHeader[i_elem].nOffset, SEEK_SET); if (hMiraMonLayer->bIsPolygon && (VFG & MM_POL_REVERSE_ARC)) // && //nStartVertice > 0) { MM_N_VERTICES_TYPE nIVertice; // Reading arcs vertices in an inverse order if (MMResizeMM_POINT2DPointer( &hMiraMonLayer->ReadFeature.pCoord, &hMiraMonLayer->ReadFeature.nMaxpCoord, nStartVertice + pArcHeader[i_elem].nElemCount * 2, // ask for twice memory to reverse 0, 0)) return 1; // Get the vertices far away from their place to be inverted later if (pArcHeader[i_elem].nElemCount != fread_function(hMiraMonLayer->ReadFeature.pCoord + nStartVertice + pArcHeader[i_elem].nElemCount, sizeof(*hMiraMonLayer->ReadFeature.pCoord), (size_t)pArcHeader[i_elem].nElemCount, pF)) { return 1; } if (hMiraMonLayer->TopHeader.bIs3d) { if (MMResizeDoublePointer( &hMiraMonLayer->ReadFeature.pZCoord, &hMiraMonLayer->ReadFeature.nMaxpZCoord, nStartVertice + pArcHeader[i_elem].nElemCount * 2, 0, 0)) return 1; // +nStartVertice MM_GetArcHeights(hMiraMonLayer->ReadFeature.pZCoord + nStartVertice + pArcHeader[i_elem].nElemCount, pF, pArcHeader[i_elem].nElemCount, pZDescription + i_elem, flag_z); // If there is a value for Z-nodata in GDAL this lines can be uncommented // and MM_GDAL_NODATA_COORD_Z can be defined /*if(!DOUBLES_DIFERENTS_DJ(punts_z[k], MM_NODATA_COORD_Z)) { MM_N_VERTICES_TYPE nIVertice; for(nIVertice=0; nIVerticeReadFeature.pZCoord[nIVertice]=MM_GDAL_NODATA_COORD_Z; } */ } // Reverse the vertices while putting on their place for (nIVertice = 0; nIVertice < pArcHeader[i_elem].nElemCount; nIVertice++) { memcpy(hMiraMonLayer->ReadFeature.pCoord + nStartVertice - ((nStartVertice > 0 && bAvoidFirst) ? 1 : 0) + nIVertice, hMiraMonLayer->ReadFeature.pCoord + nStartVertice + 2 * pArcHeader[i_elem].nElemCount - nIVertice - 1, sizeof(*hMiraMonLayer->ReadFeature.pCoord)); if (hMiraMonLayer->TopHeader.bIs3d) { memcpy(hMiraMonLayer->ReadFeature.pZCoord + nStartVertice - ((nStartVertice > 0 && bAvoidFirst) ? 1 : 0) + nIVertice, hMiraMonLayer->ReadFeature.pZCoord + nStartVertice + 2 * pArcHeader[i_elem].nElemCount - nIVertice - 1, sizeof(*hMiraMonLayer->ReadFeature.pZCoord)); } } } else { // Reading arcs vertices if (MMResizeMM_POINT2DPointer( &hMiraMonLayer->ReadFeature.pCoord, &hMiraMonLayer->ReadFeature.nMaxpCoord, nStartVertice + pArcHeader[i_elem].nElemCount, 0, 0)) return 1; if (pArcHeader[i_elem].nElemCount != fread_function(hMiraMonLayer->ReadFeature.pCoord + nStartVertice - (bAvoidFirst ? 1 : 0), sizeof(*hMiraMonLayer->ReadFeature.pCoord), (size_t)pArcHeader[i_elem].nElemCount, pF)) { return 1; } if (hMiraMonLayer->TopHeader.bIs3d) { if (MMResizeDoublePointer( &hMiraMonLayer->ReadFeature.pZCoord, &hMiraMonLayer->ReadFeature.nMaxpZCoord, nStartVertice + pArcHeader[i_elem].nElemCount, 0, 0)) return 1; // +nStartVertice MM_GetArcHeights(hMiraMonLayer->ReadFeature.pZCoord + nStartVertice - (bAvoidFirst ? 1 : 0), pF, pArcHeader[i_elem].nElemCount, pZDescription + i_elem, flag_z); // If there is a value for Z-nodata in GDAL this lines can be uncommented // and MM_GDAL_NODATA_COORD_Z can be defined /*if(!DOUBLES_DIFERENTS_DJ(punts_z[k], MM_NODATA_COORD_Z)) { MM_N_VERTICES_TYPE nIVertice; for(nIVertice=0; nIVerticeReadFeature.pZCoord[nIVertice]=MM_GDAL_NODATA_COORD_Z; } */ } } hMiraMonLayer->ReadFeature.nNumpCoord = pArcHeader[i_elem].nElemCount == 0 ? 0 : pArcHeader[i_elem].nElemCount - (bAvoidFirst ? 1 : 0); return 0; } // Reads Polygon coordinates and puts them in a buffer static int MMGetMultiPolygonCoordinates(struct MiraMonVectLayerInfo *hMiraMonLayer, MM_INTERNAL_FID i_pol, uint32_t flag_z) { struct MM_PH *pPolHeader; struct MM_AH *pArcHeader; char *pBuffer; MM_POLYGON_ARCS_COUNT nIndex; MM_BOOLEAN bAvoidFirst; MM_N_VERTICES_TYPE nNAcumulVertices = 0; // Checking if the index of the polygon is in the correct range. if (i_pol >= hMiraMonLayer->TopHeader.nElemCount) return 1; MMResetFeatureGeometry(&hMiraMonLayer->ReadFeature); MMResetFeatureRecord(&hMiraMonLayer->ReadFeature); pPolHeader = hMiraMonLayer->MMPolygon.pPolHeader + i_pol; // It's accepted not having arcs in the universal polygon if (!pPolHeader->nArcsCount) { if (i_pol == 0) return 0; else return 1; } if (MMResizeMiraMonPolygonArcs(&hMiraMonLayer->pArcs, &hMiraMonLayer->nMaxArcs, pPolHeader->nArcsCount, 0, 0)) return 1; if (MMInitFlush(&hMiraMonLayer->FlushPAL, hMiraMonLayer->MMPolygon.pF, hMiraMonLayer->MMPolygon.nPALElementSize * pPolHeader->nArcsCount, &pBuffer, pPolHeader->nOffset, 0)) { if (pBuffer) free_function(pBuffer); return 1; } hMiraMonLayer->FlushPAL.pBlockWhereToSaveOrRead = (void *)pBuffer; if (MMReadFlush(&hMiraMonLayer->FlushPAL)) { if (pBuffer) free_function(pBuffer); return 1; } hMiraMonLayer->ReadFeature.nNRings = 0; hMiraMonLayer->ReadFeature.nNumpCoord = 0; if (MMResize_MM_N_VERTICES_TYPE_Pointer( &hMiraMonLayer->ReadFeature.pNCoordRing, &hMiraMonLayer->ReadFeature.nMaxpNCoordRing, (MM_N_VERTICES_TYPE)hMiraMonLayer->ReadFeature.nNRings + 1, 10, 10)) { free_function(pBuffer); return 1; } if (MMResizeVFGPointer(&hMiraMonLayer->ReadFeature.flag_VFG, &hMiraMonLayer->ReadFeature.nMaxVFG, (MM_INTERNAL_FID)pPolHeader->nArcsCount, 0, 0)) // Perhaps more memory than needed { free_function(pBuffer); return 1; } // Preparing memory for all coordinates hMiraMonLayer->ReadFeature.pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] = 0; for (nIndex = 0; nIndex < pPolHeader->nArcsCount; nIndex++) { hMiraMonLayer->FlushPAL.SizeOfBlockToBeSaved = sizeof((hMiraMonLayer->pArcs + nIndex)->VFG); hMiraMonLayer->FlushPAL.pBlockToBeSaved = (void *)&(hMiraMonLayer->pArcs + nIndex)->VFG; if (MMReadBlockFromBuffer(&hMiraMonLayer->FlushPAL)) { if (pBuffer) free_function(pBuffer); return 1; } // Arc index if (MMReadGUInt64DependingOnVersion( hMiraMonLayer, &hMiraMonLayer->FlushPAL, &((hMiraMonLayer->pArcs + nIndex)->nIArc))) { if (pBuffer) free_function(pBuffer); return 1; } if (hMiraMonLayer->MMPolygon.MMArc.pArcHeader == nullptr) { if (pBuffer) free_function(pBuffer); return 1; } // Checking if the index of the arc is in the correct range. if ((hMiraMonLayer->pArcs + nIndex)->nIArc >= hMiraMonLayer->MMPolygon.TopArcHeader.nElemCount) { free_function(pBuffer); return 1; } pArcHeader = hMiraMonLayer->MMPolygon.MMArc.pArcHeader + (hMiraMonLayer->pArcs + nIndex)->nIArc; if (hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] > UINT64_MAX - pArcHeader->nElemCount) { free_function(pBuffer); return 1; } if (hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] > UINT64_MAX - pArcHeader->nElemCount) { free_function(pBuffer); return 1; } hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] += pArcHeader->nElemCount; } if (MMResizeMM_POINT2DPointer( &hMiraMonLayer->ReadFeature.pCoord, &hMiraMonLayer->ReadFeature.nMaxpCoord, hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings], 0, 0)) { free_function(pBuffer); return 1; } hMiraMonLayer->FlushPAL.CurrentOffset = 0; // Real work hMiraMonLayer->ReadFeature.pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] = 0; for (nIndex = 0; nIndex < pPolHeader->nArcsCount; nIndex++) { hMiraMonLayer->FlushPAL.SizeOfBlockToBeSaved = sizeof((hMiraMonLayer->pArcs + nIndex)->VFG); hMiraMonLayer->FlushPAL.pBlockToBeSaved = (void *)&(hMiraMonLayer->pArcs + nIndex)->VFG; if (MMReadBlockFromBuffer(&hMiraMonLayer->FlushPAL)) { if (pBuffer) free_function(pBuffer); return 1; } // Arc index if (MMReadGUInt64DependingOnVersion( hMiraMonLayer, &hMiraMonLayer->FlushPAL, &((hMiraMonLayer->pArcs + nIndex)->nIArc))) { if (pBuffer) free_function(pBuffer); return 1; } bAvoidFirst = FALSE; if (hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] != 0) bAvoidFirst = TRUE; // Add coordinates to hMiraMonLayer->ReadFeature.pCoord if (MMAddStringLineCoordinates(hMiraMonLayer, (hMiraMonLayer->pArcs + nIndex)->nIArc, flag_z, nNAcumulVertices, bAvoidFirst, (hMiraMonLayer->pArcs + nIndex)->VFG)) { free_function(pBuffer); return 1; } if (MMResize_MM_N_VERTICES_TYPE_Pointer( &hMiraMonLayer->ReadFeature.pNCoordRing, &hMiraMonLayer->ReadFeature.nMaxpNCoordRing, (MM_N_VERTICES_TYPE)hMiraMonLayer->ReadFeature.nNRings + 1, 10, 10)) { free_function(pBuffer); return 1; } if (hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] > UINT64_MAX - hMiraMonLayer->ReadFeature.nNumpCoord) { free_function(pBuffer); return 1; } hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] += hMiraMonLayer->ReadFeature.nNumpCoord; nNAcumulVertices += hMiraMonLayer->ReadFeature.nNumpCoord; if ((hMiraMonLayer->pArcs + nIndex)->VFG & MM_POL_END_RING) { hMiraMonLayer->ReadFeature .flag_VFG[hMiraMonLayer->ReadFeature.nNRings] = (hMiraMonLayer->pArcs + nIndex)->VFG; hMiraMonLayer->ReadFeature.nNRings++; hMiraMonLayer->ReadFeature .pNCoordRing[hMiraMonLayer->ReadFeature.nNRings] = 0; } } hMiraMonLayer->nNumArcs = pPolHeader->nArcsCount; if (pBuffer) free_function(pBuffer); return 0; } // Reads the geographical part of a MiraMon layer feature int MMGetGeoFeatureFromVector(struct MiraMonVectLayerInfo *hMiraMonLayer, MM_INTERNAL_FID i_elem) { FILE_TYPE *pF; struct MM_ZD *pZDescription; uint32_t flag_z; int num; double cz; if (hMiraMonLayer->nSelectCoordz == MM_SELECT_HIGHEST_COORDZ) flag_z = MM_STRING_HIGHEST_ALTITUDE; else if (hMiraMonLayer->nSelectCoordz == MM_SELECT_LOWEST_COORDZ) flag_z = MM_STRING_LOWEST_ALTITUDE; else flag_z = 0L; if (hMiraMonLayer->bIsPoint) { pF = hMiraMonLayer->MMPoint.pF; // Getting to the i-th element offset fseek_function(pF, hMiraMonLayer->nHeaderDiskSize + sizeof(MM_COORD_TYPE) * 2 * i_elem, SEEK_SET); // Reading the point if (MMResizeMM_POINT2DPointer(&hMiraMonLayer->ReadFeature.pCoord, &hMiraMonLayer->ReadFeature.nMaxpCoord, hMiraMonLayer->ReadFeature.nNumpCoord, 1, 1)) return 1; if (1 != fread_function(hMiraMonLayer->ReadFeature.pCoord, sizeof(MM_COORD_TYPE) * 2, 1, pF)) { return 1; } hMiraMonLayer->ReadFeature.nNRings = 1; if (MMResize_MM_N_VERTICES_TYPE_Pointer( &hMiraMonLayer->ReadFeature.pNCoordRing, &hMiraMonLayer->ReadFeature.nMaxpNCoordRing, 1, 0, 1)) return 1; hMiraMonLayer->ReadFeature.pNCoordRing[0] = 1; if (hMiraMonLayer->TopHeader.bIs3d) { pZDescription = hMiraMonLayer->MMPoint.pZSection.pZDescription + i_elem; if (pZDescription->nZCount == INT_MIN) return 1; num = MM_ARC_TOTAL_N_HEIGHTS_DISK(pZDescription->nZCount, 1); if (MMResizeDoublePointer(&hMiraMonLayer->ReadFeature.pZCoord, &hMiraMonLayer->ReadFeature.nMaxpZCoord, 1, 1, 1)) return 1; if (num == 0) hMiraMonLayer->ReadFeature.pZCoord[0] = MM_NODATA_COORD_Z; else { if (flag_z == MM_STRING_HIGHEST_ALTITUDE) // Max z cz = pZDescription->dfBBmaxz; else if (flag_z == MM_STRING_LOWEST_ALTITUDE) // Min z cz = pZDescription->dfBBminz; else { // Reading the first z coordinate fseek_function(pF, pZDescription->nOffsetZ, SEEK_SET); if ((size_t)1 != fread_function( &cz, sizeof(*hMiraMonLayer->ReadFeature.pZCoord), 1, pF)) { return 1; } } // If there is a value for Z-nodata in GDAL this lines can be uncommented // and MM_GDAL_NODATA_COORD_Z can be defined /*if(!DOUBLES_DIFERENTS_DJ(cz, MM_NODATA_COORD_Z)) hMiraMonLayer->ReadFeature.pZCoord[0]=MM_GDAL_NODATA_COORD_Z; else */ hMiraMonLayer->ReadFeature.pZCoord[0] = cz; } } return 0; } // Stringlines if (hMiraMonLayer->bIsArc && !hMiraMonLayer->bIsPolygon) { if (MMAddStringLineCoordinates(hMiraMonLayer, i_elem, flag_z, 0, FALSE, 0)) return 1; if (MMResize_MM_N_VERTICES_TYPE_Pointer( &hMiraMonLayer->ReadFeature.pNCoordRing, &hMiraMonLayer->ReadFeature.nMaxpNCoordRing, 1, 0, 1)) return 1; hMiraMonLayer->ReadFeature.pNCoordRing[0] = hMiraMonLayer->ReadFeature.nNumpCoord; return 0; } // Polygons or multipolygons if (MMGetMultiPolygonCoordinates(hMiraMonLayer, i_elem, flag_z)) return 1; return 0; } // Reads the header of a MiraMon DBF // Please read the format at this link: // https://www.miramon.cat/new_note/usa/notes/DBF_estesa.pdf int MM_ReadExtendedDBFHeader(struct MiraMonVectLayerInfo *hMiraMonLayer) { const char *pszRelFile = nullptr; struct MM_DATA_BASE_XP *pMMBDXP; const char *szDBFFileName = nullptr; // If read don't read again. It happens when Polygon reads // the database and then in initArc() it's read again. if (hMiraMonLayer->pMMBDXP) return 0; pMMBDXP = hMiraMonLayer->pMMBDXP = calloc_function(sizeof(*pMMBDXP)); if (!pMMBDXP) return 1; if (hMiraMonLayer->bIsPoint) { hMiraMonLayer->MMPoint.MMAdmDB.pMMBDXP = pMMBDXP; szDBFFileName = hMiraMonLayer->MMPoint.MMAdmDB.pszExtDBFLayerName; pszRelFile = hMiraMonLayer->MMPoint.pszREL_LayerName; } else if (hMiraMonLayer->bIsArc && !hMiraMonLayer->bIsPolygon) { hMiraMonLayer->MMArc.MMAdmDB.pMMBDXP = pMMBDXP; szDBFFileName = hMiraMonLayer->MMArc.MMAdmDB.pszExtDBFLayerName; pszRelFile = hMiraMonLayer->MMArc.pszREL_LayerName; } else if (hMiraMonLayer->bIsPolygon) { hMiraMonLayer->MMPolygon.MMAdmDB.pMMBDXP = pMMBDXP; szDBFFileName = hMiraMonLayer->MMPolygon.MMAdmDB.pszExtDBFLayerName; pszRelFile = hMiraMonLayer->MMPolygon.pszREL_LayerName; } if (MM_ReadExtendedDBFHeaderFromFile(szDBFFileName, pMMBDXP, pszRelFile)) { MMCPLError(CE_Failure, CPLE_NotSupported, "Error reading the format in the DBF file %s.", szDBFFileName); return 1; } fclose_and_nullify(&pMMBDXP->pfDataBase); return 0; } #ifdef GDAL_COMPILATION CPL_C_END // Necessary for compiling in GDAL project #endif