如何使用GDAL进行AOI裁剪

最新推荐文章于 2025-03-04 17:19:05 发布
最新推荐文章于 2025-03-04 17:19:05 发布
阅读量1.8w 收藏 19
点赞数 6
分类专栏： GIS RS GDAL 遥感 C++编程技术 GDAL交流 文章标签： GDAL AOI 裁切
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/liminlu0314/article/details/6136512
版权

    在工作中，会经常使用一个行政区的矢量图去裁剪一个遥感影像图，得到该行政区的影像图，一般的商业遥感软件都具有这个功能。今天就是用GDAL来实现这一个很实用的功能。首先用到的是GDAL中的gdalwarp，又是warp，呵呵，上一篇就是使用warp进行重采样的。

   首先需要用到gdal源码目录里面的app文件夹下的gdalwarp.cpp文件中的几个函数，大概行数是1651行，直到文件结尾，代码如下：

/************************************************************************/
/*                      GeoTransform_Transformer()                      */
/*                                                                      */
/*      Convert points from georef coordinates to pixel/line based      */
/*      on a geotransform.                                              */
/************************************************************************/

class CutlineTransformer : public OGRCoordinateTransformation
{
public:

    void         *hSrcImageTransformer;

    virtual OGRSpatialReference *GetSourceCS() { return NULL; }
    virtual OGRSpatialReference *GetTargetCS() { return NULL; }

    virtual int Transform( int nCount, 
                           double *x, double *y, double *z = NULL ) {
        int nResult;

        int *pabSuccess = (int *) CPLCalloc(sizeof(int),nCount);
        nResult = TransformEx( nCount, x, y, z, pabSuccess );
        CPLFree( pabSuccess );

        return nResult;
    }

    virtual int TransformEx( int nCount, 
                             double *x, double *y, double *z = NULL,
                             int *pabSuccess = NULL ) {
        return GDALGenImgProjTransform( hSrcImageTransformer, TRUE, 
                                        nCount, x, y, z, pabSuccess );
    }
};


/************************************************************************/
/*                            LoadCutline()                             */
/*                                                                      */
/*      Load blend cutline from OGR datasource.                         */
/************************************************************************/

static void
LoadCutline( const char *pszCutlineDSName, const char *pszCLayer, 
             const char *pszCWHERE, const char *pszCSQL, 
             void **phCutlineRet )

{
#ifndef OGR_ENABLED
    CPLError( CE_Failure, CPLE_AppDefined, 
              "Request to load a cutline failed, this build does not support OGR features./n" );
    exit( 1 );
#else // def OGR_ENABLED
    OGRRegisterAll();

/* -------------------------------------------------------------------- */
/*      Open source vector dataset.                                     */
/* -------------------------------------------------------------------- */
    OGRDataSourceH hSrcDS;

    hSrcDS = OGROpen( pszCutlineDSName, FALSE, NULL );
    if( hSrcDS == NULL )
        exit( 1 );

/* -------------------------------------------------------------------- */
/*      Get the source layer                                            */
/* -------------------------------------------------------------------- */
    OGRLayerH hLayer = NULL;

    if( pszCSQL != NULL )
        hLayer = OGR_DS_ExecuteSQL( hSrcDS, pszCSQL, NULL, NULL ); 
    else if( pszCLayer != NULL )
        hLayer = OGR_DS_GetLayerByName( hSrcDS, pszCLayer );
    else
        hLayer = OGR_DS_GetLayer( hSrcDS, 0 );

    if( hLayer == NULL )
    {
        fprintf( stderr, "Failed to identify source layer from datasource./n" );
        exit( 1 );
    }

/* -------------------------------------------------------------------- */
/*      Apply WHERE clause if there is one.                             */
/* -------------------------------------------------------------------- */
    if( pszCWHERE != NULL )
        OGR_L_SetAttributeFilter( hLayer, pszCWHERE );

/* -------------------------------------------------------------------- */
/*      Collect the geometries from this layer, and build list of       */
/*      burn values.                                                    */
/* -------------------------------------------------------------------- */
    OGRFeatureH hFeat;
    OGRGeometryH hMultiPolygon = OGR_G_CreateGeometry( wkbMultiPolygon );

    OGR_L_ResetReading( hLayer );
    
    while( (hFeat = OGR_L_GetNextFeature( hLayer )) != NULL )
    {
        OGRGeometryH hGeom = OGR_F_GetGeometryRef(hFeat);

        if( hGeom == NULL )
        {
            fprintf( stderr, "ERROR: Cutline feature without a geometry./n" );
            exit( 1 );
        }
        
        OGRwkbGeometryType eType = wkbFlatten(OGR_G_GetGeometryType( hGeom ));

        if( eType == wkbPolygon )
            OGR_G_AddGeometry( hMultiPolygon, hGeom );
        else if( eType == wkbMultiPolygon )
        {
            int iGeom;

            for( iGeom = 0; iGeom < OGR_G_GetGeometryCount( hGeom ); iGeom++ )
            {
                OGR_G_AddGeometry( hMultiPolygon, 
                                   OGR_G_GetGeometryRef(hGeom,iGeom) );
            }
        }
        else
        {
            fprintf( stderr, "ERROR: Cutline not of polygon type./n" );
            exit( 1 );
        }

        OGR_F_Destroy( hFeat );
    }

    if( OGR_G_GetGeometryCount( hMultiPolygon ) == 0 )
    {
        fprintf( stderr, "ERROR: Did not get any cutline features./n" );
        exit( 1 );
    }

/* -------------------------------------------------------------------- */
/*      Ensure the coordinate system gets set on the geometry.          */
/* -------------------------------------------------------------------- */
    OGR_G_AssignSpatialReference(
        hMultiPolygon, OGR_L_GetSpatialRef(hLayer) );

    *phCutlineRet = (void *) hMultiPolygon;

/* -------------------------------------------------------------------- */
/*      Cleanup                                                         */
/* -------------------------------------------------------------------- */
    if( pszCSQL != NULL )
        OGR_DS_ReleaseResultSet( hSrcDS, hLayer );

    OGR_DS_Destroy( hSrcDS );
#endif
}

/************************************************************************/
/*                      TransformCutlineToSource()                      */
/*                                                                      */
/*      Transform cutline from its SRS to source pixel/line coordinates.*/
/************************************************************************/
static void
TransformCutlineToSource( GDALDatasetH hSrcDS, void *hCutline,
                          char ***ppapszWarpOptions, char **papszTO_In )

{
#ifdef OGR_ENABLED
    OGRGeometryH hMultiPolygon = OGR_G_Clone( (OGRGeometryH) hCutline );
    char **papszTO = CSLDuplicate( papszTO_In );

/* -------------------------------------------------------------------- */
/*      Checkout that SRS are the same.                                 */
/* -------------------------------------------------------------------- */
    OGRSpatialReferenceH  hRasterSRS = NULL;
    const char *pszProjection = NULL;

    if( GDALGetProjectionRef( hSrcDS ) != NULL 
        && strlen(GDALGetProjectionRef( hSrcDS )) > 0 )
        pszProjection = GDALGetProjectionRef( hSrcDS );
    else if( GDALGetGCPProjection( hSrcDS ) != NULL )
        pszProjection = GDALGetGCPProjection( hSrcDS );

    if( pszProjection != NULL )
    {
        hRasterSRS = OSRNewSpatialReference(NULL);
        if( OSRImportFromWkt( hRasterSRS, (char **)&pszProjection ) != CE_None )
        {
            OSRDestroySpatialReference(hRasterSRS);
            hRasterSRS = NULL;
        }
    }

    OGRSpatialReferenceH hCutlineSRS = OGR_G_GetSpatialReference( hMultiPolygon );
    if( hRasterSRS != NULL && hCutlineSRS != NULL )
    {
        /* ok, we will reproject */
    }
    else if( hRasterSRS != NULL && hCutlineSRS == NULL )
    {
        fprintf(stderr,
                "Warning : the source raster dataset has a SRS, but the cutline features/n"
                "not.  We assume that the cutline coordinates are expressed in the destination SRS./n"
                "If not, cutline results may be incorrect./n");
    }
    else if( hRasterSRS == NULL && hCutlineSRS != NULL )
    {
        fprintf(stderr,
                "Warning : the input vector layer has a SRS, but the source raster dataset does not./n"
                "Cutline results may be incorrect./n");
    }

    if( hRasterSRS != NULL )
        OSRDestroySpatialReference(hRasterSRS);

/* -------------------------------------------------------------------- */
/*      Extract the cutline SRS WKT.                                    */
/* -------------------------------------------------------------------- */
    if( hCutlineSRS != NULL )
    {
        char *pszCutlineSRS_WKT = NULL;

        OSRExportToWkt( hCutlineSRS, &pszCutlineSRS_WKT );
        papszTO = CSLSetNameValue( papszTO, "DST_SRS", pszCutlineSRS_WKT );
        CPLFree( pszCutlineSRS_WKT );
    }

/* -------------------------------------------------------------------- */
/*      Transform the geometry to pixel/line coordinates.               */
/* -------------------------------------------------------------------- */
    CutlineTransformer oTransformer;

    /* The cutline transformer will *invert* the hSrcImageTransformer */
    /* so it will convert from the cutline SRS to the source pixel/line */
    /* coordinates */
    oTransformer.hSrcImageTransformer = 
        GDALCreateGenImgProjTransformer2( hSrcDS, NULL, papszTO );

    CSLDestroy( papszTO );

    if( oTransformer.hSrcImageTransformer == NULL )
        exit( 1 );

    OGR_G_Transform( hMultiPolygon, 
                     (OGRCoordinateTransformationH) &oTransformer );

    GDALDestroyGenImgProjTransformer( oTransformer.hSrcImageTransformer );

/* -------------------------------------------------------------------- */
/*      Convert aggregate geometry into WKT.                            */
/* -------------------------------------------------------------------- */
    char *pszWKT = NULL;

    OGR_G_ExportToWkt( hMultiPolygon, &pszWKT );
    OGR_G_DestroyGeometry( hMultiPolygon );

    *ppapszWarpOptions = CSLSetNameValue( *ppapszWarpOptions, 
                                          "CUTLINE", pszWKT );
    CPLFree( pszWKT );
#endif
}

    需要将上面的代码稍微进行改造，也就是把exit(1)之类的去掉，换成返回错误代码，具体修改过的代码就不贴了，就那么几个，还有就是记得把函数的返回值void改成int，用于获取错误代码。下面就是我写的函数代码： 

/**
* @brief AOI截图(GDAL)
* @param pszInFile        	输入文件的路径
* @param pszOutFile        	截图后输出文件的路径
* @param pszAOIFile        	AOI文件路径
* @param pszSQL            	指定AOI文件中的属性字段值来裁剪
* @param pBandInex        	指定裁剪的波段，默认为NULL，表示裁剪所有波段
* @param piBandCount    	指定裁剪波段的个数，同上一个参数同时使用
* @param iBuffer        	指定AOI文件外扩范围，默认为0，表示不外扩
* @param pszFormat        	截图后输出文件的格式
* @param pProgress        	进度条指针
* @return 成功返回
*/ 
int CutImageByAOIGDAL(const char* pszInFile, const char* pszOutFile, const char* pszAOIFile, const char* pszSQL, 
    int *pBandInex, int *piBandCount, int iBuffer, const char* pszFormat, LT_Progress *pProgress)
{
    if(pProgress != NULL)
    {
        pProgress->SetProgressCaption("AOI裁剪");
        pProgress->SetProgressTip("开始裁剪图像...");
    }

    GDALAllRegister();
    void *hCutLine = NULL;

    int iRev = LoadCutline( pszAOIFile, "", "", pszSQL, &hCutline );

    GDALDataset * pSrcDS = (GDALDataset*) GDALOpen(pszInFile, GA_ReadOnly);
    if (pSrcDS == NULL)
    {
        if (pProgress != NULL)
            pProgress->SetProgressTip("输入的栅格文件不能打开，请检查文件是否存在！");

        return RE_NOFILE;
    }
    
    int iBandCount = pSrcDS->GetRasterCount();
    const char* pszWkt = pSrcDS->GetProjectionRef();
    GDALDataType dT = pSrcDS->GetRasterBand(1)->GetRasterDataType();

    double adfGeoTransform[6] = {0};
    double newGeoTransform[6] = {0};

    pSrcDS->GetGeoTransform(adfGeoTransform);
    memcpy(newGeoTransform, adfGeoTransform, sizeof(double)*6);

    int *pSrcBand = NULL;
    int iNewBandCount = iBandCount;
    if (pBandInex != NULL && piBandCount != NULL)
    {
        int iMaxBandIndex = pBandInex[0];
        for (int i=1; i<*piBandCount; i++)
        {
            if (iMaxBandIndex < pBandInex[i])
                iMaxBandIndex = pBandInex[i];
        }

        if (iMaxBandIndex > iBandCount)
        {
            if (pProgress != NULL)
                pProgress->SetProgressTip("输入的波段序号没有指定的波段数！");

            GDALClose((GDALDatasetH) pSrcDS);
            return RE_PARAMERROR;
        }

        iNewBandCount = *piBandCount;
        pSrcBand = new int[iNewBandCount];
        for (int i=0; i<iNewBandCount; i++)
            pSrcBand[i] = pBandInex[i];
    }
    else
    {
        pSrcBand = new int[iNewBandCount];
        for (int i=0; i<iNewBandCount; i++)
            pSrcBand[i] = i+1;
    }

    OGRGeometryH hGeometry = (OGRGeometryH) hCutline;
    OGRGeometryH hMultiPoly = NULL;
    if (iBuffer > 0)
    {
        double dDistance = ABS(adfGeoTransform[1]*iBuffer);
        hMultiPoly = OGR_G_Buffer(hGeometry, dDistance, 30);
        OGR_G_AssignSpatialReference(hMultiPoly, OGR_G_GetSpatialReference(hGeometry));
        OGR_G_DestroyGeometry(hGeometry);
    }
    else
        hMultiPoly = hGeometry;

    OGREnvelope eRect;
    OGR_G_GetEnvelope(hMultiPoly, &eRect);

    int iNewWidth = 0, iNewHeigh = 0;
    int iBeginRow = 0, iBeginCol = 0;

    newGeoTransform[0] = eRect.MinX;
    newGeoTransform[3] = eRect.MaxY;

    iNewWidth = static_cast<int>((eRect.MaxX -eRect.MinX) / ABS(adfGeoTransform[1]));
    iNewHeigh = static_cast<int>((eRect.MaxY -eRect.MinY) / ABS(adfGeoTransform[5]));

    Projection2ImageRowCol(adfGeoTransform, newGeoTransform[0], newGeoTransform[3], iBeginCol, iBeginRow);
    ImageRowCol2Projection(adfGeoTransform, iBeginCol, iBeginRow, newGeoTransform[0], newGeoTransform[3]);

    GDALDriver *pDriver = GetGDALDriverManager()->GetDriverByName(pszFormat);
    if (pDriver == NULL)
    {
        if (pProgress != NULL)
            pProgress->SetProgressTip("不能创建指定类型的文件！");

        GDALClose((GDALDatasetH) pSrcDS);
        return RE_NOFILE;
    }

    GDALDataset* pDstDS = pDriver->Create(pszOutFile, iNewWidth, iNewHeigh, iNewBandCount, dT, NULL);
    if (pDstDS == NULL)
    {
        if (pProgress != NULL)
            pProgress->SetProgressTip("创建文件失败！");

        GDALClose((GDALDatasetH) pSrcDS);
        return RE_CREATEFILE;
    }

    pDstDS->SetGeoTransform(newGeoTransform);
    pDstDS->SetProjection(pszWkt);

    void *hTransformArg, *hGenImgProjArg=NULL;
    char **papszTO = NULL;
    /* -------------------------------------------------------------------- */
    /*      Create a transformation object from the source to               */
    /*      destination coordinate system.                                  */
    /* -------------------------------------------------------------------- */
    hTransformArg = hGenImgProjArg = 
        GDALCreateGenImgProjTransformer2( hSrcDS, (GDALDatasetH)pDstDS, papszTO );

    if( hTransformArg == NULL )
    {
        if (pProgress != NULL)
            pProgress->SetProgressTip("创建投影失败，请检查输入参数！");

        GDALClose((GDALDatasetH) pSrcDS);
        GDALClose((GDALDatasetH) pDstDS);

        RELEASE(pSrcBand);
        return RE_PARAMERROR;
    }

    GDALTransformerFunc pfnTransformer = GDALGenImgProjTransform;
    GDALWarpOptions *psWO = GDALCreateWarpOptions();

    psWO->papszWarpOptions = CSLDuplicate(NULL);
    psWO->eWorkingDataType = dT;
    psWO->eResampleAlg = GRA_Bilinear ;

    psWO->hSrcDS = (GDALDatasetH) pSrcDS;
    psWO->hDstDS = (GDALDatasetH) pDstDS;

    psWO->pfnTransformer = pfnTransformer;
    psWO->pTransformerArg = hTransformArg;

    psWO->pfnProgress = GDALProgress;
    psWO->pProgressArg = pProgress;

    psWO->nBandCount = iNewBandCount;
    psWO->panSrcBands = (int *) CPLMalloc(iNewBandCount*sizeof(int));
    psWO->panDstBands = (int *) CPLMalloc(iNewBandCount*sizeof(int));
    for (int i=0; i<iNewBandCount; i++)
    {
        psWO->panSrcBands[i] = pSrcBand[i];
        psWO->panDstBands[i] = i+1;
    }

    RELEASE(pSrcBand);
    
    psWO->hCutline = (void*) hMultiPoly;
    TransformCutlineToSource((GDALDatasetH) pSrcDS, (void*)hMultiPoly, &(psWO->papszWarpOptions), papszTO );


    GDALWarpOperation oWO;
    if (oWO.Initialize(psWO) != CE_None)
    {
        if(pProgress != NULL)
            pProgress->SetProgressTip("转换参数错误！");

        GDALClose((GDALDatasetH) pSrcDS);
        GDALClose((GDALDatasetH) pDstDS);

        return RE_PARAMERROR;
    }

    oWO.ChunkAndWarpImage(0, 0, iNewWidth, iNewHeigh);

    GDALDestroyGenImgProjTransformer(psWO->pTransformerArg);
    GDALDestroyWarpOptions( psWO );
    GDALClose((GDALDatasetH) pSrcDS);
    GDALClose((GDALDatasetH) pDstDS);

    if(pProgress != NULL)
        pProgress->SetProgressTip("图像裁剪完成！");

    return RE_SUCCESS;
}

裁剪后的效果图就不贴了，希望对大家有用。

………………………………………………分割线 下面更新于2013-5-2…………………………………………………………

上面的代码中有两处粘贴错误，已经修改，同时有人问代码中有两个坐标转换的函数，就是行列号和投影坐标相互转换的函数，其实很简单，就是根据六参数计算各仿射变换就OK，下面把代码贴出来。

bool Projection2ImageRowCol(double *adfGeoTransform, double dProjX, double dProjY, int &iCol, int &iRow)
{
	try
	{
		double dTemp = adfGeoTransform[1]*adfGeoTransform[5] - adfGeoTransform[2]*adfGeoTransform[4];
		double dCol = 0.0, dRow = 0.0;
		dCol = (adfGeoTransform[5]*(dProjX - adfGeoTransform[0]) - 
			adfGeoTransform[2]*(dProjY - adfGeoTransform[3])) / dTemp + 0.5;
		dRow = (adfGeoTransform[1]*(dProjY - adfGeoTransform[3]) - 
			adfGeoTransform[4]*(dProjX - adfGeoTransform[0])) / dTemp + 0.5;

		iCol = static_cast<int>(dCol);
		iRow = static_cast<int>(dRow);
		return true;
	}
	catch(...)
	{
		return false;
	}
}

bool ImageRowCol2Projection(double *adfGeoTransform, int iCol, int iRow, double &dProjX, double &dProjY)
{
	//adfGeoTransform[6]  数组adfGeoTransform保存的是仿射变换中的一些参数，分别含义见下
	//adfGeoTransform[0]  左上角x坐标 
	//adfGeoTransform[1]  东西方向分辨率
	//adfGeoTransform[2]  旋转角度, 0表示图像 "北方朝上"
	//adfGeoTransform[3]  左上角y坐标 
	//adfGeoTransform[4]  旋转角度, 0表示图像 "北方朝上"
	//adfGeoTransform[5]  南北方向分辨率

	try
	{
		dProjX = adfGeoTransform[0] + adfGeoTransform[1] * iCol + adfGeoTransform[2] * iRow;
		dProjY = adfGeoTransform[3] + adfGeoTransform[4] * iCol + adfGeoTransform[5] * iRow;
		return true;
	}
	catch(...)
	{
		return false;
	}
}

GDAL+C#实现影像拼接（TIFF格式）
哎呦喂iuv的博客
08-03 4319
目录 任务描述 任务分解 GDAL读写tif文件函数介绍及注意点 具体步骤 ①影像来源 ②确定影像范围 ③写入输入影像数据到输出影像中 任务描述 将多张有重叠部分的影像经过影像拼接，最终得到一张完整的影像，如下： ➕ ═ 任务分解 ①得到两张有重叠部分的影像； ②确定输出影像的范围（长x宽）； ③以此读取输入影像到输出影像中，保证输入影像在输出影像中的正确x、y偏移值（重叠部分默认用后一张影像值替代；无数据部分默认用黑色填充）。 GDAL读写tif文件函数介绍及注意点..
arcgis分隔图层重复出文件_shp文件自相交处理的方法
weixin_39945795的博客
11-27 1860
今天基于GDAL使用shp文件对栅格影像进行裁剪时出现了下面的问题，提示多边形自相交了Warning 1: Ring Self-intersection at or near point 112.48666420300003 34.830899357000078ERROR 1: Cutline polygon is invalid.很多人的第一反应是使用ArcGIS进行拓扑检查，或使用A...
绘制AOI批量裁切图像得到多个图斑
MCM的专栏
09-01 907
把分布不同区域的“足球场”分别裁切出图，绘制AOI批量操作如下所示： 点击图片查看清晰图 MCM智拼图-遥感影像匀色镶嵌分幅处理
python中利用gdal.warp完成矢量裁剪栅格（遥感影像掩膜、shp裁剪tif、图像裁剪）
weixin_44778507的博客
10-25 3125
在python中，实现基于shp矢量对栅格数据的裁剪或者叫掩膜。
gdal.wrap裁剪数据
张欣的专栏
03-04 234
cutlineDSName=shapefile_path, # 用于裁剪的 ShapefilecropToCutline=True, # 根据 cutline 裁剪栅格) # 裁剪区域外的像素设置为 nan。
Python GDAL/OGR Vector layer operation
豆芽菜
09-30 3221
From: https://pcjericks.github.io/py-gdalogr-cookbook/vector_layers.html 目录 Delete a file Is Ogr Installed View Auto Generated Ogr Help Get List of Ogr Drivers Alphabetically (A- Z) Is Ogr Drive...
使用arcgis中遇到的问题及解决办法整合记录【遇到问题、持续更新、有问题也可留言】
weixin_42476786的博客
07-29 7222
使用arcgis遇到的问题整合
Cannot compute bounding box of cutline
04-14 930
最近利用GDAL的gdalwarp来切图，却总是出现 Cannot compute bounding box of cutline
用nc裁剪nc
03-08
### 使用 Python 中的 GDAL 和 xarray 裁剪 NetCDF 文件 #### 方法一：基于 GDAL 的裁剪方法 可以采用一种较为简便的方式，即通过 Python 结合 GDAL 库来实现 SHP 对 NetCDF 数据的裁剪操作[^1]。 ```python from...
哨兵2号数据裁剪拼贴
最新发布
03-29
#### 使用 Python 和 GDAL 进行图像裁剪 如果偏好本地计算环境，则可采用 Python 结合 GDAL 库完成裁剪任务。以下是一个简单的例子[^3]: ```python from osgeo import gdal # 输入文件路径 input_raster = "path/to...
SARScape图像裁剪终极指南：你必须掌握的关键技术
![SARScape图像裁剪终极指南：你必须掌握的关键技术]...紧接着，本文通过实践操作部分深入分析了裁剪前的准备工作、SARScape裁剪工具的使用方法和裁剪后图像质量
【影像镶嵌与裁剪】：ERDAS遥感影像处理的艺术
遥感影像处理是地理信息系统和遥感技术领域的关键环节，它包括影像镶嵌与裁剪两大核心过程。本文首先概述了遥感影像处理的理论与实践，强调了影像镶嵌在遥感数据整合中的重要性，并介绍了影像裁剪的基本原理和工具...
如何使用GDAL进行图像镶嵌
weixin_34220834的博客
05-04 1350
相信看过和用过gdalwarp工具的同学都会对gdalwarp这个工具的强悍功能所震惊。今天主要就是用gdalwarp来进行图像镶嵌，当然这个镶嵌是比较简单的镶嵌，就是上层的图像会覆盖下层的图像。对于gdalwarp的介绍，我还要搬过来一下，见下面： The gdalwarp utility is an imagemosaicing, reproject...
Python+gdal裁剪遥感图像出现“Warning 1: TIFFReadDirectory:Sum of Photometric type-related color channels ……”
soderayer的博客
06-24 7726
运行环境：windows10、pycharm、python3.7、gdal等python使用GDAL 裁剪遥感图像的时候出现以下警告：“Warning 1: TIFFReadDirectory:Sum of Photometric type-related color channels and ExtraSamples doesn't match SamplesPerPixel. Defining non-color channels as ExtraSamples.”使用的代码段：Warp函数有很多...
python+ogr解决shp文件自相交的问题
qq_45723511的博客
03-18 3304
Warning1: RingSelf-intersectionatornearpoint112 .4866642030000334 .830899357000078 ERROR1: Cutlinepolygonisinvalid. 使用矢量shp文件剪裁栅格影像时会出现上述问题。提示shp文件自相交。搜集了大量资料，多数回答都是基于PostGIS进行处理。在下载安装过程中很复杂，本人试过多次也没有成功。这里给出一种使用python+ogr函数库的解决方法。 不废话，直接上代码。 from osgeo im
QGIS错误：裁剪矢量图形时出现“GEOS地理处理错误”的解决方法
codeuniverse的博客
10-24 1383
QGIS错误：裁剪矢量图形时出现“GEOS地理处理错误”的解决方法
python中读取栅格数据_使用rasterio处理python中的地理空间栅格数据
weixin_26713457的博客
09-01 6117
python中读取栅格数据This article is meant to provide a quick introduction into how to use the Python package Rasterio for common tasks related to geospatial raster data. This is mainly a collection of things...
分割模型nnUNet学习日记（二）：常见问题报错解决方案
Joker 007的博客
07-02 1451
分割模型nnUNet学习日记（二）：常见问题报错解决方案
评论 126
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值