using System;using System.Collections.Generic;using ponentModel;using System.Data;using System.Drawing;using System.Linq;using System.Text;using System.Windows.Forms;using System.IO;using System.Drawing.Imaging;namespace 课程作业3_bk20113268{public partial class Form1 : Form{public Form1(){InitializeComponent();}/// <summary>/// 打开图像/// </summary>/// <param name="sender"></param>/// <param name="e"></param>private void button1_Click(object sender, EventArgs e){//设置文件的类型openFileDialog1.Filter ="*.jpg,*.jpeg,*.bmp,*.gif,*.ico,*.png,*.tif,*.wmf|*.jpg;*.jpeg;*.bmp;*.gif;*.ico;*.png;*.tif;*.wmf";if (openFileDialog1.ShowDialog() == DialogResult.OK) //打开文件对话框{//根据文件的路径创建Image对象Image myImage =System.Drawing.Image.FromFile(openFileDialog1.FileName);pictureBox1.Image = myImage; //显示打开的图片pictureBox1.Height = myImage.Height;pictureBox1.Width = myImage.Width;this.button2.Enabled = true;}}/// <summary>/// 为数据表表头添加行号/// </summary>/// <param name="sender"></param>/// <param name="e"></param>private void dataGridView1_RowPostPaint(object sender, DataGridViewRowPostPaintEventArgs e){try{//添加行号SolidBrush v_SolidBrush = newSolidBrush(dataGridView1.RowHeadersDefaultCellStyle.ForeColor);int v_LineNo = 0;v_LineNo = e.RowIndex + 1;string v_Line = v_LineNo.ToString();e.Graphics.DrawString(v_Line, e.InheritedRowStyle.Font,v_SolidBrush, e.RowBounds.Location.X + 15, e.RowBounds.Location.Y + 5);}catch (Exception ex){MessageBox.Show("添加行号时发生错误，错误信息：" + ex.Message, "操作失败");}}/// <summary>/// 提取特征点/// </summary>/// <param name="sender"></param>/// <param name="e"></param>private void button2_Click(object sender, EventArgs e){Image myImage =System.Drawing.Image.FromFile(openFileDialog1.FileName);pictureBox1.Image = myImage;//显示打开的图片toolStripProgressBar1.Visible = true; //进度条可视toolStripProgressBar1.Maximum = 7; //设置进度条最大长度值toolStripProgressBar1.Value = 0; //设置进度条当前值toolStripProgressBar1.Step = 1; //设置进度条步长toolStripProgressBar1.Value += toolStripProgressBar1.Step; //进度条前进int Var_H = pictureBox1.Image.Height; //获取图象的高度int Var_W = pictureBox1.Image.Width; //获取图象的宽度Bitmap Var_bmp = (Bitmap)pictureBox1.Image; //根据图象的大小创建Bitmap对象double[,] huiduzhi = new double[Var_W, Var_H]; //用于存储各点灰度值for (int i = 0; i < Var_W; i++){for (int j = 0; j < Var_H; j++){Color tem_color = Var_bmp.GetPixel(i, j); //获取当前像素的颜色值huiduzhi[i, j] = tem_color.R * 0.299 + tem_color.G * 0.587 +tem_color.B * 0.114; //各点灰度值}}toolStripProgressBar1.Value += toolStripProgressBar1.Step;double[,] xingquzhi = new double[Var_W, Var_H]; //用于存储各点兴趣值for (int i = 2; i < Var_W - 2; i++){for (int j = 2; j < Var_H - 2; j++){double V1 = 0;for (int m = 0; m < 4; m++){V1 = V1 + Math.Pow(huiduzhi[i - 2 + m, j] - huiduzhi[i -1 + m, j], 2); //计算V1方向相邻像素灰度差平方和}double V2 = 0;for (int m = 0; m < 4; m++){V2 = V2 + Math.Pow(huiduzhi[i - 2 + m, j - 2 + m] - huiduzhi[i- 1 + m, j - 1 + m], 2); //计算V2方向相邻像素灰度差平方和}double V3 = 0;for (int m = 0; m < 4; m++){V3 = V3 + Math.Pow(huiduzhi[i, j - 2 + m] - huiduzhi[i, j- 1 + m], 2); //计算V3方向相邻像素灰度差平方和}double V4 = 0;for (int m = 0; m < 4; m++){V4 = V4 + Math.Pow(huiduzhi[i - 2 + m, j + 2 - m] - huiduzhi[i- 1 + m, j + 1 - m], 2); //计算V4方向相邻像素灰度差平方和}xingquzhi[i, j] = Math.Min(Math.Min(Math.Min(V1, V2), V3), V4);//从V1、V2、V3、V4中取最小值作为该点兴趣值}}toolStripProgressBar1.Value += toolStripProgressBar1.Step;double sum=0;for (int i = 0; i < Var_W; i++){for (int j = 0; j < Var_H; j++){sum+=xingquzhi[i,j];}}double pingjunzhi = sum / (Var_W * Var_H);toolStripProgressBar1.Value += toolStripProgressBar1.Step;int houxuan = 0; //统计候选特征点数目int c = 101; //获取准确特征点数目double zeng = 4;while ( c > 100){c = 0;double yuzhi = pingjunzhi * zeng; //设定阈值double[,] jianding = new double[Var_W, Var_H];for (int i = 0; i < Var_W; i++){for (int j = 0; j < Var_H; j++){if (xingquzhi[i, j] <= yuzhi){jianding[i, j] = 0; //选取兴趣值大于阈值的点作为特征候选点，其他点兴趣值归零}else{jianding[i, j] = 1;houxuan++;}}}int[,] tezhengzhi1 = new int[houxuan, 2]; //假定一个数组能容纳所有点皆为特征点的像素坐标矩阵int yuzhimuban = 5; //定义阈值模板int mubanbanchuang = 2;for (int i = mubanbanchuang; i < Var_W - mubanbanchuang; i = i + yuzhimuban){for (int j = mubanbanchuang; j < Var_H - mubanbanchuang; j = j + yuzhimuban){double MAX = 0; //假定5*5模板最大值起始值为第一个元素值int a = 0; //设a为最大值行int b = 0; //设b为最大值列for (int m = 0; m < yuzhimuban; m++){for (int n = 0; n < yuzhimuban; n++){if (jianding[i - mubanbanchuang + m, j - mubanbanchuang + n] == 1){if (MAX < xingquzhi[i - mubanbanchuang + m, j- mubanbanchuang + n]){MAX = xingquzhi[i - mubanbanchuang + m, j- mubanbanchuang + n]; //获取5*5模板中最大值a = i - mubanbanchuang + m; //获取最大值列b = j - mubanbanchuang + n; //获取最大值行}}else{a = 0;b = 0;}}}if ((a != 0) && (b != 0)){tezhengzhi1[c, 0] = a; //存储特征点列tezhengzhi1[c, 1] = b; //存储特征点行c++; //每有一个既不为0也不重复的最大值特征点数目加一}}}tezhengzhi = tezhengzhi1;zeng += 1;}toolStripProgressBar1.Value += toolStripProgressBar1.Step;toolStripProgressBar1.Value += toolStripProgressBar1.Step;int[,] tezhengdian2 = new int[c, 2]; //定义一个数组存储像素坐标for (int i = 0; i < c; i++){for (int j = 0; j < 2; j++){tezhengdian2[i, j] = tezhengzhi[i,j];}}tezheng = tezhengdian2;Image img = pictureBox1.Image; //将pictureBox1中图像存储入另一个变量Bitmap bmp = new Bitmap(img.Width, img.Height); //创建Bitmap对象Graphics g = Graphics.FromImage(bmp); //创建Graphics对象g.InterpolationMode =System.Drawing.Drawing2D.InterpolationMode.HighQualityBicubic; //设置高质量双三次插值法g.SmoothingMode = System.Drawing.Drawing2D.SmoothingMode.HighQuality; //设置高质量,低速度呈现平滑程度positingQuality =positingQuality.HighQuality; //设置绘制到此 Graphics 的合成图像的呈现质量g.DrawImage(img, 0, 0, img.Width, img.Height); //以img为原本重新于（0,0）点绘制g.Dispose(); //释放资源for (int i = 0; i < c; i++){Graphics myGraphics = Graphics.FromImage(bmp); //创建Graphics对象myGraphics.DrawLine(new Pen(Color.Red, 1), newPoint(tezhengdian2[i, 0], tezhengdian2[i, 1] - 5), new Point(tezhengdian2[i, 0],tezhengdian2[i, 1] + 5)); //画出竖直方向直线myGraphics.DrawLine(new Pen(Color.Red, 1), newPoint(tezhengdian2[i, 0] - 5, tezhengdian2[i, 1]), new Point(tezhengdian2[i, 0] +5, tezhengdian2[i, 1])); //画出水平方向直线myGraphics.Dispose(); //释放资源pictureBox1.Image = bmp; //显示含有“+”的图}toolStripProgressBar1.Value += toolStripProgressBar1.Step;//新建一个datatable用于保存读入的数据DataTable dt = new DataTable();//给datatable添加二个列dt.Columns.Add("列", typeof(String));dt.Columns.Add("行", typeof(String));for (int i = 0; i < c; i++){DataRow dr = dt.NewRow();dr[0] = tezhengdian2[i, 0]; //将列数据赋给表dr[1] = tezhengdian2[i, 1]; //将行数据赋给表dt.Rows.Add(dr); //将这行数据加入到datatable中}this.dataGridView1.DataSource = dt; //将datatable 绑定到datagridview上显示结果dataGridView1.AllowUserToAddRows = false;groupBox1.Text = "特征点像素坐标" + "(" + c.ToString() + "个" + ")"; toolStripProgressBar1.Visible = false; //隐藏进度条MessageBox.Show("共有" + c.ToString() + "个特征点！");}}}。
