问题标签 [r-modis]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
676 浏览

r - 使用 MODIS 包下载 MODIS hdf

我正在尝试使用 MODIS 包从 LAADS 和 LPDAAC 服务器下载 MODIS 数据。我对这些功能非常缺乏经验,因此不胜感激!

我的代码现在看起来像这样:

当我尝试使用 getHdf 函数从服务器获取 hdf 文件时,我收到以下错误:

0 投票
1 回答
81 浏览

r - MODIS R 包的 getTile("Philippines") 可能存在错误

使用 CRAN 的 R MODIS 包(1.1 版)

代码:

我收到错误:createPolygonsComment(p) 中的错误:rgeos_PolyCreateComment:孤立孔,在索引 1 处找不到包含孔的多边形#

是的,我检查了菲律宾的拼写是否正确。

0 投票
1 回答
212 浏览

r - 有没有办法将 MODIS 瓦片作为 R 中纬度/经度的函数?

我正在寻找一种方法来自动返回 MODIS 瓦片(例如 h12 v4)作为经度和纬度的函数。

我希望 MODIS 包可以使用 getTile(***) 来完成。当我使用 getTile() 时,我看到有一个不错的 GUI,但我正在寻找一种自动方法作为 long/lat 的函数。

谢谢,杰夫

0 投票
1 回答
41 浏览

r - 在 R 中使用条件

我有 3 个栅格,我想在表达式中使用它们,但我可以在 3 个栅格中找到不同的na值。例如:我可以在 2 个栅格中有一个值,但在 3 个中我有na,那么在这种情况下我不能应用我的表达式。

按照我的代码:

我的问题是这部分:!is.na(BSA[,]) & !is.na(WSA[,]) & !is.na(FDS[,])

这部分不起作用。

有人可以帮助我吗?

0 投票
0 回答
189 浏览

r - 基于经度/纬度的 MODIS L3SMI ncdf 子集

我有 ncdf 格式的时间序列 Global MODIS L3SMI 数据集,如何根据经度/纬度范围在特定位置(即印度尼西亚地区)对该数据集进行子集化。然后我想将此数据集写入与原始数据集相同的 ncdf 文件,而不更改 ncdf 结构上的任何内容(即维度和变量名称)

可以从此链接下载示例数据集https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A2018064.L3m_DAY_SST_sst_4km.nc

我尝试使用transNcdfSubsetfrom ncdf.toolspackage 但给我错误消息:

transNcdfSubset("A2018064.L3m_DAY_SST_sst_4km.nc", : 文件中的维度名称不是纬度/经度/时间。重命名!

这是我的代码:

0 投票
0 回答
51 浏览

r - 为什么“url.exist”显示 url 对于 R 中的特定 url 不存在错误?

我正在使用R 包下载 MODIS 数据。

它给出了错误

the http address does not exist! Version may be incorrect OR Server is down!

有趣的是,版本是正确的,我可以超出显示服务器未关闭的数据存储库的站点。

此 R 包试图访问的 URL 是“ https://e4ftl01.cr.usgs.gov/MOTA/MCD43A3.006/

虽然 URL 存在,但当我将其检查为

此错误的原因可能是什么?

0 投票
1 回答
313 浏览

r - 使用 R 重新投影 MODIS 数据(结果为 NA 或没有空间范围)

我将此处存储的 GLASS 反照率数据用于 2000 年前 (AVHRR) 数据,此处用于 2000 年后数据 (MODIS)。我的最终目标是创建每个月的栅格堆栈,其中包含 1982-2015 年的白色天空反照率数据。我遇到的问题是 MODIS 和 AVHRR 数据位于不同的空间参考系统中,我似乎无法将它们重新投影到同一系统中。

我使用 R 将 hdf 转换为 tif,如下所示:

以下是我尝试过的事情:

1.) 使用MODIS 重投影工具。无论出于何种原因,该工具似乎认为 MODIS .hdf 文件的子数据集只是一个图块(最左上角的图块,图块 0,0),而不是全局数据集。我的理解是 MODIS 数据是全局的(不是在瓦片中?),所以我不知道捷运为什么要这样做。

2.)在R中使用光栅包。

这将返回一个值全部为 NA 的栅格:

3.) 在 R 中使用 gdalUtils 包:

这将返回一个基本上没有空间范围的栅格。

关于为什么重新投影此 MODIS 数据如此困难的任何想法?

0 投票
1 回答
145 浏览

r - 将数据 modis EVI 转换为 r 中的日期 yy-mm-dd

我在 2000 年使用 MODIS EVI 栅格。我有 6 个年的栅格,一个月的栅格:

我想像 2000-02-09 那样转换月份,但我不知道该怎么做。

0 投票
0 回答
95 浏览

r - 使用 R 脚本下载 MODIS 图像时出错

我正在尝试使用ModisDownload()R 中 rts 包的功能下载 MODIS 图像。我使用了以下代码。

但每次相同的错误重复:

.modisHTTP(x, v = version, opt = opt) 中的错误:http 地址不存在!版本可能不正确或服务器已关闭!

0 投票
1 回答
1241 浏览

r - 在 r 中使用 HDF 数据

我正在尝试从这里下载并绘制 modis 数据http://orca.science.oregonstate.edu/1080.by.2160.monthly.hdf.vgpm.s.chl.a.sst.php

我经常使用来自海洋颜色网站的 .nc 文件,但我似乎无法以同样的方式使用 hdf 文件。可以将hdf转换为栅格吗?

sh:-c:第 0 行:寻找匹配的 ``' 时出现意外的 EOF
sh:-c:第 1 行:语法错误:文件
NULL的意外结尾

.local(.Object, ...) 中的错误:

.rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", 中的错误:
无法从此文件创建 RasterLayer 对象。(文件不存在)

我已经尝试过这里的建议将 hdf 文件读入 R 并将它们转换为 geoTIFF 栅格

并收到此错误: sh: -c: line 0: unexpected EOF while looking for matching ``' sh: -c: line 1: syntax error: unexpected end of file

我也在这里尝试过将 HDF 转换为地理参考文件(geotiff,shapefile)的建议 ,错误是: .rasterObjectFromFile(x,band = band,objecttype = "RasterLayer", 中的错误:无法从该文件创建 RasterLayer 对象。

我只想提取数据并在我拥有的其他数据旁边的线图中绘制给定位置随时间变化的生产力数据......
谢谢