在網上查了很多資料，但是有關同心圓提取切割的項目很少，我把我的思路寫下來，還沒有進行完善，因爲要一步一步得做。
剛開始做先將複雜問題進行簡化，假設我們知道同心圓的座標位置，也就是這個同心圓位於整幅圖像的正中央。那麼我們要做的：
第一步，識別兩個圓在圖像中的位置，它們的圓心相同，只是半徑不同。然後先把大圓切下來，存入image1中，再把原圖中的小圓切下來，存入image2中。
第二步，現在我們得到了兩個新圖像，image1和image2，分別是在原圖基礎上挖去了一個大圓和一個小圓。下面我們用opencv3自帶的subtract函數，對兩幅圖像進行相減，就能得到圓環也就是同心圓了。
如圖：
第三步，進行極座標變換，對於虹膜識別這個項目來說，抽象成數學模型就是將圓環展開成矩形。是通過極座標的逆變換來實現的。但是這個我轉換之後還是灰度圖像，需要進一步進行處理。
下面是代碼（代碼借鑑了網上很多人的例子，我把它根據具體問題處理了一下）：
說明：這個代碼只是初步代碼，還有好多地方需要根據具體情況進行修改。先把這部分工作記錄下來，算是個思路。
#include <opencv2/opencv.hpp>
#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
#include<opencv2/imgproc/imgproc.hpp>
#include<iostream>
using namespace cv;
using namespace std;
bool polar_to_cartesian(cv::Mat& mat_p, cv::Mat& mat_c, double rows_c, double cols_c);
int main(int argc,char* argv[])
{
Mat image = imread("/Users/oumoemoe/Downloads/timg-3.jpeg");
Mat image1(image.rows, image.cols, image.type(), Scalar(180, 120, 50));//, Scalar(180, 120, 50)
Mat image2(image.rows, image.cols, image.type(), Scalar(180, 120, 50));
Mat image3(image.rows, image.cols, image.type(), Scalar(180, 120, 50));
Mat image4(image.rows, image.cols, image.type(), Scalar(180, 120, 50));
Point center(image.cols/2,image.rows/2);
float radius = image.rows/2;
float radius1=image.rows/4;
circle(image, center, radius, Scalar(0, 200, 100), 2, 8, 0);
for (int x = 0; x < image.cols; x++)
{
for (int y = 0; y < image.rows; y++)
{
int temp = ((x - center.x) * (x - center.x) + (y - center.y) *(y - center.y));
if (temp < (radius * radius))
{
image1.at<Vec3b>(Point(x, y))[0] = image.at<Vec3b>(Point(x, y))[0];
image1.at<Vec3b>(Point(x, y))[1] = image.at<Vec3b>(Point(x, y))[1];
image1.at<Vec3b>(Point(x, y))[2] = image.at<Vec3b>(Point(x, y))[2];
}
else
{
image2.at<Vec3b>(Point(x, y))[0] = image.at<Vec3b>(Point(x, y))[0];
image2.at<Vec3b>(Point(x, y))[1] = image.at<Vec3b>(Point(x, y))[1];
image2.at<Vec3b>(Point(x, y))[2] = image.at<Vec3b>(Point(x, y))[2];
}
}
}
for (int x = 0; x < image.cols; x++)
{
for (int y = 0; y < image.rows; y++)
{
int temp = ((x - center.x) * (x - center.x) + (y - center.y) *(y - center.y));
if (temp < (radius1 * radius1))
{
image3.at<Vec3b>(Point(x, y))[0] = image.at<Vec3b>(Point(x, y))[0];
image3.at<Vec3b>(Point(x, y))[1] = image.at<Vec3b>(Point(x, y))[1];
image3.at<Vec3b>(Point(x, y))[2] = image.at<Vec3b>(Point(x, y))[2];
}
else
{
image4.at<Vec3b>(Point(x, y))[0] = image.at<Vec3b>(Point(x, y))[0];
image4.at<Vec3b>(Point(x, y))[1] = image.at<Vec3b>(Point(x, y))[1];
image4.at<Vec3b>(Point(x, y))[2] = image.at<Vec3b>(Point(x, y))[2];
}
}
}
Mat dstImage;
subtract(image2, image4, dstImage,Mat(),-1);
// imshow("image1", image1);
imshow("image2", image2);
// imshow("image3", image3);
imshow("image4", image4);
imshow("dstImage", dstImage);
Mat dstImage2;
polar_to_cartesian(dstImage, dstImage2, 2*radius-2*radius1, 2*3.14*radius);
imshow("final", dstImage2);
while (uchar(waitKey() != 'q')) {}
return 0;
}
bool polar_to_cartesian(cv::Mat& mat_p, cv::Mat& mat_c, double rows_c, double cols_c)//極座標轉換成普通座標，也就是圓轉方
{
mat_c = cv::Mat::zeros(rows_c, cols_c, CV_8UC1);
int polar_d = mat_p.cols;
double polar_r = polar_d / 2.0; // 圓圖半徑
double delta_r = polar_r / rows_c; //半徑因子
double delta_t = 2.0*3.1415 / cols_c; //角度因子
double center_polar_x = (polar_d - 1) / 2.0;
double center_polar_y = (polar_d - 1) / 2.0;
for (int i = 0; i < cols_c; i++)
{
double theta_p = i * delta_t; //方圖第i列在圓圖對應線的角度
double sin_theta = std::sin(theta_p);
double cos_theta = std::cos(theta_p);
for (int j = 0; j < rows_c; j++)
{
double temp_r = j * delta_r; //方圖第j行在圓圖上對應的半徑長度
int polar_x = (int)(center_polar_x + temp_r * cos_theta);
int polar_y = (int)(center_polar_y - temp_r * sin_theta);
mat_c.ptr<uchar>(j)[i] = mat_p.ptr<uchar>(polar_y)[polar_x];
}
}
return true;
}