using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;

namespace WindowsFormsApplication1
{
    public partial class Form1 : Form
    {
        // アフィン変換行列
        System.Drawing.Drawing2D.Matrix _matAffine;

        public Form1()
        {
            InitializeComponent();
        }

        private void Form1_Load(object sender, EventArgs e)
        {
            // アフィン変換行列
            _matAffine = new System.Drawing.Drawing2D.Matrix();

            // ３０°回転
            _matAffine.Rotate(30, System.Drawing.Drawing2D.MatrixOrder.Append);
            // ２０倍
            _matAffine.Scale(20.0f, 20.0f, System.Drawing.Drawing2D.MatrixOrder.Append);
            // 平行移動
            _matAffine.Translate(-4000f, -2300f, System.Drawing.Drawing2D.MatrixOrder.Append);

            pictureBox1.Image = new Bitmap(pictureBox1.Width, pictureBox1.Height);

            using (var bmp = new Bitmap(@"Mandrill.bmp"))
            using (var g = Graphics.FromImage(pictureBox1.Image))
            {
                g.InterpolationMode = System.Drawing.Drawing2D.InterpolationMode.NearestNeighbor;

                DrawImageLocal(g, _matAffine, bmp);
            }
        }

        private void pictureBox1_MouseMove(object sender, MouseEventArgs e)
        {
            // マウスポインタの座標（コントロールの座標）
            lblSrc.Text = String.Format("(x', y') = ({0}, {1})", e.X, e.Y);

            using (var matInvert = _matAffine.Clone())
            {
                // アフィン変換行列の逆行列を求める
                matInvert.Invert();

                // コントロールの座標
                var points = new PointF[]
                    {
                        new PointF(e.X, e.Y)
                    };

                // 元の座標(画像上の座標）を求める
                matInvert.TransformPoints(points);

                lblDst.Text = String.Format("(x, y) = ({0:#.#}, {1:0.#})", points[0].X, points[0].Y);
            }
        }

        /// <summary>
        /// アフィン変換行列に基づき画像を表示する
        /// </summary>
        /// <param name="g">グラフィックスオブジェクト</param>
        /// <param name="mat">アフィン変換行列</param>
        /// <param name="bmp">表示するBitmapオブジェクト</param>
        private void DrawImageLocal(
            Graphics g,
            System.Drawing.Drawing2D.Matrix mat,
            Bitmap bmp)
        {
            if (bmp == null) return;

                // 描画元の領域
                var srcRect = new RectangleF(-0.5f, -0.5f, bmp.Width, bmp.Height);

                // 描画先の座標(描画元に合わせる、左上、右上、左下の順)
                var points = new PointF[]
                {
                    new PointF(srcRect.Left, srcRect.Top),
                    new PointF(srcRect.Right, srcRect.Top),
                    new PointF(srcRect.Left, srcRect.Bottom),
                };
                // 描画先の座標をアフィン変換で求める（変換後の座標は上書きされる）
                mat.TransformPoints(points);

                // 描画
                g.DrawImage(
                    bmp,
                    points,
                    srcRect,
                    GraphicsUnit.Pixel
                );
        }
    }
}

var bmp = new Bitmap(@"C:\Temp\Mandrill.bmp");

// Bitmapをロック
System.Drawing.Imaging.BitmapData bmpData =
    bmp.LockBits(
        new Rectangle(0, 0, bmp.Width, bmp.Height),
        System.Drawing.Imaging.ImageLockMode.ReadWrite,
        bmp.PixelFormat
    );

// メモリの幅のバイト数を取得
var stride = Math.Abs(bmpData.Stride);
// チャンネル数を取得（8, 24, 32bitのBitmapのみを想定）
var channel = Bitmap.GetPixelFormatSize(bmp.PixelFormat) / 8;

// 元の画像データ用配列
var srcData = new byte[stride * bmpData.Height];
// 処理後の画像データ用配列
var dstData = new byte[stride * bmpData.Height];

// Bitmapデータを配列へコピー
System.Runtime.InteropServices.Marshal.Copy(
    bmpData.Scan0,
    srcData,
    0,
    srcData.Length
    );

int index;

// 移動平均処理
for (int j = 1; j < bmpData.Height - 1; j++)
{
    for (int i = channel; i < (bmpData.Width - 1) * channel; i++)
    {
        index = i + j * stride;
        dstData[index]
            = (byte)((
                srcData[index - channel - stride] + srcData[index - stride] + srcData[index + channel - stride]
              + srcData[index - channel ] + srcData[index ] + srcData[index + channel ]
              + srcData[index - channel + stride] + srcData[index + stride] + srcData[index + channel + stride]
              ) / 9);
    }
}

// 配列をBitmapデータへコピー
System.Runtime.InteropServices.Marshal.Copy(
    dstData,
    0,
    bmpData.Scan0,
    dstData.Length
);

// アンロック
bmp.UnlockBits(bmpData);

bmp.Save(@"C:\Temp\Mandrill_Smooth.bmp", System.Drawing.Imaging.ImageFormat.Bmp);