This example shows how to process data from an ENVI cube on disk using the perClass Mira Runtime API. The paths to the model and cube file should be modified to point to existing files.


The example shows how to:

  • Query the version of the acquisition library (line 21-24)
  • Initialize the acquisition plugin (line 26-30)
  • Scan for available computational devices and print their names (line 34-38)
  • Select a computational device (line 40-42)
  • Load a model from disk and print the classifier decisions (line 44-51)
  • Load the correction from disk (line 53-54)
  • Initialize acquisition (line 44-50)
  • Query frame geometry, data type, and layout, and check if it is valid (line 56-70)
  • Load the data from the scan on disk (line 72-85)
  • Enable object segmentation (line 87-89)
  • Process all lines from the cube and print information of any detect object (line 91-122)


In order to compile this example, save it as example.c, and then use the command for the operating system to compile on:


On Windows using the Microsoft Visual C/C++ compiler:

> cl example.c -I "C:\Program Files\perClass Mira\include" "C:\Program Files\perClass Mira\lib\perclass_mira_gpu.lib"


On Linux using GCC:

> gcc example.c -o example -I"/opt/perClass/Mira/include" -L"/opt/perClass/Mira/lib" -lperclass_mira_gpu

  1. #include <stdint.h>
  2. #include <stdio.h>
  3. #include <stdlib.h>

  5. #include "perclass_mira.h"

  7. #define MODEL_PATH "./model.mpl"

  9. #define CUBE_DIR "./data"
  10. #define CUBE_NAME "cube"
  11. #define CUBE_EXT ".pcf"
  12. #define CUBE_WIDTH 640
  13. #define CUBE_BANDS 200
  14. #define CUBE_HEIGHT 1550

  16. int main(int argc, const char* argv[])
  17. {
  18.    int res = MIRA_OK;
  19.    uint16_t* pBuf = NULL;

  21.    printf("Example of processing data using perClass Mira Runtime\n");

  23.    const char* pVersionStr = mira_GetVersion();
  24.    printf("Version: '%s'\n", pVersionStr);

  26.    mrkernel* pmr = mira_Init(".");
  27.    printf("Init: %s", mira_GetErrorMsg(pmr));
  28.    if (pmr == NULL) {
  29.        goto Error;
  30.    }

  32.    MIRA_CHECK(mira_RefreshDeviceList(pmr, 1, 0));

  34.    const int deviceCount = mira_GetDeviceCount(pmr);
  35.    printf("\n%d devices:\n", deviceCount);
  36.    for (int i = 0; i < deviceCount; i++) {
  37.        printf(" %d : %s\n", i, mira_GetDeviceName(pmr, i));
  38.    }

  40.    int deviceInd = 0;
  41.    MIRA_CHECK(mira_SetDevice(pmr, deviceInd));
  42.    printf("Device selected: %d '%s'\n", deviceInd, mira_GetDeviceName(pmr, deviceInd));

  44.    printf("Loading model from '%s'...\n", MODEL_PATH);
  45.    MIRA_CHECK(mira_LoadModel(pmr, MODEL_PATH));

  47.    const int decisionCount = mira_GetDecCount(pmr);
  48.    printf("\n%d model decisions:\n", decisionCount);
  49.    for (int i = 0; i < decisionCount; i++) {
  50.        printf(" %d : %s\n", i, mira_GetDecName(pmr, i));
  51.    }

  53.    printf("Loading correction...\n");
  54.    MIRA_CHECK(mira_LoadCorrection(pmr, CUBE_DIR, CUBE_NAME));

  56.    printf("Expected geometry: width=%d bands=%d lines=%d dataType=%d dataLayout=%d\n",
  57.           mira_GetInputWidth(pmr),
  58.           mira_GetInputBands(pmr),
  59.           mira_GetInputHeight(pmr),
  60.           mira_GetInputDataType(pmr),
  61.           mira_GetInputDataLayout(pmr));

  63.    if (mira_GetInputBands(pmr) != CUBE_BANDS) {
  64.        printf("The number of cube bands (%d) does not match the number of expected bands\n", CUBE_BANDS);
  65.        goto Error;
  66.    }
  67.    if (mira_GetInputHeight(pmr) != 1) {
  68.        printf("Line-scan model expected by this example\n", CUBE_BANDS);
  69.        goto Error;
  70.    }

  72.    printf("Loading data from '%s/%s'...\n", CUBE_DIR, CUBE_NAME CUBE_EXT);
  73.    FILE* pFile = fopen(CUBE_DIR "/" CUBE_NAME CUBE_EXT, "rb");
  74.    if (pFile == NULL) {
  75.        printf("Cannot open file for reading\n");
  76.        goto Error;
  77.    }

  79.    size_t frameSize = CUBE_WIDTH * CUBE_BANDS * sizeof(uint16_t);
  80.    pBuf = malloc(frameSize * CUBE_HEIGHT);

  82.    size_t frames = fread(pBuf, frameSize, CUBE_HEIGHT, pFile);
  83.    printf("Read %zu frames\n", frames);

  85.    fclose(pFile);

  87.    printf("Setting segmentation...\n");
  88.    MIRA_CHECK(mira_SetMinObjSize(pmr, 500));
  89.    MIRA_CHECK(mira_SetSegmentation(pmr, 1));

  91.    printf("Starting acquisition...\n");
  92.    MIRA_CHECK(mira_StartAcquisition(pmr));

  94.    uint16_t* ptr = pBuf;

  96.    int objectCount = 0;
  97.    int* pObjectData = NULL;

  99.    printf("Processing %zu frames:\n", frames);
  100.    for (int i = 0; i < frames; i++) {
  101.        mira_ProcessFrame(pmr, ptr);

  103.        objectCount = mira_GetObjCount(pmr);
  104.        if (objectCount > 0) {
  105.            for (int j = 0; j < objectCount; j++) {
  106.                mira_GetObjDataInt(pmr, j, &pObjectData);

  108.                printf("Object %d : %d,%d size=%d decision=%d box=[%d:%d,%d:%d]\n",
  109.                       pObjectData[MIRA_OBJECT_ID],
  110.                       pObjectData[MIRA_OBJECT_POS], pObjectData[MIRA_OBJECT_FRAME],
  111.                       pObjectData[MIRA_OBJECT_SIZE],
  112.                       pObjectData[MIRA_OBJECT_CLASS],
  113.                       pObjectData[MIRA_OBJECT_MINCOL], pObjectData[MIRA_OBJECT_MAXCOL],
  114.                       pObjectData[MIRA_OBJECT_MINFRAME], pObjectData[MIRA_OBJECT_MAXFRAME]);
  115.            }
  116.        }

  118.        ptr += frameSize;
  119.    }

  121.    printf("Stopping acquisition\n");
  122.    MIRA_CHECK(mira_StopAcquisition(pmr));

  124. Error:
  125.    if (res != MIRA_OK) {
  126.        printf("Error %d: %s", mira_GetErrorCode(pmr), mira_GetErrorMsg(pmr));
  127.    }

  129.    mira_Release(pmr);

  131.    if (pBuf != NULL) free(pBuf);

  133.    return 0;
  134. }