5
\ begingroup美元

有人知道海平面以下地区的现有数据集吗?矢量区段是理想的。出于我的目的,不需要非常详细,但如果它包括所有领域就好了(例如,高估比低估要好)。

有一个维基百科上海拔低于海平面的陆地地点列表这看起来是相当完整的,但要找到每个地点的海平面以下区域的范围不太可能是简单的。

\ endgroup美元
2
  • \ begingroup美元 在你的链接页面上,你可以找到文章中所有地方的来源。我不认为在座的任何人能对此给出更好或更详细的回答。 \ endgroup美元
    - - - - - -trond汉森
    2020年1月16日6:54
  • \ begingroup美元 我一直在用我正在做的一个朴素的地形渲染器涉水SRTM数据。SRTM覆盖纬度+60°/-56°。然而,将有几个10gb的梳理(3弧秒分辨率)。从栅格数据(SRTM或任何其他数据)沿0仰角生成多边形不是一个大问题,这取决于您的GIS系统(或一点编程练习)。困难可能在于地理参照,如果它在不同的来源上不是统一的。 \ endgroup美元
    - - - - - -user18607
    2020年1月16日11:29

1回答1

8
\ begingroup美元

您可以尝试使用全局数字高程模型构建一个。有几个免费的图像中SRTM,或ASTER GDEM.你必须寻找所有包含负值的像素。相邻的负像素会给你海平面以下的连接区域,你可以很容易地估计它们的面积,只需用一个区域的像素数乘以一个像素的面积。

我不是说这很容易。主要的挑战可能是一些区域将被分割成几个块。


这是最新情况。我花了几个小时做练习。我从ASTER全球DEM中下载了三个以Badwater Basin(美国死亡谷)为中心的瓷砖,这是维基百科列表中的一个地区。我写了一个小的(约20行)Matlab脚本,它做以下工作:

  • 二值化图像,使海平面以上的每个像素设置为0(黑色),海平面以下的每个像素设置为1(白色)。
  • 三个数组的连接,以便在下一步中只有一个图像(我手动添加了第四个,为了保持尺寸一致性,全黑色数组)。
  • 对于二值图像,Matlab中有一个有用的函数叫做regionprops它查找所有连接的区域并返回它们的属性,例如它们的面积(以像素数表示)。

结果如下:Badland_Basin_black_and_white我对它非常满意,因为它看起来非常类似于标记为“海平面以下”的黄色区域这张美国地质勘探局地图

在这里输入图像描述

现在,它如何转换成面积呢?根据Matlab,白色区域的面积为1 236 989像素。每个像素约为30x30米,海平面以下面积为1113公里$ ^ 2美元.我在网上唯一能找到的数据是1425公里$ ^ 2美元来自大英百科全书,没有来源。但对我来说,这已经足够接近我的估计了!

一些缺点:

  • 我清楚地知道哪块瓷砖包含了这个特定的海平面以下区域。如果要遍历所有的瓦片(几乎有23000个!),您必须找到一种方法来检测一个区域何时被分割到多个瓦片,以便通过连接重建该区域。看起来没那么难,但也不是小事。
  • 我没有考虑到相邻瓷砖之间有一个像素重叠的事实(它不应该改变结果)。
  • 我没有考虑像素大小实际变化的事实(分辨率为一弧秒,在赤道约为30米,但随纬度变化)。
\ endgroup美元
2
  • \ begingroup美元 是的,这看起来相当直截了当,结合上面的wiki列表,它可能是全球可行的。我想我可能不得不这样做,但我觉得值得先问问是否有什么东西已经存在了。 \ endgroup美元
    - - - - - -naught101
    2020年1月17日23:32
  • 1
    \ begingroup美元 如果你最终这样做了,请在这里分享它作为一个自我回答!:) \ endgroup美元
    - - - - - -让-玛丽•Prival
    2020年1月18日9:57

你的答案

点击“张贴您的答案”,即表示您同意我们的服务条款隐私政策而且饼干的政策

这不是你想要的答案?浏览带标签的其他问题问自己的问题