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
我想在geoTIFF中提取4个光栅子图像中心的坐标(并创建一个数据框)，并绘制最初在r2图像中拆分的2个子图像的边界，为此我尝试： library(raster)
library(rgeos)
#create GeoTIFF raster
r <- raster(ncol=100, nrow=100)
s <- stack(lapply(1:3, function(i) setValues(r, runif(ncell(r)))))
f1 <- file.path(tempdir(), "sl1.tif")
writeRaster(s,f1, form
浏览 23提问于2020-01-22得票数 0
回答已采纳
1回答
选择部分彩色地图绘制地形
、、、、
数据
.tif数据DEM是一个行政区划的高程值。我上传它
除了分区以外的区域我不想在情节上给他们看。
我的目标
用plt.pcolormesh绘制这个部门的催眠图(PS:我在阅读.tif时发现，plt.imshow()比pcolormesh__快得多。)我也不知道为什么)。
这里我展示了一个我从互联网上剪下来的例子。
我的尝试
### Using GDAL to read the .tif data
from osgeo import gdal
### Read the .tif
pathToRaster = r'./dem.tif'
raster = gda
浏览 1提问于2016-03-10得票数 3
回答已采纳
1回答
如何用不同程度的光栅制作光栅砖？
、、、、
我是R的新手，所以这个问题是非常基本的，但是我一直在与它斗争，无法找到一个有效的解决方案。我想从同一地区的一些陆地卫星图像中创建一个栅格砖。它们是以HDF-EOS格式下载的，我使用Modis重新投影工具将它们转换为.tif。
由此产生的光栅具有相同的投影，但在其范围，分辨率和来源不同。
我尝试了几种方法，概述如下：
手动定义子集范围并删除所有栅格。然后试着用子集栅格做块
重新对栅格进行采样，给它们提供相同数量的列和行。理想情况下，这将确保栅格细胞是对齐的，可以放在栅格砖中。这个选项创建了一个砖块，光栅没有价值，它们是空的。
我想知道我应该遵循什么概念来修正这个范围。创建一个空栅
浏览 2提问于2013-12-22得票数 9
回答已采纳
我试图使用来自zonal_stats包的函数rasterstats从.shp文件中的每个形状的.tif文件中获取光栅统计信息。我设法在QGIS中做到了这一点，没有任何问题，但我必须对200多个文件进行同样的操作，这将需要很长时间，所以我正在尝试Python方法。文件和复制代码都在我的中。
我的剧本是：
import rasterio
import geopandas as gpd
import numpy as np
from rasterio.plot import show
from rasterstats import zonal_stats
from rasterio.transfor
浏览 20提问于2022-09-13得票数 1
1回答
我成功地在R中编写了一个for循环，这是好的，我很高兴它能工作。但是我也想了解我所做的事情，因为我在以后的分析中也要使用循环。
我使用Raster Data (DEM)。我将它们作为光栅加载到环境中，然后在循环中使用getValues函数，因为我想做一些计算。情况如下：
list <- dir(pattern=".tif", full.names=T)
tif.files <- list()
tif.files.values <- tif.files
for (i in 1: length(list)){
tif.files[[i]] <- r
浏览 3提问于2017-11-10得票数 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
回答已采纳
我对R非常陌生，我想知道是否有一种更快的方法来在单个.pdf中保存多个栅格，而不是像下面这样手动输入每个栅格文件名：
pdf(file = "file_name.pdf", width = 11, height = 8.5)
plot(raster(file.path(dir_path_data,"nutrient_pollution_2003_impact.tif")))
plot(raster(file.path(dir_path_data,"nutrient_pollution_2004_impact.tif")))
plot(raster
浏览 3提问于2020-05-30得票数 0
回答已采纳
3回答
我是一个R新手，尤其是在空间数据方面。我正在尝试找到一种方法来高效地将多个(~600)个单波段光栅(.tif)文件导入到R中，所有这些文件都存储在同一个文件夹中。不确定这是否重要，但请注意，当在我的Mac和Windows并行虚拟机上的文件夹中查看时，每个.tif = .TIF；.tfw；.TIF.aux.xml；.TIF.ovr；.TIF.xml有以下五(5)种文件格式。无论如何，以下代码(以及我尝试过的其他类似变体)似乎不起作用：
library(sp)
library(rgdal)
library(raster)
#path to where all .tif files are loc
浏览 8提问于2018-10-11得票数 5
2回答
我已经将一个栅格文件导入到R中-它是分类的，有两个值，0和1。我需要将它存储为一个因子。
我尝试过以下几种方法：
#as.factor when importing
mydata.factor <- as.factor(raster("mydata.tif"))
#or import first then try converting
mydata <- raster(mydata.tif)
mydata.factor <- as.factor(mydata)
#or call just values
mydata <- raster(mydat
浏览 2提问于2014-08-30得票数 1
2回答
下载WorldPop人口计数光栅
、、、、
我有兴趣下载联合国调整后的刚果民主共和国2020年第二次决议(见WorldPop (见 )为GeoTIFF (.tif)文件)。我试图为多个国家做这件事--我在探索包。
devtools::install_github("wpgp/wpgpDownloadR")
library(wpgpDownloadR)
wpgpListCountryDatasets(ISO3="COD")
我能否得到一些帮助，直接下载特定国家的WorldPop人口计数.tif文件(例如，刚果民主共和国、尼日利亚、比利时、法国等)？
我有一个R代码来将这个GeoTIFF光栅(.tif文件)转
浏览 4提问于2020-11-24得票数 1
回答已采纳
我得到了一个tif，我试图用指定的投影/基准/等读入。这些是从ArcMap导出的tif，带有包含投影信息的.tif.xml和.tfw文件。在R中有没有一种方法可以通过.tif引入分配的coord.ref
在TIF中读取
r<-'example.tif'
r <- raster(r)
r
输出
class : RasterLayer
dimensions : 199, 695, 138305 (nrow, ncol, ncell)
resolution : 50000, 50000 (x, y)
extent : -17367529,
浏览 0提问于2017-06-27得票数 1
2回答
我有一个很大的tif文件，内存放不下。我想知道是否可以根据范围只读入此tif文件的一部分。
我正在使用R，到目前为止还没有找到解决方案。我知道它可以使用C#来完成，所以我想在R中应该有一种方法。
谢谢!
浏览 0提问于2018-03-29得票数 2
这个问题很模糊，但这就是我的意思。我有一个110 20的光栅文件，这是20个其他文件的组合，导致整个欧洲的.tif文件。我想使用它的二进制类型来创建3D地形，比如在this blog中，但问题是它的大小。我在gdal中尝试了一些技巧，例如： gdal_translate
-co COMPRESS=JPEG
-co TILED=YES
input.tif output.tif 但显然，JPEG压缩不适用于Float32文件。我更喜欢保存Float32格式的负值，但16位整数也可以(它排除了所有字节压缩)。到目前为止，我发现对Float32文件最好的压缩是使用PREDICTOR
浏览 81提问于2020-04-28得票数 0
回答已采纳
当我读入非常大的tiff文件时，如下所示
r = stack( 'tiffile.tif')
那么我想要一小部分。
r_part = r[1:100, 1:100, ]
我希望发现r_part是一个3D光栅堆栈。然而，它被证明是一个具有维度、宽*高和通道数的矩阵。
如何以这种方式提取3D光栅堆栈？
浏览 1提问于2018-05-29得票数 1
1回答
有没有办法将R中的*.las点云导出为正交马赛克？我使用包lidR加载了包含这些点的las文件。我想导出一个在RGB中从上方显示点云的tif，类似于正射影像的样子。这些数据是使用地面激光扫描仪收集的。 point cloud this is what I want
浏览 16提问于2021-03-28得票数 0
回答已采纳
我正在尝试使用Python.Can在Erdas中对Band1、Band2和Band3的堆栈图像进行分层。
from imagine import modeler
imagePath = "'C:\\Desktop\\Erdas_Script\\New folder'"
outputPathStack = "'C:\\Desktop\\Erdas_Script\\New folder'"
m = modeler.Model()
ri1 = m.RasterInput("BAND2.tif")
ri2 = m.Rast
浏览 0提问于2015-02-09得票数 1
问题是在将30×30m空间分辨率的栅格导入Postgis数据库时，如何选择栅格大小，以优化与复杂多边形的空间连接。
，即，在raster2pgsql中，如何选择"-t“开关来优化空间连接查询？
栅格数据表示由SRTM数据导出的30m分辨率的斜率。多边形代表布拉斯利亚人口普查的计数区(setor censitário)。一个枚举区对应于覆盖300个家庭的地区(美国人口普查轨道的1/4 )，在城市地区，可以是1至8个城市街区。
我试图获取每个枚举区的平均地形斜率，使用以下查询：
-- importing
shp2pgsql -c -s 4674:4326 -I -W LATIN1 enum
浏览 5提问于2014-06-06得票数 4
回答已采纳
1回答
我正在尝试使用R中的栅格库中的堆栈函数创建一个RasterStack对象。
library(raster)
但我对函数中使用的参数有问题。让我展示一下我在做什么：
###set working directory
setwd("myworkingdirectory")
###Upload variables
v1 <- raster("variable1.tif")
v2 <- raster("variable2.tif")
v3 <- raster("variable3.tif")
v4 <- rast
浏览 0提问于2013-06-27得票数 2
回答已采纳
1回答
我添加了一个新的变量来编码我的形状文件的单元格，我想将该文件保存为.tif，并拥有我的新变量的所有信息。我将它保存为一个新的.shp文件，然后使用shp2raster函数，但我得到了错误。 beijing10<-sf::read_sf("Landuse/Beijing2010.shp")
beijing10<-mutate(beijing10,land=ifelse(CODE%in%c(41,42),"Cultivated Land","Unused Land"))
st_write(beijing10,"beijing20
浏览 58提问于2020-08-14得票数 0
回答已采纳
1回答
我使用以下代码将一个栅格读入我的R会话：
raster <- stack("raster.tif")
现在我想创建一个简单的特征(sf)对象来表示该栅格的轮廓。我不能使用边界框，因为光栅是多部分的，所以边界框会比光栅大得多。因此，轮廓线也需要是多部分特征(sf多多边形)。
如果能帮上忙我将不胜感激。谢谢!
标记
浏览 2提问于2020-05-29得票数 0
我正在处理28个关于海洋初级生产力的HDF4文件(年度.tar文件可以在这里找到：)我的目标是做一些计算(我需要计算每个区域的浓度并获得几年的平均值，即在空间上合并所有文件)，然后将它们转换为我可以在ArcGIS中使用的地理参考文件(最好是shapefile或geotiff)。
我已经尝试了几种方法来转换为ASCII或光栅文件，然后使用gdal_translate和get_subdatasets等gdalUtils工具添加投影。但是，由于HDF4文件不是以标准命名的(与MODIS文件不同)，后者不起作用，我也无法访问子集。
下面是我用来转换为光栅的代码：
library(raster)
lib
浏览 15提问于2018-02-04得票数 1
回答已采纳
我是raster2pgsql的新手，所以请容忍我。我试图用以下代码将60 my的.tif (从高分辨率清算层项目)加载到启用postgis的数据库中：
raster2pgsql -s 5235 -C -F路径到.tif public.hrsl_lka \ psql -h localhost -U postgres -p 5432 -d项目
但是，我得到以下错误：
错误: insert_records:无法为插入语句错误分配内存: process_rasters:无法将光栅块转换为插入或复制语句错误:无法处理光栅
但是，将大约3mb的较小的.tifs加载到同一个数据库中，但是从其他来源加载
浏览 0提问于2020-12-10得票数 1
回答已采纳
2回答
我在R中使用非常大的栅格，我有一个代码，我以前已经能够用相同的大栅格运行过，但是当我试图再次运行我的代码以获得可重复性时，我得到了相同的错误“在栅格IO期间失败”。
例如，在代码的这一部分中，我试图将一个大的栅格(~ 5 GB)分成较小的块，这样我就可以在R中加快处理时间。
large_raster <- stack("D:/large_raster.tif")
tileno <- 2 # split the large raster into 2 tiles
cl <- makeCluster(3)
registerDoSNOW(cl)
foreach(i
浏览 67提问于2020-05-29得票数 0
回答已采纳
1回答
玩具数据
library(terra)
template_raster <- rast(xmin=0, xmax=10, ymin=0, ymax=10, crs="+proj=longlat +datum=WGS84 +no_defs +type=crs", resolution = 0.1)
values(template_raster) <- c(rep(NA,4000),rep(10,2000),rep(NA,4000))
polygon_to_rasterize <- as.polygons(ext(3, 6, 3, 6), crs="+
浏览 5提问于2022-09-30得票数 2
回答已采纳
2回答
我的最终目标是使用R包spdep中的函数joincounts.multi()计算连接计数。
我有一个非常大的栅格文件(n行= 19663，n行= 34073)，单元格大小为30m，它是一个二进制栅格，有两个值，0和1(以及许多NoData单元格)。它是整数格式。
首先，我需要生成一个nb对象--即邻居列表。对于这个包，有多种方法可以做到这一点。对于多边形数据，可以使用poly2nb参数。对于网格单元格，可以使用cell2nb参数。显然，后者不适用于非常大的栅格-它会使我的电脑崩溃，谷歌搜索表明其他人也有同样的问题。
编辑:我已经决定在更小的网格上做这个功能。我有所有的网格在一个文件夹中，并希望
浏览 7提问于2014-08-26得票数 1
1回答
TIF文件在R中的问题解读
、、、、
我试着用光栅读取R中的Tif文件。基本上，我需要将来自tif文件(坐标)的数据匹配到我在另一个shapefile中定义的城市的特定边界。
我的TIF文件大约是7MB。
问题是，我知道如何将这些tif文件转换为可用格式的唯一方法是首先将我的RasterLayer转换为数据格式。然而，当我尝试这样做时，由于某种原因，dataframe最终有大约3亿行，而光栅层最初只有16,000行。这导致了一切的停滞。我想知道是否有一种更简单的方法将我的TIF文件转换成可读格式为R，然后我可以在我的Over函数中使用它。我的密码在下面。任何帮助都将不胜感激。
city_lights <- raster(ti
浏览 7提问于2020-03-06得票数 1
回答已采纳
我使用raster读取栅格( GeoTIFF)，然后将其转换为表，这样我就可以进行一些分析，并对一些值进行重新分类。如下所示： r_ndvi <- raster(paste0(path, "ndvi.tif")) #read in as raster
rdf_ndvi <- as.data.frame(r_ndvi) #turn tif into df of DN 完成这些计算后，我希望将表转换回光栅格式，以便可以在QGIS中将其作为图像打开。我如何才能做到这一点呢？
浏览 127提问于2021-09-27得票数 0
回答已采纳
4回答
我有一个c#程序，它打开一个.tif映像，然后提供保存它的选项。然而，在保存图像时，质量总是会下降。
(编辑:我在保存图像时传递了一些参数，这样质量达到100 %，没有压缩，但是实际唯一颜色的数量从254到16，尽管图像属性显示了8bpp)。
(EDIT2:问题中的图像是一幅灰度图像，每像素8位-- 256种颜色/灰度--这并不是发生在每像素24位的彩色图像中，而我测试的所有颜色都保留在那里。我开始认为图像类可能只支持16度灰色)
我怎么才能避免这种情况？
下面是打开图像的代码：
public Image imageImport()
{
Stream myStream = null;
浏览 8提问于2012-11-02得票数 4
回答已采纳
我需要将形状文件转换成栅格，并且完全不知道从哪里开始。
如果有人能帮我，我会很感激的！
更新:我发现了“readOGR”-函数，但是每次我使用它时都会收到以下消息：
Warning messages:
1: In readOGR(dsn = "C:/Users/Giaco/Dropbox/Random Walk/Waterbodies", :
Dropping null geometries: 308, 309
2: In readOGR(dsn = "C:/Users/Giaco/Dropbox/Random Walk/Waterbodies", :
浏览 3提问于2016-04-15得票数 0
2回答
在R中，我遇到了一个不寻常的问题，我们试图用另一层的值填充多云像素。我可以让它和堆栈一起工作，如下所示-
library(raster)
r1 <- raster(ncols=36, nrows=18)
r1[] <- 1:ncell(r1)
r1b <- r1a <- r1
r1_stack <- stack(r1, r1a, r1b)
r2 <- setValues(r1, runif(ncell(r1)))
r2b <- r2a <- r2
r_stack <- stack(r2, r2a, r2b)
r_stack[r_stac
浏览 3提问于2017-05-31得票数 3
回答已采纳
1回答
我想导入一个.tif格式的航空图像，然后导出相同的颜色格式。我想这是一个简单的任务，但我找不到办法去做。原始文件在和Google中查看时是彩色的，但是输出是黑白的。
最终，我想要裁剪/合并多个图像，从大约6个瓷砖创建一个单一的航空图像地图。在这6块瓷砖中，我只对大约30%的面积感兴趣(水面上的碎片，加上大约20米的缓冲区)，所以我的想法是通过将我想要的部分保持为一个文件来减少图像大小，而不是进口所有6个瓷砖。我可以打开谷歌地球的一个.tif作为彩色航空图像，这是我想要的自定义边界地图/图像。
目前，我正在努力导入，然后以同样的谷歌地球可读的彩色图像格式导出单一的.tif。我正在使用raster
浏览 5提问于2021-08-27得票数 0
回答已采纳
2回答
是否有一种简单的方法可以将栅格作为矩阵直接加载到R中，而不是加载栅格，然后使用as.matrix()将其转换为矩阵。
myras <- raster("file.tif")
mymat <- as.matrix(myras)
浏览 2提问于2018-01-18得票数 1
在写入然后读取具有某些文件格式(ESRI SpatRaster标签格式、netCDF和ENVI，但不包括GeoTIFF或Erdas HFA)的.hdr时，光栅原点、分辨率和坐标会稍有变化。足以让这些参数的相等性测试返回FALSE值。下面是一个示例，先使用GeoTiff，然后使用flt。 Nrow=45; Ncol=108
r1 <- rast(nrow=Nrow,ncol=Ncol,vals=1:(Nrow*Ncol),ext=c(-180,180,-60,90))
origin(r1)
# [1] 0 0
writeRaster(r1,"test.tif",over
浏览 8提问于2021-01-21得票数 0
1回答
回答了类似的问题；但是，这个问题有点不同，我无法应用该解决方案。我已经安装了使用网站数据格式的maxent。问题是我无法预测合适的模型。m1在路径D:/maxent中的输出看起来很好。我怀疑这两个错误(下面)与rJava有关，但我不知道解决方案。请参阅我的代码如下：
> m1 <- maxent(x = d, p = id, path = "D:/maxent",
args = c("-P", "noautofeature", "nolinear", "noquadratic"
浏览 5提问于2020-02-03得票数 0
回答已采纳
我在卡雷特包里做了一个随机森林模型。我叫它rf_final。我把我的光栅变量堆叠起来，我把它叫做CqStack。我希望使用clusterR()函数来执行预测。我使用了以下代码
cores<- detectCores()
beginCluster(n=cores*0.75)
pred.map<- clusterR(CqStack,predict,args = list(model= rf_final), progress="text")
endCluster()`
我一直收到以下错误。
in clusterR(CqStack, predict, args = list
浏览 11提问于2022-03-27得票数 0
1回答
我需要将一个.tiff文件从EPSG:4326 -wgs84-地理重新投影到EPSG: 32638 (UTM Zone 38N)。当我运行下面的代码时，输出文件是一个重新投影的.tiff，大小是原始文件的一半。
# Reprojection
from osgeo import gdal
filename = r"/home/DATA/Test.tif"
input_raster = gdal.Open(filename)
output_raster = r"/home/DATA/Test_UTM38N.tif"
gdal.Warp(output_rast
浏览 5提问于2020-06-17得票数 0
1回答
我正在处理一组13个.tif栅格文件，每个116.7 MB，包含关于西非红树林分布的数据。每个文件保存该分发版本一年(2000-2012)。栅格加载到R中没有任何问题，并且绘图也相对容易，使用base plot()大约需要20秒，使用ggplot()大约需要30秒。 当我尝试对栅格进行任何类型的处理或分析时，我都会遇到问题。我试着做简单的栅格计算，从2000栅格中减去2000红树林分布栅格，以显示森林砍伐热点，但一旦我这样做了，我计算机上的内存就开始迅速消失。 我有48 GB的空闲驱动器空间，但当我开始运行光栅数学时，我开始每隔几秒钟就会丢失1 GB的存储空间。这种情况会一直持续下去，直到我
浏览 15提问于2020-12-06得票数 0
1回答
我有一个关于Netlogo上GIS可能性的一般性问题。
我在QGIS上有一个全国土壤湿度的TIF文件。我正在创建一个Netlogo模型，在这个模型中，当输入(用户给出的地理坐标)是特定值时，我需要告诉补丁‘拥有’某个水分含量的量化值。然后，模型将引用QGIS文件，并根据给定的坐标输入特定值。
我不想将整个QGIS地图覆盖到我的Netlogo模型中，这是我到目前为止所看到的。有什么办法可以解决这个问题吗？提前感谢！
浏览 10提问于2016-05-27得票数 0
2回答
只读R中的一种作物或光栅的范围
、、、、
我正在处理每日tif文件中的大量数据(数千个每日文件)。我正在分析形状文件区域中光栅的平均值，在可能的数千个形状层上重复。
我当前的.tif文件是针对整个国家的，而实际上，我只需要为每个shapefile层提供每个tif文件中的一小部分。每个形状层有不同的日期集可从中提取，即数据帧如下：
df <- data.frame(shplayer=c("shp_layerbuffer1","shp_layerbuffer2", "shp_layerbuffer3"), start=c("2000_02_18", "2004_
浏览 0提问于2018-06-13得票数 5
回答已采纳
我正在尝试使用raster2pgsql-命令将GTiff-file插入到特定的PostGIS表中。到目前为止，我已经成功地将GTiff-file插入到我所连接的PostGIS数据库中。但这将使用GTiff-file的文件名创建一个新表。我也可以在之后将栅格数据移动到目标表中，但我想有一种更有效的方法。 下面是一个示例： import psycopg2
import os
tif_path = 'test.tif'
conn = psycopg2.connect(
host = 'localhost',
port = 5432,
use
浏览 63提问于2020-04-27得票数 0
回答已采纳
1回答
我使用for循环来迭代一堆栅格(n=533)，并基于将某些值(小于353.3)转换为'NA‘来对它们进行reclassify。我现在希望有一种有效的方法来搜索重新分类的栅格列表，并删除所有'NA‘值的栅格(参见下面的输出示例)。如何做到这一点？
wfrastlist <- list.files(path = "/path/to/rasters/",
pattern='*.TIF$', all.files=TRUE, full.names=FALSE)
#generate a reclassification matrix
#in th
浏览 0提问于2018-10-17得票数 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回答
我有一个.tif格式的栅格，我试图用平面图在R中可视化。以下是我迄今所做的工作：
library(raster)
library(lattice)
data_tif <- raster("fthrt14_21.tif", RAT = TRUE)
data_tif
rat <- read.dbf("fthrt14_21.tif.vat.dbf")
rat
data_tif <- ratify(data_tif)
colnames(rat)[1] <- "ID"
levels(data_tif) <- rat
l
浏览 4提问于2019-11-20得票数 1
回答已采纳
我想在R中的一个并行循环中读取单独的tif文件。我加载了栅格软件包，并且brick在循环外工作，但在循环中尝试使用brick失败了： library(foreach)
library(doParallel)
library(raster)
cores=detectCores()
cl <- makeCluster(cores[1])
registerDoParallel(cl)
tiffile="../ESACCI-LC-L4-LC10-Map-20m-P1Y-2016-v1.0.tif"
print("step1")
bdata=brick(ti
浏览 22提问于2020-01-10得票数 1
回答已采纳
我需要创建x个文件(一个集合)，但我必须首先检查这些文件是否存在同名的文件。
例如，tiftest1.tif，tiftest2.tif，...存在，并且我必须将tiftest再次写入同一目录。我想将_x附加到文件名的末尾，其中x是一个自动递增的数字，每次我想要创建集合时。所以我可以有tiftest1_1.tif，tiftest2_1.tif，tiftest1_2.tif，tiftest2_2.tif，tiftest1_3.tif，tiftest2_3.tif等等。
这是我到目前为止所知道的：
...
DirectoryInfo root = new DirectoryInfo(fileWatc
浏览 0提问于2013-07-13得票数 0
1回答
我想生成与河流集水区的栅格数据集的地形/坡度一致的矢量箭头。
我已经在ArcGIS中创建了一个点的渔网网格，我想为设置长度的每个点创建一个箭头，这些点将遵循坡度的形状，即遵循阻力最小的路径，线将在3 x 3网格中逐渐跟随较小的数字。
我想我可以用矢量图生成矢量箭头。是否有可能实现与光栅一致的线条？
更新:我从一个点的网格中生成了大约200,000条线。我将使用R将这些转换为光栅，并将其设置为与我的坡度光栅相同的分辨率。
关于如何在坡度上分层栅格线，以便我可以让线遵循坡度的最低值，有什么想法吗？
浏览 0提问于2017-09-13得票数 0
convert 0101.jp2 -threshold 50% -type bilevel -monochrome -compress LZW ../0101.tiff
当我使用上面的命令将彩色扫描的文本页转换为黑白图像(必须是每像素一位)时，产生的图像看起来是锯齿状的。我想让它的分辨率更高，看上去更流畅。如何使用转换来做到这一点？
请注意，这样可以自动将tif图像转换为jpg格式，因此下面显示的输出图像与输出图像不一样。要在tif中获得真正的输出映像，您需要运行转换命令。
浏览 5提问于2022-11-22得票数 0
1回答
我正在服务器中使用R运行一个物种分布模型。我正在保存我的所有R环境，但是当我试图可视化和绘制保存的对象时，我得到了一个错误"Error in file(fn，"rb")：cannot open the connection In file(fn，"rb")“：警告消息: In file (fn，”rb“)：cannot open file '/localscratch/anandam.9761522.0/RtmpnCnH0y/raster/r_tmp_2020-07-13_195526_260024_04625.gri'：No to f
浏览 13提问于2020-07-23得票数 0
1回答
我有一个要添加点的栅格对象。
有没有一种方法可以给出它们的确切位置，而不依赖于图的大小？
可重现代码：
r<- raster(nrows=10, ncols=10) r <- setValues(r, 1:ncell(r)) plot(r) points(x=-50,y=20) points(x=-50,y=-90)
我不被允许发布图片。但当我将设备变小时，较低的点会移出彩色区域，或者当我将设备变大时，较低的点会移到彩色区域。
干杯
浏览 1提问于2014-11-26得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
腾讯会议活动推荐
腾讯云开发者
