raster

0热度

1回答

我想我可能发现了一个错误。我在栅格窗口中绘制栅格窗口时显示正常的栅格，但是当我尝试使用png()来保存栅格时，我在rasterVis中使用了levelplot。 library(raster) library(rasterVis) library(RColorBrewer) col<-rasterTheme(region=brewer.pal(8,"YlOrRd")) r<-raster(

1热度

1回答

[R rasterFromXYZ：X单元尺寸都是不正规的

给定一个rasterLayera， library(raster) library(rasterVis) a=raster('p190001.grd') head(a) 基于下面的说明（见下文），这也是在这个环节中， ftp://ccrp.tor.ec.gc.ca/pub/EC_data/CANGRD/ 我想出了具有以下CRS mycrs <- CRS("+pro

2热度

1回答

parLapply多个参数ř

我试图通过hargreaves方法来计算蒸发内package SPEI。这包括使用最小温度（TMIN）和上限温度（TMAX）。并行计算是我最好的选择，因为Tmin和Tmaxrasterstacks有500,000 cells and 100 layers each。 Hargreaves function以Tmin,Tmax和latitude在each grid作为输入。以下是我的第一个猜测如何

1热度

1回答

计算区域中的R

在光栅不同的土地覆盖类我怎样才能给出下面的例子中我的土地覆盖类别1,2,3,4,5的总面积： library("raster") r <- raster(nrow=10, ncol=10) r<-setValues(r,c(rep(1,20),rep(2,20),rep(3,20),rep(4,20),rep(5,20))) 一方法可能是基于这些值对栅格进行子集化处理，但是像a <- ar

0热度

1回答

在R光栅循环中找不到文件

这里我试图将CT图像的tiff文件导入R进行分析。我有大约250个tiff文件，我正试图在一个循环中进行分析。分析首先进行裁剪，然后进行阈值处理，然后根据像素强度值计算孔隙度。目前，我得到这个错误： Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", : Cannot create a RasterL

1热度

1回答

从栅格中提取沿着定制的圆形条（恒定的径向距离）的值

我有一个从.tif图像文件导入的大型栅格对象。我想在距离中心的给定径向距离处分析这个光栅的像素，以便识别我可以在图像中注意到的某种轴对称现象（下图）。为此，我想从图像中心提取与给定半径（和宽度）的圆形条相交的像素值（如下图所示）。我已经探索了几个选项来做到这一点，包括提取功能和成像器包。在imager软件包中，您可以很容易地沿着行或列提取值，但是我找不到可以按照我需要的自定义形状提取值的函数。

0热度

1回答

R如何提取海拔高于600m的每个多边形栅格的平均值？

首先，我设法提取平均光栅温度值为每个多边形，而下面的程序：您可以下载这个链接上的GIS图层： https://depots.univ-perp.fr/get?k=iTzEDSUkdyZVw2st78G ## load packages library(raster); library(rgdal) ## Read rasters ras_temp<-raster("ras_temp.

0热度

1回答

NetCDF - 转换成光栅和投影问题

我有以下NetCDF文件 - 我想转换成光栅，但有些东西是不正确的。没有给出NetCDF文件的预测，但基于我收到它的软件应该是LatLong，但可能是圆柱面积相等的区域。我尝试了这两种方法，但是我一直得到这种失真，这使得无法在正确的位置查询值。我知道网格的间距不是偶数，不确定是否会影响最终结果（这里是ArcGIS的可视化，但在R中，除非用levelplot函数绘制，否则它是相同的问题）。 libr

0热度

1回答

在R中合并2个或更多光栅堆栈R

我有2个（或更多）光栅堆栈。每个人都进行了各种处理（任何类型的东西），并且这两个堆栈保持了他们的范围和解决方案（最初，两者都具有相同的范围和分辨率）。现在我想让它们成为一个堆栈，而不是2.所以，我会有一个'merged_stack'来包含来自两个（或更多）堆栈的所有栅格。例如： #creating some data m1 = matrix(1,30,30) m2 = matrix(2,30

1热度

1回答

在％>％管道内更改栅格值

我有一个栅格，其中我感兴趣的数据（物种存在）值为1，其余为0（不存在）或NA。我想使用magrittr和dplyr，并且在管道内，我想将所有0值更改为NA，以便稍后修剪栅格并获取物种所在数据的最小和最大坐标（值= 1）然而，改变值0 NA是我通常会喜欢的东西做的事： raster[values(raster) == 0] <- NA 我有麻烦做那样的操作，“从右至左”，用magrittr（