目前用到的国产红外探测器普遍均匀度较差、且存在较多坏点，为了不影响最终的成像质量一般都会对探测器输出的图像先进行均匀性矫正和坏点去除。

基本原理

坏点去除原理很简单就是用周围的像素值来代替坏点的像素值。首先需要判断一张图像中坏点的位置，用待标定红外相机拍摄不同温度的黑体图像，坏点对温度的响应明显区别于正常的像素点，将这些点坐标标记出来。坏点替换时一般根据实际情况用3*3或者5*5或者更大的开窗的像素平均值来替换中心坏点的值。在实际操作时用均值替换的效果不如用周围像素的中值替换，因为很多红外探测器坏点喜欢集中出现，一个坏点周围可能还有坏点，用中值替换可以减少周围坏点对替换后效果的影响。

FPGA实现
坏点替换算法和中值滤波或者均值滤波算法很相似，坏点算法仅对标记为坏点的开窗计算均值或中值进行替换，而中值滤波或是均值滤波则对整幅图都计算中值或均值进行替换。所以实现坏点算法的过程和滤波算法极为相似
首先建立3*3开窗的寄存器，将图像视频流用移位寄存器ip核缓存两行，这样三行图像就实现了并行输出了，和当前行输出一起缓存三次就形成了3*3开窗。注意这里缓存行的长度包含了行空闲，如果一行长度超出了移位寄存器（shift ram ip）的最大长度就用两个寄存器。FPGA的图像处理实现基本都建立在开窗基础之上。
f373089bc0c64eab86c70028a4d3abef~noop.image?_iz=58558&from=article.jpg

将坏点标记和对应的图像像素一起缓存，中间的那行的中间像素若标记为1则代表当前开窗为坏点开窗，进行中值或均值替换，若不为1则不进行替换。
1.下面是以中值滤波为基础的坏点替换算法：
开窗建立顶层：
  2. `timescale 1ns / 1ps
  3. //
  4. // Company:
  5. // Engineer:
  6. //
  7. // Create Date: 2022/07/19 14:08:30
  8. // Design Name:
  9. // Module Name: mid_wave
  10. // Project Name:
  11. // Target Devices:
  12. // Tool Versions:
  13. // Description:
  14. //
  15. // Dependencies:
  16. //
  17. // Revision:
  18. // Revision 0.01 - File Created
  19. // Additional Comments:
  20. //
  21. //
  24. module mid_wave(
  25.    input I_clk,
  26.    input I_reset,
  27.    input I_clk_test,
  28.    input I_frame_valid,
  29.    input I_line_valid,
  30.    input[31:0] I_par_Q,
  31.    input[15:0] I_video_in,
  32.    output wire O_frame_valid,
  33.    output wire O_line_valid,
  34.    output wire[15:0] O_video_out,
  35.    output O_hblank,
  36.    output O_vblank
  37.     );
  38.    
  39. reg [18:0] S_temp0_r0;
  40. reg [18:0] S_temp0_r1;
  41. reg [18:0] S_temp0_r2;
  42. reg [18:0] S_temp1_r0;
  43. reg [18:0] S_temp1_r1;
  44. reg [18:0] S_temp1_r2;
  45. reg [18:0] S_temp2_r0;
  46. reg [18:0] S_temp2_r1;
  47. reg [18:0] S_temp2_r2;
  48. wire[18:0] S_fifo0;
  49. wire[18:0] S_fifo1;
  50. wire[18:0] S_fifo0_1;
  51. wire[18:0] S_fifo1_1;
  52. wire[18:0] S_video_in;
  53. wire[15:0] S_video_out;
  54. wire S_blind_sig;
  56. assign S_video_in = {I_frame_valid,I_line_valid,S_blind_sig,I_video_in};
  57. assign O_hblank = (~O_line_valid) & O_frame_valid;
  58. assign O_vblank = ~O_frame_valid;
  59. assign S_blind_sig = I_par_Q[28];
  62. c_shift_ram_0 c_shift_ram_0_0(
  63.     .D(S_video_in),
  64.     .CLK(I_clk),
  65.     .Q(S_fifo0_1)
  66.   );
  67.   
  68. c_shift_ram_0 c_shift_ram_0_1(
  69.     .D(S_fifo0_1),
  70.     .CLK(I_clk),
  71.     .Q(S_fifo0)
  72.   );
  73.   
  74. c_shift_ram_0 c_shift_ram_0_2(
  75.     .D(S_fifo0),
  76.     .CLK(I_clk),
  77.     .Q(S_fifo1_1)
  78.   );
  79.   
  80. c_shift_ram_0 c_shift_ram_0_3(
  81.     .D(S_fifo1_1),
  82.     .CLK(I_clk),
  83.     .Q(S_fifo1)
  84.   );
  85.   
  86. mid_data mid_data_i(
  87.     //System Interfaces
  88.     .sclk(I_clk)                 ,
  89.     .rst_n(I_reset)           ,
  90.     .mat_row1(S_video_in)     ,
  91.     .mat_row2(S_fifo0)        ,
  92.     .mat_row3(S_fifo1)        ,
  93.     .O_frame(O_frame_valid)   ,
  94.     .O_line(O_line_valid)     ,
  95.     .O_data(O_video_out)         
  96. );
  98. ///
  101. endmodule
I_par_Q为坏点标记因为是matlab计算的浮点数所以是32位也可以只用1位来表示坏点。将坏点标记和行同步帧同步一起加到数据的高3位，再一起进行缓存便于最后输出。这里是640*512 16位红外图像，每一行加上行空闲一共有700多个像素时钟周期超出移位寄存器ip最大设置，所以每一行用了2个移位寄存器ip。
中值计算
  2. `timescale 1ns / 1ps
  3. //
  4. // Company:
  5. // Engineer:
  6. //
  7. // Create Date: 2022/07/19 14:24:10
  8. // Design Name:
  9. // Module Name: mid_data
  10. // Project Name:
  11. // Target Devices:
  12. // Tool Versions:
  13. // Description:
  14. //
  15. // Dependencies:
  16. //
  17. // Revision:
  18. // Revision 0.01 - File Created
  19. // Additional Comments:
  20. //
  21. //
  24. module mid_data(
  25.     //System Interfaces
  26.     input                   sclk            ,
  27.     input                   rst_n           ,
  28.     //Communication Interfaces
  29.     input           [18:0]  mat_row1        ,
  30.     input           [18:0]  mat_row2        ,
  31.     input           [18:0]  mat_row3        ,
  32.     output          reg     O_frame         ,
  33.     output            reg        O_line            ,
  34.     output  reg     [15:0]  O_data         
  35. );
  37. //========================================================================================\
  38. //**************Define Parameter and  Internal Signals**********************************
  39. //========================================================================================/
  41. reg                 [18:0]  mat_row1_1      ;
  42. reg                 [18:0]  mat_row2_1      ;
  43. reg                 [18:0]  mat_row3_1      ;
  44. reg                 [18:0]  mat_row1_2      ;
  45. reg                 [18:0]  mat_row2_2      ;
  46. reg                 [18:0]  mat_row3_2      ;
  47.                      
  48. reg                 [18:0]  max_h1          ;
  49. reg                 [18:0]  mid_h1          ;  
  50. reg                 [18:0]  min_h1          ;
  51. reg                 [18:0]  max_h2          ;
  52. reg                 [18:0]  mid_h2          ;  
  53. reg                 [18:0]  min_h2          ;  
  54. reg                 [18:0]  max_h3          ;
  55. reg                 [18:0]  mid_h3          ;  
  56. reg                 [18:0]  min_h3          ;
  57. reg                 [18:0]  min_max         ;
  58. reg                 [18:0]  mid_mid         ;
  59. reg                 [18:0]  max_min         ;      
  60.   
  63. //========================================================================================\
  64. //**************     Main      Code        **********************************
  65. //========================================================================================/
  66. always @(posedge sclk)
  67.     begin
  68.         mat_row1_1          <=          mat_row1;
  69.         mat_row2_1          <=          mat_row2;
  70.         mat_row3_1          <=          mat_row3;
  71.         mat_row1_2          <=          mat_row1_1;
  72.         mat_row2_2          <=          mat_row2_1;
  73.         mat_row3_2          <=          mat_row3_1;
  74.     end
  75.    
  76.    
  77. always @(posedge sclk or negedge rst_n)
  78. begin
  79.     if(rst_n == 1'b0)
  80.         max_h1              <=          18'd0;        
  81.     else if(mat_row1[15:0] >= mat_row1_1[15:0] && mat_row1[15:0] >= mat_row1_2[15:0])
  82.         max_h1              <=          mat_row1;
  83.     else if(mat_row1_1[15:0] >= mat_row1[15:0] && mat_row1_1[15:0] >= mat_row1_2[15:0])
  84.         max_h1              <=          mat_row1_1;
  85.     else
  86.         max_h1              <=          mat_row1_2;
  87. end
  89. always @(posedge sclk or negedge rst_n)
  90. begin
  91.     if(rst_n == 1'b0)
  92.         mid_h1              <=          18'd0;        
  93.     else if((mat_row1[15:0] >= mat_row1_1[15:0] && mat_row1_1[15:0] >= mat_row1_2[15:0]) || (mat_row1_2[15:0] >= mat_row1_1[15:0] && mat_row1_1[15:0] >= mat_row1[15:0]))
  94.         mid_h1              <=          mat_row1_1;
  95.     else if((mat_row1_1[15:0] >= mat_row1[15:0] && mat_row1[15:0] >= mat_row1_2[15:0]) || (mat_row1_2[15:0] >= mat_row1[15:0] && mat_row1[15:0] >= mat_row1_1[15:0]))
  96.         mid_h1              <=          mat_row1;
  97.     else
  98.         mid_h1              <=          mat_row1_2;
  99. end
  100.          
  101. always @(posedge sclk or negedge rst_n)
  102. begin
  103.     if(rst_n == 1'b0)
  104.         min_h1              <=          18'd0;        
  105.     else if(mat_row1[15:0] <= mat_row1_1[15:0] && mat_row1[15:0] <= mat_row1_2[15:0])
  106.         min_h1              <=          mat_row1;
  107.     else if(mat_row1_1[15:0] <= mat_row1[15:0] && mat_row1_1[15:0] <= mat_row1_2[15:0])
  108.         min_h1              <=          mat_row1_1;
  109.     else
  110.         min_h1              <=          mat_row1_2;  
  111. end
  113. always @(posedge sclk or negedge rst_n)
  114. begin
  115.     if(rst_n == 1'b0)
  116.         max_h2              <=          18'd0;        
  117.     else if(mat_row2[15:0] >= mat_row2_1[15:0] && mat_row2[15:0] >= mat_row2_2[15:0])
  118.         max_h2              <=          mat_row2;
  119.     else if(mat_row2_1[15:0] >= mat_row2[15:0] && mat_row2_1[15:0] >= mat_row2_2[15:0])
  120.         max_h2              <=          mat_row2_1;
  121.     else
  122.         max_h2              <=          mat_row2_2;
  123. end
  125. always @(posedge sclk or negedge rst_n)
  126. begin
  127.     if(rst_n == 1'b0)
  128.         mid_h2              <=          18'd0;        
  129.     else if((mat_row2[15:0] >= mat_row2_1[15:0] && mat_row2_1[15:0] >= mat_row2_2[15:0]) || (mat_row2_2[15:0] >= mat_row2_1[15:0] && mat_row2_1[15:0] >= mat_row2[15:0]))
  130.         mid_h2              <=          mat_row2_1;
  131.     else if((mat_row2_1[15:0] >= mat_row2[15:0] && mat_row2[15:0] >= mat_row2_2[15:0]) || (mat_row2_2[15:0] >= mat_row2[15:0] && mat_row2[15:0] >= mat_row2_1[15:0]))
  132.         mid_h2              <=          mat_row2;
  133.     else
  134.         mid_h2              <=          mat_row2_2;
  135. end
  136.          
  137. always @(posedge sclk or negedge rst_n)
  138. begin
  139.     if(rst_n == 1'b0)
  140.         min_h2              <=          18'd0;        
  141.     else if(mat_row2[15:0] <= mat_row2_1[15:0] && mat_row2[15:0] <= mat_row2_2[15:0])
  142.         min_h2              <=          mat_row2;
  143.     else if(mat_row2_1[15:0] <= mat_row2[15:0] && mat_row2_1[15:0] <= mat_row2_2[15:0])
  144.         min_h2              <=          mat_row2_1;
  145.     else
  146.         min_h2              <=          mat_row2_2;  
  147. end
  149. always @(posedge sclk or negedge rst_n)
  150. begin
  151.     if(rst_n == 1'b0)
  152.         max_h3              <=          18'd0;        
  153.     else if(mat_row3[15:0] >= mat_row3_1[15:0] && mat_row3[15:0] >= mat_row3_2[15:0])
  154.         max_h3              <=          mat_row3;
  155.     else if(mat_row3_1[15:0] >= mat_row3[15:0] && mat_row3_1[15:0] >= mat_row3_2[15:0])
  156.         max_h3              <=          mat_row3_1;
  157.     else
  158.         max_h3              <=          mat_row3_2;
  159. end
  161. always @(posedge sclk or negedge rst_n)
  162. begin
  163.     if(rst_n == 1'b0)
  164.         mid_h3              <=          18'd0;        
  165.     else if((mat_row3[15:0] >= mat_row3_1[15:0] && mat_row3_1[15:0] >= mat_row3_2[15:0]) || (mat_row3_2[15:0] >= mat_row3_1[15:0] && mat_row3_1[15:0] >= mat_row3[15:0]))
  166.         mid_h3              <=          mat_row3_1;
  167.     else if((mat_row3_1[15:0] >= mat_row3[15:0] && mat_row3[15:0] >= mat_row3_2[15:0]) || (mat_row3_2[15:0] >= mat_row3[15:0] && mat_row3[15:0] >= mat_row3_1[15:0]))
  168.         mid_h3              <=          mat_row3;
  169.     else
  170.         mid_h3              <=          mat_row3_2;
  171. end
  173.          
  174. always @(posedge sclk or negedge rst_n)
  175. begin
  176.     if(rst_n == 1'b0)
  177.         min_h3              <=          18'd0;        
  178.     else if(mat_row3[15:0] <= mat_row3_1[15:0] && mat_row3[15:0] <= mat_row3_2[15:0])
  179.         min_h3              <=          mat_row3;
  180.     else if(mat_row3_1[15:0] <= mat_row3[15:0] && mat_row3_1[15:0] <= mat_row3_2[15:0])
  181.         min_h3              <=          mat_row3_1;
  182.     else
  183.         min_h3              <=          mat_row3_2;
  184. end
  187. always @(posedge sclk or negedge rst_n)
  188. begin
  189.     if(rst_n == 1'b0)
  190.         min_max             <=          18'd0;
  191.     else if(max_h1[15:0] <= max_h2[15:0] && max_h1[15:0] <= max_h3[15:0])
  192.         min_max             <=          max_h1;
  193.     else if(max_h2[15:0] <= max_h1[15:0] && max_h2[15:0] <= max_h3[15:0])
  194.         min_max             <=          max_h2;
  195.     else
  196.         min_max             <=          max_h3;
  197. end
  198.          
  199.    
  200. always @(posedge sclk or negedge rst_n)
  201. begin
  202.     if(rst_n == 1'b0)
  203.         mid_mid             <=          18'd0;
  204.     else if((mid_h1[15:0] >= mid_h2[15:0] && mid_h2[15:0] >= mid_h3[15:0]) || (mid_h3[15:0] >= mid_h2[15:0] && mid_h2[15:0] >= mid_h1[15:0]))
  205.         mid_mid             <=          mid_h2;
  206.     else if((mid_h2[15:0] >= mid_h1[15:0] && mid_h1[15:0] >= mid_h3[15:0]) || (mid_h3[15:0] >= mid_h1[15:0] && mid_h1[15:0] >= mid_h2[15:0]))
  207.         mid_mid             <=          mid_h1;
  208.     else
  209.         mid_mid             <=          mid_h3;
  210. end
  212. always @(posedge sclk or negedge rst_n)
  213. begin
  214.     if(rst_n == 1'b0)
  215.         max_min             <=          18'd0;
  216.     else if(min_h1[15:0] <= min_h2[15:0] && min_h1[15:0] <= min_h3[15:0])
  217.         max_min             <=          min_h1;
  218.     else if(min_h2[15:0] <= min_h1[15:0] && min_h2[15:0] <= min_h3[15:0])
  219.         max_min             <=          min_h2;
  220.     else
  221.         max_min             <=          min_h3;
  222. end
  223.            
  224. always @(posedge sclk or negedge rst_n)
  225. begin
  226.     if(rst_n == 1'b0)
  227.       begin
  228.         O_frame             <=          1'd0;
  229.         O_line                <=            1'd0;
  230.         O_data              <=          16'd0;
  231.       end
  232.     else if((mid_mid[15:0] >= min_max[15:0] && min_max[15:0] >= max_min[15:0]) || (max_min[15:0] >= min_max[15:0] && min_max[15:0] >= mid_mid[15:0]))
  233.       begin
  234.         O_frame             <=          mat_row2_1[18];
  235.         O_line                <=            mat_row2_1[17];
  236.         if(mat_row2_1[16] == 1)
  237.             O_data              <=          min_max[15:0];
  238.         else
  239.             O_data              <=          mat_row2_1[15:0];
  240.       end
  241.     else if((min_max[15:0] >= mid_mid[15:0] && mid_mid[15:0] >= max_min[15:0]) || (max_min[15:0] >= mid_mid[15:0] && mid_mid[15:0] >= min_max[15:0]))
  242.       begin
  243.         O_frame             <=          mat_row2_1[18];
  244.         O_line                <=            mat_row2_1[17];
  245.         if(mat_row2_1[16] == 1)
  246.             O_data              <=          mid_mid[15:0];
  247.         else
  248.             O_data              <=          mat_row2_1[15:0];
  249.       end
  250.     else
  251.       begin
  252.         O_frame             <=          mat_row2_1[18];
  253.         O_line              <=          mat_row2_1[17];
  254.         if(mat_row2_1[16] == 1)
  255.             O_data              <=          max_min[15:0];
  256.         else
  257.             O_data              <=          mat_row2_1[15:0];
  258.       end
  259. end
  260.          
  261.         
  262. endmodule
以上两个模块构成了完整的坏点替换算法，可以用带有坏点的模拟图像输入至模块顶层即可输出替换后的结果，I_frame_valid为输入帧同步高有效，I_line_valid为输入行同步高有效，I_video_in为像素数据，I_clk为像素同步时钟输入这几个信号就可以得到输出结果。可以防止试试。
2.基于均值滤波的坏点替换模块
  2. `timescale 1ns / 1ps
  3. //
  4. // Company:
  5. // Engineer:
  6. //
  7. // Create Date: 2022/01/06 10:48:23
  8. // Design Name:
  9. // Module Name: blind_point
  10. // Project Name:
  11. // Target Devices:
  12. // Tool Versions:
  13. // Description:
  14. //参考均值滤波法，使用两个深度为1280的19位fifo级联，延时前两行，每行输出端使用两个寄存器暂存前两个图像数据，形成3*3寄存器矩阵
  15. //对3*3寄存器取均值后直接输出，可将行有效帧有效和盲元标志并入像素数据高三位，直接输入，计算结果直接输出不必做状态机，不必管输入的帧有效和行有效
  16. // Dependencies:
  17. //
  18. // Revision:
  19. // Revision 0.01 - File Created
  20. // Additional Comments:
  21. //
  22. //
  25. module blind_point(
  26.    input I_clk,
  27.    input I_reset,
  28.    input I_clk_test,
  29.    input I_frame_valid,
  30.    input I_line_valid,
  31.    input[15:0] I_video_in,
  32.    input[31:0] I_par_Q,
  33.    output wire O_frame_valid,
  34.    output wire O_line_valid,
  35.    output wire[15:0] O_video_out,
  36.    output O_hblank,
  37.    output O_vblank
  38.     );
  39.    
  40. reg [18:0] S_temp0_r0;
  41. reg [18:0] S_temp0_r1;
  42. reg [18:0] S_temp0_r2;
  43. reg [18:0] S_temp1_r0;
  44. reg [18:0] S_temp1_r1;
  45. reg [18:0] S_temp1_r2;
  46. reg [18:0] S_temp2_r0;
  47. reg [18:0] S_temp2_r1;
  48. reg [18:0] S_temp2_r2;
  49. wire S_blind_sig;
  50. wire[18:0] S_fifo0;
  51. wire[18:0] S_fifo1;
  52. wire [18:0] S_fifo0_1;
  53. wire [18:0] S_fifo1_1;
  54. wire [18:0] S_add;
  55. wire [18:0] S_video_in;
  56. wire [15:0] S_result;
  57. wire S_frame_valid, S_line_valid;
  58. wire [15:0] S_video_out;
  59. assign S_blind_sig = I_par_Q[28];
  60. assign S_video_in = {I_frame_valid,I_line_valid,S_blind_sig,I_video_in};
  61. assign O_hblank = (~O_line_valid) & O_frame_valid;
  62. assign O_vblank = ~O_frame_valid;
  65. c_shift_ram_0 c_shift_ram_0_0(
  66.     .D(S_video_in),
  67.     .CLK(I_clk),
  68.     .Q(S_fifo0_1)
  69.   );
  70.   
  71. c_shift_ram_0 c_shift_ram_0_1(
  72.     .D(S_fifo0_1),
  73.     .CLK(I_clk),
  74.     .Q(S_fifo0)
  75.   );
  76.   
  77. c_shift_ram_0 c_shift_ram_0_2(
  78.     .D(S_fifo0),
  79.     .CLK(I_clk),
  80.     .Q(S_fifo1_1)
  81.   );
  82.   
  83. c_shift_ram_0 c_shift_ram_0_3(
  84.     .D(S_fifo1_1),
  85.     .CLK(I_clk),
  86.     .Q(S_fifo1)
  87.   );
  90. ///
  91. always @(posedge I_clk or negedge I_reset)
  92. begin
  93. if(!I_reset)
  94.   begin
  95.    S_temp0_r2 <= 0;
  96.    S_temp0_r1 <= 0;
  97.    S_temp0_r0 <= 0;
  98.   end
  99. else
  100.   begin
  101.    S_temp0_r2 <= S_video_in;
  102.    S_temp0_r1 <= S_temp0_r2;
  103.    S_temp0_r0 <= S_temp0_r1;
  104.   end
  105. end
  107. always @(posedge I_clk or negedge I_reset)
  108. begin
  109. if(!I_reset)
  110.   begin
  111.    S_temp1_r2 <= 0;
  112.    S_temp1_r1 <= 0;
  113.    S_temp1_r0 <= 0;
  114.   end
  115. else
  116.   begin
  117.    S_temp1_r2 <= S_fifo0;
  118.    S_temp1_r1 <= S_temp1_r2;
  119.    S_temp1_r0 <= S_temp1_r1;
  120.   end
  121. end
  123. always @(posedge I_clk or negedge I_reset)
  124. begin
  125. if(!I_reset)
  126.   begin
  127.    S_temp2_r2 <= 0;
  128.    S_temp2_r1 <= 0;
  129.    S_temp2_r0 <= 0;
  130.   end
  131. else
  132.   begin
  133.    S_temp2_r2 <= S_fifo1;
  134.    S_temp2_r1 <= S_temp2_r2;
  135.    S_temp2_r0 <= S_temp2_r1;
  136.   end
  137. end
  140. assign S_add = (S_temp0_r2[15:0] + S_temp0_r1[15:0] + S_temp0_r0[15:0])
  141.                 + (S_temp1_r2[15:0] + S_temp1_r0[15:0])
  142.                 + (S_temp2_r2[15:0] + S_temp2_r1[15:0] + S_temp2_r0[15:0]);
  143. assign S_result =  S_add/8 ;              
  144. assign O_video_out = (S_temp1_r1[16]) ? S_result : S_temp1_r1[15:0];
  145. assign O_frame_valid = S_temp1_r1[18];
  146. assign O_line_valid = S_temp1_r1[17];
  149. endmodule
这个是用均值进行替换的3*3开窗坏点替换，可以直接输入图像试试。

来源：沃爱单片机