掩码并从vrt文件中提取单元格值？

作者: 不丶离
发布时间: 2024-07-02 11:46:11 (10天前)
转自：

2 条回复

0#
回复此人
敢嫁就敢娶 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 数据似乎在Pseudo Mercator投影（EPSG 3857）中。因此，您应该指定范围 <code> projWin </code> 在那个坐标系中，或者添加 <code> projWinSRS </code> 如果你想在不同的坐标系中提供它们。 </p> <P> 另外，如果你想 <code> gdal.Translate </code> 要输出到VRT文件，您应该添加 <code> format='VRT </code> 。因为在你的代码片段输出到默认文件格式，这是GeoTIFF。 </p> <P> 当我假设您的坐标是WGS84（EPSG 4326）时，它定义了斯里兰卡南部海洋上的一个小区域。鉴于数据的性质，这没有多大意义。 </p> <P> 如果要读取坐标给出的数组，可以使用： </p> <pre> <code> invrt = 'GHS_BUILT_LDSMT_GLOBE_R2015B_3857_38_v1_0.vrt' outfile = '/vsimem/tmpfile' ds = gdal.Translate(outfile, invrt, projWin=[80.439, 5.341, 81.048, 4.686], projWinSRS='EPSG:4326') data = ds.ReadAsArray() ds = None gdal.Unlink(outfile) </code> </pre> <P> 绘制的数组看起来像： <a href="https://i.stack.imgur.com/l9tcA.png" rel="nofollow noreferrer"> <img src =“https://i.stack.imgur.com/l9tcA.png”alt =“在此处输入图片说明”/> </A> </p> </DIV>

编辑

登录后才能参与评论