{"id":1062376,"date":"2024-12-31T15:51:08","date_gmt":"2024-12-31T07:51:08","guid":{"rendered":"https:\/\/docs.pingcode.com\/ask\/ask-ask\/1062376.html"},"modified":"2024-12-31T15:51:11","modified_gmt":"2024-12-31T07:51:11","slug":"python%e5%a6%82%e4%bd%95%e6%a0%b9%e6%8d%ae%e5%9d%90%e6%a0%87%e7%82%b9%e8%af%bb%e5%8f%96%e6%a0%85%e6%a0%bc%e5%80%bc","status":"publish","type":"post","link":"https:\/\/docs.pingcode.com\/ask\/ask-ask\/1062376.html","title":{"rendered":"python\u5982\u4f55\u6839\u636e\u5750\u6807\u70b9\u8bfb\u53d6\u6805\u683c\u503c"},"content":{"rendered":"<p style=\"text-align:center;\" ><img decoding=\"async\" src=\"https:\/\/cdn-docs.pingcode.com\/wp-content\/uploads\/2024\/12\/1e04a216-3f67-4b1b-bafc-af03c8e599c6.webp?x-oss-process=image\/auto-orient,1\/format,webp\" alt=\"python\u5982\u4f55\u6839\u636e\u5750\u6807\u70b9\u8bfb\u53d6\u6805\u683c\u503c\" \/><\/p>\n<p><p> <strong>\u5982\u4f55\u7528Python\u6839\u636e\u5750\u6807\u70b9\u8bfb\u53d6\u6805\u683c\u503c<\/strong><\/p>\n<\/p>\n<p><p><strong>\u4f7f\u7528GDAL\u5e93\u3001\u4f7f\u7528Rasterio\u5e93\u3001\u5904\u7406\u5750\u6807\u7cfb\u8f6c\u6362<\/strong><\/p>\n<\/p>\n<p><p>\u5728Python\u4e2d\uff0c\u53ef\u4ee5\u4f7f\u7528GDAL\u548cRasterio\u5e93\u6839\u636e\u5750\u6807\u70b9\u8bfb\u53d6\u6805\u683c\u503c\u3002<strong>GDAL\u5e93\u5f3a\u5927\u4e14\u5e7f\u6cdb\u652f\u6301\u4e0d\u540c\u7684\u6805\u683c\u6570\u636e\u683c\u5f0f\u3001Rasterio\u5e93\u6613\u4e8e\u4f7f\u7528\u4e14\u66f4\u73b0\u4ee3\u5316\u3001\u5904\u7406\u5750\u6807\u7cfb\u8f6c\u6362\u53ef\u4ee5\u786e\u4fdd\u5750\u6807\u70b9\u4e0e\u6805\u683c\u6570\u636e\u7684\u5750\u6807\u7cfb\u4e00\u81f4<\/strong>\u3002\u63a5\u4e0b\u6765\uff0c\u6211\u4eec\u5c06\u8be6\u7ec6\u4ecb\u7ecd\u8fd9\u4e09\u79cd\u65b9\u6cd5\u3002<\/p>\n<\/p>\n<p><h3>\u4e00\u3001\u4f7f\u7528GDAL\u5e93<\/h3>\n<\/p>\n<p><p>GDAL\uff08Geospatial Data Abstraction Library\uff09\u662f\u4e00\u4e2a\u5f3a\u5927\u7684\u5e93\uff0c\u652f\u6301\u591a\u79cd\u6805\u683c\u548c\u77e2\u91cf\u6570\u636e\u683c\u5f0f\u3002\u4f7f\u7528GDAL\u5e93\u53ef\u4ee5\u65b9\u4fbf\u5730\u8bfb\u53d6\u6805\u683c\u6570\u636e\u5e76\u83b7\u53d6\u6307\u5b9a\u5750\u6807\u70b9\u7684\u6805\u683c\u503c\u3002<\/p>\n<\/p>\n<p><h4>1\u3001\u5b89\u88c5GDAL\u5e93<\/h4>\n<\/p>\n<p><p>\u9996\u5148\uff0c\u786e\u4fdd\u4f60\u5df2\u7ecf\u5b89\u88c5\u4e86GDAL\u5e93\u3002\u4f60\u53ef\u4ee5\u4f7f\u7528\u4ee5\u4e0b\u547d\u4ee4\u5b89\u88c5\uff1a<\/p>\n<\/p>\n<p><pre><code class=\"language-bash\">pip install gdal<\/p>\n<p><\/code><\/pre>\n<\/p>\n<p><h4>2\u3001\u8bfb\u53d6\u6805\u683c\u6570\u636e<\/h4>\n<\/p>\n<p><p>\u4f7f\u7528GDAL\u5e93\u8bfb\u53d6\u6805\u683c\u6570\u636e\uff0c\u5e76\u6839\u636e\u6307\u5b9a\u7684\u5750\u6807\u70b9\u83b7\u53d6\u6805\u683c\u503c\u3002\u4ee5\u4e0b\u662f\u4e00\u4e2a\u793a\u4f8b\u4ee3\u7801\uff1a<\/p>\n<\/p>\n<p><pre><code class=\"language-python\">from osgeo import gdal<\/p>\n<p>def get_raster_value(raster_path, lon, lat):<\/p>\n<p>    dataset = gdal.Open(raster_path)<\/p>\n<p>    geotransform = dataset.GetGeoTransform()<\/p>\n<p>    # \u83b7\u53d6\u6805\u683c\u7684\u5143\u4fe1\u606f<\/p>\n<p>    origin_x = geotransform[0]<\/p>\n<p>    origin_y = geotransform[3]<\/p>\n<p>    pixel_width = geotransform[1]<\/p>\n<p>    pixel_height = geotransform[5]<\/p>\n<p>    # \u5c06\u5750\u6807\u8f6c\u6362\u4e3a\u6805\u683c\u50cf\u7d20\u5750\u6807<\/p>\n<p>    x_offset = int((lon - origin_x) \/ pixel_width)<\/p>\n<p>    y_offset = int((lat - origin_y) \/ pixel_height)<\/p>\n<p>    band = dataset.GetRasterBand(1)<\/p>\n<p>    raster_value = band.ReadAsArray(x_offset, y_offset, 1, 1)<\/p>\n<p>    return raster_value[0][0]<\/p>\n<h2><strong>\u793a\u4f8b\u4f7f\u7528<\/strong><\/h2>\n<p>raster_path = &#39;path\/to\/your\/raster.tif&#39;<\/p>\n<p>lon = 120.0<\/p>\n<p>lat = 30.0<\/p>\n<p>value = get_raster_value(raster_path, lon, lat)<\/p>\n<p>print(f&#39;The raster value at ({lon}, {lat}) is {value}&#39;)<\/p>\n<p><\/code><\/pre>\n<\/p>\n<p><p>\u5728\u4e0a\u9762\u7684\u4ee3\u7801\u4e2d\uff0c\u6211\u4eec\u9996\u5148\u4f7f\u7528<code>gdal.Open<\/code>\u51fd\u6570\u6253\u5f00\u6805\u683c\u6587\u4ef6\uff0c\u7136\u540e\u83b7\u53d6\u6805\u683c\u7684\u5730\u7406\u8f6c\u6362\u4fe1\u606f\uff08GeoTransform\uff09\u3002\u63a5\u7740\uff0c\u6211\u4eec\u5c06\u6307\u5b9a\u7684\u7ecf\u7eac\u5ea6\u5750\u6807\u8f6c\u6362\u4e3a\u6805\u683c\u50cf\u7d20\u5750\u6807\uff0c\u5e76\u4f7f\u7528<code>ReadAsArray<\/code>\u51fd\u6570\u8bfb\u53d6\u6805\u683c\u503c\u3002<\/p>\n<\/p>\n<p><h3>\u4e8c\u3001\u4f7f\u7528Rasterio\u5e93<\/h3>\n<\/p>\n<p><p>Rasterio\u662f\u4e00\u4e2a\u7528\u4e8e\u8bfb\u53d6\u548c\u5199\u5165\u6805\u683c\u6570\u636e\u7684\u73b0\u4ee3\u5316\u5e93\uff0c\u6613\u4e8e\u4f7f\u7528\u4e14\u529f\u80fd\u5f3a\u5927\u3002\u4f7f\u7528Rasterio\u5e93\u4e5f\u53ef\u4ee5\u65b9\u4fbf\u5730\u8bfb\u53d6\u6805\u683c\u6570\u636e\u5e76\u83b7\u53d6\u6307\u5b9a\u5750\u6807\u70b9\u7684\u6805\u683c\u503c\u3002<\/p>\n<\/p>\n<p><h4>1\u3001\u5b89\u88c5Rasterio\u5e93<\/h4>\n<\/p>\n<p><p>\u9996\u5148\uff0c\u786e\u4fdd\u4f60\u5df2\u7ecf\u5b89\u88c5\u4e86Rasterio\u5e93\u3002\u4f60\u53ef\u4ee5\u4f7f\u7528\u4ee5\u4e0b\u547d\u4ee4\u5b89\u88c5\uff1a<\/p>\n<\/p>\n<p><pre><code class=\"language-bash\">pip install rasterio<\/p>\n<p><\/code><\/pre>\n<\/p>\n<p><h4>2\u3001\u8bfb\u53d6\u6805\u683c\u6570\u636e<\/h4>\n<\/p>\n<p><p>\u4f7f\u7528Rasterio\u5e93\u8bfb\u53d6\u6805\u683c\u6570\u636e\uff0c\u5e76\u6839\u636e\u6307\u5b9a\u7684\u5750\u6807\u70b9\u83b7\u53d6\u6805\u683c\u503c\u3002\u4ee5\u4e0b\u662f\u4e00\u4e2a\u793a\u4f8b\u4ee3\u7801\uff1a<\/p>\n<\/p>\n<p><pre><code class=\"language-python\">import rasterio<\/p>\n<p>from rasterio.transform import from_origin<\/p>\n<p>def get_raster_value(raster_path, lon, lat):<\/p>\n<p>    with rasterio.open(raster_path) as dataset:<\/p>\n<p>        # \u5c06\u5750\u6807\u8f6c\u6362\u4e3a\u6805\u683c\u50cf\u7d20\u5750\u6807<\/p>\n<p>        row, col = dataset.index(lon, lat)<\/p>\n<p>        # \u8bfb\u53d6\u6805\u683c\u503c<\/p>\n<p>        raster_value = dataset.read(1)[row, col]<\/p>\n<p>    return raster_value<\/p>\n<h2><strong>\u793a\u4f8b\u4f7f\u7528<\/strong><\/h2>\n<p>raster_path = &#39;path\/to\/your\/raster.tif&#39;<\/p>\n<p>lon = 120.0<\/p>\n<p>lat = 30.0<\/p>\n<p>value = get_raster_value(raster_path, lon, lat)<\/p>\n<p>print(f&#39;The raster value at ({lon}, {lat}) is {value}&#39;)<\/p>\n<p><\/code><\/pre>\n<\/p>\n<p><p>\u5728\u4e0a\u9762\u7684\u4ee3\u7801\u4e2d\uff0c\u6211\u4eec\u4f7f\u7528<code>rasterio.open<\/code>\u51fd\u6570\u6253\u5f00\u6805\u683c\u6587\u4ef6\uff0c\u5e76\u4f7f\u7528<code>dataset.index<\/code>\u51fd\u6570\u5c06\u6307\u5b9a\u7684\u7ecf\u7eac\u5ea6\u5750\u6807\u8f6c\u6362\u4e3a\u6805\u683c\u50cf\u7d20\u5750\u6807\u3002\u63a5\u7740\uff0c\u6211\u4eec\u4f7f\u7528<code>dataset.read<\/code>\u51fd\u6570\u8bfb\u53d6\u6805\u683c\u503c\u3002<\/p>\n<\/p>\n<p><h3>\u4e09\u3001\u5904\u7406\u5750\u6807\u7cfb\u8f6c\u6362<\/h3>\n<\/p>\n<p><p>\u5728\u5b9e\u9645\u5e94\u7528\u4e2d\uff0c\u6805\u683c\u6570\u636e\u548c\u5750\u6807\u70b9\u53ef\u80fd\u4f7f\u7528\u4e0d\u540c\u7684\u5750\u6807\u7cfb\u3002\u4e3a\u4e86\u786e\u4fdd\u5750\u6807\u70b9\u4e0e\u6805\u683c\u6570\u636e\u7684\u5750\u6807\u7cfb\u4e00\u81f4\uff0c\u6211\u4eec\u9700\u8981\u8fdb\u884c\u5750\u6807\u7cfb\u8f6c\u6362\u3002\u53ef\u4ee5\u4f7f\u7528Pyproj\u5e93\u8fdb\u884c\u5750\u6807\u7cfb\u8f6c\u6362\u3002<\/p>\n<\/p>\n<p><h4>1\u3001\u5b89\u88c5Pyproj\u5e93<\/h4>\n<\/p>\n<p><p>\u9996\u5148\uff0c\u786e\u4fdd\u4f60\u5df2\u7ecf\u5b89\u88c5\u4e86Pyproj\u5e93\u3002\u4f60\u53ef\u4ee5\u4f7f\u7528\u4ee5\u4e0b\u547d\u4ee4\u5b89\u88c5\uff1a<\/p>\n<\/p>\n<p><pre><code class=\"language-bash\">pip install pyproj<\/p>\n<p><\/code><\/pre>\n<\/p>\n<p><h4>2\u3001\u8fdb\u884c\u5750\u6807\u7cfb\u8f6c\u6362<\/h4>\n<\/p>\n<p><p>\u4f7f\u7528Pyproj\u5e93\u8fdb\u884c\u5750\u6807\u7cfb\u8f6c\u6362\uff0c\u5e76\u7ed3\u5408Rasterio\u5e93\u83b7\u53d6\u6805\u683c\u503c\u3002\u4ee5\u4e0b\u662f\u4e00\u4e2a\u793a\u4f8b\u4ee3\u7801\uff1a<\/p>\n<\/p>\n<p><pre><code class=\"language-python\">import rasterio<\/p>\n<p>from pyproj import Transformer<\/p>\n<p>def get_raster_value_with_projection(raster_path, lon, lat, src_crs=&#39;EPSG:4326&#39;, dst_crs=&#39;EPSG:3857&#39;):<\/p>\n<p>    # \u521b\u5efa\u5750\u6807\u7cfb\u8f6c\u6362\u5668<\/p>\n<p>    transformer = Transformer.from_crs(src_crs, dst_crs, always_xy=True)<\/p>\n<p>    # \u8f6c\u6362\u5750\u6807\u7cfb<\/p>\n<p>    x, y = transformer.transform(lon, lat)<\/p>\n<p>    with rasterio.open(raster_path) as dataset:<\/p>\n<p>        # \u5c06\u8f6c\u6362\u540e\u7684\u5750\u6807\u8f6c\u6362\u4e3a\u6805\u683c\u50cf\u7d20\u5750\u6807<\/p>\n<p>        row, col = dataset.index(x, y)<\/p>\n<p>        # \u8bfb\u53d6\u6805\u683c\u503c<\/p>\n<p>        raster_value = dataset.read(1)[row, col]<\/p>\n<p>    return raster_value<\/p>\n<h2><strong>\u793a\u4f8b\u4f7f\u7528<\/strong><\/h2>\n<p>raster_path = &#39;path\/to\/your\/raster.tif&#39;<\/p>\n<p>lon = 120.0<\/p>\n<p>lat = 30.0<\/p>\n<p>value = get_raster_value_with_projection(raster_path, lon, lat)<\/p>\n<p>print(f&#39;The raster value at ({lon}, {lat}) is {value}&#39;)<\/p>\n<p><\/code><\/pre>\n<\/p>\n<p><p>\u5728\u4e0a\u9762\u7684\u4ee3\u7801\u4e2d\uff0c\u6211\u4eec\u4f7f\u7528Pyproj\u5e93\u7684<code>Transformer<\/code>\u7c7b\u521b\u5efa\u4e86\u4e00\u4e2a\u5750\u6807\u7cfb\u8f6c\u6362\u5668\uff0c\u5c06\u6307\u5b9a\u7684\u7ecf\u7eac\u5ea6\u5750\u6807\u4ece\u6e90\u5750\u6807\u7cfb\uff08src_crs\uff09\u8f6c\u6362\u4e3a\u76ee\u6807\u5750\u6807\u7cfb\uff08dst_crs\uff09\u3002\u7136\u540e\uff0c\u6211\u4eec\u4f7f\u7528\u8f6c\u6362\u540e\u7684\u5750\u6807\u8bfb\u53d6\u6805\u683c\u503c\u3002<\/p>\n<\/p>\n<p><h3>\u56db\u3001\u603b\u7ed3<\/h3>\n<\/p>\n<p><p>\u901a\u8fc7\u672c\u6587\u7684\u4ecb\u7ecd\uff0c\u6211\u4eec\u5b66\u4e60\u4e86\u5982\u4f55\u4f7f\u7528Python\u6839\u636e\u5750\u6807\u70b9\u8bfb\u53d6\u6805\u683c\u503c\u3002\u5177\u4f53\u6765\u8bf4\uff0c\u6211\u4eec\u4ecb\u7ecd\u4e86\u4f7f\u7528GDAL\u5e93\u548cRasterio\u5e93\u7684\u65b9\u6cd5\uff0c\u5e76\u8ba8\u8bba\u4e86\u5982\u4f55\u5904\u7406\u5750\u6807\u7cfb\u8f6c\u6362\u3002<strong>GDAL\u5e93\u5f3a\u5927\u4e14\u5e7f\u6cdb\u652f\u6301\u4e0d\u540c\u7684\u6805\u683c\u6570\u636e\u683c\u5f0f\u3001Rasterio\u5e93\u6613\u4e8e\u4f7f\u7528\u4e14\u66f4\u73b0\u4ee3\u5316\u3001\u5904\u7406\u5750\u6807\u7cfb\u8f6c\u6362\u53ef\u4ee5\u786e\u4fdd\u5750\u6807\u70b9\u4e0e\u6805\u683c\u6570\u636e\u7684\u5750\u6807\u7cfb\u4e00\u81f4<\/strong>\u3002\u5e0c\u671b\u672c\u6587\u5bf9\u4f60\u6709\u6240\u5e2e\u52a9\uff0c\u80fd\u591f\u5728\u5b9e\u9645\u5e94\u7528\u4e2d\u66f4\u597d\u5730\u5904\u7406\u6805\u683c\u6570\u636e\u3002<\/p>\n<\/p>\n<h2><strong>\u76f8\u5173\u95ee\u7b54FAQs\uff1a<\/strong><\/h2>\n<p> <strong>\u5982\u4f55\u4f7f\u7528Python\u8bfb\u53d6\u6805\u683c\u6570\u636e\u4e2d\u7684\u7279\u5b9a\u5750\u6807\u7684\u503c\uff1f<\/strong><br \/>\u5728Python\u4e2d\uff0c\u53ef\u4ee5\u4f7f\u7528\u591a\u4e2a\u5e93\u6765\u8bfb\u53d6\u6805\u683c\u6570\u636e\uff0c\u4f8b\u5982Rasterio\u548cGDAL\u3002\u4f7f\u7528Rasterio\uff0c\u4f60\u53ef\u4ee5\u6253\u5f00\u6805\u683c\u6587\u4ef6\uff0c\u5e76\u901a\u8fc7\u63d0\u4f9b\u7684\u5750\u6807\u8f6c\u6362\u4e3a\u884c\u5217\u7d22\u5f15\uff0c\u4ece\u800c\u83b7\u53d6\u76f8\u5e94\u7684\u6805\u683c\u503c\u3002\u4ee3\u7801\u793a\u4f8b\u5982\u4e0b\uff1a<\/p>\n<pre><code class=\"language-python\">import rasterio\n\nwith rasterio.open(&#39;your_raster_file.tif&#39;) as src:\n    row, col = src.index(x_coordinate, y_coordinate)\n    raster_value = src.read(1)[row, col]\n<\/code><\/pre>\n<p><strong>Python\u4e2d\u6709\u54ea\u4e9b\u5e93\u53ef\u4ee5\u5904\u7406\u6805\u683c\u6570\u636e\uff1f<\/strong><br \/>Python\u4e2d\u5904\u7406\u6805\u683c\u6570\u636e\u7684\u4e3b\u8981\u5e93\u5305\u62ecRasterio\u3001GDAL\u3001NumPy\u548cMatplotlib\u3002\u5176\u4e2d\uff0cRasterio\u548cGDAL\u662f\u6700\u5e38\u7528\u7684\u6805\u683c\u6570\u636e\u8bfb\u53d6\u548c\u5199\u5165\u5de5\u5177\uff0cNumPy\u63d0\u4f9b\u4e86\u5f3a\u5927\u7684\u6570\u7ec4\u64cd\u4f5c\u529f\u80fd\uff0c\u800cMatplotlib\u5219\u53ef\u4ee5\u7528\u6765\u53ef\u89c6\u5316\u6805\u683c\u6570\u636e\u3002<\/p>\n<p><strong>\u5982\u4f55\u5904\u7406\u8bfb\u53d6\u6805\u683c\u503c\u65f6\u7684\u5750\u6807\u7cfb\u7edf\u4e0d\u540c\u95ee\u9898\uff1f<\/strong><br \/>\u5728\u5904\u7406\u6805\u683c\u6570\u636e\u65f6\uff0c\u4e0d\u540c\u7684\u5750\u6807\u7cfb\u7edf\u53ef\u80fd\u4f1a\u5bfc\u81f4\u8bfb\u53d6\u7684\u503c\u4e0d\u51c6\u786e\u3002\u4e3a\u4e86\u907f\u514d\u8fd9\u79cd\u60c5\u51b5\uff0c\u5efa\u8bae\u5728\u8bfb\u53d6\u6805\u683c\u6570\u636e\u524d\uff0c\u5148\u4e86\u89e3\u6805\u683c\u7684\u5750\u6807\u53c2\u8003\u7cfb\u7edf\uff0c\u5e76\u4f7f\u7528\u76f8\u5e94\u7684\u8f6c\u6362\u5de5\u5177\uff08\u5982pyproj\u5e93\uff09\u5c06\u76ee\u6807\u5750\u6807\u8f6c\u6362\u4e3a\u6805\u683c\u7684\u5750\u6807\u7cfb\u7edf\uff0c\u786e\u4fdd\u8bfb\u53d6\u503c\u7684\u51c6\u786e\u6027\u3002<\/p>\n","protected":false},"excerpt":{"rendered":"\u5982\u4f55\u7528Python\u6839\u636e\u5750\u6807\u70b9\u8bfb\u53d6\u6805\u683c\u503c \u4f7f\u7528GDAL\u5e93\u3001\u4f7f\u7528Rasterio\u5e93\u3001\u5904\u7406\u5750\u6807\u7cfb\u8f6c\u6362 \u5728Python [&hellip;]","protected":false},"author":3,"featured_media":1062381,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[37],"tags":[],"acf":[],"_links":{"self":[{"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/posts\/1062376"}],"collection":[{"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/comments?post=1062376"}],"version-history":[{"count":"1","href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/posts\/1062376\/revisions"}],"predecessor-version":[{"id":1062382,"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/posts\/1062376\/revisions\/1062382"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/media\/1062381"}],"wp:attachment":[{"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/media?parent=1062376"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/categories?post=1062376"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/docs.pingcode.com\/wp-json\/wp\/v2\/tags?post=1062376"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}