用IDL写中值滤波的代码

以下是用IDL语言实现中值滤波的代码：

; 定义中值滤波函数
FUNCTION MedianFilter, image, kernel_size
  ; 获取图像的行数和列数
  dims = SIZE(image)
  rows = dims[1]
  cols = dims[2]
  
  ; 计算中值滤波核的半径
  radius = (kernel_size - 1) / 2
  
  ; 初始化输出图像
  output = FLTARR(rows, cols)
  
  ; 遍历图像中每个像素
  FOR i=radius, rows-1-radius DO BEGIN
    FOR j=radius, cols-1-radius DO BEGIN
      ; 获取当前像素周围的像素值
      values = []
      FOR k=-radius, radius DO BEGIN
        FOR l=-radius, radius DO BEGIN
          value = image[i+k, j+l]
          values = [values, value]
        ENDFOR
      ENDFOR
      
      ; 对周围像素值进行排序并取中值
      sorted_values = SORT(values)
      median_value = sorted_values[(kernel_size^2-1)/2]
      
      ; 将中值赋给输出图像的当前像素
      output[i, j] = median_value
    ENDFOR
  ENDFOR
  
  ; 返回输出图像
  RETURN, output
END

其中，`image`是输入的图像，`kernel_size`是中值滤波核的大小。函数中，首先计算中值滤波核的半径，然后遍历输入图像中的每个像素，对其周围的像素值进行排序并取中值，最后将中值赋给输出图像的当前像素。函数返回输出图像。