在 R 中，似乎 raster 包不区分 GeoTIFF 的正旋转和负旋转。我有一种感觉，那就是 R 忽略了旋转矩阵中的负号。我不太懂得深入挖掘 raster 源代码来验证，但我确实创建了一个可重现的示例来演示这个问题:读取 R 徽标并保存为 GeoTIFF。library(raster) b...

我想在R Studio 中制作一个类似于 Arc Map 中下面的图表： 我尝试了以下代码：# data processing library(ggplot2) # spatial library(raster) library(rasterVis) library(rgdal) # te...

我正在进行一项物理实验，使用C语言参与杨氏干涉实验，我编写了一个程序，将大量像素打印到文件中：for (i=0; i < width*width; i++) { fwrite(hue(raster_matrix[i]), 1, 3, file); } 当给定一个值[0..255]...

我怎样从时间序列中提取日期？ 这是一个时间序列：x = seq (1, 768) myts <- ts(x, start=1982, frequency=24) 最初我需要创建一个向量来保存rts函数的日期/时间数据，观测开始于1982年，每月进行2次测量，一直持续到2013年。

我在R中有几个高分辨率的栅格图层，正在处理中。对于一些分析而言，细节级别太高了，我希望通过降低分辨率来加快速度。 坐标系统是UTM，单位为米。分辨率是30,30（x,y）。因此，分辨率似乎是30米。 请问有人可以建议我如何将分辨率改为120米？我已经阅读了resample()和projec...

我一直在尝试找到一种在R中合并多个光栅图像的高效方法。这些是来自南部乞力马扎罗地区的相邻ASTER场景，我的目标是将它们组合在一起以获得一个大图像。 这就是我目前所拥有的（对象“ast14dmo”代表了一个RasterLayer对象的列表）：# Loop through single AST...

我正在使用 {raster} 对基于非规则形状文件（亚马逊生物群落）剪切（或裁剪）栅格，但输出总是具有矩形范围。但是，我需要输出与形状文件完全相同的几何图形。有什么提示吗？谢谢。library(raster) library(rgdal) myshp <- readOGR("Amazo...

我有一个数据框，其中指定了笛卡尔坐标（x，y）的值（l），如下面的最小工作示例。set.seed(2013) df <- data.frame( x = rep( 0:1, each=2 ), y = rep( 0:1, 2), ...

我想创建一些GIS图表，想知道R语言能否用于这个任务。以下是一些类似于我想制作的图表示例： 一个美国的温度图（或等高线图），根据各州GDP的大小决定颜色（或高度）。因此，州界将在生成的图中产生不连续性。 一个美国的温度图，其中使用海拔高度进行数据分析。在这种情况下，结果图应该在各州边界之间...

我有一个由三个栅格图像组成的栅格堆叠stk，这是一个简单的示例：# set up a raster stack with three layers > library(raster) > r <- raster(nrows=10,ncols=10) > r[] <...