ndvi最大值合成_envi ndvi最大值合成 - CSDN
精华内容
参与话题
  • ENVI 最大值合成(MVC)

    千次阅读 2019-06-21 22:52:02
    返回两个波段中的大 function bmax,b1,b2 b=(b1 ge b2) * b1 + (b2 gt b1) * b2 return,b end 转载于:https://www.cnblogs.com/EggKiller/archive/2008/09/27/1300066.html

    返回两个波段中的大值

    function bmax,b1,b2

        b=(b1 ge b2) * b1 + (b2 gt b1) * b2

        return,b

    end

    转载于:https://www.cnblogs.com/EggKiller/archive/2008/09/27/1300066.html

    展开全文
  • python_MOD13A3_NDVI 最大合成

    千次阅读 2018-09-18 16:57:52
    如图,输入的文件路径下存放NDVI数据和 EVI数据,因此需要用到正则表达式匹配文件名。 import arcpy from arcpy import env from arcpy.sa import * import re arcpy.CheckOutExtension("spatial") #...

      如图,输入的文件路径下存放NDVI数据和 EVI数据,因此需要用到正则表达式匹配文件名。

    import arcpy
    from arcpy import env
    from arcpy.sa import *
    import re
    
    arcpy.CheckOutExtension("spatial") #检查模块许可
    
    # Set environment settings
    env.workspace = r"G:\MOD13A3_NDVI\Extract_Tarim"  #NDVI数据存放目录
    outpath=r'G:\MOD13A3_NDVI\MaxMonth_NDVI\2002_MaxNDVI.tif'  # 年最大合成NDVI数据文件名
    
    list_raster=[] #空LIST 用于存放全部用于合成的月NDVI
    re_info=r'[2][0][0][2][.a-zA-Z0-9_]{16}[N][D][V][I].tif$' #正则表达式匹配字符串  用于匹配2002年12个月NDVI数据,文件后缀为tif
    
    for file1 in arcpy.ListFiles("*.tif"): #循环列出数据存放目录.tif文件
        s=file1
        if re.findall(re_info,s):   #使用re.findall 匹配2002年的月数据,用于合成2002年最大合成NDVI
            print s
            list_raster.append(s)  #将符合文件命名要求的2002年月数据名称放入列表
    print list_raster
    
    ### Execute CellStatistics
    
    outCellStatistics = CellStatistics(list_raster, "MAXIMUM", "DATA") # 合成2002年12个月的NDVI数据,方法为 MAXIMUM , 无数据的地方以有数据的最大值为准DATA
    
    ### Save the output 
    outCellStatistics.save(outpath)  #保存数据
    print'success'
    

    输出文件示例

     

    展开全文
  • 我这是处理一年的,要是多年的,需要对数据的名称进行统一格式,方便循环语句的编写 ...str="H:\\1999monthNDVI\\"+"NDVI" a1 = str+"01.tif" a2 = str+"02.tif" a3 = str+"03.tif" a4 = str+"04.tif" a5 = s...

    我这是处理一年的,要是多年的,需要对数据的名称进行统一格式,方便循环语句的编写

    import arcpy
    arcpy.CheckOutExtension("spatial")
    str="H:\\1999monthNDVI\\"+"NDVI"
    a1 = str+"01.tif"
    a2 = str+"02.tif"
    a3 = str+"03.tif"
    a4 = str+"04.tif"
    a5 = str+"05.tif"
    a6 = str+"06.tif"
    a7 = str+"07.tif"
    a8 = str+"08.tif"
    a9 = str+"09.tif"
    a10 = str+"10.tif"
    a11 = str+"11.tif"
    a12 = str+"12.tif"
    out = "H:\\1999monthNDVI\\"+"NDVImax_1999.tif"
    arcpy.gp.CellStatistics_sa((a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12), out, "MAXIMUM", "DATA")
    print(" has done")
     

    展开全文
  • Arcpy 实现NDVI最大值合成/平均值合成

    千次阅读 2019-04-10 21:52:19
    当研究NDVI产品数据时,需要日产品月最大值合成/月平均值合成,或者月产品数据季度最大值合成/平均值合成,ArcGIS很方便的实现这些功能,Arcpy辅助更加便捷的处理批量数据。 首先介绍一下Arcpy函数CellStatistics...

    当研究NDVI产品数据时,需要日产品月最大值合成/月平均值合成,或者月产品数据季度最大值合成/平均值合成,ArcGIS很方便的实现这些功能,Arcpy辅助更加便捷的处理批量数据。

    首先介绍一下Arcpy函数CellStatistics(像元统计):

    根据多个栅格数据计算每个像元的统计数据。

    可用的统计数据有:众数、最大值、均值、中位数、最小值、少数、范围、标准差、总和及变异度。

    用法

    • 使用此工具时输入栅格的顺序不重要。

    • 将多波段栅格指定为输入栅格或常量值(Python 中的 in_rasters_or_constants)之一时,将使用所有波段。

      要处理一系列来自多波段栅格的波段,可先用波段合成工具创建由这些特定波段组成的新栅格数据集,并在输入栅格或常量值(Python 中的 in_rasters_or_constants)列表中使用该结果。

    • 对于最大值、最小值、均值、中位数、众数、少数及总和这些统计数据类型,如果将单个栅格数据用作输入,则输出像元值将与输入像元值相同。对于范围和 STD这些统计数据类型,输出像元值全部为 0;对于变异度这种统计数据,输出像元值均为 1。

    python源代码:

    outCellStatistics = CellStatistics([inRaster01, inRaster02, inRaster03], "RANGE", "NODATA")
    outCellStatistics.save("C:/sapyexamples/output/cellstats")

    至关重要(加在关键代码之前):

    arcpy.CheckOutExtension("Spatial")

    效果展示:

     

    平均值合成:

    展开全文
  • ArcGIS Python实现Modis NDVI批量化月最大合成

    万次阅读 多人点赞 2014-08-29 14:46:19
    import arcpy arcpy.CheckOutExtension("spatial") Sname1="F:\\Modis_16\\1Moasic\\" Sname2=".1_km_16_days_NDVI.tif" out = "F:\\Modis_16\\2MVC\\ddwwq" for i in range(2004,2014): for j in range(1,13):
  • 目的:利用ArcGIS自带工具Cell Statistics实现多个栅格图层最大值合成NDVI年或月的最大值。 工具:ArcToolbox>>Spatial Analyst Tools>>Local>>Cell Statistics 问题:NDVI栅格影像值为-1~1之间的float浮点型,...
  • 一年中的12个月份的月最大合成(mvc)数据放在“F:\\Vegetation Change\\Data\\GIMMS Data\\1MVC\\",数据名称格式为mvc_198801,mvc_198802........mvc_198812。处理年份为1981-2006,代码为: import arcpy arcpy....
  • 接下来把tif图像转成我们习惯使用的文本格式,即经度,纬度,ndvi。共三列。当然也可以根据个人偏好制作为nc文件等(制作nc文件可参见我的博客)。 第一步:把tif文件转为ascii格式的文本文件。所谓ascii格式的文本...
  • GIMMS NDVI3g V1.0发布

    千次阅读 2018-09-12 11:34:12
    GIMMS NDVI3g V1.0发布。 如需要Geotiff格式文件,请联系作者,qq:350150279 中国区域的tif格式数据分发地址: https://pan.baidu.com/disk/home#/all?vmode=list&path=%2FGeoData%2FNDVI3g_V1.0。   ...
  • 基于Sentinel-2(哨兵2)遥感卫星的植被覆盖度提取

    千次阅读 多人点赞 2019-08-31 14:57:05
    基于Sentinel-2(哨兵2)遥感卫星的植被覆盖度提取 一、Sentinel-2简单介绍 哨兵-2系列卫星的主要载荷多光谱成像仪拥有13个谱段,从可见光到近红外再到短波红外,空间分辨率从10m到60m,能实现全所未有的陆地海洋...
  • 04-SNAP处理Sentinel-2 L2A级数据(二)

    千次阅读 多人点赞 2020-04-23 17:30:08
    04-SNAP处理Sentinel-2 L2A级数据(二)前言数据集波段指数计数植被指数计算纹理特征提取主成分变换波段叠加结语 (原创文章,转载请注明来源!) 前言 03篇中谈到了利用SNAP对Sentinel-2进行的一些预处理,但是...
  • Python实现NDVI的MVC算法

    千次阅读 2018-11-28 00:07:45
    因为所下数据为NDVI旬数据,合并为月数据时采用国际通用最大值合成法(MVC)将旬数据合并成月数据,因数据较多,遂考虑用代码解决。主要分为遍历文件,读取制定名称的栅格文件,以及ArcGIS工具箱中的Raster ...
  • 本文为作者硕士学位论文遥感物候研究数据处理过程总结。 目录 一、数据下载 ...因文章研究的需要和数据时间段的限制,我从EcoCast网站上下载到的GIMMS 3G NDVI数据序列是30年(1982-2013年)。时间分辨率...
  • 像元二分模型是一种简单实用的遥感估算模型,它假设一个像元的地表由有植被覆盖部分地表与无植被覆盖部分地表组成,而遥感传感器观测到的光谱信息也由这2个组分因子线性加权合成,各因子的权重是各自的面积在像元中...
  • 像元二分模型将像元的植被覆盖结构分为纯像元与混合像元两类,纯像元植被完全覆盖覆盖度为1,混合想要由植被与废止被构成,其NDVINDVINDVI是植被覆盖部分的NDVINDVINDVI与非植被覆盖部分NDVINDVINDVI的线性加权...
  • 基于改进的 CASA 模型反演 NPP 数据为地理空间数据云提供的 TM 影像,其具体信息如图 1 所示。该实验中所使用的 TM 数据已进行过 QUAC 快速大气校正,来消除大气和光照等因素对地物反射的影响。...
  • 时空分析

    千次阅读 2017-08-02 22:50:22
    1.MODIS NDVI数据产品 1.1 数据介绍 MODIS产品的详细介绍:https://wenku.baidu.com/view/f41dff553186bceb18e8bb59.html 我采用的数据为02-12年的生长季MODIS ndvi数据 参考标准的数据说明: 以及具体编号说明...
  • 我们可以选取2018年5月初华东地区MODIS中的MOD11A2和MOD13A2的16天合成LST和NDVI产品数据,下载地址:MODIS数据下载 网站下载数据需要注册,在此声明,下载数据不需要“科学上网”,但注册需要谷歌Gmail邮箱,注册...
  • 第一部分:计算2015-2020年(每年8月份合成)的NDVI均值 首先建立一个空的列表用来存储影像 var col=ee.List([]); 筛选2015-2020年五年的影像集(每年选择8月份的影像合成) for(var i=2015;i<2020;i++){ var ...
1 2 3 4 5
收藏数 100
精华内容 40
关键字:

ndvi最大值合成