使用Planet Explorer QGIS插件查看洪水事件


完成时间:根据您对QGIS和卫星图像的熟悉程度,本教程预计需要30分钟到2小时不等。188bet手机滚球

它安装了时间

首先是首先,您需要购买一个星球订阅这里或注册14天免费试用这里。可以安装行星浏览器插件,但未使用,直到您使用Planet帐户登录。此外,某些功能,例如图像下载可能会限制试用帐户。

其次,如果您还没有QGIS,请通过下载并安装它下载QGIS.。您很可能希望使用QGIS独立安装程序而不是OSGeo4W网络安装程序。Planet的QGIS插件旨在用于3.6+的QGIS版本。请记住,您可以在同一台计算机上同时安装不同版本的不同版本,因此如果您尚未升级,则最好获取最新版本,而不必担心破坏旧项目或需要卸载旧版本。

安装后,启动QGIS桌面(您可以选择用草版本或没有)。在QGIS菜单栏中,选择插件下拉,然后选择全部如果它尚未选择,并在“行星浏览器”中键入搜索…盒子。在列表中选择它并按下安装插件。将添加一个新的按钮在工具栏中。您的QGIS可能看起来与下面的一个不同,具体取决于版本,哪个面板和工具栏是打开的,但是一般过程应该是相同的。

注意:如果是输入凭据弹出窗口,输入您的行星登录密码两次,然后按好的

请按行星探险家现在显示在QGIS中的按钮工具栏部分。请注意,您的按钮可能会在不同的地方出现,具体取决于您过去的工具栏配置。

按下按钮后Planet Explorer面板将出现,位于QGIS窗口的右侧。请注意,如果您已经在QGIS窗口的右侧打开了一个或多个面板,则从以前的QGIS会话,Planet Explorer面板将出现在右下角,在它们下面。

使用电子邮件地址密码输入框用于输入您的凭据(如果您已经拥有它们)。如果你没有,但想注册免费试用14天,请使用报名!链接,注册,然后返回此处并输入您的凭据。经过身份验证,Planet Explorer面板会是这样的:

如果你曾经关闭Planet Explorer面板并希望回到它,只需点击看法在菜单栏中,然后面板,然后勾选旁边的复选框行星探险家

这是疯狂的快点

有时人们想要跳进。如果这听起来像你,请使用本节作为30秒的快速启动,然后继续进行教程。否则,跳到教程。

  1. 添加数据集到QGIS
  2. 点击范围>当前可见程度在里面Planet Explorer面板
  3. 点击搜索
  4. 收购日期,检查其中一个复选框。在它旁边,按森伯斯特设置按钮选择添加预览图层到地图

这是教学时间

2019年11月初,佛蒙特州密斯quoi河发生洪水。让我们通过探索一个受影响的地方来看看那里发生了什么:佛蒙特州的里奇福德。我们将看到发生在2019年11188bet手机滚球月1日前后的洪水之前和期间的卫星图像。我们从开放地图中提取了佛蒙特州北部的一些村庄点并将它们准备在一个名为geojson的文件中村庄点.Geojson.,位于此目录中。下载该文件并将其拖放到QGIS图层面板

使用>画矩形按钮Planet Explorer面板在代表Richford的点绘制一个粗糙的矩形,这是该数据集中的最北端。为此,请确保在单击并将矩形拖动到QGIS Map窗口时按住鼠标按钮。注意:我们使用QGIS的QuickMapservices插件将底图添加到项目中。这是可选的。

地图自动放大到该位置,并在其周围放置一条蓝色虚线。这表示搜索区域,也称为兴趣区域(AOI)。接下来,因为我们需要来自特定时间段的图像,所以使用日期范围(UTC)输入框提供两个日期:2019/10/23和2019/11/03。这些是在过滤器面板的一部分。你可能需要滚动过滤器部分看看日期范围(UTC)输入部分,具体取决于您的屏幕空间有多少。然后单击搜索

一些PlanetScope 4波段场景出现在收购日期部分。将鼠标悬停在这些场景列表中,并注意到地图上出现蓝线。这些显示您在该日期的卫星捕获的图像的位置。根据搜索矩形的初始大小,您可能希望缩小一下以便查看整个行星剖面图像占地面积。

每个场景旁边的图像缩略图还将为您提供场景中云盖的一般数量的感觉。虽然您的可用场景列表可能看起来不同,但根据您画的矩形的大小,在我们的情况下,我们列表中的第一个和第三场景看起来很混乱。

要在大地图视图中看到两个不那么多云的图像,请单击每个右侧的Blue Sunburst设置按钮,然后单击添加预览图层到地图。请密切关注以下屏幕截图,以便您可以尝试向您的地图添加相同的两个场景,因为您的列表可能与我们的名单不同。具体来说,我们正在加入2019年11月24日和2019年10月24日的场景。

注意QGIS图层面板现在列出这两个场景作为图像预览。你可以在图层面板然后平移和缩放地图来检查它们。

2019年11月2日,图片预览打开了地图,但2019年10月24日的图片预览关闭了。你现在可以在地图上看到“万圣节”风暴期间的一些洪水:

2019年11月02日图像相当明确。但是,如果您点击10月24日图像,您可能会注意到它看起来过混浊。自从2014年10月24日的场景将成为我们的图片,它可能是可以扩大我们的日期搜索稍微早一点的时间框架,看看我们是否可以得到一个更好的图像。

改变第一个日期过滤器第2019/10/10/20节并按搜索再次。使用新的搜索结果将在结果区域中出现一个新选项卡。将2019年10月21日的场景添加到地图上的预览,就像我们使用其他两个图像一样。这个主要是无云,并将作为我们的工作图像。关闭第一个结果选项卡,以便只需开放一个结果“选项卡:

它现在是时候下订单来获得图像,因为我们只是在现在看预览。订购图像将可用于下载,我们可以在地理灯接地中运行更高级的分析。展开11月2日场景清单,以便看到我们可以订购的三个图像。使用“100A”后缀单击图像旁边的复选框。

展开10月21日场景列表,并使用“0F32”后缀展开图像列表。单击列表旁边的复选框20191021_134736_0F32。现在点击这一点命令按钮打开订单对话框。为订单添加名称并检查剪辑项目到AOI将它们自动剪辑到我们的位置。点击下订单

请注意,对话框的底部现在显示订单请求。点击提供以打开的链接订单监控对话框

这将显示您之前订单的状态以及此新订单。你可以随时到达订单监控对话框同样,如果您关闭它,通过用户菜单右下角的按钮行星探险家面板所示:

订单状态完成后下载当你按下按钮时,按钮将被激活刷新

您还将通过电子邮件通知订单已准备就绪,如果您愿意,您可以直接从提供的电子邮件链接下载。

如果我们早期探索的预览场景仍然在您的QGIS项目中,通过右键单击它们并选择删除它们去掉。这将使项目更容易工作。

下载zip文件订单监控将两个最大的<.tif>文件拖放到QGIS图层面板。注意,解压缩下载后,数据将位于名为files的子文件夹中。这些文件被称为:

20191021 _134736_0f32_3b_analyticms_clip.tif 20191102 _152019_100a_3b_analyticms_clip.tif

请注意,第一个是10月21日日期的场景,从文件名开始时显而易见,第二个是11月2日日期的场景。

这些文件中的每一个都有4个频段。您可以看到每个像素通过使用来包含四个值识别特征按钮查看其中一些值。PlanetScope波段排序:蓝色,绿色,红色,近红外。记住这一点识别特征工具工作,你必须有一个层选择在面板(或者你可以右键单击地图,选择你正在识别的层)。

我们希望利用归一化差异水指数(NDWI)来描绘两种场景之间的水的变化程度,其用带2(对应于绿色)和带4(对应于近红外)来计算。NWDI的公式是:

NDWI =(绿色-近红外)/(绿色+近红外)

对于PlanetScope 4波段图像,公式可以表示为:

NDWI =(带2-带4)/(频段2 +带4)

(在本教程中,我们使用绿色和NIR频带方法而不是近红外和短波红外方法计算的两种方法。)我们希望为我们的两个图像中的每一个计算NDWI。在QGIS中,点击栅格在菜单栏中选择光栅计算器......所有光栅乐队都将显示在光栅的乐队部分。使用公式,双击相应的乐队并填写10月21日场景的相应数学符号,直到您有一个如此如此。请记住,您可以双击栅格频段光栅的乐队列表填充光栅计算器表达式输入框:

但不要压好的然而!这是棘手的部分:这些行星扫描场景值是辐射值,我们需要将被称为反射值的值。为此,我们必须将我们的数据乘以一定的系数,我们下载数据集的元数据中找到的元数据。返回到10月21日数据集随附的元数据XML文件,并在文本编辑器中打开它(在Windows上,您将右键单击并选择打开>记事本)并滚动到文件的底部。这是将列出场景的系数的位置。以下是如何与适当的系数一起看的屏幕截图。请注意我们想要的系数分别在频带号2和频段编号4部分中列出。

这幅图像的两个系数是:

波段2:0.0000438963652969波段4:0.0000747815337148

将这些乘法器应用于计算中的频段,我们将具有以下表达式,该表达式是格式化的,以便您只需复制/粘贴“如果您想要。(请注意,在某些情况下,复制/粘贴将导致罗斯特计算器看起来奇怪的格式,但请放心,粘贴工作,计算也将仍然工作。)本质上,你将用这个取代之前在计算器里的计算,因为它有我们需要的系数。

(( “20191021_134736_0f32_3B_AnalyticMS_clip @ 2” * 0.0000438963652969)( “20191021_134736_0f32_3B_AnalyticMS_clip @ 4” * 0.0000747815337148))/(( “20191021_134736_0f32_3B_AnalyticMS_clip @ 2” * 0.0000438963652969)+( “20191021_134736_0f32_3B_AnalyticMS_clip @ 4” * 0.0000747815337148))

控件中为输出提供一个名称输出层在右上角的输入框栅格计算器对话框。我的建议是确保名字里有这个场景的日期。然后单击好的。它会对图像中的每个像素进行计算。

为11月2日场景执行整个过程,记住从其元数据文件中获取系数,因为它们将是不同的。

当完成它将看起来像这样,记住你的结果数字将略有不同,因为你的AOI将略有不同。

McFeeters(1996)报告说,当它们代表土壤和植被时,NDWI导致消极或零像素值,但是当它们代表水时,结果将具有正值。在我们的情况下,在目视检查时,数据似乎可以将自身归因于略低较低的阈值-0.1和更高。要与现有文献争论并继续委托的最低错误,我们仍将使用零作为我们的较低阈值,即使可能无法捕获一些水。

因此,我们有兴趣仅在两个结果层中可视化的值,该值大于零。打开层样式面板如果它尚未使用使用看法在QGIS的菜单栏中下拉,然后下拉面板>图层样式。对于这两层,请遵循相同的程序。在里面层样式面板将第一个下拉到刚刚计算的图层中的一个更改为(可能已经是这种情况),然后将第二个下拉缩短到Singleband伪彩色。改变最小值输入框为0.彩色坡道应该设置为布鲁斯如果它尚未。该项目应该看起来像这样:

控件中“0”值旁边的颜色矩形层样式面板,改变不透明度输入框为0%并按好的。这使得所有的零值透明。看看相应场景图层的结果。

为了我们的示范目的,这似乎为我们提供了这种情况的充分和保守的水量近似。做与另一层相同。现在,我们在暴风事件之后,我们可以在暴风雨事件之前对水中的变化进行可衡量的描绘。

引文:

关键词:水体特征,归一化差水指数(NDWI),遥感,DOI: 10.1080/01431169608948714

请注意,在许多情况下,洪水事件与高云事件相对应。因此,不可能对每个洪水事件都执行相同类型的分析。此外,这种分析是不真实的,没有估计的错误被试图。本教程仅用于演示目的,在应用于实际场景之前,必须对这种计算进行适当的使用。