如何叠加前哨2波段和前哨2波段的NDVI波段以改进qgis中的植被分类
在QGIS中，要叠加前哨2波段和前哨2波段的NDVI波段以改进植被分类，可以按照以下步骤进行操作：
- 首先，确保已经加载了前哨2波段和前哨2波段的NDVI波段的遥感影像数据。
- 打开QGIS软件，点击菜单栏中的"Raster"，选择"Raster Calculator"。
- 在弹出的对话框中，输入表达式来叠加NDVI波段。表达式的格式为："(前哨2波段 + 前哨2波段的NDVI波段) / 2"。这样可以将两个波段的像素值相加后再除以2，得到叠加后的NDVI波段。
- 点击"OK"按钮，QGIS会生成一个新的图层，即叠加后的NDVI波段。
- 接下来，可以使用该叠加后的NDVI波段进行植被分类。在QGIS中，可以使用分类算法（如最大似然分类、支持向量机等）对图像进行分类。具体操作可以参考QGIS的相关文档或教程。
需要注意的是，叠加前哨2波段和前哨2波段的NDVI波段只是一种改进植被分类的方法之一，具体效果还需要根据实际情况进行调整和评估。
推荐的腾讯云相关产品：腾讯云地理信息系统（GIS），该产品提供了丰富的地理信息处理和分析功能，可用于处理和分析遥感影像数据。产品介绍链接地址：https://cloud.tencent.com/product/gis
相关·内容
我正在尝试使用前哨2图像进行植被分类，我知道我也可以应用ndvi来改进我的分类。现在我已经从前哨2创建了ndvi，并且我已经将其堆叠在一起，结果光栅图像是黑色的。我想知道是不是应该这样，或者有什么我需要做的。我正在使用qgis
浏览 16提问于2019-12-03得票数 0
1回答
我正在与谷歌地球引擎工作，我试图提取/过滤(剪辑)像素在一个波段使用另一个图像(波段)。我计算了NDVI并创建了一个阈值，它呈现了一个NDVI > 0.3的图像，但是我希望提取可见的NIR波段中的相应像素。
这是代码片段。
var s2 = ee.ImageCollection('COPERNICUS/S2');
var s2_filtered = s2.filterDate('2017-01-01', '2017-12-31')
.filterBounds(geometry) //custom Geom
浏览 0提问于2018-03-15得票数 0
回答已采纳
1回答
我想开发一种工作流程，通过遥感检测土豆植物中的一种特定疾病。
我在一架在植物上空5米高度飞行的无人机上安装了一台多光谱相机，从而获得了土豆田的图像。
多光谱相机有5个波段，即:蓝、绿、红、近红外和RedEdge。
我已经将所有波段的DN(原始数字)值转换为反射率值。
我首先训练支持向量机从植物中分割土壤，然后应用SAVI(土壤调整植被指数)来细化植物中的土壤分割。
现在，我想应用NDVI(归一化差异植被指数)来确定植物的健康程度。
这是要遵循的正确方法吗？NDVI是否适用于仅在5米高拍摄的图像？或者还有更好的方法吗？
致以最良好的问候。
浏览 18提问于2018-02-15得票数 1
1回答
我有两组卫星图像。一个是一组经典的卫星图像，另一个是一组红外卫星图像。我正在尝试用检测黄色多边形内和红色区域外的植被。
可见图像
红外图像
根据图像记录，红外图像的颜色光谱发生了变化:红外波段占据红色波段，红色波段占据绿色波段。
为了计算NDVI图像，我执行以下操作：
# Images are in BGR color space.
ndvi = (img_irc[:,:,2].astype(int) - img_visible[:,:,2].astype(int))/(img_irc[:,:,2] + img_visible[:,:,2] + 0.001)
然
浏览 3提问于2021-09-18得票数 0
我试图对NOOA NDVI应用一条规则，这个已经被分类的区域只显示茂密的植被区域，但是它说ndvi不是一种功能，任何帮助都是非常感谢的。以下是代码：
//collecting data using ndvi
var good = ee.ImageCollection('NOAA/CDR/AVHRR/NDVI/V5')
.filter(ee.Filter.date('2018-05-01', '2018-06-01'))
//.median()
浏览 3提问于2022-01-21得票数 0
1回答
我的数据是这样的。我使用melt函数将数据排列如下
Legend variable value
1 Grassland NDVI 0.139
2 Grassland NDVI 0.285
3 Grassland NDVI 0.134
4 Grassland NDVI 0.243
5 Grassland NDVI 0.113
6
浏览 0提问于2016-12-14得票数 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
回答已采纳
1回答
我正在使用陆地卫星数据进行一些地形校正，使用Hantson & Chuvieco (2010)提出的方法，该方法涉及一种基于特定NDVI值的覆盖分离方法，以区分植被和非植被地区。
因此，我从一个多波段栅格对象变成两个多波段栅格对象，它们被分开处理，然后合并形成原始的、经过地形校正的场景。
问题是我在场景中有一些NA值，所以如果我不去掉它们，我会得到一个空白图像。我使用sum函数并将na.rm选项定义为TRUE来解决此问题，但结果是一个单波段栅格对象，而不是堆栈或栅格砖对象。
这是我的代码：
#topographic correction
VEGB1TC<- topocorr(VE
浏览 0提问于2016-07-07得票数 1
1回答
我正在致力于将多处理集成到一些数字图像处理工作流程中。以下脚本1)将图像波段转换为numpy数组，2)计算归一化差分植被指数(NDVI)，3)将numpy数组转换回光栅并写入磁盘。该脚本旨在了解使用多处理提高速度的想法。第一部分通过迭代工作空间，在每个栅格上执行处理，并写入磁盘(总时间=2分钟)来正常工作。第二个多处理部分无限期地“挂起”，并且不产生任何输出。在脚本的多进程部分，我哪里出了问题？
import arcpy, os, time
from multiprocessing import Pool
arcpy.env.workspace = r'C:\temp\tiles&#
浏览 1提问于2015-04-18得票数 1
1回答
我有一个如下所示的数据集：
Distance Mean SD Median VI Vegetation.Index Direction X X.1 X.2 X.3
1 10m 0.525 0.082 0.530 NDVI NDVI Whole Landscape NA NA NA NA
2 25m 0.517 0.085 0.523 NDVI NDVI Whole Landscape NA NA NA NA
3 50m 0.509 0.086 0.514 NDV
浏览 11提问于2016-11-14得票数 1
回答已采纳
我想使用监督分类来分类一个有明确时间模式的模式。例如，识别针叶林中落叶树的林分。NDVI会在落叶林中以一种很容易检测到的规律改变超时时间。我假设有一种简单的方法可以将时态数据集压缩成一幅图像，这样就可以在分类算法中使用图像中的波段。也许使用.map(....)
下面是一些代码来构建答案：
var startDate = '2016-05-01';
var endDate = '2016-09-01';
var lng = -122.3424; var lat = 37.9344; //SF
var region = ee.Geometry.Point(lng,
浏览 0提问于2018-04-27得票数 0
回答已采纳
1回答
我有三个时间序列数据。
归一化差异植被指数(NDVI)
降水
温度 所有这些都有自己的单位。现在我想找出NDVI与降水、NDVI和温度之间的相似性/相关性。基本上，我的目的是找出"NDVI与降水或温度更相关“。我应该将降水和温度正常化到NDVI值吗？
浏览 0提问于2019-02-06得票数 1
回答已采纳
3回答
我有一个数据框架与评估地点，计算不同的植被指数和不同的日期。我需要将所有信息连接到一个新的data.frame中，其中嵌入在植被指数中的日期信息包含在输出数据框架中的一个单独的列中。
我的数据框架有以下结构：
df.16 <- data.frame(ID=c("a","b","c"),
SUGAR=c(152232.92, 117937.06, 72080.81),
EVI_20160616_re=c(0.51, 0.59, 0.37), # The date
浏览 3提问于2020-07-04得票数 2
回答已采纳
1回答
我有两个ndvi图像具有值-1到1作为最小/最大值...我想从第一个图像中减去第二个图像，反之亦然，但当我这样做时，它给了我0 value...which的光栅，我认为它应该是这样的，因为我们从另一个波段中减去一个波段或从another....kindly中减去一个图像帮助我在哪里获得wrong.....here是我用来查看变化的代码...
b3<-raster(stackIm1,3) #band 3 of first input image
b4<-raster(stackIm1,4) #band 4 of first input image
ndvi1&
浏览 0提问于2013-02-21得票数 0
我正在处理一些tif文件，我必须绘制基于.tif文件的温度和蔬菜索引之间的依赖关系。那只是菲伊。现在我的编程问题。我使用python2.7 (x64)。我有一个大的ndarray表单NumPy库，包含温度值和第二个值(相同大小)和植被idex。mergedmask是我的掩码(和其他数组一样大小)，其中虚值意味着它是有效数据。
maxTS = np.amax(toa[mergedmask==False])
minTS = np.amin(toa[mergedmask==False])
maxVI = np.amax(ndvi1[mergedmask==False])
minVi = np
浏览 1提问于2015-06-01得票数 0
回答已采纳
我正试着在有很多云的地区(印尼云林)的谷歌地球引擎中制作一个无云的陆地卫星合成图。以前，我成功地实现了这一点，方法是制作一个最绿色的像素合成，其中我使用具有最高NDVI值的像素，以确保在合成图像中使用非云像素。
//Filter landsat 8 image collection by date, area
var collection = landsat
.filterBounds(bounds)
.filterDate(2016-08-01, 2016-10-31);
// Sort from least to most cloudy and get first (least clou
浏览 0提问于2018-09-27得票数 0
如果我的问题表达得不好，我很抱歉，但基本上我是在使用Python和GDAL来识别TIF文件中的植物(只使用RGB波段)。我查看了图像中的波段值，并查看了一些随机采样点，发现每个波段的植被倾向于落在特定的范围内(例如: 30-71蓝色，125-175绿色，等等)。所以我是这样写出来的，“数组”是一个包含三个波段的数组： vegetation = (arrays[0] > 67 & arrays[0] < 133) & (arrays[1] > 125 & arrays[1] < 175) & (arrays[2] > 30 &
浏览 19提问于2020-04-01得票数 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
回答已采纳
1回答
我感兴趣的是在现有的图像中添加一个高斯形状的物体，类似于附加的图像。我想将对象添加到的基本图像是8位无符号图像，其值范围为0-255。附图中明亮的物体实际上是一棵由归一化差异植被指数(NDVI)数据表示的树。所附的脚本是我到目前为止所拥有的。我如何添加一个高斯形状的物体(即一棵树)，其值从110-155到现有的NDVI图像？
，可以与此脚本一起使用来计算NDVI。
file = 'F:\path\to\fourband\image.tif';
[I R] = geotiffread(file);
outputdir = 'F:\path\to\output\
浏览 3提问于2014-07-05得票数 0
回答已采纳
我在一个与植被探测相关的项目中工作，我遵循Matlab的这个教程来计算NDVI。
但我犯了一个错误，我不知道为什么会这样。这是代码：
% Create a figure with a 1-by-2 aspect ratio
h = figure;
p = h.Position;
h.Position = [p(1,1:3),p(3)/2];
subplot(1,2,1)
% Create the scatter plot
plot(red, NIR, '+b')
hold on
plot(red(q(:)), NIR(q(:)), 'g+')
ax = g
浏览 1提问于2015-12-18得票数 0
回答已采纳
1回答
图像的数据类型是什么？
、、、
我推导了归一化差异植被指数( NDVI )，它是( NIR -R)/(NIR+R)的比值，其中近红外带是近红外波段，R是红带。这个索引范围从-1到1。所以我编写了一个pyopencl代码，下面是我所做的和观察到的。
Python代码：
import pyopencl as cl
import cv2
from PIL import Image
import numpy as np
from time import time
import matplotlib.pyplot as plt
#get kernel file
def getKernel():
kernel = open(
浏览 3提问于2020-06-02得票数 0
回答已采纳
我需要使用Sentinel 2图像进行分类，为此，我需要执行光谱可分性分析，以选择要使用的最佳波段和植被指数。因此，我需要计算训练站点的平均值和标准差。我试着使用这段代码，但结果没有用
// Get the Mean of the bands of the image for the polygons of the Vegetation class
var MeanTraining = Image.reduceRegions({
collection: Vegetation, // Vegetation is a FeatureCollection of polygons
r
浏览 5提问于2020-08-25得票数 1
1回答
我如何通过改编@Rodrigo E. Principe编写的这段漂亮的代码来提取集合上的索引植被点： Extract pixel values by points and convert to a table in Google Earth Engine 我尝试提取所有值，因此只有NDVI或EVI可以正常工作。
浏览 18提问于2019-10-18得票数 0
我需要从4波段图像中提取3个波段。我正在使用一个名为NumpyArrayToRaster()的函数，它只接受最多3个波段的图像。如何使其适用于4个波段的图像？
这是我现在的代码-
import arcpy
arcpy.CheckOutExtension("Spatial")
from PIL import Image
import matplotlib.pyplot as plt
import numpy as np
from skimage import io
from skimage.segmentation import quickshift
arcpy.env.ov
浏览 6提问于2015-10-26得票数 3
2回答
我使用scikit-learn库对卫星图像执行监督分类(支持向量机分类器)。我的主要问题是如何训练我的SVM分类器。我在youtube上看过很多视频，也读过一些关于如何在scikit-learn中训练支持向量机模型的教程。我看过的所有教程，都使用了著名的Iris数据集。为了在scikit-learn中进行有监督的支持向量机分类，我们需要有标签。对于Iris数据集，我们有一个Iris.target，它是我们试图预测的标签('setosa'，'versicolor'，'virginica')。通过阅读scikit-learn文档，培训过程非常简单。
浏览 0提问于2017-04-10得票数 5
回答已采纳
我的工作是遥感(陆地卫星5，30米分辨率图像)，我想使用陆地卫星波段的反射率和植被指数作为X预测因子，并使用一组基于地面的生物量测量作为Y输出来估计林木生物量。我在yaImpute package (R)中构建了kNN模型，并在gbm (R)软件包中使用了相同的数据集成boosting模型。
现在我想画一张光栅地图。函数AsciiGridImpute() (用于yai对象)和AsciiGridPredict() (用于其他对象，如GBM)必须提供帮助。但我得到了一个错误。
首先，我有带模型的yai-object和带观察值和预测值的impute对象：
test_yai <- yai(x =
浏览 4提问于2019-06-11得票数 0
1回答
谷歌地球引擎脚本
、、、
我需要从谷歌地球引擎将Geotiff图像文件转换为PNG，这可能吗？
我创建了以下代码，它生成一个geotiff以导出到Google，以便稍后下载并使用Qgis中的PNG映像。在声明参数、图像和变量之后，我将向您展示以下要点：
Map.addLayer(NDWI_table,{color:'red'}, 'vectors');
// Export result to Google Drive
Export.table.toDrive({
collection: NDWI_table,
description: 'NDWI_table'
浏览 10提问于2022-07-22得票数 0
我有两张光栅图像，一张来自第4波段，末尾是B4，另一张来自第5波段，最后是B5。我想将B5光栅子集为800x600，然后显示它并将其保存为GeoTiff。然后，我想计算NDVI (我假设我需要B4和B5来完成这个任务，但不确定)。然后，我想显示B5光栅的NDVI子集。显示它并将其保存为GeoTiff。
我如何创建像TIFF光栅图像的800 x 600像素子集？我也想拿那个TIFF，为这个子集生成一个NDVI映像。
注:我正在处理一张陆地卫星图像。图像在文件标题的末尾有B5。
我到目前为止所做的事：
import rasterio
from rasterio.windows import Win
浏览 3提问于2018-11-18得票数 2
我应该将哪些声明与逻辑函数/索引操作结合起来，这样Cython才能完成繁重的任务？
我有两个大小相同的大光栅。第一个数组包含植被索引值，第二个数组包含字段ID。目标是按田间平均植被指数。这两个数组都有讨厌的nodata值(-9999)，我想忽略它们。
目前，这个函数需要超过60秒的时间来执行，通常我不会介意，但我将处理可能的数百个图像。即使是30秒的改善也是很重要的。因此，我一直在探索Cython，以帮助加快速度。我一直在使用作为向导。
test_cy.pyx代码：
import numpy as np
cimport numpy as np
cimport cython
@cython.b
浏览 0提问于2018-02-17得票数 0
回答已采纳
1回答
利用栅格和数据的代数
、、、
我想用2个遥感植被指数(VIs)对多个月的多个林分进行植被健康预测。之前，我使用for()循环来迭代多波段光栅的列表，并使用给定的方程计算每个光栅(月)的两个VIs。然后我使用raster::extract()提取每个支架对应的像素。然而，我现在想在我对植被健康的预测中加入一些额外的变量，并且很难用相同的方法来集成它们，因为它们只是数据帧中的列，而不是栅格。我愿意用不同的方法来做这件事，我只是想不出任何办法。
示例：
#Part 1: Loading libraries and creating some sample data
library(sf)
library(raster)
lib
浏览 5提问于2022-06-07得票数 0
1回答
R和多列图形
、、
我就像R的新手，我在大学学习时就用过它，但它离我很远…
我有一个包含4列的表: vine ID，以及3列NDVI (一种植被指数)在3个日期的值。
ID 09052017 25052017 16062017
1 233 244 238
2 225 234 247
3 224 231 245
4 124 115 124
我知道如何读取我的表，如何使用它创建变量，如何选择列或行，如何创建plot(x,y)。
我的目标是为每个ID表示一条包含3个NDVI值的线，并在同一图形窗口中显示所有这些值
但
浏览 0提问于2017-07-11得票数 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
我有一个NDVI(Normalized植被指数)光栅对象，现在我想在相同的latitude.Like上计算像素值，我想知道30°N latittude.How中的像素值，我能用R做这个吗？
但我发现在“栅格”或"terra“包中没有任何功能可以做到这一点。
浏览 5提问于2022-11-08得票数 0
回答已采纳
我有一个包含NDVI数据和森林类型的数据集，我的NDVI列充满了NaN，我尝试平均2个NDVI列(植被指数从0到1)，然后按森林类型计算平均值。我的方法是创建第三列，然后计算平均值。类似于： data_full %>%
mutate(sum = sum(meanNDVIyear11,meanNDVIyear12, na.rm=TRUE)) %>%
mutate(sum = sum/2) %>%
group_by(OriginalForest) %>%
summarise(meanNDVI = mean(sum,na.rm=TRUE))
浏览 7提问于2020-04-20得票数 1
回答已采纳
我刚刚开始使用QGIS。我想自动计算，例如，
NDMI=(NIR-SWIR)/(NIR+SWIR)=(B08-B11)/(B08+B11)
适用于Sentinel-2和其他使用非RGB/近红外波段的索引，当使用SCP插件在QGIS中下载一批Sen2或陆地卫星-8图像时。
为此，我需要定义一个通用公式与通用名称的频带上传作为一个TXT文件在“频带计算”，将工作在所有的图像。例如，我知道Sen2 B08的特殊名称是"#NIR#"，但"#SWIR#“不适用于B11，可能是因为B11只是Sen2中可用的SWIR波段之一。
在SCP文档和互联网上，我找不到包含Sen2和Lands
浏览 2提问于2018-06-06得票数 1
1回答
我和R3.4.0一起工作。
我有多个单波段光栅，收集在一个文件夹，我想把它放在一起，在一个多波段光栅。按照其他帖子的指示，我写了以下简单的几行：
s2_list <- list.files(path="C:/Users/LAURA/Documents/S2", pattern=".tif$")
s2_stack <- stack(s2_list)
此过程给出了以下错误：
.local中的错误(.Object，.)：
.rasterObjectFromFile中的错误(x，band = band，objecttype = " RasterLa
浏览 0提问于2017-06-08得票数 1
1回答
我希望从陆地卫星数据中获得NDVI的月平均值，然后生成折线图。
我已经成功地对Landsat8进行了云掩蔽，并添加了NDVI波段。
使用另一篇文章中建议的代码，我能够在5年内每月创建一个NDVI值为1的图像集合。
然而，现在当我尝试用图表表示这个图像集合时，没有"system:time_start“可用，因此我无法生成时间序列。
任何帮助都将不胜感激。也许我必须用一种不同的方法来计算月平均NDVI值。
我的代码如下。
浏览 12提问于2019-11-15得票数 0
我试图从多光谱照片中计算多光谱植被指数，使用Python。
到目前为止，我已经成功地计算了NDVI，这是很容易实现的，因为网上有很多例子。
#Read the image file
image_file = "my_tif_image.tif"
with rasterio.open(image_file) as src:
band_red = src.read(3)
with rasterio.open(image_file) as src:
band_nir = src.read(4)
# Allow division by zero
nump
浏览 7提问于2020-11-20得票数 1
回答已采纳
从像'NDVI_TS.tif'这样的文件中的geotiff (NDVI的时间序列)堆栈中，我想得到单个的带名。例如：'Band 086: 20190803T004719'。例如，当我将堆栈加载到QGIS中时，我可以看到这一点。我需要追溯名称中的日期，在上面的示例中是2019-08-03。但是，我无法找到从Python访问它的方法。我可以通过索引访问乐队，但这并不能帮助我找到它们来自哪一天。
from osgeo import gdal
NDVI = gdal.Open('NDVI_TS.tif', gdal.GA_ReadOnly)
#for
浏览 1提问于2022-03-11得票数 0
回答已采纳
2回答
我正在使用HDF卫星数据来提取波段，并从中得出不同的植被指数。hdf数据中的每个波段都是灰度格式，是一幅灰度图像。经过HDF数据处理后，我可以使用色彩映射表转换为颜色(我使用jet作为色彩映射表)。我的疑问是如何在使用imwrite时将灰度图像转换为彩色映射。如何在imwrite中使用colourmap。我尝试了很多次，但输出只是全蓝色，这破坏了输出图像。请帮我做这件事。
浏览 0提问于2014-09-29得票数 0
我想在灰度tiff上使用NDVI索引实现渲染转换。我点击这个链接：- 什么类型的数据适用于它？
我尝试了Geoserver手册中提供的sld，但似乎对我的数据不起作用，但在应用该样式时，它不起作用
<NamedLayer>
<Name>cite:NDVI</Name>
<UserStyle>
<Title>NDVI</Title>
<FeatureTypeStyle>
<Transformation>
<ogc:Fun
浏览 11提问于2019-10-17得票数 0
我想使用工具箱sen2r来处理R中的前哨2 L2A数据。我已经手动下载了.SAFE格式的图像。 我已经使用s2_translate()将.SAFE格式转换为geotif： in_dir <- "D:/data/s2"
out_dir <-"D:/s2_geotifs"
## translate .safe to geotif
s2_example <- file.path(
in_dir,
"S2B_MSIL2A_20200525T104619_N0214_R051_T31UFT_20200525T133932.SAFE&
浏览 27提问于2020-08-18得票数 2
回答已采纳
1回答
我想用一张栅格表做一些计算。因此，我写了一个循环，到目前为止，我只使用光栅就可以正常工作了。现在我用栅格和数字vavlues进行计算，这会导致叠加函数中的一个错误。我的代码是：
列表文件
NDVI <- list.files(path =
"C:/Users/Felix/Desktop/Bachelorarbeit/Daten/Datenverarbeitung_R/L8_NDVI",
pattern = 'NDVI.tif$', full.names=T)
定义函数
fun_PV <- function(x,y,z)
浏览 0提问于2018-09-18得票数 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
我有一个18年的月平均NDVI值的时间序列。我想估计一下时间点，它证明了植被季节的开始或结束。我已经尝试过使用bfast，但这并没有给我带来我想要的结果。每一年都应该有很多断点。
有没有其他的R-包可以帮助你？
浏览 0提问于2019-11-18得票数 1
1回答
对于包含在该特定分类中的每个定制帖子，我必须为我的自定义分类(类别)添加简单的选项。
分类学名称：类型
分类法： Blues，Jazz，Rock
自定义帖子名称：波段
波段：第一波段、第二波段第三及第四波段
因此，基本上来看 Jazz 分类法选项，应该会给我Jazz乐队的选项。
带2 {text_field}
带四{text_field}
我不能用高级自定义字段来实现这一点，所以我想我需要写它？还是有一个插件来救我的命？如果没有，我应该从哪里开始寻找呢？
浏览 4提问于2014-03-27得票数 1
回答已采纳
1回答
R中的栅格重分类
、、
我有一个光栅图像，它包含NDVI的值在-0.2到0.6之间。
NDVI
#class : RasterLayer
#dimensions : 706, 953, 672818 (nrow, ncol, ncell)
#resolution : 30, 30 (x, y)
#extent : 368085, 396675, 2038125, 2059305 (xmin, xmax, ymin, ymax)
#crs : +proj=utm +zone=43 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +tow
浏览 1提问于2020-06-09得票数 0
回答已采纳
2回答
光栅重采样非常慢
、、、
)我想重采样哨兵-2图像，以便两个波段的分辨率适合。我使用的是raster：：raster(raster1，raster2，method="ngb")，但是执行过程非常长。它至少运行了15分钟，到目前为止我还没有得到任何输出。我不确定这件事是否会在某个时候结束。
我用的是哨兵2号的08波段(10米分辨率)和12波段(20米分辨率)。
用QGIS重采样同一光栅需要一分钟左右。我想问一下，这是否是此函数的常见问题，以及是否有任何解决办法或其他函数/包可用于此。
谢谢:-)
浏览 5提问于2021-08-18得票数 0
回答已采纳
我有一个栅格堆栈的列表，它是针对不同波段的时间序列(300+层)。时间序列是不规则的，因此我想创建基于中位数ndvi的年度组合。因此，从一年的可用图像中，选择具有(接近)中值ndvi的像素。对于其他波段，我想创建基于ndvi组合的年度组合。我试图制作一个掩码，其中每个像素都有用于中值ndvi合成的图像索引的值。我将应用于其他乐队，所以我有“相同”的年度组合为每个乐队每年。
不幸的是，我被困在做面具上了。我创建了一些虚拟数据，并以某种方式返回两个索引器(一个值为1-3，另一个为2-4)，而我希望返回一个索引器(值为1-4)。另外，我的函数无法处理NA值，并且将na.rm添加到calc函数并不能解
浏览 1提问于2018-11-15得票数 0
回答已采纳
2回答
我已经从非洲的MODIS下载了所有的植被数据，现在我正在尝试从图像中创建马赛克，并将它们保存为geo tif。在我尝试在MacPro 6核、12个线程上并行运行它之前，我完全可以做到这一点。代码运行，但它只使用1%的核心可用性，并且永远需要完成进程。
我真的需要帮助，我有超过70 lab的MODIS颗粒我需要转换成马赛克，如果我不使用我在实验室获得的所有计算机能力，这将永远需要。
以下是并行处理的代码部分：
有两个部分并行运行，1.选择NDVI；2.要创建马赛克，将它们投影并保存为.tif。
# register the cluster with 10 cores
registerDoPar
浏览 0提问于2015-07-23得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
