如何使用R计算月度NDVI栅格

文章/答案/技术大牛

发布

1回答

r、raster

NDVI RasterStack 我有一个为期15天的1981-2015 NDVI RasterStack。我需要用15d的数据计算每月的NDVI。我想知道如何将相同名称的MM栅格的平均值计算到新的1981-2015月度栅格堆栈中感谢您的帮助！非常感谢。 names XYYYY.MM.DD

浏览 10提问于2020-10-24得票数 0

1回答

我使用raster读取栅格( GeoTIFF)，然后将其转换为表，这样我就可以进行一些分析，并对一些值进行重新分类。如下所示： r_ndvi <- raster(paste0(path, "ndvi.tif")) #read in as raster rdf_ndvi <- as.data.frame(r_ndvi)#turn tif into df of DN 完成这些计算后，我希望将表转换回光<em

浏览 127提问于2021-09-27得票数 0

回答已采纳

1回答

叠加()函数中的错误，用光栅和R中的数值计算

r、loops、overlay、raster

我想用一张栅格表做一些计算。因此，我写了一个循环，到目前为止，我只使用光栅就可以正常工作了。现在我用栅格和数字vavlues进行计算，这会导致叠加函数中的一个错误。我的代码是：NDVI <- list.files(path = "C:/Users/Felix/Desktop/Bachelorarbeit/Daten/Datenverarbeitung_R/L8_N

浏览 0提问于2018-09-18得票数 0

回答已采纳

1回答

是否有一种方法可以将numpy.where()用于以NaN值作为无数据的光栅数据？

python、python-3.x、numpy、nan、rasterio

我有一个栅格数据，其中包括NaN值作为无数据。我想从它中计算新的栅格，比如如果raster==0做statement1，如果raster==1做statement2，如果光栅在0和1之间做statement3，否则不会更改值。如何使用numpy.where()函数完成此操作？/NDVI_output/NDVI.tiff") min_value = ndvi_s = np.nanmin(ndvi) #<e

浏览 3提问于2020-04-26得票数 2

回答已采纳

1回答

在R中创建一个循环，根据存储文件的变量的名称来命名(保存)文件

r、gis、raster

我的代码非常简单，有一堆计算得到的栅格文件存储在变量:NDWI，NDWI，VSI等，每个计算都是基于一个卫星图像，其中有一个名称与日期在它。我会找到一种方法来提取日期并将其存储在一个变量中。本质上，我想要的是我的代码的一部分，它遍历所有创建的文件，将它们保存为我笔记本电脑上特定路径下的栅格，格式为(“变量名”_“日期”.tif。)(NDVI+NDWII) writeRaster(NDWI, "C:/Users/Perrin/Desktop/INRA&#x

浏览 70提问于2021-09-15得票数 2

回答已采纳

1回答

如何为每个栅格单元选择一个点？

r、sample、r-raster

我有一个高度聚集的点shapefile ("search_effort.shp")和一个NDVI光栅(分辨率m: 30.94948，30.77829)。我想通过选择每个栅格单元格一个点来设置我的形状，并创建一个新的search_effort search_effort.shp文件。我使用的是R版本4.0.3 我想我本可以使用‘gridsample’包(在'raster‘v1.3-1中)，但它已经从CRAN存储库中删除了，我不希望使用归档版本。在R</

浏览 13提问于2020-11-05得票数 0

回答已采纳

1回答

如何计算R在光栅层中相同纬度上的pixiel值？

r、raster

我有一个NDVI(Normalized植被指数)光栅对象，现在我想在相同的latitude.Like上计算像素值，我想知道30°N latittude.How中的像素值，我能用R做这个吗？但我发现在“栅格”或"terra“包中没有任何功能可以做到这一点。

浏览 5提问于2022-11-08得票数 0

回答已采纳

2回答

terra extract()给出R中的NA值

r、terra

我正在尝试使用terra::extract()来获取地块中的平均NDVI值。我注意到来自terra::extract()的数据帧给出了NA值(实际上，在4519个观察值中，我有2494个NA值)。我的NDVI光栅的空间分辨率约为0.5英尺(0.15米)。下面是我的代码。我想知道使用哪个更好，exactextractr::exact_extract()还是terra::extract()。任何帮助/指导都是很棒的！)) epa_ndvi_sr<- terra::

浏览 1提问于2021-03-01得票数 0

1回答

基于R的GIMMS 3g.v1处理

r、date、process、stack、r-raster

我想用R下载和处理GIMMS NDVI 3g.v1集合，使用gimms包，我成功地下载了完整的集合，但是我不知道如何在文件上分配日期。我的最终结果是，我希望成为一个具有良好质量(零标志值)的NDVI的栅格堆栈。因此，每月2次，每年。我真的很感谢你的帮助。my_dir <- "my_directory_here"my_ndvi <- sta

浏览 0提问于2019-07-01得票数 1

回答已采纳

1回答

如何用不同程度的光栅制作光栅砖？

r、stack、crop、raster、resampling

我是R的新手，所以这个问题是非常基本的，但是我一直在与它斗争，无法找到一个有效的解决方案。我想从同一地区的一些陆地卫星图像中创建一个栅格砖。它们是以HDF-EOS格式下载的，我使用Modis重新投影工具将它们转换为.tif。重新对栅格进行采样，给它们提供相同数量的列和行。理想情况下，这将确保栅格细胞是对齐的，可以放在

浏览 2提问于2013-12-22得票数 9

回答已采纳

1回答

R中的栅格重分类

r、raster、r-raster

我有一个光栅图像，它包含NDVI的值在-0.2到0.6之间。NDVI #dimensions : 706, 953, 672818 (nrow, ncol, ncell)其中lsat8_pv是另一个光栅，它是以下面的格式计算<

浏览 1提问于2020-06-09得票数 0

回答已采纳

1回答

Python多处理脚本无限期挂起

python、image、multiprocessing

以下脚本1)将图像波段转换为numpy数组，2)计算归一化差分植被指数(NDVI)，3)将numpy数组转换回光栅并写入磁盘。该脚本旨在了解使用多处理提高速度的想法。第一部分通过迭代工作空间，在每个栅格上执行处理，并写入磁盘(总时间=2分钟)来正常工作。第二个多处理部分无限期地“挂起”，并且不产生任何输出。在脚本的多进程部分，我哪里出了问题？import arcpy, os, time arcpy.env.workspace = r&#x

浏览 1提问于2015-04-18得票数 1

1回答

如何从关联栅格中查找p值

r、geospatial、correlation

我想计算p值并将其与相关光栅图合并，但我不知道如何计算光栅p值并将其与相关光栅图进行边缘处理。我的R编程代码如下(y13,y14,y15,y16,y17,y18,y19,y20,y21,y22,y23,y24,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12)z2 <- r.stack[[13:24]] z3 <- stack(z1,z

浏览 1提问于2016-05-30得票数 0

1回答

如何从R中的栅格堆栈中提取时间序列

r、time-series、r-raster

我已经从MODIS数据的NDVI层创建了RasterStack。现在我想从这些数据的不同位置提取时间序列数据，这样我就可以使用BFAST/greenbrown包来估计趋势和断点。2018.12.27",# , SDSstring = "1000000000000000000000") NDVI_files_path <- "/med

浏览 11提问于2019-07-07得票数 0

5回答

栅格包包含所有硬盘驱动器

r、raster

我正在处理栅格(modis ndvi图像)的时间序列，以计算该序列的平均值和st.deviation。每个年度系列由23个ndvi.tif图像组成，每个图像大小为508Mb，因此总共需要处理11 to的图像。以下是一年的脚本。我必须在几年内重复这个过程。library(raster)filesndvi <- list.files(, pattern="NDVI.tif",full.names=TRUE) filesetnd

浏览 0提问于2014-08-21得票数 11

1回答

使用过多内存从{raster}包中提取

r、memory、extract、raster

我一直在使用raster包中的extract函数，使用shapefiles定义的区域从光栅文件中提取数据。但是，我遇到了这个进程现在所需的内存量问题。栅格文件较大(约1.6 The )shp <- mclapply(list.files(pattern="*.shp",full.names=TRUE), readShapePoly,mc.cores=6)mc<- functio

浏览 2提问于2013-03-29得票数 3

回答已采纳

2回答

如何使用空间掩码限制光栅处理范围？

image、matlab、image-processing、shapefile、arcpy

下面是我正在使用的一些(可重复的)样本数据： file = 'C:\path\to\doi1m2011_41111h4nw_usda.tif';outputdir = 'C:\output\' NIR = im2single, R, 'GeoKeyDirecto

浏览 9提问于2014-02-24得票数 6

回答已采纳

1回答

提高气候数据大栅格叠加计算速度

我有一个很大的栅格堆栈(4396663296个元素，13.3Mb)，其中包含1920年至2015年阿拉斯加和加拿大的月度降水量数据，缩小到2公里。我想计算观察期(1920-2015)的年平均降水量，但我发现我的代码在R中停滞不前。有没有一种更快更有效的方式来执行下面的代码？library(raster) r <- stack(list.files(pattern="tif$", full.names=F

浏览 11提问于2018-09-01得票数 0

1回答

如何应用从[-1，1]到[0,255]的线性拉伸？

image、matlab、image-processing、stretch

我有一个脚本，它从4波段的地理信息中计算出，并将NDVI写成一个地理信息。我正在研究1米空间分辨率的航空图像。输出NDVI测井值范围为-1 ~1，为32位浮点.如何应用线性拉伸，以便在中输出范围为0- 255。file = 'D:\create_ndvi\3711101_se1.tif';outputdir = 'D:\create_ndvi\' NIR= im2singl

浏览 4提问于2014-02-24得票数 0

回答已采纳

2回答

使用像素值和种群大小分割栅格图像

r、raster

我通常使用ENVI来完成这一过程，但是我正在尝试将所有内容都保留在R中，以编写更多的栅格处理脚本。在ENVI中，我使用的工具称为“分割图像”。这些参数包括像素阈值、最小种群和考虑的邻居数量。有没有办法在R中做到这一点？这似乎是一个相对简单的分析，但我似乎找不到一个可以做到这一点的函数。谢谢你！

浏览 22提问于2019-02-19得票数 1

回答已采纳

点击加载更多