如何在r中创建栅格的随机点(1000)
在R中创建栅格的随机点（1000），可以使用sp和raster包来实现。以下是一个完善且全面的答案：
在R中创建栅格的随机点（1000）可以通过以下步骤实现：
- 首先，确保已安装
sp和raster包。如果没有安装，可以使用以下命令进行安装：
install.packages("sp")
install.packages("raster")- 导入所需的包：
library(sp)
library(raster)- 创建一个空的栅格对象，指定栅格的行数、列数和分辨率：
r <- raster(nrow=100, ncol=100, xmn=0, xmx=10, ymn=0, ymx=10)这将创建一个100x100的栅格对象，范围从0到10。
- 使用
sampleRandom()函数在栅格上生成随机点：
random_points <- sampleRandom(r, 1000, sp=TRUE)这将在栅格上生成1000个随机点，并将其存储在random_points变量中。
- 可以使用
plot()函数将随机点绘制在栅格上进行可视化：
plot(r)
points(random_points, col="red")这将在栅格上绘制随机点，并将其标记为红色。
至于栅格的随机点的应用场景，它可以用于模拟自然资源分布、环境研究、地理信息系统等领域。例如，在生态学研究中，可以使用栅格的随机点来模拟动物的分布情况，进而进行种群数量估计和生态系统评估。
推荐的腾讯云相关产品是腾讯云地理信息系统（GIS），它提供了丰富的地理信息处理和分析功能，可以与R语言进行集成。您可以通过以下链接了解更多关于腾讯云GIS的信息：
希望以上信息对您有所帮助！
相关·内容
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 3
回答已采纳
1回答
如何快速可视化R?中的大型矩阵
我有时使用大数值矩阵(例如3000 x 3000)，快速可视化它们是一个非常有用的质量控制步骤。在我以前选择的语言Matlab中，这是非常简单和快速的。例如，显示1000x1000矩阵需要0.5秒：
rand_matrix = rand(1000,1000);
tic
imagesc(rand_matrix)
toc
>> Elapsed time is 0.463903 seconds.
我希望R中有相同的幂，但不幸的是，可视化矩阵在R中看起来非常慢。例如，使用image.plot()显示相同的随机矩阵需要超过10秒：
require(tictoc
浏览 0提问于2018-05-16得票数 1
回答已采纳
我想要实现一个天气层，它使用GeoServer在地图中显示从-40到40度的温度。我有一个SQL Server表，它看起来类似于：
create table temperature (id int identity primary key, geom geometry, temperature float)
Geom包含一个适用于给定温度测量的点。我在GeoServer中将表添加为图层，但我仍然无法将它呈现出来。现在，我得到的最好的是一个大的红色正方形:D，它似乎呈现为一个纯点，而不是光栅。
如何为光栅的良好显示创建正确的绑定等？
是在SLD级别还是在其他地方完成的？
浏览 1提问于2022-03-14得票数 0
回答已采纳
我尝试在raster包中使用distanceFromPoints函数，如下所示：
distanceFromPoints(object,xy,...)
其中，object是栅格，xy是x和y坐标的矩阵
现在，如果我的栅格有1000个像元，xy表示一个点，我就会得到1000个值，表示xy和每个栅格像元之间的距离。我的问题是当xy有多个坐标时，例如，10个点。函数描述表明xy可以是多个点，但是当我使用多个XY点运行这个函数时，我仍然只得到1000个值，而我期望XY中的每个坐标都有1000个值。这是怎么回事？
谢谢!
浏览 2提问于2016-12-28得票数 1
我已经编写了SAFT (合成孔径聚焦技术)的代码，并成功地绘制了2D图像。现在，我在3D中绘制它有一个问题，因为我不知道我应该使用什么类型的函数。我想创建3D图像，如下面的图像链接所示(右侧的图像具有实体体积)，但我对他们使用的功能一无所知。我所知道的就是左边的图片是拼接出来的。
如果有人能帮我解决这个问题，我将不胜感激。提前谢谢你！
这是我在3D数组中使用的示例代码。请注意，我只使用了强度的随机值，这里是V，因为实际值存储在excel文件中。我在使用等值面时遇到错误，因为我的数组不是使用meshgrid创建的。
r=linspace(50,500,101)';
o=zero
浏览 1提问于2014-05-03得票数 1
1回答
我正在尝试使用raster包中的predict函数来预测一堆栅格的负二项模型。我需要包括一个偏移量来规范化我的count变量。我曾尝试使用一种方法使其工作，但失败了，该方法将偏移项包括在模型中，如下所示：
condor.glm <- glm.nb(y_count ~ logsafefood + logpigharvest +
logintdist + houseden + pubforest + pubrange +
privforest + privrange + offset(log(off
浏览 0提问于2016-02-27得票数 3
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
我需要在gstat库中创建一个用于空间插值的遮罩网格。具体地说，我有不同的随机分布的采样点，我需要创建包围这些点的最小凸多边形。然后，我必须创建一个空间网格，它应该被计算的外壳裁剪，以便将插值限制在这个多边形的范围内。如果有人能给我解释一下详细的程序并提供一些例子，我将不胜感激。提前谢谢你。
浏览 0提问于2012-10-18得票数 2
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
回答已采纳
我想提取位于光栅图像中心的像素的纬度和经度信息。如果光栅图像不是精确的正方形(例如，它不是50x50)，我想得到大约中心像素的经度/经度。以下是一些可重现的数据： library(raster)
r = raster(ncol=25, nrow=25, xmn=-1000, xmx=1000, ymn=-100, ymx=900, vals = 1:(25*25))
crs(r) = crs(raster())
浏览 14提问于2021-04-19得票数 0
回答已采纳
1回答
使用空间多边形改变栅格值
、、、、
要在SpatialPoints下更改光栅值，只需使用[即可。
r <- raster(system.file("external/test.grd", package="raster"))
rTp <- rasterToPoints(r, spatial = T)
set.seed(666)
rTpS <- rTp[sample(1:length(rTp), 500),]
plot(r)
plot(rTpS, add = TRUE, pch = ".")
现在，空间点下面的栅格值会发生如下变化：
r1 <- r
r1[r
浏览 4提问于2016-07-27得票数 8
回答已采纳
1回答
我正在使用QGIS，并试图在栅格层中每个像素的中心生成均匀的空间点。我能够生成均匀分布的点，但无法使它们足够靠近或位于每个像素的中心。每个像素应该是1.85m，并且在UTM中的投影CRS上。我没有访问ARIGIS的权限。在QGIS中有没有办法做到这一点？甚至是在R中？我尝试生成一个网格以及生成点，但默认距离是0.0001，当我试图变小时，它会自动恢复。考虑到我在UTM，我也不确定如何解释这个默认值。
浏览 6提问于2018-01-31得票数 1
1回答
将栅格与圆相交
、、、、
我有以下问题。我有来自CH州(弗里堡)的数据，我想将3个圆与在该地区创建的网格单元相交。以下是生成附加地图和多边形的一些代码：
library(raster)
centroid.x <- c(566052, 576768, 560652)
centroid.y <- c(155501, 185501, 165325)`
centroids <- data.frame(centroid.x = centroid.x, centroid.y = centroid.y)
radius <- 1000
n = 1000
pts = seq(0, 2 * pi, leng
浏览 0提问于2017-06-27得票数 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
回答已采纳
我正在使用webots在一个模拟机器人模型上测试我的算法。我想在不平坦的地面上测试我的机器人，以证明其健壮性。我在webots中找到了一个名为unevenTerrain的地板对象，它可以基于珀林噪声随机生成不平坦的地形。还好，但我想我可能需要一个比这个物体更不平坦的地面。我如何在Webots中获得更不平坦的地面，或者有什么方法可以创建我自己的不平坦的地形？
浏览 130提问于2020-04-12得票数 1
我有一个大小为n_X * n_Y的离散网格，其中中心单元位于(0,0)。在这个网格中，在(x_target, y_target)上有一个相对于中心的随机点，它可以随着时间的推移改变它的位置。
我想根据每个单元格与随机点的接近程度，给每个单元格一个介于0和1之间的权重。最远的点应该是0，最近的点应该是1。中间单元应该基于它们的距离在两者之间进行插值。
系统：
除了知道随机点的位置外，你知道我如何估计每个单元格的权重吗?假设最远的点是0，最近的点是1。谢谢。
浏览 0提问于2021-09-27得票数 0
5回答
我正在尝试写一个简单的光线跟踪器作为一个爱好项目，现在一切都很好，除了我根本不能让柔和阴影工作。我对软阴影的想法是，光源被认为是有位置和半径的。为了在这个灯光上做阴影测试，我取主要光线击中场景中对象的点，并向光源投射n-数量的光线，其中每条新光线都有一个到每个轴的随机分量，其中随机分量在半径和-radius之间变化。
如果这样的光线击中场景中的一个对象，我会递增一个命中计数器(如果一条光线击中多个对象，它仍然只递增一个)。如果它在没有碰撞的情况下到达光源，我将主光线的交点到光源中心的距离添加到一个变量中。
当取了n个样本后，我计算了碰撞的光线的比率，并将光线的颜色乘以这个比率(因此，颜色为10
浏览 1提问于2009-02-05得票数 3
回答已采纳
1回答
我的以下问题基于@jbaum在这篇文章中提出的解决方案：
为了重现示例，我有一个到最近海岸线的距离的栅格数据集：
library(rasterVis); library(raster); library(maptools)
data(wrld_simpl)
# Create a raster template for rasterizing the polys.
r <- raster(xmn=-180, xmx=180, ymn=-90, ymx=90, res=1)
# Rasterize and set land pixels to NA
r2 <- rasterize(
浏览 14提问于2017-08-01得票数 6
回答已采纳
我有一个包含三列的文件。这三个都有不同的值。要将它绘制在具有第三列颜色渐变的平滑表面上，我应该做什么？前两列是伪随机分布的。最后一列也是如此。
数据文件如下所示：
8.4295190 0.3860565 0.3706621
-2.9886350 -0.1156874 -0.1314160
8.4375611 0.2617630 0.3710158
8.4092863 0.3195774 0.3697725
8.4237288 0.3930579 0.3704075
-1.1439280 -0.7286996 -0.091929
浏览 7提问于2014-08-16得票数 1
1回答
我正在做一个猜词游戏，我已经创建了一个文本文件，每行都是一个单词( 1000行)。 我不想把整个文本文件放到内存中。我只想读取文本文件中的任意一行，并将其放入字符串中。 我做了一个快速搜索，每个人都将整个文件读取到内存中，并从中获得随机行，但我想读取字符串形式的随机行。
浏览 24提问于2019-08-31得票数 1
回答已采纳
1回答
我们如何决定延长时间w.r.t。线程
当我运行1200秒的24000线程时，它工作得很好，但是当我减少到1000秒时，我得到了随机错误(503)，但是当我减少到500-800时，我得到了50%的错误等等。线程和爬升比有什么标准吗？
提前谢谢。
浏览 4提问于2014-08-07得票数 1
回答已采纳
在GIS extension中，我将一个以橙色表示的栅格导入到Netlogo中，如下图所示。我的目标是在用蓝色表示的栅格的每个块中随机放置一只海龟(即在9个块中)。蓝色光栅是橙色光栅的一部分。
下面是我将海龟随机放入块中的代码
let number 1
ask n-of number patches with [ (max-pxcor - ((x-increment + 1) * (max-pxcor / 3))) <= pxcor and pxcor <= (max-pxcor - (x-increment * (max-pxcor / 3))) and (max-pycor -
浏览 1提问于2014-10-09得票数 0
我已经从包含R中某一景观的高度值的光栅中创建了以下过渡对象：
wd <- "C:/Users/LG/Dropbox/Random Walk"
setwd(wd)
library(gdistance)
library(raster)
r <- raster("altitude.tif")
alt
浏览 0提问于2016-04-11得票数 0
1回答
我试着在银行账户之间建立随机交易。我创建了以下查询：
//Create transactions
CALL apoc.periodic.iterate("
match (a:BANK_ACCOUNT)
WITH apoc.coll.randomItem(collect(a)) as sender
return sender", "
MATCH (b:BANK_ACCOUNT)
WHERE NOT sender = b
WITH apoc.coll.randomItem(collect(b)) as receiver
MERGE (sender)-[r:HAS_TRANS
浏览 5提问于2020-10-21得票数 1
回答已采纳
1回答
我在栅格化shapefile以在0.5*0.5网格上生成点时遇到了问题。shapefile表示全球珊瑚礁的风险级别(低-0，中-100，高-1000，极高-1500)到综合威胁的分类。
我从另一个运行良好的示例中提取了代码，但是当我对我的数据尝试它时，我从plot函数中得不到任何东西。有关shapefile和我的代码的链接，请参见以下内容：
# Read shapefile into R
library(rgdal)
library(raster)
int.threat.2030 <- readOGR(dsn = "Global_Threats/Integrated_
浏览 7提问于2016-01-30得票数 3
回答已采纳
我正在优化我的粒子渲染器与GL_POINTS一起工作，现在我需要调整点的大小，使用顶点着色器中的gl_PointSize来缩放粒子，从顶点着色器的正确数量。这是我现在拥有的顶点着色器：
#version 330 core
layout (location = 0) in vec3 position;
layout (location = 1) in uint uv;
uniform mat4 projection;
uniform mat4 view;
void main(){
gl_PointSize = 10; // No difference with gl_PointSi
浏览 6提问于2021-09-05得票数 4
回答已采纳
3回答
我正在制作一个当按钮被按下时旋转的物体，我想让它旋转360度，然后继续旋转一个随机数，这样每次它都会落在不同的地方。这是我的xml文件(称为动画)中的内容，它完美地旋转了360。
<rotate
android:fromDegrees = "0"
android:toDegrees = "360"
android:pivotX = "50%"
android:pivotY = "50%"
android:startOffset = "0"
androi
浏览 0提问于2012-11-08得票数 1
回答已采纳
1回答
栅格中的越界颜色比例图
、、、、
我有一个光栅，我想用image()绘制它。到目前为止，我用自己的颜色比例col绘制了光栅。
## read the libraries
library(raster)
library(fields)
library(grDevices)
##random raster object
set.seed(1)
r <- raster(ncol=5, nrow=5)
r[] <- rnorm(n=ncell(r),mean=2)
col = colorRampPalette(c("darkred","red","lightskyblue
浏览 7提问于2018-08-16得票数 1
我在selenium中创建了一个随机的电子邮件方法，但我想知道如何在Selenium web驱动程序中为相同的方法创建一个永久方法，这样可以节省测试时间。 driver.findElement(By.id("email")).click();
Random randomEmail = new Random();
int randomInt = randomEmail.nextInt(1000);
driver.findElement(By.id("email")).sendKeys("username" + ra
浏览 17提问于2021-11-23得票数 0
回答已采纳
2回答
我有大量的数据要放入一个数据库，这是一个GPS点的列表。
我想遍历这个数据库并创建一个“热点”表，其中在一定大小的区域(无论是正方形区域还是圆形区域--我不需要精确地说)中有大量的数据库点。
有没有人能推荐一些对我有帮助的现有算法？
提前感谢！
r3mo
浏览 0提问于2011-01-26得票数 1
回答已采纳
我正在处理R (raster软件包)中的MODIS栅格数据，我想知道为什么当将两个栅格文件相乘时，栅格的大小会变得7000倍。在我的例子中，一个栅格是一个简单的MODIS栅格(介于0和255 (整数)之间的值)，另一个是一种模板(0和1(整数))。所以我不明白为什么栅格要这么大？当尝试用writeRaster保存栅格时，这是一个很大的问题，因为我的R会话总是崩溃。
我的代码如下所示：
a <- list.files('All_MODIS_files_of_one_year.tif')
for (day in 1:length(a)){
ra <- raster(a[
浏览 0提问于2018-06-16得票数 2
我正在尝试在R中创建物种分布模型。我已经在ArcMap中创建了栅格层，并将它们导入到R中。它们不能堆叠，除非范围完全相同，并且它们都具有相同的行数和列数。
然而，当我改变这些因素来成功地堆叠它们时，它们会丢失所有的值，并且我堆叠的数据帧只填充了NAs。
有谁知道我如何改变栅格图层的范围和分辨率，以便可以成功地堆叠它们--这样我就可以将环境信息附加到存在点。
干杯
浏览 2提问于2015-08-27得票数 3
我有一个巨大的光栅(> 20 GB)，并试图将其转换为一个txt文件。有没有办法在并行R计算中做到这一点？我无法将RAM加载到内存中，因为它太小了。
我的电脑有8个内核和16 Gb的RAM
最好的
浏览 0提问于2015-09-02得票数 0
嗨，我正在做我的第一个随机漫步程序。我只是想出了随机游走程序本身。它是一个简单的一维随机行走，有1000步。到目前为止，我的代码如下所示：
stepslim = 1000;
rnew(1) = 0
r=0;
%Now we will set up the for loop
for isteps = 2:1:stepslim;
if rand <= 0.5 %The if-else statement will tell the walker
step = -1; %which direction it steps i
浏览 4提问于2017-11-30得票数 1
回答已采纳
2回答
打印较大的栅格文件会生成白色边框
、、、、
我目前使用的是256x256像素的ASCII矩阵。我正确地将其导入到R中，并对其进行了光栅化，得到了我所期望的值(即，正确的x和y边界以及最小和最大"z“值)。然而，当我绘制它时，我得到了一个空白的光栅，就像矩阵中的每个值都是零一样。 我试着创建另一个5x5矩阵的文件，没有任何问题。我是不是遗漏了什么？ 文件和屏幕截图如下：my 256x256 raster https://gofile.io/d/JGApXI ascii矩阵链接
浏览 25提问于2021-05-07得票数 0
回答已采纳
1回答
我正在尝试使用R中的ggplot2一次显示多个图。
数据中有200多个城市，但我只创建了一个随机样本数据，如下所示：
#sample data
df <- data.frame(
city = c('paris', 'tokyo', 'seoul', 'shanghai', 'berlin', 'rome')
, week = c(41,42)
, type = c('A','B','C')
, count =
浏览 2提问于2016-11-02得票数 0
我有一列1000个随机数，我想计算一个窗口长度为10的滚动和，在Excel或R中，对于前10个观察，滚动和仍然应该用部分窗口长度来计算。不过，我没有在Excel或R中找到滚动和功能。
我目前在R中的尝试是使用循环。然而，循环工作正常，效率很低。请看下面我的尝试。
vector = as.vector(matrix(1,1000,nrow = 1))
for (i in 1:1000) {
if (i<11) {
rolling.sum[i]=sum(vector[1:i])
} else {
rolling.sum[i] = sum(vector[(i-10):
浏览 4提问于2015-08-10得票数 0
回答已采纳
1回答
从以下向量中：
depth = [0:0.1:20];
我如何创建一个新的向量，它的数据与'depth‘相同，但每0.5取一次值，结果为：
d2 = depth(1:5:end);
因为我的数据实际上不是以0.1的间隔间隔，而是以随机的间隔间隔，所以我不能应用上面显示的简单方法。实现这一目标的最佳方法是什么？
浏览 0提问于2012-11-06得票数 1
回答已采纳
1回答
我正在做Asteroids (Atari)在C。我为小行星设计了ADT，所有的ADT都与1颗小行星相关。
我需要在一个随机的位置创建一个小行星，但是在(0, y)或(x, 0)中是随机的
我有我的函数asteroid_create()
asteroid_t *asteroid_create(float x, float y, float radio){
asteroid_t *ast = malloc(sizeof(asteroid_t));
ast->x_pos = x;
ast->y_pos = y;
ast->radio = radi
浏览 0提问于2019-07-25得票数 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
回答已采纳
1回答
栅格转多边形-零值
、、、、
我正在尝试将单元大小为50x50的栅格转换为多边形，但没有成功。
光栅文件有一些值，其中大多数是零值。当我使用工具raster to polygon转换光栅时，结果是所有值都等于零的大多边形，而不是每个值为零的单元的多边形。
有人知道如何为每个零值创建一个独立的多边形吗？
我附上了一张我的结果的图片。
浏览 8提问于2018-07-30得票数 0
1回答
我在试着理解一个愚蠢的句子：
这个函数将进行一次实验。它将产生npoint随机有理数对(x，y)，形成npoint随机二维点。这些点的x和y坐标应保持在区间-r，r中。
因此，在这个函数中，npoint是整数，值为1000，可以是不同的值。函数得到数字a和b，度thera1和theta2，还有r，也就是r(θ) = ae^bθ。那么，谁能解释一下，这个函数到底要做什么呢？
PS。它并不是真正的bug，问题是我无法理解它可能是因为我的英语不是那么好，对不起，但它不是英文网站的代码:)
浏览 1提问于2013-12-12得票数 0
回答已采纳
1回答
我尝试用(a+(b-a)*R)产生离散的均匀随机数，其中R是随机样本，它是由线性同余法产生的。但我仍然对D统一数的创建有疑问。请给我离散均匀随机数的完美公式？
浏览 2提问于2014-04-07得票数 0
在bash中生成我使用的1和0随机数。
#!/bin/bash
for i in `seq 1 100000`; do
let rand=$RANDOM%2
echo $i $rand >> r-test.dat
done
同时使用，
集randm = awk -v min=1000 -v max=10000 'BEGIN{srand(); print int(min+rand()*(max-min+1))}'
产生随机数，例如在csh中产生1000到10000之间的随机数。现在，我的问题是如何修改这个命令，或者使用一种更有效的方法在csh中随机
浏览 1提问于2018-05-24得票数 1
回答已采纳
我想让writeRaster编写我在R中构建的RAT (栅格属性表)。
我运行的是R3.0.1，raster2.1-49和rgdal 0.8-10。
我的输入栅格如下所示：
r <-raster("F:/test.img")
class : RasterLayer
dimensions : 3, 3, 9 (nrow, ncol, ncell)
resolution : 30, 30 (x, y)
extent : 347325, 347415, 4301655, 4301745 (xmin, xmax, ymin, ymax)
coo
浏览 0提问于2014-05-24得票数 9
1回答
我已经通过以下代码生成了1000次模拟。
ar1 = 0.4857
ar2 = 0.0173
ma1 = -0.8054
r0 <- 0.002937432 #mean of margin.logrtn
e0 <- 0.002976873 #mean of fit.margin.logrtn$res
r.sims <- matrix(rep(0,1000*52),nrow=52, ncol=1000)
e.sims <- matrix( rnorm(52*1000, mean = 0.002976873, sd = 0.1056021), 52, 1000)
浏览 3提问于2013-10-15得票数 0
回答已采纳
我正在尝试使用matplotlib绘制球体周围磁场的流线，它确实工作得很好。然而，生成的图像不是对称的，但它应该是(我认为)对称的。
这是用于生成图像的代码。请原谅我的长度，但我认为这会比仅仅发布一个不能工作的代码片段更好。而且，它也不是很有pythonic风格；这是因为我将它从Matlab转换而来，这比我预想的要容易。
from __future__ import division
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Circle
def cart2sphe
浏览 3提问于2013-04-27得票数 13
回答已采纳
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
我的补丁包含诸如标高之类的属性：
set mnt gis:load-dataset "F:/StageM2/Modelisation/Modele/mnt.asc"
gis:apply-raster mnt alt
gis:set-transformation (list 567887.504252 573503.504252 6183200.86463 6187628.86463) (list min-pxcor max-pxcor min-pycor max-pycor)
gis:set-world-envelope gis:envelope-of mnt
海龟是从森
浏览 5提问于2014-07-15得票数 2
我需要创建3个这样的整数：
for i in range(1000):
x = random.randint(0, 4)
g = random.randint(0, 4)
z = random.randint(0, 4)
如果x，g，z不相等，我试过if x != r and z and g:，但没有工作，它只给了我3个随机整数，我的代码也是这样。
浏览 2提问于2021-05-06得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
