Python提取多边形内外NDVI值方法
本文介绍了如何使用Python编程提取栅格图像（如NDVI图像）中，多个多边形内部以及多边形外部的平均NDVI值。主要依赖rasterio和fiona库，通过加载矢量多边形数据，裁剪栅格图像，并计算裁剪区域的平均值，从而实现NDVI值的提取和分析。本文提供代码示例，帮助读者理解和应用相关技术。
本文介绍了如何使用 Python 编程提取栅格图像（如 NDVI 图像）中，多个多边形内部以及多边形外部的平均 NDVI 值。主要依赖 rasterio 和 fiona 库，通过加载矢量多边形数据，裁剪栅格图像，并计算裁剪区域的平均值，从而实现 NDVI 值的提取和分析。本文提供代码示例，帮助读者理解和应用相关技术。
正文
本文将指导您如何使用 Python 提取栅格数据（例如 NDVI 图像）中多边形内部和外部的平均值。 我们将使用 rasterio 和 fiona 库来完成此任务。 rasterio 用于读取和写入栅格数据，而 fiona 用于读取矢量数据（例如 shapefile）。
准备工作
在开始之前，请确保已安装以下 Python 库：
- rasterio
- fiona
- numpy
您可以使用 pip 安装这些库：
pip install rasterio fiona numpy
代码实现
以下代码演示了如何提取多边形内部和外部的平均 NDVI 值：
import rasterio
import fiona
import rasterio.mask
import numpy as np
# 1. 定义输入文件路径
shapefile_path = "path/to/your/shapefile.shp" # 替换为你的 shapefile 文件路径
raster_path = "path/to/your/ndvi.tif" # 替换为你的 NDVI 栅格文件路径
# 2. 加载 shapefile 中的多边形
with fiona.open(shapefile_path, "r") as sf:
shapes = [feature["geometry"] for feature in sf]
# 3. 读取栅格数据
with rasterio.open(raster_path) as src:
# 3.1 提取多边形内部的 NDVI 值
out_image, out_transform = rasterio.mask.mask(src, shapes, crop=True)
# 清理无效值（例如 NaN 或 NoData 值）
masked_array = np.ma.masked_invalid(out_image)
NDVI_mean_inside = np.mean(masked_array)
# 3.2 提取多边形外部的 NDVI 值
out_image_outside, out_transform_outside = rasterio.mask.mask(src, shapes, crop=True, invert=True)
# 清理无效值（例如 NaN 或 NoData 值）
masked_array_outside = np.ma.masked_invalid(out_image_outside)
NDVI_mean_outside = np.mean(masked_array_outside)
# 4. 打印结果
print(f"多边形内部平均 NDVI 值: {NDVI_mean_inside}")
print(f"多边形外部平均 NDVI 值: {NDVI_mean_outside}")代码解释
- 导入必要的库: 导入 rasterio，fiona 和 numpy 库。
- 定义文件路径: 指定 shapefile 和栅格文件的路径。请务必替换为实际的文件路径。
- 加载 shapefile: 使用 fiona.open() 函数打开 shapefile，并提取所有多边形的几何信息。
- 读取栅格数据并提取 NDVI 值:
- 使用 rasterio.open() 函数打开栅格文件。
- 使用 rasterio.mask.mask() 函数裁剪栅格数据，提取多边形内部的 NDVI 值。 crop=True 参数表示裁剪结果仅包含多边形范围内的像素。
- 使用 np.mean() 函数计算裁剪区域的平均 NDVI 值。
- 使用 rasterio.mask.mask() 函数裁剪栅格数据，提取多边形外部的 NDVI 值。invert=True 参数表示裁剪区域为多边形外部。
- 使用 np.mean() 函数计算裁剪区域的平均 NDVI 值。
- 打印结果: 将计算得到的平均 NDVI 值打印到控制台。
注意事项
- 确保 shapefile 和栅格文件位于指定的路径，并且具有正确的格式。
- 如果栅格数据包含无效值（例如 NaN 或 NoData 值），请在使用 np.mean() 函数之前使用 np.ma.masked_invalid() 函数进行屏蔽，避免影响计算结果。
- rasterio.mask.mask() 函数返回的是一个 numpy 数组，你可以根据需要对其进行进一步的处理和分析。
- 如果 shapefile 包含多个多边形，代码将提取所有多边形内部和外部的 NDVI 值，并计算平均值。如果需要分别提取每个多边形的 NDVI 值，请修改代码，循环处理每个多边形。
- rasterio.mask 允许设置 all_touched=True 参数，以便包含与多边形边界相交的像素。默认情况下，仅包含完全位于多边形内的像素。
总结
本文介绍了使用 Python 和 rasterio、fiona 库提取栅格数据中多边形内部和外部平均 NDVI 值的方法。 通过加载 shapefile，裁剪栅格数据，并计算裁剪区域的平均值，您可以轻松地提取和分析特定区域的 NDVI 值。 这种方法可以应用于各种遥感应用，例如土地覆盖分类，植被监测和环境评估。
Postman是一套面向API全生命周期的开发与测试平台，可在Mac上发送和调试多种协议请求，管理集合、环境及认证信息，并完成脚本测试、模拟服务、监控、文档发布和团队协作。它适合开发、测试及运维人员使用，提供个人免费方案。
Docker Desktop 是 Docker 官方为 Mac 提供的容器开发环境，集成 Docker Engine、命令行工具、Compose、Build、Kubernetes 与图形化管理界面。
Android Studio 是 Google 面向安卓应用开发提供的官方集成开发环境，基于 IntelliJ IDEA 构建。它整合智能代码编辑、Gradle 构建、界面预览、设备模拟、性能分析、调试测试及 APK 分析等工具。
Xcode 是苹果面向 Mac 提供的官方集成开发环境，将代码编辑、项目构建、SwiftUI 实时预览、设备模拟、自动化测试、断点调试、性能分析和应用分发整合在同一套工作流中，适合开发与维护 iPhone、iPad、Mac、Apple Wa
JetBrains WebStorm 是面向 JavaScript、TypeScript 和现代 Web 技术的专业集成开发环境，可智能理解项目结构，提供代码补全、重构、错误检查、运行调试、测试及 Git 工具。
PhpStorm 是 JetBrains 推出的专业 PHP 集成开发环境，可在 Mac 上完成代码编写、智能检查、重构、调试、测试和数据库管理。它支持主流 PHP 框架、Composer、Git、Docker 与远程解释器。
JetBrains PyCharm 是面向 Python 开发者的集成开发环境，提供智能补全、代码检查、重构、调试、测试、版本控制及虚拟环境管理功能。Mac 版兼容 Intel 与 Apple 芯片，核心功能可免费使用。
IntelliJ IDEA 是 JetBrains 推出的专业集成开发环境，重点服务于 Java、Kotlin 与 JVM 项目开发。软件集成代码补全、静态检查、重构、调试、测试、版本控制及构建工具，并可借助插件扩展语言和框架支持。
Visual Studio Code 是微软提供的免费代码编辑器，支持 macOS，并为 Intel 与 Apple 芯片 Mac 提供适配版本。软件集代码编辑、调试、Git 源代码管理、集成终端和扩展市场于一体。
REAPER是Cockos开发的数字音频工作站，提供多轨音频与MIDI录制、剪辑、处理、混音和母带制作工具。Mac版兼容Intel与Apple芯片，支持AU、VST、VST3、CLAP等插件格式，并提供高度可定制的工作流程。
Adobe Audition 是面向播客、影视后期、音乐制作及内容创作者的专业音频工作站，提供波形精修、多轨录制与混音、频谱编辑、降噪修复、响度匹配和效果处理工具，并可与 Adobe Premiere 协同完成视频声音制作。
Pro Tools 是 Avid 推出的专业数字音频工作站，覆盖多轨录音、无损剪辑、MIDI 编曲、插件处理、自动化混音及声音与画面同步等工作。它拥有 Intro、Artist、Studio 和 Ultimate 等方案。
PixPix 是面向品牌、卖家和内容团队的在线 AI 图片视频创作平台，提供图视频 Agent、无限画布、图片生成、视频生成、自然语言编辑及商品套图、产品精修和带货视频等电商工具。
Postman 是面向 API 全生命周期的开发与协作平台。Windows 桌面版可用于构建和发送 HTTP、GraphQL、gRPC、WebSocket、MQTT 等请求，并提供环境变量、集合、脚本测试、模拟服务、监控、文档及团队工作区等
Docker Desktop 为 Windows 开发者提供集成式容器环境，内置 Docker Engine、命令行工具、Compose、Build 与 Kubernetes 等组件，可在本机完成镜像构建、容器运行、日志查看和资源管理。
Android Studio 是 Google 面向 Android 应用开发提供的官方集成开发环境，以 IntelliJ IDEA 为基础，集成智能代码编辑、Gradle 构建、设备模拟、界面预览、调试测试及性能分析工具，
WebStorm 是 JetBrains 面向 JavaScript、TypeScript 和现代 Web 技术打造的集成开发环境。它提供项目级代码分析、智能补全、安全重构、导航搜索，并集成调试、测试、终端、Git 与常用前端工具，
PhpStorm 是 JetBrains 面向 PHP 与现代 Web 开发打造的专业集成开发环境，提供代码补全、静态检查、重构、调试、测试、数据库、版本控制及前端语言支持。Windows 用户可直接下载安装程序，
PyCharm 是 JetBrains 推出的 Python 集成开发环境，提供智能代码补全、错误检查、重构、调试、测试、Git、终端及 Jupyter Notebook 支持。统一版安装后可体验 30 天 Pro 功能，
JetBrains IntelliJ IDEA 是以 Java、Kotlin 和 JVM 项目为核心的专业集成开发环境，提供智能补全、代码检查、安全重构、调试、测试、版本控制及构建工具集成。统一发行版的核心功能可免费用于商业和非商业项目，
Visual Studio Code 是微软推出的免费代码编辑器，Windows 版集成智能补全、代码导航、调试、终端和版本控制，并可通过扩展添加语言、格式化及远程开发能力。它启动较快、配置灵活，
REAPER 是 Cockos 推出的数字音频工作站，适用于录音、MIDI 编曲、音频编辑、混音、母带制作及声音设计。软件安装包小巧，支持多轨录制、灵活路由、自动化、第三方插件和脚本扩展，并提供完整功能的 60 天免费评估期。
Adobe Audition 是面向播客、影视后期、广播与音乐制作的专业音频工作站，提供波形编辑、多轨混音、频谱分析、录音、降噪修复和效果处理等功能。Windows 版通过 Adobe Creative Cloud 安装和管理，
Pro Tools 是 Avid 推出的专业数字音频工作站，提供多轨录音、音频剪辑、MIDI编曲、插件处理、自动化混音及影视声音制作能力。它适合音乐人、录音师、制作人和音频后期团队使用，并提供免费入门版、订阅方案及永久授权选择。