GDAL 2.1.0工具开发之gdal_translate（图像裁剪、缩放、拉伸、赋坐标投影、格式转换等）

最新推荐文章于 2024-04-07 22:14:02 发布
最新推荐文章于 2024-04-07 22:14:02 发布
阅读量1.4w 收藏 8
点赞数
分类专栏： gdal 文章标签： gdal 图像裁剪 格式转换
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zhanghaoping/article/details/51645547
版权

为便于在自己的程序中使用GDAL的示例代码提供的那些强大功能，计划建立一个动态库，在动态库中模拟实现gdal的demo的功能。

初步想法是建立一个类CGdalTools，类的成员即是这些demo的名字。实现与demo一样的功能。

下面的例子是以gdal_translate为例。

找到gdal_translate_bin.cpp（该文件已经相当简单，加上各种注释仅300行出头），复制其中的main函数内容到一个单独的gdal_translate函数中，函数参数与main函数同，返回值与其核心函数GDALTranslate的返回值同，为GDALDatasetH，如下：

GDALDatasetH Gdal_Translate(int argc, char** argv) ;

复制过来的代码不同直接使用，需要做一些修改：

  1. 继续复制GDALTranslateOptionsForBinaryNew和GDALTranslateOptionsForBinaryFree
  2. Gdal驱动的注册改为在CGdalTools类的构建函数中进行，而销毁则放在析构函数中。
  3. 删除EarlySetConfigOptions函数调用。
  4. exit函数全部改为return null;
  5. 认为不必要的地方，可以直接删除。比如输出使用说明、帮助信息、版本信息等地方。
  6. 有些输出Usage函数的地方，改为return NULL;
  7. 函数最后删除GDALClose(hOutDS); return nRetCode;改为return hOutDS;
  8. 找不到GDALTranslateOptionsForBinary类，经检查，该类其实定义于gdal_utils_priv.h中，该定义主要是用来演示DEMO的使用而添加，并不是GDAL的核心功能，因此，GDAL的编译并没有将它导出。为便于开发，最简单的方法是将该文件复制到gdal开发版本的include目录下，然后在所有结构体的定义前加CPL_DLL即可。

修改后的代码已经可以运行了，其输入参数与在命令行中调用gdal_translate完全相同。代码示例如下，需要注意的是输出的GDALDatasetH须在外面关闭，也就是调用在上面的第7步中删除的GDALClose(hOutDS);

GDALDatasetH CGdalTools::Gdal_Translate(int argc, char** argv)
{
    GDALDatasetH    hDataset, hOutDS;
    int bUsageError;

    argc = GDALGeneralCmdLineProcessor(argc, &argv, 0);
    if (argc < 1)
        return NULL;

    if (CPLGetConfigOption("GDAL_MAX_DATASET_POOL_SIZE", NULL) == NULL)
    {
        CPLSetConfigOption("GDAL_MAX_DATASET_POOL_SIZE", "450");
    }

    GDALTranslateOptionsForBinary* psOptionsForBinary = GDALTranslateOptionsForBinaryNew();
    GDALTranslateOptions *psOptions = GDALTranslateOptionsNew(argv + 1, psOptionsForBinary);
    CSLDestroy(argv);

    if (psOptions == NULL)
    {
        return NULL;
    }

    if (psOptionsForBinary->pszSource == NULL)
    {
        return NULL;
    }

    if (psOptionsForBinary->pszDest == NULL)
    {
        return NULL;
    }

    if (strcmp(psOptionsForBinary->pszDest, "/vsistdout/") == 0)
    {
        psOptionsForBinary->bQuiet = TRUE;
    }

    if (!(psOptionsForBinary->bQuiet))
    {
        GDALTranslateOptionsSetProgress(psOptions, GDALTermProgress, NULL);
    }

    if (!psOptionsForBinary->bQuiet && !psOptionsForBinary->bFormatExplicitlySet)
        CheckExtensionConsistency(psOptionsForBinary->pszDest, psOptionsForBinary->pszFormat);

    /* -------------------------------------------------------------------- */
    /*      Attempt to open source file.                                    */
    /* -------------------------------------------------------------------- */

    hDataset = GDALOpenEx(psOptionsForBinary->pszSource, GDAL_OF_RASTER | GDAL_OF_VERBOSE_ERROR, NULL,
        (const char* const*)psOptionsForBinary->papszOpenOptions, NULL);

    if (hDataset == NULL)
    {
        return NULL;
    }

    /* -------------------------------------------------------------------- */
    /*      Handle subdatasets.                                             */
    /* -------------------------------------------------------------------- */
    if (!psOptionsForBinary->bCopySubDatasets
        && CSLCount(GDALGetMetadata(hDataset, "SUBDATASETS")) > 0
        && GDALGetRasterCount(hDataset) == 0)
    {
        fprintf(stderr,
            "Input file contains subdatasets. Please, select one of them for reading.\n");
        GDALClose(hDataset);
        return NULL;
    }

    if (CSLCount(GDALGetMetadata(hDataset, "SUBDATASETS")) > 0
        && psOptionsForBinary->bCopySubDatasets)
    {
        char **papszSubdatasets = GDALGetMetadata(hDataset, "SUBDATASETS");
        char *pszSubDest = (char *)CPLMalloc(strlen(psOptionsForBinary->pszDest) + 32);
        int i;

        CPLString osPath = CPLGetPath(psOptionsForBinary->pszDest);
        CPLString osBasename = CPLGetBasename(psOptionsForBinary->pszDest);
        CPLString osExtension = CPLGetExtension(psOptionsForBinary->pszDest);
        CPLString osTemp;

        const char* pszFormat = NULL;
        if (CSLCount(papszSubdatasets) / 2 < 10)
        {
            pszFormat = "%s_%d";
        }
        else if (CSLCount(papszSubdatasets) / 2 < 100)
        {
            pszFormat = "%s_%002d";
        }
        else
        {
            pszFormat = "%s_%003d";
        }

        const char* pszDest = pszSubDest;

        for (i = 0; papszSubdatasets[i] != NULL; i += 2)
        {
            char* pszSource = CPLStrdup(strstr(papszSubdatasets[i], "=") + 1);
            osTemp = CPLSPrintf(pszFormat, osBasename.c_str(), i / 2 + 1);
            osTemp = CPLFormFilename(osPath, osTemp, osExtension);
            strcpy(pszSubDest, osTemp.c_str());
            hDataset = GDALOpenEx(pszSource, GDAL_OF_RASTER, NULL,
                (const char* const*)psOptionsForBinary->papszOpenOptions, NULL);
            CPLFree(pszSource);
            if (!psOptionsForBinary->bQuiet)
                printf("Input file size is %d, %d\n", GDALGetRasterXSize(hDataset), GDALGetRasterYSize(hDataset));
            hOutDS = GDALTranslate(pszDest, hDataset, psOptions, &bUsageError);
            if (bUsageError == TRUE)
                return NULL;
            if (hOutDS == NULL)
                break;
            GDALClose(hOutDS);
        }

        GDALClose(hDataset);
        GDALTranslateOptionsFree(psOptions);
        GDALTranslateOptionsForBinaryFree(psOptionsForBinary);
        CPLFree(pszSubDest);

        GDALDestroyDriverManager();
        return 0;

    }

    if (!psOptionsForBinary->bQuiet)
        printf("Input file size is %d, %d\n", GDALGetRasterXSize(hDataset), GDALGetRasterYSize(hDataset));

    hOutDS = GDALTranslate(psOptionsForBinary->pszDest, hDataset, psOptions, &bUsageError);
    if (bUsageError == TRUE)
        return NULL;
    int nRetCode = (hOutDS) ? 0 : 1;

    /* Close hOutDS before hDataset for the -f VRT case */
    //GDALClose(hOutDS);
    GDALClose(hDataset);
    GDALTranslateOptionsFree(psOptions);
    GDALTranslateOptionsForBinaryFree(psOptionsForBinary);

    return hOutDS;
}

继续对该函数进行改造。上面的代码中，变化最大的当属第8步，需要变更GDAL的源文件，gdal_utils_priv.h文件。实际上，GDALTranslateOptionsForBinaryNew只是用作GDALTranslateOptionsNew的构建参数，也就是在GDALTranslateOptionsNew中将某些参数输出给GDALTranslateOptionsForBinaryNew，具体内容可以参看GDALTranslateOptionsForBinaryNew的定义，在gdal_utils_priv.h的第52行，具体内容不展开。需要注意的是在GDALTranslateOptionsNew是如何添加GDALTranslateOptionsForBinaryNew信息的。

在gdal_translate_lib.cpp中，第1674行，能看到GDALTranslateOptionsNew函数。里面与GDALTranslateOptionsForBinaryNew相关的内容如下：

    if( EQUAL(papszArgv[i],"-of") && i < argc-1 )
    {
        ++i;
        CPLFree(psOptions->pszFormat);
        psOptions->pszFormat = CPLStrdup(papszArgv[i]);
        if( psOptionsForBinary )
        {
            psOptionsForBinary->bFormatExplicitlySet = TRUE;
        }
    }
    else if( EQUAL(papszArgv[i],"-q") || EQUAL(papszArgv[i],"-quiet") )
    {
        if( psOptionsForBinary )
            psOptionsForBinary->bQuiet = TRUE;
    }
    else if( EQUAL(papszArgv[i],"-sds")  )
    {
        if( psOptionsForBinary )
            psOptionsForBinary->bCopySubDatasets = TRUE;
    }
    else if( EQUAL(papszArgv[i], "-oo") && i+1 < argc )
    {
        i++;
        if( psOptionsForBinary )
        {
            psOptionsForBinary->papszOpenOptions =
                CSLAddString( psOptionsForBinary->papszOpenOptions,
                                            papszArgv[i] );
        }
    }
    else if( !bGotSourceFilename )
    {
        bGotSourceFilename = true;
        if( psOptionsForBinary )
            psOptionsForBinary->pszSource = CPLStrdup(papszArgv[i]);
    }
    else if( !bGotDestFilename )
    {
        bGotDestFilename = true;
        if( psOptionsForBinary )
            psOptionsForBinary->pszDest = CPLStrdup(papszArgv[i]);
    }
    if( psOptionsForBinary )
    {
        psOptionsForBinary->pszFormat = CPLStrdup(psOptions->pszFormat);
    }

回到上文，不进行第8步，也就是不变动gdal_utils_priv.h文件，也就是在函数外自行构建并管理GDALTranslateOptionsForBinaryNew里的信息。因此，上面的8步中，第1步也不需要。为保证输入文件与输出文件信息，变更函数声明如下：

GDALDatasetH Gdal_Translate(char ** argv, char* srcFile, char* dstFile);

为保证程序调用时传入GDALDatasetH也可以完成相关工作，这里加一个重载函数

GDALDatasetH Gdal_Translate(char ** argv, GDALDatasetH hDataset, char* dstFile);

前者详细如下：

GDALDatasetH CGdalTools::Gdal_Translate(char** argv, char* srcFile, char* dstFile)
{
    if (srcFile == NULL)
    {
        return NULL;
    }   

    char** papszOpenOptions=nullptr;

    int argc = CSLCount(argv);
    for (int i = 0; i < argc; i++)
    {

        if (EQUAL(argv[i], "-oo") && i + 1 < argc)
        {
            i++;
            papszOpenOptions = CSLAddString(papszOpenOptions, argv[i]);
        }
    }

    GDALDatasetH    hDataset;

    hDataset = GDALOpenEx(srcFile, GDAL_OF_RASTER | GDAL_OF_VERBOSE_ERROR, NULL,
        (const char* const*)papszOpenOptions, NULL);

    if (hDataset == NULL)
    {
        return NULL;
    }

    GDALDatasetH hOutDs = Gdal_Translate(argv, hDataset, dstFile);

    GDALClose(hDataset);

    return hOutDs;
}

后者详细如下：

GDALDatasetH CGdalTools::Gdal_Translate(char** argv, GDALDatasetH hDataset, char* dstFile)
{
    if (hDataset == NULL)
    {
        return NULL;
    }

    if (dstFile == NULL)
    {
        return NULL;
    }

    bool bQuiet = false, bCopySubDatasets = false, bFormatExplicitlySet = false;
    char* pszFormat = NULL;
    char** papszOpenOptions=nullptr;

    int argc = CSLCount(argv);
    for (int i = 0; i < argc; i++)
    {
        if (EQUAL(argv[i], "-of") && i < argc - 1)
        {
            ++i;
            pszFormat = CPLStrdup(argv[i]);
            bFormatExplicitlySet = TRUE;
        }
        else if (EQUAL(argv[i], "-q") || EQUAL(argv[i], "-quiet"))
        {
            bQuiet = TRUE;
        }
        else if (EQUAL(argv[i], "-sds"))
        {
            bCopySubDatasets = TRUE;
        }       
        else if (EQUAL(argv[i], "-oo") && i + 1 < argc)
        {
            i++;
            papszOpenOptions = CSLAddString(papszOpenOptions, argv[i]);
        }
    }

    GDALDatasetH    hOutDS;
    int bUsageError;

    argc = GDALGeneralCmdLineProcessor(argc, &argv, 0);
    if (argc < 1)
        return NULL;

    if (CPLGetConfigOption("GDAL_MAX_DATASET_POOL_SIZE", NULL) == NULL)
    {
        CPLSetConfigOption("GDAL_MAX_DATASET_POOL_SIZE", "450");
    }

    GDALTranslateOptions *psOptions = GDALTranslateOptionsNew(argv + 1, NULL);
    CSLDestroy(argv);

    if (strcmp(dstFile, "/vsistdout/") == 0)
    {
        bQuiet = TRUE;
    }

    if (!bQuiet)
    {
        GDALTranslateOptionsSetProgress(psOptions, GDALTermProgress, NULL);
    }

    if (!bQuiet && !bFormatExplicitlySet)
        CheckExtensionConsistency(dstFile, pszFormat);


    /* -------------------------------------------------------------------- */
    /*      Handle subdatasets.                                             */
    /* -------------------------------------------------------------------- */
    if (!bCopySubDatasets
        && CSLCount(GDALGetMetadata(hDataset, "SUBDATASETS")) > 0
        && GDALGetRasterCount(hDataset) == 0)
    {
        fprintf(stderr,
            "Input file contains subdatasets. Please, select one of them for reading.\n");
        return NULL;
    }

    if (CSLCount(GDALGetMetadata(hDataset, "SUBDATASETS")) > 0
        && bCopySubDatasets)
    {
        char **papszSubdatasets = GDALGetMetadata(hDataset, "SUBDATASETS");
        char *pszSubDest = (char *)CPLMalloc(strlen(dstFile) + 32);
        int i;

        CPLString osPath = CPLGetPath(dstFile);
        CPLString osBasename = CPLGetBasename(dstFile);
        CPLString osExtension = CPLGetExtension(dstFile);
        CPLString osTemp;

        const char* pszFormat = NULL;
        if (CSLCount(papszSubdatasets) / 2 < 10)
        {
            pszFormat = "%s_%d";
        }
        else if (CSLCount(papszSubdatasets) / 2 < 100)
        {
            pszFormat = "%s_%002d";
        }
        else
        {
            pszFormat = "%s_%003d";
        }

        const char* pszDest = pszSubDest;

        for (i = 0; papszSubdatasets[i] != NULL; i += 2)
        {
            char* pszSource = CPLStrdup(strstr(papszSubdatasets[i], "=") + 1);
            osTemp = CPLSPrintf(pszFormat, osBasename.c_str(), i / 2 + 1);
            osTemp = CPLFormFilename(osPath, osTemp, osExtension);
            strcpy(pszSubDest, osTemp.c_str());
            hDataset = GDALOpenEx(pszSource, GDAL_OF_RASTER, NULL,
                (const char* const*)papszOpenOptions, NULL);
            CPLFree(pszSource);
            if (!bQuiet)
                printf("Input file size is %d, %d\n", GDALGetRasterXSize(hDataset), GDALGetRasterYSize(hDataset));
            hOutDS = GDALTranslate(pszDest, hDataset, psOptions, &bUsageError);
            if (bUsageError == TRUE)
                return NULL;
            if (hOutDS == NULL)
                break;
            GDALClose(hOutDS);
        }

        GDALTranslateOptionsFree(psOptions);

        CPLFree(pszSubDest);

        return NULL;

    }

    if (!bQuiet)
        printf("Input file size is %d, %d\n", GDALGetRasterXSize(hDataset), GDALGetRasterYSize(hDataset));

    hOutDS = GDALTranslate(dstFile, hDataset, psOptions, &bUsageError);
    if (bUsageError == TRUE)
    {
        GDALClose(hOutDS);
        hOutDS = NULL;
    }

    GDALTranslateOptionsFree(psOptions);

    return hOutDS;
}

现在编写一个主程序，对上述代码进行测试，

    char** pszInfo = NULL;
    pszInfo = CSLAddString(pszInfo, "change_format");
    pszInfo = CSLAddString(pszInfo, "-of");
    pszInfo = CSLAddString(pszInfo, "JPEG");
    gt.Gdal_Translate( pszInfo, "G:\\sss.tif", "G:\\sss.jpg");

经测试，程序运行界面如下

这里写图片描述

杂货铺-zzz
关注
  • 0
    点赞
  • 踩
  • 8
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
lambda-gdal_translate:在AWS Lambda中运行gdal_translate二进制文件以处理Amazon S3中的大型图像集合
05-15
这是基于Hector Castro的，他在其中展示了如何包装gdalinfo以在AWS Lambda上运行。 有关于AWS LAMBDA运行任意可执行文件的概述。 lambda-gdal_translate 该脚本允许您使用执行环境运行实用程序。 通常，它允许您运行批处理操作，其中的一行可能看起来像这样， gdal_translate -b 1 -b 2 -b 3 -of GTiff -outsize 50P% -co tiled=yes -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 ' -co PHOTOMETRIC=YCBCR -co COMPRESS=JPEG -co JPEG_QUALITY= ' 85 ' input.tif output.tif 但可以通过无服务器方式使用AWS Lambda。 通常的想法是，使用AWS Lambda可以
GMT画图：gdal_translate转换画图所需的“.grd”格式
weixin_42215013的博客
07-12 3262
平常使用GMT画图，如果需要用到栅格数据的话，则需要用“.grd”的数据，早前都是使用matlab将二进制的数据转换成xyz的三列的文本，再用gmt 内置的 xyz2grd进行转换，这样的转换不仅效率低，而且极其占用硬盘空间。 而GDAL工具中的gdal_translate工具则很好地解决了这个问题，只需要将数据转换成tif格式，再利用gdal_translate就可以把栅格数据转换成GMT所需的...
GDAL APP
长歌行
04-05 1590
原文地址：GDAL附带工具作者：xxzxlx 【本文转载于：http://wiki.woodpecker.org.cn/moin/lilin/gdal-utils，版权属于原创作者】 http://blog.163.com/lionyue@126/blog/static/10793071200710120467784/   gdal是一个库，所以它附带的程序本身是很容易
使用GDAL进行简单的坐标系转换
最新发布
GIS探索之路
04-07 380
使用python GDAL进行简单的坐标系转换，暂时不考虑不同基准坐标系转换的精度问题。
gdal应用-gdal_translate
弗里曼的小伙伴
03-31 9834
gdal应用-gdal_translate [email protected] 2015年1月16日   1 gdal_translate:转换工具。包括裁剪，格式转换等。可能会进行重采样、取子集、变换分辨率等。 注意：只能使用nearest重采样方法。 gdal_translate[--help-general]        [-ot{Byte/Int16/UInt16/UInt3
GDAL转换数据格式C++实现（gdal_translate）
Aldenw的博客
08-05 1146
转载至：https://headfirst-gdal.readthedocs.io/en/latest/gdal-tool-code.html#gdalinfo 添加头文件： #include "gdal_priv.h" #include "cpl_string.h" #include "gdal_utils.h" 函数实现 /** * @param pszDest 输出路径 * @param hSrcDataset 输入数据集 * @param psOptionsIn 命令行参数
python gdal 选择通道将tif格式影像转换为png
qq_43509463的博客
11-02 2898
def RGBNir_tiftoRGB_PNG(): options = gdal.TranslateOptions(format='PNG', bandList=[4, 3, 2]) gdal.Translate('RGB.png', 'GF.tif', options=options) 其中bandlist为原tif中用于生成png的波段数组，如缺省该参数，则默认为[1, 2, 3] 以GF-2影像为例，选择波段[4 , 3, 2]，生成的PNG如下 如由其他要求可自行前往源代.
gdal_translate用法
msl201467的博客
12-25 5770
1、geotiff转换为GMT所需的grd文件 gdal_translate -of GMT -a_nodata NAN input.tif output.tif 2、有时候geotiff不止一个波段，需要转成单波段的geotiff才能进一步转化为GMT的grd文件 gdal_translate -b 1 input.tif output.tif 3、更改投影系统，比如把横轴墨卡托投影转换为带有经纬度的WGS84坐标系 gdalwrap -s_srs def -t_srs def input.tif out
python GDAL遥感影像16位转8位——gdal_translate实现
llc的博客
03-09 2031
例如，此选项可用于将16位（0至65535）数据缩放为8位（0到255）。前两个描述输入图像的最小值和最大值，后一个描述输出栅格的相同值。（其中bn是频段号），可以多次使用该选项来应用于特定波段。，则值将重新缩放为8位等效值（0到255）。也被省略，该函数将自动计算输入栅格的值范围。文档中提到，通过使用。
gdal.Translate和gdal.Wrap支持的creationOptions（命令行中的-co参数）
llc的博客
04-07 3079
在使用gdal_translate和gdalwarp时，使用-co参数时有时候不知道支持哪些可选项，接下来研究一下： $ gdal_translate.exe --format GTiff Format Details: Short Name: GTiff Long Name: GeoTIFF Supports: Raster Extensions: tif tiff Mime...
GDAL对图像文件格式的转换
m0_48011081的博客
03-21 4560
我们常常在图像处理过程中遇到不同软件或程序要求输入的图像格式不同（有些程序或软件支持的数据格式不是常用的Tiff，Img等数据格式），因此需要对不同的数据格式相互进行转换。 我这里以GTiff(.tif)数据转换为PCRaster(.map)数据为例。 首先需要安装GDAL，我这里是在Anaconda上直接安装了基于Python的GDAL，可以在下面网站自行下载，https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 例如下面对应的就是Pyth...
GDAL的C++编译和调用
05-03
对GDAL进行编译并且使用的方法，使用GDAL源码版本GDAL1.9.2，目标是将GDAL编译成支持Win64位的VC++应用程序的类库。
gdal裁剪tif_tif遥感图像_GDAL裁剪_遥感影像_GDAL_gdal批量裁剪_
10-01
基于gdal实现批量裁剪tif图像，对于遥感影像同样适用。
gdal.rar_C 遥感_gdal库_gdal开发_读取遥感图像_遥感
07-15
使用GDAL库读取遥感图像，开发环境为Visual Studio
grb_read.zip_gdal tiff_gdal库格式转换_gdal读取文件_grb文件
09-23
程序基于C++开源库GDAL读取了GRB格式文件，并将其转换为tiff格式
GDAL学习笔记.zip_GDAL_模型格式转换_矢量转换
07-14
包括读取、写入、转换、处理各种栅格数据格式（有些特定的格式对一些操作如写入等不支持）。它使用了一个单一的抽象数据模型就支持了大多数的 栅格数据（GIS 对栅格，矢量，3D 数据模型的抽象能力实在令人叹服）。...
gdal_opengl.rar_GDAL_gdal tiff_opengl gdal_tiff_读取tiff图像
09-20
同时使用gdal和opengl实现tiff图像的读取和现实
GDAL常用命令——遥感图像处理
m0_45074715的博客
08-07 1147
GDAL常用命令——遥感图像处理 最近在做一些遥感图像处理的工作，经常会调用一些gdal的命令对遥感影像进行处理，下面我将对一些常用的gdal命令进行简介，以便今后的学习与使用 gdal 栅格图像处理常用命令 1、 gdalinfo 查询栅格数据集信息， 常用参数：-json 以json格式输出数据； 例子：gdalinfo XXX.tif 2、gdal_translate 该命令可以在不同格式之间转换栅格数据，如数据类型转换以及设置无效值等 常用参数： -ot 强制输出图像带具有驱动程序支持的特
PythonGDAL给图像设置投影坐标系
weixin_40501429的博客
03-22 7498
网站上下载的部分遥感数据只有地理坐标，没有投影坐标，本文章尝试给“只有地理坐标的图片”添加投影坐标。 原始数据：山东地区ASTER GDEM30m数据，该数据只有WGS84地理坐标 目标数据：给该数据添加UTM-50N 1.确定你的目标投影坐标信息，就可以通过下述代码实现。 （1）首先你得知道想添加的目标投影信息，可以通过网上查询的方式。 此网站查的比较费劲，我反正没在这找到：https://developers.arcgis.com/javascript/3/jshelp/pcs.htm， 该
为什么python解释器下有GDAL库还会报错gdal_translate'
06-06
即使你已经在Python解释器下成功导入了GDAL库，但是在运行`gdal_translate`命令时仍然会报错，这是因为`gdal_translate`是一个独立的命令行工具，它并不是GDAL库的一部分，而是GDAL提供的一个工具。这个工具通常是在GDAL的安装目录下的/bin目录内，也就是GDAL的可执行文件路径。 在Python中调用`gdal_translate`命令时，Python解释器会尝试在系统路径中查找该命令，而不是在GDAL的可执行文件路径中查找。因此，如果系统路径中没有包含GDAL的可执行文件路径，Python解释器就无法找到`gdal_translate`命令，从而报错。 解决这个问题的方法有两种： 1. 在系统路径中添加GDAL的可执行文件路径 可以将GDAL的可执行文件路径添加到系统路径中，这样Python解释器就可以找到`gdal_translate`命令了。在Windows系统中，可以按以下步骤操作： - 打开“控制面板” -> “系统和安全” -> “系统” -> “高级系统设置” -> “环境变量”。 - 在“系统变量”中找到“Path”变量，点击“编辑”。 - 在变量值的最后添加GDAL的可执行文件路径，例如：“C:\Program Files\GDAL\bin”。 - 确认保存并关闭所有窗口。重新启动Python解释器，再次尝试运行`gdal_translate`命令。 2. 在Python程序中指定GDAL的可执行文件路径 可以在Python程序中手动指定`gdal_translate`命令的路径，例如： ```python import subprocess gdal_translate_path = r'C:\Program Files\GDAL\bin\gdal_translate.exe' subprocess.call([gdal_translate_path, 'input_file.tif', 'output_file.tif']) ``` 这样就可以保证Python程序能够正确地调用`gdal_translate`命令。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值