我有一张地图的图像,我希望在iOS上的内置地图上进行地理参考和叠加.我正在以编程方式使用GDAL.

第1步(确定) – 地理参考(工作正常)

目前我正在计算三个地面控制点的GeoTransform(6个系数),以便对图像进行地理参考,它给出了正确的6个系数.

第2步(问题) – Warp图像获取GeoTransform以获取新的变换图像

图像变得颠倒了！这是因为目标坐标系统(Apple在MapKit中自己的坐标系统)具有反转的y轴,当你向南时,它会增加.

题

如何让GDAL正确地扭曲图像(同时给我一个正确的GeoTransform来配合它)？

我试过的

在变形之前,我在原始GeoTransform中更改了值5/6.这给出了图像的正确扭曲,但新的GeoTransform是错误的.

当前代码

- (WarpResultC*)warpImageWithGeoTransform:(NSArray<NSNumber*>*)geoTransformArray sourceFile:(Nsstring*)inFilepath destinationFile:(Nsstring*)outFilepath
{
    GDALAllRegister();

    GDALDriverH hDriver;
    GDALDataType eDT;
    GDALDatasetH hDstDS;
    GDALDatasetH hSrcDS;

    // Open the source file.
    hSrcDS = GDALOpen( inFilepath.UTF8String,GA_ReadOnly );
    CPLAssert( hSrcDS != NULL );

    // Set the GeoTransform on the source image
    // HERE IS WHERE I NEED NEGATIVE VALUES OF 4 & 5 TO GET A PROPER IMAGE
    double geoTransform[] = { geoTransformArray[0].doubleValue,geoTransformArray[1].doubleValue,geoTransformArray[2].doubleValue,geoTransformArray[3].doubleValue,-geoTransformArray[4].doubleValue,-geoTransformArray[5].doubleValue };
    GDALSetGeoTransform(hSrcDS,geoTransform);

    // Create output with same datatype as first input band.
    eDT = GDALGetRasterDataType(GDALGetRasterBand(hSrcDS,1));

    // Get output driver (GeoTIFF format)
    hDriver = GDALGetDriverByName( "GTiff" );
    CPLAssert( hDriver != NULL );

    // Create a transformer that maps from source pixel/line coordinates
    // to destination georeferenced coordinates (not destination
    // pixel line).  We do that by omitting the destination dataset
    // handle (setting it to NULL).
    void *hTransformArg = GDALCreateGenImgProjTransformer( hSrcDS,NULL,FALSE,1 );
    CPLAssert( hTransformArg != NULL );

    // Get approximate output georeferenced bounds and resolution for file.
    double adfDstGeoTransform[6];
    int nPixels=0,nLines=0;
    CPLErr eErr = GDALSuggestedWarpOutput( hSrcDS,GDALGenImgProjTransform,hTransformArg,adfDstGeoTransform,&nPixels,&nLines );
    CPLAssert( eErr == CE_None );
    GDALDestroyGenImgProjTransformer( hTransformArg );

    // Create the output file.
    hDstDS = GDALCreate( hDriver,outFilepath.UTF8String,nPixels,nLines,4,eDT,NULL );
    CPLAssert( hDstDS != NULL );

    // Write out the projection deFinition.
    GDALSetGeoTransform( hDstDS,adfDstGeoTransform );

    // copy the color table,if required.
    GDALColorTableH hCT = GDALGetRasterColorTable( GDALGetRasterBand(hSrcDS,1) );
    if( hCT != NULL )
        GDALSetRasterColorTable( GDALGetRasterBand(hDstDS,1),hCT );

    // Setup warp options.
    GDALWarpOptions *psWarpOptions = GDALCreateWarpOptions();
    psWarpOptions->hSrcDS = hSrcDS;
    psWarpOptions->hDstDS = hDstDS;

    /* -------------------------------------------------------------------- */
    /*      Do we have a source alpha band?                                 */
    /* -------------------------------------------------------------------- */
    bool enableSrcAlpha = GDALGetRasterColorInterpretation( GDALGetRasterBand(hSrcDS,GDALGetRasterCount(hSrcDS) )) == GCI_AlphaBand;
    if(enableSrcAlpha) { printf( "Using band %d of source image as alpha.\n",GDALGetRasterCount(hSrcDS) ); }

    /* -------------------------------------------------------------------- */
    /*      Setup band mapping.                                             */
    /* -------------------------------------------------------------------- */
    if(enableSrcAlpha)
        psWarpOptions->nBandCount = GDALGetRasterCount(hSrcDS) - 1;
    else
        psWarpOptions->nBandCount = GDALGetRasterCount(hSrcDS);

    psWarpOptions->panSrcBands = (int *) CPLMalloc(psWarpOptions->nBandCount*sizeof(int));
    psWarpOptions->panDstBands = (int *) CPLMalloc(psWarpOptions->nBandCount*sizeof(int));

    for( int i = 0; i < psWarpOptions->nBandCount; i++ )
    {
        psWarpOptions->panSrcBands[i] = i+1;
        psWarpOptions->panDstBands[i] = i+1;
    }

    /* -------------------------------------------------------------------- */
    /*      Setup alpha bands used if any.                                  */
    /* -------------------------------------------------------------------- */
    if( enableSrcAlpha )
        psWarpOptions->nSrcAlphaBand = GDALGetRasterCount(hSrcDS);

    psWarpOptions->nDstAlphaBand = GDALGetRasterCount(hDstDS);

    psWarpOptions->pfnProgress = GDALTermProgress;

    // Establish reprojection transformer.
    psWarpOptions->pTransformerArg = GDALCreateGenImgProjTransformer( hSrcDS,hDstDS,0.0,1 );
    psWarpOptions->pfnTransformer = GDALGenImgProjTransform;

    // Initialize and execute the warp operation.
    GDALWarpOperation oOperation;
    oOperation.Initialize( psWarpOptions );
    CPLErr warpError = oOperation.ChunkAndWarpImage( 0,GDALGetRasterXSize( hDstDS ),GDALGetRasterYSize( hDstDS ) );
    CPLAssert( warpError == CE_None );

    GDALDestroyGenImgProjTransformer( psWarpOptions->pTransformerArg );
    GDALDestroyWarpOptions( psWarpOptions );
    GDALClose( hDstDS );
    GDALClose( hSrcDS );

    WarpResultC* warpResultC = [WarpResultC new];
    warpResultC.geoTransformValues = @[@(adfDstGeoTransform[0]),@(adfDstGeoTransform[1]),@(adfDstGeoTransform[2]),@(adfDstGeoTransform[3]),@(adfDstGeoTransform[4]),@(adfDstGeoTransform[5])];
    warpResultC.newX = nPixels;
    warpResultC.newY = nLines;

    return warpResultC;
}

您可以使用 GDALCreateGenImgProjTransformer执行此操作.从文档：
创建一个从源像素/线坐标映射到目标地理坐标坐标(不是目标像素线)的变换器.我们通过省略目标数据集句柄(将其设置为NULL)来做到这一点.

从gdal_alg.h开始的任务的其他相关信息：
创建图像到图像变换器.

此函数创建一个转换对象,该对象从一个图像上的像素/线坐标映射到另一个图像上的像素/线坐标.图像可能在不同的坐标系中被地理参考,并且可以使用GCP在它们的像素/线坐标和地理参考坐标之间进行映射(与应该使用它们的地理变换的默认假设相反).

此变换器可能执行三个连接转换.

第一阶段是从源图像像素/线坐标到源图像地理参考坐标,并且可以使用地理变换来完成,或者如果没有使用从GCP导出的多项式模型来定义.如果使用GCP,则使用GDALGCPTransform()完成此阶段.

第二阶段是将源坐标系中的投影更改为目标坐标系,假设它们不同.这是使用GDALReprojectionTransform()在内部完成的.

第三阶段是从目标图像地理配准坐标转换为目标图像坐标.这是使用目标图像地理转换完成的,如果不可用,则使用从GCP派生的多项式模型.如果使用GCP,则使用GDALGCPTransform()完成此阶段.如果在创建转换时hDstDS为NULL,则跳过此阶段.

TL; DR：苹果已经在 Java中删除了本机功能支持.特别是触控板手势不再工作.我写了一个使用JNI完成工作的库.看下面我的答案.

我正在寻找几个小时,使Apple Java Extensions在JDK 7或8上工作.我一直在努力尝试所有的想法.互联网上的所有线程关于这个简单地指责JDK.我认为这应该是可以解决的,因为它在JDK 6中有效.为了清楚起见,我可以使其在JDK 6中工作,但是我不希望这个遗留依赖关系,因为不是所有使用我的应用程序的人都有JDK 6或者知道安装吧,让我们说.目标是使用可安装的每个版本的Java(大于JDK 6)的操作系统X功能(菜单栏,手势,处理程序等)运行的双击可用JAR文件.

自从OS X El Capitan的开发者预览版本开始,99％的普通计算机用户就不可能安装Java 6,所以这对我来说已经不再是一个选择.我终于决定编写一个本地库,通过JNI通知Java手势.该项目运行良好,由 GitHub托管.

看看这里：
https://github.com/mcourteaux/MultiTouch-Gestures-Java