Kt_Irc_deal.c 62.3 KB
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <pthread.h>
#include <unistd.h>

#include "JZsdkLib.h"
#include "version_choose.h"
#include "BaseConfig.h"
#include "MediaProc/MediaParm.h"
#include <math.h>
#include "../../StreamProc/StreamProc.h"
#include "../../ImageProc/PseudoColor/PseudoColor.h"
#include "../../ImageProc/BadPixelProc/BadPixelProc.h"
#include "./Kt_Irc_deal.h"

#define SPC_DATA_DIR "/root/spc_save_data.raw"

extern int JZsdk_RawBadPixelCheck_Flag;  //原始盲点记录标志
extern int JZsdk_GrayBadPixelCheck_flag; //灰度图坏点纠正标志
extern int JZsdk_SPCCheck_flag; // 单点纠正标志
extern int JZsdk_GasEnhancementCorrect_flag ; //气体增强纠正标志


extern int g_CameraRegionBox_FirstPixel[2];  
extern int g_CameraRegionBox_SecondPixel[2]; 
extern int g_CameraRegionBox_ClickTime;

int g_temp_GasValueMax = 160;
int g_temp_GasValueMin = 30;

#define BAD_PIXEL_MODE JZ_FLAGCODE_OFF

static T_JZsdkReturnCode Kt_Irc_TempCalibrate_Read(int flag, U16_t *Image, unsigned int maxLen);
static T_JZsdkReturnCode Kt_Irc_TempCalibrate_Save(int flag, U16_t *Image, unsigned int maxLen);
static T_JZsdkReturnCode Kt_Irc_ReadAndSave_threshold(int flag, int order,unsigned int *num);


/**********************************************
 * 
 *  直方均衡, 提高图像的对比度,以降低过曝和曝光不足问题
 *  从直方图来看,即将 密集的直方图向两边扩散,从而变疏和均匀
 *  https://blog.csdn.net/qq_43416206/article/details/138680036
 * 
 * 
 *  输入源数据 直方图数据 输出承载数组 配置信息
 * *****************************************************/
static T_JZsdkReturnCode StraightBalancing_U16(U16_t *rawData, unsigned int *StraightDraw, U16_t *outData, struct Kt_Irc_deal *dealInfo )
{
    double result=0;

    for (int i = 0; i < dealInfo->PixelNum; i++)
    {
        result = 0.0;
        for (int temp = 0; temp < rawData[i]; temp++)
        {
            result += StraightDraw[temp];
        }

        result = (result*dealInfo->ExpectedMax) / dealInfo->PixelNum;
        outData[i] =  (U16_t)((int)result); 
    }

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}


/**********************************************
 * 
 *  直方拉伸
 *  https://blog.csdn.net/qq_43416206/article/details/138680036
 *  输入源数据 直方图数据 输出承载数组 配置信息
 * *****************************************************/
static T_JZsdkReturnCode StraightStretching(U16_t *rawData, unsigned int *StraightDraw, U16_t *outData, struct Kt_Irc_deal *dealInfo )
{
    int LeftDrop = dealInfo->LeftDrop;
    int RightDrop = dealInfo->RightDrop;
    int hist_cnt = 0;
    int find_min = 0, find_max = 0;
    int LowIndex = 0, HighIndex = 0;
    for (int i = 0; i < dealInfo->ExpectedMax; i++)
    {
        hist_cnt += StraightDraw[i];
        
        if ((hist_cnt >= LeftDrop) && (!find_min))
        {
            
            LowIndex = i;
            //JZSDK_LOG_INFO("lowindex:%d", LowIndex);
            find_min = 1;
        }
        if((hist_cnt >= (dealInfo->PixelNum - RightDrop) )&&(!find_max))
        {
            HighIndex = i;
            //JZSDK_LOG_INFO("HighIndex:%d", HighIndex);
            find_max = 1;
        }
    }

    if (HighIndex == LowIndex)
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }
    
    for (int i = 0; i < dealInfo->PixelNum; i++)
    {
        if (rawData[i]<= LowIndex) 
        {
            outData[i]= 0;
        }
        else if (rawData[i]>= HighIndex) 
        {
            outData[i]= dealInfo->ExpectedMax;
        }
        else
        { 
            outData[i]= ((rawData[i]- LowIndex)*dealInfo->ExpectedMax) /(HighIndex- LowIndex);
        }
    }
    
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/*计算直方图均衡化并转换为8位灰度值
	    U16_t in_str 输入的数据
                int in_str_len 输入的数据长度
    灰度0为黑色,灰度255为白色
*/
static T_JZsdkReturnCode JZsdk_vKT(U16_t *in_str, U8_t **out_str, int *out_str_len, struct Kt_Irc_deal *dealInfo)
{
    //无图像
    if (in_str == NULL || dealInfo == NULL) 
    {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }  

    // 分配输出结果的内存空间
    *out_str = (unsigned char*)malloc(sizeof(unsigned char) * dealInfo->PixelNum);
    if (*out_str == NULL) {
        // 如果内存分配失败,打印错误信息并返回空指针
        JZSDK_LOG_ERROR("Error: Memory allocation failed\n");
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

/**********************************************
 * 
 *  获取图像的直方图
 * 
 * 
 * *****************************************************/

    // 初始化变量
    unsigned int ThisPixelGray = 0;     //当前处理像素点的灰度值
    long SumAllLegalGray = 0;       //存储的总灰度值

    // 分配当前灰度值出现次数的内存空间,并初始化为0
    //用于记录每个灰度值在输入向量中出现的次数。具体来说,它是一个无符号整型数组,
    //其索引表示灰度值,数组中的值表示该灰度值在输入向量中出现的次数。
    //在遍历输入向量时,将每个像素的灰度值记录在相应的 ThisGrayAppendTime 数组中。
    //这是直方图均衡化算法中需要用到的一个重要步骤,用于统计每个灰度值的出现次数 
    //为保证每一个灰度值都被包括,所以用 dealInfo->ExpectedMax 写长度
    unsigned int* ThisGrayAppendTime = (unsigned int*)calloc(dealInfo->ExpectedMax + 1, sizeof(unsigned int));
    if (ThisGrayAppendTime == NULL) {
        JZSDK_LOG_ERROR("Error: Memory allocation failed\n");
        free(*out_str);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    // 遍历输入向量,统计各灰度值出现次数和总灰度和
    for (int i = 0; i < dealInfo->PixelNum; ++i) 
    {
        //获取本像素点的灰度值,且避免超过最大灰度值
        ThisPixelGray = in_str[i];
        if (ThisPixelGray > dealInfo->ExpectedMax) 
        {
            ThisPixelGray = dealInfo->ExpectedMax;
        }

        //直方图的 x坐标ThisPixelGray 的y值 即次数+1
        ThisGrayAppendTime[ThisPixelGray]++;

        //计算总灰度值
        SumAllLegalGray += ThisPixelGray;
    }

/*
    注:图像对比度:直方图分布越均匀,图像的对比度越高
        图像亮度:  直方图分布越靠右(最大值),则亮度越高
*/

/********************
 * 
 * 直方拉伸
 * 
 * 
 * *****************/
    U16_t StrDraw[dealInfo->PixelNum];
    StraightStretching(in_str, ThisGrayAppendTime, StrDraw, dealInfo);

/****************************
 * 
 *  将14位图,右移为8位图
 * 
 * 
 * *****************************/
   
    // 将均衡化后的灰度值转换为8位,并存储在输出向量中
    //直方图均衡化后的结果进行操作,将每个值转换为8位灰度值。具体操作是将in_str[i]中的值右移(InBits - 8)位
    //然后将结果转换为unsigned char类型,存储到输出向量对应的位置上。这样操作的目的是将原本可能不是8位的灰度值转换为8位灰度值,以便后续的处理和存储。
    for (int i = 0; i < dealInfo->PixelNum; i ++) 
    {
        (*out_str)[i] = (unsigned char)(StrDraw[i] >> (dealInfo->ImgDataBits - 8));
        //(*out_str)[i] = (unsigned char)(in_str[i] >> (dealInfo->ImgDataBits - 8));
    }

    *out_str_len = dealInfo->PixelNum; 
   
   
    // 释放动态分配的内存空间
    free(ThisGrayAppendTime);

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/*计算直方图均衡化并转换为8位灰度值
	    U16_t in_str 输入的数据
                int in_str_len 输入的数据长度
    灰度0为黑色,灰度255为白色
*/
static T_JZsdkReturnCode Kt_Irc_Histogram_vKT(U16_t *in_str, U8_t **out_str, int *out_str_len, struct Kt_Irc_deal *dealInfo)
{
    //无图像
    if (in_str == NULL || dealInfo == NULL) 
    {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }  

    // 分配输出结果的内存空间
    *out_str = (unsigned char*)malloc(sizeof(unsigned char) * dealInfo->PixelNum);
    if (*out_str == NULL) {
        // 如果内存分配失败,打印错误信息并返回空指针
        JZSDK_LOG_ERROR("Error: Memory allocation failed\n");
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

/**********************************************
 * 
 *  获取图像的直方图
 * 
 * 
 * *****************************************************/

    // 初始化变量
    unsigned int ThisPixelGray = 0;     //当前处理像素点的灰度值
    long SumAllLegalGray = 0;       //存储的总灰度值

    // 分配当前灰度值出现次数的内存空间,并初始化为0
    //用于记录每个灰度值在输入向量中出现的次数。具体来说,它是一个无符号整型数组,
    //其索引表示灰度值,数组中的值表示该灰度值在输入向量中出现的次数。
    //在遍历输入向量时,将每个像素的灰度值记录在相应的 ThisGrayAppendTime 数组中。
    //这是直方图均衡化算法中需要用到的一个重要步骤,用于统计每个灰度值的出现次数 
    //为保证每一个灰度值都被包括,所以用 dealInfo->ExpectedMax 写长度
    unsigned int* ThisGrayAppendTime = (unsigned int*)calloc(dealInfo->ExpectedMax + 1, sizeof(unsigned int));
    if (ThisGrayAppendTime == NULL) {
        JZSDK_LOG_ERROR("Error: Memory allocation failed\n");
        free(*out_str);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    // 遍历输入向量,统计各灰度值出现次数和总灰度和
    for (int i = 0; i < dealInfo->PixelNum; ++i) 
    {
        //获取本像素点的灰度值,且避免超过最大灰度值
        ThisPixelGray = in_str[i];
        if (ThisPixelGray > dealInfo->ExpectedMax) 
        {
            ThisPixelGray = dealInfo->ExpectedMax;
        }

        //直方图的 x坐标ThisPixelGray 的y值 即次数+1
        ThisGrayAppendTime[ThisPixelGray]++;

        //计算总灰度值
        SumAllLegalGray += ThisPixelGray;
    }

/*
    注:图像对比度:直方图分布越均匀,图像的对比度越高
        图像亮度:  直方图分布越靠右(最大值),则亮度越高
*/

/**********************************************************************************************************/

    // 初始化边界值和阈值增益
    unsigned int LowLimit = 0; //用于指定直方图均衡化的最低灰度值范围,通常为0。
    unsigned int HighLimit = 0; //用于指定直方图均衡化的最高灰度值范围,通常为最大灰度值。
    unsigned int VMax = dealInfo->ExpectedMax; //表示输入图像中可能的最大灰度值。
    unsigned int VMin = 0; //通常为0,表示输入图像中可能的最小灰度值。
    unsigned int ThresholdLimits = dealInfo->Gain * dealInfo->ImgDataBits; //阈值参数,是一个根据增益和位深度计算得出的阈值,用于在直方图均衡化时进行限制,确保输出图像的对比度合适。
    double AvgGray = SumAllLegalGray / (dealInfo->PixelNum * 1.0); //平均灰度值,表示输入图像中所有合法灰度值的平均值,用于直方图均衡化算法中计算每个像素点的增益值。

    // 计算灰度直方图的边界值,找到合适的灰度范围用于直方图均衡化
    for (int i = 0; i < dealInfo->ExpectedMax; ++i) 
    {
        if (LowLimit < dealInfo->LeftDrop) {
            LowLimit += ThisGrayAppendTime[i];
            VMin = i;
        }
        if (HighLimit < (dealInfo->PixelNum - dealInfo->RightDrop)) 
        {
            HighLimit += ThisGrayAppendTime[i];
            VMax = i;
        }
    }

    //(ThresholdLimits / 12) 是一个阈值的限制
    //根据增益和位深度计算得出的阈值,用于在直方图均衡化时进行限制,确保输出图像的对比度合适

    //将 VMax 增加 (ThresholdLimits / 12) 的目的是为了保证直方图均衡化后的最高灰度值范围与阈值限制的要求相符,
    //从而确保输出图像的对比度得到合适的控制,并且不会超出预期的范围。
    if (VMax < dealInfo->ExpectedMax - (ThresholdLimits / 12)) 
    {
        VMax += (ThresholdLimits / 12); //最大范围右移 1/12个阈值
    }

    //ThresholdDiff 阈值差异,计算阈值和灰度均值之间的差值
    //ThresholdDiff 被计算为 VMax - VMin + (ThresholdLimits / 6.0)
    //代表最高和最低灰度值范围之间的差距,再加上阈值的六分之一
    //ThresholdDiff 用于衡量最大最小灰度范围的宽度,以确定是否需要对像素进行进一步处理。
    double ThresholdDiff = VMax - VMin + (ThresholdLimits / 6.0);

    int LowThreshold = 0;
    if (AvgGray > (ThresholdDiff / 2.0)) 
    {
        //a. 如果 AvgGray 大于 ThresholdDiff / 2.0,则将 LowThreshold 的值设为 AvgGray - (ThresholdDiff / 2.0) 
        //的无符号整数部分。让 LowThreshold 在保证图像对比度的同时,尽可能地靠近平均灰度值 AvgGray。
        LowThreshold = (unsigned int)(AvgGray - (ThresholdDiff / 2.0));
    } 
    else 
    {
        //如果 AvgGray 小于等于 ThresholdDiff / 2.0,则将 LowThreshold 的值设为零。
        //避免 LowThreshold 取负值,确保不会对图像的过滤效果产生不良影响。
        LowThreshold = 0;
    }

    //重新计算阈值差值
    ThresholdDiff = VMax - LowThreshold;

    //避免阈值差值低于阈值限制,
    if (ThresholdDiff < ThresholdLimits) 
    {
        ThresholdDiff = ThresholdLimits;
        
        //根据新的阈值差值,重新计算 低阈值    
        LowThreshold = (unsigned int)(AvgGray - ThresholdDiff / 2.0);
    
        //最大灰度值与低阈值加上阈值限制之间的差异,如果仍然满足条件
        if (VMax > LowThreshold + ThresholdLimits) 
        {
            //重新计算阈值差值
            ThresholdDiff = VMax - LowThreshold;
        }
    }

    // 对输入向量进行直方图均衡化处理,并转换为8位灰度值
    for (int i = 0; i < dealInfo->PixelNum; ++i) 
    {
        //如果该像素点的灰度值小于最小阈值,使其等于最小阈值
        //确保直方图均衡化的值不会低于设定的阈值,从而保证图像的对比度和亮度的均衡。
        if (in_str[i] < LowThreshold) 
        {
            in_str[i] = LowThreshold;
        }

        //原始直方图值in_str[i]减去低阈值LowThreshold
        //然后乘以一个缩放系数(dealInfo->ExpectedMax / (ThresholdDiff * 1.0)),将结果赋值回in_str[i]。
        //目的将直方图的灰度值映射到更广的范围,以增强图像的对比度和细节。
        in_str[i] = (unsigned int)((in_str[i] - LowThreshold) * (dealInfo->ExpectedMax / (ThresholdDiff * 1.0)));

        //避免超出灰度最大值
        if (in_str[i] > dealInfo->ExpectedMax) 
        {
            in_str[i] = dealInfo->ExpectedMax;
        }
    }

    // 将均衡化后的灰度值转换为8位,并存储在输出向量中
    //直方图均衡化后的结果进行操作,将每个值转换为8位灰度值。具体操作是将in_str[i]中的值右移(InBits - 8)位
    //然后将结果转换为unsigned char类型,存储到输出向量对应的位置上。这样操作的目的是将原本可能不是8位的灰度值转换为8位灰度值,以便后续的处理和存储。
    for (int i = 0; i < dealInfo->PixelNum; i += 4) {
        (*out_str)[i]     = (unsigned char)(in_str[i] >> (dealInfo->ImgDataBits - 8));
        (*out_str)[i + 1] = (unsigned char)(in_str[i + 1] >> (dealInfo->ImgDataBits - 8));
        (*out_str)[i + 2] = (unsigned char)(in_str[i + 2] >> (dealInfo->ImgDataBits - 8));
        (*out_str)[i + 3] = (unsigned char)(in_str[i + 3] >> (dealInfo->ImgDataBits - 8));
    }

    *out_str_len = dealInfo->PixelNum;

    // 释放动态分配的内存空间
    free(ThisGrayAppendTime);

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}


/*3. (组成RGB数组) static vector<unsigned char> supportFun_RGB(const vector<unsigned char> In) {}

// 	输入: vector<unsigned char>	即存有 8 位灰度值的容器
// 	输出: vector<unsigned char> 即可被转换的带有 RGB 值的容器,长度 245760 (81920 * 3)

// 	将输入数组的每个元素,扩展三倍,如:

// 	示例输入: 10 30 60 90 ...
// 	示例输出: 10 10 10 30 30 30 60 60 60 90 90 90 ...

// 	后续元素以此类推。
*/ 
T_JZsdkReturnCode Kt_Irc_GrayTo_RGB(U8_t *in_str, U8_t **out_str, unsigned int *out_str_len, struct Kt_Irc_deal *dealInfo) 
{
    *out_str_len = (dealInfo->PixelNum * 3);
    *out_str = (unsigned char*)malloc(*out_str_len * sizeof(unsigned char));
    if (out_str == NULL)
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
    
    //这个会导致y轴倒置, 需要转码后才正常
    for (int i = 0; i < *out_str_len; i += 3) 
    {
        (*out_str)[i] = in_str[i / 3];
        (*out_str)[i + 1] = in_str[i / 3];
        (*out_str)[i + 2] = in_str[i / 3];
    }

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/******
 * 
 *  盲元矫正
 *  输入变量,画面数据,处理地址
 * 
 * *****/
T_JZsdkReturnCode Kt_Irc_CutBlindPixel(U16_t *ImageData, struct Kt_Irc_deal *dealInfo) 
{  
    if (ImageData == NULL || dealInfo == NULL) 
    {  
        // 如果输入数据为空,则返回指针错误
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }  
   
    // 遍历盲元数组
    for (int i = 0; i < (dealInfo->BlindNum)*2; i+=2)
    {
        int PosX = dealInfo->Blind[i]; //盲元在图像上的x轴坐标
        int PosY = dealInfo->Blind[i+1]; //盲元在图像上的y轴坐标
        int Pos  = PosY * dealInfo->Width + PosX; //盲元在图像一维数组中的位置

        //检查盲元的3*3的位置是否处于图像内
        if( (PosX > 0) && (PosX < (dealInfo->Width-1)) && (PosY > 0) && (PosY < (dealInfo->Height - 1)) ) 
        {
            // 计算3x3邻域的平均值
            int Avged = (   ImageData[Pos - dealInfo->Width - 1] + ImageData[Pos - dealInfo->Width] + ImageData[Pos - dealInfo->Width + 1] +
                            ImageData[Pos - 1]                   +                                  + ImageData[Pos + 1]+
                            ImageData[Pos + dealInfo->Width - 1] + ImageData[Pos + dealInfo->Width] + ImageData[Pos + dealInfo->Width + 1]) / 8;  

            // 如果得出的平均值小于纠正范围Extern,并且位置在图像内部有足够边距,则计算5x5邻域的平均值 
            if ( (fabs(ImageData[Pos] - Avged) < dealInfo->BlindExtern) && (PosX > 1) && (PosX < (dealInfo->Width - 2)) && (PosY > 1) && (PosY < (dealInfo->Height - 2)) ) 
            {  
                int Sum = 0;  
                // 计算5x5邻域的和  
                for (int y = -2; y <= 2; y++) {  
                    for (int x = -2; x <= 2; x++) {  
                        int index = (PosY + y) * dealInfo->Width + (PosX + x);  
                        Sum += ImageData[index];  
                    }  
                }  

                //除去中间块
                Sum -= ImageData[Pos];

                // 计算5x5邻域的平均值  
                Avged = Sum / 24;  
            }  
        }
    }
    
    // 返回处理后的图像数据  
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;  
}  

/******
 * 
 *  增加盲元
 *  输入变量,画面数据,左上角坐标,选取区域的右下角坐标,处理地址
 *  视频画面
 * *****/
T_JZsdkReturnCode Kt_Irc_AddBlind(U16_t *ImageData,struct Kt_Irc_deal *dealInfo, int topLeft_X, int topLeft_Y, int bottomRight_X, int bottomRight_Y, int blindIsMax) 
{  
    // 检查红外图像数据是否为空
    if (ImageData == NULL || dealInfo == NULL) 
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
  
    // 检查两点是否重合,至少是面,不是线
    if (topLeft_X == bottomRight_X || topLeft_Y == bottomRight_Y) 
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }
  
    // 计算矩形区域的大小
    int areaWidth = bottomRight_X - topLeft_X + 1;
    int areaHeight = bottomRight_Y - topLeft_Y + 1;

    // 分配内存以存储矩形区域的数据
    U16_t *area = (U16_t *)malloc(areaWidth * areaHeight * sizeof(U16_t));
    if (area == NULL) 
    { 
        // 内存分配失败处理
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
  
    // 从画面数据中提取矩形区域的数据
    for (int y = 0; y < areaHeight; y++) 
    {
        memcpy(area + y * areaWidth, ImageData + (topLeft_Y + y) * dealInfo->Width + topLeft_X, areaWidth * sizeof(U16_t));
    }

    // 初始化盲元点的位置和灰度值
    // 更新盲元点数组
    int *newBlindPoints = (int *)realloc(dealInfo->Blind, (dealInfo->BlindNum + 1) * sizeof(int) * 2);
    if (newBlindPoints != NULL) 
    {  
        dealInfo->Blind = newBlindPoints;  
    } else {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    dealInfo->Blind[dealInfo->BlindNum*2] = 0; //x坐标
    dealInfo->Blind[dealInfo->BlindNum*2+1] = 0; //y坐标
    
    // 遍历矩形区域找到最大或最小灰度值的点
    int blindGray = area[0]; //比较用的基准灰度值
    for (int i = 0; i < areaWidth * areaHeight; i++) 
    {
        int gray = area[i];
        if ((gray > blindGray && blindIsMax) || (gray < blindGray && !blindIsMax)) 
        {
            blindGray = gray;
            dealInfo->Blind[dealInfo->BlindNum*2] = i % areaWidth + topLeft_X;
            dealInfo->Blind[dealInfo->BlindNum*2+1] = i / areaWidth + topLeft_Y;
        }
    }
    dealInfo->BlindNum++; //增加一个盲元点

    // 释放矩形区域数据的内存
    free(area);
}

/******
 * 
 *  删除某个盲元
 *  输入变量,处理地址
 *  
 * *****/
T_JZsdkReturnCode Kt_Irc_RemoveBlind(struct Kt_Irc_deal *dealInfo, int BlindNumindex) 
{  
    // 检查红外图像数据是否为空
    if (dealInfo == NULL) 
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
  
    // 检查是否存在盲元
    if (dealInfo->BlindNum == 0) 
    {
        JZSDK_LOG_ERROR("无盲元,不可删除");
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }

    // 检查删除的盲元是否处理数组内
    if (dealInfo->BlindNum < BlindNumindex || BlindNumindex < 1)
    {
        JZSDK_LOG_ERROR("删除盲元不处于盲元范围内");
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }
    
  
    // 将要删除元素之后的元素向前移动一个位置  
    for (unsigned int i = BlindNumindex; i < dealInfo->BlindNum - 1; i++) {  
        dealInfo->Blind[i * 2] = dealInfo->Blind[(i + 1) * 2];  
        dealInfo->Blind[i * 2 + 1] = dealInfo->Blind[(i + 1) * 2 + 1];  
    }  
  
    // 减少盲元数量  
    dealInfo->BlindNum--;  
  
    // 如果需要,重新分配内存以缩小数组大小(注意:这可能会导致性能下降,因为频繁的内存分配/释放)  
    // 如果不重新分配,那么Blind数组将包含一些未使用的空间  
    if (dealInfo->BlindNum > 0) 
    {  
        // 重新分配内存  
        unsigned int newSize = dealInfo->BlindNum * 2 * sizeof(unsigned int);  
        unsigned int *newBlind = (unsigned int *)realloc(dealInfo->Blind, newSize);  
        if (newBlind != NULL) 
        {  
            dealInfo->Blind = newBlind;  
        } else {  
            // realloc失败,处理错误(例如,可能不需要处理,因为原始内存仍然有效)  
            return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
        }  
    } 
    else 
    {  
        // 如果没有盲元了,完全释放内存  
        free(dealInfo->Blind);  
        dealInfo->Blind = NULL;  
    }  

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/******
 * 
 *  删除所有盲元
 *  输入变量,处理地址
 *  
 * *****/
T_JZsdkReturnCode Kt_Irc_RemoveALLBlind(struct Kt_Irc_deal *dealInfo) 
{  
    // 释放盲元数组占用的内存  
    if (dealInfo->Blind != NULL) 
    {  
        free(dealInfo->Blind);  
    }  
  
    // 重置相关变量  
    dealInfo->Blind = NULL;  
    dealInfo->BlindNum = 0;  
}

/******
 * 
 *  增加屏蔽或标色区域  
 *  输入变量,处理地址
 *  
 * *****/
static T_JZsdkReturnCode Kt_Irc_AddShielding_or_ColorCoded_Areas(U16_t *ImageData,struct Kt_Irc_deal *dealInfo, int topLeft_X, int topLeft_Y, int bottomRight_X, int bottomRight_Y, int IsEnable) 
{  
    // 计算区域内的像素数  
    int AreaPixelNum = abs((bottomRight_X - topLeft_X) * (bottomRight_Y - topLeft_Y));  
      
    // 遍历区域内的每个像素并进行处理  
    for (int y = topLeft_Y; y < bottomRight_Y; y++) 
    {  
        for (int x = topLeft_X; x < bottomRight_X; x++) 
        {  
            int index = y * dealInfo->Width + x; // 计算一维数组中的索引  
              
            if (IsEnable) 
            {  
                // 如果是标色(IsEnable为true),将像素值从1改为0  
                dealInfo->DGCE_Area[index] = 0;   
            } 
            else 
            {  
                // 如果是屏蔽(IsEnable为false),将像素值从0改为1  
                dealInfo->DGCE_Area[index] = 1;   
            }  
        }  
    }  
}

/******
 * 
 *  低温循环标定
 *  手动划定区域屏蔽
 *  
 * *****/
static T_JZsdkReturnCode Kt_Irc_LowT_CycleCalibration(U16_t *ImageData,struct Kt_Irc_deal *dealInfo) 
{
    if (ImageData == NULL || dealInfo == NULL)
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    //检测到开始获取标定温度,清空标定数组,并赋予第一组基准值
    if (dealInfo->LowT_flag == JZ_FLAGCODE_ON)
    {
        JZSDK_LOG_DEBUG("开始低温标定");
        memset(dealInfo->LowT_NineFrameAdd, 0, dealInfo->PixelNum);
        memcpy(dealInfo->LowT_NineFrameAdd, ImageData, dealInfo->PixelNum * sizeof(U16_t));
        dealInfo->LowT_flag++;
    }
    else if (dealInfo->LowT_flag == JZ_FLAGCODE_OFF)
    {
        //不是,你这不是都没打开吗
        return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
    }
    else
    {
        //否则累加存储数组
        for (int i = 0; i < (dealInfo->Width * dealInfo->Height); i++)
        {
            dealInfo->LowT_NineFrameAdd[i] += ImageData[i];
        }
        dealInfo->LowT_flag++;
    }
    
    if (dealInfo->LowT_flag > 9) //已经存储了9组数据
    {
        for (int i = 0; i < (dealInfo->Width * dealInfo->Height); i++)
        {
            dealInfo->LowT_NineFrame_Avg[i] = dealInfo->LowT_NineFrameAdd[i] / 9;
        }   

        //存放到本地
        Kt_Irc_TempCalibrate_Save(1, dealInfo->LowT_NineFrame_Avg, dealInfo->PixelNum);

        dealInfo->LowT_flag = JZ_FLAGCODE_OFF; //将标志位归位

        JZSDK_LOG_DEBUG("低温数据记录完成");
    }
    
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/******
 * 
 *  高温循环标定
 *  手动划定区域屏蔽
 *  
 * *****/
static T_JZsdkReturnCode Kt_Irc_HighT_CycleCalibration(U16_t *ImageData, struct Kt_Irc_deal *dealInfo) 
{
    if (ImageData == NULL || dealInfo == NULL)
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    //检测到开始获取标定温度,清空标定数组,并赋予第一组基准值
    if (dealInfo->HighT_flag == JZ_FLAGCODE_ON)
    {
        JZSDK_LOG_DEBUG("开始高温标定");
        memset(dealInfo->HighT_NineFrameAdd, 0, dealInfo->PixelNum);
        memcpy(dealInfo->HighT_NineFrameAdd, ImageData, dealInfo->PixelNum * sizeof(U16_t));
        dealInfo->HighT_flag++;
    }
    else if (dealInfo->HighT_flag == JZ_FLAGCODE_OFF)
    {
        //不是,你这不是都没打开吗
        return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
    }
    else
    {
        //否则累加存储数组
        for (int i = 0; i < (dealInfo->Width * dealInfo->Height); i++)
        {
            dealInfo->HighT_NineFrameAdd[i] += ImageData[i];
        }
        dealInfo->HighT_flag++;
    }
    
    if (dealInfo->HighT_flag > 9) //已经存储了9组数据
    {
        for (int i = 0; i < (dealInfo->Width * dealInfo->Height); i++)
        {
            dealInfo->HighT_NineFrame_Avg[i] = dealInfo->HighT_NineFrameAdd[i] / 9;
        }   

        Kt_Irc_TempCalibrate_Save(0, dealInfo->HighT_NineFrame_Avg, dealInfo->PixelNum);

        dealInfo->HighT_flag = JZ_FLAGCODE_OFF; //将标志位归位

        JZSDK_LOG_DEBUG("高温数据记录完成");
    }
    
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/******
 * 
 *  两点矫正
 *  
 *  
 * *****/
T_JZsdkReturnCode Kt_Irc_TPC(U16_t *ImageData,struct Kt_Irc_deal *dealInfo) 
{
    //无图像
    if (ImageData == NULL || dealInfo == NULL) 
    {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    } 
  
    // 如果tpc的斜率和tpc的截距都为0,重新计算斜率与斜距
    int allZeroSlope_flag = 1, allZeroDiff_flag = 1;  
    for (int i = 0; i < dealInfo->PixelNum; i++) 
    {  
        if (dealInfo->TPC_Slope[i] != 0) 
        {  
            allZeroSlope_flag = 0;  
        }  
        if (dealInfo->TPC_Diff[i] != 0) {  
            allZeroDiff_flag = 0;  
        }  
    }  
  
    //如果出现tpc的斜率和tpc的截距都为0,重新计算斜率与斜距
    if (allZeroSlope_flag && allZeroDiff_flag) 
    {  
        double AvgSingleFrame_LowT = 0, AvgSingleFrame_HighT = 0;  

        for (int i = 0; i < dealInfo->PixelNum; i++) 
        {  
            AvgSingleFrame_LowT += dealInfo->LowT_NineFrame_Avg[i];  
            AvgSingleFrame_HighT += dealInfo->HighT_NineFrame_Avg[i];  
        }  

        AvgSingleFrame_LowT = AvgSingleFrame_LowT / dealInfo->PixelNum;  
        AvgSingleFrame_HighT = AvgSingleFrame_HighT / dealInfo->PixelNum;  
  
        for (int i = 0; i < dealInfo->PixelNum; i++) 
        {  
            if (dealInfo->HighT_NineFrame_Avg[i] > dealInfo->LowT_NineFrame_Avg[i]) 
            {  
                dealInfo->TPC_Slope[i] = (AvgSingleFrame_HighT - AvgSingleFrame_LowT) / (dealInfo->HighT_NineFrame_Avg[i] - dealInfo->LowT_NineFrame_Avg[i]);  
            } 
            else 
            {  
                dealInfo->TPC_Slope[i] = 0;  
            }  
            dealInfo->TPC_Diff[i] = AvgSingleFrame_LowT - dealInfo->TPC_Slope[i] * dealInfo->LowT_NineFrame_Avg[i];  
        }  
    }  
  
    // 应用两点校正公式  
    for (int i = 0; i < dealInfo->PixelNum; i++) 
    {  
        ImageData[i] = (int)(dealInfo->TPC_Slope[i] * ImageData[i] + dealInfo->TPC_Diff[i]);  
        if (ImageData[i] < 0) 
        {  
            ImageData[i] = 0;  
        } 
        else if (ImageData[i] > dealInfo->ExpectedMax) 
        {  
            ImageData[i] = dealInfo->ExpectedMax;  
        }  
    }  
  
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS; // 返回校正后的灰度数组  
}


/****************
 * 
 *  单点纠正的画面获取/保存
 * 
 * ***************/
static T_JZsdkReturnCode Kt_Irc_SPC_FrameSaveOrRead(struct Kt_Irc_deal *dealInfo, int order, U16_t *MarkData)
{
    FILE  *fp = NULL;
    //order为0是获取
    if (order == 0)
    {
        fp = fopen(SPC_DATA_DIR, "rb"); 
        if (fp == NULL)
        {
            JZSDK_LOG_ERROR("spc 数据获取失败1");
            return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
        }

        // 从文件读取到数组  
        int read = fread(MarkData, sizeof(U16_t), dealInfo->PixelNum, fp);  
        if (read != dealInfo->PixelNum) { 
            if (feof(fp)) {  
                // 如果到达文件末尾且没有读取到足够的数据,可能是文件不完整  
                JZSDK_LOG_ERROR("Reached end of file before reading all data.\n");  
            } 
            else 
            { 
                JZSDK_LOG_ERROR("Failed to read data from file");  
            }
            fclose(fp);  
            JZSDK_LOG_ERROR("spc 数据获取失败2");
            return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
        }  
    }
    else if (order == 1)
    {
        // 打开文件以写入  
        fp = fopen(SPC_DATA_DIR, "wb");  
        if (fp == NULL) {  
            JZSDK_LOG_ERROR("Failed to open file for writing");  
            return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;  
        }  
    
        // 写入数据到文件  
        int write = fwrite(MarkData, sizeof(uint16_t), dealInfo->PixelNum, fp);
        if (write != dealInfo->PixelNum) {  
            JZSDK_LOG_ERROR("Failed to write data to file");  
            fclose(fp);  
            return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;  
        }  
    }
    else
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }
        
    // 关闭文件  
    fclose(fp); 
    fp = NULL;
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/*************************
 * 
 *  单点纠正数值纠正
 * 
 * *************************/
static T_JZsdkReturnCode Kt_Irc_SPC_Param(struct Kt_Irc_deal *dealInfo, U16_t *MarkData)
{
    // 如果图像数据和Markdata完全相等,则计算校正参数  
    double Avg = 0.0;  
    for (int i = 0; i < dealInfo->PixelNum; i++) 
    {  
        Avg += MarkData[i]; 
    }  

    Avg /= dealInfo->PixelNum;  

    for (int i = 0; i < dealInfo->PixelNum; i++) 
    {  
        if (MarkData[i] <= 1) 
        {  
            // 处理除数为零的情况  
            // 可以设置一个默认值,比如斜率设为1(不进行斜率校正),或者设置为一个特定的错误码  
            dealInfo->SPC_Slope[i] = Avg / 1.0 ;  
            dealInfo->SPC_Diff[i] = Avg - MarkData[i];   
        } else {  
            dealInfo->SPC_Slope[i] = Avg / MarkData[i];  
            dealInfo->SPC_Diff[i] = Avg - MarkData[i];  
        }  
    }  

    JZSDK_LOG_INFO("重新计算IPC参数 Avg%f",Avg);
}

/******
 * 
 *  单点矫正
 *  
 * *****/
T_JZsdkReturnCode Kt_Irc_SPC(U16_t *ImageData,struct Kt_Irc_deal *dealInfo) 
{
    //无图像
    if (ImageData == NULL || dealInfo == NULL) 
    {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }  

    // 应用校正参数 到图像数据
    for (int i = 0; i < dealInfo->PixelNum; i++) 
    {    
        double temp = ( ((dealInfo->SPC_Slope[i] * ImageData[i]) + (dealInfo->SPC_Diff[i] + ImageData[i])) / 2);
        if (temp <= dealInfo->ExpectedMax && temp >= 0)
        {
            ImageData[i] = (int)temp;
        }
        if (ImageData[i] < 0) 
        {  
            ImageData[i] = 0;  
        } 
        else if (ImageData[i] > dealInfo->ExpectedMax) 
        {  
            ImageData[i] = dealInfo->ExpectedMax;  
        }  
    }  

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/******
 * 
 *  动态气体色彩增强
 *  手动划定区域屏蔽
 *  Imagedata rgb过后的图像
 *  rawdata 未经过spc处理的原数据
 * *****/
T_JZsdkReturnCode Kt_Irc_DynamicGasesColorEnhance(U8_t *ImageData, U16_t *RawData, struct Kt_Irc_deal *dealInfo) 
{  

    //气体增强前,需要先经过一次spc单点矫正,然后,再vkt直方图处理 imageData 是已经处理得到的rgb图

    //无图像
    if (ImageData == NULL || dealInfo == NULL) 
    {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }  

    // 遍历每个像素  
    for (int i = 0; i < dealInfo->PixelNum; i++) 
    {  
        int index = i * 3; // RGB值在一维数组中的索引  

        // 检查像素是否在指定区域内,并且差值在指定范围内  
        if (dealInfo->DGCE_Area[i] == 0
                && ( fabs(RawData[i] - dealInfo->DGCE_Mark_Data[i]) >= g_temp_GasValueMin )
                && (fabs(RawData[i] - dealInfo->DGCE_Mark_Data[i]) <= g_temp_GasValueMax) ) 
        {  
            // 根据标色类型设置RGB值  
            switch (dealInfo->DGCE_ColorType) {  
                case 0: // 红色  
                    ImageData[index + 0] = 255;  
                    ImageData[index + 1] = 0;  
                    ImageData[index + 2] = 0;  
                    break;  
                case 1: // 绿色  
                    ImageData[index + 0] = 0;  
                    ImageData[index + 1] = 255;  
                    ImageData[index + 2] = 0;
                case 2: // 蓝色
                    ImageData[index + 0] = 0;  
                    ImageData[index + 1] = 0;  
                    ImageData[index + 2] = 255;
                    break;
            }
        }
        
    }

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}



/********************
 * 
 *  Kt_Irc_deal 结构体的初始化
 * 
 * 
 * ******************/
T_JZsdkReturnCode Kt_Irc_deal_Init(struct Kt_Irc_deal **index, int height, int width)
{
    T_JZsdkReturnCode ret;

    struct Kt_Irc_deal *IrcDealCfg = NULL;

    //为Kt_Irc_deal 结构体注册内存
    IrcDealCfg = (struct Kt_Irc_deal *)malloc(sizeof(struct Kt_Irc_deal));  
    if (IrcDealCfg == NULL) {  
        JZSDK_LOG_ERROR("Kt_Irc_deal 结构体的初始化内存分配失败");  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }

    IrcDealCfg->ImgDataBits = 14; //图像为14位图像
    IrcDealCfg->Height = height;
    IrcDealCfg->Width  = width;
    IrcDealCfg->PixelNum = IrcDealCfg->Height * IrcDealCfg->Width;
    IrcDealCfg->ExpectedMax = (1 << IrcDealCfg->ImgDataBits) - 1;

    IrcDealCfg->DGCE_CurSubMarkedMin = 60;
    IrcDealCfg->DGCE_CurSubMarkedMax = 240; //240

    IrcDealCfg->Gain = 11;		//对比度,这里设为 11
    IrcDealCfg->LeftDrop = 5000;	//低灰度舍弃点数,这里设为 5000
    IrcDealCfg->RightDrop = 5000;	//高灰度舍弃点数,这里设为 5000

    IrcDealCfg->OutputPixelColorMode = 0;
    IrcDealCfg->First_DGCE_flag = JZ_FLAGCODE_ON;
    IrcDealCfg->SPC_flag = JZ_FLAGCODE_ON;

    IrcDealCfg->LowT_flag = JZ_FLAGCODE_OFF;
    IrcDealCfg->HighT_flag = JZ_FLAGCODE_OFF;
    IrcDealCfg->TPC_flag = JZ_FLAGCODE_OFF;

    IrcDealCfg->TPC_Diff = (double *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(double));  
    if (IrcDealCfg->TPC_Diff == NULL)
    {
        JZSDK_LOG_ERROR("TPC_Diff注册失败");
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
    
    IrcDealCfg->TPC_Slope = (double *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(double));  
    if (IrcDealCfg->TPC_Slope == NULL)
    {
        JZSDK_LOG_ERROR("TPC_Slope注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->HighT_NineFrameAdd = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->HighT_NineFrameAdd == NULL)
    {
        JZSDK_LOG_ERROR("HighT_NineFrameAdd注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->LowT_NineFrameAdd = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->LowT_NineFrameAdd == NULL)
    {
        JZSDK_LOG_ERROR("LowT_NineFrameAdd注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
    
    IrcDealCfg->HighT_NineFrame_Avg = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->HighT_NineFrame_Avg == NULL)
    {
        JZSDK_LOG_ERROR("HighT_NineFrame_Avg注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
    
    IrcDealCfg->LowT_NineFrame_Avg = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->LowT_NineFrame_Avg == NULL)
    {
        JZSDK_LOG_ERROR("LowT_NineFrame_Avg注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg->HighT_NineFrame_Avg);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->SPC_Diff = (double *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(double));  
    if (IrcDealCfg->SPC_Diff == NULL)
    {
        JZSDK_LOG_ERROR("SPC_Diff注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg->HighT_NineFrame_Avg);
        free(IrcDealCfg->LowT_NineFrame_Avg);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->SPC_Slope = (double *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(double));  
    if (IrcDealCfg->SPC_Slope == NULL)
    {
        JZSDK_LOG_ERROR("SPC_Slope注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg->HighT_NineFrame_Avg);
        free(IrcDealCfg->LowT_NineFrame_Avg);
        free(IrcDealCfg->SPC_Diff);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->DGCE_Area = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->DGCE_Area == NULL)
    {
        JZSDK_LOG_ERROR("DGCE_Area注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg->HighT_NineFrame_Avg);
        free(IrcDealCfg->LowT_NineFrame_Avg);
        free(IrcDealCfg->SPC_Diff);
        free(IrcDealCfg->SPC_Slope);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->DGCE_Mark_Data = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->DGCE_Mark_Data == NULL)
    {
        JZSDK_LOG_ERROR("DGCE_Mark_Data注册失败");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg->HighT_NineFrame_Avg);
        free(IrcDealCfg->LowT_NineFrame_Avg);
        free(IrcDealCfg->SPC_Diff);
        free(IrcDealCfg->SPC_Slope);
        free(IrcDealCfg->DGCE_Area);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    IrcDealCfg->SPC_Mark_Data = (U16_t *)malloc(IrcDealCfg->Height * IrcDealCfg->Width * sizeof(U16_t));  
    if (IrcDealCfg->SPC_Mark_Data == NULL)
    {
        JZSDK_LOG_ERROR("SPC_Mark_Data");
        free(IrcDealCfg->TPC_Diff);
        free(IrcDealCfg->TPC_Slope);
        free(IrcDealCfg->HighT_NineFrameAdd);
        free(IrcDealCfg->LowT_NineFrameAdd);
        free(IrcDealCfg->HighT_NineFrame_Avg);
        free(IrcDealCfg->LowT_NineFrame_Avg);
        free(IrcDealCfg->SPC_Diff);
        free(IrcDealCfg->SPC_Slope);
        free(IrcDealCfg->DGCE_Area);
        free(IrcDealCfg->DGCE_Mark_Data);
        free(IrcDealCfg);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    //尝试往markdata中填入数据
    ret = Kt_Irc_SPC_FrameSaveOrRead(IrcDealCfg, 0, IrcDealCfg->SPC_Mark_Data);
    if (ret != JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS)
    {
        IrcDealCfg->FirstSPC_flag = JZ_FLAGCODE_ON;
    }
    else
    {
        JZSDK_LOG_INFO("存在本地的spc数据");

        //计算一次spc参数
        Kt_Irc_SPC_Param(IrcDealCfg, IrcDealCfg->SPC_Mark_Data);
    }

    //加载坏点与灰坏点纠正
    IrcDealCfg->BadPixel = (unsigned int *)malloc(sizeof(unsigned int) * IrcDealCfg->Width * IrcDealCfg->Height);
    IrcDealCfg->GrayPixel = (unsigned int *)malloc(sizeof(unsigned int) * IrcDealCfg->Width * IrcDealCfg->Height);
    JZsdk_RawGetBadPixel_U16(IrcDealCfg->BadPixel, &IrcDealCfg->BadPixelNum, IrcDealCfg->Height*IrcDealCfg->Width);
    JZsdk_GrayGetBadPixel_U8(IrcDealCfg->GrayPixel, &IrcDealCfg->GrayPixelNum, IrcDealCfg->Height*IrcDealCfg->Width);
    IrcDealCfg->BadPixelExtern = 1000;
    IrcDealCfg->GrayPixelExtern = 100;

    //加载高低温
    ret = Kt_Irc_TempCalibrate_Read(0, IrcDealCfg->HighT_NineFrame_Avg, IrcDealCfg->Width * IrcDealCfg->Height);
    if (ret != JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS)
    {
        JZSDK_LOG_ERROR("没有高温数据");
    }
    ret = Kt_Irc_TempCalibrate_Read(1, IrcDealCfg->LowT_NineFrame_Avg, IrcDealCfg->Width * IrcDealCfg->Height);
    if (ret != JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS)
    {
        JZSDK_LOG_ERROR("没有低温数据");
    }

    //加载气体阈值
    Kt_Irc_ReadAndSave_threshold(0, JZ_FLAGCODE_GET, &g_temp_GasValueMax);
    Kt_Irc_ReadAndSave_threshold(1, JZ_FLAGCODE_GET, &g_temp_GasValueMin);

    *index = IrcDealCfg;

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

T_JZsdkReturnCode Kt_Irc_WriteGraphical(U8_t *rgb_data, struct Kt_Irc_deal *dealInfo)
{
    //不画
    if (g_CameraRegionBox_ClickTime == 0)
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
    }

    //画十字
    if (g_CameraRegionBox_ClickTime == 1)
    {
        Stream_rgb888_WriteCross(rgb_data, 
                            dealInfo->Width, dealInfo->Height, 
                            g_CameraRegionBox_FirstPixel[0],
                            g_CameraRegionBox_FirstPixel[1],
                            255, 215, 0, 20, 20);

        return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
    }
    
    //画十字
    if (g_CameraRegionBox_ClickTime == 2)
    {
        Stream_rgb888_WriteRectangle(rgb_data,
                                dealInfo->Width, dealInfo->Height,
                                g_CameraRegionBox_FirstPixel[0],
                                g_CameraRegionBox_FirstPixel[1],
                                g_CameraRegionBox_SecondPixel[0],
                                g_CameraRegionBox_SecondPixel[1], 
                                255, 215, 0, 1);

        return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
    }
    
    return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
}

//将14位灰度图数据转换为原始码流
T_JZsdkReturnCode Kt_Irc_14bitGrayData_to_RawData(U8_t *data, int data_len, U8_t **raw_data, int *raw_data_len, struct Kt_Irc_deal *dealInfo, int rawType)
{
    T_JZsdkReturnCode ret;

    //无图像
    if (data == NULL || dealInfo == NULL) 
    {  
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;  
    }  

    U16_t u16_RawData[dealInfo->PixelNum];     //原始16位码流
    U16_t u16_CorrentData[dealInfo->PixelNum];
    unsigned int u16_data_len; //16位码流的数据长度

    //1、合成像素,u8转换合并成u16
    JZsdk_Merge_U8_to_U16_byReverse(data, data_len,u16_RawData, &u16_data_len);
    if (u16_data_len != dealInfo->PixelNum)
    {
        JZSDK_LOG_ERROR("像素合成失败");
    }
    
    //2、数据检查
    for (int i = 0; i < dealInfo->PixelNum; i++)
    {
        if (u16_RawData[i] > dealInfo->ExpectedMax)
        {
            u16_RawData[i] = dealInfo->ExpectedMax;
        }
    }   

    //3、如果开启了原始流坏点寻找
    if (JZsdk_RawBadPixelCheck_Flag == JZ_FLAGCODE_ON)
    {
        JZsdk_RawCheckisBadPixel_U16(u16_RawData, dealInfo->Width, dealInfo->Height, dealInfo->BadPixel, &dealInfo->BadPixelNum);
        JZsdk_RawBadPixelCheck_Flag = JZ_FLAGCODE_OFF;
    }
    //盲元纠正
    JZsdk_CutBadPixel_U16(u16_RawData, dealInfo->Width, dealInfo->Height, dealInfo->BadPixel, dealInfo->BadPixelNum,dealInfo->BadPixelExtern, 1);


/***************************************************************************
 * 
 *  图像标定记录
 * 
 * 
 * *******************************************************************/
    //3、如果开启了低温循环标定
    if (dealInfo->LowT_flag != JZ_FLAGCODE_OFF)
    {
        Kt_Irc_LowT_CycleCalibration(u16_RawData, dealInfo);
    }
    
    //4、如果开启了高温循环标定
    if (dealInfo->HighT_flag != JZ_FLAGCODE_OFF)
    {
        Kt_Irc_HighT_CycleCalibration(u16_RawData, dealInfo);
    }


    //单点校正spc如果未进行过数据的标定,则进行一次
    if (dealInfo->FirstSPC_flag == JZ_FLAGCODE_ON)
    {
        JZSDK_LOG_INFO("未进行过数据标定");
        memcpy(dealInfo->SPC_Mark_Data, u16_RawData, dealInfo->PixelNum * sizeof(U16_t));
        Kt_Irc_SPC_FrameSaveOrRead(dealInfo, 1, dealInfo->SPC_Mark_Data);
        Kt_Irc_SPC_Param(dealInfo, dealInfo->SPC_Mark_Data);
        dealInfo->FirstSPC_flag = JZ_FLAGCODE_OFF;
    }
    //如果开启spc重置标定
    if (JZsdk_SPCCheck_flag == JZ_FLAGCODE_ON)
    {
        JZSDK_LOG_INFO("SPC重置标定");
        memcpy(dealInfo->SPC_Mark_Data, u16_RawData, dealInfo->PixelNum * sizeof(U16_t));
        Kt_Irc_SPC_FrameSaveOrRead(dealInfo, 1, dealInfo->SPC_Mark_Data);
        Kt_Irc_SPC_Param(dealInfo, dealInfo->SPC_Mark_Data);
        JZsdk_SPCCheck_flag = JZ_FLAGCODE_OFF;
    }

    //标记气体增图像图像
    //如果开启气体增强时,未标记过
    if (dealInfo->OutputPixelColorMode == 2 && dealInfo->First_DGCE_flag == JZ_FLAGCODE_ON)
    {
        memcpy(dealInfo->DGCE_Mark_Data, u16_RawData, dealInfo->PixelNum * sizeof(U16_t) );
        dealInfo->First_DGCE_flag = JZ_FLAGCODE_OFF;
    }
    //如果开启气体增强重置标定
    if (JZsdk_GasEnhancementCorrect_flag == JZ_FLAGCODE_ON)
    {
        JZSDK_LOG_INFO("气体增强重置标定");
        memcpy(dealInfo->DGCE_Mark_Data, u16_RawData, dealInfo->PixelNum * sizeof(U16_t));
        JZsdk_GasEnhancementCorrect_flag = JZ_FLAGCODE_OFF;
    }

/***************************************************************************
 * 
 *  校正
 * 
 * 
 * *******************************************************************/
    //5、校正
    //先复制一份纠正用数据
    memcpy(u16_CorrentData, u16_RawData, dealInfo->PixelNum * sizeof(U16_t));


    //如果打开了两点校正
    if (dealInfo->TPC_flag == JZ_FLAGCODE_ON)
    {
        Kt_Irc_TPC(u16_CorrentData, dealInfo);
    }

    
    //如果打开了单点校正
    if (dealInfo->SPC_flag == JZ_FLAGCODE_ON)
    {
        Kt_Irc_SPC(u16_CorrentData, dealInfo);
    }

    //7、图像输出模式
    U8_t *GrayImage = NULL; //灰度图数组
    unsigned int GrayImageLen = 0; //灰度图长度
    U8_t *RgbImage = NULL; //rgb888的图像数组
    unsigned int RgbImageLen = 0; //rgb数组的长度

    switch (dealInfo->OutputPixelColorMode)
    {
        case 0: //默认输出模式
        {
            //u16转换为灰度图
            //Kt_Irc_Histogram_vKT(u16_CorrentData, &GrayImage, &GrayImageLen, dealInfo);
            JZsdk_vKT(u16_CorrentData, &GrayImage, &GrayImageLen, dealInfo);
            
            if (JZsdk_GrayBadPixelCheck_flag == JZ_FLAGCODE_ON)
            {
                JZsdk_GrayCheckisBadPixel_U8(GrayImage, dealInfo->Width, dealInfo->Height, dealInfo->GrayPixel, &dealInfo->GrayPixelNum);
                JZsdk_GrayBadPixelCheck_flag = JZ_FLAGCODE_OFF;
            }
                
            //JZsdk_CutBadPixel_U8(GrayImage, dealInfo->Width, dealInfo->Height, dealInfo->GrayPixel, dealInfo->GrayPixelNum,dealInfo->GrayPixelExtern, 1);
            
            //灰度图转rgb888
            Kt_Irc_GrayTo_RGB(GrayImage, &RgbImage, &RgbImageLen,dealInfo);
        }
        break;

        case 1: //伪彩输出模式
        {
            //u16转换为灰度图
            //Kt_Irc_Histogram_vKT(u16_CorrentData, &GrayImage, &GrayImageLen, dealInfo);
            JZsdk_vKT(u16_CorrentData, &GrayImage, &GrayImageLen, dealInfo);

            if (JZsdk_GrayBadPixelCheck_flag == JZ_FLAGCODE_ON)
            {
                JZsdk_GrayCheckisBadPixel_U8(GrayImage, dealInfo->Width, dealInfo->Height, dealInfo->GrayPixel, &dealInfo->GrayPixelNum);
                JZsdk_GrayBadPixelCheck_flag = JZ_FLAGCODE_OFF;
            }

            //JZsdk_CutBadPixel_U8(GrayImage, dealInfo->Width, dealInfo->Height, dealInfo->GrayPixel, dealInfo->GrayPixelNum,dealInfo->GrayPixelExtern, 1);

            //灰度图转伪彩rgb888
            ret = PseudoColor_Gray2Rgb(GrayImage, &RgbImage, &RgbImageLen, dealInfo->PixelNum);
            if (ret != JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS)
            {
                Kt_Irc_GrayTo_RGB(GrayImage, &RgbImage, &RgbImageLen,dealInfo);
            }
            
        }
        break;

        case 2: //气体色彩增强输出模式
        {            
            //进行前需要自动打开单点纠正模式
            //u16转换为灰度图
            //Kt_Irc_Histogram_vKT(u16_CorrentData, &GrayImage, &GrayImageLen, dealInfo);
            JZsdk_vKT(u16_CorrentData, &GrayImage, &GrayImageLen, dealInfo); 

            // if (JZsdk_GrayBadPixelCheck_flag == JZ_FLAGCODE_ON)
            // {
            //     JZsdk_GrayCheckisBadPixel_U8(GrayImage, dealInfo->Width, dealInfo->Height, dealInfo->GrayPixel, &dealInfo->GrayPixelNum);
            //     JZsdk_GrayBadPixelCheck_flag = JZ_FLAGCODE_OFF;
            // }
            // JZsdk_CutBadPixel_U8(GrayImage, dealInfo->Width, dealInfo->Height, dealInfo->GrayPixel, dealInfo->GrayPixelNum,dealInfo->GrayPixelExtern, 1);
            
            //转为rgb
            Kt_Irc_GrayTo_RGB(GrayImage, &RgbImage, &RgbImageLen,dealInfo);

            //灰度图转气体增强rgb888
            Kt_Irc_DynamicGasesColorEnhance(RgbImage, u16_RawData, dealInfo);
        }
        break;
    
        default:
            return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
        break;
    }

    // 在rgb图上画图形
    Kt_Irc_WriteGraphical(RgbImage, dealInfo);

/*****************
 * 
 * 
 *  如果码流类型要求为yuv420p
 * 
 * **********************/
    if (rawType == 0)
    {
        //rgb888转yuv
        U8_t *yuv_frame = (U8_t *)malloc(dealInfo->Width*dealInfo->Height*3/2 ); 
        if (yuv_frame == NULL)
        {
            printf("内存注册失败\n");
            if (GrayImage != NULL)
            {
                free(GrayImage);
                GrayImage = NULL;
            }
            if (RgbImage != NULL)
            {
                free(RgbImage);
                RgbImage = NULL;
            }
            return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
        }

        // 调用转换函数 
        Stream_rgb888_to_yuv420p(RgbImage,dealInfo->Width,dealInfo->Height,yuv_frame);
        
        // 释放临时缓冲区,保留yuv_frame  
        if (GrayImage != NULL)
        {
            free(GrayImage);
            GrayImage = NULL;
        }
        if (RgbImage != NULL)
        {
            free(RgbImage);
            RgbImage = NULL;
        }

        *raw_data = yuv_frame;
        *raw_data_len = dealInfo->Width*dealInfo->Height*3/2;
    }

/*****************
 * 
 * 
 *  如果码流类型要求为rgb888
 * 
 * **********************/
    if (rawType == 1)
    {
        // 释放临时缓冲区
        if (GrayImage != NULL)
        {
            free(GrayImage);
            GrayImage = NULL;
        }

        *raw_data = RgbImage;
        *raw_data_len = dealInfo->Width*dealInfo->Height*3;
    }
    
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}


#define TEMP_CALIBRATE_SAVE_DIR "/root/TempCalibrate"
#define TEMP_CALIBRATE_HIGH_PIXEL_DIR "tempCalibrateHigh.txt"
#define TEMP_CALIBRATE_LOW_PIXEL_DIR "tempCalibrateLow.txt"

/**********************************
 * 
 *  高低温参数的标定值读取
 *  flag == 0 为高温标定
 *  flag == 1 为低温标定
 * ********************************/
static T_JZsdkReturnCode Kt_Irc_TempCalibrate_Read(int flag, U16_t *Image, unsigned int maxLen)
{
    unsigned char str[256];
    if (flag == 0)
    {
        snprintf(str, sizeof(str), "%s/%s",TEMP_CALIBRATE_SAVE_DIR, TEMP_CALIBRATE_HIGH_PIXEL_DIR);
    }
    else if (flag == 1)
    {
        snprintf(str, sizeof(str), "%s/%s",TEMP_CALIBRATE_SAVE_DIR, TEMP_CALIBRATE_LOW_PIXEL_DIR);
    }
    else
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }
    
    JZSDK_LOG_INFO("get_to%s",str);
    FILE  *Save_fp = fopen(str, "r");
    if (Save_fp == NULL)
    {
        JZSDK_LOG_ERROR("高低温参数的标定值打开失败"); 
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    if(fread(Image, sizeof(U16_t), maxLen,Save_fp) != maxLen)
    {
        JZSDK_LOG_ERROR("高低温参数的标定值读取失败"); 
        fclose(Save_fp);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    fclose(Save_fp);
}

/**********************************
 * 
 *  高低温参数的标定值存储
 * 
 * 
 * ********************************/
static T_JZsdkReturnCode Kt_Irc_TempCalibrate_Save(int flag, U16_t *Image, unsigned int maxLen)
{
    T_JZsdkReturnCode ret;
    //2、保存
    //检查文件夹
    ret = JZsdk_check_directory_exists_posix(TEMP_CALIBRATE_SAVE_DIR); 
    if (ret != JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS)
    {

        JZSDK_LOG_ERROR("不存在高低温参数文件夹,建立一个");
        ret = JZsdk_Osal_Mkdir(TEMP_CALIBRATE_SAVE_DIR);
        if (ret == JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE)
        {
            JZSDK_LOG_ERROR("高低温参数文件夹建立失败");         
            return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
        }
    }
    
    //写入文件
    unsigned char str[256];
    if (flag == 0)
    {
        snprintf(str, sizeof(str), "%s/%s",TEMP_CALIBRATE_SAVE_DIR, TEMP_CALIBRATE_HIGH_PIXEL_DIR);
    }
    else if (flag == 1)
    {
        snprintf(str, sizeof(str), "%s/%s",TEMP_CALIBRATE_SAVE_DIR, TEMP_CALIBRATE_LOW_PIXEL_DIR);
    }
    else
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }

    JZSDK_LOG_INFO("save_to%s",str);
    FILE  *Save_fp = fopen(str, "w");
    if (Save_fp == NULL)
    {
        JZSDK_LOG_ERROR("高低温文件建立失败"); 
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }

    if(fwrite(Image, sizeof(U16_t), maxLen,Save_fp) != maxLen)
    {
        JZSDK_LOG_ERROR("高低温参数的标定值保存失败"); 
        fclose(Save_fp);
        return JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
    }
  
    fclose(Save_fp); // 关闭文件  
  
    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS; // 返回成功代码  

}



/***********************
 * 
 *  气体阈值
 *  flag 0指高阈值 1低阈值
 *  
 * ************************/
#define HIGH_THRESHOLD_DIR "/root/highThreshold.txt"
#define LOW_THRESHOLD_DIR "/root/lowThreshold.txt"

static T_JZsdkReturnCode Kt_Irc_ReadAndSave_threshold(int flag, int order,unsigned int *num)
{
    //写入文件	
    FILE *fp = NULL;

    if (flag == 0)
    {
        if (order == JZ_FLAGCODE_SET)
        {      
            fp = fopen(HIGH_THRESHOLD_DIR, "wb+");
            if (fp == NULL) {  
                JZSDK_LOG_ERROR("Failed to open highThreshold file for writing.\r\n");  
                return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE; // 替换为合适的错误码  
            }  
            
            int ret = fwrite((unsigned int *)num, sizeof(unsigned int), 1,fp);
            if (ret != sizeof(unsigned int)) 
            {
                JZSDK_LOG_ERROR("Write highThreshold error.\r\n");
            }
            else{
                JZSDK_LOG_INFO("Write highThreshold=%d\n",*num);
            }

            fclose(fp);
        }
        else if (order == JZ_FLAGCODE_GET)
        {
            *num = g_temp_GasValueMax;
            fp = fopen(HIGH_THRESHOLD_DIR, "rb+");
            if (fp == NULL) 
            {
                fp = fopen(HIGH_THRESHOLD_DIR, "wb+");
                if (fp == NULL) 
                {
                    JZSDK_LOG_ERROR("threshold创建失败");
                    JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
                }
                *num = g_temp_GasValueMax;
            } 
            else 
            {
                int ret = fseek(fp, 0, SEEK_SET);
                if (ret != 0) {
                    printf("Seek log count file error, ret: %d.\r\n", ret);
                }

                ret = fread((unsigned int *)num, sizeof(unsigned int), 1,fp);
                if (ret != 1) 
                {
                    printf("Read highThreshold error.\r\n");
                }
                else{
                    printf("Read highThreshold=%d\n",*num);
                }
            }
            fclose(fp);    
        }
    }
    else if (flag == 1)
    {
        if (order == JZ_FLAGCODE_SET)
        {      
            fp = fopen(LOW_THRESHOLD_DIR, "wb+");
            if (fp == NULL) {  
                JZSDK_LOG_ERROR("Failed to open lowThreshold file for writing.\r\n");  
                return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE; // 替换为合适的错误码  
            }  
            
            int ret = fwrite((unsigned int *)num, sizeof(unsigned int), 1,fp);
            if (ret != sizeof(unsigned int)) 
            {
                JZSDK_LOG_ERROR("Write lowThreshold error.\r\n");
            }
            else{
                JZSDK_LOG_INFO("Write lowThreshold=%d\n",*num);
            }

            fclose(fp);
        }
        else if (order == JZ_FLAGCODE_GET)
        {
            *num = g_temp_GasValueMin;
            fp = fopen(LOW_THRESHOLD_DIR, "rb+");
            if (fp == NULL) 
            {
                fp = fopen(LOW_THRESHOLD_DIR, "wb+");
                if (fp == NULL) 
                {
                    JZSDK_LOG_ERROR("threshold创建失败");
                    JZ_ERROR_SYSTEM_MODULE_CODE_INVALID_PARAMETER;
                }
                *num = g_temp_GasValueMin;
            } 
            else 
            {
                int ret = fseek(fp, 0, SEEK_SET);
                if (ret != 0) {
                    printf("Seek log count file error, ret: %d.\r\n", ret);
                }

                ret = fread((unsigned int *)num, sizeof(unsigned int), 1, fp);
                if (ret != 1) 
                {
                    printf("Read lowThreshold error.\r\n");
                }
                else{
                    printf("Read lowThreshold=%d\n",*num);
                }
            }
            fclose(fp);    
        }
    }
    else
    {
        return JZ_ERROR_SYSTEM_MODULE_CODE_FAILURE;
    }

    return JZ_ERROR_SYSTEM_MODULE_CODE_SUCCESS;
}

/***********************
 * 
 *  气体阈值
 *  flag 0指高阈值 1低阈值
 *  
 * ************************/
T_JZsdkReturnCode Set_GasThreshold(int flag, int num)
{
    int temp = 0;
    if (flag == 0)
    {
        temp = g_temp_GasValueMax + num;
        if (temp <= 0)
        {
            temp = 0;
        }
        g_temp_GasValueMax = temp;
        Kt_Irc_ReadAndSave_threshold(flag, JZ_FLAGCODE_SET, &temp);
    }
    else if (flag == 1)
    {
        temp = g_temp_GasValueMin + num;
        if (temp <= 0)
        {
            temp = 0;
        }
        g_temp_GasValueMin = temp;
        Kt_Irc_ReadAndSave_threshold(flag, JZ_FLAGCODE_SET, &temp);
    }
    
}