{"id":839,"date":"2023-04-19T06:50:12","date_gmt":"2023-04-19T06:50:12","guid":{"rendered":"https:\/\/www.gislxz.com\/?p=839"},"modified":"2023-04-19T07:33:34","modified_gmt":"2023-04-19T07:33:34","slug":"zonalstatistics","status":"publish","type":"post","link":"https:\/\/www.gislxz.com\/index.php\/2023\/04\/19\/zonalstatistics\/","title":{"rendered":"C++ GDAL\u5f00\u53d1 \u5b9e\u73b0\u5206\u533a\u7edf\u8ba1\u529f\u80fd"},"content":{"rendered":"\n

\u4ee3\u7801\u4ed3\u5e93\uff1aZonalStatistics: GDAL\u5206\u533a\u7edf\u8ba1\u5b9e\u73b0 (gitee.com)<\/a><\/p>\n\n\n\n

\u80cc\u666f<\/h2>\n\n\n\n

\u5b9e\u73b0\u8fd9\u4e2a\u7b97\u6cd5\u7684\u60f3\u6cd5\u840c\u751f\u4e8e\u53bb\u5e74\uff0c\u6211\u7684\u5bfc\u5e08\u6307\u5bfc\u6211\u53bb\u534f\u52a9\u53c2\u4e0e\u91cd\u5e86\u519c\u4e1a\u9879\u76ee\uff0c\u4efb\u52a1\u662f\u66f4\u65b0\u91cd\u5e86\u5e02\u6240\u6709\u5730\u5757\u4e00\u5e74\u7684NDVI\u6570\u636e\uff0c\u8fd9\u4e2a\u4efb\u52a1\u770b\u4f3c\u7b80\u5355\u4f46\u5de5\u4f5c\u91cf\u6781\u5927\uff0c\u9664\u53bb\u6570\u636e\u51c6\u5907\u5916\uff0c\u8be5\u4efb\u52a1\u6700\u6838\u5fc3\u7684\u6b65\u9aa4\u5c31\u662f\u4f7f\u7528\u5206\u533a\u7edf\u8ba1\u5c06\u6805\u683c\u683c\u5f0f\u7684NDVI\u6570\u636e\u8d4b\u503c\u5230\u4ee5\u591a\u8fb9\u5f62\u4e3a\u7ec4\u7ec7\u5f62\u5f0f\u7684\u5730\u5757\u6570\u636e\u4e2d\u3002\u7136\u800c\u5c31\u662f\u8fd9\u6837\u7b80\u5355\u7684\u5206\u533a\u7edf\u8ba1\u529f\u80fd\u5374\u8017\u8d39\u4e86\u6211\u4eec\u5927\u91cf\u65f6\u95f4\uff0cArcGIS\u548cQGIS\u8f6f\u4ef6\u63d0\u4f9b\u7684\u5206\u533a\u7edf\u8ba1\u529f\u80fd\u5e76\u4e0d\u8db3\u4ee5\u652f\u6301\u6211\u4eec\u9ad8\u6548\u7cbe\u51c6\u7684\u5b8c\u6210\u8fd9\u9879\u4efb\u52a1\uff0c\u53ea\u80fd\u9000\u800c\u6c42\u5176\u6b21\u9009\u62e9\u4e86\u4e2d\u5fc3\u70b9\u8d4b\u503c\u6cd5\u5b8c\u6210\u4e86NDVI\u8d4b\u503c\u7684\u4efb\u52a1\uff0c\u4f46\u8fd9\u6837\u4e0d\u5229\u4e8e\u6570\u636e\u4e0d\u786e\u5b9a\u6027\u7684\u63a7\u5236\uff0c\u8fdb\u800c\u5f71\u54cd\u540e\u7eed\u8ba1\u7b97\u4e0e\u9884\u6d4b\uff0c\u9020\u6210\u8bef\u5dee\u7684\u4f20\u9012\u548c\u653e\u5927\u3002<\/p>\n\n\n\n

\u73b0\u6709GIS\u8f6f\u4ef6\u5206\u533a\u7edf\u8ba1\u529f\u80fd\u7684\u4e0d\u8db3<\/h2>\n\n\n\n

Arcgis\uff1a\u56e0\u4e3a\u6211\u53ea\u6709arcgis10.8\u7684\u7248\u672c\uff0c\u4e0d\u592a\u6e05\u695a\u65b0\u7684Arcgis Pro 3\u7684\u7b97\u6cd5\u5b9e\u73b0\u600e\u4e48\u6837\uff0c\u4f4610.8\u7248\u672c\u4e2d\u5206\u533a\u7edf\u8ba1\u7b97\u6cd5\u592a\u6162\u4e86\uff0c\u6bd4QGIS\u6162\u6570\u5341\u500d\uff0c\u8f93\u51fa\u7ed3\u679c\u662f\u4e00\u4e2a\u8868\u683c\uff0c\u8fd8\u8981\u518d\u8fdb\u884c\u5b57\u6bb5\u8fde\u63a5\uff0c\u4e0d\u592a\u884c\u3002<\/p>\n\n\n\n

QGIS\uff1aQGIS\u5904\u7406\u5012\u633a\u5feb\u7684\uff0c\u95ee\u9898\u662f\u4e0d\u80fd\u5199\u5165\u539f\u6587\u4ef6\uff0c\u53ea\u80fd\u8f93\u51fa\u4e00\u4e2a\u65b0\u77e2\u91cf\u6587\u4ef6\uff0c\u90a3\u6211\u4e00\u6b21\u8fd0\u7b9712\u4e2a\u6708\u7684\u6805\u683c\u56fe\u50cf\uff0c\u6ca1\u6cd5\u6279\u5904\u7406\u8fd0\u884c\u4e86\u3002<\/p>\n\n\n\n

\u601d\u6765\u60f3\u53bb\u8fd8\u662f\u81ea\u5df1\u91cd\u5199\u4e00\u4e2a<\/p>\n\n\n\n

\u4f7f\u7528GDALWarp\u6765\u5b9e\u73b0<\/h2>\n\n\n\n

\u4e00\u5f00\u59cb\u5012\u662f\u60f3\u7528\u73b0\u6210\u7684API\u6765\u5b9e\u73b0\u7684\uff0cGDALWarp\u4e2d\u6709\u56fe\u50cf\u88c1\u5207\u7684API\u53ef\u4ee5\u62ff\u6765\u7528\u3002\u5206\u533a\u7edf\u8ba1\u672c\u8d28\u4e0a\u5c31\u662f\u7edf\u8ba1\u6bcf\u4e00\u4e2a\u591a\u8fb9\u5f62\u5185\u90e8\u50cf\u5143\u70b9\u7684\u503c\uff0c\u5176\u5b9e\u4e5f\u5c31\u662f\u8981\u5bf9\u56fe\u50cf\u8fdb\u884c\u88c1\u5207\u3002\u53ea\u4e0d\u8fc7\u6211\u4eec\u4e0d\u9700\u8981\u5c06\u88c1\u5207\u540e\u7684\u56fe\u50cf\u4fdd\u5b58\u4e3a\u65b0\u6587\u4ef6\uff0c\u53ea\u7528\u628a\u50cf\u5143\u503c\u5b58\u50a8\u7684\u4e8c\u7ef4\u6570\u7ec4\u8f93\u5165\u5230\u5185\u5b58\u91cc\u5728\u5bf9\u6570\u7ec4\u8fdb\u884c\u8fd0\u7b97\u5c31\u597d\u4e86\u3002\u9057\u61be\u7684\u662fGDAL\u5b66\u4e60\u8d44\u6599\u592a\u5c11\uff0c\u674e\u6c11\u5f55\u5148\u751f\u5199\u7684\u300aGDAL\u6e90\u7801\u5256\u6790\u4e0e\u5f00\u53d1\u6307\u5357\u300b\u4ecb\u7ecd\u4e86Warp\u7684\u57fa\u672c\u7528\u6cd5\uff0c\u4f46\u662f\u5e76\u6ca1\u8bf4\u600e\u4e48\u8f93\u5165\u5230\u5185\u5b58\uff0c\u53ea\u8bf4\u4e86\u8f93\u5165\u5230\u5185\u5b58\u662f\u7528WarpRegionToBuffer\u8fd9\u4e2a\u6267\u884c\u8bed\u53e5\u3002\u6211\u4e5f\u662f\u8fdb\u884c\u4e86\u4e00\u756a\u5c1d\u8bd5\uff0c\u8fd9\u4e2aWarp\u6ce8\u610f\u4e8b\u9879\u4e3b\u8981\u6709\u4e24\u70b9\uff0c\u4e00\u662f\u8f93\u5165\u7684AOI\u8303\u56f4\u9700\u8981\u662f\u591a\u8fb9\u5f62\u6216\u591a\u591a\u8fb9\u5f62\uff0c\u5176\u5750\u6807\u5e94\u8be5\u662f\u50cf\u7d20\u884c\u5217\u53f7\u800c\u4e0d\u662f\u5730\u7406\u5750\u6807\uff0c\u4e14\u8981\u6c42\u662fWKTGeo\u683c\u5f0f\uff0c\u597d\u5728OGRGeometry\u7c7b\u662f\u6709exportToWkt\u7684\u51fd\u6570\u3002\u7b2c\u4e8c\u70b9\u5c31\u662f\u4f7f\u7528Warp\u505a\u88c1\u5207\u8fd8\u9700\u8981\u8bbe\u7f6e\u5750\u6807\u8f6c\u6362\u7684\u53c2\u6570\uff0c\u4e66\u4e2d\u7ed9\u7684\u793a\u4f8b\u662f\u7528\u7684GDALCreateGenImgProjTransformer2\u8fd9\u4e2a\u51fd\u6570\uff0c\u8f93\u5165\u8f93\u5165\u548c\u8f93\u51fa\u56fe\u50cf\u5bf9\u5e94\u7684GDALDataset\u5c31\u884c\u4e86\uff0c\u4f46\u6211\u4eec\u4e0d\u9700\u8981\u8f93\u51fa\u5230\u6587\u4ef6\uff0c\u81ea\u7136\u4e5f\u6ca1\u6709\u521b\u5efa\u8f93\u51fa\u6805\u683c\u56fe\u50cf\u7684GDALDataset\uff0c\u5982\u679c\u521b\u5efa\u52bf\u5fc5\u4f1a\u6781\u5927\u5f71\u54cd\u5904\u7406\u901f\u5ea6\uff0c\u4e8e\u662f\u6211\u53bb\u67e5\u4e86\u5176\u4ed6\u51e0\u4e2aGDALCreateGenImgProjTransformer\u51fd\u6570\uff0c\u53d1\u73b0\u53ef\u4ee5\u8fd9\u6837\u66ff\u4ee3\u3002<\/p>\n\n\n\n

\/\/\u4e66\u4e2d\u793a\u4f8b\u4ee3\u7801\nvoid *hTransformArg=GDALCreateGenImgProjTransformer2((GDALDatasetH)pSrcDS, (GDALDatasetH)pDstDS,NULL);\n\/\/\u53ea\u521b\u5efa\u8f93\u51fa\u56fe\u50cf\u7684\u5730\u7406\u8f6c\u6362\u516d\u53c2\u6570\uff0c\u8c03\u7528\u4e24\u6b21\u539f\u56fe\u50cf\u7684OGRSpatialReference\nvoid* hTransformArg = GDALCreateGenImgProjTransformer4((OGRSpatialReferenceH*)MyRasterdata->GetSpatialRef(), pSrcGeoTransform, (OGRSpatialReferenceH*)MyRasterdata->GetSpatialRef(), pDstGeoTransform, NULL);<\/code><\/pre>\n\n\n\n
\u8c03\u7528GDALWarp\u8fdb\u884c\u88c1\u526a\u7684\u5b8c\u6574\u51fd\u6570\u4ee3\u7801<\/p>\n\n\n\n
\/\/\u8c03\u7528GDALWarp\u8fdb\u884c\u88c1\u526a\uff08\u7ed3\u679c\u8f93\u51fa\u81f3\u5185\u5b58\uff0c\u76f8\u5f53\u4e8e\u63a9\u819c\u63d0\u53d6\uff09\ntemplate <typename T>\nT* ZonalStatistics<T>::imageCutByAoi(const OGRGeometry& poGeometry, int& size) {\n\tGDALAllRegister();\n\n\t\/\/Geometry\u8f6c\u884c\u5217\u53f7WKTPolygon\n\tOGRGeometry* pixelGeometry = NULL;\n\t\/\/OGRPolygon* Pixelpolygon = (OGRPolygon*)OGRGeometryFactory::createGeometry(wkbPolygon);\n\tMyParcelManager->PLtoXY(poGeometry, &pixelGeometry);\n\n\t\/\/\u83b7\u53d6\u8f6c\u6362\u53c2\u6570\n\tdouble pSrcGeoTransform[6] = { 0 };\n\tdouble pDstGeoTransform[6] = { 0 };\n\tMyRasterdata->GetGeoTransform(pSrcGeoTransform);\t\t\t\/\/\u56fe\u50cf\u7684\u4eff\u5c04\u53d8\u6362\u4fe1\u606f\n\tmemcpy(pDstGeoTransform, pSrcGeoTransform, sizeof(double) * 6);\n\n\t\/\/\u83b7\u53d6\u56fe\u50cf\u5c3a\u5bf8\n\tint iBandCount = MyRasterdata->GetRasterCount();\n\tint iSrcWidth = MyRasterdata->GetRasterXSize();\n\tint iSrcHeight = MyRasterdata->GetRasterYSize();\n\n\t\/\/OGRGeometry\u8f6cWKT\n\tstd::string temp = pixelGeometry->exportToWkt();\n\tchar* WKTGeo = new char[temp.size() + 1];\n\tstrcpy(WKTGeo, temp.c_str());\n\t\/\/std::cout << WKTGeo << std::endl;\n\tOGREnvelope poRect;\n\tpixelGeometry->getEnvelope(&poRect);\n\tint iDstWidth = static_cast<int>(poRect.MaxX - poRect.MinX);\n\tint iDstHeight = static_cast<int>(poRect.MaxY - poRect.MinY);\n\n\t\/\/ \u8bbe\u7f6e\u8f93\u51fa\u56fe\u50cf\u7684\u5de6\u4e0a\u89d2\u5750\u6807\n\tGDALApplyGeoTransform(pSrcGeoTransform, poRect.MinX, poRect.MinY, (&pDstGeoTransform[0]), &(pDstGeoTransform[3]));\n\n\t\/\/\u6784\u9020\u5750\u6807\u8f6c\u6362\u5173\u7cfb\n\tvoid* hTransformArg = GDALCreateGenImgProjTransformer4((OGRSpatialReferenceH*)MyRasterdata->GetSpatialRef(), pSrcGeoTransform, (OGRSpatialReferenceH*)MyRasterdata->GetSpatialRef(), pDstGeoTransform, NULL);\n\tGDALTransformerFunc pfnTransformer = GDALGenImgProjTransform;\n\n\tT* pBuf = (T*)CPLMalloc(sizeof(float) * iDstWidth * iDstHeight);\n\tsize = iDstWidth * iDstHeight;\n\n\t\/\/\u6784\u9020GDALWarp\u7684\u53d8\u6362\u9009\u9879\n\tGDALWarpOptions* psWO = GDALCreateWarpOptions();\n\tpsWO->papszWarpOptions = CSLDuplicate(NULL);\n\tpsWO->eWorkingDataType = uDT;\n\tpsWO->eResampleAlg = GRA_NearestNeighbour;\n\t\/\/\u6d4b\u8bd5\u53d1\u73b0\u4e0d\u80fd\u6b63\u5e38\u8bbe\u7f6enodata\u503c\uff0c\u8fd8\u672a\u89e3\u51b3\n\tdouble* nodata = new double{ static_cast<double>(NoDataValue) };\n\tpsWO->padfDstNoDataReal = nodata;\n\n\tpsWO->hSrcDS = (GDALDatasetH)(MyRasterdata.get());\n\tpsWO->pfnTransformer = pfnTransformer;\n\tpsWO->pTransformerArg = hTransformArg;\n\tpsWO->nBandCount = 1;\n\tpsWO->papszWarpOptions = CSLSetNameValue(psWO->papszWarpOptions, \"CUTLINE\", WKTGeo);\n\tpsWO->panSrcBands = (int*)CPLMalloc(psWO->nBandCount * sizeof(int));\n\tpsWO->panSrcBands[0] = 1;\n\t\/\/\u521b\u5efaGDALWarp\u6267\u884c\u5bf9\u8c61\uff0c\u5e76\u4f7f\u7528GDALWarpOptions\n\tGDALWarpOperation oWO;\n\toWO.Initialize(psWO);\n\n\n\t\/\/\u6267\u884c\u5904\u7406(\u8f93\u51fa\u5230\u5185\u5b58)\n\t\/\/oWO.ChunkAndWarpImage(0, 0, iDstWidth, iDstHeight);\n\tif (oWO.WarpRegionToBuffer(0, 0, iDstWidth, iDstHeight, pBuf, uDT)) {\n\t\tstd::cout << \"cut failed\" << std::endl;\n\t\tdelete[] WKTGeo;\n\t\tCPLFree(pBuf);\n\t\t\/\/ \u91ca\u653e\u8d44\u6e90\u548c\u5173\u95ed\u6587\u4ef6\n\t\tGDALDestroyGenImgProjTransformer(psWO->pTransformerArg);\n\t\tGDALDestroyWarpOptions(psWO);\n\t\texit(0);\n\t}\n\telse {\n\t\tdelete[] WKTGeo;\n\t\t\/\/ \u91ca\u653e\u8d44\u6e90\u548c\u5173\u95ed\u6587\u4ef6\n\t\tGDALDestroyGenImgProjTransformer(psWO->pTransformerArg);\n\t\tGDALDestroyWarpOptions(psWO);\n\t}\n\treturn pBuf;\n}<\/code><\/pre>\n\n\n\n
\u6d4b\u8bd5\u540e\u53d1\u73b0\u8c03\u7528GDALWarp\u6765\u5b9e\u73b0\u5b58\u5728\u4ee5\u4e0b\u51e0\u70b9\u4e0d\u8db3\uff1a<\/p>\n\n\n\n
1\u5904\u7406\u901f\u5ea6\u6bd4QGIS\u8fd8\u662f\u6162\u56db\u5230\u4e94\u500d<\/p>\n\n\n\n
2\u4e0d\u80fd\u7cbe\u786e\u6355\u83b7NODATA\u503c<\/p>\n\n\n\n
3 \u9762\u79ef\u8fc7\u5c0f\u7684\u5730\u5757\u53d1\u751f\u9519\u8bef<\/p>\n\n\n\n
\u4e8e\u662f\u6211\u53ea\u80fd\u5c1d\u8bd5\u4ece\u5934\u5b9e\u73b0\uff0c\u81ea\u5df1\u6765\u5199\u88c1\u5207\u529f\u80fd\u3002<\/p>\n\n\n\n
\u4ece\u5934\u5f00\u59cb\u5e95\u5c42\u5b9e\u73b0<\/h2>\n\n\n\n
\u7b97\u6cd5\u8bbe\u8ba1\u601d\u8def<\/h3>\n\n\n\n
\u5728\u5224\u65ad\u591a\u8fb9\u5f62\u4e0e\u6805\u683c\u56fe\u50cf\u76f8\u4ea4\u8303\u56f4\u4e2d\uff0c\u6211\u91c7\u7528\u4e86\u4e2d\u5fc3\u70b9\u5224\u65ad\u6cd5\uff0c\u53ca\u50cf\u5143\u4e2d\u5fc3\u70b9\u88ab\u591a\u8fb9\u5f62\u6240\u5305\u542b\uff0c\u5373\u8ba4\u4e3a\u591a\u8fb9\u5f62\u4e0e\u6805\u683c\u56fe\u50cf\u7684\u76f8\u4ea4\u8303\u56f4\u5305\u542b\u8be5\u50cf\u5143\u3002<\/p>\n\n\n\n
\"\"<\/figure><\/div>\n\n\n\n
\u8ba1\u7b97\u591a\u8fb9\u5f62\u4e0e\u6805\u683c\u56fe\u50cf\u76f8\u4ea4\u8303\u56f4\u6b65\u9aa4\uff1a<\/p>\n\n\n\n
1\u901a\u8fc7\u6805\u683c\u56fe\u50cf\u884c\u5217\u53f7\u8f6c\u5730\u7406\u5750\u6807\u7684\u8f6c\u6362\u516d\u53c2\u6570\u53cd\u63a8\u5730\u7406\u5750\u6807\u8f6c\u884c\u5217\u53f7\u7684\u516d\u53c2\u6570\u3002\uff08\u9700\u6ce8\u610f\u8f6c\u6362\u516d\u53c2\u6570\u4e2d\u7684\u8d77\u59cb\u70b9\u4e3a\u5de6\u4e0a\u89d2\u50cf\u7d20\u70b9\u7684\u5de6\u4e0a\u89d2\u4f4d\u7f6e\u800c\u975e\u4e2d\u5fc3\u70b9\uff0c\u56e0\u6b64\u53cd\u7b97\u65f6\u9700\u8981\u5c06\u8d77\u59cb\u70b9\u7684\u5750\u6807\u5206\u522b\u52a0\u4e0a\u5bf9\u5e94\u65b9\u5411\u50cf\u7d20\u5bbd\u5ea6\u7684\u4e00\u534a\uff09<\/p>\n\n\n\n
2\u5bf9\u4efb\u4e00\u591a\u8fb9\u5f62\u904d\u5386\u5176\u5916\u73af\uff08ExteriorRing\uff09\uff0c\u5c06\u5f97\u5230\u7684\u70b9\u5750\u6807\u8f6c\u4e3a\u884c\u5217\u53f7\uff0c\u901a\u8fc7\u524d\u540e\u4e24\u4e2a\u70b9\u8fde\u6210\u7684\u76f4\u7ebf\u8ba1\u7b97\u8de8\u8d8a\u884c\u53f7\u65f6\u4e0e\u5bf9\u5e94\u76f4\u7ebf\u7684\u4ea4\u70b9\u3002<\/p>\n\n\n\n
\"\"<\/figure><\/div>\n\n\n\n
3\u5bf9\u6bcf\u4e00\u884c\u7684\u4ea4\u70b9\u8fdb\u884c\u6392\u5e8f\uff0c\u4ece\u5c0f\u5230\u5927\u5373\u4ece\u5de6\u5230\u53f3\uff0c\u8fd9\u6837\u5305\u542b\u5728\u4fe9\u70b9\u95f4\u7684\u50cf\u7d20\u70b9\u5fc5\u7136\u90fd\u5728\u591a\u8fb9\u5f62\u5185\u3002\u5bf9\u4e8e\u51f9\u591a\u8fb9\u5f62\uff0c\u591a\u591a\u8fb9\u5f62\uff0c\u5185\u73af\u7b49\u60c5\u51b5\u8be5\u7b97\u6cd5\u90fd\u9002\u7528\u3002<\/p>\n\n\n\n
\u5728\u83b7\u5f97\u591a\u8fb9\u5f62\u5305\u542b\u7684\u50cf\u5143\u70b9\u540e\u5373\u53ef\u8bfb\u53d6\u6805\u683c\u5c06\u5bf9\u5e94\u50cf\u7d20\u70b9\u7684\u6570\u503c\u8d4b\u7ed9\u591a\u8fb9\u5f62\uff0c\u518d\u8fdb\u884c\u7edf\u8ba1\u8ba1\u7b97\u3002\u8fd9\u91cc\u6211\u6ce8\u610f\u5230\u5982\u679c\u904d\u5386\u591a\u8fb9\u5f62\u4ee5\u5176\u56db\u81f3\u8303\u56f4\u53bb\u8bfb\u53d6\u6805\u683c\u5e76\u8d4b\u503c\u4f1a\u9020\u6210\u6781\u5927\u7684\u786c\u76d8IO\u5f00\u9500\uff0c\u56e0\u800c\u5c06\u7b56\u7565\u8c03\u6574\u6210\u9010\u884c\u8bfb\u53d6\u6805\u683c\u56fe\u50cf\uff0c\u5e76\u5bf9\u8be5\u884c\u6240\u6709\u76f8\u4ea4\u7684\u591a\u8fb9\u5f62\u8fdb\u884c\u8d4b\u503c\u5904\u7406\u3002<\/p>\n\n\n\n
\u56de\u5fc6NDVI\u8d4b\u503c\u7684\u4efb\u52a1\uff0c\u6bcf\u6708\u5bf9\u5e94\u7684\u4e00\u5e45NDVI\u56fe\u50cf\u7684\u5c3a\u5bf8\uff0c\u89d2\u5ea6\u548c\u8303\u56f4\u53ef\u4ee5\u8ba4\u4e3a\u662f\u76f8\u540c\u7684\uff0c\u8fd9\u5c31\u610f\u5473\u7740\u6bcf\u4e2a\u591a\u8fb9\u5f62\u8ba1\u7b97\u7684\u76f8\u4ea4\uff08\u63a9\u819c\uff09\u8303\u56f4\u662f\u53ef\u4ee5\u91cd\u590d\u5229\u7528\u7684\uff0c\u56e0\u6b64\u6211\u53c8\u589e\u52a0\u4e86\u4e00\u4e2a\u6279\u5904\u7406\u52a0\u901f\u7684\u7b56\u7565\uff0c\u5f53\u8f93\u5165\u591a\u5e45\u56fe\u50cf\u65f6\u53ef\u4ee5\u4f7f\u7528\u6279\u5904\u7406\u52a0\u901f\uff0c\u5728\u5904\u7406\u5b8c\u6bcf\u4e2a\u5730\u5757\u65f6\u4e0d\u91ca\u653e\u76f8\u4ea4\u8303\u56f4\u548c\u8bb0\u5f55\u50cf\u5143\u503c\u7f13\u51b2\u533a\u800c\u53ea\u662f\u91cd\u7f6e\u6709\u6548\u50cf\u5143\u8ba1\u6570\uff0c\u7528\u66f4\u591a\u7684\u5185\u5b58\u5360\u7528\u6362\u53d6\u5904\u7406\u65f6\u95f4\u4e0a\u7684\u51cf\u5c11\u3002<\/p>\n\n\n\n
\u7b97\u6cd5\u7ed3\u6784\u8bbe\u8ba1<\/h3>\n\n\n\n
\"\"<\/figure><\/div>\n\n\n\n
\u8bf4\u660e\uff1a<\/p>\n\n\n\n
1\u6574\u4e2a\u4ee3\u7801\u4f7f\u7528Entrance\u51fd\u6570\u4f5c\u4e3a\u5165\u53e3\uff0c\u6253\u5f00\u7b2c\u4e00\u5e45\u6805\u683c\u56fe\u50cf\u4ee5\u83b7\u53d6\u6805\u683c\u6570\u636e\u7c7b\u578b\u4ece\u800c\u521b\u5efa\u76f8\u5e94ZonalStatistics\uff08\u6a21\u677f\u7c7b\uff09\uff0c\u5e76\u628a\u6253\u5f00\u6805\u683c\u56fe\u50cf\u6240\u7528\u7684GDALDatasetUniquePtr\u8f6c\u7ed9\uff08std::move\uff09ZonalStatistics\u3002<\/p>\n\n\n\n
2 ZonalStatistics\u7c7b\u4e2d\u5177\u6709\u4e00\u4e2aParcelManager\u6307\u9488\uff0c\u5728\u6784\u9020\u51fd\u6570\u4e2d\u4f1a\u6784\u9020\u4e00\u4e2a\u5176\u72ec\u6709\u7684ParcelManager\u3002ZonalStatistics\u7c7b\u5305\u542b\u4e86\u6805\u683c\u56fe\u50cf\u6587\u4ef6\u7684\u6307\u9488\uff0c\u9700\u8981\u6dfb\u52a0\u7684\u5c5e\u6027\u4fe1\u606f\uff0c\u5176\u6210\u5458\u51fd\u6570void poiStatistics()\u548cvoid MyStatistics()\u5206\u522b\u5b9e\u73b0\u4e86\u501f\u52a9GDALWrap\u88c1\u526a\u56fe\u50cf\u548c\u5e95\u5c42\u6784\u5efa\u7684\u5206\u533a\u7edf\u8ba1\u7b97\u6cd5\u3002<\/p>\n\n\n\n
3 ParcelManager\u4f5c\u4e3a\u7ba1\u7406ParcelUnit\u7684\u7c7b\uff0c\u5176\u5173\u952e\u6210\u5458\u5728\u4e8e\u5177\u6709\u4e24\u4e2astd::forward_list<ParcelUnit<T>*>**\u6210\u5458\u5206\u522b\u7528\u4e8e\u8bb0\u5f55\u6805\u683c\u56fe\u50cf\u67d0\u4e00\u884c\u9700\u8981\u542f\u52a8\u548c\u7ed3\u7b97\u7684\u5730\u5757\u5bf9\u8c61\uff0c\u5c06\u5176\u52a0\u5165\u5904\u7406\u961f\u5217\u6216\u4ece\u5904\u7406\u961f\u5217\u4e2d\u5220\u9664\u5e76\u8fdb\u884c\u7ed3\u7b97\uff08\u6279\u5904\u7406\u52a0\u901f\u4e2d\u53ea\u91cd\u7f6e\u8ba1\u6570\uff09\u3002\u5904\u7406\u961f\u5217\u4f7f\u7528std::unordered_set<ParcelUnit<T>*>*\u7ed3\u6784\u5b9e\u73b0\u65b9\u4fbf\u63d2\u5165\u548c\u5220\u9664\u3002<\/p>\n\n\n\n
4 ParcelUnit\u4f5c\u4e3a\u5730\u5757\u5355\u5143\uff0c\u4f7f\u7528\u81ea\u5b9a\u4e49\u7ed3\u6784\u4f53MyExtent\u8bb0\u5f55\u51e0\u4f55\u4f53\u8303\u56f4\uff0c\u4ee5\u4fbf\u63d2\u5165ParcelManager\u7684\u542f\u52a8\u548c\u7ed3\u7b97\u961f\u5217\uff1b\u4f7f\u7528std::vector<double>**\u8bb0\u5f55\u76f8\u4ea4\uff08\u63a9\u819c\uff09\u8303\u56f4\uff1b\u4f7f\u7528double\u6570\u7ec4\u8bb0\u5f55\u5305\u542b\u7684\u50cf\u5143\u503c\u3002\u540c\u65f6\u5176\u5177\u6709\u4e00\u4e2a\u4e2d\u5fc3\u70b9\u5750\u6807\u548c\u5bf9\u5e94\u50cf\u5143\u503c\uff08\u5982\u679c\u5728\u6805\u683c\u56fe\u50cf\u8303\u56f4\u5185\uff09\u4ee5\u4fbf\u5904\u7406\u5730\u5757\u8fc7\u5c0f\u4e0d\u5305\u542b\u4efb\u4e00\u50cf\u5143\u4e2d\u5fc3\u70b9\u7684\u60c5\u51b5\u3002<\/p>\n\n\n\n
\u7ed3\u679c\u5bf9\u6bd4<\/h3>\n\n\n\n
\u5b9e\u9a8c\u6570\u636e\u4e3a<\/p>\n\n\n\n
\u00b7\u91cd\u5e86\u5e022020\u5e7410\u6708\u53ca2020\u5e7411\u6708NDVI\u5408\u6210\u56fe\u50cf\uff08\u5355\u6ce2\u6bb5\uff0c\u5c3a\u5bf854615\u00d744970\uff0c\u5355\u5e45\u6570\u636e\u5927\u5c0f9.89GB\uff09<\/p>\n\n\n\n
\u00b7\u91cd\u5e86\u5e02\u7da6\u6c5f\u533a\u8015\u5730\u5730\u5757\u6570\u636e\uff0c\u5171761496\u4e2a\u51e0\u4f55\u4f53\uff0c\u5305\u542b\u51f9\u591a\u8fb9\u5f62\uff0c\u591a\u591a\u8fb9\u5f62\uff0c\u5185\u73af\u7b49\u5404\u79cd\u60c5\u51b5\u3002<\/p>\n\n\n\n
\u5355\u4f4d\uff08\u79d2\/s\uff09<\/td>\u5206\u9879<\/td>QGIS<\/td>\u5206\u9879<\/td>\u672c\u7b97\u6cd5<\/td><\/tr>
\u5904\u7406\u65f6\u95f4\uff08\u5355\u5e45\u56fe\u50cf\uff09<\/td>\u603b\u8ba1<\/td>107s<\/td>\u603b\u8ba1<\/td>116s<\/td><\/tr>
 <\/td>\u8ba1\u7b97\u5f97\u5230\u4e34\u65f6\u56fe\u5c42<\/td>64s<\/td>\u5411\u539f\u6587\u4ef6\u5c5e\u6027\u8868\u4e2d\u6dfb\u52a0\u5217<\/td>67s<\/td><\/tr>
 <\/td>\u5199\u5165\u786c\u76d8<\/td>43s<\/td>\u8ba1\u7b97\u5e76\u5199\u5165<\/td>49s<\/td><\/tr>
\u5904\u7406\u7b2c\u4e8c\u5e45\u56fe\u8c61<\/td>\u603b\u8ba1<\/td>121s<\/td>\u603b\u8ba1\uff08\u6279\u5904\u7406\u52a0\u901f\uff09<\/td>97s<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
1\u672c\u7b97\u6cd5\u4e0eQGIS\u5206\u533a\u7edf\u8ba1\u7b97\u6cd5\u8ba1\u7b97\u7ed3\u679c\u57fa\u672c\u76f8\u540c\uff0cQGIS\u5bf9\u5c0f\u5730\u5757\u6709\u7279\u6b8a\u5904\u7406\u7b56\u7565\uff0c\u56e0\u800c\u9020\u6210\u4e2a\u522b\u5730\u5757\u7ed3\u679c\u5b58\u5728\u4e9b\u8bb8\u5dee\u5f02<\/p>\n\n\n\n
2\u5728\u5904\u7406\u5355\u5e45\u56fe\u50cf\u4e2d\u672c\u7b97\u6cd5\u4e0eQGIS\u5904\u7406\u901f\u5ea6\u57fa\u672c\u4e00\u81f4\uff0c\u5728\u6279\u91cf\u5904\u7406\u4e2d\uff0cQGIS\u56e0\u4e3a\u5c5e\u6027\u8868\u4f53\u79ef\u7684\u589e\u957f\u800c\u51cf\u901f\uff0c\u672c\u7b97\u6cd5\u5728\u5411\u5c5e\u6027\u8868\u4e2d\u6dfb\u52a0\u5217\u6b65\u9aa4\u4e2d\u8017\u65f6\u6ca1\u6709\u660e\u663e\u53d8\u5316\u4e14\u7edf\u8ba1\u8ba1\u7b97\u6b65\u9aa4\u7531\u4e8e\u4f18\u5316\u7b56\u7565\u6709\u660e\u663e\u63d0\u901f\u3002<\/p>\n","protected":false},"excerpt":{"rendered":"
C++ GDAL\u5b9e\u73b0\u5206\u533a\u7edf\u8ba1<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[25],"tags":[],"class_list":["post-839","post","type-post","status-publish","format-standard","hentry","category-gis"],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/posts\/839","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/comments?post=839"}],"version-history":[{"count":2,"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/posts\/839\/revisions"}],"predecessor-version":[{"id":844,"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/posts\/839\/revisions\/844"}],"wp:attachment":[{"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/media?parent=839"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/categories?post=839"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.gislxz.com\/index.php\/wp-json\/wp\/v2\/tags?post=839"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}