R:可以仅从栅格层的非NA区域采样随机点吗？
A: 在栅格数据处理中，非NA区域指的是不包含无效或缺失值的区域。采样随机点是一种常用的方法，用于从栅格层的非NA区域中选择一些点进行分析或表示。
优势：
- 采样随机点可以减少数据量，提高计算效率。
- 通过采样随机点，可以更好地代表整个非NA区域的特征，避免过多的重复信息。
应用场景：
- 地理信息系统（GIS）中，通过采样随机点可以获取栅格数据中非NA区域的空间分布特征。
- 环境科学中，可以利用采样随机点来研究栅格数据中非NA区域的环境变量分布情况。
- 农业领域中，可以利用采样随机点来分析栅格数据中非NA区域的土壤质量、植被分布等。
推荐的腾讯云相关产品和产品介绍链接地址：
腾讯云地理信息服务（Tencent Cloud GIS）是一项提供地理信息处理和分析的云服务，可以用于处理栅格数据和采样随机点的需求。具体产品介绍和功能可以参考腾讯云官方文档：腾讯云地理信息服务
请注意，以上答案仅供参考，具体的产品选择和使用应根据实际需求和情况进行评估和决策。
相关·内容
1回答
我有几个图层的堆叠栅格数据集，但是，我想计算每个单元格的总和，用于不同的图层选择，并最终生成一个新的图层，谁有一些好的建议，使用计算或叠加或一些其他栅格计算在R？我可以通过循环进行计算，但当我有很多层时，它会消耗很多次，也会使用很多存储空间，我的脚本如下所示：
## library(raster)
make_calc <- function(rr, start, end) {
rr <- as.array(rr)
start <- as.array(start)
end <- as.array(end)
dms <- dim(raster)
t
浏览 1提问于2018-04-12得票数 0
1回答
从网格像元内的栅格中提取随机点
、、、、
我想要从每个网格单元内的光栅的随机坐标中提取非NA值。
栅格示例
library(raster)
r <- raster(ncol = 10, nrow = 10, xmx = -80, xmn = -150, ymn = 20, ymx = 60)
values(r) <- runif(ncell(r))
网格的一个示例
grid <- raster(extent(r))
res(grid) <- 15
proj4string(grid)<- proj4string(r)
gridpolygon <- rasterToPolygons(grid)
plo
浏览 12提问于2017-12-27得票数 2
回答已采纳
1回答
R中的坐标值
、、
我是R的新手，正在尝试学习它是如何用于一些物种分布建模的。我正在尝试使用maptools包中包含的wrld_simpl地图创建意大利的蒙版。当我运行代码时，当我试图创建spxy时，我得到"Error in .local(obj，...)：NA values in create“。我怀疑这个问题可能根源于光栅化步骤，但我不确定……我做错了什么？
对代码转储表示歉意
# Create library
library(raster)
library(dismo)
library(sf)
library(maptools)
library(rgdal)
library(sp)
library(
浏览 0提问于2020-09-22得票数 0
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 3
回答已采纳
2回答
我想在我的shapefile中添加一个z字段。Z值是排水沟中心的高程。我不知道我该怎么做。质心的x和y值在表中。我有这个地区的dem。 ?
浏览 14提问于2019-04-18得票数 0
因此，我有一条常见的错误消息，其中图像的行数和列数不匹配(在交叉引用中)。我已经概括了我的一个图像，然后使用展开使分辨率再次匹配。然而，在这个过程中，我丢失了一些列(这并不困扰我)，但是，我不知道该怎么做才能使我的两个图像再次保持相同的大小。有人能帮我吗？
非常感谢你。
浏览 2提问于2014-05-29得票数 0
我想在两个分类光栅的每一个可能组合下计算面积。我使用以下代码
library(terra)
#First create two rasters
r1 <- r2 <- rast(nrow=100, ncol=100)
#Assign random cell values
set.seed(123)
values(r1) <- runif(ncell(r1), min=0, max=1)
values(r2) <- runif(ncell(r2), min=0, max=1)
# classify the values into two groups
m_r1 &l
浏览 7提问于2022-04-23得票数 0
回答已采纳
我有11个大光栅文件，里面有大约314578920个细胞。这些细胞中有些含有NA。如果其他栅格在特定的单元格中有价值，我想将NA替换为零。但是所有栅格中含有NA的细胞应该是相同的。例如，请参见下面的图像：
我创造了五个栅格(上排)来说明我的问题(但实际上我有11个栅格)。输出光栅应该像底部的行光栅。
将所有NA替换为零(即使用r[is.na(r[])] <- 0)使得光栅太大，内存无法处理。或者，省略NA并不能达到我的目的。任何关于如何解决这一问题的想法都将受到高度赞赏。
样本代码：
library(raster)
r <- raster(nrow=5, ncol=5)
浏览 0提问于2018-01-18得票数 0
回答已采纳
2回答
我试图给一些站点数据赋值，这些数据就在我有天气数据的区域之外。我试图根据最近的单元格值提取，如果可能的话，在40公里内提取一个单元格值。
我的光栅看起来是这样的：
class(r)
class : RasterBrick
dimensions : 201, 464, 93264, 23376 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : -40.5, 75.5, 25.25, 75.5 (xmin, xmax, ymin, ymax)
coord. ref. : +p
浏览 3提问于2014-10-30得票数 2
回答已采纳
我正在处理一个由20个层组成的数据集，这些层堆叠在一个RasterBrick中(源自一个数组)。我已经查看了用'calc‘和'cellStats’计算的层的总和。我已经使用calc来计算总值的总和，并使用cellStats来查看每层值的平均值(对于时间序列很有用)。然而，当我对每一层的平均值求和时，它是另一个计算和的值的一半。造成这种差异的原因是什么？我忽略了什么？ 代码如下所示： testarray <- runif(54214776,0,1)
# Although testarray should contain a raster of 127x147 with
浏览 125提问于2020-04-20得票数 1
回答已采纳
我已经编写了SAFT (合成孔径聚焦技术)的代码，并成功地绘制了2D图像。现在，我在3D中绘制它有一个问题，因为我不知道我应该使用什么类型的函数。我想创建3D图像，如下面的图像链接所示(右侧的图像具有实体体积)，但我对他们使用的功能一无所知。我所知道的就是左边的图片是拼接出来的。
如果有人能帮我解决这个问题，我将不胜感激。提前谢谢你！
这是我在3D数组中使用的示例代码。请注意，我只使用了强度的随机值，这里是V，因为实际值存储在excel文件中。我在使用等值面时遇到错误，因为我的数组不是使用meshgrid创建的。
r=linspace(50,500,101)';
o=zero
浏览 1提问于2014-05-03得票数 1
6回答
我正在尝试实现。我不能理解一件事:我检查的所有实现都选择图像上的随机点来计算粒子逃逸的路径。他们为什么要这么做？为什么不遍历所有像素呢？
随机点的作用是什么？点越多，图片越好，所以我认为所有像素都是最好的图片--我错了吗？
根据我的测试数据：
处理400x400的图片。所以如果我遍历的话，需要160000个像素来迭代。
使用随机抽样，图片只有在一百万个点之后才开始成形。好的结果显示了大约10亿个随机点，这需要几个小时的计算。
浏览 3提问于2009-09-29得票数 15
回答已采纳
1回答
我试图在R中的特定列中用随机小数替换NAs，但是R生成具有相同尾随分数的随机小数，并且只是在十进制之前更改部分。以下是我尝试过的方法：
df_LT$ATC[is.na(df_LT$ATC)] <- sample(seq(10.2354897,23.78954214), size=sum(is.na(df_LT$ATC)), replace=T)
dplyr
df_LT <- df_LT %>%mutate_at(vars(df_LT$ATC), ~replace_na(., sample(10.2354897:23.78954214, size=sum(is.na(ATC
浏览 0提问于2019-10-30得票数 0
回答已采纳
2回答
如何检查光栅和多边形是否重叠？
、、、、
我有一个栅格，我想检查有值的单元格是否与我感兴趣的区域重叠，这是一个空间多边形文件。
有什么快速的方法吗？
我的光栅看起来是这样的：
我试过的是：
r <- raster("C:/user/test.jp2")
studyarea <- readOGR("C:/user/boundary.shp")
dummi <- as(extent(r), "SpatialPolygons")
if (gContainsProperly(studyarea, dummi)) {
print("Raster extent
浏览 3提问于2020-09-23得票数 1
回答已采纳
2回答
我正在使用rgdal包在R中编写一个图像分类脚本。所讨论的光栅是一个具有28个通道的PCIDSK文件:一个训练数据通道、一个验证数据通道和26个光谱数据通道。目标是填充一个数据帧，该数据帧包含训练数据通道中非零的每个像素的值，以及26个频带中的相关谱值。
在Python/Numpy中，我可以轻松地将整个图像的所有波段导入到多维数组中，然而，由于内存限制，R中唯一的选择似乎是逐块导入此数据块，这非常慢：
library(rgdal)
raster = "image.pix"
training_band = 2
validation_band = 1
BlockWidth = 5
浏览 2提问于2010-11-16得票数 2
回答已采纳
1回答
我的光栅文件具有相同的分辨率和范围，但NA的数量不同。我想统一它们之间的NA数。如果一个单元格在所有光栅文件中都不是NA，是否可以通过将其视为非NA来实现？
下面是一个示例：
library(raster)
library(terra)
f <- system.file("external/test.grd", package="raster")
r1 <- raster(f)
r2 <- calc(r1, fun=function(x){ x[x < 500] <- NA; return(x)} )
r1 <- ca
浏览 4提问于2021-05-19得票数 1
回答已采纳
1回答
在R中按属性提取栅格
、、、
我一直在尝试创建一个新的栅格对象，它只包含现有栅格中的几个值。我使用的是这里找到的类栅格：。
class : RasterLayer dimensions : 14902, 19161, 285537222 (nrow, ncol, ncell)
resolution : 0.002349, 0.002349 (x, y)
extent : 110, 155.0092, -45.0048, -9.999999 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs
浏览 0提问于2012-03-03得票数 0
回答已采纳
我是一个MATLAB的业余用户，我想有一个4D颜色图，其中的颜色代表第四维度。我尝试使用scatter3，如下所示：
fid = fopen('abc.dat','r');
datacell = textscan(fid, '%f%f%f%f');
fclose(fid);
all_data = cell2mat(datacell); %converted into a matrix
boxx = all_data(:,1);
boxy = all_data(:,2);
boxz = all_data(:,3);
disden = all
浏览 0提问于2015-10-05得票数 3
我正在计算一个新的栅格(输出ras)基于2个栅格(输入ras)和一个“阶层”光栅。Stratum光栅值(1至4)指偏置和权重数据中的行。地层值'4‘被用来填充地层栅格中的任何'NA’，否则该函数将崩溃。需要输入以下内容。
# load library
library(raster)
# reproducing the bias and weight data.frames
bias <- data.frame(
ras_1 = c(56,-7,-30,0),
ras_2 = c(29,18,-52,0),
ras_3 = c(44,4,-15,0)
)
rowna
浏览 3提问于2017-04-13得票数 0
回答已采纳
我最近访问了一个网站
这给了我任何地球位置的每个点的海拔数据。现在，我想在C#中构建一个桌面应用程序，使用google或其他方式，在这个应用程序中，用户将搜索该位置并显示该位置的图像，用户可以在该位置放大或缩小并选择所需的区域。选择后，该工具将获取选定区域内各点的高程数据，并在该点上创建相应的颜色，从而使最高点具有白色，最低点将具有黑色，然后将数据转换为图像格式。因此，如果我能做到这一点，我基本上可以在几秒钟内创建高度图。
有人能告诉我如何做到这一点吗？我不知道如何获得选定区域内每个点的高程数据，以及如何从这些数据中计算和创建彩色图像。
谢谢。
浏览 3提问于2014-04-09得票数 0
回答已采纳
我正在使用R中具有以下结构的dataframe：
summary(data)
X Y Z type
Min. : 2.94 Min. :22.48 Min. :1.000 Length:225176
1st Qu.: 22.39 1st Qu.:48.17 1st Qu.:1.000 Class :character
Median : 29.64 Median :52.19 Median :1.000 Mode
浏览 3提问于2022-07-14得票数 0
1回答
如果一个程序在一个分辨率A的显示器上显示一个像素的X，Y，我能准确地预测在什么坐标下，同一像素将显示在分辨率B？
更多信息
2种显示分辨率如下：
A->1366x768
B->1600 x 900
将每个方向的最大分辨率除以后产生：
X方向标度因子= 1600/1366 = 1.171303075
Y方向标度系数= 900/768 = 1.171875
例如，假设显示A上唯一的红色像素发生在像素(1,1)处。如果我仅仅使用这些因子来放大，那么在显示B上，红色像素将显示在像素(1.171303075，1.171875)上。我不知道该如何解释，因为我习惯于把像素看作整数值。如果我知道屏
浏览 1提问于2014-08-08得票数 3
我有一个大小为n_X * n_Y的离散网格，其中中心单元位于(0,0)。在这个网格中，在(x_target, y_target)上有一个相对于中心的随机点，它可以随着时间的推移改变它的位置。
我想根据每个单元格与随机点的接近程度，给每个单元格一个介于0和1之间的权重。最远的点应该是0，最近的点应该是1。中间单元应该基于它们的距离在两者之间进行插值。
系统：
除了知道随机点的位置外，你知道我如何估计每个单元格的权重吗?假设最远的点是0，最近的点是1。谢谢。
浏览 0提问于2021-09-27得票数 0
1回答
我正在使用QGIS，并试图在栅格层中每个像素的中心生成均匀的空间点。我能够生成均匀分布的点，但无法使它们足够靠近或位于每个像素的中心。每个像素应该是1.85m，并且在UTM中的投影CRS上。我没有访问ARIGIS的权限。在QGIS中有没有办法做到这一点？甚至是在R中？我尝试生成一个网格以及生成点，但默认距离是0.0001，当我试图变小时，它会自动恢复。考虑到我在UTM，我也不确定如何解释这个默认值。
浏览 6提问于2018-01-31得票数 1
1回答
我正在尝试学习多重重要性抽样(MIS)，我需要根据每种抽样技术的PDF来计算样本的权重。我理解来自BRDF的PDF，以及基于光的半径命中特定球面光的PDF。但是对于网格发射器，我如何计算命中发射器的PDF？只是发射三角形的立体角度吗？
浏览 1提问于2016-09-03得票数 0
1回答
我有问题要把这个空间聚合数据转换成光栅。当我进行转换时，光栅有NA作为它的值。如下所示：
DL3
[1]
class : SpatialPolygonsDataFrame
features : 126
extent : -15.04001, 46.1036, 3.759985, 31.71804 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +towgs84=0,0,0 +ellps=WGS84
variables : 1
names : LFRP
m
浏览 0提问于2015-07-09得票数 2
回答已采纳
4回答
分区统计QGIS
、、、
是否有任何用于分区统计工具(ArcGIS Spatial Analyst)的开源替代工具？什么是最好的工具(我可以在脚本中使用)来计算栅格文件的统计数据？
浏览 2提问于2011-04-27得票数 5
2回答
我需要使用空间数据框点从光栅中提取高度值。
通过两个对象(光栅和矢量)，我刚刚获得了NA。
> lat41_42_sp
class : RasterLayer
dimensions : 1201, 10801, 12972001 (nrow, ncol, ncell)
resolution : 0.0008333333, 0.0008333333 (x, y)
extent : -9.000417, 0.0004166633, 40.99958, 42.00042 (xmin, xmax, ymin, ymax)
crs : +proj=longla
浏览 0提问于2019-06-19得票数 0
2回答
我试图准备的环境层(在特定的亚洲地区)，用于Maxent模型。但是，我在最后一行中遇到了一些错误消息：
library(sp)
library(maptools)
library(rworldmap)
library(dismo)
# A specified range of Asia area that suitable for special species
tsta <- read.csv('CM10_Kop_Shp_V1.2/Asiaclip/Asiaclipt.csv',as.is=TRUE)[https://drive.google.com/file/d/0B
浏览 0提问于2015-12-25得票数 1
回答已采纳
1回答
我有一个NetCDF文件的全球海洋学(OmegaA)数据，在相对粗糙的空间分辨率，33个深度。我还有一个分辨率更高的全球测深光栅。我的目标是从OmegaA文件中获取海底NetCDF数据，使用测深数据来确定所需的深度。到目前为止我的代码；
library(raster)
library(rgdal)
library(ncdf4)
# Aragonite data. Defaults to CRS WGS84
ncin <- nc_open("C:/..../GLODAPv2.2016b.OmegaA.nc")
ncin.depth <- nc
浏览 2提问于2018-01-04得票数 1
回答已采纳
1回答
我有一堆植被趋势栅格，在1982-2013年间每年都有一层。此外，我在矢量中有一些数据，我想将其与栅格进行比较。但是，这个矢量数据并不完整，有些年份没有值。
有没有办法只用我的vetor数据有值的年份对应的栅格做一个新的堆栈？我有180个向量，所以如果我可以避免为每个向量手动执行此操作，我会很高兴。
我不知道怎么做，所以如果你有任何建议，请分享。
下面是我的数据的一个示例：
#yearly vegetation trend residuals
res <- list.files(pattern="*res.tif")
s <- stack(res)
> s
浏览 0提问于2017-04-26得票数 2
我正在试图提取关于未被烧毁的地区内不同植被类型的百分比的信息。然而，我为公园提供的数据包括在NAD_83投影中作为栅格文件被烧毁的区域，以及WGS_84投影中作为多边形层的植被覆盖。本质上，我试图抹去被烧毁的地区和植被层之间的重叠，只看那些没有被烧毁的地区的植被覆盖类型。植被被分解成12类标签，基于下层/灌木覆盖类型，我想要一个百分比的每种植被类型发生在未被烧毁的地区。有人能帮到这个转变吗？
我可以访问ArcGIS桌面10.8.2与高级许可证。
我试着把栅格转换成有效的点，但我不知道如何将点转换成多边形，使之与原始光栅的大小/形状完美匹配。除了尝试将光栅转换为具有相同坐标系的多边形外，我不知道
浏览 11提问于2022-11-14得票数 1
1回答
最近我遇到了一篇关于使用反卷积网络进行语义分割的论文：用于语义分割的学习反褶积网络。
网络的基本结构如下：
📷
最终的目标是生成一个概率图。我很难搞清楚如何实现反褶积层。文件中说：
非池层的输出是一个扩大的，但稀疏的激活映射。反褶积层将通过使用多个学习滤波器的卷积类运算解池获得的稀疏激活紧密化。然而，与将滤波器窗口内的多个输入激活连接到单个激活的卷积层相反，反卷积层将单个输入激活与多个输出相关联。反褶积层的输出是一幅放大密集的活化图。我们裁剪扩大的激活映射的边界，以保持输出映射的大小与前一个解池层的大小相同。反卷积层中的学习滤波器对应于重构输入对象形状的基。因此，类似于卷积网络，采用一种分层的
浏览 0提问于2015-11-24得票数 5
回答已采纳
3回答
lat long
7.16 124.21
8.6 123.35
8.43 124.28
8.15 125.08
考虑这些坐标，这些坐标对应于测量降雨量数据的气象站。
R中的gstat包简介使用meuse数据集。在本教程中的某一点：，这些人在这行代码中使用了"meuse.grid“：
data("meuse.grid")
我没有这样的文件，我不知道如何创建它，我可以使用这些坐标创建一个吗？或者至少给我一些讨论如何为自定义区域创建自定义网格的资料(即不使用GADM的行政边界)。
可能用错了字眼，甚至不知道这个问题对精明的人来说是否有意义。尽管如此，我
浏览 4提问于2017-04-16得票数 9
回答已采纳
在下面的栅格中，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回答
如果我有一个二值图像和一个不规则的凸多边形，我如何计算它们是否相交？多边形的坐标以图像的形式描述。
从碰撞检测或填充算法的角度来看，我对此有一些想法，但我认为这两种方法都不是最优的。我确信有一个经过试验和测试的方法可以解决这个问题，但是想不到关键字。
下面是我的意思的一个例子：
在这种情况下，它应该返回true。
浏览 0提问于2012-03-10得票数 1
回答已采纳
这是我第一次在gstat中使用协同克里格法。我的问题是，当感兴趣的变量和辅助变量不是在相同的位置测量时，我不确定如何准备提供给协同克里格法的数据帧。
包含感兴趣变量z_ar的数据框为voi_gs
> head(voi_gs)
x y z_ar
8974 312216.6 530439.8 49.03470
8283 312084.6 530559.8 57.15355
5057 311772.6 530883.8 49.13453
1551 311976.6 531207.8 44.73679
10116 312168.6 5301
浏览 6提问于2018-04-26得票数 0
1回答
我有个问题。我必须建立一个程序，用蒙特卡罗方法计算一个立方体的体积。立方体在XYZ轴上的起点为(0,0,0)。也可以有一些球体。球体可以是立方体，那么立方体的体积就是:=立方体体积。- sphere vol.此球体也可以位于立方体外部。我知道这种方法是如何工作的，但我在定义随机点的间隔时遇到了问题(我只知道用这种方法计算积分)。参数是:a-正方体边的长度，点-球体中心的(X，Y，Z)点，r -radius。例如a= 1，point1 = (0,0,0)，r1 = 0.5，point2 = (1.25,1.25,1.25) r2 = 0.25。在这种情况下，如何确定随机点的区间？
浏览 3提问于2015-03-15得票数 0
我正在尝试在R中创建物种分布模型。我已经在ArcMap中创建了栅格层，并将它们导入到R中。它们不能堆叠，除非范围完全相同，并且它们都具有相同的行数和列数。
然而，当我改变这些因素来成功地堆叠它们时，它们会丢失所有的值，并且我堆叠的数据帧只填充了NAs。
有谁知道我如何改变栅格图层的范围和分辨率，以便可以成功地堆叠它们--这样我就可以将环境信息附加到存在点。
干杯
浏览 2提问于2015-08-27得票数 3
1回答
最新的OpenGL规范(4.6core)第14.5.1章(基本行段栅格化)(将在上找到)指定如何对别名行(由两个浮点端点指定)进行栅格化。
理想情况下，GL使用“菱形退出”规则来确定通过光栅化线段生成的片段。对于以窗口坐标x_f和y_f为中心的每个片段f，定义一个菱形区域，即四个半平面的交点：
本质上，从p_a开始到以p_b结尾的线段产生的片段f与R_f相交，除非p_b包含在R_f中。
然而，考虑到这一情况：
当正在绘制的直线通过上述公式所描述的两颗相邻钻石的角时，任何碎片都不应被光栅化。表达式|x - x_f| + |y - y_f|的最小值是1/2，1/
浏览 3提问于2017-12-28得票数 4
回答已采纳
2回答
我已经为这件事困惑了好几天了。我有一个栅格，我需要根据sf对象中的多边形属性进行编辑。这两个对象具有相同的范围。 填充了9的栅格 r <- raster( nrow=10, ncol=10, xmn=0, xmx=10, ymn=0, ymx=10 )
values(r) <- 9
> r
class : RasterLayer
dimensions : 10, 10, 100 (nrow, ncol, ncell)
resolution : 1, 1 (x, y)
extent : 0, 10, 0, 10 (xmin, xmax, ymin,
浏览 20提问于2020-01-11得票数 2
回答已采纳
2回答
我在R工作。 我有一个数据框列，里面有母的或公的，还有一些NAs。 现在，我想为这一列中的NA值随机分配女性或男性。我不希望所有的NA都是男性或女性，但每个NA都随机分配或。 我该怎么做？ 最好的，corkinabottle
浏览 34提问于2020-12-29得票数 0
回答已采纳
1回答
查找重叠散点图和下采样散点图
、、、、
我有大量的数据，这些数据由点(x，y)组成，需要用散点图表示。 为了找到相同点之间的重叠并消除看不见的差异，我编写了以下代码： void DownSampler::makeDownSample(QVector<double> keys,QVector<double> values, QVector<int> pixmapIdx, QSize resolution, sViewSize view)
{
m_calculating = true;
QTime now = QTime::currentTime();
QMap<QPair<double
浏览 27提问于2021-01-20得票数 0
回答已采纳
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
回答已采纳
图片比文字好，所以请看一看
我拥有的是
一个RasterLayer对象(这里填充了随机值，只是为了说明，实际值并不重要)
一个包含大量多边形的SpatialPolygons对象
您可以使用以下代码重新创建用于图像的示例数据：
library(sp)
library(raster)
library(rgeos)
# create example raster
r <- raster(nrows=10, ncol=15, xmn=0, ymn=0)
values(r) <- sample(x=1:1000, size=150)
# create example
浏览 5提问于2016-11-17得票数 5
回答已采纳
我正在尝试使用ArcPy脚本将多个光栅文件转换为NetCDF文件。当我运行下面的脚本时，我得到以下错误消息：
Message File Name Line Position
Traceback
<module> <module1> 19
RasterToNetCDF C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\md.py 253
ExecuteError: Failed to execute. Parameters a
浏览 69提问于2017-01-19得票数 0
回答已采纳
2回答
我正在评估在2D自上而下的程序游戏中生成岩石的“块状”或“脉”的方法。显然，最常见的建议是使用各种类型的噪音，但我觉得噪音(单独)无法达到我想要的效果。
我使用了voronoi和perlin噪声的组合，因为它更容易生成不同噪声的“区域”。比起单独使用“高度图”，我可以更容易地把这些区域变成岩石或其他东西的“脉”--高度图往往更像漩涡/云图。
📷
有了一个高阈值和一些四舍五入，我可以把它变成一个对静脉有效的面罩：
📷
改变间距和比例就像玩数学/阈值一样简单。
困难在于选择哪个“岩石”(或“矿石”)将填补每个脉。这主要是因为我的岩石是加权的。稀有岩脉不应常见。
我目前把最初的voronoi细胞值映
浏览 0提问于2017-11-09得票数 5
我正在处理关于树木物种和生境的数据，两者都以矢量格式存在(点和多边形)。每个点发生在多边形的特定网格单元内，该多边形具有描述该网格单元的主要土地使用的属性。
我的目标是使用树形点数据创建另一个数据集。理想情况下，该数据集将考虑到当前的生境类型(网格单元格)，并在特定距离内搜索与该类型匹配的所有单元格。最后，一个新的点将被随机放置在与原始类型匹配的网格单元中。
这在QGIS中是可能的吗？
浏览 2提问于2017-10-03得票数 0
1回答
我有一个站点ID和舞台高度的点形状文件。我想要创建一个栅格，其中每个单元的阶段高度值(以米为单位)最近的原地站，该细胞。
我想让这个光栅和另一个光栅相匹配。因此，如果我能够输入我创建的栅格(下面描述的数据集3)和我的点shapefile (1)，我想要它。
数据集：
1)具有河三角洲阶高的点形状文件
( 2)河三角洲范围的形状
3) NA's代表陆地的三角洲栅格图(如果需要的话也可以是0)，1是水，两个数据集10米分辨率和30米分辨率。
我有一个概念上的问题，就是我拥有的小溪的数量。例如(下图所示)，站1(蓝色圈)在技术上更接近黑色x区域，而不是站2(红圈)，但红色的舞台高度值更能代表
浏览 1提问于2017-01-03得票数 0
2回答
计算漫反射(朗伯反射)
、、、、
我需要帮助找出如何计算(朗伯)漫反射-即当光线击中表面，然后以随机方向反弹。因此，如果我有一个来自光源的光向量L，一个三维点X和一个点X处的法向量N，我如何计算随机反射的光线？
在我读的一本书中，他们说使用这个等式：
Wd =(θ，φ)= cos-1(sqrt( rand1 )，2*pi* rand2 )，其中Wd是反射射线，rand1和rand2是0,1范围内的随机数。
这对我不起作用。我使用了这个方程，然后将球面坐标转换为笛卡尔坐标，但这使得光线总是在相同的几个方向上反射，而不管法线的方向。
如有任何帮助，我们不胜感激！
浏览 6提问于2013-01-23得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
