tiff文件读取

转自：http://blog.youkuaiyun.com/zhouxuguang236/article/details/7846615

以下是VC下读取TIFF文件的代码

[cpp]  view plain copy

1. char* szFileName = "K:\\地图\\fujian-DEM\\fujian1.tif";  
2.     TIFF* tiff = TIFFOpen(szFileName, "r");//打开Tiff文件，得到指针，以后所有的操作都通过指针进行  
3.   
4.     int nTotalFrame = TIFFNumberOfDirectories(tiff);    //得到图像的总帧数  
5.   
6.     //TIFFSetDirectory(tiff,0);  
7.     //我们打开第一幅图，也就是第0帧，如果是第1帧，第二个参数写1，由此类推。因为Windows下图像基本  
8.     //操作都是以BMP格式进行，我们读出该帧并转成BMP格式。  
9.   
10.     char *dtitle;  
11.     TIFFGetField(tiff,TIFFTAG_PAGENAME,&dtitle);  
12.     //得到该帧的名字，存放在dtitle中。  
13.   
14.     int width,height;  
15.     TIFFGetField(tiff, TIFFTAG_IMAGEWIDTH, &width);  //得到宽度  
16.     TIFFGetField(tiff, TIFFTAG_IMAGELENGTH, &height);//得到高度  
17.   
18.     float resolution = max(width,height);  
19.   
20.     uint16 bitspersample = 1;  
21.     uint16 samplesperpixel = 1;  
22.   
23.     TIFFGetField(tiff, TIFFTAG_SAMPLESPERPIXEL, &samplesperpixel);  
24.     //每个像素占多少机器字，24位图samplesperpixel应该等于3。  
25.     TIFFGetField(tiff, TIFFTAG_BITSPERSAMPLE, &bitspersample);  
26.     //每一个机器字长，这里应为8。  
27.   
28.     uint16 bitsperpixel = bitspersample * samplesperpixel;  
29.     //算出每个像素占多少bit，24位图，值为24  
30.     DWORD dwBytePerLine = (width*bitsperpixel+31)/32 *4;  
31.     //由上面几个参数算出图像每行所占字节(BYTE)数。  
32.   
33.   
34.     DWORD64 dwLeng = height*dwBytePerLine;//在内存里存放这帧图像数据所需要的长度  
35.     BYTE* pData = new BYTE[dwLeng];    //为存放数据分配内存空间  
36.   
37.   
38.     uint32* raster;          
39.     uint32 *row;  
40.     raster = (uint32*)malloc(width * height * sizeof (uint32));  
41.     TIFFReadRGBAImage(tiff, width, height, (uint32*)pData, 1);   
42.     //以上几行读出该帧数据,保存到raster中。  
43.   
44.     row = &raster[0];  
45.     LPBYTE bits2 = pData;  
46.     for (int y = 0; y < height; y++)   
47.     {  
48.   
49.         LPBYTE bits = bits2;  
50.         for (int x = 0; x < width; x++)   
51.         {  
52.             *bits++ = (BYTE)TIFFGetB(row[x]);  
53.             *bits++ = (BYTE)TIFFGetG(row[x]);  
54.             *bits++ = (BYTE)TIFFGetR(row[x]);  
55.         }  
56.         row += width;  
57.         bits2 += dwBytePerLine;  
58.     }  
59.     _TIFFfree(raster);  
60.   
61.     //因为Tif的数据存放顺序和Windows下的BMP相反，上面这几句进行转换。  
62.     //转换结束后，数据存在pData里，释放raster所用内存。  
63.   
64.   
65.   
66.     LPBITMAPINFO pInfo = new BITMAPINFO;  
67.     pInfo->bmiHeader.biSize        = sizeof(BITMAPINFOHEADER);  
68.     pInfo->bmiHeader.biWidth        = width;  
69.     pInfo->bmiHeader.biHeight        = width;  
70.     pInfo->bmiHeader.biCompression    = BI_RGB;  
71.   
72.     pInfo->bmiHeader.biClrUsed        = 0;  
73.     pInfo->bmiHeader.biClrImportant    = 0;  
74.     pInfo->bmiHeader.biPlanes        = 1;  
75.     pInfo->bmiHeader.biBitCount = 24;  
76.     pInfo->bmiHeader.biSizeImage        = dwLeng;  
77.   
78.     float xres,yres;  
79.     uint16 res_unit;   
80.     //解析度单位：如是英寸，厘米  
81.     TIFFGetFieldDefaulted(tiff, TIFFTAG_RESOLUTIONUNIT, &res_unit);  
82.   
83.     if(TIFFGetField(tiff, TIFFTAG_XRESOLUTION, &xres) == 0)  
84.     {  
85.         pInfo->bmiHeader.biXPelsPerMeter = 0;  
86.     }  
87.     else  
88.     {  
89.         if(res_unit == 2)    //英寸  
90.         {  
91.             pInfo->bmiHeader.biXPelsPerMeter = xres * 10000 / 254;  
92.         }  
93.         else if(res_unit == 3)    //厘米  
94.         {  
95.             pInfo->bmiHeader.biXPelsPerMeter = xres * 100;  
96.         }  
97.         else  
98.         {  
99.             pInfo->bmiHeader.biXPelsPerMeter = 0;  
100.         }  
101.     }  
102.     //得到该帧TIFF横向解析度，并计算出m_pInfo->bmiHeader.biXPelsPerMeter  
103.   
104.     if(TIFFGetField(tiff, TIFFTAG_YRESOLUTION, &yres) == 0)  
105.     {  
106.         pInfo->bmiHeader.biYPelsPerMeter = 0;  
107.     }  
108.     else  
109.     {  
110.         if(res_unit == 2)    //英寸  
111.         {  
112.             pInfo->bmiHeader.biYPelsPerMeter = yres * 10000 / 254;  
113.         }  
114.         else if(res_unit == 3)    //厘米  
115.         {  
116.             pInfo->bmiHeader.biYPelsPerMeter = yres * 100;  
117.         }  
118.         else  
119.         {  
120.             pInfo->bmiHeader.biYPelsPerMeter = 0;  
121.         }  
122.     }  
123.     //得到该帧TIFF纵向解析度，并计算出m_pInfo->bmiHeader.biYPelsPerMeter  
124.   
125.   
126.     BITMAPFILEHEADER bmheader;  
127.     bmheader.bfType=0x4d42;  
128.     bmheader.bfSize=0;  
129.     bmheader.bfReserved1=0;  
130.     bmheader.bfReserved2=0;  
131.     bmheader.bfOffBits=54;  
132.     //这几句是生成bmp文件的头结构  
133.   
134.     CFile bmpFile;  
135.     bmpFile.Open(_T("c://test.bmp"),CFile::modeCreate|CFile::modeWrite);  
136.     bmpFile.Write(&bmheader,sizeof(BITMAPFILEHEADER));  
137.     bmpFile.Write(&(pInfo->bmiHeader),sizeof(BITMAPINFOHEADER));  
138.     bmpFile.Write(pData,dwLeng);  
139.     bmpFile.Close();  
140.   
141.     //这里，把该帧TIFF保存到了C盘的test.bmp中，可以用看图软件打开浏览一下。  
142.   
143.     //记得释放内存空间  
144.     delete pInfo;  
145.     pInfo = NULL;  
146.     delete pData;  
147.     pData = NULL;  
148.     //如果想直接显示，就不需要释放，调用StretchDIBits在客户区的DC上就可以显示了。  
149.   
150.   
151.     //如果再打开其他帧的话，从TIFFSetDirectory开始循环运行，比如取下一帧就是  
152.     TIFFSetDirectory(tiff,1);  
153.     //记得保存时另换一个bmp文件名。  
154.     //最后，对这个TIFF文件全部操作结束，记得调用  
155.     TIFFClose(tiff);  

下面的代码是用GDAL打开的

[cpp]  view plain copy

1. char* szFileName = "K:\\地图\\fujian-DEM\\fujian1.tif";  
2.     GDALDataset *poDataset;   //GDAL数据集  
3.     GDALAllRegister();  
4.   
5.     poDataset = (GDALDataset*)GDALOpen(szFileName,GA_ReadOnly);  
6.     if( poDataset == NULL )  
7.     {  
8.         AfxMessageBox(_T("文件打开失败!!!"));  
9.         return;  
10.     }   
11.   
12.     GDALRasterBand *poBand;   //遥感的一个波段  
13.     int nBandCount = poDataset->GetRasterCount();  
14.     poBand = poDataset->GetRasterBand(1);   //和数组下标有点不同  
15.   
16.     //获得图像显示窗口的尺寸  
17.     GetClientRect(&m_ViewRect);  
18.   
19.     int nImgSizeX = poDataset->GetRasterXSize();  
20.     int nImgSizeY = poDataset->GetRasterYSize();  
21.   
22.     double adfGeoTransform[6];  
23.     poDataset->GetGeoTransform( adfGeoTransform );  
24.   
25.     double right  = adfGeoTransform[0] + nImgSizeX*adfGeoTransform[1];  
26.     double bottom  = adfGeoTransform[3] + nImgSizeY*adfGeoTransform[5];  
27.   
28.     int nBufferSizeX,nBufferSizeY;  
29.   
30.     nBufferSizeX = nImgSizeX;  
31.     nBufferSizeY = nImgSizeY;  
32.   
33.     int nScrrenWidth = m_ViewRect.Width();  
34.     int nScrrenHeight= m_ViewRect.Height();  
35.   
36.     BYTE *pafScanblock1,*TempLock1;  
37.     pafScanblock1 = (BYTE *) CPLMalloc((nScrrenWidth)*(nScrrenHeight));  
38.     TempLock1 = pafScanblock1;  
39.   
40.     poBand->RasterIO( GF_Read, 0, 0,nBufferSizeX,nBufferSizeY,   
41.         pafScanblock1,nScrrenWidth,nScrrenHeight, GDT_Byte,0, 0 );  
42.   
43.     //在View逐点显示图像  
44.     DWORD dwBytes = (nScrrenWidth * 24) / 8;  
45.     while(((DWORD) dwBytes) % 4)   
46.     {  
47.         dwBytes++;  
48.     }  
49.   
50.     BYTE *szBuffer = new BYTE[nScrrenHeight*dwBytes];  
51.     memset(szBuffer,0,nScrrenHeight*dwBytes);  
52.     BYTE *pTemp = szBuffer;  
53.     CClientDC dc(this);  
54.     int nIndex = 0;  
55.     for (int i=0;i<nScrrenHeight;i++)  
56.     {  
57.         for (int j=0;j<nScrrenWidth;j++)  
58.         {  
59.   
60.             BYTE dn1 = *pafScanblock1;  
61.   
62.             memcpy(szBuffer,(char*)(&dn1),1);  
63.             szBuffer += 1;  
64.   
65.             pafScanblock1 ++;  
66.   
67.         }  
68.   
69.         szBuffer = pTemp+dwBytes*i;  
70.   
71.     }  
72.     CPLFree(TempLock1);  
73.   
74.     BITMAPINFOHEADER bmiHdr;  
75.     BITMAPINFO MapInfo;  
76.     memset(&bmiHdr, 0, sizeof(BITMAPINFOHEADER));  
77.     bmiHdr.biBitCount = 3*8;  
78.     bmiHdr.biClrImportant = 0;  
79.     bmiHdr.biClrUsed = 0;  
80.     bmiHdr.biCompression = BI_RGB;  
81.     bmiHdr.biHeight = -nScrrenHeight;  
82.     bmiHdr.biPlanes = 1;  
83.     bmiHdr.biSize = sizeof(BITMAPINFOHEADER);  
84.     bmiHdr.biSizeImage = 0;  
85.     bmiHdr.biWidth = nScrrenWidth;  
86.   
87.     bmiHdr.biXPelsPerMeter = 0;  
88.     bmiHdr.biYPelsPerMeter = 0;  
89.   
90.     MapInfo.bmiHeader = bmiHdr;  
91.     MapInfo.bmiColors[0].rgbBlue = 0;  
92.     MapInfo.bmiColors[0].rgbGreen = 0;  
93.     MapInfo.bmiColors[0].rgbRed = 0;  
94.     MapInfo.bmiColors[0].rgbReserved = 0;  
95.   
96.     dc.SetStretchBltMode(MAXSTRETCHBLTMODE);  
97.     ::StretchDIBits(dc.GetSafeHdc(), 0, 0, nScrrenWidth, nScrrenHeight,  
98.         0, 0, bmiHdr.biWidth, -bmiHdr.biHeight,  
99.         pTemp, (LPBITMAPINFO)(&MapInfo), DIB_RGB_COLORS, SRCCOPY);  
100.   
101.     GDALClose(poDataset);  
102.     delete []pTemp;