地球自然特性在1:10 million-scale 5类型的光栅文件来满足您的带宽和重点内容。该版本是10 million-scale栅格数据提供:高分辨率文件在21600 x 10800像素和低分辨率为16200 x 8100。10米的光栅文件登记准确矢量数据。嵌入式光栅内容包括:土地覆盖、阴影救济,海水和流域湖泊。
在GIS中,栅格属性里有关于栅格自身的信息,背景(nodata value)对于识别一张图像的边界像元尤为重要,我们目的只要把每行每列中的第一次出现不是nodata的像元和最后一次出现nodata的前一个像元就可以了。
对于栅格,可以用ArcPy中的RasterToNumpyArray函数将将栅格转成numpy数组,然后就可以按照所想读取出每行列中首尾像元。
以下是部分代码提取边界像元的核心算法,其实是很简单的一个思路(假设0是nodata value)。
a=[[0 for col in r