Программа фильтрации шумов
    
    Задание.
    Создать программу, осуществляющую фильтрацию шумов в растровых изображениях методами усредняющего, порогового и медианного фильтров.
    Программа написана на языке Object Pascal и выполняется в среде Win32. Общий вид программы показан на рис. № 1.
    
    Рис. № 2. Общий вид главного окна программы c загруженным исходным изображением.
    
    1. Усредняющий фильтp.
    Алгоритм работы усредняющего фильтра заключается в замене значения яркости в
    текущем пикселе на среднюю яркость, вычисленную по его 8 окрестностям, включая и сам элемент. Этот фильтр является самым простым. К недостаткам его можно отнести сглаживание ступенчатых и пилообразных функций. Кроме того пиксели, имеющее существенно отличное значение яркости и являющимися шумовыми могут вносить значительный вклад в результат обработки.
    
    Реализация фильтра представлена в виде процедуры:
    Procedure AverageFilter(Value:Integer);
    
    Данная процедура осуществляет алгоритм усредняющего фильтра применительно к объекту TBitmap. В него предварительно должна быть загружено изображение (Рис. № 2). Результат работы усредняющего фильтра можно увидеть на рис. № 3. Параметр Value - порог при котором производятся манипуляции с пикселом.
    
    Procedure TMainForm.AverageFilter;
    var
    PrevisionLine:pByteArray;
    CurrentLine:pByteArray;
    NextLine:pByteArray;
    I,J:Integer;
    Summ:Integer;
    begin
    if Image1.Picture.Bitmap.PixelFormat = pf8bit then
    begin
    for I := 0 to Image1.Picture.Bitmap.Height - 1 do
    begin
    CurrentLine := Image1.Picture.Bitmap.ScanLine[I];
    for J := 0 to Image1.Picture.Bitmap.Width - 1 do
    begin
    Summ := 0;
    
    if I > 0 then
    begin
    PrevisionLine := Image1.Picture.Bitmap.ScanLine[I - 1];
    if J > 0 then
    begin
    Summ := Summ + PrevisionLine^[J - 1];
    end;
    
    Summ := Summ + PrevisionLine^[J];
    
    if J + 1 < Image1.Picture.Bitmap.Width then
    begin
    Summ := Summ + PrevisionLine^[J + 1];
    end;
    end;
    
    if J > 0 then
    begin
    Summ := Summ + CurrentLine^[J - 1];
    end;
    
    Summ := Summ + CurrentLine^[J];
    
    if J + 1 < Image1.Picture.Bitmap.Width then
    begin
    Summ := Summ + CurrentLine^[J + 1];
    end;
    
    if I + 1 < Image1.Picture.Bitmap.Height then
    begin
    NextLine := Image1.Picture.Bitmap.ScanLine[I + 1];
    if J > 0 then
    begin
    Summ := Summ + NextLine^[J - 1];
    end;
    
    Summ := Summ + NextLine^[J];
    
    if J + 1 < Image1.Picture.Bitmap.Width then
    begin
    Summ := Summ + NextLine^[J + 1];
    end;
    end;
    
    if (Summ div 9) = 0) and (VertB < Image1.Picture.Bitmap.Height) then
    BoxCurrentLine := Image1.Picture.Bitmap.ScanLine[VertB];
    
    for HorB := (Hor - (Value div 2)) to (Hor + (Value div 2)) do
    begin
    
    if (HorB >= 0) and (VertB >= 0) and
    (HorB < Image1.Picture.Bitmap.Width) and
    (VertB < Image1.Picture.Bitmap.Height) then
    PixelArray[Counter] := BoxCurrentLine^[HorB]
    else
    PixelArray[Counter] := 0;
    
    Inc(Counter);
    end;
    end;
    
    // Сортируем массив
    for VertB := 0 to Value*Value - 1 do
    begin
    for HorB := VertB to Value*Value - 1 do
    begin
    if PixelArray[VertB] > PixelArray[HorB] then
    begin
    Temp := PixelArray[VertB];
    PixelArray[VertB] := PixelArray[HorB];
    PixelArray[HorB] := Temp;
    end;
    end;
    end;
    
    // Берем то что посередине и присваиваем текущему пикселю
    CurrentLine^[Hor] := PixelArray[((Value*Value) div 2) + 1];
    end;
    
    end;
    
    Image1.Visible := False;
    Image1.Visible := True;
    N4.Enabled := True;
    end
    else
    MessageBox(Handle,'Такой формат файла пока не поддерживается...',
    'Слабоват я пока...',MB_OK or MB_ICONSTOP or MB_APPLMODAL);
    end;
    end;
    
    Результат работы фильтра можно увидеть на рис. ............