2

我想将一个HDF从UTM（WGS84）重新投影到正弦（WGS84），所以我尝试用GDALAutoCreateWarpedVRT它来完成它。代码如下：

hSrcDS = (GDALDataset*)GDALOpen("HJ1ACCD1.hdf", GA_ReadOnly);
const char *pszSrcWKT = NULL;
char*  pszDstWKT = NULL;
//pszSrcWKT = ProjectionStr;
pszSrcWKT=GDALGetProjectionRef(hSrcDS);
CPLAssert( pszSrcWKT != NULL &&strlen(pszSrcWKT) > 0 );

OGRSpatialReference oSRS;
oSRS.SetSinusoidal(0,0,0);
oSRS.SetWellKnownGeogCS("WGS84");
oSRS.exportToWkt(&pszDstWKT );

GDALWarpOptions*psWarpOptions = GDALCreateWarpOptions();  
psWarpOptions->dfWarpMemoryLimit=500*1024*1024;

hDstDS=(GDALDataset*)(GDALDataset*)GDALAutoCreateWarpedVRT(hSrcDS,pszSrcWKT,pszDstWKT,GRA_Bilinear ,20,psWarpOptions);
GDALDriver *poDriverTiff;
poDriverTiff=GetGDALDriverManager()->GetDriverByName("GTIFF");
poDriverTiff->CreateCopy("d:\\toto.tif",(GDALDataset*)hDstDS,false,NULL,NULL,NULL);

当我设置oSRS.SetSinusoidal(0,0,0)时，结果看起来不错，但分辨率翻了一番（从 30 到 60）。太奇怪了。

4

1 回答 1

1

根据API 文档GDALAutoCreateWarpedVRT：

GDALSuggestedWarpOutput ()函数用于确定输出虚拟文件的边界和分辨率，该文件应该足够大以包含所有输入图像

还有一个GDALSuggestedWarpOutput2() 函数可以帮助建议一组类似要求的输出文件大小。

于 2013-10-06T20:57:56.543 回答