{"id":358961,"date":"2024-05-21T00:22:38","date_gmt":"2024-05-21T00:22:38","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=358961"},"modified":"-0001-11-30T00:00:00","modified_gmt":"2024-05-21T00:22:38","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=358961","title":{"rendered":"\u0421\u0447\u0438\u0442\u0430\u0435\u043c \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u0431\u044b\u0441\u0442\u0440\u0435\u0435 numpy<\/span>"},"content":{"rendered":"
<\/div>\n
\n
\n
\n
\n

\u0423\u0432\u0430\u0436\u0430\u0435\u043c\u044b\u0435 \u043a\u043e\u043b\u043b\u0435\u0433\u0438! \u0412\u0430\u0448\u0435\u043c\u0443 \u0432\u043d\u0438\u043c\u0430\u043d\u0438\u044e \u043f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u043d\u0435\u0431\u043e\u043b\u044c\u0448\u043e\u0439 «\u044d\u0442\u044e\u0434 \u0432\u044b\u0445\u043e\u0434\u043d\u043e\u0433\u043e \u0434\u043d\u044f», \u043f\u043e\u0441\u0432\u044f\u0449\u0451\u043d\u043d\u044b\u0439 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e, \u0441\u043a\u0430\u0436\u0435\u043c \u0442\u0430\u043a, \u043d\u0435\u0442\u0440\u0430\u0434\u0438\u0446\u0438\u043e\u043d\u043d\u043e\u043c\u0443 \u0441\u043f\u043e\u0441\u043e\u0431\u0443 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u043c\u0430\u0441\u0441\u0438\u0432\u0430 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441 \u043f\u043b\u0430\u0432\u0430\u044e\u0449\u0435\u0439 \u0442\u043e\u0447\u043a\u043e\u0439. \u0412\u043a\u0440\u0430\u0442\u0446\u0435 \u2014 \u043c\u044b \u0441\u0434\u0435\u043b\u0430\u0435\u043c \u044d\u0442\u043e \u0432 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043f\u0440\u043e\u0445\u043e\u0434\u043e\u0432 \u043f\u043e \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c\u0443 \u043c\u0430\u0441\u0441\u0438\u0432\u0443 (\u0434\u0432\u0430 \u0434\u043b\u044f \u043e\u0434\u0438\u043d\u0430\u0440\u043d\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438 \u0438\u043b\u0438 \u0447\u0435\u0442\u044b\u0440\u0435 \u0434\u043b\u044f \u0434\u0432\u043e\u0439\u043d\u043e\u0439), \u0432\u044b\u0447\u0438\u0441\u043b\u044f\u044f \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u043f\u043e \u0441\u043b\u043e\u0432\u0430\u043c, \u043d\u0430\u0447\u0438\u043d\u0430\u044f \u0441 \u0431\u043e\u043b\u0435\u0435 \u0437\u043d\u0430\u0447\u0430\u0449\u0438\u0445, \u043f\u043e\u043b\u044c\u0437\u0443\u044f\u0441\u044c \u043f\u0440\u0438 \u044d\u0442\u043e\u043c \u0442\u043e\u043b\u044c\u043a\u043e \u0446\u0435\u043b\u043e\u0447\u0438\u0441\u043b\u0435\u043d\u043d\u043e\u0439 \u0430\u0440\u0438\u0444\u043c\u0435\u0442\u0438\u043a\u043e\u0439, \u0447\u0442\u043e \u0434\u0430\u0441\u0442 \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u043e\u0441\u0442\u044c \u0432 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u043b\u0443\u0447\u0430\u044f\u0445 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043e\u0431\u043e\u0433\u043d\u0430\u0442\u044c \u043f\u043e \u0441\u043a\u043e\u0440\u043e\u0441\u0442\u0438 «\u0442\u0440\u0430\u0434\u0438\u0446\u0438\u043e\u043d\u043d\u044b\u0435» \u043a\u043b\u0430\u0441\u0441\u0438\u0447\u0435\u0441\u043a\u0438\u0435 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u044b. \u0412\u043e\u0437\u043c\u043e\u0436\u043d\u043e \u0434\u0430\u043d\u043d\u0430\u044f «\u0437\u0430\u0440\u0438\u0441\u043e\u0432\u043a\u0430» \u043a\u0430\u043a \u0438\u0434\u0435\u044f \u043e\u043a\u0430\u0436\u0435\u0442\u0441\u044f \u043a\u043e\u043c\u0443-\u043d\u0438\u0431\u0443\u0434\u044c \u043f\u043e\u043b\u0435\u0437\u043d\u0430.<\/p>\n

\u041f\u0440\u0435\u0436\u0434\u0435 \u0447\u0435\u043c \u043f\u0440\u043e\u0434\u043e\u043b\u0436\u0438\u0442\u044c \u2014 \u043c\u0430\u043b\u0435\u043d\u044c\u043a\u0438\u0439 \u0434\u0438\u0441\u043a\u043b\u0435\u0439\u043c\u0435\u0440: \u044f \u043d\u0438\u043a\u0430\u043a \u043d\u0435 \u0441\u0432\u044f\u0437\u0430\u043d \u043d\u0438 \u0441 \u043a\u0430\u043a\u043e\u0439 \u0438\u0437 \u0443\u043f\u043e\u043c\u0438\u043d\u0430\u044e\u0449\u0438\u0445\u0441\u044f \u043d\u0438\u0436\u0435 \u043a\u043e\u043c\u043f\u0430\u043d\u0438\u0439, \u043f\u043e\u0441\u0442 \u043d\u0435 \u043d\u043e\u0441\u0438\u0442 \u0440\u0435\u043a\u043b\u0430\u043c\u043d\u043e\u0433\u043e \u0445\u0430\u0440\u0430\u043a\u0442\u0435\u0440\u0430, \u043a\u043e\u0434 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 \u043d\u0430 \u0441\u0432\u043e\u0439 \u0441\u0442\u0440\u0430\u0445 \u0438 \u0440\u0438\u0441\u043a.<\/p>\n

\u0410 \u043d\u0430\u0447\u0430\u043b\u0430\u0441\u044c \u044d\u0442\u0430 \u0438\u0441\u0442\u043e\u0440\u0438\u044f \u0441 \u0440\u044f\u0434\u043e\u0432\u043e\u0439 \u0432 \u043e\u0431\u0449\u0435\u043c-\u0442\u043e \u0437\u0430\u0434\u0430\u0447\u043a\u0438 \u2014 \u043c\u043d\u0435 \u043d\u0443\u0436\u043d\u043e \u0431\u044b\u043b\u043e \u043d\u0430\u0439\u0442\u0438 \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u0440\u0435\u043d\u0442\u0433\u0435\u043d\u043e\u0432\u0441\u043a\u043e\u0439 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438 (\u044d\u0442\u043e \u043e\u0431\u044b\u0447\u043d\u043e 12-14 \u0431\u0438\u0442 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435) \u0432 \u0440\u0430\u043c\u043a\u0430\u0445 \u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 Flat Field Correction<\/a>. \u041c\u043e\u044f \u043e\u0441\u043d\u043e\u0432\u043d\u0430\u044f \u0440\u0430\u0431\u043e\u0447\u0430\u044f \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 — NI Vision Development Module, \u0438 \u0442\u0430\u043c \u0435\u0441\u0442\u044c \u0438 \u0441\u0440\u0435\u0434\u043d\u0435\u0435 \u0438 \u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u043e\u0435 \u043e\u0442\u043a\u043b\u043e\u043d\u0435\u043d\u0438\u0435, \u0438 \u043c\u043d\u043e\u0433\u043e \u0447\u0435\u0433\u043e \u0434\u0440\u0443\u0433\u043e\u0433\u043e, \u0430 \u0432\u043e\u0442 \u043f\u0440\u0438\u043c\u0438\u0442\u0438\u0432\u043d\u043e\u0439 \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u043d\u0435 \u043e\u043a\u0430\u0437\u0430\u043b\u043e\u0441\u044c. \u041d\u0430 «\u043f\u043e\u0434\u0445\u0432\u0430\u0442\u0435» \u044f \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e \u0442\u0430\u043a\u0436\u0435 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0438 OpenCV \u043b\u0438\u0431\u043e Intel IPP, \u043d\u043e \u0438 \u0442\u0430\u043c \u043d\u0430\u0432\u0441\u043a\u0438\u0434\u043a\u0443 \u0442\u0430\u043a\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0432 \u0433\u043e\u0442\u043e\u0432\u043e\u043c \u0432\u0438\u0434\u0435 \u044f \u043d\u0435 \u043e\u0431\u043d\u0430\u0440\u0443\u0436\u0438\u043b (\u0430 \u043c\u043e\u0436\u0435\u0442 \u043f\u043b\u043e\u0445\u043e \u0438\u0441\u043a\u0430\u043b). \u041d\u0443 \u043d\u0435 \u0431\u0435\u0434\u0430 \u2014 \u0440\u0430\u0441\u0447\u0435\u0445\u043b\u0438\u043c \u043a\u043e\u043c\u043f\u0438\u043b\u044f\u0442\u043e\u0440 \u0434\u0430 \u043d\u0430\u043f\u0438\u0448\u0435\u043c \u0441\u0430\u043c\u0438. \u0414\u0430, \u0435\u0441\u043b\u0438 \u0432\u044b \u0437\u043d\u0430\u0435\u0442\u0435, \u043a\u0430\u043a \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u0439 \u043c\u0435\u0442\u043e\u0434, \u0442\u043e \u043c\u043e\u0436\u0435\u0442\u0435 \u0441\u0440\u0430\u0437\u0443 \u043f\u0440\u043e\u043c\u0430\u0442\u044b\u0432\u0430\u0442\u044c \u043f\u0440\u0438\u043c\u0435\u0440\u043d\u043e \u0434\u043e \u0441\u0435\u0440\u0435\u0434\u0438\u043d\u044b, \u043b\u0438\u0431\u043e \u043f\u0440\u044f\u043c\u043e \u0434\u043e \u043a\u043e\u043d\u0446\u0430, \u0442\u0430\u043c \u0433\u0434\u0435 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b \u0437\u0430\u043c\u0435\u0440\u0430 \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438.<\/p>\n

\u041f\u0440\u0438 \u0440\u0430\u0431\u043e\u0442\u0435 \u043f\u043e \u0448\u0435\u0441\u0442\u043d\u0430\u0434\u0446\u0430\u0442\u0438\u0431\u0438\u0442\u043d\u043e\u043c\u0443 \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u043e\u0431\u044b\u0447\u043d\u043e \u043f\u0440\u0438\u043d\u044f\u0442\u043e \u0432\u044b\u0447\u0438\u0441\u043b\u044f\u0442\u044c, \u043f\u043e\u043b\u044c\u0437\u0443\u044f\u0441\u044c \u043d\u0435\u0441\u043b\u043e\u0436\u043d\u044b\u043c \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u043c \u043c\u0435\u0442\u043e\u0434\u043e\u043c. \u041c\u0435\u0434\u0438\u0430\u043d\u0430 \u2014 \u044d\u0442\u043e \u0432\u0435\u0434\u044c \u0442\u0430\u043a\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435, \u043a\u043e\u0442\u043e\u0440\u0430\u044f \u0434\u0435\u043b\u0438\u0442 \u043c\u0430\u0441\u0441\u0438\u0432 (\u0435\u0441\u043b\u0438 \u0435\u0433\u043e \u043e\u0442\u0441\u043e\u0440\u0442\u0438\u0440\u043e\u0432\u0430\u0442\u044c, \u0440\u0430\u0437\u0443\u043c\u0435\u0435\u0442\u0441\u044f) \u043f\u043e\u043f\u043e\u043b\u0430\u043c, \u0442\u0430\u043a \u0447\u0442\u043e \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u0430 \u0432\u0441\u0435\u0445 \u043e\u0441\u0442\u0430\u043b\u044c\u043d\u044b\u0445 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043f\u043e \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u0435 \u0431\u0443\u0434\u0435\u0442 \u0441\u043b\u0435\u0432\u0430 \u043e\u0442 \u043c\u0435\u0434\u0438\u0430\u043d\u044b, \u0430 \u0434\u0440\u0443\u0433\u0430\u044f \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u0430 — \u0441\u043f\u0440\u0430\u0432\u0430. \u041f\u043e\u044d\u0442\u043e\u043c\u0443 \u0437\u0430 \u043e\u0434\u0438\u043d \u043f\u0440\u043e\u0445\u043e\u0434 \u043f\u043e \u0432\u0441\u0435\u043c\u0443 \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u0441\u0442\u0440\u043e\u044f\u0442 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443, \u0430 \u0434\u0430\u043b\u044c\u0448\u0435 \u043f\u0435\u0440\u0435\u043c\u0435\u0449\u0430\u044e\u0442\u0441\u044f \u043f\u043e \u043d\u0435\u0439, \u043e\u0441\u0442\u0430\u043d\u0430\u0432\u043b\u0438\u0432\u0430\u044f\u0441\u044c \u0432 \u0442\u043e\u0442 \u043c\u043e\u043c\u0435\u043d\u0442, \u043a\u043e\u0433\u0434\u0430 \u0447\u0438\u0441\u043b\u043e \u043e\u0442\u0441\u0447\u0451\u0442\u043e\u0432 \u0434\u043e\u0441\u0442\u0438\u0433\u043d\u0435\u0442 \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u044b \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0430 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432: <\/p>\n

<\/figure>\n

\u041f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u044f \u0443\u0436\u0435 \u043e\u0442\u043d\u043e\u0441\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043d\u0435\u043c\u043e\u043b\u043e\u0434 \u0438 \u0432 \u043c\u0435\u0440\u0443 \u043b\u0435\u043d\u0438\u0432, \u0442\u043e \u0434\u043b\u044f \u0437\u0430\u0442\u0440\u0430\u0432\u043a\u0438 \u044f \u0432\u043e\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e\u0441\u044c ChatGPT. \u041d\u0443 \u0434\u0430\u0432\u0430\u0439, \u0433\u043e\u0440\u0448\u043e\u0447\u0435\u043a, \u0432\u0430\u0440\u0438! \u041f\u043e\u0434\u043d\u0430\u0442\u0443\u0436\u0438\u0432\u0448\u0438\u0441\u044c, \u043e\u043d\u043e \u043c\u043d\u0435 \u0432\u044b\u0434\u0430\u043b\u043e \u0432\u043e\u0442 \u0442\u0430\u043a\u043e\u0439 \u0448\u0435\u0434\u0435\u0432\u0440, \u043f\u0440\u0438\u0432\u0435\u0434\u0443 \u0435\u0433\u043e \u043a\u0430\u043a \u0435\u0441\u0442\u044c:<\/p>\n

unsigned short computeMedian(unsigned short* arr, int size) {     \/\/ Create histogram array     int maxVal = 65536; \/\/ Maximum value for unsigned short     int* histogram = (int*)calloc(maxVal, sizeof(int));      \/\/ Populate histogram     for (int i = 0; i < size; i++) {         histogram[arr[i]]++;     }      \/\/ Find the median value     int medianPos = (size + 1) \/ 2;     unsigned short median = 0;     int count = 0;      for (int i = 0; i < maxVal; i++) {         count += histogram[i];         if (count >= medianPos) {             median = i;             break;         }     }      \/\/ Free memory     free(histogram);      return median; } <\/code><\/pre>\n
\u041d\u0430 \u0441\u0430\u043c\u043e\u043c \u0434\u0435\u043b\u0435 \u043d\u0435\u043f\u043b\u043e\u0445\u043e, \u043e\u043d\u043e \u0434\u0430\u0436\u0435 \u0437\u043d\u0430\u0435\u0442, \u0447\u0442\u043e calloc \u0438\u043d\u0438\u0446\u0438\u0430\u043b\u0438\u0437\u0438\u0440\u0443\u0435\u0442 \u043f\u0430\u043c\u044f\u0442\u044c, \u043d\u043e \u0435\u0441\u0442\u044c \u043d\u044e\u0430\u043d\u0441 \u2014 \u043c\u043d\u0435 \u0434\u0430\u0436\u0435 \u0437\u0430\u043f\u0443\u0441\u043a\u0430\u0442\u044c \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u043d\u0435 \u043d\u0443\u0436\u043d\u043e, \u0447\u0442\u043e\u0431\u044b \u043f\u043e\u043d\u044f\u0442\u044c, \u0447\u0442\u043e \u044d\u0442\u043e \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c, \u0442\u043e\u043b\u044c\u043a\u043e \u0435\u0441\u043b\u0438 \u0432 \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u043d\u0435\u0447\u0451\u0442\u043d\u043e\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432. \u0410 \u0432 \u0447\u0451\u0442\u043d\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435 \u043c\u044b \u0431\u0435\u0440\u0451\u043c \u0441\u0440\u0435\u0434\u043d\u0435\u0435 \u0430\u0440\u0438\u0444\u043c\u0435\u0442\u0438\u0447\u0435\u0441\u043a\u043e\u0435 \u043e\u0442 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u043b\u0435\u0432\u0430 \u0438 \u0441\u043f\u0440\u0430\u0432\u0430 \u043e\u0442 \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u0438 \u043c\u043e\u0439 \u0432\u043d\u0443\u0442\u0440\u0435\u043d\u043d\u0438\u0439 \u043f\u0435\u0440\u0444\u0435\u043a\u0446\u0438\u043e\u043d\u0438\u0441\u0442 \u043f\u0440\u043e\u0441\u0438\u0442, \u0447\u0442\u043e\u0431\u044b \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0442\u043e\u0447\u043d\u043e \u0441\u043e\u0432\u043f\u0430\u0434\u0430\u043b\u043e \u0441 \u0440\u0435\u0444\u0435\u0440\u0435\u043d\u0441\u043d\u044b\u043c.  \u041f\u043e\u043f\u0440\u043e\u0441\u0438\u043c \u043f\u043e\u043f\u0440\u0430\u0432\u0438\u0442\u044c:<\/p>\n
double computeMedian(unsigned short* arr, int size) {     \/\/ Create histogram array     int maxVal = 65536; \/\/ Maximum value for unsigned short     int* histogram = (int*)calloc(maxVal, sizeof(int));      \/\/ Populate histogram     for (int i = 0; i < size; i++) {         histogram[arr[i]]++;     }      \/\/ Find the median position     int medianPos = (size + 1) \/ 2;     int count = 0;     int medianVal = 0;      \/\/ Iterate over histogram to find the median value     for (int i = 0; i < maxVal; i++) {         count += histogram[i];         if (count >= medianPos) {             medianVal = i;             break;         }     }      \/\/ If the size is odd, return the median value directly     if (size % 2 != 0) {         return (double)medianVal;     }      \/\/ If the size is even, find the next value to compute the median     int nextVal = medianVal + 1;     while (histogram[nextVal] == 0) {         nextVal++;     }      \/\/ Compute the median value for even-sized array     return (medianVal + nextVal) \/ 2.0; } <\/code><\/pre>\n
\u0423\u0436\u0435 \u043b\u0443\u0447\u0448\u0435, \u043d\u043e \u0432\u043e-\u043f\u0435\u0440\u0432\u044b\u0445, \u043e\u043d\u043e \u043f\u043e\u0447\u0435\u043c\u0443-\u0442\u043e \u043f\u043e\u0442\u0435\u0440\u044f\u043b\u043e \u043e\u0441\u0432\u043e\u0431\u043e\u0436\u0434\u0435\u043d\u0438\u0435 \u0432\u044b\u0434\u0435\u043b\u0435\u043d\u043d\u043e\u0439 \u043f\u043e\u0434 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u043f\u0430\u043c\u044f\u0442\u0438, \u0430 \u0432\u043e-\u0432\u0442\u043e\u0440\u044b\u0445, \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u0431\u0443\u0434\u0435\u0442 \u043f\u0440\u043e\u0444\u0443\u043a\u0438\u0432\u0430\u0442\u044c \u0438\u0441\u0442\u0438\u043d\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u0435\u0434\u0438\u0430\u043d\u044b, \u0435\u0441\u043b\u0438 \u0432 \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u0431\u0443\u0434\u0443\u0442 \u043f\u043e\u0432\u0442\u043e\u0440\u044f\u044e\u0449\u0438\u0435\u0441\u044f \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u044b. \u0414\u0430\u043b\u044c\u043d\u0435\u0439\u0448\u0438\u0435 \u043f\u0440\u043e\u0441\u044c\u0431\u044b \u0438\u0441\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u0438\u0439 \u0431\u0443\u0434\u0443\u0442 \u0440\u0430\u0437\u0434\u0443\u0432\u0430\u0442\u044c \u043a\u043e\u0434 \u0435\u0449\u0451 \u0441\u0438\u043b\u044c\u043d\u0435\u0435 (\u0430 \u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u0438\u0441\u0442\u043e\u0432 \u0441\u043a\u043e\u043f\u0438\u043f\u0430\u0441\u0442\u044f\u0442 \u0435\u0433\u043e \u043f\u0440\u044f\u043c\u043e \u0432 \u0442\u0430\u043a\u043e\u043c \u0432\u0438\u0434\u0435?). \u041f\u043e\u0436\u0430\u043b\u0443\u0439 \u0431\u0435\u0437 \u0440\u0430\u0431\u043e\u0442\u044b ChatGPT \u043d\u0430\u0441 \u043f\u043e\u043a\u0430 \u0447\u0442\u043e \u043d\u0435 \u043e\u0441\u0442\u0430\u0432\u0438\u0442. \u0417\u0430\u0441\u0443\u0447\u0438\u0432\u0430\u0435\u043c \u0440\u0443\u043a\u0430\u0432\u0430 \u0438 \u043f\u0440\u0430\u0432\u0438\u043c \u043a\u043e\u0434 \u0440\u0443\u0447\u043a\u0430\u043c\u0438, \u043d\u0443 \u043a\u0430\u043a-\u0442\u043e \u0432\u043e\u0442 \u0442\u0430\u043a \u0447\u0442\u043e \u043b\u0438:<\/p>\n
double computeMedian(unsigned short* arr, size_t size) { size_t* histogram; double Median; if (!arr || !size) return -1; if (!( histogram = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t)))) return -2;  for (int i = 0; i < size; i++) histogram[arr[i]]++;    \/\/ Populate histogram size_t medianPos = (size + 1) \/ 2, count = 0, medianVal; \/\/ Find the median pos for (medianVal = 0; medianVal <= USHRT_MAX; medianVal++) { \/\/ Iterate over histo count += histogram[medianVal]; if (count >= medianPos) break; }  if (!(size % 2 || count > medianPos)) { \/\/ If the size is even, find the next size_t nextVal = medianVal; while (!histogram[++nextVal]); Median = (medianVal + nextVal) \/ 2.0; \/\/ Middle } else Median = (double)medianVal;  free(histogram); return Median; } <\/code><\/pre>\n
\u041d\u0430 \u043a\u043e\u0434 \u0440\u0435\u0432\u044c\u044e, \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e, \u043c\u043e\u0433\u0443\u0442 \u0432\u043e\u0437\u043d\u0438\u043a\u043d\u0443\u0442\u044c \u043a\u043e\u0435-\u043a\u0430\u043a\u0438\u0435 \u0432\u043e\u043f\u0440\u043e\u0441\u044b, \u043d\u043e \u0432\u0441\u0435\u0433\u0434\u0430 \u043c\u043e\u0436\u043d\u043e \u0441\u0432\u0430\u043b\u0438\u0442\u044c \u043d\u0430 \u0438\u0441\u043a\u0443\u0441\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0439 \u0438\u043d\u0442\u0435\u043b\u043b\u0435\u043a\u0442, \u043e\u043f\u044f\u0442\u044c \u0436\u0435 \u0441\u0442\u0430\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0439 \u0430\u043d\u0430\u043b\u0438\u0437\u0430\u0442\u043e\u0440 \u043c\u043e\u043b\u0447\u0438\u0442 \u043a\u0430\u043a \u043f\u0430\u0440\u0442\u0438\u0437\u0430\u043d.<\/p>\n
32 \u0431\u0438\u0442\u0430<\/h3>\n
\u0420\u044f\u0434\u043e\u0432\u0430\u044f \u0432 \u043e\u0431\u0449\u0435\u043c-\u0442\u043e \u0437\u0430\u0434\u0430\u0447\u043a\u0430 \u043d\u0435 \u0441\u0442\u043e\u0438\u043b\u0430 \u0431\u044b \u043f\u043e\u0441\u0442\u0430 \u043d\u0430 \u0425\u0430\u0431\u0440\u0435,  \u0434\u0430 \u0438 \u0435\u0433\u043e \u0431\u044b \u0438 \u043d\u0435 \u043f\u043e\u044f\u0432\u0438\u043b\u043e\u0441\u044c, \u0435\u0441\u043b\u0438 \u0431\u044b \u0432\u043e \u0432\u0440\u0435\u043c\u044f \u0435\u0436\u0435\u0432\u0435\u0447\u0435\u0440\u043d\u0435\u0439 \u043f\u0440\u043e\u0433\u0443\u043b\u043a\u0438 \u0441 \u0441\u043e\u0431\u0430\u043a\u0435\u043d\u044b\u043c \u043c\u043e\u0439 \u0432\u043d\u0443\u0442\u0440\u0435\u043d\u043d\u0438\u0439 \u0433\u043e\u043b\u043e\u0441 \u0435\u0445\u0438\u0434\u043d\u043e \u0438 \u0432\u043a\u0440\u0430\u0434\u0447\u0438\u0432\u043e \u043d\u0435 \u043f\u043e\u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043e\u0432\u0430\u043b\u0441\u044f: «\u2014 \u0430 \u0432\u043e\u0442 \u0435\u0441\u043b\u0438, \u043f\u0440\u0435\u0434\u043f\u043e\u043b\u043e\u0436\u0438\u043c, \u0443 \u0442\u0435\u0431\u044f \u043c\u0430\u0441\u0441\u0438\u0432 \u043d\u0435 16-\u0438, \u0430 32-\u0445 \u0431\u0438\u0442\u043d\u044b\u0439, \u0442\u043e \u043a\u0430\u043a \u0442\u043e\u0433\u0434\u0430?»<\/em>. \u0412\u0435\u0434\u044c \u043d\u0430 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0432 \u0447\u0435\u0442\u044b\u0440\u0435 \u0441 \u043b\u0438\u0448\u043d\u0438\u043c \u043c\u0438\u043b\u043b\u0438\u0430\u0440\u0434\u043e\u0432 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043d\u0438\u043a\u0430\u043a\u043e\u0439 \u043f\u0430\u043c\u044f\u0442\u0438 \u043d\u0435 \u0445\u0432\u0430\u0442\u0438\u0442! \u041f\u043e\u0440\u0430\u0437\u043c\u044b\u0448\u043b\u044f\u0432 \u043d\u0435\u043c\u043d\u043e\u0433\u043e, \u044f \u0435\u0449\u0451 \u0434\u043e \u043a\u043e\u043d\u0446\u0430 \u043f\u0440\u043e\u0433\u0443\u043b\u043a\u0438 \u0438\u0437\u043c\u044b\u0441\u043b\u0438\u043b \u043f\u0440\u043e\u0441\u0442\u043e\u0439 \u043d\u0435\u0437\u0430\u0442\u0435\u0439\u043b\u0438\u0432\u044b\u0439 \u043c\u0435\u0442\u043e\u0434 \u2014 \u0432\u043d\u0430\u0447\u0430\u043b\u0435 \u043d\u0430\u0434\u043e \u043f\u0440\u043e\u0439\u0442\u0438 \u043f\u043e \u0441\u0442\u0430\u0440\u0448\u0438\u043c \u0441\u043b\u043e\u0432\u0430\u043c, \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u0432 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0432 65 \u043a\u0438\u043b\u043e\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043a\u0430\u043a \u0438 \u0432\u044b\u0448\u0435, \u0430 \u0437\u0430\u0442\u0435\u043c, \u043e\u0442\u0442\u0430\u043b\u043a\u0438\u0432\u0430\u044f\u0441\u044c \u043e\u0442 \u043d\u0430\u0439\u0434\u0435\u043d\u043d\u043e\u0433\u043e \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f \u043f\u0440\u043e\u0439\u0442\u0438 \u043f\u043e \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u0435\u0449\u0451 \u0440\u0430\u0437, \u0441\u0442\u0440\u043e\u044f \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0443\u0436\u0435 \u043f\u043e \u043c\u043b\u0430\u0434\u0448\u0435\u043c\u0443 \u0441\u043b\u043e\u0432\u0443 \u0438 \u0443\u0442\u043e\u0447\u043d\u044f\u044f \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u0435\u0434\u0438\u0430\u043d\u044b. \u0412 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435 \u043c\u0435\u0442\u043e\u0434 \u043d\u0435 \u043d\u043e\u0432 \u2014 \u043c\u044b \u0444\u0430\u043a\u0442\u0438\u0447\u0435\u0441\u043a\u0438 \u0440\u0430\u0437\u0431\u0438\u0432\u0430\u0435\u043c \u0432\u0435\u0441\u044c \u0434\u0438\u0430\u043f\u0430\u0437\u043e\u043d \u043d\u0430 65536 \u043e\u0442\u0440\u0435\u0437\u043a\u043e\u0432, \u0438 \u043d\u0430\u0439\u0434\u044f \u043f\u043e\u0434\u0445\u043e\u0434\u044f\u0449\u0438\u0439 \u043e\u0442\u0440\u0435\u0437\u043e\u043a, \u0438\u0449\u0435\u043c \u0442\u0430\u043c \u043e\u043a\u043e\u043d\u0447\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435. \u041c\u0430\u043b\u0435\u043d\u044c\u043a\u0438\u0439 \u043d\u044e\u0430\u043d\u0441 \u0435\u0441\u0442\u044c \u0441 \u0447\u0430\u0441\u0442\u043d\u044b\u043c \u0441\u043b\u0443\u0447\u0430\u0435\u043c \u043a\u043e\u0433\u0434\u0430 \u0432 \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u0447\u0451\u0442\u043d\u043e\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0438 \u043d\u0430\u0448\u0430 \u043c\u0435\u0434\u0438\u0430\u043d\u0430 \u043b\u0435\u0436\u0438\u0442 \u043c\u0435\u0436\u0434\u0443 \u0440\u0430\u0437\u043d\u044b\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438 \u0441\u0442\u0430\u0440\u0448\u0435\u0433\u043e \u0441\u043b\u043e\u0432\u0430. \u041d\u043e \u0435\u0441\u043b\u0438 \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u0438\u0442\u044c \u0432 \u0443\u043c\u0435 \u043e\u0442\u0441\u043e\u0440\u0442\u0438\u0440\u043e\u0432\u0430\u043d\u043d\u044b\u0439 \u043c\u0430\u0441\u0441\u0438\u0432, \u0442\u043e \u0441\u0442\u0430\u043d\u0435\u0442 \u044f\u0441\u043d\u043e, \u0447\u0442\u043e \u043e\u0442\u0434\u0435\u043b\u044c\u043d\u044b\u0435 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u044b \u0442\u0443\u0442 \u0443\u0436\u0435 \u043d\u0435 \u043d\u0443\u0436\u043d\u044b, \u043f\u0440\u043e\u0441\u0442\u043e \u0441\u043b\u0435\u0432\u0430 (\u0434\u043b\u044f \u043c\u0435\u043d\u044c\u0448\u0435\u0433\u043e \u0441\u0442\u0430\u0440\u0448\u0435\u0433\u043e \u0431\u0430\u0439\u0442\u0430) \u0441\u043b\u0435\u0434\u0443\u0435\u0442 \u0432\u0437\u044f\u0442\u044c \u043c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u043b\u0430\u0434\u0448\u0435\u0433\u043e \u0431\u0430\u0439\u0442\u0430, \u0430 \u0441\u043f\u0440\u0430\u0432\u0430 (\u0434\u043b\u044f \u0431\u043e\u043b\u044c\u0448\u0435\u0433\u043e) \u2014 \u043c\u0438\u043d\u0438\u043c\u0430\u043b\u044c\u043d\u043e\u0435, \u0438 \u043f\u043e\u0442\u043e\u043c \u0443\u0436\u0435 \u043d\u0430\u0439\u0442\u0438 \u0441\u0440\u0435\u0434\u043d\u0435\u0435, \u0432\u043e\u0442 \u0438 \u0432\u0441\u0451: <\/p>\n
<\/figure>\n
 \u042f \u043d\u0430\u0431\u0440\u043e\u0441\u0430\u043b \u0431\u044b\u043b\u043e \u043a\u043e\u0434 \u043d\u0430 LabVIEW, \u043d\u043e \u043f\u043e\u0436\u0430\u043b\u0443\u0439 \u043d\u0435 \u0441\u0442\u0430\u043d\u0443 \u0432\u0430\u0441 \u043c\u0443\u0447\u0438\u0442\u044c, \u0445\u043e\u0442\u044c \u043c\u044b \u0438 \u0432 \u0445\u0430\u0431\u0435 «\u041d\u0435\u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0435». <\/p>\n
\u041d\u0430 \u0421\u0438 \u043e\u043d\u043e \u0431\u0443\u0434\u0435\u0442 \u0442\u0430\u043a \u0436\u0435, \u043a\u0430\u043a \u0438 \u0432 \u0441\u043b\u0443\u0447\u0430\u0435 16 \u0431\u0438\u0442, \u0432\u043e\u0442 \u0442\u043e\u043b\u044c\u043a\u043e \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u043c\u044b \u0441\u0442\u0440\u043e\u0438\u043c \u043f\u043e \u0441\u0442\u0430\u0440\u0448\u0438\u043c \u0431\u0438\u0442\u0430\u043c, \u0437\u0430\u0434\u0432\u0438\u0433\u0430\u044f \u0438\u0445 \u043d\u0430\u043f\u0440\u0430\u0432\u043e:<\/p>\n
LIBMED_API int MedianU32(unsigned int* ptr, size_t Length, double* med) {     \/\/...     histo = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t));     for (int i = 0; i < Length; i++) histo[(ptr[i] & 0xFFFF0000) >> 16]++;     \/\/... <\/code><\/pre>\n
\u0414\u043e \u044d\u0442\u043e\u0433\u043e \u043c\u043e\u043c\u0435\u043d\u0442\u0430 \u0432\u0441\u0451 \u0440\u043e\u0432\u043d\u043e \u043a\u0430\u043a \u0438 \u0440\u0430\u043d\u044c\u0448\u0435:<\/p>\n
uint64_t medianPos = (Length + 1) \/ 2, count = 0;  for (int i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram to find the med count += histo[i]; if (count >= medianPos) { medianValMSW = i; break; } } if (!(Length % 2 || count > medianPos)) { \/\/ If the size is even, find the next <\/code><\/pre>\n
\u0415\u0441\u043b\u0438 \u0443\u0441\u043b\u043e\u0432\u0438\u0435 \u0441\u0440\u0430\u0431\u043e\u0442\u0430\u043b\u043e, \u0442\u043e \u043d\u0430 \u0432\u0442\u043e\u0440\u043e\u043c \u043f\u0440\u043e\u0445\u043e\u0434\u0435 \u043f\u0440\u043e\u0441\u0442\u043e \u0441\u043c\u043e\u0442\u0440\u0438\u043c \u043c\u0438\u043d\u0438\u043c\u0443\u043c \u0438 \u043c\u0430\u043a\u0441\u0438\u043c\u0443\u043c:<\/p>\n
if (!(Length % 2 || count > medianPos)) { \/\/ If the size is even, find the next nextVal = medianValMSW; while (!histo[++nextVal]);  for (max = 0, min = 65535, i = 0; i < Length; i++) { indexMSW = (ptr[i] & 0xFFFF0000) >> 16; if (medianValMSW == indexMSW && (ptr[i] & 0x0000FFFF) > max) max = (ptr[i] & 0x0000FFFF); if (indexMSW == nextVal && (ptr[i] & 0x0000FFFF) < min) min = (ptr[i] & 0x0000FFFF); }  medianLeft = (medianValMSW << 16) | max; medianRight = (nextVal << 16) | min; *med = ((double)medianLeft + (double)medianRight) \/ 2.0; <\/code><\/pre>\n
\u0410 \u0435\u0441\u043b\u0438 \u0443\u0441\u043b\u043e\u0432\u0438\u0435 \u043d\u0435 \u0441\u0440\u0430\u0431\u043e\u0442\u0430\u043b\u043e, \u0442\u043e \u0441\u0431\u0440\u0430\u0441\u044b\u0432\u0430\u0435\u043c \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0438 \u0434\u0435\u043b\u0430\u0435\u043c \u0432\u0442\u043e\u0440\u043e\u0439 \u043f\u0440\u043e\u0445\u043e\u0434 \u043f\u043e \u043c\u043b\u0430\u0434\u0448\u0435\u043c\u0443 \u0441\u043b\u043e\u0432\u0443 \u0432\u043e\u0442 \u0442\u0430\u043a, \u0437\u0430\u043d\u043e\u0441\u044f \u0432 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0442\u043e\u043b\u044c\u043a\u043e \u0442\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f \u0443 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u0442\u0430\u0440\u0448\u0435\u0435 \u0441\u043b\u043e\u0432\u043e \u0440\u0430\u0432\u043d\u043e \u043d\u0430\u0439\u0434\u0435\u043d\u043d\u043e\u043c\u0443 \u043d\u0430 \u043f\u0435\u0440\u0432\u043e\u043c \u043f\u0440\u043e\u0445\u043e\u0434\u0435 \u0438 \u043f\u043e\u043f\u0443\u0442\u043d\u043e \u0441\u0447\u0438\u0442\u0430\u044f \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0441\u043b\u0435\u0432\u0430:<\/p>\n
ZeroMemory(histo, (USHRT_MAX + 1) * sizeof(size_t)); med_ = medianValMSW << 16; for (count = 0, i = 0; i < Length; i++) { indexMSW = (ptr[i] & 0xFFFF0000) >> 16; if (medianValMSW == indexMSW) histo[(ptr[i] & 0x0000FFFF)]++; if (ptr[i] < med_) count++; }  \/\/ Find the median position; medianValLSW = 0; for (int i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram  count += histo[i]; if (count >= medianPos) { medianValLSW = i; break; } } <\/code><\/pre>\n
\u041d\u0443 \u0438 \u0437\u0434\u0435\u0441\u044c \u043d\u0430\u0441 \u043c\u043e\u0436\u0435\u0442 \u043d\u0430\u0441\u0442\u0438\u0433\u043d\u0443\u0442\u044c \u0442\u0430\u043a\u0430\u044f \u0436\u0435 \u0441\u0438\u0442\u0443\u0430\u0446\u0438\u044f, \u043a\u043e\u0433\u0434\u0430 \u043b\u0438\u0431\u043e \u043c\u0435\u0434\u0438\u0430\u043d\u0430 \u0441\u043e\u0431\u0438\u0440\u0430\u0435\u0442\u0441\u044f \u0438\u0437 \u0434\u0432\u0443\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439, \u043b\u0438\u0431\u043e \u0438\u0437 \u043e\u0434\u043d\u043e\u0433\u043e:<\/p>\n
if (!(Length % 2 || count > medianPos)) { nextVal = medianValLSW; while (!histo[++nextVal]);  medianLeft = (medianValMSW << 16) | medianValLSW; medianRight = (medianValMSW << 16) | nextVal; *med = ((double)medianLeft + (double)medianRight) \/ 2.0; } else { medianValI32 = (medianValMSW << 16) | medianValLSW; *med = (double)medianValI32; } <\/code><\/pre>\n
\u0412\u043e\u0442, \u0441\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e \u0438 \u0432\u0441\u0451. \u041f\u043e\u043d\u044f\u0442\u043d\u043e, \u0447\u0442\u043e \u043c\u0435\u0442\u043e\u0434 \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u0438 \u0434\u043b\u044f 64-\u0445 \u0431\u0438\u0442\u043d\u044b\u0445 \u043c\u0430\u0441\u0441\u0438\u0432\u043e\u0432, \u043f\u0440\u043e\u0441\u0442\u043e \u043f\u0440\u043e\u0439\u0442\u0438 \u043f\u043e \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u043f\u0440\u0438\u0434\u0451\u0442\u0441\u044f \u0447\u0435\u0442\u044b\u0440\u0435 \u0440\u0430\u0437\u0430, \u043f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u0443 \u043d\u0430\u0441 \u0447\u0435\u0442\u044b\u0440\u0435 \u0441\u043b\u043e\u0432\u0430 \u043f\u043e \u0434\u0432\u0430 \u0431\u0430\u0439\u0442\u0430, \u043f\u0440\u0438\u0447\u0451\u043c \u0447\u0430\u0441\u0442\u043d\u044b\u0439 \u0441\u043b\u0443\u0447\u0430\u0439 \u0441 \u0440\u0430\u0437\u043d\u044b\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438 \u0441\u043b\u0435\u0432\u0430 \u0438 \u0441\u043f\u0440\u0430\u0432\u0430 \u043f\u0440\u0438 \u0447\u0451\u0442\u043d\u043e\u043c \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0435 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043c\u043e\u0436\u0435\u0442 \u0432\u043e\u0437\u043d\u0438\u043a\u043d\u0443\u0442\u044c \u043d\u0430 \u043b\u044e\u0431\u043e\u043c \u043f\u0440\u043e\u0433\u043e\u043d\u0435.<\/p>\n
\n\u043f\u043e\u043b\u043d\u044b\u0439 \u043b\u0438\u0441\u0442\u0438\u043d\u0433 MedianU32<\/summary>\n
\n
 LIBMED_API int MedianU32(unsigned int* ptr, size_t Length, double* med) { size_t* histo; int i; uint64_t medianLeft, medianRight, indexMSW, nextVal, med_; uint64_t medianValI32, medianValLSW, medianValMSW = 0, min, max;  if (!ptr || !Length) return -1; histo = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t)); if (!histo) return -2;  \/\/MSW HISTO: for (int i = 0; i < Length; i++) histo[(ptr[i] & 0xFFFF0000) >> 16]++;  \/\/ Find the median position; uint64_t medianPos = (Length + 1) \/ 2, count = 0;  for (int i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram to find the med count += histo[i]; if (count >= medianPos) { medianValMSW = i; break; } } if (!(Length % 2 || count > medianPos)) { \/\/ If the size is even, find the next nextVal = medianValMSW; while (!histo[++nextVal]);  for (max = 0, min = 65535, i = 0; i < Length; i++) { indexMSW = (ptr[i] & 0xFFFF0000) >> 16; if (medianValMSW == indexMSW && (ptr[i] & 0x0000FFFF) > max)  max = (ptr[i] & 0x0000FFFF); if (indexMSW == nextVal && (ptr[i] & 0x0000FFFF) < min) min = (ptr[i] & 0x0000FFFF); }  medianLeft = (medianValMSW << 16) | max; medianRight = (nextVal << 16) | min; *med = ((double)medianLeft + (double)medianRight) \/ 2.0; } else { \/\/LSW HISTO ZeroMemory(histo, (USHRT_MAX + 1) * sizeof(size_t)); med_ = medianValMSW << 16; for (count = 0, i = 0; i < Length; i++) { indexMSW = (ptr[i] & 0xFFFF0000) >> 16; if (medianValMSW == indexMSW) histo[(ptr[i] & 0x0000FFFF)]++; if (ptr[i] < med_) count++; }  \/\/ Find the median position; medianValLSW = 0; for (int i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram  count += histo[i]; if (count >= medianPos) { medianValLSW = i; break; } }  if (!(Length % 2 || count > medianPos)) { nextVal = medianValLSW; while (!histo[++nextVal]);  medianLeft = (medianValMSW << 16) | medianValLSW; medianRight = (medianValMSW << 16) | nextVal; *med = ((double)medianLeft + (double)medianRight) \/ 2.0; } else { medianValI32 = (medianValMSW << 16) | medianValLSW; *med = (double)medianValI32; } } free(histo); return 0; } \/\/MedianU32 <\/code><\/pre>\n<\/div>\n<\/details>\n
Double<\/h3>\n
\u041d\u043e \u0432\u043d\u0443\u0442\u0440\u0435\u043d\u043d\u0438\u0439 \u0433\u043e\u043b\u043e\u0441 \u043d\u0435 \u043e\u0442\u0441\u0442\u0430\u0432\u0430\u043b \u2014 \u0432\u0435\u0434\u044c \u0435\u0441\u0442\u044c \u0436\u0435 \u0435\u0449\u0451 \u0447\u0438\u0441\u043b\u0430 \u0441 \u043f\u043b\u0430\u0432\u0430\u044e\u0449\u0435\u0439 \u0442\u043e\u0447\u043a\u043e\u0439, \u0442\u0430\u043a \u0447\u0442\u043e \u043a\u043e\u0434 \u044d\u0442\u043e\u0442 \u0434\u043b\u044f \u0431\u0435\u0437\u0437\u043d\u0430\u043a\u043e\u0432\u044b\u0445 32-\u0445 \u0431\u0438\u0442\u043d\u044b\u0445 \u0446\u0435\u043b\u044b\u0445  \u044f \u0434\u043e\u0442\u043e\u0448\u043d\u043e \u043d\u0435 \u043f\u0440\u043e\u0432\u0435\u0440\u044f\u043b, \u043f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u0443 \u043c\u0435\u043d\u044f \u0437\u0430\u0447\u0435\u0441\u0430\u043b\u0438\u0441\u044c \u0440\u0443\u043a\u0438 \u043f\u0440\u043e\u0434\u0435\u043b\u0430\u0442\u044c \u0442\u0430\u043a\u043e\u0439 \u0436\u0435 \u0444\u043e\u043a\u0443\u0441 \u0434\u043b\u044f \u0447\u0438\u0441\u0435\u043b \u0441 \u043f\u043b\u0430\u0432\u0430\u044e\u0449\u0435\u0439 \u0437\u0430\u043f\u044f\u0442\u043e\u0439, \u043f\u0440\u0438\u0447\u0451\u043c \u0434\u0432\u043e\u0439\u043d\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438 \u0438 \u0441\u043e \u0437\u043d\u0430\u043a\u043e\u043c, \u0432\u0435\u0434\u044c \u0435\u0441\u043b\u0438 \u0432\u0437\u0433\u043b\u044f\u043d\u0443\u0442\u044c \u043d\u0430 \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u0438\u0435 IEEE 754, \u043a\u043e\u0442\u043e\u0440\u044b\u043c \u043c\u044b \u043f\u043e\u043b\u044c\u0437\u0443\u0435\u043c\u0441\u044f, \u0442\u043e \u043b\u0435\u0433\u043a\u043e \u0437\u0430\u043c\u0435\u0442\u0438\u0442\u044c, \u0447\u0442\u043e \u0432 \u043f\u043e\u043b\u043e\u0436\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0439 \u043e\u0431\u043b\u0430\u0441\u0442\u0438 \u043c\u043e\u043d\u043e\u0442\u043e\u043d\u043d\u043e \u0432\u043e\u0437\u0440\u0430\u0441\u0442\u0430\u044e\u0449\u0438\u043c \u0447\u0438\u0441\u043b\u0430\u043c \u043d\u0430 \u0448\u043a\u0430\u043b\u0435 \u0441 \u043f\u043b\u0430\u0432\u0430\u044e\u0449\u0435\u0439 \u0442\u043e\u0447\u043a\u043e\u0439 \u0441\u043e\u043e\u0442\u0432\u0435\u0442\u0441\u0442\u0432\u0443\u044e\u0442 \u043c\u043e\u043d\u043e\u0442\u043e\u043d\u043d\u043e \u0432\u043e\u0437\u0440\u0430\u0441\u0442\u0430\u044e\u0449\u0438\u0435 \u0446\u0435\u043b\u044b\u0435 \u0447\u0438\u0441\u043b\u0430 \u0432 \u0442\u043e\u043c \u0436\u0435 \u0431\u0438\u0442\u043e\u0432\u043e\u043c \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u0438\u0438, \u0442\u043e \u0435\u0441\u0442\u044c \u0435\u0441\u043b\u0438 \u043f\u0440\u043e\u0441\u0442\u043e \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u0441 \u044d\u0442\u0438\u043c \u043c\u0430\u0441\u0441\u0438\u0432\u043e\u043c \u043a\u0430\u043a \u0441 \u043c\u0430\u0441\u0441\u0438\u0432\u043e\u043c \u0446\u0435\u043b\u044b\u0445 (\u0442\u043e\u0439 \u0436\u0435 \u0440\u0430\u0437\u0440\u044f\u0434\u043d\u043e\u0441\u0442\u0438), \u0430 \u043f\u043e\u0442\u043e\u043c \u043f\u0435\u0440\u0435\u0432\u0435\u0441\u0442\u0438 \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u0438\u0435 \u043e\u0431\u0440\u0430\u0442\u043d\u043e, \u0442\u043e \u043c\u044b \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u0438\u0441\u0442\u0438\u043d\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u0435\u0434\u0438\u0430\u043d\u044b.<\/p>\n
\u0412\u043e\u0442 \u0441\u043c\u043e\u0442\u0440\u0438\u0442\u0435 \u043d\u0430 \u043f\u0440\u0438\u043c\u0435\u0440\u0435 \u043e\u0434\u0438\u043d\u0430\u0440\u043d\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438:<\/p>\n
\u041d\u0443\u043b\u044c \u043e\u043d \u0442\u0430\u043a \u0438 \u0431\u0443\u0434\u0435\u0442 \u043d\u0443\u043b\u044c. \u0415\u0441\u043b\u0438 \u043c\u044b \u0432\u044b\u0441\u0442\u0430\u0432\u0438\u043c \u043f\u0435\u0440\u0432\u044b\u0439 \u0431\u0438\u0442, \u0442\u043e \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u043d\u0430\u0438\u043c\u0435\u043d\u044c\u0448\u0435\u0435 \u043f\u043e\u0437\u0438\u0442\u0438\u0432\u043d\u043e\u0435 \u0447\u0438\u0441\u043b\u043e: 1(0000 0001)->1,4012984643248171E-45 \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0435\u0435 \u0432\u043e\u0442 \u0442\u0430\u043a\u043e\u0435 2(0000 0002)->2,8025969286496341E-45 ... 65535(0000 FFFF)->9,1834094859526886E-41 65536(0001 0000)->9,1835496157991211E-41 ... \u0412\u043e\u0442 \u043d\u0430 \u044d\u0442\u043e\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0438 (\u0432\u044b\u0441\u0442\u0430\u0432\u0438\u0432 23-\u0439 \u0431\u0438\u0442) \u043c\u044b \u0434\u043e\u0441\u0442\u0438\u0433\u0430\u0435\u043c \u0442\u043e\u0447\u043d\u043e\u0439 \u0435\u0434\u0438\u043d\u0438\u0446\u044b: 1065353216(3F80 0000)->1,000000000000000000 ... \u0438 \u0434\u0430\u043b\u0435\u0435 \u0434\u043e \u043c\u0430\u043a\u0441\u0438\u043c\u0443\u043c\u0430: 2139095039(7F7F FFFF)->3,4028234663852886E+38 2139095040(7F80 0000)->Inf 2139095041(7F80 0001)->NaN ... \u0435\u0434\u0435\u043c \u0434\u0430\u043b\u044c\u0448\u0435 2147483647(7FFF FFFF)->NaN \u0410 \u0432\u043e\u0442 \u043f\u0440\u0438 \u0441\u0442\u0430\u0440\u0448\u0435\u043c \u0431\u0438\u0442\u0435 \u043c\u044b \u0443\u0445\u043e\u0434\u0438\u043c \u0432 \u043d\u0435\u0433\u0430\u0442\u0438\u0432\u043d\u044b\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f: 2147483648(8000 0000)->-0 2147483649(8000 0001)->-1,4012984643248171E-45 2147483650(8000 0002)->-2,8025969286496341E-45 3212836864(BF80 0000)->-1 4286578687(FF7F FFFF)->-3,4028234663852886E+38 4286578688(FF80 0000)->-Inf 4286578689(FF80 0001)->NaN ... <\/code><\/pre>\n
\u041e\u0441\u0442\u0430\u0451\u0442\u0441\u044f \u043b\u0438\u0448\u044c \u043f\u043e\u0441\u043b\u0435\u0434\u043d\u0438\u0439 \u0432\u043e\u043f\u0440\u043e\u0441 \u2014 \u043a\u0430\u043a \u043e\u0431\u043e\u0439\u0442\u0438\u0441\u044c \u0441 \u043d\u0435\u0433\u0430\u0442\u0438\u0432\u043d\u044b\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438? \u0415\u0441\u043b\u0438 \u043c\u044b \u0437\u0430\u0440\u0430\u043d\u0435\u0435 \u0437\u043d\u0430\u0435\u043c, \u0447\u0442\u043e \u043d\u0435\u0433\u0430\u0442\u0438\u0432\u043d\u044b\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0432 \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u043d\u0435 \u0431\u0443\u0434\u0435\u0442, \u0442\u043e \u043d\u0438\u043a\u0430\u043a\u0438\u0445 \u0434\u043e\u043f\u043e\u043b\u043d\u0438\u0442\u0435\u043b\u044c\u043d\u044b\u0445 \u0442\u0435\u043b\u043e\u0434\u0432\u0438\u0436\u0435\u043d\u0438\u0439 \u0434\u0435\u043b\u0430\u0442\u044c \u043d\u0435 \u043d\u0443\u0436\u043d\u043e \u2014 \u0432\u0441\u0451 \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u043a\u0430\u043a \u0441 \u0431\u0435\u0437\u0437\u043d\u0430\u043a\u043e\u0432\u044b\u043c\u0438 \u0446\u0435\u043b\u044b\u043c\u0438. \u0410 \u0432\u043e\u0442 \u0435\u0441\u043b\u0438 \u043d\u0435\u0442, \u0442\u043e \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u043e \u0441\u0431\u0440\u043e\u0441\u0438\u0442\u044c \u0441\u0442\u0430\u0440\u0448\u0438\u0439 \u0431\u0438\u0442 \u0443 \u043d\u0435\u0433\u0430\u0442\u0438\u0432\u043d\u044b\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439, \u0438 \u0443\u0441\u0442\u0430\u043d\u043e\u0432\u0438\u0442\u044c \u0435\u0433\u043e \u0443 \u043f\u043e\u0437\u0438\u0442\u0438\u0432\u043d\u044b\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 (\u0447\u0442\u043e \u0441\u0434\u0432\u0438\u043d\u0435\u0442 \u0438\u0445 \u0432\u043f\u0440\u0430\u0432\u043e \u043f\u043e \u0431\u0435\u0437\u0437\u043d\u0430\u043a\u043e\u0432\u043e\u0439 \u0448\u043a\u0430\u043b\u0435), \u043d\u043e \u043f\u0440\u0438 \u044d\u0442\u043e\u043c \u0442\u0430\u043a\u0436\u0435 «\u043f\u0435\u0440\u0435\u0432\u0435\u0440\u043d\u0443\u0442\u044c» \u0432\u0441\u044e \u043f\u043e\u0441\u043b\u0435\u0434\u043e\u0432\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u044c \u043d\u0435\u0433\u0430\u0442\u0438\u0432\u043d\u044b\u0445, \u0447\u0442\u043e\u0431\u044b \u043c\u0435\u043d\u044c\u0448\u0438\u0435 \u0447\u0438\u0441\u043b\u0430 \u0440\u0430\u0441\u043f\u043e\u043b\u0430\u0433\u0430\u043b\u0438\u0441\u044c \u043b\u0435\u0432\u0435\u0435 \u0431\u043e\u043b\u044c\u0448\u0438\u0445. \u0412\u043e\u0442 \u0447\u0442\u043e \u044f \u0438\u043c\u0435\u044e \u0432\u0432\u0438\u0434\u0443: <\/p>\n
<\/figure>\n
 \u0412\u043e\u0442 \u0442\u0435\u043f\u0435\u0440\u044c \u0432\u043e \u0432\u0441\u0451\u043c \u0434\u0438\u0430\u043f\u0430\u0437\u043e\u043d\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u0442\u0430\u0440\u0448\u0435\u0435 \u0441\u043b\u043e\u0432\u043e \u0431\u0443\u0434\u0435\u0442 \u043c\u043e\u043d\u043e\u0442\u043e\u043d\u043d\u043e \u0440\u0430\u0441\u0442\u0438 \u0438 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0430 \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c (\u0432 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435 \u0435\u0441\u0442\u044c \u0440\u0430\u0437\u043d\u044b\u0435 \u0441\u043f\u043e\u0441\u043e\u0431\u044b \u0434\u043e\u0441\u0442\u0438\u0447\u044c \u044d\u0442\u043e\u0433\u043e \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u0430, \u043c\u043d\u0435 \u043f\u043e\u043a\u0430\u0437\u0430\u043b\u043e\u0441\u044c \u043f\u0440\u043e\u0449\u0435 \u0432\u043e\u0442 \u0442\u0430\u043a). <\/p>\n
\u0418\u0442\u0430\u043a, \u0431\u0443\u0434\u0435\u043c \u043f\u0438\u043b\u0438\u0442\u044c \u0432\u043e\u0442 \u0442\u0430\u043a\u0443\u044e \u0444\u0443\u043d\u043a\u0446\u0438\u044e:<\/p>\n
LIBMED_API int MedianDBL(double* ptr, size_t Length, double* med) <\/code><\/pre>\n
\u0413\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0430 \u0431\u0443\u0434\u0435\u0442 \u043f\u043e-\u043f\u0440\u0435\u0436\u043d\u0435\u043c\u0443 64\u041a \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432:<\/p>\n
size_t* histo = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t)); <\/code><\/pre>\n
\u0421\u0430\u043c\u043e\u0435 \u0433\u043b\u0430\u0432\u043d\u043e\u0435 — \u0431\u0435\u0440\u0451\u043c \u0443\u043a\u0430\u0437\u0430\u0442\u0435\u043b\u044c \u043d\u0430 \u0432\u0445\u043e\u0434\u043d\u043e\u0439 \u043c\u0430\u0441\u0441\u0438\u0432 \u0438 \u0433\u043e\u0432\u043e\u0440\u0438\u043c, \u0447\u0442\u043e \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u0431\u0443\u0434\u0435\u043c \u0441 \u0446\u0435\u043b\u044b\u043c\u0438:<\/p>\n
uint64_t* ptr_U64 = (uint64_t*)ptr; \/\/Cast <\/code><\/pre>\n
\u0427\u0442\u043e\u0431\u044b \u0443 \u043c\u0435\u043d\u044f \u043d\u0435 \u0440\u044f\u0431\u0438\u043b\u043e \u0432 \u0433\u043b\u0430\u0437\u0430\u0445 \u043e\u0442 \u043d\u0443\u043b\u0435\u0439, \u0441\u0434\u0435\u043b\u0430\u044e \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e «\u0443\u043f\u0440\u043e\u0449\u0430\u044e\u0449\u0438\u0445» \u043f\u043e\u0434\u043c\u0435\u043d, \u0434\u0430 \u043f\u0440\u043e\u0441\u0442\u044f\u0442 \u043c\u0435\u043d\u044f \u0432\u0441\u0435, \u043a\u0442\u043e \u0431\u0443\u0434\u0435\u0442 \u043f\u0440\u0430\u0432\u0438\u0442\u044c \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u043f\u043e\u0441\u043b\u0435 \u043c\u0435\u043d\u044f:<\/p>\n
#define xFFFF 0xFFFFFFFFFFFFFFFF #define x8000 0x8000000000000000 #define x0000 0x0000000000000000 #define xF000 0xFFFF000000000000 #define x0F00 0x0000FFFF00000000 #define x00F0 0x00000000FFFF0000 #define x000F 0x000000000000FFFF  \/\/0 - LSW; 3 - MSW #defineW3 (ptr_U64[i] & 0xFFFF000000000000) #defineW2 (ptr_U64[i] & 0x0000FFFF00000000) #defineW1 (ptr_U64[i] & 0x00000000FFFF0000) #defineW0 (ptr_U64[i] & 0x000000000000FFFF) <\/code><\/pre>\n
\u0415\u0449\u0451 \u043c\u043d\u0435 \u043f\u043e\u043d\u0430\u0434\u043e\u0431\u044f\u0442\u0441\u044f \u043f\u0440\u043e\u0441\u0442\u0435\u043d\u044c\u043a\u0438\u0435 \u044d\u043d\u043a\u043e\u0434\u0435\u0440 \u0438 \u0434\u0435\u043a\u043e\u0434\u0435\u0440:<\/p>\n
#define ENCODE(var) if (!((var) & x8000)) (var) ^= x8000; else (var) ^= xFFFF #define DECODE(var) if (((var) & x8000)) (var) ^= x8000; else (var) ^= xFFFF <\/code><\/pre>\n
\u042d\u043d\u043a\u043e\u0434\u0435\u0440 \u0431\u0443\u0434\u0435\u0442 \u0438\u043d\u0432\u0435\u0440\u0442\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u0441\u0442\u0430\u0440\u0448\u0438\u0439 \u0431\u0438\u0442, \u0435\u0441\u043b\u0438 \u043e\u043d \u043d\u0435 \u0443\u0441\u0442\u0430\u043d\u043e\u0432\u043b\u0435\u043d, \u0432 \u043f\u0440\u043e\u0442\u0438\u0432\u043d\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435 \u0438\u043d\u0432\u0435\u0440\u0442\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u0432\u043e\u043e\u0431\u0449\u0435 \u0432\u0441\u0451. \u0414\u0435\u043a\u043e\u0434\u0435\u0440 \u0431\u0443\u0434\u0435\u0442 \u0432\u043e\u0437\u0432\u0440\u0430\u0449\u0430\u0442\u044c \u0432\u0441\u0451 \u043e\u0431\u0440\u0430\u0442\u043d\u043e.<\/p>\n
\u0414\u0435\u043b\u0430\u0435\u043c \u043f\u0435\u0440\u0432\u044b\u0439 \u043f\u0440\u043e\u0445\u043e\u0434 — \u0433\u043e\u0442\u043e\u0432\u0438\u043c \u0434\u0430\u043d\u043d\u044b\u0435 \u0438 \u0441\u043e\u0431\u0438\u0440\u0430\u0435\u043c \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0438\u0437 \u0441\u0442\u0430\u0440\u0448\u0435\u0433\u043e \u0441\u043b\u043e\u0432\u0430:<\/p>\n
for (i = 0; i < Length; i++) { ENCODE(ptr_U64[i]); histo[W3 >> 48]++; } <\/code><\/pre>\n
\u041c\u0435\u043d\u044f\u0442\u044c \u0432\u0445\u043e\u0434\u043d\u044b\u0435 \u0434\u0430\u043d\u043d\u044b\u0435 \u0432 \u043e\u0431\u0449\u0435\u043c \u2014 \u043d\u0443 \u0442\u0430\u043a \u0441\u0435\u0431\u0435 \u0438\u0434\u0435\u044f, \u043d\u043e \u043c\u044b \u0432 \u043a\u043e\u043d\u0446\u0435 \u043f\u043e\u0442\u0438\u0445\u043e\u043d\u044c\u043a\u0443 \u0432\u0435\u0440\u043d\u0451\u043c \u0432\u0441\u0451 \u043e\u0431\u0440\u0430\u0442\u043d\u043e, \u0430\u0432\u043e\u0441\u044c \u043d\u0438\u043a\u0442\u043e \u043d\u0435 \u0437\u0430\u043c\u0435\u0442\u0438\u0442 (\u0434\u043e \u0442\u0435\u0445 \u043f\u043e\u0440 \u043f\u043e\u043a\u0430 \u043d\u0438\u043a\u0442\u043e \u043d\u0435 \u0431\u0443\u0434\u0435\u0442 \u0447\u0438\u0442\u0430\u0442\u044c \u044d\u0442\u043e\u0442 \u043c\u0430\u0441\u0441\u0438\u0432 \u0432\u043e \u0432\u0440\u0435\u043c\u044f \u0432\u044b\u043f\u043e\u043b\u043d\u0435\u043d\u0438\u044f). \u0423 numpy, \u043a\u0441\u0442\u0430\u0442\u0438, \u044f \u0432\u0438\u0434\u0435\u043b \u0444\u043b\u0430\u0433 overwrite_input \u2014 \u043c\u043e\u0436\u043d\u043e \u0438\u043b\u0438 \u043d\u0435\u0442 \u043c\u043e\u0434\u0438\u0444\u0438\u0446\u0438\u0440\u043e\u0432\u0430\u0442\u044c \u0432\u0445\u043e\u0434\u044f\u0449\u0438\u0439 \u043c\u0430\u0441\u0441\u0438\u0432.<\/p>\n
\u0414\u0435\u043b\u0430\u0435\u043c \u043f\u0435\u0440\u0432\u044b\u0439 \u043f\u0440\u043e\u0445\u043e\u0434 \u043f\u043e \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0435, \u043e\u0441\u0442\u0430\u043d\u0430\u0432\u043b\u0438\u0432\u0430\u0435\u043c\u0441\u044f, \u043a\u043e\u0433\u0434\u0430 \u043d\u0430\u0431\u0435\u0440\u0451\u043c \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u0443:<\/p>\n
medianPos = (Length + 1) \/ 2; for (count = 0, i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram count += histo[i]; if (count >= medianPos) { medValW3 = i; break; } \/\/ >> single Median at MSW } <\/code><\/pre>\n
\u0412 \u0441\u043b\u0443\u0447\u0430\u0435 \u0447\u0451\u0442\u043d\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0430 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0438 \u0442\u043e\u0447\u043d\u043e\u0433\u043e \u0441\u043e\u0432\u043f\u0430\u0434\u0435\u043d\u0438\u044f \u043f\u0440\u043e\u043c\u0430\u0442\u044b\u0432\u0430\u0435\u043c \u0434\u043e \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0435\u0433\u043e \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f:<\/p>\n
if (!(Length & 1) && count == medianPos) { \/\/W3 nextVal = medValW3; while (!histo[++nextVal]); medValW3 <<= 48; nextVal <<= 48; <\/code><\/pre>\n
\u0422\u0435\u043f\u0435\u0440\u044c \u0432\u0442\u043e\u0440\u043e\u0439 \u043f\u0440\u043e\u0445\u043e\u0434: \u0438\u0449\u0435\u043c \u043c\u0438\u043d\u0438\u043c\u0443\u043c \u0438 \u043c\u0430\u043a\u0441\u0438\u043c\u0443\u043c \u0432\u043e \u0432\u0442\u043e\u0440\u043e\u043c \u0441\u043b\u043e\u0432\u0435, \u0442\u0430\u043c, \u0433\u0434\u0435 \u0441\u0442\u0430\u0440\u0448\u0435\u0435 \u0441\u043e\u0432\u043f\u0430\u0434\u0430\u0435\u0442 \u0441 \u043d\u0430\u0439\u0434\u0435\u043d\u043d\u044b\u043c\u0438 \u0432\u044b\u0448\u0435:<\/p>\n
for (max2 = 0, min2 = _UI64_MAX, i = 0; i < Length; i++) { \/\/1.2 pass if (medValW3 == W3 && W2 > max2) max2 = W2; if (W3 == nextVal && W2 < min2) min2 = W2; } <\/code><\/pre>\n
\u041f\u043e\u0432\u0442\u043e\u0440\u044f\u0435\u043c \u0434\u043b\u044f \u0434\u0432\u0443\u0445 \u043e\u0441\u0442\u0430\u0432\u0448\u0438\u0445\u0441\u044f \u0441\u043b\u043e\u0432<\/p>\n
for (max1 = 0, min1 = _UI64_MAX, i = 0; i < Length; i++) { \/\/1.3 pass if (medValW3 == W3 && max2 == W2 && W1 > max1) max1 = W1; if (W3 == nextVal && min2 == W2 && W1 < min1) min1 = W1; }  for (max0 = 0, min0 = _UI64_MAX, i = 0; i < Length; i++) { \/\/1.4 pass, final if (medValW3 == W3 && max2 == W2 && max1 == W1 && W0 > max0) max0 = W0; if (W3 == nextVal && min2 == W2 && min1 == W1 && W0 < min0) min0 = W0; DECODE(ptr_U64[i]); } \/\/ 1.2 - 1.4 <\/code><\/pre>\n
\u0412 \u043a\u043e\u043d\u0446\u0435 \u043d\u0435 \u0437\u0430\u0431\u044b\u0432\u0430\u0435\u043c \u0432\u0435\u0440\u043d\u0443\u0442\u044c \u043c\u0430\u0441\u0441\u0438\u0432 \u0432 \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u0435 \u0441\u043e\u0441\u0442\u043e\u044f\u043d\u0438\u0435.<\/p>\n
\u041d\u0443 \u0438 «\u0441\u043e\u0431\u0438\u0440\u0430\u0435\u043c» \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u0438\u0437 \u0434\u0432\u0443\u0445 \u0447\u0430\u0441\u0442\u0435\u0439:<\/p>\n
medianLeft = medValW3 | max2 | max1 | max0; medianRight = nextVal | min2 | min1 | min0; DECODE(medianLeft); DECODE(medianRight); *med = (*(double*)&medianLeft + *(double*)&medianRight) \/ 2.0; <\/code><\/pre>\n
\u0410 \u0432 \u0441\u043b\u0443\u0447\u0430\u0435 \u043d\u0435\u0447\u0451\u0442\u043d\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0430 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0441\u043e\u0431\u0438\u0440\u0430\u0435\u043c \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0443\u044e \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443:<\/p>\n
memset(histo, 0, (USHRT_MAX + 1) * sizeof(size_t)); med_ = (medValW3 <<= 48); for (count = 0, i = 0; i < Length; i++) { \/\/2.2 pass if (medValW3 == W3) histo[W2 >> 32]++; if (ptr_U64[i] < med_) count++; } \/\/ 2.2  for (int i = 0; i <= USHRT_MAX; i++) { count += histo[i]; \/\/continue if (count >= medianPos) { medValW2 = i; break; } } <\/code><\/pre>\n
\u041d\u0443 \u0438 \u0442\u0430\u043a \u0434\u0430\u043b\u0435\u0435, \u0442\u0435\u043f\u0435\u0440\u044c \u0441\u043e\u0431\u0438\u0440\u0430\u0435\u043c \u0432\u0441\u0451 \u0432\u043c\u0435\u0441\u0442\u0435 \u043c\u0435\u0442\u043e\u0434\u043e\u043c \u043e\u0442\u0447\u0430\u044f\u043d\u043d\u043e\u0433\u043e \u043a\u043e\u043f\u0438\u043f\u0430\u0441\u0442\u0438\u043d\u0433\u0430 (\u044f \u043e\u0442\u043d\u044e\u0434\u044c \u043d\u0435 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u043e \u0434\u043e\u0431\u0430\u0432\u0438\u043b \u0445\u0430\u0431 «\u043d\u0435\u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0435»):<\/p>\n
\n\u043f\u043e\u043b\u043d\u044b\u0439 \u043b\u0438\u0441\u0442\u0438\u043d\u0433 MedianDBL<\/summary>\n
\n
#define xFFFF 0xFFFFFFFFFFFFFFFF #define x8000 0x8000000000000000 #define x0000 0x0000000000000000 #define xF000 0xFFFF000000000000 #define x0F00 0x0000FFFF00000000 #define x00F0 0x00000000FFFF0000 #define x000F 0x000000000000FFFF  \/\/0 - LSW; 3 - MSW #defineW3 (ptr_U64[i] & 0xFFFF000000000000) #defineW2 (ptr_U64[i] & 0x0000FFFF00000000) #defineW1 (ptr_U64[i] & 0x00000000FFFF0000) #defineW0 (ptr_U64[i] & 0x000000000000FFFF)  #define ENCODE(var) if (!((var) & x8000)) (var) ^= x8000; else (var) ^= xFFFF #define DECODE(var) if (((var) & x8000)) (var) ^= x8000; else (var) ^= xFFFF  LIBMED_API int MedianDBL(double* ptr, size_t Length, double* med) { if (!ptr || !Length) return -1; size_t* histo = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t)); if (!histo) return -2; uint64_t count, i, medianPos; uint64_t min0, min1, min2, max0, max1, max2; uint64_t medValW0 = 0, medValW1 = 0, medValW2 = 0, medValW3 = 0;  uint64_t medianLeft, medianRight, medVal, nextVal, med_;   uint64_t* ptr_U64 = (uint64_t*)ptr; \/\/Cast  for (i = 0; i < Length; i++) { ENCODE(ptr_U64[i]); histo[W3 >> 48]++; } \/\/ Find the initial median position; medianPos = (Length + 1) \/ 2; for (count = 0, i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram count += histo[i]; if (count >= medianPos) { medValW3 = i; break; } \/\/ >> single Median at MSW }  if (!(Length & 1) && count == medianPos) { \/\/W3 nextVal = medValW3; while (!histo[++nextVal]); medValW3 <<= 48; nextVal <<= 48;  for (max2 = 0, min2 = _UI64_MAX, i = 0; i < Length; i++) { \/\/1.2 pass if (medValW3 == W3 && W2 > max2) max2 = W2; if (W3 == nextVal && W2 < min2) min2 = W2; }  for (max1 = 0, min1 = _UI64_MAX, i = 0; i < Length; i++) { \/\/1.3 pass if (medValW3 == W3 && max2 == W2 && W1 > max1) max1 = W1; if (W3 == nextVal && min2 == W2 && W1 < min1) min1 = W1; }  for (max0 = 0, min0 = _UI64_MAX, i = 0; i < Length; i++) { \/\/1.4 pass, final if (medValW3 == W3 && max2 == W2 && max1 == W1 && W0 > max0) max0 = W0; if (W3 == nextVal && min2 == W2 && min1 == W1 && W0 < min0) min0 = W0; DECODE(ptr_U64[i]); } \/\/ 1.2 - 1.4  medianLeft = medValW3 | max2 | max1 | max0; DECODE(medianLeft); medianRight = nextVal | min2 | min1 | min0; DECODE(medianRight); *med = (*(double*)&medianLeft + *(double*)&medianRight) \/ 2.0; } else { memset(histo, 0, (USHRT_MAX + 1) * sizeof(size_t)); med_ = (medValW3 <<= 48); for (count = 0, i = 0; i < Length; i++) { \/\/2.2 pass if (medValW3 == W3) histo[W2 >> 32]++; if (ptr_U64[i] < med_) count++; } \/\/ 2.2  for (int i = 0; i <= USHRT_MAX; i++) { count += histo[i]; \/\/continue if (count >= medianPos) { medValW2 = i; break; } }  if (!(Length & 1) && count == medianPos) { \/\/ W2 nextVal = medValW2; while (!histo[++nextVal]); medValW2 <<= 32; nextVal <<= 32;  for (max1 = 0, min1 = _UI64_MAX, i = 0; i < Length; i++) { \/\/2.3 pass if (medValW3 == W3 && medValW2 == W2 && W1 > max1) max1 = W1; if (W3 == medValW3 && nextVal == W2 && W1 < min1) min1 = W1; }  for (max0 = 0, min0 = _UI64_MAX, i = 0; i < Length; i++) { \/\/2.4 pass, final if (medValW3 == W3 && medValW2 == W2 && max1 == W1 && W0 > max0) max0 = W0; if (W3 == medValW3 && nextVal == W2 && min1 == W1 && W0 < min0) min0 = W0; DECODE(ptr_U64[i]); } \/\/2.3-2.4  medianLeft = medValW3 | medValW2 | max1 | max0; DECODE(medianLeft); medianRight = medValW3 | nextVal | min1 | min0; DECODE(medianRight); *med = (*(double*)&medianLeft + *(double*)&medianRight) \/ 2.0; } else { memset(histo, 0, (USHRT_MAX + 1) * sizeof(size_t)); med_ = medValW3 | (medValW2 <<= 32); for (count = 0, i = 0; i < Length; i++) { \/\/3.3 pass if (medValW3 == W3 && medValW2 == W2) histo[W1 >> 16]++; if (ptr_U64[i] < med_) count++; } \/\/3.3  for (int i = 0; i <= USHRT_MAX; i++) { count += histo[i]; if (count >= medianPos) { medValW1 = i; break; } }  if (!(Length & 1) && count == medianPos) { \/\/W1 nextVal = medValW1; while (!histo[++nextVal]); medValW1 <<= 16; nextVal <<= 16;  for (max0 = 0, min0 = _UI64_MAX, i = 0; i < Length; i++) { \/\/3.4 pass, final if ((medValW3 == W3) && (medValW2 == W2) && (medValW1 == W1) && W0 > max0) max0 = W0; if ((W3 == medValW3) && (medValW2 == W2) && (nextVal == W1) && W0 < min0) min0 = W0; DECODE(ptr_U64[i]); } \/\/3.4  medianLeft = medValW3 | medValW2 | medValW1 | max0; DECODE(medianLeft); medianRight = medValW3 | medValW2 | nextVal | min0; DECODE(medianRight); *med = (*(double*)&medianLeft + *(double*)&medianRight) \/ 2.0; } else { memset(histo, 0, (USHRT_MAX + 1) * sizeof(size_t)); med_ = medValW3 | medValW2 | (medValW1 <<= 16); for (count = 0, i = 0; i < Length; i++) {\/\/!! Spaziergang !! if (medValW3 == W3 && medValW2 == W2 && medValW1 == W1) histo[W0]++; if (ptr_U64[i] < med_) count++; DECODE(ptr_U64[i]); \/\/restore back } \/\/4.4, final  for (int i = 0; i <= USHRT_MAX; i++) { count += histo[i]; if (count >= medianPos) { medValW0 = i; break; } }  if (!(Length & 1) && count == medianPos) { nextVal = medValW0; while (!histo[++nextVal]); medianLeft = medValW3 | medValW2 | medValW1 | medValW0; DECODE(medianLeft); medianRight = medValW3 | medValW2 | medValW1 | nextVal; DECODE(medianRight); *med = (*(double*)&medianLeft + *(double*)&medianRight) \/ 2.0; } else { medVal = medValW3 | medValW2 | medValW1 | medValW0; DECODE(medVal); *med = *(double*)&medVal; } } } } free(histo); return 0; } \/\/MedianDBL  <\/code><\/pre>\n<\/div>\n<\/details>\n
 \u0412\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u0434\u043e\u0432\u043e\u043b\u044c\u043d\u043e \u0441\u0443\u0440\u043e\u0432\u043e, \u043d\u043e \u0435\u0441\u043b\u0438 \u0441\u0440\u0430\u0432\u043d\u0438\u0442\u044c \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u0441 «\u043a\u043b\u0430\u0441\u0441\u0438\u043a\u043e\u0439» \u0442\u0438\u043f\u0430 QuickSelect \u0438\u0437 «\u0427\u0438\u0441\u043b\u0435\u043d\u043d\u044b\u0445 \u043c\u0435\u0442\u043e\u0434\u043e\u0432», \u0442\u043e \u0432 \u043f\u0440\u0435\u0434\u043b\u043e\u0436\u0435\u043d\u043d\u043e\u043c \u043f\u043e\u0434\u0445\u043e\u0434\u0435 \u043c\u043d\u043e\u0433\u043e\u0435 \u0434\u0430\u0436\u0435 \u043f\u0440\u043e\u0449\u0435. \u0418 \u0434\u0430, \u044f \u0443\u0447\u0438\u043b \u0421\u0438 \u043f\u043e \u041a\u0435\u0440\u043d\u0438\u0433\u0430\u043d\u0443 \u0438 \u0420\u0438\u0442\u0447\u0438. <\/p>\n
\u041f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u044c<\/h3>\n
\u041d\u0443 \u0432\u043e\u0442, \u0442\u0435\u043f\u0435\u0440\u044c \u043c\u043e\u0436\u043d\u043e \u0438 \u043f\u043e\u0433\u043e\u043d\u044f\u0442\u044c \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u0438 \u0441\u0440\u0430\u0432\u043d\u0438\u0442\u044c \u0441 \u0434\u0440\u0443\u0433\u0438\u043c\u0438 \u0440\u0435\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u044f\u043c\u0438.<\/p>\n
\u0423 \u043c\u0435\u043d\u044f \u0432 \u0440\u0443\u043a\u0430\u0445 \u0435\u0441\u0442\u044c \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0434\u0438\u0430\u043d \u0438\u0437 \u0434\u0432\u0443\u0445 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0438 — \u043e\u0434\u043d\u0430 \u044d\u0442\u043e Numpy Median<\/a>, \u0430 \u0432\u0442\u043e\u0440\u0430\u044f — NI Median<\/a>.<\/p>\n
\u0412\u043e\u0442 \u0441 \u043d\u0438\u043c\u0438 \u043c\u044b \u0438 \u0431\u0443\u0434\u0435\u043c \u0441\u043e\u0440\u0435\u0432\u043d\u043e\u0432\u0430\u0442\u044c\u0441\u044f.<\/p>\n
\u041d\u0430\u0447\u043d\u0443 \u044f \u0441 \u041f\u0438\u0442\u043e\u043d\u0430, \u043f\u043e\u0436\u0430\u043b\u0443\u0439, \u0432 \u0441\u0438\u043b\u0443 \u0437\u043d\u0430\u0447\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u0431\u043e\u043b\u044c\u0448\u0435\u0439 \u0440\u0430\u0441\u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0451\u043d\u043d\u043e\u0441\u0442\u0438 (\u043f\u0440\u043e\u0434\u0443\u043a\u0442\u044b NI — \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e\u0442\u0441\u044f \u043d\u0435 \u0442\u0430\u043a \u0447\u0430\u0441\u0442\u043e). \u042f \u043e\u0444\u043e\u0440\u043c\u043b\u044e \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u0432\u044b\u0448\u0435 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u0447\u043a\u043e\u0439 libMed.dll \u0438 \u0432\u044b\u0437\u043e\u0432\u0443 \u0438\u0437 \u041f\u0438\u0442\u043e\u043d\u0430 \u0432\u043e\u0442 \u0442\u0430\u043a (\u044f \u043d\u0438 \u0440\u0430\u0437\u0443 \u043d\u0435 \u043f\u0440\u043e\u0444\u0438 \u0432 \u041f\u0438\u0442\u043e\u043d\u0435, \u0438 \u0432\u044b\u0437\u044b\u0432\u0430\u044e DLL \u0438\u0437 \u043d\u0435\u0433\u043e \u0432 \u043f\u0435\u0440\u0432\u044b\u0439 \u0440\u0430\u0437 \u0432 \u0436\u0438\u0437\u043d\u0438, \u0442\u0430\u043a \u0447\u0442\u043e \u043f\u043e\u043b\u0443\u0447\u0438\u043b\u043e\u0441\u044c \u043a\u0430\u043a \u043f\u043e\u043b\u0443\u0447\u0438\u043b\u043e\u0441\u044c):<\/p>\n
import os from random import random from time import time from ctypes import CDLL, POINTER, byref, c_size_t, c_double  mylib = CDLL(os.path.dirname(os.path.abspath(__file__)) + os.path.sep + \"libMed.dll\") ND_POINTER_1 = np.ctypeslib.ndpointer(dtype=np.float64, ndim=1, flags=\"C\") mylib.fnMedian.argtypes = [ND_POINTER_1, c_size_t, POINTER( c_double )] mylib.fnMedian.restype = np.int32 # Test Data T1M = np.array([random() for _ in range(1024*1024)]) MyRes = c_double( 0.0 )  start = time() mylib.fnMedian(T1M, T1M.size, byref(MyRes)) print (\"My 1M:\", (time() - start)*1000, \"ms; res = \", MyRes) #... <\/code><\/pre>\n
(\u0433\u0434\u0435-\u0442\u043e \u0432\u044b \u0431\u0443\u0434\u0435\u0442\u0435 \u0432\u0441\u0442\u0440\u0435\u0447\u0430\u0442\u044c \u0438\u043c\u044f fnMedian, \u0430 \u043d\u0435 MedianDBL, \u044d\u0442\u043e \u044f \u043f\u0440\u043e\u0441\u0442\u043e \u0434\u0435\u043a\u043e\u0440\u0430\u0446\u0438\u044e \u0438\u043c\u0451\u043d \u0432\u044b\u043a\u043b\u044e\u0447\u0438\u043b, \u0434\u0430 \u0442\u0430\u043a \u043e\u043d\u043e \u0438 \u043e\u0441\u0442\u0430\u043b\u043e\u0441\u044c).<\/p>\n
\u042f \u0432 \u043a\u0443\u0440\u0441\u0435 \u043d\u0435\u0432\u044b\u0441\u043e\u043a\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438 time(), \u043d\u043e \u043c\u043d\u0435 \u0438 \u043d\u0435 \u043d\u0443\u0436\u043d\u043e \u0433\u043d\u0430\u0442\u044c\u0441\u044f \u0437\u0430 \u043c\u0438\u043a\u0440\u043e\u0441\u0435\u043a\u0443\u043d\u0434\u0430\u043c\u0438, \u0432 \u043a\u043e\u043d\u0446\u0435 \u043a\u043e\u043d\u0446\u043e\u0432 \u044f \u043c\u043e\u0433\u0443 \u043f\u0440\u043e\u043a\u0440\u0443\u0442\u0438\u0442\u044c \u043a\u043e\u0434 \u0434\u0435\u0441\u044f\u0442\u043e\u043a \u0440\u0430\u0437 \u0434\u0430 \u043f\u043e\u043b\u0443\u0447\u0438\u0442\u044c \u0441\u0440\u0435\u0434\u043d\u0435\u0435. \u042f \u043f\u0440\u043e\u0432\u0435\u0440\u044e \u043d\u0430 \u0442\u0440\u0451\u0445 \u0441\u0430\u043c\u044b\u0445 «\u0445\u043e\u0434\u043e\u0432\u044b\u0445» (\u0432 \u043c\u043e\u0451\u043c \u0447\u0430\u0441\u0442\u043d\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435) \u0440\u0430\u0437\u043c\u0435\u0440\u0430\u0445 \u043a\u0430\u0440\u0442\u0438\u043d\u043e\u043a \u2014 1024\u04451024, 2048\u04452048 \u0438 4096\u04454096, \u043d\u0443 \u0442\u043e \u0435\u0441\u0442\u044c 1, 4 \u0438 16 \u043c\u0438\u043b\u043b\u0438\u043e\u043d\u043e\u0432 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0432 \u043c\u0430\u0441\u0441\u0438\u0432\u0435.<\/p>\n
\n\u043f\u043e\u043b\u043d\u044b\u0439 \u043b\u0438\u0441\u0442\u0438\u043d\u0433 \u0442\u0435\u0441\u0442\u0430 \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438<\/summary>\n
\n
import os import numpy as np from random import random from time import time from ctypes import CDLL, POINTER, byref, c_size_t, c_double  mylib = CDLL(os.path.dirname(os.path.abspath(__file__)) + os.path.sep + \"libMed.dll\") ND_POINTER_1 = np.ctypeslib.ndpointer(dtype=np.float64, ndim=1, flags=\"C\") mylib.fastMedian.argtypes = [ND_POINTER_1, c_size_t, POINTER( c_double )] mylib.fastMedian.restype = np.int32 # Test Data T1M = np.array([random() for _ in range(1024*1024+1)]) T4M = np.array([random() for _ in range(2048*2048+1)]) T16M = np.array([random() for _ in range(4096*4096+1)]) MyRes = c_double( 0.0 )  start = time() for x in range(10):mylib.fastMedian(T1M, T1M.size, byref(MyRes)) print (\"My 1M:\", (time() - start)*100, \"ms; res = \", MyRes)                    start = time()           for x in range(10):NumpyRes = np.median(T1M) print (\"np 1M:\", (time() - start)*100, \"ms; res = \", NumpyRes) if (MyRes == NumpyRes): print(\"passed\")  start = time() for x in range(10):mylib.fastMedian(T4M, T4M.size, byref(MyRes)) print (\"My 4M:\", (time() - start)*100, \"ms; res = \", MyRes)  start = time()           for x in range(10):NumpyRes = np.median(T4M) print (\"np 4M:\", (time() - start)*100, \"ms; res = \", NumpyRes) if (MyRes == NumpyRes): print(\"passed\")  start = time() for x in range(10):mylib.fastMedian(T16M, T16M.size, byref(MyRes)) print (\"My 16M:\", (time() - start)*100, \"ms; res = \", MyRes)  start = time()           for x in range(10):NumpyRes = np.median(T16M) print (\"np 16M:\", (time() - start)*100, \"ms; res = \", NumpyRes) if (MyRes == NumpyRes): print(\"passed\")   <\/code><\/pre>\n<\/div>\n<\/details>\n
\u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b (My — \u044d\u0442\u043e \u043c\u043e\u0439 \u043c\u0435\u0442\u043e\u0434, \u0430 np \u2014 \u044d\u0442\u043e numpy):<\/p>\n
>python fastMedTest.py My 1M: 9.675335884094238 ms; res =  c_double(0.49913991096863186) np 1M: 17.79000759124756 ms; res =  0.49913991096863186 passed My 4M: 38.81254196166992 ms; res =  c_double(0.5001439678114372) np 4M: 69.83757019042969 ms; res =  0.5001439678114372 passed My 16M: 152.37784385681152 ms; res =  c_double(0.5000523968034343) np 16M: 259.2304468154907 ms; res =  0.5000523968034343 passed <\/code><\/pre>\n
\u0412 \u043c\u043e\u0451\u043c Precision M6700, \u043d\u0430 \u043a\u043e\u0442\u043e\u0440\u043e\u043c \u044f \u0443\u043f\u0440\u0430\u0436\u043d\u044f\u044e\u0441\u044c, \u043a\u0430\u043c\u0443\u0448\u0435\u043a i7-3740QM, \u043f\u043e\u044d\u0442\u043e\u043c\u0443 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b \u0442\u0430\u043a \u0441\u0435\u0431\u0435. \u0410 \u043d\u043e\u0443\u0442, \u043a\u0441\u0442\u0430\u0442\u0438, \u0441\u0447\u0438\u0442\u0430\u044e \u043e\u0434\u043d\u0438\u043c \u0438\u0437 \u043b\u0443\u0447\u0448\u0438\u0445 \u0432\u0441\u0435\u0445 \u0432\u0440\u0435\u043c\u0451\u043d (\u0432\u043f\u0435\u0440\u0435\u0434\u0438 \u0442\u043e\u043b\u044c\u043a\u043e \u043e\u0442\u0434\u0435\u043b\u044c\u043d\u044b\u0435 ThinkPad\u044b), \u0435\u0441\u043b\u0438 \u0431 \u043d\u0435 \u0432\u0435\u0441 \u043f\u043e\u0434 \u0447\u0435\u0442\u044b\u0440\u0435 \u043a\u0438\u043b\u043e.<\/p>\n
\u0415\u0441\u043b\u0438 \u043f\u0435\u0440\u0435\u0441\u0435\u0441\u0442\u044c \u043d\u0430 i7-7700 \u043f\u043e\u0441\u0432\u0435\u0436\u0435\u0435, \u0442\u043e \u0441\u0442\u0430\u043d\u0435\u0442 \u0447\u0443\u0442\u044c \u0432\u0435\u0441\u0435\u043b\u0435\u0435:<\/p>\n
>python fastMedTest.py  My 1M: 4.303789138793945 ms; res =  c_double(0.5003872509309095) np 1M: 10.99863052368164 ms; res =  0.5003872509309095 passed My 4M: 18.0020809173584 ms; res =  c_double(0.4999627738635442) np 4M: 46.71292304992676 ms; res =  0.4999627738635442 passed My 16M: 78.09562683105469 ms; res =  c_double(0.5001147029114446) np 16M: 167.00820922851562 ms; res =  0.5001147029114446 passed <\/code><\/pre>\n
\u041d\u043e \u0432 \u043b\u044e\u0431\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435 \u0432 \u043d\u0430\u043b\u0438\u0447\u0438\u0438 \u043f\u0440\u0430\u043a\u0442\u0438\u0447\u0435\u0441\u043a\u0438 \u0434\u0432\u0443\u043a\u0440\u0430\u0442\u043d\u044b\u0439 \u0432\u044b\u0438\u0433\u0440\u044b\u0448 \u043f\u043e \u0432\u0440\u0435\u043c\u0435\u043d\u0438, \u043f\u0440\u0438 \u044d\u0442\u043e\u043c \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u0432 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438 \u0441\u043e\u0432\u043f\u0430\u0434\u0430\u0435\u0442. \u0412 \u0441\u043b\u0443\u0447\u0430\u0435, \u043a\u043e\u0433\u0434\u0430 \u0434\u0430\u043d\u043d\u044b\u0435 \u0431\u0443\u0434\u0443\u0442 \u0441\u043e\u0434\u0435\u0440\u0436\u0430\u0442\u044c \u0442\u043e\u043b\u044c\u043a\u043e \u043f\u043e\u043b\u043e\u0436\u0438\u0442\u0435\u043b\u044c\u043d\u044b\u0435 \u0447\u0438\u0441\u043b\u0430, \u0430 \u0442\u0435\u043c \u0431\u043e\u043b\u0435\u0435 \u0432 \u0441\u043b\u0443\u0447\u0430\u0435 \u043e\u0434\u0438\u043d\u0430\u0440\u043d\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438 (\u043a\u043e\u0433\u0434\u0430 \u043c\u044b \u043d\u0430\u0439\u0434\u0451\u043c \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u0432\u0441\u0435\u0433\u043e \u0437\u0430 \u0434\u0432\u0430 \u043f\u0440\u043e\u0445\u043e\u0434\u0430) \u0432\u044b\u0438\u0433\u0440\u044b\u0448 \u0431\u0443\u0434\u0443\u0442 \u0435\u0449\u0451 \u0437\u043d\u0430\u0447\u0438\u0442\u0435\u043b\u044c\u043d\u0435\u0435. \u042d\u0442\u0443 \u0436\u0435 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0443 \u044f \u043c\u043e\u0433\u0443 \u0438 \u0438\u0437 LabVIEW \u0432\u044b\u0437\u0432\u0430\u0442\u044c \u0438 \u043f\u0440\u043e\u0432\u0435\u0440\u0438\u0442\u044c \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u0438 \u0432\u0440\u0435\u043c\u044f: <\/p>\n
<\/figure>\n
 \u0411\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 \u043e\u0442 NI \u0434\u0430\u0451\u0442 \u043f\u0440\u0438\u043c\u0435\u0440\u043d\u043e \u043f\u043e\u0445\u043e\u0436\u0438\u0435 \u0441 numpy \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b. \u0412 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435 \u043d\u0430\u0447\u0438\u043d\u0430\u044f \u0441 LabVIEW 2019 \u043d\u0430\u043c \u0437\u0430\u0432\u0435\u0437\u043b\u0438 \u043f\u043e\u0434\u0434\u0435\u0440\u0436\u043a\u0443 \u041f\u0438\u0442\u043e\u043d\u0430 \u043f\u0440\u044f\u043c\u043e \u0432 LabVIEW \u0438 \u044f \u043c\u043e\u0433\u0443 \u0432\u044b\u0437\u0432\u0430\u0442\u044c \u041c\u0435\u0434\u0438\u0430\u043d\u0443 \u0438\u0437 numpy \u0432\u043e\u0442 \u0442\u0430\u043a: <\/p>\n
<\/figure>\n
 \u041d\u043e \u0434\u043b\u044f \u043e\u0446\u0435\u043d\u043a\u0438 \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438 \u044d\u0442\u043e \u043d\u0435 \u043e\u0447\u0435\u043d\u044c \u043f\u043e\u0434\u0445\u043e\u0434\u0438\u0442, \u043f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u044f \u0432\u0438\u0436\u0443 \u043e\u0447\u0435\u043d\u044c \u0441\u0438\u043b\u044c\u043d\u044b\u0435 \u0437\u0430\u0434\u0435\u0440\u0436\u043a\u0438 \u043f\u0440\u0438 \u043f\u0440\u043e\u0431\u0440\u043e\u0441\u0435 LabVIEW \u043c\u0430\u0441\u0441\u0438\u0432\u0430 \u0432 \u0441\u043a\u0440\u0438\u043f\u0442 \u043d\u0430 \u041f\u0438\u0442\u043e\u043d\u0435. <\/p>\n
\u0412 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435 \u0442\u0430\u043c \u0435\u0441\u0442\u044c \u0435\u0449\u0451 \u0447\u0442\u043e \u043f\u043e\u043e\u043f\u0442\u0438\u043c\u0438\u0437\u0438\u0440\u043e\u0432\u0430\u0442\u044c, \u0447\u0430\u0441\u0442\u044c \u043c\u043e\u0436\u043d\u043e \u043f\u0435\u0440\u0435\u043f\u0438\u0441\u0430\u0442\u044c \u043d\u0430 AVX2\/AVX512, \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u043e, \u0418\u043d\u0442\u0435\u043b\u043e\u0432\u0441\u043a\u0438\u0439 \u043a\u043e\u043c\u043f\u0438\u043b\u044f\u0442\u043e\u0440 \u0441\u043c\u043e\u0436\u0435\u0442 \u0434\u0430\u0442\u044c \u0447\u0443\u0442\u044c \u0431\u043e\u043b\u044c\u0448\u0443\u044e \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u044c. \u0420\u0430\u0437\u0443\u043c\u0435\u0435\u0442\u0441\u044f, \u0432\u044b\u0438\u0433\u0440\u044b\u0448 \u044d\u0442\u043e\u0442 \u043c\u0435\u0442\u043e\u0434 \u0434\u0430\u0441\u0442 \u043d\u0430\u0447\u0438\u043d\u0430\u044f \u0441 \u043e\u043f\u0440\u0435\u0434\u0435\u043b\u0451\u043d\u043d\u043e\u0433\u043e \u0440\u0430\u0437\u043c\u0435\u0440\u0430 \u043c\u0430\u0441\u0441\u0438\u0432\u0430, \u0442\u0430\u043a \u043a\u0430\u043a \u043f\u0440\u043e\u0433\u0443\u043b\u043a\u0438 \u043f\u043e \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0430\u043c \u0441\u0430\u043c\u0438 \u043f\u043e \u0441\u0435\u0431\u0435 \u043d\u0435 \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u044b, \u0438 \u0434\u043e 10000-15000 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043e\u043d \u0431\u0443\u0434\u0435\u0442 \u043c\u0435\u0434\u043b\u0435\u043d\u043d\u0435\u0435, \u043a\u0440\u043e\u043c\u0435 \u0442\u043e\u0433\u043e \u0434\u043e\u0432\u043e\u043b\u044c\u043d\u043e \u0441\u0438\u043b\u044c\u043d\u043e \u0437\u0430\u0432\u0438\u0441\u0435\u0442\u044c \u043e\u0442 \u0432\u0445\u043e\u0434\u043d\u044b\u0445 \u0434\u0430\u043d\u043d\u044b\u0445. \u041d\u043e \u0432 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435 \u0443\u0436\u0435 \u0441 \u0440\u0430\u0437\u043c\u0435\u0440\u043e\u0432 \u0432 \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u0443 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u044b, \u0433\u0434\u0435-\u0442\u043e \u0441 25000 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 — \u043f\u0440\u0438\u0440\u043e\u0441\u0442 \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438 \u0437\u0430\u043c\u0435\u0442\u0435\u043d: <\/p>\n
<\/figure>\n
 (QS — \u044d\u0442\u043e Quick Select \u0438\u0437 «\u0427\u0438\u0441\u043b\u0435\u043d\u043d\u044b\u0445 \u041c\u0435\u0442\u043e\u0434\u043e\u0432»). \u0410 \u044d\u0442\u043e \u0433\u0440\u0430\u0444\u0438\u043a \u0434\u043e \u043c\u0438\u043b\u043b\u0438\u043e\u043d\u0430 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432. \u0412\u0438\u0434\u043d\u043e. \u0447\u0442\u043e NI \u0438 Numpy \u0431\u0435\u0433\u0443\u0442 \u043f\u043e\u0447\u0442\u0438 \u0432\u0440\u043e\u0432\u0435\u043d\u044c (numpy \u0447\u0443\u0442\u044c \u043b\u0443\u0447\u0448\u0435), \u0430 \u043f\u0440\u0435\u0434\u043b\u043e\u0436\u0435\u043d\u043d\u044b\u0439 \u043c\u0435\u0442\u043e\u0434 \u0443\u0432\u0435\u0440\u0435\u043d\u043d\u043e \u043b\u0438\u0434\u0438\u0440\u0443\u0435\u0442 (\u0447\u0435\u043c \u0431\u043e\u043b\u044c\u0448\u0435 \u043c\u0430\u0441\u0441\u0438\u0432, \u0442\u0435\u043c \u043c\u0435\u043d\u044c\u0448\u0435 \u0432\u043a\u043b\u0430\u0434 \u043e\u0442 \u043f\u0440\u043e\u0431\u0435\u0436\u0435\u043a \u043f\u043e \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0430\u043c): <\/p>\n
<\/figure>\n
\u0412\u043e\u0442 \u0441\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e \u0438 \u0432\u0441\u0451. \u041a\u043e\u0434 (\u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043b\u0430\u0441\u044c \u0421\u0442\u0443\u0434\u0438\u044f 2022 \u0432\u0435\u0440\u0441\u0438\u0438 17.7.6) \u043b\u0435\u0436\u0438\u0442 \u043d\u0430 \u0433\u0438\u0442\u0445\u0430\u0431\u0435<\/a>, \u0435\u0441\u043b\u0438 \u043a\u043e\u043c\u0443 \u043d\u0443\u0436\u043d\u043e.<\/p>\n
\u0415\u0441\u043b\u0438 \u044f \u0433\u0434\u0435-\u043b\u0438\u0431\u043e \u043e\u0448\u0438\u0431\u0441\u044f, \u043b\u0438\u0431\u043e \u0442\u0430\u043a\u0438 \u0438\u0437\u043e\u0431\u0440\u0451\u043b \u0432\u0435\u043b\u043e\u0441\u0438\u043f\u0435\u0434 \u2014 \u043d\u0435 \u0441\u043e\u0447\u0442\u0438\u0442\u0435 \u0437\u0430 \u0442\u0440\u0443\u0434, \u043d\u0430\u043f\u0438\u0448\u0438\u0442\u0435, \u043f\u043e\u0436\u0430\u043b\u0443\u0439\u0441\u0442\u0430 \u0432 \u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0430\u0445.<\/p>\n
\u0412\u0441\u0435\u043c \u0434\u043e\u0431\u0440\u0430!<\/p>\n<\/div>\n<\/div>\n<\/div>\n
<\/div>\n

 \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438  https:\/\/habr.com\/ru\/articles\/771010\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"
<\/div>\n
\n
\n
\n
\n
\u0423\u0432\u0430\u0436\u0430\u0435\u043c\u044b\u0435 \u043a\u043e\u043b\u043b\u0435\u0433\u0438! \u0412\u0430\u0448\u0435\u043c\u0443 \u0432\u043d\u0438\u043c\u0430\u043d\u0438\u044e \u043f\u0440\u0435\u0434\u043b\u0430\u0433\u0430\u0435\u0442\u0441\u044f \u043d\u0435\u0431\u043e\u043b\u044c\u0448\u043e\u0439 «\u044d\u0442\u044e\u0434 \u0432\u044b\u0445\u043e\u0434\u043d\u043e\u0433\u043e \u0434\u043d\u044f», \u043f\u043e\u0441\u0432\u044f\u0449\u0451\u043d\u043d\u044b\u0439 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e, \u0441\u043a\u0430\u0436\u0435\u043c \u0442\u0430\u043a, \u043d\u0435\u0442\u0440\u0430\u0434\u0438\u0446\u0438\u043e\u043d\u043d\u043e\u043c\u0443 \u0441\u043f\u043e\u0441\u043e\u0431\u0443 \u0432\u044b\u0447\u0438\u0441\u043b\u0435\u043d\u0438\u044f \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u043c\u0430\u0441\u0441\u0438\u0432\u0430 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441 \u043f\u043b\u0430\u0432\u0430\u044e\u0449\u0435\u0439 \u0442\u043e\u0447\u043a\u043e\u0439. \u0412\u043a\u0440\u0430\u0442\u0446\u0435 \u2014 \u043c\u044b \u0441\u0434\u0435\u043b\u0430\u0435\u043c \u044d\u0442\u043e \u0432 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043f\u0440\u043e\u0445\u043e\u0434\u043e\u0432 \u043f\u043e \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c\u0443 \u043c\u0430\u0441\u0441\u0438\u0432\u0443 (\u0434\u0432\u0430 \u0434\u043b\u044f \u043e\u0434\u0438\u043d\u0430\u0440\u043d\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u0438 \u0438\u043b\u0438 \u0447\u0435\u0442\u044b\u0440\u0435 \u0434\u043b\u044f \u0434\u0432\u043e\u0439\u043d\u043e\u0439), \u0432\u044b\u0447\u0438\u0441\u043b\u044f\u044f \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u043f\u043e \u0441\u043b\u043e\u0432\u0430\u043c, \u043d\u0430\u0447\u0438\u043d\u0430\u044f \u0441 \u0431\u043e\u043b\u0435\u0435 \u0437\u043d\u0430\u0447\u0430\u0449\u0438\u0445, \u043f\u043e\u043b\u044c\u0437\u0443\u044f\u0441\u044c \u043f\u0440\u0438 \u044d\u0442\u043e\u043c \u0442\u043e\u043b\u044c\u043a\u043e \u0446\u0435\u043b\u043e\u0447\u0438\u0441\u043b\u0435\u043d\u043d\u043e\u0439 \u0430\u0440\u0438\u0444\u043c\u0435\u0442\u0438\u043a\u043e\u0439, \u0447\u0442\u043e \u0434\u0430\u0441\u0442 \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u043e\u0441\u0442\u044c \u0432 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u043b\u0443\u0447\u0430\u044f\u0445 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043e\u0431\u043e\u0433\u043d\u0430\u0442\u044c \u043f\u043e \u0441\u043a\u043e\u0440\u043e\u0441\u0442\u0438 «\u0442\u0440\u0430\u0434\u0438\u0446\u0438\u043e\u043d\u043d\u044b\u0435» \u043a\u043b\u0430\u0441\u0441\u0438\u0447\u0435\u0441\u043a\u0438\u0435 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u044b. \u0412\u043e\u0437\u043c\u043e\u0436\u043d\u043e \u0434\u0430\u043d\u043d\u0430\u044f «\u0437\u0430\u0440\u0438\u0441\u043e\u0432\u043a\u0430» \u043a\u0430\u043a \u0438\u0434\u0435\u044f \u043e\u043a\u0430\u0436\u0435\u0442\u0441\u044f \u043a\u043e\u043c\u0443-\u043d\u0438\u0431\u0443\u0434\u044c \u043f\u043e\u043b\u0435\u0437\u043d\u0430.<\/p>\n
\u041f\u0440\u0435\u0436\u0434\u0435 \u0447\u0435\u043c \u043f\u0440\u043e\u0434\u043e\u043b\u0436\u0438\u0442\u044c \u2014 \u043c\u0430\u043b\u0435\u043d\u044c\u043a\u0438\u0439 \u0434\u0438\u0441\u043a\u043b\u0435\u0439\u043c\u0435\u0440: \u044f \u043d\u0438\u043a\u0430\u043a \u043d\u0435 \u0441\u0432\u044f\u0437\u0430\u043d \u043d\u0438 \u0441 \u043a\u0430\u043a\u043e\u0439 \u0438\u0437 \u0443\u043f\u043e\u043c\u0438\u043d\u0430\u044e\u0449\u0438\u0445\u0441\u044f \u043d\u0438\u0436\u0435 \u043a\u043e\u043c\u043f\u0430\u043d\u0438\u0439, \u043f\u043e\u0441\u0442 \u043d\u0435 \u043d\u043e\u0441\u0438\u0442 \u0440\u0435\u043a\u043b\u0430\u043c\u043d\u043e\u0433\u043e \u0445\u0430\u0440\u0430\u043a\u0442\u0435\u0440\u0430, \u043a\u043e\u0434 \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0439\u0442\u0435 \u043d\u0430 \u0441\u0432\u043e\u0439 \u0441\u0442\u0440\u0430\u0445 \u0438 \u0440\u0438\u0441\u043a.<\/p>\n
\u0410 \u043d\u0430\u0447\u0430\u043b\u0430\u0441\u044c \u044d\u0442\u0430 \u0438\u0441\u0442\u043e\u0440\u0438\u044f \u0441 \u0440\u044f\u0434\u043e\u0432\u043e\u0439 \u0432 \u043e\u0431\u0449\u0435\u043c-\u0442\u043e \u0437\u0430\u0434\u0430\u0447\u043a\u0438 \u2014 \u043c\u043d\u0435 \u043d\u0443\u0436\u043d\u043e \u0431\u044b\u043b\u043e \u043d\u0430\u0439\u0442\u0438 \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u0440\u0435\u043d\u0442\u0433\u0435\u043d\u043e\u0432\u0441\u043a\u043e\u0439 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438 (\u044d\u0442\u043e \u043e\u0431\u044b\u0447\u043d\u043e 12-14 \u0431\u0438\u0442 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435) \u0432 \u0440\u0430\u043c\u043a\u0430\u0445 \u043f\u0440\u043e\u0432\u0435\u0440\u043a\u0438 \u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 Flat Field Correction<\/a>. \u041c\u043e\u044f \u043e\u0441\u043d\u043e\u0432\u043d\u0430\u044f \u0440\u0430\u0431\u043e\u0447\u0430\u044f \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0430 — NI Vision Development Module, \u0438 \u0442\u0430\u043c \u0435\u0441\u0442\u044c \u0438 \u0441\u0440\u0435\u0434\u043d\u0435\u0435 \u0438 \u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u043e\u0435 \u043e\u0442\u043a\u043b\u043e\u043d\u0435\u043d\u0438\u0435, \u0438 \u043c\u043d\u043e\u0433\u043e \u0447\u0435\u0433\u043e \u0434\u0440\u0443\u0433\u043e\u0433\u043e, \u0430 \u0432\u043e\u0442 \u043f\u0440\u0438\u043c\u0438\u0442\u0438\u0432\u043d\u043e\u0439 \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u043d\u0435 \u043e\u043a\u0430\u0437\u0430\u043b\u043e\u0441\u044c. \u041d\u0430 «\u043f\u043e\u0434\u0445\u0432\u0430\u0442\u0435» \u044f \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e \u0442\u0430\u043a\u0436\u0435 \u0431\u0438\u0431\u043b\u0438\u043e\u0442\u0435\u043a\u0438 OpenCV \u043b\u0438\u0431\u043e Intel IPP, \u043d\u043e \u0438 \u0442\u0430\u043c \u043d\u0430\u0432\u0441\u043a\u0438\u0434\u043a\u0443 \u0442\u0430\u043a\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 \u0432 \u0433\u043e\u0442\u043e\u0432\u043e\u043c \u0432\u0438\u0434\u0435 \u044f \u043d\u0435 \u043e\u0431\u043d\u0430\u0440\u0443\u0436\u0438\u043b (\u0430 \u043c\u043e\u0436\u0435\u0442 \u043f\u043b\u043e\u0445\u043e \u0438\u0441\u043a\u0430\u043b). \u041d\u0443 \u043d\u0435 \u0431\u0435\u0434\u0430 \u2014 \u0440\u0430\u0441\u0447\u0435\u0445\u043b\u0438\u043c \u043a\u043e\u043c\u043f\u0438\u043b\u044f\u0442\u043e\u0440 \u0434\u0430 \u043d\u0430\u043f\u0438\u0448\u0435\u043c \u0441\u0430\u043c\u0438. \u0414\u0430, \u0435\u0441\u043b\u0438 \u0432\u044b \u0437\u043d\u0430\u0435\u0442\u0435, \u043a\u0430\u043a \u0440\u0430\u0431\u043e\u0442\u0430\u0435\u0442 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u0439 \u043c\u0435\u0442\u043e\u0434, \u0442\u043e \u043c\u043e\u0436\u0435\u0442\u0435 \u0441\u0440\u0430\u0437\u0443 \u043f\u0440\u043e\u043c\u0430\u0442\u044b\u0432\u0430\u0442\u044c \u043f\u0440\u0438\u043c\u0435\u0440\u043d\u043e \u0434\u043e \u0441\u0435\u0440\u0435\u0434\u0438\u043d\u044b, \u043b\u0438\u0431\u043e \u043f\u0440\u044f\u043c\u043e \u0434\u043e \u043a\u043e\u043d\u0446\u0430, \u0442\u0430\u043c \u0433\u0434\u0435 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u044b \u0437\u0430\u043c\u0435\u0440\u0430 \u043f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438.<\/p>\n
\u041f\u0440\u0438 \u0440\u0430\u0431\u043e\u0442\u0435 \u043f\u043e \u0448\u0435\u0441\u0442\u043d\u0430\u0434\u0446\u0430\u0442\u0438\u0431\u0438\u0442\u043d\u043e\u043c\u0443 \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u043c\u0435\u0434\u0438\u0430\u043d\u0443 \u043e\u0431\u044b\u0447\u043d\u043e \u043f\u0440\u0438\u043d\u044f\u0442\u043e \u0432\u044b\u0447\u0438\u0441\u043b\u044f\u0442\u044c, \u043f\u043e\u043b\u044c\u0437\u0443\u044f\u0441\u044c \u043d\u0435\u0441\u043b\u043e\u0436\u043d\u044b\u043c \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u043d\u044b\u043c \u043c\u0435\u0442\u043e\u0434\u043e\u043c. \u041c\u0435\u0434\u0438\u0430\u043d\u0430 \u2014 \u044d\u0442\u043e \u0432\u0435\u0434\u044c \u0442\u0430\u043a\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435, \u043a\u043e\u0442\u043e\u0440\u0430\u044f \u0434\u0435\u043b\u0438\u0442 \u043c\u0430\u0441\u0441\u0438\u0432 (\u0435\u0441\u043b\u0438 \u0435\u0433\u043e \u043e\u0442\u0441\u043e\u0440\u0442\u0438\u0440\u043e\u0432\u0430\u0442\u044c, \u0440\u0430\u0437\u0443\u043c\u0435\u0435\u0442\u0441\u044f) \u043f\u043e\u043f\u043e\u043b\u0430\u043c, \u0442\u0430\u043a \u0447\u0442\u043e \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u0430 \u0432\u0441\u0435\u0445 \u043e\u0441\u0442\u0430\u043b\u044c\u043d\u044b\u0445 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043f\u043e \u0432\u0435\u043b\u0438\u0447\u0438\u043d\u0435 \u0431\u0443\u0434\u0435\u0442 \u0441\u043b\u0435\u0432\u0430 \u043e\u0442 \u043c\u0435\u0434\u0438\u0430\u043d\u044b, \u0430 \u0434\u0440\u0443\u0433\u0430\u044f \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u0430 — \u0441\u043f\u0440\u0430\u0432\u0430. \u041f\u043e\u044d\u0442\u043e\u043c\u0443 \u0437\u0430 \u043e\u0434\u0438\u043d \u043f\u0440\u043e\u0445\u043e\u0434 \u043f\u043e \u0432\u0441\u0435\u043c\u0443 \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u0441\u0442\u0440\u043e\u044f\u0442 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443, \u0430 \u0434\u0430\u043b\u044c\u0448\u0435 \u043f\u0435\u0440\u0435\u043c\u0435\u0449\u0430\u044e\u0442\u0441\u044f \u043f\u043e \u043d\u0435\u0439, \u043e\u0441\u0442\u0430\u043d\u0430\u0432\u043b\u0438\u0432\u0430\u044f\u0441\u044c \u0432 \u0442\u043e\u0442 \u043c\u043e\u043c\u0435\u043d\u0442, \u043a\u043e\u0433\u0434\u0430 \u0447\u0438\u0441\u043b\u043e \u043e\u0442\u0441\u0447\u0451\u0442\u043e\u0432 \u0434\u043e\u0441\u0442\u0438\u0433\u043d\u0435\u0442 \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u044b \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u0433\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0430 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432: <\/p>\n
<\/figure>\n
\u041f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u044f \u0443\u0436\u0435 \u043e\u0442\u043d\u043e\u0441\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u043d\u0435\u043c\u043e\u043b\u043e\u0434 \u0438 \u0432 \u043c\u0435\u0440\u0443 \u043b\u0435\u043d\u0438\u0432, \u0442\u043e \u0434\u043b\u044f \u0437\u0430\u0442\u0440\u0430\u0432\u043a\u0438 \u044f \u0432\u043e\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e\u0441\u044c ChatGPT. \u041d\u0443 \u0434\u0430\u0432\u0430\u0439, \u0433\u043e\u0440\u0448\u043e\u0447\u0435\u043a, \u0432\u0430\u0440\u0438! \u041f\u043e\u0434\u043d\u0430\u0442\u0443\u0436\u0438\u0432\u0448\u0438\u0441\u044c, \u043e\u043d\u043e \u043c\u043d\u0435 \u0432\u044b\u0434\u0430\u043b\u043e \u0432\u043e\u0442 \u0442\u0430\u043a\u043e\u0439 \u0448\u0435\u0434\u0435\u0432\u0440, \u043f\u0440\u0438\u0432\u0435\u0434\u0443 \u0435\u0433\u043e \u043a\u0430\u043a \u0435\u0441\u0442\u044c:<\/p>\n
unsigned short computeMedian(unsigned short* arr, int size) {     \/\/ Create histogram array     int maxVal = 65536; \/\/ Maximum value for unsigned short     int* histogram = (int*)calloc(maxVal, sizeof(int));      \/\/ Populate histogram     for (int i = 0; i < size; i++) {         histogram[arr[i]]++;     }      \/\/ Find the median value     int medianPos = (size + 1) \/ 2;     unsigned short median = 0;     int count = 0;      for (int i = 0; i < maxVal; i++) {         count += histogram[i];         if (count >= medianPos) {             median = i;             break;         }     }      \/\/ Free memory     free(histogram);      return median; } <\/code><\/pre>\n
\u041d\u0430 \u0441\u0430\u043c\u043e\u043c \u0434\u0435\u043b\u0435 \u043d\u0435\u043f\u043b\u043e\u0445\u043e, \u043e\u043d\u043e \u0434\u0430\u0436\u0435 \u0437\u043d\u0430\u0435\u0442, \u0447\u0442\u043e calloc \u0438\u043d\u0438\u0446\u0438\u0430\u043b\u0438\u0437\u0438\u0440\u0443\u0435\u0442 \u043f\u0430\u043c\u044f\u0442\u044c, \u043d\u043e \u0435\u0441\u0442\u044c \u043d\u044e\u0430\u043d\u0441 \u2014 \u043c\u043d\u0435 \u0434\u0430\u0436\u0435 \u0437\u0430\u043f\u0443\u0441\u043a\u0430\u0442\u044c \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u043d\u0435 \u043d\u0443\u0436\u043d\u043e, \u0447\u0442\u043e\u0431\u044b \u043f\u043e\u043d\u044f\u0442\u044c, \u0447\u0442\u043e \u044d\u0442\u043e \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c, \u0442\u043e\u043b\u044c\u043a\u043e \u0435\u0441\u043b\u0438 \u0432 \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u043d\u0435\u0447\u0451\u0442\u043d\u043e\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432. \u0410 \u0432 \u0447\u0451\u0442\u043d\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435 \u043c\u044b \u0431\u0435\u0440\u0451\u043c \u0441\u0440\u0435\u0434\u043d\u0435\u0435 \u0430\u0440\u0438\u0444\u043c\u0435\u0442\u0438\u0447\u0435\u0441\u043a\u043e\u0435 \u043e\u0442 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439 \u0441\u043b\u0435\u0432\u0430 \u0438 \u0441\u043f\u0440\u0430\u0432\u0430 \u043e\u0442 \u043c\u0435\u0434\u0438\u0430\u043d\u044b \u0438 \u043c\u043e\u0439 \u0432\u043d\u0443\u0442\u0440\u0435\u043d\u043d\u0438\u0439 \u043f\u0435\u0440\u0444\u0435\u043a\u0446\u0438\u043e\u043d\u0438\u0441\u0442 \u043f\u0440\u043e\u0441\u0438\u0442, \u0447\u0442\u043e\u0431\u044b \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0442\u043e\u0447\u043d\u043e \u0441\u043e\u0432\u043f\u0430\u0434\u0430\u043b\u043e \u0441 \u0440\u0435\u0444\u0435\u0440\u0435\u043d\u0441\u043d\u044b\u043c.  \u041f\u043e\u043f\u0440\u043e\u0441\u0438\u043c \u043f\u043e\u043f\u0440\u0430\u0432\u0438\u0442\u044c:<\/p>\n
double computeMedian(unsigned short* arr, int size) {     \/\/ Create histogram array     int maxVal = 65536; \/\/ Maximum value for unsigned short     int* histogram = (int*)calloc(maxVal, sizeof(int));      \/\/ Populate histogram     for (int i = 0; i < size; i++) {         histogram[arr[i]]++;     }      \/\/ Find the median position     int medianPos = (size + 1) \/ 2;     int count = 0;     int medianVal = 0;      \/\/ Iterate over histogram to find the median value     for (int i = 0; i < maxVal; i++) {         count += histogram[i];         if (count >= medianPos) {             medianVal = i;             break;         }     }      \/\/ If the size is odd, return the median value directly     if (size % 2 != 0) {         return (double)medianVal;     }      \/\/ If the size is even, find the next value to compute the median     int nextVal = medianVal + 1;     while (histogram[nextVal] == 0) {         nextVal++;     }      \/\/ Compute the median value for even-sized array     return (medianVal + nextVal) \/ 2.0; } <\/code><\/pre>\n
\u0423\u0436\u0435 \u043b\u0443\u0447\u0448\u0435, \u043d\u043e \u0432\u043e-\u043f\u0435\u0440\u0432\u044b\u0445, \u043e\u043d\u043e \u043f\u043e\u0447\u0435\u043c\u0443-\u0442\u043e \u043f\u043e\u0442\u0435\u0440\u044f\u043b\u043e \u043e\u0441\u0432\u043e\u0431\u043e\u0436\u0434\u0435\u043d\u0438\u0435 \u0432\u044b\u0434\u0435\u043b\u0435\u043d\u043d\u043e\u0439 \u043f\u043e\u0434 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u043f\u0430\u043c\u044f\u0442\u0438, \u0430 \u0432\u043e-\u0432\u0442\u043e\u0440\u044b\u0445, \u044d\u0442\u043e\u0442 \u043a\u043e\u0434 \u0431\u0443\u0434\u0435\u0442 \u043f\u0440\u043e\u0444\u0443\u043a\u0438\u0432\u0430\u0442\u044c \u0438\u0441\u0442\u0438\u043d\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u0435\u0434\u0438\u0430\u043d\u044b, \u0435\u0441\u043b\u0438 \u0432 \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u0431\u0443\u0434\u0443\u0442 \u043f\u043e\u0432\u0442\u043e\u0440\u044f\u044e\u0449\u0438\u0435\u0441\u044f \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u044b. \u0414\u0430\u043b\u044c\u043d\u0435\u0439\u0448\u0438\u0435 \u043f\u0440\u043e\u0441\u044c\u0431\u044b \u0438\u0441\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u0438\u0439 \u0431\u0443\u0434\u0443\u0442 \u0440\u0430\u0437\u0434\u0443\u0432\u0430\u0442\u044c \u043a\u043e\u0434 \u0435\u0449\u0451 \u0441\u0438\u043b\u044c\u043d\u0435\u0435 (\u0430 \u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u0438\u0441\u0442\u043e\u0432 \u0441\u043a\u043e\u043f\u0438\u043f\u0430\u0441\u0442\u044f\u0442 \u0435\u0433\u043e \u043f\u0440\u044f\u043c\u043e \u0432 \u0442\u0430\u043a\u043e\u043c \u0432\u0438\u0434\u0435?). \u041f\u043e\u0436\u0430\u043b\u0443\u0439 \u0431\u0435\u0437 \u0440\u0430\u0431\u043e\u0442\u044b ChatGPT \u043d\u0430\u0441 \u043f\u043e\u043a\u0430 \u0447\u0442\u043e \u043d\u0435 \u043e\u0441\u0442\u0430\u0432\u0438\u0442. \u0417\u0430\u0441\u0443\u0447\u0438\u0432\u0430\u0435\u043c \u0440\u0443\u043a\u0430\u0432\u0430 \u0438 \u043f\u0440\u0430\u0432\u0438\u043c \u043a\u043e\u0434 \u0440\u0443\u0447\u043a\u0430\u043c\u0438, \u043d\u0443 \u043a\u0430\u043a-\u0442\u043e \u0432\u043e\u0442 \u0442\u0430\u043a \u0447\u0442\u043e \u043b\u0438:<\/p>\n
double computeMedian(unsigned short* arr, size_t size) { size_t* histogram; double Median; if (!arr || !size) return -1; if (!( histogram = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t)))) return -2;  for (int i = 0; i < size; i++) histogram[arr[i]]++;    \/\/ Populate histogram size_t medianPos = (size + 1) \/ 2, count = 0, medianVal; \/\/ Find the median pos for (medianVal = 0; medianVal <= USHRT_MAX; medianVal++) { \/\/ Iterate over histo count += histogram[medianVal]; if (count >= medianPos) break; }  if (!(size % 2 || count > medianPos)) { \/\/ If the size is even, find the next size_t nextVal = medianVal; while (!histogram[++nextVal]); Median = (medianVal + nextVal) \/ 2.0; \/\/ Middle } else Median = (double)medianVal;  free(histogram); return Median; } <\/code><\/pre>\n
\u041d\u0430 \u043a\u043e\u0434 \u0440\u0435\u0432\u044c\u044e, \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e, \u043c\u043e\u0433\u0443\u0442 \u0432\u043e\u0437\u043d\u0438\u043a\u043d\u0443\u0442\u044c \u043a\u043e\u0435-\u043a\u0430\u043a\u0438\u0435 \u0432\u043e\u043f\u0440\u043e\u0441\u044b, \u043d\u043e \u0432\u0441\u0435\u0433\u0434\u0430 \u043c\u043e\u0436\u043d\u043e \u0441\u0432\u0430\u043b\u0438\u0442\u044c \u043d\u0430 \u0438\u0441\u043a\u0443\u0441\u0441\u0442\u0432\u0435\u043d\u043d\u044b\u0439 \u0438\u043d\u0442\u0435\u043b\u043b\u0435\u043a\u0442, \u043e\u043f\u044f\u0442\u044c \u0436\u0435 \u0441\u0442\u0430\u0442\u0438\u0447\u0435\u0441\u043a\u0438\u0439 \u0430\u043d\u0430\u043b\u0438\u0437\u0430\u0442\u043e\u0440 \u043c\u043e\u043b\u0447\u0438\u0442 \u043a\u0430\u043a \u043f\u0430\u0440\u0442\u0438\u0437\u0430\u043d.<\/p>\n
32 \u0431\u0438\u0442\u0430<\/h3>\n
\u0420\u044f\u0434\u043e\u0432\u0430\u044f \u0432 \u043e\u0431\u0449\u0435\u043c-\u0442\u043e \u0437\u0430\u0434\u0430\u0447\u043a\u0430 \u043d\u0435 \u0441\u0442\u043e\u0438\u043b\u0430 \u0431\u044b \u043f\u043e\u0441\u0442\u0430 \u043d\u0430 \u0425\u0430\u0431\u0440\u0435,  \u0434\u0430 \u0438 \u0435\u0433\u043e \u0431\u044b \u0438 \u043d\u0435 \u043f\u043e\u044f\u0432\u0438\u043b\u043e\u0441\u044c, \u0435\u0441\u043b\u0438 \u0431\u044b \u0432\u043e \u0432\u0440\u0435\u043c\u044f \u0435\u0436\u0435\u0432\u0435\u0447\u0435\u0440\u043d\u0435\u0439 \u043f\u0440\u043e\u0433\u0443\u043b\u043a\u0438 \u0441 \u0441\u043e\u0431\u0430\u043a\u0435\u043d\u044b\u043c \u043c\u043e\u0439 \u0432\u043d\u0443\u0442\u0440\u0435\u043d\u043d\u0438\u0439 \u0433\u043e\u043b\u043e\u0441 \u0435\u0445\u0438\u0434\u043d\u043e \u0438 \u0432\u043a\u0440\u0430\u0434\u0447\u0438\u0432\u043e \u043d\u0435 \u043f\u043e\u0438\u043d\u0442\u0435\u0440\u0435\u0441\u043e\u0432\u0430\u043b\u0441\u044f: «\u2014 \u0430 \u0432\u043e\u0442 \u0435\u0441\u043b\u0438, \u043f\u0440\u0435\u0434\u043f\u043e\u043b\u043e\u0436\u0438\u043c, \u0443 \u0442\u0435\u0431\u044f \u043c\u0430\u0441\u0441\u0438\u0432 \u043d\u0435 16-\u0438, \u0430 32-\u0445 \u0431\u0438\u0442\u043d\u044b\u0439, \u0442\u043e \u043a\u0430\u043a \u0442\u043e\u0433\u0434\u0430?»<\/em>. \u0412\u0435\u0434\u044c \u043d\u0430 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0432 \u0447\u0435\u0442\u044b\u0440\u0435 \u0441 \u043b\u0438\u0448\u043d\u0438\u043c \u043c\u0438\u043b\u043b\u0438\u0430\u0440\u0434\u043e\u0432 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043d\u0438\u043a\u0430\u043a\u043e\u0439 \u043f\u0430\u043c\u044f\u0442\u0438 \u043d\u0435 \u0445\u0432\u0430\u0442\u0438\u0442! \u041f\u043e\u0440\u0430\u0437\u043c\u044b\u0448\u043b\u044f\u0432 \u043d\u0435\u043c\u043d\u043e\u0433\u043e, \u044f \u0435\u0449\u0451 \u0434\u043e \u043a\u043e\u043d\u0446\u0430 \u043f\u0440\u043e\u0433\u0443\u043b\u043a\u0438 \u0438\u0437\u043c\u044b\u0441\u043b\u0438\u043b \u043f\u0440\u043e\u0441\u0442\u043e\u0439 \u043d\u0435\u0437\u0430\u0442\u0435\u0439\u043b\u0438\u0432\u044b\u0439 \u043c\u0435\u0442\u043e\u0434 \u2014 \u0432\u043d\u0430\u0447\u0430\u043b\u0435 \u043d\u0430\u0434\u043e \u043f\u0440\u043e\u0439\u0442\u0438 \u043f\u043e \u0441\u0442\u0430\u0440\u0448\u0438\u043c \u0441\u043b\u043e\u0432\u0430\u043c, \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u0432 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0432 65 \u043a\u0438\u043b\u043e\u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043a\u0430\u043a \u0438 \u0432\u044b\u0448\u0435, \u0430 \u0437\u0430\u0442\u0435\u043c, \u043e\u0442\u0442\u0430\u043b\u043a\u0438\u0432\u0430\u044f\u0441\u044c \u043e\u0442 \u043d\u0430\u0439\u0434\u0435\u043d\u043d\u043e\u0433\u043e \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f \u043f\u0440\u043e\u0439\u0442\u0438 \u043f\u043e \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u0435\u0449\u0451 \u0440\u0430\u0437, \u0441\u0442\u0440\u043e\u044f \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0443\u0436\u0435 \u043f\u043e \u043c\u043b\u0430\u0434\u0448\u0435\u043c\u0443 \u0441\u043b\u043e\u0432\u0443 \u0438 \u0443\u0442\u043e\u0447\u043d\u044f\u044f \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u0435\u0434\u0438\u0430\u043d\u044b. \u0412 \u043f\u0440\u0438\u043d\u0446\u0438\u043f\u0435 \u043c\u0435\u0442\u043e\u0434 \u043d\u0435 \u043d\u043e\u0432 \u2014 \u043c\u044b \u0444\u0430\u043a\u0442\u0438\u0447\u0435\u0441\u043a\u0438 \u0440\u0430\u0437\u0431\u0438\u0432\u0430\u0435\u043c \u0432\u0435\u0441\u044c \u0434\u0438\u0430\u043f\u0430\u0437\u043e\u043d \u043d\u0430 65536 \u043e\u0442\u0440\u0435\u0437\u043a\u043e\u0432, \u0438 \u043d\u0430\u0439\u0434\u044f \u043f\u043e\u0434\u0445\u043e\u0434\u044f\u0449\u0438\u0439 \u043e\u0442\u0440\u0435\u0437\u043e\u043a, \u0438\u0449\u0435\u043c \u0442\u0430\u043c \u043e\u043a\u043e\u043d\u0447\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435. \u041c\u0430\u043b\u0435\u043d\u044c\u043a\u0438\u0439 \u043d\u044e\u0430\u043d\u0441 \u0435\u0441\u0442\u044c \u0441 \u0447\u0430\u0441\u0442\u043d\u044b\u043c \u0441\u043b\u0443\u0447\u0430\u0435\u043c \u043a\u043e\u0433\u0434\u0430 \u0432 \u0438\u0441\u0445\u043e\u0434\u043d\u043e\u043c \u043c\u0430\u0441\u0441\u0438\u0432\u0435 \u0447\u0451\u0442\u043d\u043e\u0435 \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0438 \u043d\u0430\u0448\u0430 \u043c\u0435\u0434\u0438\u0430\u043d\u0430 \u043b\u0435\u0436\u0438\u0442 \u043c\u0435\u0436\u0434\u0443 \u0440\u0430\u0437\u043d\u044b\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438 \u0441\u0442\u0430\u0440\u0448\u0435\u0433\u043e \u0441\u043b\u043e\u0432\u0430. \u041d\u043e \u0435\u0441\u043b\u0438 \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u0438\u0442\u044c \u0432 \u0443\u043c\u0435 \u043e\u0442\u0441\u043e\u0440\u0442\u0438\u0440\u043e\u0432\u0430\u043d\u043d\u044b\u0439 \u043c\u0430\u0441\u0441\u0438\u0432, \u0442\u043e \u0441\u0442\u0430\u043d\u0435\u0442 \u044f\u0441\u043d\u043e, \u0447\u0442\u043e \u043e\u0442\u0434\u0435\u043b\u044c\u043d\u044b\u0435 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u044b \u0442\u0443\u0442 \u0443\u0436\u0435 \u043d\u0435 \u043d\u0443\u0436\u043d\u044b, \u043f\u0440\u043e\u0441\u0442\u043e \u0441\u043b\u0435\u0432\u0430 (\u0434\u043b\u044f \u043c\u0435\u043d\u044c\u0448\u0435\u0433\u043e \u0441\u0442\u0430\u0440\u0448\u0435\u0433\u043e \u0431\u0430\u0439\u0442\u0430) \u0441\u043b\u0435\u0434\u0443\u0435\u0442 \u0432\u0437\u044f\u0442\u044c \u043c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043c\u043b\u0430\u0434\u0448\u0435\u0433\u043e \u0431\u0430\u0439\u0442\u0430, \u0430 \u0441\u043f\u0440\u0430\u0432\u0430 (\u0434\u043b\u044f \u0431\u043e\u043b\u044c\u0448\u0435\u0433\u043e) \u2014 \u043c\u0438\u043d\u0438\u043c\u0430\u043b\u044c\u043d\u043e\u0435, \u0438 \u043f\u043e\u0442\u043e\u043c \u0443\u0436\u0435 \u043d\u0430\u0439\u0442\u0438 \u0441\u0440\u0435\u0434\u043d\u0435\u0435, \u0432\u043e\u0442 \u0438 \u0432\u0441\u0451: <\/p>\n
<\/figure>\n
 \u042f \u043d\u0430\u0431\u0440\u043e\u0441\u0430\u043b \u0431\u044b\u043b\u043e \u043a\u043e\u0434 \u043d\u0430 LabVIEW, \u043d\u043e \u043f\u043e\u0436\u0430\u043b\u0443\u0439 \u043d\u0435 \u0441\u0442\u0430\u043d\u0443 \u0432\u0430\u0441 \u043c\u0443\u0447\u0438\u0442\u044c, \u0445\u043e\u0442\u044c \u043c\u044b \u0438 \u0432 \u0445\u0430\u0431\u0435 «\u041d\u0435\u043d\u043e\u0440\u043c\u0430\u043b\u044c\u043d\u043e\u0435 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u0435». <\/p>\n
\u041d\u0430 \u0421\u0438 \u043e\u043d\u043e \u0431\u0443\u0434\u0435\u0442 \u0442\u0430\u043a \u0436\u0435, \u043a\u0430\u043a \u0438 \u0432 \u0441\u043b\u0443\u0447\u0430\u0435 16 \u0431\u0438\u0442, \u0432\u043e\u0442 \u0442\u043e\u043b\u044c\u043a\u043e \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u043c\u044b \u0441\u0442\u0440\u043e\u0438\u043c \u043f\u043e \u0441\u0442\u0430\u0440\u0448\u0438\u043c \u0431\u0438\u0442\u0430\u043c, \u0437\u0430\u0434\u0432\u0438\u0433\u0430\u044f \u0438\u0445 \u043d\u0430\u043f\u0440\u0430\u0432\u043e:<\/p>\n
LIBMED_API int MedianU32(unsigned int* ptr, size_t Length, double* med) {     \/\/...     histo = (size_t*)calloc(USHRT_MAX + 1, sizeof(size_t));     for (int i = 0; i < Length; i++) histo[(ptr[i] & 0xFFFF0000) >> 16]++;     \/\/... <\/code><\/pre>\n
\u0414\u043e \u044d\u0442\u043e\u0433\u043e \u043c\u043e\u043c\u0435\u043d\u0442\u0430 \u0432\u0441\u0451 \u0440\u043e\u0432\u043d\u043e \u043a\u0430\u043a \u0438 \u0440\u0430\u043d\u044c\u0448\u0435:<\/p>\n
uint64_t medianPos = (Length + 1) \/ 2, count = 0;  for (int i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram to find the med count += histo[i]; if (count >= medianPos) { medianValMSW = i; break; } } if (!(Length % 2 || count > medianPos)) { \/\/ If the size is even, find the next <\/code><\/pre>\n
\u0415\u0441\u043b\u0438 \u0443\u0441\u043b\u043e\u0432\u0438\u0435 \u0441\u0440\u0430\u0431\u043e\u0442\u0430\u043b\u043e, \u0442\u043e \u043d\u0430 \u0432\u0442\u043e\u0440\u043e\u043c \u043f\u0440\u043e\u0445\u043e\u0434\u0435 \u043f\u0440\u043e\u0441\u0442\u043e \u0441\u043c\u043e\u0442\u0440\u0438\u043c \u043c\u0438\u043d\u0438\u043c\u0443\u043c \u0438 \u043c\u0430\u043a\u0441\u0438\u043c\u0443\u043c:<\/p>\n
if (!(Length % 2 || count > medianPos)) { \/\/ If the size is even, find the next nextVal = medianValMSW; while (!histo[++nextVal]);  for (max = 0, min = 65535, i = 0; i < Length; i++) { indexMSW = (ptr[i] & 0xFFFF0000) >> 16; if (medianValMSW == indexMSW && (ptr[i] & 0x0000FFFF) > max) max = (ptr[i] & 0x0000FFFF); if (indexMSW == nextVal && (ptr[i] & 0x0000FFFF) < min) min = (ptr[i] & 0x0000FFFF); }  medianLeft = (medianValMSW << 16) | max; medianRight = (nextVal << 16) | min; *med = ((double)medianLeft + (double)medianRight) \/ 2.0; <\/code><\/pre>\n
\u0410 \u0435\u0441\u043b\u0438 \u0443\u0441\u043b\u043e\u0432\u0438\u0435 \u043d\u0435 \u0441\u0440\u0430\u0431\u043e\u0442\u0430\u043b\u043e, \u0442\u043e \u0441\u0431\u0440\u0430\u0441\u044b\u0432\u0430\u0435\u043c \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0438 \u0434\u0435\u043b\u0430\u0435\u043c \u0432\u0442\u043e\u0440\u043e\u0439 \u043f\u0440\u043e\u0445\u043e\u0434 \u043f\u043e \u043c\u043b\u0430\u0434\u0448\u0435\u043c\u0443 \u0441\u043b\u043e\u0432\u0443 \u0432\u043e\u0442 \u0442\u0430\u043a, \u0437\u0430\u043d\u043e\u0441\u044f \u0432 \u0433\u0438\u0441\u0442\u043e\u0433\u0440\u0430\u043c\u043c\u0443 \u0442\u043e\u043b\u044c\u043a\u043e \u0442\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f \u0443 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u0442\u0430\u0440\u0448\u0435\u0435 \u0441\u043b\u043e\u0432\u043e \u0440\u0430\u0432\u043d\u043e \u043d\u0430\u0439\u0434\u0435\u043d\u043d\u043e\u043c\u0443 \u043d\u0430 \u043f\u0435\u0440\u0432\u043e\u043c \u043f\u0440\u043e\u0445\u043e\u0434\u0435 \u0438 \u043f\u043e\u043f\u0443\u0442\u043d\u043e \u0441\u0447\u0438\u0442\u0430\u044f \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u0441\u043b\u0435\u0432\u0430:<\/p>\n
ZeroMemory(histo, (USHRT_MAX + 1) * sizeof(size_t)); med_ = medianValMSW << 16; for (count = 0, i = 0; i < Length; i++) { indexMSW = (ptr[i] & 0xFFFF0000) >> 16; if (medianValMSW == indexMSW) histo[(ptr[i] & 0x0000FFFF)]++; if (ptr[i] < med_) count++; }  \/\/ Find the median position; medianValLSW = 0; for (int i = 0; i <= USHRT_MAX; i++) { \/\/ Iterate over histogram  count += histo[i]; if (count >= medianPos) { medianValLSW = i; break; } } <\/code><\/pre>\n
\u041d\u0443 \u0438 \u0437\u0434\u0435\u0441\u044c \u043d\u0430\u0441 \u043c\u043e\u0436\u0435\u0442 \u043d\u0430\u0441\u0442\u0438\u0433\u043d\u0443\u0442\u044c \u0442\u0430\u043a\u0430\u044f \u0436\u0435 \u0441\u0438\u0442\u0443\u0430\u0446\u0438\u044f, \u043a\u043e\u0433\u0434\u0430 \u043b\u0438\u0431\u043e \u043c\u0435\u0434\u0438\u0430\u043d\u0430 \u0441\u043e\u0431\u0438\u0440\u0430\u0435\u0442\u0441\u044f \u0438\u0437 \u0434\u0432\u0443\u0445 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0439, \u043b\u0438\u0431\u043e \u0438\u0437 \u043e\u0434\u043d\u043e\u0433\u043e:<\/p>\n
if (!(Length % 2 || count > medianPos)) { nextVal = medianValLSW; while (!histo[++nextVal]);  medianLeft = (medianValMSW << 16) | medianValLSW; medianRight = (medianValMSW << 16) | nextVal; *med = ((double)medianLeft + (double)medianRight) \/ 2.0; } else { medianValI32 = (medianValMSW << 16) | medianValLSW; *med = (double)medianValI32; } <\/code><\/pre>\n
\u0412\u043e\u0442, \u0441\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e \u0438 \u0432\u0441\u0451. \u041f\u043e\u043d\u044f\u0442\u043d\u043e, \u0447\u0442\u043e \u043c\u0435\u0442\u043e\u0434 \u0431\u0443\u0434\u0435\u0442 \u0440\u0430\u0431\u043e\u0442\u0430\u0442\u044c \u0438 \u0434\u043b\u044f 64-\u0445 \u0431\u0438\u0442\u043d\u044b\u0445 \u043c\u0430\u0441\u0441\u0438\u0432\u043e\u0432, \u043f\u0440\u043e\u0441\u0442\u043e \u043f\u0440\u043e\u0439\u0442\u0438 \u043f\u043e \u043c\u0430\u0441\u0441\u0438\u0432\u0443 \u043f\u0440\u0438\u0434\u0451\u0442\u0441\u044f \u0447\u0435\u0442\u044b\u0440\u0435 \u0440\u0430\u0437\u0430, \u043f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u0443 \u043d\u0430\u0441 \u0447\u0435\u0442\u044b\u0440\u0435 \u0441\u043b\u043e\u0432\u0430 \u043f\u043e \u0434\u0432\u0430 \u0431\u0430\u0439\u0442\u0430, \u043f\u0440\u0438\u0447\u0451\u043c \u0447\u0430\u0441\u0442\u043d\u044b\u0439 \u0441\u043b\u0443\u0447\u0430\u0439 \u0441 \u0440\u0430\u0437\u043d\u044b\u043c\u0438 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438 \u0441\u043b\u0435\u0432\u0430 \u0438 \u0441\u043f\u0440\u0430\u0432\u0430 \u043f\u0440\u0438 \u0447\u0451\u0442\u043d\u043e\u043c \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u0435 \u044d\u043b\u0435\u043c\u0435\u043d\u0442\u043e\u0432 \u043c\u043e\u0436\u0435\u0442 \u0432\u043e\u0437\u043d\u0438\u043a\u043d\u0443\u0442\u044c \u043d\u0430 \u043b\u044e\u0431\u043e\u043c \u043f\u0440\u043e\u0433\u043e\u043d\u0435.<\/p>\n
\n\u043f\u043e\u043b\u043d\u044b\u0439 \u043b\u0438\u0441\u0442\u0438\u043d\u0433 MedianU32<\/summary>\n
\n
 LIBMED_API int MedianU32(unsigned int* ptr, size_t Length, double* med) { size_t* histo; int i; uint64_t medianLeft, medianRight, indexMSW, nextVal, med_;<\/code><\/pre>\n<\/div>\n<\/details>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[],"tags":[],"class_list":["post-358961","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/358961","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=358961"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/358961\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=358961"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=358961"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=358961"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}