Opencv系列5_opencv命名方式读取指定路径DCM并批量另存为JPG图像

快来打我* 2022-09-03 03:26 156阅读 0赞

# **实例5：opencv命名方式读取指定路径DCM并批量另存为JPG图像** #

#include <iostream>
    #include <io.h>  // 当中含有_finddata_t
    #include <string>
    #include <vector>
    
    #include <opencv2/core/core.hpp>
    #include <opencv2/highgui/highgui.hpp>
    #include <opencv2/imgproc/imgproc.hpp>
    
    #include <vtkSmartPointer.h>
    #include <vtkImageViewer2.h>
    #include <vtkImageCast.h>
    #include <vtkDICOMImageReader.h>
    #include <vtkRenderWindow.h>
    #include <vtkRenderWindowInteractor.h>
    #include <vtkRenderer.h>
    #include <vtkImageData.h>
    #include <vtkCoordinate.h>
    
    #define NUMBER 72//(CT文件数量为NUMBER-1)
    
    using namespace std;
    using namespace cv;
    
    // 读入一个CT图，返回它的像素矩阵，使用OpenCV的Mat类型返回
    // VTK读取DICOM图像并将像素值转给OpenCV的Mat对象
    void dicomread(string inputFilename, Mat &img,vtkSmartPointer<vtkDICOMImageReader> &reader)
    {
        img.create(512,512,CV_32SC1); 
    
    	vtkSmartPointer<vtkImageCast> imageCast = 
    		vtkSmartPointer<vtkImageCast>::New();
    
        reader->SetFileName(inputFilename.c_str());
    	
        reader->Update();
    
    	imageCast->SetInputConnection(reader->GetOutputPort());
    	imageCast->SetOutputScalarTypeToInt();
    	imageCast->Update();
    
    	// 图像的基本信息
    	int dims[3];
    	reader->GetOutput()->GetDimensions(dims);
    	
    	//图像的像素值
    	for(int k=0;k<dims[2];k++)
    	{
    	    for(int j=0;j<dims[1];j++)
    		{
    		    for(int i=0;i<dims[0];i++)
    			{
    				//转换数据类型，使用imagecast转到double（或float）
    			    int* pixel = 
    					(int*)(imageCast->GetOutput()->GetScalarPointer(i,j,k)); // 第i列第j行的像素值
    				img.at<int>(j,i) = int(*pixel); // 第j行第i列的像素值
    			}
    		}
    	}
    }
    //可视化DICOM图像
    void showdicom(Mat I)
    {
    	double maxx=0,minn=0;
        double *max = &maxx;
    	double *min = &minn;
    	I.convertTo(I,CV_64FC1);
    	minMaxIdx(I,min,max);
    	for(int i=0;i<I.rows;i++)
    	{
    	    for(int j=0;j<I.cols;j++)
    		{
    		    I.at<double>(i,j) = 255*(I.at<double>(i,j)-minn)*1/(maxx-minn);
    		}
    	}
    	
    	minMaxIdx(I,min,max);
    	for(int i=0;i<I.rows;i++)
    	{
    		for(int j=0;j<I.cols;j++)
    			I.at<double>(i,j) = (I.at<double>(i,j)-minn)*1/(maxx-minn);
    	}
    	//cout << I <<endl;  //I为图像的像素值
    	imshow("DICOM Image",I);
    	waitKey(0);
    }
    
    int main()
    {
    	vector<Mat> images;
    	char ad[128] = { 0 };
    	for (int j = 1;j < NUMBER;j+=10) {
    		string filename = format("F:\\Software\\VTK8.2.0\\vtk1\\helloVtk\\CT\\1066244_20110617_CT_102_215_%03d.dcm", j);
    		Mat I1, G1;//I1为原图，G1为转化图
    		vtkSmartPointer<vtkDICOMImageReader> reader =
    			vtkSmartPointer<vtkDICOMImageReader>::New();
    		// 读入dicom图,并将像素值传递给Mat对象I1
    		dicomread(filename, I1, reader);
    		flip(I1, I1, 0);
    		//显示得到的Mat对象I1的信息*(单通道 大小512*512)
    		cout << I1.channels() << "  " << I1.size() << endl;
    		showdicom(I1);
    		//保存图像
    		sprintf_s(ad, "F:\\Software\\VTK8.2.0\\vtk1\\helloVtk\\JPG_images\\%d.jpg", j);
    		imwrite(ad, I1);
    	}
    }

每按下空格，进行下一张图像读取和显示：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2][]

程序运行完后，在F:\\Software\\VTK8.2.0\\vtk1\\helloVtk\\JPG\_images路径下便可以查看另存的JPG图像：

![20210421103213452.png][]

**[本例程配套素材见源码整理文章下载（点击进入）][Link 1]**

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70]: https://img-blog.csdnimg.cn/20210420144151660.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1]: https://img-blog.csdnimg.cn/20210420144227393.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2]: /images/20220829/09b40b3d8e264b70b9616306c0850dd2.png
[20210421103213452.png]: /images/20220829/c1a6de8afc554104950804495646ae86.png
[Link 1]: https://blog.csdn.net/qq_32809093/article/details/115689942