将多层 RasterBrick 保存到硬盘

2023-12-31

我有一个代表地形图的多层 RasterBrick,我想将其以 grd 或 tif 格式保存到硬盘,以便其他人以后可以使用。

这是 RasterBrick:

class       : RasterBrick 
dimensions  : 2400, 4200, 10080000, 3  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 480000, 9e+05, 62000, 302000  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : in memory
names       :  layer.1,  layer.2,  layer.3 
min values  :   2.8725,   2.8725,   2.8725 
max values  : 254.5175, 254.5175, 254.5175 

我尝试使用以下命令保存它:

outfile <- writeRaster(brick, filename='grid.tif', format="GTiff", overwrite=TRUE)  

和这个:

outfile <- writeRaster(m, filename='grid.grd', format="raster", overwrite=TRUE)  

但是 tif 文件已损坏,并且 grd 对象仅包含一层,当我使用它读回时,它不会被识别为多层 RasterBrickraster().

目的是利用地形图作为专题图的背景。


尝试这个:

outfile <- writeRaster(brick, filename='grid.tif', format="GTiff", overwrite=TRUE,options=c("INTERLEAVE=BAND","COMPRESS=LZW"))
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

将多层 RasterBrick 保存到硬盘 的相关文章

  • 如何减小 R Plot 中图例的大小,同时仍使其可读?

    我试图用 R 中的两个 y 轴绘制多年来的一些数据 但是 每当我尝试包含图例时 图例就会主导我的绘图 当我使用其他地方建议的解决方案时 例如keyword和 或使用cex论据 在另一篇文章中建议here https stackoverflo
  • 将 data.frame 的列中的值替换为另一个 data.frame 中的值

    我的情况是 我有一个数据框 其中有一列填充了整数 1 到 6 我想用更具描述性的标签替换这些整数 这些标签在另一个充当 键 的数据框中提供 V1 V2 1 1 LABEL1 2 2 LABEL2 3 3 LABEL3 4 4 LABEL4
  • 计算级别内的值

    我在 R 中生成了一组级别cut 例如假设 0 到 1 之间的小数值 分为 0 1 个区间 gt frac lt cut c 0 1 breaks 10 gt levels frac 1 0 001 0 1 0 1 0 2 0 2 0 3
  • 单击并按住 R 中的按钮闪亮?

    我希望能够通过单击 R 闪亮按钮来更改参数的值 所以我需要按钮 一个用于增加值 一个用于减少值 我想在按住按钮的同时保持值以一定的速度减少 增加 通过释放按钮的点击 动作应该停止 到目前为止我还没有找到这个选项actionButtons在
  • 在 for 循环中绘制的多个 ggplot2 绘图的网格

    作为一个新的 ggplot2 用户 我对可能性的数量感到有点迷失 并且很难在网上找到我认为简单问题的简单答案 我想在同一张纸上显示 ggplot2 的多个图 但知道这些图来自 for 循环 以下示例无法编译 仅用于说明 for i in c
  • profvis() 何时以及为何显示“源不可用”?

    我经常分析 R 代码 并大量使用 profvis 对于某些函数 浏览器窗口的上半部分会显示源代码 有时则不会 我不知道什么时候会出现这种情况 对我来说这似乎是随机的 有谁知道 profvis 何时以及为什么无法在顶部窗口中显示代码 发生这种
  • 有效地将环境从内部功能转移到全局环境

    我有一个在其中创建环境的函数 我希望将该环境分配给全局环境 目前我通过将环境分配给来做到这一点globalenv 作为最后一步 如下 funfun lt function inc 1 dataEnv lt new env dataEnv d
  • 正则表达式挑选括号之间的一些文本[重复]

    这个问题在这里已经有答案了 可能的重复 提取 R 中所有括号内的信息 正则表达式 https stackoverflow com questions 8613237 extract info inside all parenthesis i
  • glm() 模型的交叉验证

    我正在尝试对我之前在 R 中构建的一些 glm 模型进行 10 倍交叉验证 我对cv glm 函数在boot包 尽管我已经阅读了很多帮助文件 当我提供以下公式时 library boot cv glm data glmfit K 10 这里
  • readRDS() 加载额外的包

    什么情况下会出现readRDS R 中的函数尝试加载包 命名空间 我很惊讶地在新的 R 会话中看到以下内容 gt loadedNamespaces 1 base datasets graphics grDevices methods sta
  • 在ggplot2中创建部分虚线

    我正在 R 中创建一个图 并且需要创建一条线 其中某些值是投影 投影用虚线表示 这是代码 df data frame date c rep 2008 2013 by 1 value c 303 407 538 696 881 1094 gg
  • 为什么我在 ddply 中看到“错误:length(rows) == 1 is not TRUE”?

    我有一个数据框 比如工资单 例如 payroll lt read table text AgencyName Rate PayBasis Status NumRate HousingAuthority 26 843 00 Annual Fu
  • R 中带有边缘箱线图的直方图

    如何使直方图中的 X 轴与边缘箱线图匹配 data lt rnorm 1000 nf lt layout mat matrix c 1 2 2 1 byrow TRUE height c 1 3 layout show nf par mar
  • 使用 R 中的 tidyverse 重新调整因子和重新排序因子

    我想使用这些功能重新调平 and 重新排序 在我的数据框中 我了解重新调整级别的工作原理 但我不明白为什么我在 data frame 中看不到级别的变化 例如 假设我有鸢尾花数据集 library tidyverse head iris g
  • 如何使用r中的dplyr在特定位置插入空白行

    我想在数据框中的特定位置插入空白行 我的数据框是这样的 dat lt data frame group c rep A 1 rep B 4 rep C 2 rep D 2 group 1 A 2 B 3 B 4 B 5 B 6 C 7 C
  • 如何在 R 中查找平衡面板数据(又名,如何查找面板中的哪些条目在给定窗口内完整)

    我有来自 Compustat 的大量数据 我向其中添加了一些手工收集的数据 认真地从一堆旧书中手工收集 但我不想手工收集整个面板 只想随机选择一个子集 为了找到更大的集合 我从中随机选择 我想从 Compustat 的平衡面板开始 我看到p
  • 在 R Shiny 中显示/隐藏整个框元素

    我目前正在尝试找到一种方法来隐藏 显示 R Shiny 中的整个 box 元素 以及里面的所有内容 我想创建一个可能的按钮 它允许用户展开特定框 然后使用相同 甚至不同 的按钮隐藏它 我不想使用条件面板 因为我的应用程序非常大并且会产生一些
  • 绘制带有颜色渐变的geom_segment线? (或者还有另一种方法来强调开始与结束吗?)

    我的数据框中有大量行 100 000 的两组纬度和经度变量 我正在尝试绘制一个连接这两组坐标的图 即 从纬度1 经度1 to 纬度2 经度2 使用 geom segment 使用非常低的 alpha 使线条透明 因为线条太多 我想强调这些线
  • 一次评论多个对象

    假设您有一个变量 a 到 j 的列表 for x in 1 10 assign letters x x 您将如何评论那些最近创建的对象 我尝试过类似的事情 for x in 1 10 comment get letters x lt pas
  • R Shiny - 使用 DataTable 移动列名称

    我有一个非常复杂的闪亮代码 其中有几个面板和这些面板内的几个表格 启动应用程序时 列名称与列值正确对齐 但是 一旦我更改应用程序表格下的页码 列名称就会移动到左侧 而值仍保留在中间 如何强制应用程序使列名称与列值对齐 一个可重现的例子 li

随机推荐

  • 如何使用 sqlalchemy 在 Docker 中访问 postgresql?

    我正在尝试使用 SQLAlchemy 与 dockerized PostgreSQL 服务器交互 就像是 engine create engine postgresql user user password localhost 5432 d
  • 如何在RTL中使用时钟门控?

    我正在对一些时钟进行门控latch以及我设计中的逻辑 我在综合和布局布线方面没有太多经验 在 RTL 中实现时钟门控的正确方法是什么 示例1 always comb begin gated clk clk latch update en e
  • 如何使用 connect w/ Redux 从 this.props 获得简单的调度?

    我有一个简单的 React 组件 可以连接 映射一个简单的数组 状态 为了避免引用商店的上下文 我想要一种直接从道具获取 调度 的方法 我见过其他人使用这种方法 但由于某种原因无法使用它 以下是我当前使用的每个 npm 依赖项的版本 rea
  • 在表单提交时,从 javascript 发送带有表单值的 mailto

    我有一个表单 当提交表单 输入类型 提交 时 我想使用预先填充的电子邮件消息打开客户端默认邮件浏览器 因此 当用户单击提交时 需要发生两件事 打开电子邮件并提交表格 另外 如何使用表单中输入的值来预填充电子邮件 我是 javascript
  • 找不到创建 iOS 应用程序开发证书的选项 [关闭]

    Closed 这个问题不符合堆栈溢出指南 help closed questions 目前不接受答案 请参阅此图片 iOS 应用程序开发选项已禁用 我如何启用它 您需要成为团队代理或管理员才能创建开发证书 Source https deve
  • C++ 错误:与“operator=”不匹配

    给数组赋值时遇到问题 我创建了一个名为Treasury 我创建了另一个类TradingBook我想包含一个全局数组Treasury可以从所有方法访问TradingBook 这是我的 TradingBook 和 Treasury 的头文件 c
  • 如何将数据传输到另一个页面

    海 我有一个包含 150 个控件的 ASP NET 页面 我想将这些控件的数据传输到另一个 ASP NET 页面 什么方法最适合这项任务 控制数量可能会增加 提前致谢 有很多方法 使用查询字符串 可能不适用于您的情况 仅适用于传输少量数据
  • 引用 IBM.Data.DB2 中的代码使得该程序集对我的解决方案的其余部分不可用

    我有一个包含三个程序集的 C 控制台应用程序 Main Common and Utilities 在一个文件中Main集会 Main cs 我有一行 using Utilities 在目录内Common程序集 我有DLLIBM Data D
  • 从 pmdarima 导入 auto_arima 时:错误:无法从 'scipy.misc' 导入名称 'factorial'

    我有 python 3 7 1 和 scipy 版本 1 3 0 我在调用 auto arima 时遇到错误 无法从 scipy misc 导入名称 factorial 只是这个基本的导入导致了问题 从 pmdarima arima 导入
  • Angular2 CLI:为什么“--prod”的包大小小于“--prod --aot”?

    我正在为一个项目使用最新的 angular cli beta 18 我知道 cli 仍处于非常早期的阶段 但我很困惑为什么我的最终包大小实际上更小without AoT 当我跑步时ng build prod 为 1 08 MB Asset
  • 将命令发送到现有 MATLAB IDE 命令窗口

    打开 MATLAB IDE 后 我们是否仍然可以从另一个进程写入其命令窗口 在 Windows 中 我们可以使用 MATLAB Engine API 附加到现有的 MATLAB 但这在 Linux Mac 中不起作用 Linux 下我们有什
  • 在 AngularJS 的 ionic 框架中使用 ui-router

    我正在开发一个使用离子框架的应用程序 这又使用了 ui router 目前 我有一个非常基本的两页应用程序 然而 它会扩大到更大 此时 当我从第一个视图转换到第二个视图时出现错误 错误说 TypeError Cannot read prop
  • 如何解决 Media type=multipart/form-data 错误未找到 MessageBodyWriter

    如何为简单的多表单帖子配置提供程序 任何建议 指示将不胜感激 堆栈跟踪 org glassfish jersey message internal MessageBodyProviderNotFoundException MessageBo
  • 从 ControlTemplate 中的样式 DataTrigger 进行模板绑定

    在下面的 XAML 中 我使用带边框的矩形作为 ToggleButton 的模板 我希望 BorderBrush 具有不同的颜色以反映 ToggleButton IsChecked 的变化值 不幸的是 我在这里尝试在 DataTrigger
  • pcap 纳秒 Python

    有没有办法使用现有的 python 库从 pcap 中获取纳秒 我有一个纳秒 pcap 文件 可以与 Wireshark 配合使用 但与 Python 配合使用效果很好库甚至不会导入该文件 此功能确实存在于 c libpcap 中 请参阅
  • YouTube v3 API 对于有效请求返回状态代码 400

    目前 我正在尝试使用 OAuth 对 YouTube API 执行看似简单的请求 我可以使用评论端点成功请求评论列表 但是 当尝试删除评论或设置其审核状态时 就会出现此问题 在审核状态更新的情况下 我将向以下地址执行 POST 当然使用有效
  • Swing 渲染器在列中显示完整文本会创建无限循环

    我想渲染一个摆动表 根据列中显示的文本量来调整行的高度 以便所有文本都显示在单元格中 问题是 当我设置行的高度时 会创建无限循环 我该如何解决呢 这是我的渲染器 final class RenderTextAreaMessage exten
  • 会话超时时调用 Ajax

    所有 如果在会话超时时进行 ajax 调用 我会尝试重定向到登录页面 这是我到目前为止所做的 为所有操作定义一个操作过滤器 public class AuthenticateFilterAttribute ActionFilterAttri
  • 为什么不建议@import链接css?

    Why import不建议链接css 什么 import不能做哪些可以做 import 与浏览器有任何不兼容吗 What in 使其可取 还有有没有什么专业 import哪个有用 我们应该总是 总是使用 这是一篇文章的链接 该文章讨论了性能
  • 将多层 RasterBrick 保存到硬盘

    我有一个代表地形图的多层 RasterBrick 我想将其以 grd 或 tif 格式保存到硬盘 以便其他人以后可以使用 这是 RasterBrick class RasterBrick dimensions 2400 4200 10080