在R中使用if语句从栅格创建新栅格
的过程如下：
- 首先，确保已经安装并加载了所需的R包，例如raster和rgdal。
- 使用raster函数读取原始栅格数据，并将其存储在一个变量中，例如raster_data。
library(raster)
raster_data <- raster("path_to_raster_file")- 创建一个新的空白栅格，其行列数、范围和分辨率与原始栅格相同。
new_raster <- raster(nrows=ncell(raster_data), ncols=ncell(raster_data),
extent(raster_data), res(raster_data))- 使用if语句遍历原始栅格的每个像素，并根据特定条件设置新栅格的像素值。
for (i in 1:ncell(raster_data)) {
if (raster_data[i] > threshold) {
new_raster[i] <- 1 # 设置新栅格像素值为1
} else {
new_raster[i] <- 0 # 设置新栅格像素值为0
}
}在上述代码中，threshold是一个阈值，用于确定是否将像素值设置为1。
- 可以选择将新栅格保存为栅格文件，以便后续使用。
writeRaster(new_raster, "path_to_output_raster_file", format="GTiff")这样，你就可以使用if语句从栅格创建新栅格了。
请注意，以上答案中没有提及任何特定的云计算品牌商，因为问题与云计算无关。如果你需要了解与云计算相关的内容，请提供相关问题。
相关·内容
2回答
我有一个栅格层，在当前显示SST，但只显示我知道适合物种栖息地的像素。所有其他像素都设置为NA。 现在我有了相同的栅格层，但显示了2050年的预测SST。在此光栅图层中，所有像素都有一个值(陆地表面上的NA除外)。 现在我想做以下事情:我想在第二个栅格层中搜索第一个栅格层中给出的值，并将所有其他值设置为0。因此，结果应该是一个只显示像素的光栅，这些像素的值也可以在第一个光栅中找到。 我认为正确的函数应该是if else function。 有没有办法<
浏览 51提问于2020-04-23得票数 0
回答已采纳
1回答
我需要通过使用z <- rast(x, vals = y)分配值来创建一些新的栅格。我希望新的光栅(z)在与旧光栅(x)相同的位置上有NA值，但是在非NA单元中有新的赋值值(S) (y)。，什么是创建一个新栅格的最有效的方法？
我没有看到任何参数来处理rast()文档中列出的NA值，所以是否需要执行两个步骤，首先将值分配给所有单元格，然后将NAs从旧的<em
浏览 5提问于2022-04-14得票数 0
回答已采纳
2回答
R的栅格功能在磁盘上的内存映射栅格上提供了一个方便的句柄，允许操作非常大的数据集。aggregate函数和包中的其他函数接受一个输入光栅，并创建一个分辨率较小的新光栅。新返回的栅格是驻留在内存中，还是使用新的内存映射创建了磁盘上的临时副本？
浏览 23提问于2019-01-23得票数 1
回答已采纳
我有两个栅格，大小相同，包含来自相同位置的数据，但数据类型不同(一个栅格具有坡度数据，另一个栅格具有纵横比数据)。我希望能够一次查看一个方面的坡度数据，所以我尝试创建一个设置(可能是if/else语句？)其中我说“如果(纵横比条件)在一个栅格中得到满足，那么坡度数据将从另一个栅格中的相同像素中提取。indices of the raster that face north to a new one
nort
浏览 27提问于2019-06-21得票数 1
回答已采纳
我正在使用R，并使用以下库: raster、sp、rgeos、rgdalraster::writeRaster(ras8b,
filename="filepath", format="GTiff" ,但是我使用的代码将栅格保存为64 BitsPerPixel<
浏览 2提问于2017-11-03得票数 0
2回答
条件命令从栅格计算器传输到R
、、、、
到目前为止，我在ArcGIS的光栅计算器中使用了以下语句：这将计算一个新的栅格，该栅格仅包含Land_use等于但我想在R中自动化这一点，因为我必须做很多次。x <- overlay(r2, r1, fun=function(x,y){ y[x!值的栅格，其中<em
浏览 5提问于2016-09-21得票数 3
我已经使用stars包加载了一个栅格，然后在栅格上创建了瓦片。现在，我根据这些瓦片对栅格进行子集。tiles[[i]] <- st_bbox(c(xmin=x0,ymin=y0,xmax=x1,ymax=y1),crs=st_crs(r))
crop_tiles[[i]] <- r[tiles[[i]]] 这里的r是使用read_stars("filename.tif")加载的光栅。现
浏览 13提问于2020-10-16得票数 0
2回答
R:光栅层的覆盖xy坐标
、、、、
. : NA names : janW1 我使用指定的公式计算了拉特和长和弦3.5 1616.5 06 46.99401 -14.90345 5.5 1616.5 0
我如何覆盖或创建一个空间范围为
浏览 2提问于2016-03-15得票数 5
在使用raster软件包时，resample或projectRaster的双线性插值是否总是采用最近的四个质心，或者它是否考虑了单元大小的差异？在的关于重采样中双线性插值的解释中，它说使用了四个最近的质心。在这里的例子中，单元大小是相似的。= 392, replace = TRUE))
此处创建的是要projectRaster到的栅格，然后在projectRaster之后执行<e
浏览 14提问于2018-08-17得票数 3
回答已采纳
1回答
我有几个图层的堆叠栅格数据集，但是，我想计算每个单元格的总和，用于不同的图层选择，并最终生成一个新的图层，谁有一些好的建议，使用计算或叠加或一些其他栅格计算在R？我可以通过循环进行计算，但当我有很多层时，它会消耗很多次，也会使用很多存储空间，我的脚本如下所示：make_calc <- function(rr, start,
浏览 1提问于2018-04-12得票数 0
1回答
我在不同的表格上有几个按时间分隔的栅格。例如，表1表示2016-09-21的光栅，表2表示2016-09-22的光栅，依此类推。我想把这些栅格上的信息汇总一周。我可以使用ST_Value从每个波段获取值，使用avg函数，然后创建一个包含agg信息的新栅格，但这将是大量的工作。
有没有更好的方法呢？
浏览 3提问于2016-09-09得票数 0
3回答
我有两个栅格如下set.seed(2) set.seed(2) 现在我想创建一个基于多种条件的新光栅条件就像当a>0.5时，<e
浏览 6提问于2020-11-10得票数 0
回答已采纳
1回答
我在R中运行物种分布模型，并希望为非洲大陆创建可变栅格，而不是岛屿。我只能找到非洲及其岛屿的shapefiles，而不仅仅是大陆。
1)哪里可以只下载大陆的shapefile？2)如果没有shapefile，我想手动从我的栅格中删除孤岛。有办法做到这一点吗f.e.是否将某些坐标之间的部分栅格设置为NA？
浏览 16提问于2020-05-27得票数 0
回答已采纳
1回答
我想使用线性插值为每一天创建新的层。在这种情况下，前两个光栅属于Feb和29 days月份，第二个栅格属于March和31 days。我想知道如何创建每日光栅对象，通过考虑月份的天数(2月为29个栅格，3月为31个栅格)来填充时间段。谢谢!library(raster)feb15 <- setValues(r1, runif(ncell(r1))
浏览 3提问于2015-10-26得票数 1
回答已采纳
2回答
调用writeRaster overwrite=TRUE时，新的光栅值保持不变。我最初编写了错误的栅格对象，然后更正了代码，并以相同的文件名编写了一个新的栅格。写入文件的属性表中的值与原始文件中的值相同，即使我正在写入的栅格对象在R中查看时具有正确的属性。
解决方法是为新光栅指定一个不同的名称(或手动删除旧光栅)。R 3.0.0，Windows 7 64-b
浏览 0提问于2013-06-01得票数 2
回答已采纳
我有一个非常大的光栅文件，其中包含的单元格的值从1到25不等。我想要创建25个单独的栅格，其中只包括从原始单元格与相同的值。我已经查看了r的“光栅”包，并认为RasterFromCells可能是我的答案，但它需要单元格号，而且由于我原来的光栅有171,681,006个单元格，为每个单元格提供每个值的小区号将是非常耗时的。在r中是否有一个函数，我可以用来自动创建一个新的栅格，其中包括值为1的所有单元格
浏览 0提问于2019-01-03得票数 0
回答已采纳
1回答
我在R中有两个栅格对象，其中216个栅格或单元具有相同的经度和纬度坐标，一个栅格具有每个栅格的日最大风速值，另一个栅格具有表示单个栅格中的人口计数的数字。第一个栅格有183个层，由下面的信息表示。50.625, 53.625 (xmin, xmax, ymin, ymax)source : IPSL-IPSL-CM5A-MR_
浏览 23提问于2021-10-29得票数 2
回答已采纳
1回答
在R中按属性提取栅格
、、、
我一直在尝试创建一个新的栅格对象，它只包含现有栅格中的几个值。我使用的是这里找到的类栅格：。# but the file is to bigs <- r[r>30 & r<33] # but the file is to bigrc <- reclass(r, c(-Inf,30,NA, 31,32, 1, 33,Inf,NA))
这会生
浏览 0提问于2012-03-03得票数 0
回答已采纳
在R中，但是我很难让它对每个栅格应用一个简单的条件语句(所有栅格层都是相同的范围，分辨率，并且是共同注册的)。我希望将所有小于每个栅格最大值的95%的像素值设置为NA。下面是我的代码： #Get max value from each raster layerrstack[rstack < (r
浏览 12提问于2020-03-26得票数 0
回答已采纳
2回答
我在R中使用非常大的栅格，我有一个代码，我以前已经能够用相同的大栅格运行过，但是当我试图再次运行我的代码以获得可重复性时，我得到了相同的错误“在栅格IO期间失败”。例如，在代码的这一部分中，我试图将一个大的栅格(~ 5 GB)分成较小的块，这样我就可以在R中加快处理时间。getRasterData(object@file@con, offset = offs, region.dim
浏览 67提问于2020-05-29得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
