GDAL矢量栅格化

CString COGRVToR::Run()<br /> {<br /> GDALAllRegister();<br /> OGRRegisterAll();<br /> CPLSetErrorHandler(CPLLoggingHandler);//在刚开始的GDAL调用时<br /> CString strErr;<br /> OGRDataSource *poDS;<br /> poDS = OGRSFDriverRegistrar::Open(m_strInputV, FALSE );<br /> if( poDS == NULL )<br /> {<br /> strErr.Format(“Open %s failed./n”,m_strInputV);<br /> return strErr;<br /> // exit( 1 );<br /> }<br /> OGRLayer* poLayer;<br /> poLayer = poDS->GetLayer(0);<br /> const char * pszOutFileName =m_strOutputR;<br /> char * pszProjection;<br /> OGRSpatialReference * poSRS = poLayer->GetSpatialRef();<br /> OGREnvelope psExent;<br /> poLayer->GetExtent(&psExent);<br /> int m_nImageWidth =1024;// podataset->GetRasterXSize();<br /> int m_nImageHeight =1024;// podataset->GetRasterYSize();<br /> if ( !poSRS )<br /> {<br /> strErr.Format(“无法打开%s地理信息”,m_strInputV);<br /> m_nImageWidth = psExent.MaxX;<br /> m_nImageHeight = psExent.MaxY;<br /> CString str;<br /> str.Format(“%d ,%d” , m_nImageWidth,m_nImageHeight);<br /> str = str+strErr;<br /> AfxMessageBox(str);<br /> // AfxMessageBox(strErr);<br /> // exit(1);<br /> // return strErr;<br /> }<br /> else<br /> {<br /> poSRS->exportToWkt( &pszProjection );</p> <p> }</p> <p> const char *pszFormat = “GTiff”;<br /> GDALDriver *poDriver;<br /> poDriver = GetGDALDriverManager()->GetDriverByName(pszFormat);<br /> // if( poDriver == NULL )<br /> // exit( 1 );<br /> GDALDataset *poNewDS = poDriver->Create( pszOutFileName, m_nImageWidth,m_nImageHeight,1, GDT_Float32 , NULL );<br /> if ( !poSRS )<br /> {<br /> // strErr.Format(“因为无地理信息,%s无法栅格化”,m_strInputV);</p> <p> // AfxMessageBox(strErr);<br /> // return strErr;<br /> //fot jy<br /> double adfGeoTransform[6];<br /> adfGeoTransform[0] = psExent.MinX;<br /> adfGeoTransform[1] = (psExent.MaxX – psExent.MinX)/m_nImageWidth;<br /> adfGeoTransform[2] = 0;<br /> adfGeoTransform[3] = psExent.MaxY;<br /> adfGeoTransform[4] = 0;<br /> adfGeoTransform[5] = (psExent.MinY – psExent.MaxY)/m_nImageHeight;<br /> GDALSetGeoTransform( poNewDS, adfGeoTransform );<br /> // OGRSpatialReference * poOGRSR;<br /> // poOGRSR->SetUTM(50);<br /> // poSRS->exportToWkt(&pszProjection);<br /> // poNewDS->SetProjection(pszProjection);<br /> //for jy</p> <p> }<br /> else<br /> {<br /> double adfGeoTransform[6];<br /> adfGeoTransform[0] = psExent.MinX;<br /> adfGeoTransform[1] = (psExent.MaxX – psExent.MinX)/m_nImageWidth;<br /> adfGeoTransform[2] = 0;<br /> adfGeoTransform[3] = psExent.MaxY;<br /> adfGeoTransform[4] = 0;<br /> adfGeoTransform[5] = (psExent.MinY – psExent.MaxY)/m_nImageHeight;<br /> GDALSetGeoTransform( poNewDS, adfGeoTransform );<br /> poNewDS->SetProjection(pszProjection); </p> <p> }<br /> // GDALDatasetH poDstDs;<br /> // poDstDs = (GDALDatasetH)poNewDS;//GDALOpen(pszOutFileName,GA_Update);<br /> int * pnbandlist;<br /> pnbandlist = new int [1];<br /> pnbandlist[0]=1;<br /> double *dburnValues;<br /> dburnValues = new double[1];<br /> dburnValues[0]=255;<br /> OGRLayerH * player;<br /> player = new OGRLayerH[1];<br /> player[0] = (OGRLayerH)poLayer;<br /> // poLayer->GetInfo()<br /> char **papszOptions = NULL;<br /> // papszOptions = CSLSetNameValue( papszOptions, “ALL_TOUCHED”, “TRUE” );<br /> // papszOptions = CSLSetNameValue( papszOptions, “ATTRIBUTE”, “Band1Mean” );<br /> papszOptions = CSLSetNameValue( papszOptions, “CHUNKSIZE”, “1” );<br /> papszOptions = CSLSetNameValue( papszOptions, “ATTRIBUTE”, m_strPara ); </p> <p>// AfxMessageBox(“start”);<br /> // CPLErr err= GDALRasterizeLayers(poDstDs,1,pnbandlist,1,player,NULL,NULL,dburnValues,papszOptions,NULL,NULL);<br /> // void * pTransformArg;<br /> // pTransformArg = GDALCreateApproxTransformer();<br /> void * pTransformArg=NULL;<br /> void * m_hGenTransformArg=NULL;<br /> m_hGenTransformArg = GDALCreateGenImgProjTransformer( NULL,<br /> pszProjection,<br /> (GDALDatasetH)poNewDS,<br /> poNewDS->GetProjectionRef(),<br /> FALSE, 1000.0, 3 );<br /> pTransformArg = GDALCreateApproxTransformer( GDALGenImgProjTransform,m_hGenTransformArg, 0.125 );</p> <p> CPLErr err= GDALRasterizeLayers((GDALDatasetH)poNewDS,1,pnbandlist,1,player,GDALGenImgProjTransform,m_hGenTransformArg,NULL,papszOptions,GDALTermProgresscjs,”矢量栅格化”);<br /> // AfxMessageBox(“end”);<br /> // CPLErr err= GDALRasterizeLayers(po<br /> // DstDs,1,pnbandlist,1,player,NULL,NULL,NULL,papszOptions,NULL,NULL);<br /> GDALDestroyGenImgProjTransformer(m_hGenTransformArg);<br /> GDALDestroyApproxTransformer(pTransformArg);<br /> GDALClose(poNewDS);<br /> OGR_DS_Destroy( poDS );<br /> // GDALClose( poDstDs );<br /> delete[]player;<br /> delete[]pnbandlist;<br /> delete[]dburnValues;<br /> strErr = “栅格化成功”;<br /> return strErr;<br /> }

转载自：https://blog.csdn.net/haojknm/article/details/5854317