raster:修改resample()函数中的方法
raster: 修改resample()函数中的方法
在云计算领域中，raster（光栅）是一种用于表示和处理图像、地理空间数据等的数据结构。它由像素组成的网格，每个像素都包含了图像或地理空间数据的值。在处理和分析大规模的光栅数据时，常常需要对其进行重采样（resample）操作，以改变其分辨率或投影等属性。
resample()函数是用于执行重采样操作的函数，它可以根据指定的方法对光栅数据进行修改。在修改resample()函数中的方法时，可以考虑以下几个方面：
- 重采样方法：根据具体的需求，可以选择不同的重采样方法，例如最近邻插值、双线性插值、三次样条插值等。不同的方法会对重采样结果产生不同的影响，需要根据具体情况进行选择。
- 分辨率调整：重采样操作可以改变光栅数据的分辨率，即像素的大小。可以根据需要调整分辨率，例如增加分辨率以提高图像质量，或减小分辨率以降低数据大小。
- 投影变换：在进行地理空间数据的重采样时，可能需要进行投影变换以适应不同的坐标系统或地图投影。可以在resample()函数中添加投影变换的功能，以便在重采样过程中进行坐标转换。
- 边缘处理：在重采样过程中，可能会出现边缘效应，即重采样后的图像边缘出现锯齿或模糊等问题。可以通过在resample()函数中添加边缘处理算法来解决这些问题，例如使用卷积核进行平滑处理。
- 并行计算：对于大规模的光栅数据，重采样操作可能会非常耗时。可以考虑使用并行计算的方法来加速重采样过程，例如使用多线程或分布式计算框架。
腾讯云提供了一系列与光栅数据处理相关的产品和服务，例如腾讯云地理信息服务（Tencent Cloud GIS），可以用于处理和分析地理空间数据。具体的产品介绍和文档可以在腾讯云官网上找到。
请注意，以上答案仅供参考，具体的修改方法和推荐产品需要根据实际需求和情况进行选择。
相关·内容
1回答
我想将两个栅格相乘(在R中)，但它们并不完全匹配。可以对栅格进行倍增的选项有哪些？ > veg_g
class : RasterLayer
dimensions : 2416, 2807, 6781712 (nrow, ncol, ncell)
resolution : 0.25, 0.25 (x, y)
extent : 109067.2, 109769, 533016.5, 533620.5 (xmin, xmax, ymin, ymax)
crs : +proj=sterea +lat_0=52.1561605555556 +lon_0=5
浏览 38提问于2021-07-31得票数 0
回答已采纳
我有一个50米分辨率的光栅。我正试着跑到30米。我试过使用raster::disaggregate函数，但是分辨率没有改变到30米。 library(raster)
a = raster("file.tif")
b = disaggregate(a, fact=30/50, method='bilinear') 如果你有其他的解决方案，这将是非常有帮助的。提前谢谢。
浏览 36提问于2021-10-11得票数 0
回答已采纳
我有三个光栅。Raster1是一种土地覆盖类型的土地覆盖文件。Raster2和raster3是显示变量'NPP‘的光栅。正如你所看到的，每个光栅都有不同的范围和分辨率。我想知道根据raster1的地形，在栅格2和栅格3中的NPP是多少。但是，如何才能使所有栅格达到相同的范围和分辨率，并根据raster1中的landcover类在raster2和raster3中找到NPP？ (我如何知道应该为所有栅格选择哪种分辨率？) > raster1
class : RasterLayer
dimensions : 2803, 5303, 14864309 (n
浏览 21提问于2020-03-16得票数 0
回答已采纳
1回答
我有两个同一天的栅格，但有不同的线条。我想把它们结合起来，但我意识到检索算法可能是不同的。这两个光栅的尺寸相同。请问用R做这件事最简单的方法是什么？我会在一个列表上运行它。
library(raster)
A <- raster(nrows=108, ncols=21, xmn=0, xmx=10)
A[] <- 1:ncell(A)
xy <- matrix(rnorm(ncell(A)),108,21)
B<- raster(xy)
## Induce NAs in raster B:
B[sample(1:ncell(B), 1000)] <- NA
浏览 0提问于2016-02-10得票数 1
1回答
我有分辨率为0.5度的光栅数据。我想用最近邻插值将它们插值到0.125度和1度。
如何在R中实现这一点？以下是一些示例数据：
library(raster)
ras <- raster(res=0.5)
#get polygons for Iceland
iceland <- getData('GADM', country="Iceland", level=1)
cr <- crop(ras, iceland, snap="out")
values(cr) <- 1:ncell(cr)
fr <- mask(cr,
浏览 3提问于2020-07-20得票数 1
回答已采纳
1回答
用for循环重采样多个光栅
、、、、
我试图在列表中的一些栅格上循环一个重采样函数，以便它们的尺寸、范围和分辨率与列表中的第一个栅格相匹配。但我一直收到以下错误
(函数(类、fdef、mtable)中的错误:无法找到用于签名‘list’、“RasterLayer”的函数“重采样”的继承方法
下面是我的代码：
library(raster)
rastlist <- list.files(path = ".", pattern='.tif$',
all.files=T, full.names=F)
allrasters <- lapply(rastlist, raster)
nraster
浏览 1提问于2021-04-09得票数 1
2回答
光栅重采样非常慢
、、、
)我想重采样哨兵-2图像，以便两个波段的分辨率适合。我使用的是raster：：raster(raster1，raster2，method="ngb")，但是执行过程非常长。它至少运行了15分钟，到目前为止我还没有得到任何输出。我不确定这件事是否会在某个时候结束。
我用的是哨兵2号的08波段(10米分辨率)和12波段(20米分辨率)。
用QGIS重采样同一光栅需要一分钟左右。我想问一下，这是否是此函数的常见问题，以及是否有任何解决办法或其他函数/包可用于此。
谢谢:-)
浏览 5提问于2021-08-18得票数 0
回答已采纳
1回答
我有两个栅格，我想使一个空间范围到另一个。那就把它保存成一个新的光栅。我用了下面的代码。然而，我无法将2013年的图像保存为新的空间范围作为新的栅格。如有任何指导，将不胜感激。
raster_2013 <- raster("avgt2013.tif")
extent(raster_2013)
class : Extent
xmin : 112.91
xmax : 153.64
ymin : -43.75
ymax : -9
> res(raster_2013)
[1] 0.01 0.01
浏览 2提问于2017-07-31得票数 2
2回答
用R语言将哨兵-2波段降到10米
、、、、
我试图用R语言中的使用来计算。
NDRE = (nir-re)/(nir+re)的公式
nir-近InfraRed (Band8)
re - RedEdge (Band5)
我的守则：
library(raster)
library(RStoolbox)
re_path <- "D:/R/T43PHS_20190223T050811_B05.jp2"
nir_band <- "D:/R/T43PHS_20190223T050811_B08.jp2"
re <- raster(re_path)
nir <- raster(nir_band)
浏览 9提问于2021-03-15得票数 2
我有两个不同的空间分辨率和投影光栅图像。我想重采样多个粗图像(0.25度)，以更好的分辨率(4km)的基础上，更精细的分辨率，是4公里。
两者都在不同的文件夹中。我试着用R写代码，但我错了。
>library(raster)
>setwd("D:/mtech project/data/DEC 16/RH/1-DEC-16/extracted/")
>inFiles1 <- list.files(pattern="*.tif")
> nFiles1 <- length(inFiles1)
> setwd(
浏览 0提问于2018-11-10得票数 1
回答已采纳
1回答
我正在学习R tool，我也学习了一个例子，但是尽管我做的一切都一样，但在the plot raster过程中我得到了一个错误。
以下是代码：
library(raster)
library(rgdal)
myRaster1 <- raster(nrow=4, ncol=4)
myRaster1[]<- 1:ncell(myRaster1)
myRaster2=raster(nrow=8, ncol=8)
resample(myRaster1, myRaster2, method='bilinear')
plot(myRaster2, main="Rast
浏览 3提问于2017-04-06得票数 1
回答已采纳
1回答
首先，我成功地提取了每个多边形的平均栅格温度值，使用以下程序：
您可以在此链接上下载GIS层：
## load packages
library(raster); library(rgdal)
## Read rasters
ras_temp<-raster("ras_temp.tif")
plot(ras_temp)
ras_alti<-raster("ras_alti.tif")
## read polygon
polygon <- readOGR(dsn = getwd(), layer = "polygon
浏览 10提问于2017-05-03得票数 1
回答已采纳
2回答
我正在简化我的问题，所以希望这是有意义的。
我有三个正在处理的栅格。
library(raster)
LC <- raster("finalRRaster2018.tif")
nativeAET <- raster('kohala_aet_mm_ann_int.tif')
nonnativeAET <- raster('pred_AET.tif`)
栅格LC是具有landcover属性值的栅格，1表示原生landcover，0表示非原生landcover。
nativeAET和nonnativeAET分别是本地树种和nonnative树
浏览 1提问于2018-04-21得票数 0
2回答
我在尝试叠加多波段光栅。它们是通过欧空局哨兵1预处理工具SNAP创建的。每个tif文件都是逐层传输的。
我装载了两个栅格，并试图堆叠它们：
rs1 <- raster::stack("example/rs1.tif")
rs2 <- raster::stack("example/rs2.tif")
rsstack <- stack(rs1,rs2)
然后，我得到以下错误消息：
误差在compareRaster(x)中的不同程度
栅格的范围几乎是相同的：
> rs1
class : RasterBrick
dimensions
浏览 5提问于2020-08-13得票数 0
回答已采纳
1回答
我需要你的小帮助重采样多个文件具有相同的名称，但在不同的文件夹使用R。
我在reclass1中有一个文件夹，其图像名称的范围为1:172。我在reclass2中有另一个文件夹，其图像名称的范围为1:172。我需要用reclass2中的1.tif重采样reclass 1中的1.tif，并对所有172个文件进行类似的采样。我可以进行单次重采样，但不能在循环中完成所有操作。
我试过了：
setwd("D:/archived_121316_before_files/Reclass2")
rlist2 = list.files(getwd(), pattern="tif$
浏览 0提问于2017-02-12得票数 0
我想知道在raster包中对resample()函数做一些修改的可能性。首先，在"bilinear"方法中，默认情况下是分配最近的四个单元格的加权平均值，而我想改变为不同的最近单元格的数量，是可能的吗？其次，是否也可以创建mean方法来计算最近的n个单元的算术平均值？例如，在第一种情况下25个单元格：resample (myraster, myresolution, window=matrix (nrow=5,nol=5), method="bilinear")，在第二种情况下：resample (myraster, myresolution, window=
浏览 22提问于2020-03-23得票数 0
1回答
我有两个光栅:一个("log")，我想用来遮蔽另一个("roi")中的像素。
> log
class : RasterLayer
dimensions : 11957, 13462, 160965134 (nrow, ncol, ncell)
resolution : 30, 30 (x, y)
extent : 621465, 1025325, 8464135, 8822845 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=21 +south +datum=
浏览 0提问于2017-02-07得票数 0
回答已采纳
2回答
我在使用R中的DEM文件时遇到了问题，该文件可以作为ESRI网格adf文件(例如，包含各种adf文件的文件夹，如hdr.adf、w001001.adf等)。
我可以使用光栅包来读取和绘制文件：
dem<-raster("w001001.adf")
plot(dem)
结果的情节看起来很好，我可以做的事情，如裁剪文件。但是，我想从以下位置重新投影该文件：
+proj=laea +lat_0=-100 +lon_0=6370997 +x_0=45 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,
浏览 1提问于2020-02-14得票数 1
1回答
我正在写一个脚本，将采取任何三个栅格，并将它们裁剪到最小可能的程度。所有这三个光栅都将具有相同的分辨率和投影。但是，最小程度的裁剪不会更改这三个光栅的范围。我尝试过setExtent，同样的事情也发生了。如果有人能给出建议，我将不胜感激。以下是示例代码：
library(raster)
#Projection of all three rasters
newproj<- "+proj=utm +zone=4 +datum=WGS84 +units=m +no_defs +ellps=WGS84
+towgs84=0,0,0"
#Create three raste
浏览 0提问于2018-03-01得票数 1
1回答
我有一个具有不同分辨率的DTM和DSM。以下是每个栅格层的摘要。
> raster_dsm
class : RasterLayer
dimensions : 2001, 2501, 5004501 (nrow, ncol, ncell)
resolution : 0.5, 0.5 (x, y)
extent : -112500.2, -111249.8, 388999.8, 390000.2 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=16.333333333333
浏览 0提问于2019-05-06得票数 0
回答已采纳
1回答
我用两个光栅工作，每个光栅都有不同的分辨率。我想知道是否有一种更有效的方法来匹配更粗的光栅分辨率和更精细的光栅分辨率。现在，我正在使用掩码函数来节省一些时间，剪辑到正确的程度，并更改分辨率：
library(raster)
#the raster template with the desired resolution
r <- raster(extent(-180, 180, -64, 84), res=0.04166667)
# set some pixels to values, others to NA
r <- setValues(r, sample(c
浏览 6提问于2016-02-25得票数 2
回答已采纳
我目前正在处理一些来自瑞士的100x100米分辨率的数据。我想将数据从基于国家网格的千米方聚合到1x1千米。X= 500000，y= 130000)，与我的RasterLayer对象的原始范围不同。我提供了基于原始范围(=r)的RasterLayer和基于瑞士国家网格单元(=r.agg)的新范围的RasterLayer的一些代码，并在r上绘制了r.agg (后者，我将其涂成红色以使边界更明显)。 library(raster)
ncol <- 3677 # same as original RasterLayer
nrow <- 2261 # same as original R
浏览 33提问于2019-10-19得票数 1
回答已采纳
我有点问题，把.csv文件转换为光栅在R.我的.csv文件包含坐标(long和lat)、radius (在deg中)和站点类型。我能够将坐标转换为光栅，并能够使用st_buffer()绘制圆圈，但我面临两个问题：
我不能把圆圈变成光栅..。我试过使用rasterize()和fasterize()，但都没有工作，我得到的只是一个空光栅层，我似乎无法根据站点类型对坐标和圆圈进行分类。
知道我做错了什么吗？我该如何分类我的圆圈呢？提前谢谢你！
下面是我使用的代码：
> head(sp_csv_data)
Longitude Latitude Radius Site_Type
1 -1
浏览 1提问于2021-04-14得票数 0
回答已采纳
1回答
光栅图像似乎是用R的传单移动的。
、、、、
我想用R中的传单包来绘制一些空间数据，但是所生成的光栅图像似乎与参考网格相比发生了变化。我怀疑地图投影问题，但我不是这方面的专家，因此，任何帮助将不胜感激。
下面是绘制地图的最低限度代码：
library(leaflet)
library(sp)
library(raster)
set.seed(111)
# create dummy data -rectangular grid with random values
m = 10
n = 10
x = seq(45,48,length.out = m)
y = seq(15,18,length.out = n)
X = matrix(re
浏览 2提问于2015-11-16得票数 3
回答已采纳
1回答
有人能帮我把下面的代码写得更有效率吗？由于数据集非常大，每次运行提取部分(栅格、多边形)都非常缓慢。我认为它们可以结合在一起，所以只能完成一次，但我很难做到这一点。
result$mean <- extract(raster, polygons, weights=T, na.rm=F, fun=mean)
result$cells <- unlist(lapply(extract(raster, polygons, weights=T, na.rm=F),
FUN = function(x)(length(x[,1]))))
浏览 1提问于2018-05-25得票数 1
2回答
重采样光栅
、、、
我试图重新采样一个高分辨率(25米)和分类数据(1到13)的森林覆盖光栅到一个新的低分辨率(约1公里)的RasterLayer。我的想法是将森林覆盖数据与其他低分辨率光栅数据结合起来：
我尝试了raster::resample()，但是由于数据是绝对的，我丢失了很多信息：
摘要(as.factor(df$loss_year_mosaic_30m))0 1 2 3 4 5 6 7 8 10 11 12 12 3777691 65 101 50 151 159 295 291 134 102 126 104 91
正如你所看到的，新的光栅有所需的分辨率，但也有很多零。我认为这是正常的，因为我在
浏览 5提问于2016-06-21得票数 6
回答已采纳
2回答
我有一个全世界的光栅文件，它的详细信息是：
class : RasterLayer
dimensions : 18000, 43200, 777600000 (nrow, ncol, ncell)
resolution : 0.008333333, 0.008333333 (x, y)
extent : -180, 180, -60, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs
data source
浏览 3提问于2016-09-07得票数 0
1回答
我想使用terra包中的warp函数，但我的terra库中似乎没有该函数。这是一个bug吗？或者它可能与系统gdal的安装有关？我还没有在R之外安装gdal。
编辑:经过更多的变通之后，我猜函数terra::warp已经是较新的版本了，重命名为terra::resample。因此，可能只是帮助中存在不一致之处。
R version 4.0.2 (2020-06-22)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)
Matrix products: defaul
浏览 0提问于2020-08-30得票数 0
2回答
我想调整一个光栅的尺寸(在这种情况下放大)，以匹配另一个光栅的尺寸。使用disaggregate I可以放大光栅，但调整因子参数只接受整数，因此输出更接近目标尺寸，但不匹配目标尺寸。下面的示例通过尝试调整ras2大小以匹配ras1维度来说明：
> require(raster)
> ras1 = raster(volcano)
> ras2 = aggregate(raster(apply(volcano, 1, rev)), fact=2, FUN=mean) # aggregated & rotated
> extent(ras1) = extent(ras
浏览 5提问于2014-08-26得票数 0
回答已采纳
2回答
我正试图创建这个网格，“将人口网格划分为两个部分：(A)检测到照明的人口，(B)没有检测到照明的人口”。我有两个光栅数据，一个population.tif和一个nighttimelights.tif。基本上，我想满足两个条件：
人口>1，夜间灯> 20，这是A段人口>1，夜视灯< 20，这是段B，我想在R中创建条件语句，我到目前为止尝试的是：
library(raster)
ntl = raster("path/ntl_atprk.tif")
pop = raster("path/pop.tif"
浏览 0提问于2022-08-07得票数 0
回答已采纳
1回答
所以我试着做了一些我在这里读到的关于多个栅格的例子。我在r上使用随机森林来制作一些预测地图。到目前为止，我的地形协变量起作用了(相对坡度位置，saga湿度指数...)。我有三个新的图像添加红色，绿色和透明。在我执行了栅格堆栈命令后，它给出了这个错误"Error in compareRaster(x)：different extent“，下面是我使用的栅格和代码的数据信息。 library(raster)
catch <- raster("Catchment Slope.tif")
EM <- raster("EM_data_April_14.tif
浏览 42提问于2021-05-08得票数 0
1回答
如何用不同程度的光栅制作光栅砖？
、、、、
我是R的新手，所以这个问题是非常基本的，但是我一直在与它斗争，无法找到一个有效的解决方案。我想从同一地区的一些陆地卫星图像中创建一个栅格砖。它们是以HDF-EOS格式下载的，我使用Modis重新投影工具将它们转换为.tif。
由此产生的光栅具有相同的投影，但在其范围，分辨率和来源不同。
我尝试了几种方法，概述如下：
手动定义子集范围并删除所有栅格。然后试着用子集栅格做块
重新对栅格进行采样，给它们提供相同数量的列和行。理想情况下，这将确保栅格细胞是对齐的，可以放在栅格砖中。这个选项创建了一个砖块，光栅没有价值，它们是空的。
我想知道我应该遵循什么概念来修正这个范围。创建一个空栅
浏览 2提问于2013-12-22得票数 9
回答已采纳
我有一个名为AOT的文件夹，其图像名称的范围为1:744。我有另一个名为FOG的文件夹，其中有一个图像名称- '01DEC2016_0000_L2C_FOG_FOG.tif‘。我需要用雾中的'01DEC2016_0000_L2C_FOG_FOG.tif‘重采样AOT中的1.tif，并且类似于AOT文件夹中的所有744个文件。我正在做重采样与循环的帮助，但一个错误来了。 我试过了： > require(raster)
#reading files from folder AOT
> flname1 <- list.files(path = 'D:/mt
浏览 51提问于2019-02-07得票数 0
回答已采纳
1回答
我正在使用R中的asc文件，并尝试在其他人编写的脚本上通过一些GIS处理来运行它。我目前正在处理的代码块是 asc <- (raster(inputRas)) asc.extent <- extent(asc) head <- read.ascii.grid.header(inputRas) writeRaster(asc, "outputRas", format="SAGA", overwrite=TRUE) 它一直工作到最后一行。它在这里失败了，因为我的光栅是在直角坐标中，这是SAGA不喜欢的。我花了很多时间搜索和尝试不同的功能，但没有
浏览 13提问于2020-10-03得票数 0
我有两个不同范围的栅格，但是当我在R中使用重采样时，输出是一个空白的raster......How，我能解决这个问题吗？请参阅下面的一些代码，used....the重采样栅格参数的值中带有NA……
show(alt)
class : RasterLayer
dimensions : 1800, 4320, 7776000 (nrow, ncol, ncell)
resolution : 0.08333334, 0.08333334 (x, y)
extent : -180, 180, -60, 90.00001 (xmin, xmax, ymin, ymax
浏览 4提问于2018-05-18得票数 1
1回答
具有任意函数的非整数因子聚合光栅
、、、、
我想把人口栅格加起来1.5倍，把细胞的值加起来。
虽然aggregate()允许我在聚合时对值进行求和，但它的factor参数只接受整数值。projectRaster()和resample()允许我精确地调整分辨率，但(据我所知)我仅限于预先打包的双线性插值和最近邻计算方法。
是否有一种方法可以通过非整数因子聚合光栅并指定聚合时要使用的函数？
library(raster)
set.seed(10)
proj <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
r <- raster(resolution =
浏览 1提问于2017-03-12得票数 6
回答已采纳
我处理Sentinel2图像，并试图重新整理它们。
我尝试了以下代码：
import os, fnmatch
INPUT_FOLDER = "/d/afavro/Bureau/test_resampling/original"
OUTPUT_FOLDER = "/d/afavro/Bureau/test_resampling/resampling_10m"
def findRasters (path, filter):
for root, dirs, files in os.walk(path):
for fi
浏览 1提问于2019-04-24得票数 2
回答已采纳
1回答
我想将高分辨率的栅格重新采样到较粗的分辨率，但以这样的方式，即较粗的网格单元保留单元的最大值。
由于R的raster包中的resample函数中没有fun参数，因此我组合了一个简单的自定义函数：
resampleCustom <- function(r1, r2) {
resRatio <- as.integer(res(r2) / res(r1))
ret <- aggregate(r1, fact = resRatio, fun = max)
if (!compareRaster(ret, r2, stopiffalse = FALSE)) {
浏览 10提问于2018-01-23得票数 2
回答已采纳
1回答
我正在尝试计算R中DEM的地形湿度指数(TWI)。
对于TWI，我想使用dynatopmodel包中的upslope.area函数。运行函数时，我得到以下错误："Raster具有不同的x和y单元分辨率“。但是当查看我的栅格的属性时，x和y像元大小都是1。有人遇到过类似的问题并找到了解决方案吗？
下面是我的代码：
twi <- upslope.area(marion_dem, log-TRUE, atb=TRUE, deg=0.1, fill.sinks=TRUE)
Raster properties:
class : RasterLayer
dimensions
浏览 20提问于2019-10-25得票数 3
1回答
我有一个包含空白(NAs)的网格，我想使用插值来填充它。我的网格在x和y维上显示自相关，所以我想尝试双线性插值。我发现的大多数解决方案都集中在“过采样”(为了增加样本数量/网格大小而进行插值)，但我不想/不需要更改网格大小。我只想用内插法来填充NAs。其他可能的解决方案，或是，而不是双线性插值，或简单的。
我发现，使用光栅包，我可以输入包含NAs的网格(作为光栅)，并使用“重新采样”命令输出相同大小的网格。然而，结果看起来像最近邻插值，而不是双线性插值。
我是不是遗漏了一些东西，有办法用光栅包做双线性插值？还是有更好的方法来做双线性插值，简单地填充NAs？
library(raster)
#
浏览 0提问于2018-06-16得票数 8
回答已采纳
2回答
我只想在y轴上有平面图边距。默认情况下，边距出现在x轴和y轴上。下面是一个示例：
library (raster)
library(rasterVis)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
levelplot(r, margin=T)
浏览 1提问于2015-08-27得票数 2
回答已采纳
我是c++和im的新手，他们试图编写一个简单的c++包装器来集成这个第三方c++库；这个库附带了糟糕的文档。
要集成的函数(这是该功能附带的所有文档)：
virtual ImageCoord Raster::groundToImage(const XyzCoord & groundPt,
double desiredPrecision=0.001,
double* achievedPrecision=null) const
函数文本:此方法将给定的groundPt (x、y、z以米为单位)转换为返回的图像坐标(行，样本在完整
浏览 7提问于2022-08-19得票数 -2
1回答
我的问题是如何计算R中两个光栅之间的角度？下面是一个可复制的示例和我得到的错误：
angle <-function(vectora, vectorb){
theta <- acos( sum(vectora*vectorb) / ( sqrt(sum(vectora * vectora)) * sqrt(sum(vectorb * vectorb)) ))
}
raster1 <- raster(ncols=150, nrows=150, xmn=0)
raster1[] <- rnorm(150*150,ncell(1))
raster2 <-
浏览 4提问于2013-10-25得票数 1
回答已采纳
2回答
我有一组全球定位系统数据，我试图使用“内核平滑”包中的bkde2d来平滑内核。我在'ks‘包中使用了Hpi带宽估计器来确定我的带宽，但是当我运行内核平滑并将结果列表转换为光栅时，结果产品似乎具有不同的x，y分辨率，因此无法作为ascii导出。当作为GRD文件导出时，也不可能将这个光栅读入GIS工具，因为它似乎是损坏的，大概是因为分辨率不同。
下面是我运行的一些示例代码。我的数据是在UTM30，WGS84中预测的：
bnd=Hpi(x=cbind(GPS$lon, GPS$lat))
coord <- cbind(GPS$lon, GPS$lat)
est <- bkde2
浏览 2提问于2016-03-02得票数 1
回答已采纳
我使用一个函数来根据ID计算每个单元格的行串长度，并将其存储在一个列表中，将列表中的每个元素转换为一个RasterLayer，并将该列表转换为一个RasterStack，平均所有层并得到一个栅格。
#function
# build_length_raster <- function(one_df) {
intersect_list <- by(
one_df ,
one_df$sub_id,
function(subid_df) sf::st_intersection(grid2, subid_df) %>%
dplyr::
浏览 3提问于2022-04-21得票数 0
回答已采纳
1回答
如何在ggmap上叠加透明栅格？
、、、、
我正在尝试找到一种方法来添加alpha函数，以便在ggmap上覆盖raster layer。这里有一些关于的建议，但我还是想使用raster file。是否有任何方法可以添加alpha并使raster图层成为transparent
library(raster)
library(ggmap)
location <- get_map(location = c(lon =22, lat =40), zoom = 6, maptype = "hybrid")
# get the raster layer
tmin <- raster::getData('world
浏览 3提问于2018-02-24得票数 0
我试着在我的光栅砖上做一些有条件的逻辑。在下面的代码中，myFun1正确地生成raster.out1。但是myFun2在试图生成raster.out2时会生成一个错误。错误信息是
错误(测试)：对“哪个”不符合逻辑的参数
除了使用值为5的变量之外，这两个函数看起来是相同的。我明显漏掉了一些东西。
library(raster)
raster.in <- raster(nrows=100, ncols=100)
raster.in[] <- runif(ncell(raster.in), min = -10, max = 10)
const1 <- 5
myFun1 &l
浏览 7提问于2020-04-30得票数 2
回答已采纳
我有一个具有大量属性的raster对象，我希望用一个属性来绘制R中的空间数据，并对其进行颜色编码。我还没有想出如何使用特定属性的信息来实现这一点。到目前为止，我已经成功地使用factorValues()提取了所选择的属性，但我现在无法确定如何将这些信息合并到plot()函数中。我尝试使用光栅包文档中提到的ratify()和level()函数，但我不明白如何将简化的在线示例用于具有多个属性的光栅。
任何关于如何实现这一目标的建议都将不胜感激。
# read in shapefile
shp = readOGR(".", "grid")
#convert to
浏览 0提问于2015-02-19得票数 3
回答已采纳
在下面的栅格中，NA值增加了。
library(raster)
filename <- system.file("external/test.grd", package="raster")
r <- raster(filename)
r[r<300] <- NA
summary(getValues(r))
是否有可能只填入NA细胞？我一直在使用，但如下所示，NA值仍然保留在最终产品中。
fill.na <- function(x, i=5) {
if( is.na(x)[i] ) {
return( round(mean
浏览 1提问于2017-08-11得票数 6
回答已采纳
1回答
光栅的合并(镶嵌)改变分辨率
、、、、
我正在使用我开发的R脚本合并两个MODIS块，这些是产品：
因此，我从相同的日期(2019180)打开两个产品(瓷砖h15v05和h16v05)，然后打开每个h16v05并将它们合并在一起(从h15v05到h16v05的00h等等)
对这两种产品的总体(使用合并选项)的可视化：
紫色方格是分隔两块瓷砖的分界线的位置。
使用我的代码，我获得了一个具有不同分辨率(以及不同的min/max值)的像素的图，我不明白为什么：
我怀疑所取得的结果是由于：
1-由正弦波型向长端型WGS84型转变；
2-使用重采样(方法ngb)与马赛克。
我的代码很广泛，但是下面是其中的一些部分：
# O
浏览 5提问于2021-07-23得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
