【OPENCV学习】GDAL读取的数据 转为OPENCV类型 并进行最简单的区域生长
2017-03-03 16:37
573 查看
由于遥感影像通常是TIFF数据,OPENCV读取TIFF会出问题,所以采用了先用GDAL读取数据以及影像信息,再转换为OPENCV的Mat类型,这样结合了GDAL强大的支持多种数据格式的能力,又能方便地调用OPENCV的算法函数。
实验结果图如下图所示:
得到的实验结果图上面有很多点点,是因为原图上就有很多噪声,哈哈
代码中有出处的有:
http://blog.csdn.net/robin__chou/article/details/50071313
但这个博客只写了区域生长的代码,是用于医学影像的
还有Mat、cvMat之间的转换是参考:
http://www.cnblogs.com/Key-Ky/p/4150531.html
#include <iostream> #include <opencv2/core/core.hpp> #include <opencv2/highgui/highgui.hpp> #include <gdal.h> #include <gdal_priv.h> #include <highgui.h> using namespace cv; Mat RegionGrow(Mat src,Point2i pt,int th) { GDALAllRegister(); Point2i ptGrowing;//待生长点位置 int nGrowLable=0;//标记是否生长过 int nSrcValue=0;//生长起点灰度值 int nCurValue=0;//当前生长点灰度值 Mat matDst=Mat::zeros(src.size(),CV_8UC1);//创建一个空白区域,填充为黑色 //生长方向顺序数据 int DIR[8][2]={{-1,-1},{0,-1},{1,-1},{1,0},{1,1},{0,1},{-1,1},{-1,0}}; Vector<Point2i> vcGrowPt;//生长点线 vcGrowPt.push_back(pt);//将生长点压入栈中 matDst.at<uchar>(pt.y,pt.x)=255;//标记生长点 nSrcValue=src.at<uchar>(pt.y,pt.x); while(!vcGrowPt.empty())//生长栈不为空则生长 { pt=vcGrowPt.back();//取出一个生长点 vcGrowPt.pop_back(); //八个方向进行生长 for(int i=0;i<9;++i) { ptGrowing.x=pt.x+DIR[i][0]; ptGrowing.y=pt.y+DIR[i][1]; //检查是否是边缘点 if(ptGrowing.x<0||ptGrowing.y<0||ptGrowing.x>(src.cols-1)||(ptGrowing.y>src.rows-1)) { continue; } nGrowLable=matDst.at<uchar>(ptGrowing.y,ptGrowing.x);//待生长点的灰度值 if(nGrowLable==0)//如果标记点还没有生长 { nCurValue=src.at<uchar>(ptGrowing.y,ptGrowing.x); if(abs(nSrcValue-nCurValue)<th) { matDst.at<uchar>(ptGrowing.y,ptGrowing.x)=255; vcGrowPt.push_back(ptGrowing); } } } } return matDst.clone(); } int main() { //注册驱动 GDALAllRegister(); GDALDataset *pData=NULL; pData=(GDALDataset*)GDALOpen("F:\\ETM.tif",GA_ReadOnly); int XSize=pData->GetRasterXSize(); int YSize=pData->GetRasterYSize(); int XBuff=XSize; int YBuff=YSize; int nBand=pData->GetRasterCount(); float *pBuffer; pBuffer=(float*)CPLMalloc(sizeof(float)*(XBuff)*(YBuff)*(nBand)); pData->RasterIO(GF_Read,0,0,XSize,YSize,pBuffer,XBuff,YBuff,GDT_UInt16,nBand,NULL,1,0,1); CvMat src; cvInitMatHeader(&src,YBuff,XBuff,CV_8UC1,pBuffer); Mat src_mat=Mat(&src,true); Point pt= Point(715,632); int th=70; src_mat=RegionGrow(src_mat,pt,th); imwrite("F:\\test2.bmp",src_mat); waitKey(0); }
实验结果图如下图所示:
得到的实验结果图上面有很多点点,是因为原图上就有很多噪声,哈哈
代码中有出处的有:
http://blog.csdn.net/robin__chou/article/details/50071313
但这个博客只写了区域生长的代码,是用于医学影像的
还有Mat、cvMat之间的转换是参考:
http://www.cnblogs.com/Key-Ky/p/4150531.html
相关文章推荐
- 计算机视觉与图像处理学习笔记(三)opencv的基本数据类型与简单图像处理函数
- Python学习日志 - 2.变量和简单数据类型(习题解答)
- 学习opencv——读取视频并进行canny检测
- 学习opencv——从摄像机读取视频流并进行边缘检测
- opencv学习-2-0初探-基本类型数据
- 【C++】学习笔记二——简单数据类型
- 格式化程序尝试对消息反序列化时引发异常: 尝试对参数 http://tempuri.org/ 进行反序列化时出错: GetLzdtArticleResult。InnerException 消息是“反序列化对象 属于类型 lzdt.DTO.Dtolzdt[] 时出现错误。读取 XML 数据时,超出最大
- BingMap读取shape(arcgis)格式数据(将shape文件转成Sql server2008的geometry类型,再进行读取)
- 基于Qt与OpenCV的图片读取、保存,并利用Canny进行简单的边缘处理结果显示
- C语言学习笔记(二)--从最简单的数据类型开始
- XML学习笔记4——XSD简单数据类型
- OpenCV学习五:图片矩形区域的读取
- opencv学习-2-0初探-基本类型数据
- VS2013MFC对话框工程学习笔记七 - C++数据类型 - 整形 的简单学习
- iOS学习笔记3—NSKeyedUnarchiver进行本地数据的读取|存储
- TensorFlow学习笔记2——数据类型及简单运算
- opencv 2.x学习笔记(三)基本数据类型(一)
- opencv 2.x学习笔记(四)基本数据类型(二)
- 【OpenCV学习笔记】2.1OpenCV基本数据类型
- C++学习笔记——2.简单的变量和数据类型