根据标准差拉伸波段创建新栅格(.tif)，适用于数据堆栈，但不能写入新文件Python
根据标准差拉伸波段创建新栅格（.tif），适用于数据堆栈，但不能写入新文件。
根据标准差拉伸波段是一种数据处理技术，用于增强图像的对比度和可视化效果。标准差是一种衡量数据分布广度的统计量，而波段是指图像中的不同光谱通道。
在Python中，可以使用开源库GDAL（Geospatial Data Abstraction Library）进行栅格数据处理。下面是一个示例代码，展示了如何使用GDAL实现根据标准差拉伸波段并创建新的栅格数据：
import gdal
import numpy as np
# 打开原始栅格数据
dataset = gdal.Open('input.tif')
# 读取波段数据
band = dataset.GetRasterBand(1)
band_data = band.ReadAsArray()
# 计算标准差
std_dev = np.std(band_data)
# 计算拉伸范围
min_value = np.mean(band_data) - 2 * std_dev
max_value = np.mean(band_data) + 2 * std_dev
# 根据标准差拉伸波段
stretched_band_data = (band_data - min_value) / (max_value - min_value) * 255
# 创建新的栅格数据
driver = gdal.GetDriverByName('GTiff')
output_dataset = driver.Create('output.tif', dataset.RasterXSize, dataset.RasterYSize, 1, gdal.GDT_Byte)
output_dataset.SetProjection(dataset.GetProjection())
output_dataset.SetGeoTransform(dataset.GetGeoTransform())
# 写入拉伸后的波段数据
output_band = output_dataset.GetRasterBand(1)
output_band.WriteArray(stretched_band_data)
# 释放资源
dataset = None
output_dataset = None在以上示例代码中，input.tif表示原始栅格数据文件名，output.tif表示输出的新栅格数据文件名。通过计算标准差和拉伸范围，将波段数据进行线性拉伸，并将结果写入新的栅格数据文件。
需要注意的是，以上示例只处理了单波段数据。如果有多个波段，可以循环处理或者进行其他处理方式。
对于数据堆栈，可以根据具体需求选择适当的栅格数据处理方法，并使用GDAL或其他相关库进行实现。
相关·内容
1回答
我有一个函数，将多波段图像读入为光栅砖块对象，迭代各波段进行各种计算，然后将光栅写成新的.tif。所有这些都很好，但是新图像文件的文件大小大约要大4倍(我假设是因为原始图像有4个波段)。我想知道writeRaster()函数中是否有我不知道的参数，或者是否有其他方法可以确保输出图像与输入文件大小基本相同。
原始文件大小为134 MB；输出范围从471到530 MB左右，取决于格式。
简化代码：
library(rgdal)
library(raster)
path = "/Volumes/ENVI Standard Files/"
img = "qb_tile.img
浏览 11提问于2017-02-04得票数 7
如果标题不清楚，我很抱歉，我是python的新手，我的词汇量有限。 我想要做的是对.tif栅格中的每个波段应用标准差拉伸，然后通过使用GDAL (Python)堆叠这些波段来创建一个新的栅格(.tif)。 我可以用不同的波段组合创建新的假彩色光栅并保存它们，我可以使用dstack (第一块代码)在python中创建我想要的图像，但我无法将该图像保存为经过地理校正的.tif文件。 因此，要使用dstack创建拉伸图像，我的代码如下所示： import os
import numpy as np
import matplotlib.pyplot as plt
import math
from o
浏览 30提问于2020-01-20得票数 0
我想要crop一个多波段光栅(4波段)由空间多边形(在SpatialPolygonsDataFrame)。当我在QGIS中显示原始栅格和裁剪后的栅格时，我发现裁剪后的栅格与原始栅格的颜色不同。下面是我的代码：
library(raster)
mosaic_shp <- shapefile("mo_clipper.shp")
mosaic <- brick('orthomosaic.tif')
mosaic_sub <- crop(mosaic, extent(mosaic_shp))
writeRaster(mosaic_sub, &#
浏览 0提问于2018-11-20得票数 0
1回答
我希望使用R提取地理参考光栅图像的每个像素的颜色值。我需要这些像素的颜色稍后才能像前面的问题(参见)中提到的那样绘制一个tif (或geotif)。由于带有数据波段的光栅图像不能隐式地使用带值为它们分配颜色(无法表示模式填充)，所以我需要那些像素颜色。我已经知道如何访问在向量中列出了所有可能的256个颜色的colortable。但是，它们不是以此形式映射的。下面是我用来加载光栅图像和提取未映射的colortable的代码
raster1 <- raster(paste(
workingDir, "/HUEK200_Durchlaessigkeit001_proj00
浏览 0提问于2013-10-06得票数 0
回答已采纳
我有两个包含大约15000个.tif文件的文件夹。第一个文件夹中的每个文件都是一个带有5个波段的栅格，命名为AA_"number“，意思是它看起来像
AA_1.tif，AA_2.tif，…，AA_15000.tif。
第二个文件夹中的每个文件都是一个栅格，带有两个名为BB_"number“的带，看起来如下
BB_1.tif，BB_2.tif，…，BB_15000.tif。
我的目标是在文件夹AA的第一个文件中添加第一个波段1-3，从BB文件夹中的第一个文件中添加第一个波段，创建一个4波段光栅，并制作15000个4波段光栅。在QGIS中做了一些研究和测试之后，我相信GDAL的工具
浏览 15提问于2022-09-06得票数 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
我需要知道如何创建虚拟栅格。我迭代了一个文件夹，其中包含一些二进制栅格(1，0)。我使用numpy.concatenate将这些栅格添加到一个numpy数组中。然后，我想创建一个虚拟栅格，使用连接的栅格数作为此栅格将拥有的波段数。不过，我收到了以下消息：
---------------------------------------------------------------------------
IndexError Traceback (most recent call last)
<ipython-input-93
浏览 43提问于2019-07-12得票数 0
回答已采纳
1回答
读取多个栅格tif作为数据
、、、、
我想读取几个光栅tif文件，并获得一个数据作为输出。每个光栅有两个波段。
我想得到这样的输出：
我尝试读取所有的tif文件，但我不知道如何将它们堆栈为dataframe：
import numpy as np
import glob
import rasterio as rio
final = []
tif_files = glob.glob(os.path.join(Repo_img, r'*.tif'))
for f in tif_files:
im = rio.open(f).read(1)
imarray = np.array(im)
浏览 6提问于2022-04-08得票数 1
回答已采纳
我读取了我的高光谱(.raw)文件，并将三个波段组合为"gai_out_r“，然后输出如下：
writeRaster(gai_out_r,filepath,format="GTiff")
最后我得到了gai_out_r.tif
但是，为什么Win10不能将这个小的tif显示为我从envi输出的图片--将图像保存为--tif
Win10显示的两个tiffs如下：
浏览 46提问于2019-01-06得票数 0
我有一个脚本，它从4波段的地理信息中计算出，并将NDVI写成一个地理信息。我正在研究1米空间分辨率的航空图像。输出NDVI测井值范围为-1 ~1，为32位浮点.如何应用线性拉伸，以便在中输出范围为0- 255。
file = 'D:\create_ndvi\3711101_se1.tif';
[I R] = geotiffread(file);
outputdir = 'D:\create_ndvi\'
NIR = im2single(I(:,:,4));
red = im2single(I(:,:,1));
ndvi = (NIR - red) ./ (N
浏览 4提问于2014-02-24得票数 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
1回答
R中光栅处理的内存问题
、、、、
我已经编写了一个脚本来自动处理光栅数据，但是它在完成一半之后就失败了。循环在目录unzipped中查找包含.tif文件的子目录，并将那些在/base_layers中有一些基本映射的子目录堆栈起来。它一直未能做到：
矩阵中的错误(unlist(Ini)，ncol = 2，byrow = TRUE)：'data‘必须是向量类型的，是'NULL’
和
In writeBin(v，x@file@con，size = x@file@dsize)：写入连接的问题
和
错误.local(.Object，.)：GDAL 3:可用磁盘空间为875151360字节，而至少需要24739624
浏览 5提问于2019-11-08得票数 0
回答已采纳
我正在处理R中的Landsat 8场景，以计算NDVI并运行土地覆盖分类算法。我对光栅包中的writeRaster函数有问题，特别是在磁盘上编写光栅堆栈时。
我开始加载landsat 8场景的12个波段，并将它们作为层叠加在光栅堆栈中。由于它们是作为16位图像交付的，所以所有层的值范围从0到65535不等。在磁盘上写入光栅堆栈之后，当我从R环境中的磁盘重新加载新创建的文件时，所有层的值范围都与原始值不同。我不知道为什么，我在网上也找不到任何解决办法。
这是代码：
library(raster)
# Load the individual bands of the Landsat scene.
浏览 3提问于2017-08-31得票数 1
回答已采纳
我无法在R-工作室中使用abline()绘制贴合线。
我使用analysis运行在Windows8.1机器上，并在LULC地图上下载了Landsat图像进行温度分析。利用Erdas Imagine中的第6波段，在ArcGIS中生成和导入LST图像。利用6个多边形，以网格格式生成6幅LST光栅图像。默认情况下，ArcGIS生成32位深度数据类型的栅格。我在这里的问题是，在将网格转换为ascii文件之前，是否有必要先将栅格转换为8位光栅网格中的第一个栅格，如果是，那么为什么？如果不是，那么在R中处理转换的最佳方法是什么?我自己使用ArcGIS将光栅转换为ascii，但我无法在其中显示任何数据，而不
浏览 3提问于2014-09-11得票数 0
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
回答已采纳
我有这段代码，我的目标是计算我的栅格的罪过的0.8次方。
import os
os.chdir('D:/NOA/Soil_Erosion/test_Project/Workspace/Input_Data_LS_Factor')
import rasterio
import math
data = rasterio.open('Slope_degrees_clipped.tif')
band = data.read(1) # array of float32 with size (3297,2537)
w = band.shape[0]
print(w)
浏览 8提问于2021-05-04得票数 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
我试图使用来自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
我试图从landsat8 geotiff文件中读取像素值和坐标，但在将文件作为数组读取后，它只对我的geotiff文件显示了零。我使用以下代码来读取该文件：
from osgeo import gdal
gtif = gdal.Open( "FileName.TIF" )
band = gtif.GetRasterBand(1)
print(band.ReadAsArray())
"FileName.TIF“是陆地卫星的波段图像之一。
我还想遍历整个图像并获得坐标，所以请建议是否有任何简单的方法来遍历图像。以下是我正在使用的图像的链接:：
https://drive.g
浏览 1提问于2018-04-06得票数 0
我有一个" tiff“文件，它是由陆地卫星7，波段3拍摄的。我想要得到它的地理信息。但是，当我尝试在图像上执行"geotiffread“函数时，它会说：
?? Error using ==> geotiffinfo>readinfo
Multiple images exist in the file and the HEIGHT or WIDTH values are not equal.
Error in ==> geotiffinfo at 215
info = readinfo(filename);
浏览 2提问于2012-04-03得票数 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
2回答
我正在尝试将成堆的图像写入TIFF文件。图像大小为256*256像素，每个堆栈包含1000幅图像。编写其中一个文件大约需要4分钟--所以我的代码很可能有问题。
我正在做的事情是：
void Tiff_WriterSplit(float data[], int PicNum, int WIDTH, int LENGTH, char PATH[]) {
int i;
int rows = LENGTH, columns = WIDTH;
TIFF* tif;
if (PicNum == 0)
tif = TIFFOpen(PATH, &#
浏览 5提问于2015-10-27得票数 3
回答已采纳
4回答
如何将点映射到扭曲的网格上
、、、、
假设你有一个点的集合，这些点在笛卡尔坐标系上有坐标。
您想要绘制另一个点，并且您知道它在同一笛卡尔坐标系中的坐标。
但是，您在其上绘制的图形与原始图形不同。想象一下，将原始平面打印在橡皮板上，然后以不对称的方式(没有重叠或任何复杂的东西)在某些地方拉伸，在另一些地方收缩。
()
您知道每组点的拉伸和未拉伸坐标，但不知道底层的拉伸函数。您知道新点的未拉伸坐标。
如何根据附近点的拉伸位置来估计在拉伸坐标中绘制新点的位置？它不需要精确，因为您不能从一组重新映射的点确定实际的拉伸函数，除非您有更多信息。
其他可能的关键字:扭曲扭曲网格网格平面坐标取消扭曲
浏览 0提问于2009-10-21得票数 11
回答已采纳
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
回答已采纳
当我读入非常大的tiff文件时，如下所示
r = stack( 'tiffile.tif')
那么我想要一小部分。
r_part = r[1:100, 1:100, ]
我希望发现r_part是一个3D光栅堆栈。然而，它被证明是一个具有维度、宽*高和通道数的矩阵。
如何以这种方式提取3D光栅堆栈？
浏览 1提问于2018-05-29得票数 1
我想堆叠和写一些Landsat波段/tiff文件，以ENVI格式的BIP交织。然而，结果总是以BSQ的形式出现，即使我将带宽顺序更改为BIP。
下面是我的代码：
library(raster)
library(rgdal)
library(gdalUtils)
inbands <-list.files(pattern= "*.tif")
stk<-stack(inbands[2], inbands[3], inbands[4])
writeRaster(stk, "BIP_test", format="ENVI", bandorde
浏览 7提问于2015-01-07得票数 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
回答已采纳
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
回答已采纳
我需要抓取光栅图像的每个像素值(.tif，单波段，像素值作为高程值)，并将其与另一个图像进行比较，看看像素值是否相同。我尝试过gdalcompare.py，但这只给出了文件名、文件类型、文件大小等一般差异。
我只能访问免费软件，如果能找到一种方法来做到这一点，那就太棒了，因为我在谷歌上的搜索一直是徒劳的
浏览 0提问于2018-08-06得票数 0
我想使用由FAA 提供的航空图表TIFF文件。
python模块rasterio似乎是我的正确工具，我已经看过几个教程和参考文档。他们都为我工作，但不是我想要使用的数据。我似乎无法使它在颜色上正确地显示。
主要的区别似乎是，我的数据是单波段的(如dataset.count所报告的)，而所有的教程都有三个或更多的波段。
例如，下面的代码应该显示图表：
import rasterio
from rasterio.plot import show
dataset = rasterio.open ('Green_Bay_101/Green Bay SEC 101.tif')
show(
浏览 0提问于2020-10-28得票数 0
1回答
下面是一个名为的光栅文件。
使用R，我可以得到这样的颜色表：
my.raster <- raster( "/path/to/singleband-raster.tif" )
colortable( my.raster )
它将给出如下所示的RGB值：
[1] "#FFFFFF" "#969696" "#FA00FF" "#F200FF" "#E600FF" # shortened for brevity
这是可以的，但我想得到的是RGBA (红色、绿色、蓝色、alpha)值，而不仅仅是RGB
浏览 8提问于2020-10-26得票数 1
回答已采纳
_rat = subprocess.check_output('gdalinfo -json ' + dataset_uri, shell=True)
我想提取.tif文件的栅格属性表。在上面的命令中，我能够将信息放入_rat，但不确定如何从_rat中提取rat部分。有什么建议吗？
浏览 15提问于2017-02-04得票数 0
回答已采纳
4回答
我有一个栅格堆栈stk，其中包含R中的三个栅格图像。下面是一个简单的示例
# set up a raster stack with three layers
> library(raster)
> r <- raster(nrows=10,ncols=10)
> r[] <- rnorm(100)
> stk <- stack(r,r,r)
# layer names are set by default
> names(stk)
[1] "layer.1" "layer.2" "layer.3"
浏览 1提问于2014-11-06得票数 23
1回答
我有几个文件夹，每个文件夹包含几个Landsat 8波段。我想遍历每个文件夹，并选择相同的4个特定波段。每个文件夹中的带按相同的顺序排列，因此可以根据它们的列表值进行查询。
这是我的代码，但不起作用。
library(raster)
library(rgdal)
setwd("Z:/TasseledCap")
folders <- list.dirs()
for(f in 1:length(folders)){
rasterlist<-list.files(folders[f], full.names=FALSE)
rasterlist
B <-
浏览 3提问于2017-02-03得票数 0
1回答
我正在通过一个未标记的预测函数处理一个大型数据堆栈(pm)。处理后，我将使用paroutPred、SE、for和as作为模板，将预测的、SE的、下的和上层的数据粘在一起，并将这些栅格叠加起来进行绘图。
我现在的代码在下面。foreach循环似乎运行得很好，并且我得到了必要的变量。说到底，compareRaster(rs，pm)的结果是正确的。只有不同的值。
在执行plot(rs)之后，四个栅格绘制，但它们都是顶起的，如下所示：。
我还没搞清楚到底出了什么问题。我在一个for循环中按顺序运行它(3天处理时间)，以确认光栅输出是错误的。
有人对我的问题有任何见解吗？我们非常感谢所有的帮助。
par
浏览 4提问于2021-02-07得票数 1
回答已采纳
1回答
我想同时遍历两个目录。两个目录都包含光栅图像，一个是初始图像，另一个是这些图像的掩码.
我想在循环之外访问它们，因为我必须做一些计算，但也保留我的文件夹结构。
我想，这些循环不一定要相互嵌套，但是我如何访问循环之外的文件呢？在外部创建一个空堆栈/光栅，并用循环的堆栈/光栅的值填充它，这似乎不是诀窍。
到目前为止我得到了这个：
library(raster)
files1 = list.files(path2, pattern = "*.tif", full.names = TRUE)
files2 = list.files(path4, pattern = "*_cm
浏览 0提问于2018-03-19得票数 0
回答已采纳
我有4个子文件夹，包含5个具有连续值的栅格。因此，构建一个带有"for“函数的循环以：
列出了这些光栅文件，在每个文件夹中堆叠这些文件，即4个光栅堆栈对象(包含5个光栅)，I为转换二进制光栅中的连续光栅设置了一个按钮，最后我使用wirte.raster函数编写了二进制光栅。
我的问题在第4步。尽管我在writeRaster函数中使用了参数“writeRaster= T”，但是保存的光栅是一个带有5个二进制光栅的光栅堆栈。我要把它写在栅格，文件，每个波段
如果有人给我任何见解，我真的很感激
setwd("Vole_raw_mean_Present/")
sub <
浏览 1提问于2020-04-15得票数 0
回答已采纳
1回答
注意:当创建一个从未被LibTiff.Net库更改过的新tiff映像时，我可以使用LibTiff.Net库添加/更改标记，而不会出现任何问题(第一次)。
在我使用LibTiff.Net库修改了tiff标记之后，当我再次尝试修改它们时，图像就会损坏。标签看起来很好，但图像本身显示为黑色图像(损坏)。
这是我当前在我的类中使用的标记扩展器...
private static Tiff.TiffExtendProc m_parentExtender;
private static List<TIFFTAGS_TAG> m_lTagsToWrite;
private static voi
浏览 0提问于2012-08-12得票数 1
回答已采纳
3回答
我正在尝试通过find复制由REGEX表达式选择的图像。下面的命令会运行，但不会复制文件。如何将循环变量(例如读取文件行)传递给find -regex参数？
cat hit_images_regex_list.txt | while read line;
do
find . -regex "${line}" -exec cp --parents {} /destination_dir \;
done
我从中提取REGEX表达式的hit_images_regex_list.txt文件如下所示：
".*B - 12.*tif"
".*D - 09.*ti
浏览 29提问于2020-02-13得票数 0
我正在处理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
1回答
我使用R重命名了330个文件，我在.csv文件中保存了原始名称和新名称。我使用了一个脚本，它不会出现错误，但不会更改名称。
下面是一个新名称的示例：(Df1)
D:\Modis_EVI\Original\EVI_Smoothed\ MODIS_EVI_20010101.tif
D:\Modis_EVI\Original\EVI_Smoothed\ MODIS_EVI_20010117.tif
D:\Modis_EVI\Original\EVI_Smoothed\ MODIS_EVI_20010201.tif
以及原始名称的示例：(Df2)
D:\Modis_EVI\Original\EVI_S
浏览 9提问于2016-06-24得票数 1
回答已采纳
1回答
我有两个JP2格式的光栅文件。我需要将两者结合起来，并针对波段执行计算。在.NET和C#中有什么方法可以做到这一点吗？我看到的大多数关于执行此操作的参考资料都使用了python中的Gdal的calc函数。
我已经尝试使用Gdal.Core和Gdal.Core.WindowsRuntime，但我没有看到Calculate调用的任何包装器。以前有没有人尝试过这样做，如果有，你是如何进行调用的，或者你使用的是什么库？
谢谢,
浏览 14提问于2020-03-10得票数 0
我有一个randomForest模型，我试图预测回到研究领域。我能够在没有错误的情况下运行预测，我可以看到它会产生光栅，但是当我尝试运行writeRaster()或plot()时，我得到了错误Error in .local(.object, ...)，而没有任何更多的信息。下面是代码(注意，envStack是环境预测器的光栅堆栈，rf.full是randomForest模型对象)：
library(here)
library(raster)
library(randomForest)
setwd(here("species", model_species, "outpu
浏览 0提问于2020-02-17得票数 3
1回答
我是r的新手，所以我还在摸索基础知识。我正在尝试获取多个shapefile的列表，并使用它们来裁剪几个栅格堆栈。然后使用裁剪的光栅堆栈，并从它们创建直方图。我知道我必须列出shapefiles的列表，我已经想好了怎么做。我还想出了如何单独裁剪光栅和制作直方图。不幸的是，这个循环把我搞糊涂了。有人有什么建议吗？
这是我到目前为止所知道的：
library('rgdal')
library('raster')
library('landsat')
library('ggplot2')
library('rasterVis
浏览 0提问于2016-01-11得票数 0
我有几百个Landsat8场景，每个有12个乐队。--我把它们都存储在一个文件夹中。，现在我尝试把它们堆在一起，以便更容易地对索引进行批处理。这是im使用的代码：
#Create List of files (whole path is used)
L8files = list.files(path =
"C:/Users/Felix/Desktop/Bachelorarbeit/Daten/R_Test/Only_TIF",
full.names=TRUE)
#Write Names of Files in List
getprefi
浏览 0提问于2018-07-28得票数 0
回答已采纳
1回答
我正在尝试将多个.tif堆栈(每个堆栈包含40幅图像)合并到一个tiff堆栈中。我更愿意使用python来完成这个任务。到目前为止，我尝试过的是(请记住，我编写代码的经验不多，如果我错过了一些显而易见的东西，很抱歉)：
import numpy as np
from skimage import io
im1 = io.imread('filename1.ome.tif')
for i in range(2,10):
im = io.imread('filename'+str(i)+'.ome.tif')
im1 = np.co
浏览 1提问于2017-11-08得票数 1
回答已采纳
由于数字的内部表示的经典问题，我认为在写入外部文件时会修改栅格值：
library(terra)
r <- rast(ncol=2, nrow=2, vals=c(5.3, 7.1, 3, 1.2))
sprintf("%.8f",values(r))
# [1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"
t <- writeRaster(r, "test.tif", overwrite=TRUE)
sprintf(&
浏览 0提问于2020-12-19得票数 1
我正在尝试寻找一个apply()类型的函数，它可以运行一个在两个数组而不是一个数组上操作的函数。
有点像：
apply(X1 = doy_stack, X2 = snow_stack, MARGIN = 2, FUN = r_part(a, b))
数据是来自Landsat的波段数组堆栈，这些数组使用rbind堆叠在一起。每一行都包含来自单个平铺的数据，最后，我需要对此堆栈中的每一列(像素)数据应用一个函数。一个这样的堆栈包含每个像素上是否有雪，另一个堆栈包含该行的日期。我想在每个像素上运行一个分类器(rpart)，让它识别每个像素一年中的无雪日。
我现在做的事情非常愚蠢：mapply(pas
浏览 1提问于2011-02-23得票数 5
回答已采纳
1回答
我有一个ndvi光栅。我必须使用阈值conditon.When来提取值，我绘制了提取的栅格的图形，x，y坐标不会太长。 我已经尝试了ggplot()，也使用了rastertopoints和rasterfromXYZ来绘制栅格，但都是徒劳的。代码如下。无法绘制s。xy协作者不会太长。 library(raster)
library(rgdal)
raster <- ("D:/Project/ndvisoybean.tif")
plot(ndvi)
ndvi[is.na(ndvi)] <-0
s <- ndvi[ndvi@data@values<maxVa
浏览 39提问于2018-12-29得票数 0
回答已采纳
我从谷歌地球引擎下载了哨兵2堆栈。对于出口，我使用了分辨率(scale: 10)，而有些波段应该具有20m分辨率(SWIR)或60m分辨率(例如Cirrus)：
Export.image.toDrive({
image: exportImgSentinel,
description: 'Sentinel',
folder: 'Sentinel',
fileNamePrefix: 'Sentinel',
region: region,
scale: 10,
crs: 'EPSG:4326'
}
浏览 6提问于2019-10-23得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
