基于opencv检测斑马鱼重心

// stdafx.h : 标准系统包含文件的包含文件，

// 或是经常使用但不常更改的

// 特定于项目的包含文件

//

#pragma once

#include "targetver.h"

#include <stdio.h>

#include <tchar.h>

// TODO: 在此处引用程序需要的其他头文件

#include<iostream>

#include<string>

#include<opencv2/core/core.hpp>

#include<opencv2/imgproc/imgproc.hpp>

#include<opencv2/highgui/highgui.hpp>

using namespace std;

using namespace cv;

// zebrafish.cpp : 定义控制台应用程序的入口点。

//

#include "stdafx.h"

int _tmain(int argc, _TCHAR* argv[])

{

Mat img1,img2,img3,img4;

img1=imread("E:\\visual studio 2010\\zebrafish\\zebrafish.jpg");//读入RGB图像

imshow("zebrafish_rgb",img1);//显示RGB图像

cvtColor(img1,img2,7);//RGB to GRAY

imshow("zebrafish_gray",img2);//显示灰度图像

threshold(img2,img3,157,255,0);//灰度图像二值化

imshow("zebrafish_bin",img3);//显示二值图像

dilate(img3,img4, Mat(3,3,CV_8U), Point(-1,-1), 1);//做膨胀处理

imshow("zebrafish_dilate",img4);

double area = 0;
//面积

int max = 0;
//最大面积

int count = 0;
//面积符合要求的区域计数

vector<vector<Point> > contours;

vector<Vec4i> hierarchy;

Mat img5=255-img4;

imshow("zebrafish_reve",img5);//图像反转

// 寻找轮廓

findContours( img5, contours, hierarchy, CV_RETR_TREE, CV_CHAIN_APPROX_SIMPLE, Point(0, 0) );

// 绘出轮廓

Mat img6= Mat::zeros( img5.size(), CV_8UC1 );

for( int i = 0; i<contours.size(); i++ )

{

area = contourArea(contours[i]);//依次求每个区域的面积

if (area>max)

{

max=area;

count=i;

}

}

drawContours( img6, contours, count,255,-1,8);

//imshow( "Contours", img6 );

int s=img6.at<unsigned char>(0,0);

long x=0;

long y=0;

int num=0;

//计算连通域的重心

int w=img6.cols;

int h=img6.rows;

int m,n;

for(m=0;m<h;m++)

{

for(n=0;n<w;n++)

{

if(img6.at<unsigned char>(m,n)==255)

{

x=x+m;

y=y+n;

num=num+1;

}

}

}

x=x/num;

y=y/num;

img6.row(x)=0;

img6.col(y)=0;

//在img6上绘出坐标点

Point point=Point(y,x);

char st[4];

int arr[2];

arr[0]=y;

arr[1]=x;

sprintf(st,"(%d,%d)",arr[0],arr[1]);

//String str(st);

putText(img6,st,point,1,1.0,1,1,8,false);

imshow( "Contours", img6 );

waitKey(0);

return 0;

}