我需要在循环中堆叠一些栅格,如下所示:
for(month in 1:12){
.
.
.
"some algorithm spiting out a raster called 'sm_esa'"
sm_esa_stack<-stack(sm_esa)
}
最后,我想创建一个有12层(每层一个月)的堆栈。但我的最后一行显然会覆盖每一个新的栅格,而不是堆叠起来。有什么提示吗?
有没有办法在R中的时间序列NDVI图像上运行savitzky golay滤波器?我已经尝试过使用‘signal’包中给出的以下代码;sg <- sgolayfilt(timeseries,3,5)。但是它返回下面的错误;
Error in if (all(is.na(x))) return(x) :
argument is not interpretable as logical
这里的文件"timeseries“是一个堆叠的光栅NDVI图像。有人能在这方面帮我吗?谢谢你的好意帮助。
我有一个被河网屏蔽的二进制形式的栅格层。同时计算了N0。在R(频率(栅格))和QGIS中使用r.report的像素,我在两者中找到了相同的数字。但在面积计算中,以平方为单位。km我发现(tapply( area (raster),raster[],sum))和QGIS计算的面积R有差异。然而,我遇到的主要问题是,为什么面积计算与像素数量不一致?栅格的分辨率为30秒(约1 km *1 km),因此像素的数量必须大约等于以平方千米为单位的面积。栅格具有地理坐标系OGC:CRS84 -WGS84 (CRS84) -地理,并且为.grd形式。我还将其投影到用于QGIS的UTM,这会略微增加面积,但不会
我有一个为我的任务定义正确的函数的问题。为了过滤图像,我设置了一个过滤器矩阵eucdis
library(rgdal)
library(raster)
refm=matrix(1,nrow=11,ncol=11)
M = dim(refm)[1]
N = dim(refm)[2]
eucdis = matrix(NaN, nrow=11, ncol=11)
for (i in -5:5){
for (j in -5:5){
eucdis[i+6,j+6] = 2*(sqrt(sum(abs(0-i)^2+abs(0-j)^2))) #euclidean
我花了很多时间寻找这个问题的答案,但什么也找不到,所以我在这里
我有一个完整的ascii光栅文件对应的气温和露点温度的某个地区超过744个小时的时间步长。(所以我有744个空气温度和744个露点文件,对应31天的一个月)。每个文件只有大约45 kB。
我想把它们堆叠在一起,这样我就可以对它们进行一些分析,我还想把它们的单位从K转换为deg F。
文件名为air Tair1.txt、Tair2.txt、...Tair744.txt和Eair1.txt,Eair2.txt,...Eair744.txt.
使用raster包,我可以轻松地将所有文件加载为光栅:
for (i in 1:744) {
如果以前有人问过这个问题,我很抱歉,但我无法找到解决我的问题的方法,我相信这是一个相当简单的问题。我有一个单一的源栅格数据集,其中包含100到500之间的连续浮点值。我想要做的是以50为增量循环遍历此源栅格,以导出/创建小于增量的所有值的新栅格数据集。例如,我有以下R代码(使用光栅库)来指定光栅和识别增量。我想开发一种方法来自动创建9个输出栅格数据集,这些数据集小于或等于每个增量的值。我好像到不了那里。有人能帮上忙吗?蒂娅!
#Trying to iteratively create new raster datasets
#Based on increments of Source Ras
我正在处理R (raster软件包)中的MODIS栅格数据,我想知道为什么当将两个栅格文件相乘时,栅格的大小会变得7000倍。在我的例子中,一个栅格是一个简单的MODIS栅格(介于0和255 (整数)之间的值),另一个是一种模板(0和1(整数))。所以我不明白为什么栅格要这么大?当尝试用writeRaster保存栅格时,这是一个很大的问题,因为我的R会话总是崩溃。
我的代码如下所示:
a <- list.files('All_MODIS_files_of_one_year.tif')
for (day in 1:length(a)){
ra <- raster(a[
当前,对于同名文件,当您按顺序上传它们时,它们会被堆叠在另一个文件上。我们能控制基于元数据属性的版本堆叠吗?
因此,下面这样的序列将创建同一个文档的3个版本
filename "test.txt" with metadata attribute "x" - attempt 1 - creates version 1
filename "test.txt" with metadata attribute "x" attempt 2 - creates version 2 of above
filename "test.t
我试图使用来自zonal_stats包的函数rasterstats从.shp文件中的每个形状的.tif文件中获取光栅统计信息。我设法在QGIS中做到了这一点,没有任何问题,但我必须对200多个文件进行同样的操作,这将需要很长时间,所以我正在尝试Python方法。文件和复制代码都在我的中。
我的剧本是:
import rasterio
import geopandas as gpd
import numpy as np
from rasterio.plot import show
from rasterstats import zonal_stats
from rasterio.transfor
当在datashader中绘制一组数据时,如果X轴有离散的数字和欠采样,就会在可以看到背景的列之间留下空白。
我一直试图通过设置一个更大的点大小或者使用动态传递函数来解决这个问题。没有运气-很可能我只是不知道正确的方式应用这些。
下面是重现我的意思的示例代码:
import pandas as pd
import numpy as np
import datashader as ds, colorcet
import holoviews as hv
from holoviews.operation.datashader import datashade
from holoviews impo