GDAL能正确运行的代码：使用不同ZOOM读取瓦片数据，转换为PNG

我用AI搜索了很久总算完成了这个代码。工作思路无人机在空中飞行需要匹配地面影像确认自己的真实位置。瓦片数据使用QGIS下载瓦片地图-CSDN博客代码能编译运行能正确产生不同ZOOM的图片。#include stdio.h #include stdlib.h #include gdal.h #include gdal_utils.h #define MAPTILE_FILE wuxi.mbtiles #define POINT_LONGITUDE 120.337 #define POINT_LATITUDE 31.505 #define AREA_SIZE 0.010 static GDALDatasetH openMapFile(const int zoom) { #if 1 char param1[32]; sprintf(param1, ZOOM_LEVEL%d, zoom); const char* openOpts[] {param1, NULL}; // 3. 使用GDALOpenEx打开数据集 GDALDatasetH hSrcDS GDALOpenEx( MAPTILE_FILE, // 文件名 GDAL_OF_RASTER | GDAL_OF_READONLY, // 只读方式打开栅格数据 NULL, // 不限制驱动 openOpts, // 传入我们设置的打开选项 NULL // 无关联文件 ); #else GDALDatasetH hSrcDS GDALOpen(MAPTILE_FILE, GA_ReadOnly); #endif return hSrcDS; } static int convertToPng(const char* file, const GDALDatasetH hSrcDS, const GDALTranslateOptions* psOptions) { // 6. 执行裁剪和格式转换 GDALDatasetH hDstDS GDALTranslate(file, hSrcDS, psOptions, NULL); if (hDstDS) { printf(成功将 %s 裁剪并保存为 %s\n, MAPTILE_FILE, file); return -1; } else { printf(转换失败\n); } return 0; } static GDALTranslateOptions* initTransOptions() { // 2, 将边界值转换为字符串并加入到参数数组中 char szMinLon[32], szMaxLat[32], szMaxLon[32], szMinLat[32]; snprintf(szMinLon, sizeof(szMinLon), %.6f, (POINT_LONGITUDE-AREA_SIZE/2)); snprintf(szMaxLat, sizeof(szMaxLat), %.6f, (POINT_LATITUDE AREA_SIZE/2)); snprintf(szMaxLon, sizeof(szMaxLon), %.6f, (POINT_LONGITUDEAREA_SIZE/2)); snprintf(szMinLat, sizeof(szMinLat), %.6f, (POINT_LATITUDE -AREA_SIZE/2)); // 3. 构建 参数数组 // 重新构建包含边界值的参数数组 // 示例范围围绕北京的区域 (经度: 116.3 ~ 116.5, 纬度: 39.8 ~ 40.0) // 注意-projwin 的顺序是左上角 (minLon, maxLat) 到右下角 (maxLon, minLat) const char* const apszArgvFull[] { -projwin, szMinLon, szMaxLat, szMaxLon, szMinLat, -projwin_srs, EPSG:4326, -of, PNG, NULL }; // 4. 创建 GDALTranslateOptions 对象 GDALTranslateOptions* psOptions GDALTranslateOptionsNew((char**)apszArgvFull, NULL); if (!psOptions) { printf(解析 GDALTranslateOptions 参数失败\n); return NULL; } return psOptions; } int main() { // 1. 注册所有 GDAL 驱动 GDALAllRegister(); GDALTranslateOptions* psOptions initTransOptions(); if (!psOptions) { return -1; } for (int zoom8; zoom16; zoom) { // 5. 打开源数据集 (MBTiles) GDALDatasetH hSrcDS openMapFile(zoom); if (hSrcDS NULL) { printf(无法打开源文件: %s(zoom%d)\n, MAPTILE_FILE, zoom); continue; } char file[32]; sprintf(file, zoom%d.png, zoom); convertToPng(file, hSrcDS, psOptions); GDALClose(hSrcDS); } // 7. 清理资源 GDALTranslateOptionsFree(psOptions); GDALClearMemoryCaches(); GDALDestroy(); // Why core dump ? return 0; }编译脚本export LD_LIBRARY_PATH${LD_LIBRARY_PATH}:${HOME}/map_gdal/lib:${HOME}/map_proj/lib export PKG_CONFIG_PATH${PKG_CONFIG_PATH}:${HOME}/map_gdal/lib/pkgconfig rm core.* reset SRC_FILEtest2.cpp EXE_FILEtest if [ -f ${EXE_FILE} ]; then rm ${EXE_FILE} fi g ${SRC_FILE} \ -o ${EXE_FILE} \ $(pkg-config --cflags --libs gdal) ./${EXE_FILE}