在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语句从栅格创建新栅格了。

请注意，以上答案中没有提及任何特定的云计算品牌商，因为问题与云计算无关。如果你需要了解与云计算相关的内容，请提供相关问题。

相关·内容

C++ GDAL创建大量栅格文件并分别写入数据

，如果大家需要生成别的格式的话可以修改这里；auto x : my_file表示从我们前期已经获取到的需要处理的栅格文件列表中遍历（虽然我们这里是需要建立新的栅格文件，但由于我这里新的栅格文件的命名规则是与原有的栅格文件一致的...，所以就还是从原有的文件列表中遍历），my_file就是前期已经获取到的需要处理的栅格文件列表，具体获取方法可以参考文章C++遍历文件夹从而获取指定格式或名称的文件。 ...在GDAL库中，如果我们想用自己的数据生成栅格文件，首先需要基于CreateCopy()函数新建一个栅格文件，随后通过RasterIO()函数写入数据。...完成以上全部配置后，即可依据CreateCopy()函数进行新的栅格文件的创建。 ...至此，我们仅仅是完成了GDAL库中栅格文件的创建，但此时还没有将数据导入进去，因此在资源管理器中也是看不到具体的新的栅格文件的。

3002 0

PostGIS空间数据库简明教程

图片要使用“geometry”数据类型创建一个表，我们可以运行下面的语句：CREATE TABLE building (id UUID PRIMARY KEY,geom geometry);这将创建一个表...图片栅格数据支持包含在一个单独的 postgis 扩展中，在我们可以使用之前需要安装它：CREATE EXTENSION IF NOT EXISTS postgis_raster;然后我们可以使用栅格类型创建一个表...为了解决这个问题，空间索引使用 R-Tree（“Rectangle”中的“R”）结构，它构建了一个矩形树，其中每个子节点矩形都包含在父节点矩形中。...相同的语法可以应用于栅格列，但在这种情况下，我们在栅格图像周围索引边界框，因此该语句需要包含 ST_ConvexHull 函数。...，因为 PostGIS 需要将新对象插入 R-Tree 索引。

2.9K3 0

executeupdate mysql_使用Mysql中的executeUpdate在SQL语句中创建表

catch (SQLException e) { // TODO Auto-generated catch block e.printStackTrace(); } } executeUpdate创建...DB并使用他的前两个方法工作正常。

5.5K2 0

arcpy怎么用_python arcpy

r”c:\data” “\”在编码中需转义，前面加r,直接按照原文意思输出。...arcpy在解决以上的差异性问题时，提供的工具有：使用游标访问矢量数据单个属性项数据，使用numpy数组遍历栅格数据单个栅格值。矢量属性访问游标是包含从要素类或表中获取的一行或多行数据的内存对象。...向表中插入 25 个新行。...String where_clause 用于限制在游标中返回的行的可选表达式。有关 where 子句和 SQL 语句的详细信息，请参阅在查询表达式中使用的元素的 SQL 参考。...然后，将会创建一个新的栅格数据。

2.3K2 0

遥感影像依据分幅号的8天合成：Python实现

如下图所示，可以看到在GF1WFV3传感器对应的2022年346天遥感影像结束后，新的GF1WFV4传感器对应的遥感影像又是从2022年的开头开始的。...在之前的文章中，我们介绍过在GEE中计算每1个8天时间间隔内遥感影像数据平均值的方法；而这一次我们将基于Python，将每1个8天时间间隔内遥感影像拼接起来。本文所用到的代码如下。...接下来，基于ListRasters("*", "tiff")获取待拼接的所有.tif格式栅格文件，并将其排序后存储在image_list列表中；image_dict是一个字典，用于存储栅格影像按日期和分幅号进行分组的结果...）；随后，使用CreateRasterDataset_management()函数创建输出栅格数据集，命名规则为年份+天数+分幅号。...同时，打印信息，指示栅格数据集创建完成。最后，即可使用Mosaic_management()将影像文件列表拼接为一个栅格数据集，命名规则同上；同时，打印信息，指示栅格数据集拼接完成。

1151 0

C++中GDAL：用像素均为0的栅格填补时序遥感数据中缺少的时相

首先，我们需要基于文件夹中遥感影像文件的文件名称特征，遍历生成文件名列表。在这里，我们使用两个嵌套的for循环，生成所有可能的栅格图像文件名，并将这些文件名保存在all_file_path向量中。...接下来，我们使用GDALOpen函数，从2018001这一天开始，通过循环打开对应名字的文件，直到找到文件夹中第一个实际存在的栅格图像文件（poDataset_actual），并获取其栅格图像的行列数（...对于不存在的栅格图像文件，使用GDALDriver创建一个新的数据集（poDataset），并将其中的像元值设置为0。如果栅格图像文件已经存在，则跳过不处理。...其中，在对缺失的栅格图像加以生成时，我们首先使用GetGDALDriverManager()->GetDriverByName函数获取GDAL驱动程序对象，然后使用CreateCopy函数创建新的栅格图像...随后，我们用0填充新创建的栅格图像，并使用RasterIO函数对栅格图像的像元进行写入操作。最后，在上述处理完成后，使用GDALClose函数关闭数据集，并输出新创建的栅格图像的文件名。

2113 0

【GEE】4、 Google 地球引擎中的数据导入和导出

3将您自己的数据带入 Earth Engine 在本练习中，我们将讨论如何将您自己的数据移动到 GEE、从数据集中提取值以及从 GEE 中导出这些值。...3.3定义天气变量在本课中，我们使用 Google Earth Engine 作为将遥感数据（即我们的栅格）与我们的点位置相关联的方法。...我们可以使用这些信息来确定哪些数据与美洲狮在特定日期的位置有关。重要提示：在图像集合中有许多图像，我们将创建具有大量波段的单个图像。...虽然我们可以在 GEE 中更多地使用这些数据，但很容易将它们导入 R 或 Excel。有几个选项可以定义导出数据的最终位置。一般来说，将数据保存到 Google Drive 帐户是一个安全的选择。...从 GEE 导出特征时出现的用户定义参数示例。 3.5.2导出栅格在处理所有这些空间数据时，您可能已经意识到，在美洲狮上收集数据的时间段内显示中值的栅格可能是非常有用的信息。

9422 1

交通时空大数据如何分析，我写了本书！

在我们的日常活动中，手机会记录下我们到访过的地点；在使用城市公交IC卡、共享单车等服务时，服务供应商会知道这些出行需求产生的时间与地点；公交车与出租车的定位信息，也可以告诉我们城市交通状态的具体情况。...特定处理方法：针对各类特定数据提供相应处理方法，如从出租车GPS数据中提取订单起讫点，从手机信令数据中识别居住地与工作地，从地铁网络GIS数据构建网络拓扑结构并计算最短路径等。...TransBigData可以通过pip或者conda安装，在命令提示符中运行下面代码即可安装： pip install -U transbigdata 安装完成后，在Python中运行如下代码即可导入...在使用预处理方法时需要传入数据表中重要信息列所对应的列名，代码如下： # 数据预处理#剔除研究范围外的数据，计算原理是在方法中先栅格化后栅格匹配研究范围后实现对应。...中，我们可以对出租车数据使用简单的代码在jupyter notebook中快速进行交互可视化。

2.1K1 0

原创丨利用ArcGIS制作多年栅格渐变

● 托管方式：如果由地理数据库托管栅格目录，那么栅格目录中的栅格数据集将以物理方式存储在地理数据库中。如果某行或某栅格从目录中删除，则该行或栅格也将从地理数据库中删除。...创建好之后的成果如下图： STEP 02 找到原始的栅格数据，批量建立金字塔打开“工具箱”，找到“数据管理工具”，点击“栅格”，打开“栅格属性”；在弹出的界面中，点击文件夹图标，添加我们所有的栅格数据到下图的红色框线内...STEP 04 加载输出完毕的栅格目录进ArcMAP，并插入新的数据框，设置数据驱动工具看下图，刚刚添加进来的数据是这样的，是一个矩形面，并没有我们想用的栅格数据显示在工作界面中，出现这种问题的原因呢...我们接下来呢，插入一个新的数据框，并把我们的栅格目录复制粘贴到新的数框中。图中显示的是新插入的数据框内新粘贴的栅格目录。...，下拉选择从栅格目录，在弹出来的计算统计数据面板中，自己选择一种计算方式，然后计算添加；把RGB合成与色彩映射表两个渲方式移除，最后点击确定，查看成果。

1.9K2 0

GeoTrellis整体介绍

，元数据转换成JSON GeoTrellis解决了三个核心的问题 1.创建可扩展的，高性能的地理信息处理WEB服务 2.创建分布式的地理信息处理服务，用来处理海量数据集 3.完成并行化地理信息处理操作...，以利用多核架构的优势 GeoTrellis可以将数据(Tiff) 从本地，HDFS，S3中导入到本地，HDFS，Accumulo，HBASE，CASSANDRA，S3等，可选方式很多，而且是通过Spark...GeoTrellis的第一步就是要将数据切片(无论是存储在内存还是持久化)，然而即使能力再大，在实际工作中也难以处理以下几种需求：全球(大范围) 高分辨率遥感影像数据，数据量在TB级局部地区数据更新...多波段瓦片读取 val multiTile = tileReader.reader[SpatialKey, MultibandTile](LayerId(name, zoom)).read(key) 从多波段中获取单个波段...Http.Bind(service, host, port) 需要使用以下语句系统遍自动的在host和相应的port上发起服务。具体路由信息需要在service类中定义。

2781 0

ArcPy切分大量遥感影像为多个网格区域的方法

现有同一区域内的多张栅格遥感影像，我们希望对于每一景栅格遥感影像而言，都将其分割为多个小矩形；其中，分割后的每一个小矩形就是一个新的小的栅格文件。知道了需求，我们便可以开始具体的实践操作。...其中，有2种实现这一需求的方法——首先第一种方法，我们可以创建一个渔网的矢量文件，并依据这一渔网矢量文件对每一景栅格加以分割，如下图所示；其中，关于渔网矢量文件的创建，大家可以参考文章ArcGIS创建渔网并批量获得指定大小的网格矢量...，关于接下来用渔网中的格网去裁剪栅格图像，大家可以参考文章依据矢量要素批量分割栅格为多部分：ArcPy。 ...其中，我们前面也提到了，本文代码与上述文章代码不一致的地方在于，这里是在SplitRaster_management()函数中手动划定每一个小矩形的具体大小，而上述文章则是用一个已经得到的矢量要素文件去分割栅格...此外，其他的参数就和前述文章比较一致了（当然也要记得将前述文章中关于矢量要素的几个参数设置为空白）。在 IDLE (Python GUI) 中运行代码。

1981 0

我用MATLAB撸了一个2D LiDAR SLAM

*思路为：在预测的下一位姿上做一些细小的调整(对x、y、theta做细小调整)；对于某一次调整后的预测下一位姿，利用下一位姿的扫描数据，构建下一位姿的栅格地图；以下一位姿的栅格地图与当前位姿的栅格地图的重合度作为目标函数...如果我们在最后一步执行了 mini更新，我们将更新局部点集图和局部栅格地图（粗略） % 1....localMap 栅格单元尺寸对应实际长度以pixelSize 创建占用栅格地图 gridMap2 = OccuGrid(localMap, pixelSize/2);%从点集localMap...栅格单元尺寸对应实际长度以pixelSize/2 创建占用栅格地图 end % 2....& points(:,2) < maxY; %从全局地图中提取到的当前扫描点 localMap = points(isAround, :); (7)OccuGrid.m % 从点集创建占用栅格地图

1.8K4 0

ArcGIS空间分析笔记（汤国安）

关系类的创建是在源类的主键和目标类的外键之间创建的。主键：是储存能够唯一标识表中的每一个对象的字段。外键：记录有源表主键信息的字段。...使用向前标注可以从源类找到目标类使用向后标注，可以从目标类找到源类注释类注释是用于储存描述性文本信息的专门要素类，和储存在地图文档中的标注不同，注释类储存在地理数据库中。...创建新的几何网络需要创建在要素数据集下最好在内容列表中把数据导入定义投影指按照地图信息源原有的投影方式，为数据添加投影信息。...数据提取数据提取是从已有的数据中，根据属性表内容选择符合条件的数据，构成新的数据层可以通过设置SQL表达式进行条件选择泰森多边形用途——定性分析、统计分析、邻近分析网络中流动的管线...; 使用TrackingAnalyst回放管理器回放时间数据，可使用不同的速度进行正向和反向数据回放; 通过创建数据时钟来分析时间数据中存在的模式; 针对时间数据创建和应用操作;

3.3K2 0

遥感数据机器学习的准备工作：python将栅格数据提取至EXCEL

大部分我们处理的降水、气温等栅格数据的格式是nc形式，需要我们将他转换成栅格数据并导入至Arcgis中，进行下一步操作。...，下面我们通过python代码实现快速的栅格属性导出至EXCEL表中。...:\CN\RasterT_tif1.shp' #点数据文件 InputRasterFolder = r'E:\CN\train' #放栅格数据的文件夹 # 设置Excel编码 file = xlwt.Workbook...#excel表的第2列为纬度 for i in range(0, len(tif_files)): #遍历tif sheet1.write(0, i + 2, filename) #在表格第一行设置列名...#读取影像 data = band.ReadAsArray(xOffset, yOffset, 1, 1) #读出从(

1.3K2 0

Python识别文件名中的字段从而分类、归档栅格文件到不同文件夹

本文介绍基于Python语言，针对一个文件夹下的大量栅格遥感影像文件，基于其各自的文件名，分别创建指定名称的新文件夹，并将对应的栅格遥感影像文件复制到不同的新文件夹下的方法。 ...其中，如上图中紫色框所示，每一景遥感影像文件的文件名称中，都有一个表示其编号的字段；我们希望基于这一编号字段，将带有相同编号字段的栅格遥感影像文件，以及其对应的辅助信息文件，都复制到一个结果文件夹中；这个结果文件夹如下图所示...例如，我们希望将所有文件名称中带有15字段的栅格遥感影像文件及其辅助信息文件，都复制到结果文件夹中名称为15的子文件夹中，以此类推。知道了具体需求，我们即可开始代码的撰写。...但是这里需要注意，由于我们需要将每一个文件都放入结果文件夹中的子文件夹，因此需要首先判断当前子文件夹有没有被建立；如果没有创建的话，我们需要创建一下这个子文件夹。...这一步骤通过if判断语句及其内部的代码即可实现——通过os.path.exists()函数判断是否存在指定的子文件夹，如果不存在的话就新建这一子文件夹。随后，就可以开始文件的复制工作了。

1471 0

R语言计算大量栅格图像平均值、标准差

本文介绍基于R语言中的raster包，批量读取多张栅格图像，对多个栅格图像计算平均值、标准差，并将所得新的栅格结果图像保存的方法。 ...在文章R语言raster包读取栅格遥感影像中，我们介绍了基于R语言raster包，对单张或多张栅格图像加以平均值、标准差计算的方法；但这一篇文章中的标准差计算方法仅仅可以对一张栅格图像的全部像元加以计算...本文就介绍另一种方法，可以对多个时相的大量栅格影像加以逐像元平均值、标准差的计算，从而使得最终的结果是一景表示各个像元在全部时相的图像中的平均值或标准差的图像。 ...当然，前述提到的文章R语言raster包读取栅格遥感影像中的方法也是可以对多个栅格图像计算平均值的。...tif_sd <- calc(tif_file_all, fun = sd) plot(tif_sd) 此外，上述代码在calc()函数运行时，若某一空间位置上的像元在多张栅格遥感影像中，存在至少一个无效值

5762 0

大栅格数据如何更快运算

背景介绍这两周我在使用python进行大量的栅格数据的运算，在运算过程中遇到了数据量超级大但算力不足的问题。通过这两周的探索，也慢慢找到了一些加快栅格数据计算的方法，和读者分享。...问题与解决方法（1）数据量超过电脑内存，使用分块运算在计算栅格数据时，是把数据放到内存中进行计算，如果栅格计算数量巨大，会爆内存。分块方法就是采用横纵切割原始栅格，最后再将数据拼接起来。...栅格的运算一般使用的是numpy模块，然后将数据转为数组array放到内存中计算。但如果你的栅格数据过大，就需要用到mmap_array，这是一个内存映射数组，可以保存到硬盘中。...这种数组的数据存储在磁盘上的一个文件中，而不是直接存储在内存中。numpy.memmap 的主要优点是，它允许您处理比可用内存更大的数据集，因为数据只在需要时才从磁盘加载到内存中。..._mmap.close() （3）数据读取写入速度慢，在固态硬盘中运行刚开始我使用了机械硬盘作为数据的运行盘和数据保存盘。

3322 0

R语言批量拼接遥感影像

本文介绍基于R语言中的raster包，遍历文件夹，读取文件夹下的大量栅格遥感影像，并逐一对每一景栅格图像加以拼接、融合，使得全部栅格遥感影像拼接为完整的一景图像的方法。 ...其中，本文是用R语言来进行操作的；如果希望基于Python语言实现类似的批量拼接、镶嵌操作，大家可以参考Python中arcpy栅格创建与多景遥感影像批量拼接与ArcPy依据成像时间分别批量拼接不同时相的遥感影像这两篇文章...raster包；关于这一包的配置，大家可以参考R语言raster包读取栅格遥感影像。...max与TRUE，二者分别是指重叠区域的像元以2景遥感影像中的最大值像元为准，以及在计算时不考虑NoData值的影响。 ...其实，这里的merge()函数与前述的mosaic()函数功能大致一样，但merge()函数在处理重叠区域时，默认选择位于顶层的遥感影像的像元数值，就没有mosaic()函数中的这么多计算方法选择了。

2591 0

ArcGIS PRO基础教程（一）

前言大家好，我是南南好久不见，这次准备写几个PRO的入门教程，让大家了解一下PRO中的GIS操做特点和新的功能另外由于公众号改版，大家可以加个给我星标，不然可能会收不到消息操作要求 1.面积为...实验数据操作步骤创建工程，模板选地图就可以了（注：在arcgis pro中创建工程可以看作在arcgis中创建mxd文档，不同的是arcgispro可以选择创建二维，三维等等，同时会给创建的工程额外创建一个工作目录...，并添加在线地图作为底图）添加数据还是和arcgis一样，链接文件夹，然后再添加数据根据操作要求【不能选在有耕地、园地内】选择出合适的土地利用地区，然后导出一个新的shp图层（注：pro中的选择语句变成中文的了...，方便我们使用，不像以前sql，还要注意符号啥的）根据操作要求【坡度小于15度，高程在以下1930】需要使用DEM进行处理分析，根据等高线图CONTOUR创建TIN(注：坐标系与CONTOUR相同...打开栅格转面工具根据操作要求【距离水源地在300米以内】对水源地图层water建立缓冲区对水系缓冲区，坡度，高程，土地利用图层相交根据操作要求【面积为50-80亩】（注：1亩=666.67

1.5K3 0

求栅格序列每个像元的变化趋势和对应P值

在处理的过程中，我们是求每个像元在时间维度上的变化趋势，类似下图：引用自arcgis网站也就是说我们对上图中的每一个条柱时间序列求趋势即可。有了思路，就非常简单了，我们直接上代码。...数据创建这里说的数据创建就是把我们的多个栅格序列组成上面类似的时空立方体（这里仅指时空栅格）。.../slope.tif') as src: show(src) 到这里就完成了每个像元的线性趋势计算，不过上面的代码只保存了趋势值，并没有保存R方和p值，读者根据代码改一下即可。...因为在计算趋势的时候，如果你不关心截距，那么年份是从0-35还是1980-2015，你算出来的趋势值（也就是下面公式中的a，x是年份）都是一样的，那么就不必要多浪费那点算力了： Y = ax + b...总结处理栅格序列的时候，元信息一般不变，所以可以利用某一个原始数据的元信息作为模版，方便保存处理后的结果；对于栅格数据的值，就是一个数组而已，巧用numpy的函数可以实现很多我们需要的功能；在能简化算法的时候

2.7K4 0

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云