[GDAL]2.读取栅格和矢量数据

MFC项目要选择多字节字符集。

GDAL读取栅格：

void DlgImageInfo::OnBnClickedSearch()
{
    // TODO: 在此添加控件通知处理程序代码
    CFileDialog dlg(TRUE,NULL,NULL,OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,_T("IMG影像文件(*.img)|*.img|TIFF影像文件(*.tiff)|*.tiff||"),AfxGetMainWnd());
    CString str;
    if (dlg.DoModal()==IDOK)
    {
        str=dlg.GetPathName();
        m_FileName.SetWindowText(str);

        GDALAllRegister();  
        CPLSetConfigOption("GDAL_FILENAME_IS_UTF8","NO");
        const char* pszFile =(LPCTSTR)str;
        GDALDataset* poDataset;  
        //使用只读方式打开图像   
        poDataset = (GDALDataset*)GDALOpen(pszFile,GA_ReadOnly); 
        if( poDataset == NULL )  
        {  
            CString ex;
            ex.AppendFormat("File: %s不能打开！\n",str);
            listImageInfo.AddString(ex);  //            
        }  
        CString sDriver;
        sDriver.AppendFormat( "Driver: %s/%s\n",
            poDataset->GetDriver()->GetDescription(), 
            poDataset->GetDriver()->GetMetadataItem( GDAL_DMD_LONGNAME ) );

        listImageInfo.AddString(sDriver);         
        
        CString sSize;
        sSize.AppendFormat("Size is %dx%dx%d\n",poDataset->GetRasterXSize(),poDataset->GetRasterYSize(),poDataset->GetRasterCount());
        listImageInfo.AddString(sSize);
        //输出图像的投影信息 
        if( poDataset->GetProjectionRef()  != NULL ) 
        {
            CString sPrj;
            sPrj.Format("Projection is `%s'\n", poDataset->GetProjectionRef());
            listImageInfo.AddString(sPrj)  ; 
        }
        //输出图像的坐标和分辨率信息 
        double adfGeoTransform[6]; 
        if( poDataset->GetGeoTransform( adfGeoTransform ) == CE_None ) 
        { 
            CString sOrgin;
            sOrgin.Format("Origin = (%.6f,%.6f)\n", 
                adfGeoTransform[0], adfGeoTransform[3]);
            CString sPixel;
            sPixel.Format("Pixel Size = (%.6f,%.6f)\n", 
                adfGeoTransform[1], adfGeoTransform[5]);
            listImageInfo.AddString(sOrgin)  ; 
            listImageInfo.AddString(sPixel)  ; 
        }
        GDALClose((GDALDatasetH )poDataset); 
    }
}

 

OGR读取矢量数据：

void VectorInfoDlg::OnBnClickedSearchVector()
{
    CFileDialog dlg(TRUE,NULL,NULL,OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,_T("ESRI Shapefile(*.shp)|*.shp|MapInfo File(*.mid *.mif *.tab)|*.mid *.mif *tab|Arc/Info .E00 (*.e00)|*.e00|AutoCAD DXF(*.dxf)|*.dxf|Comma Separated Value (.csv)|*.csv|GML(*.gml)|*.gml|KML(*.kml)|*.kml|所有文件|*.*||"),AfxGetMainWnd());
    CString strFilePath;
    if (dlg.DoModal()==IDOK)
    {
        strFilePath=dlg.GetPathName();
        m_FileName.SetWindowText(strFilePath);
        const char* pszFile =(LPCTSTR)strFilePath;

        CString strFilename=dlg.GetFileName();
        strFilename.Delete(strFilename.GetLength()-4,4);
        const char* pFileName =(LPCTSTR)strFilename;
        OGRRegisterAll();
        CPLSetConfigOption("GDAL_FILENAME_IS_UTF8","NO");
        OGRDataSource  *poDS;

        poDS = OGRSFDriverRegistrar::Open(pszFile, FALSE );
        if( poDS == NULL )
        {
            AfxMessageBox("Open failed.\n");
            return;
        }

        OGRLayer  *poLayer;

        poLayer = poDS->GetLayerByName( pFileName );
        OGRFeature *poFeature;
        poLayer->ResetReading();
        while( (poFeature = poLayer->GetNextFeature()) != NULL )
        {
            OGRFeatureDefn *poFDefn = poLayer->GetLayerDefn();
            int iField;
            CString strFieldDefn;
            strFieldDefn.AppendFormat("共有%d个属性字段:",poFDefn->GetFieldCount());
            for( iField = 0; iField < poFDefn->GetFieldCount(); iField++ )
            {
                OGRFieldDefn *poFieldDefn = poFDefn->GetFieldDefn( iField );
                
                if( poFieldDefn->GetType() == OFTInteger )
                    strFieldDefn.AppendFormat( "%d,", poFeature->GetFieldAsInteger( iField ) );
                else if( poFieldDefn->GetType() == OFTReal )
                    strFieldDefn.AppendFormat( "%.3f,", poFeature->GetFieldAsDouble(iField) );
                else if( poFieldDefn->GetType() == OFTString )
                    strFieldDefn.AppendFormat( "%s,", poFeature->GetFieldAsString(iField) );
                else
                    strFieldDefn.AppendFormat( "%s,", poFeature->GetFieldAsString(iField) );
            }
            m_VectorInfo.AddString(strFieldDefn);
            OGRGeometry *poGeometry;
            CString strGeometry;
            poGeometry = poFeature->GetGeometryRef();
            if( poGeometry != NULL 
                && wkbFlatten(poGeometry->getGeometryType()) == wkbPoint )
            {
                OGRPoint *poPoint = (OGRPoint *) poGeometry;

                strGeometry.AppendFormat( "%.3f,%3.f\n", poPoint->getX(), poPoint->getY() );
            }
            else
            {
                strGeometry.AppendFormat( "no point geometry\n" );
            }  
            m_VectorInfo.AddString(strFieldDefn);
            OGRFeature::DestroyFeature( poFeature );
        }

        OGRDataSource::DestroyDataSource( poDS );//释放数据源
    }
}

注意：

如果VS编译的是64位程序，那么需要的是64位的GDAL的DLL，不然会报找不到函数依赖项的错误。

#include "..\\GDAL\\include\\gdal_priv.h"
#include "..\\GDAL\\include\\gdal_priv.h"
#include "..\\GDAL\\include\\ogr_spatialref.h"
#include "..\\GDAL\\include\\ogr_attrind.h"
#include "..\\GDAL\\include\\ogrsf_frmts.h"
#include "..\\GDAL\\include\\gdalwarper.h"
#include "..\\GDAL\\include\\gdal.h"
#include "..\\GDAL\\include\\cpl_string.h"
#pragma comment(lib,"..\\GDAL\\lib\\gdal_i.lib")