- 博客(13)
- 收藏
- 关注
RSS订阅
原创 python活用gdal库进行批量重投影、重采样、裁剪
还在为arcgis批量处理烦恼吗？快来看看我的python gdal 批量处理方案吧，仅需10s解决arcgis1h的工作量，解决你超多时间。
2025-01-24 11:57:21
1360
原创 日分辨率的干旱指标SPEI计算，干旱事件识别
本文给出了SPEI在日、8day、月分辨率下的计算，允许读者自由选择：1. 计算潜在蒸散发的方式；2. 需要的分辨率，即日、8day、月三种；3. 不同的累计水平衡的时间，即不同的时间尺度；4.SPEI的拟合分布及拟合分布的参数的优化方法。
2024-09-15 11:57:06
4595
14
原创 基于Landsat 1-8提取点或区域的覆盖度
植被覆盖度（fractional vegetation coverage）是一些生态水文过程模型的重要输入参数，本文的主要目的是提取年均覆盖度的时间序列，以约束生态水文模型。
2024-09-13 10:05:52
1139
原创 Python进行nc文件、tif文件、ndarray数据的相互转换
netcdf文件，一般简称nc文件，包含经纬度、时间三个维度的文件格式。tif文件，常用的geotiff文件，包含经纬度和波段数三个维度的信息。作为学术上公开数据集中常用的两种数据保存格式，tif和nc文件相互转换的熟练进行，是增加工作效率的保障之一。下面的代码主要分为4个部分：1. 将numpy数组格式的ndarray数据保存为nc文件；2. 将numpy数组格式的ndarray数据保存为tif文件；3. 直接将nc文件转化为tif文件；4. 直接将tif文件保存为nc文件。
2025-01-24 11:09:05
968
原创 基于MCD43A4产品提取8day的NDVI、EVI，并批量下载到Drive
在之前我的一篇SCI文章中，用到了MCD43A4数据去计算kNDVI和NIRv指数，因为这两种指数比传统的NDVI更能反映真实的GPP变化(Badgley et al., 2017;本文介绍的是NDVI、EVI的提取方法，如果需要还可根据我文章中的公式计算得到kNDVI和NIRv。使用数据NIR:这个产品是daily分辨率的，需要进行一定的处理以排除噪声，我这里是取了8天的平均以反映草地生产力迅速变化的特征。对于森林，8day也是一个很好的指标，16day可能也行，因为森林一般变化没这么快。
2024-09-13 14:48:15
1310
2
原创 基于Landsat和MODIS进行多波段融合（反射率、地表温度）
landsat去云后，缺失数据比较多，希望用MODIS数据进行填补，更准确的术语是：数据融合（data fusion）。具体来说是用，MCD43A4的反射率、地表温度去插补Landsat 8 的反射率和地表温度，示例的时间选为2017年的5-7月。
2024-09-13 14:23:42
2594
2
原创 基于MODIS提取某个点的覆盖度
植被覆盖度（fractional vegetation coverage）是一些生态水文过程模型的重要输入参数。前一篇博客介绍了基于Landsat提取覆盖度的方法，本文是为了对覆盖度的结果进行多个产品之间的验证。
2024-09-13 10:17:51
503
原创 NDVI 数据产品和常见传感器说明（Modis,TM,AVHRR,SPOT）
1，modis (moderate resolution imaging spectroradimeter)Modis是由nasa发射的，共有36个通道，从可见光到热红外波段，在1999年发射的terra， 2002年5月发射Aqua， Terra每天经过2次，分别为早上10：30和22：30， aqua每天同样两次，1：30am和13：30Modis主要提供以下变量每日地表反射率，雪盖，NDVI、EVI资料，地表温度和大气每日变量2. TM (Landset Thematic Mapper)
2021-10-30 11:00:17
9374
原创 HDF官方文档阅读
1. HDF介绍1.1 HDF：分层数据格式，用于存储和分发科学数据的自我描述、多对象文件格式。1.2 6种基本数据类型： 1）光栅数据 raster image 2）调色板 palette 3）科学数据集 scientific data set 4）注释 annotation 5）虚拟数据 vdata 6）虚拟组 group1.3 HDF4与HDF5区别 HDF4...
2021-10-30 10:31:16
514
原创 HDFView使用
General Object InfoData Type: HDF4 SDS在MODIS 1B数据中主要有两种数据对象SDS和VdataCompression：Gzip（level = 5）压缩比(1~9),越小压缩效果越差,但是越大处理越慢Fill Value：填充值Show Data with Options Spreadsheet：表格 Image:Select palett...
2021-10-30 10:30:48
1336
原创 netCDF4阅读笔记
建立、打开、关闭netCDF4文件mode = 'w' or 'r' or 'a'，w write r read a alter.from netCDF4 import Datasetrootgrp = Dataset("test.nc", mode = "w", format="NETCDF4")群组>>> rootgrp = Dataset("test.nc", "a")>>> fcstgrp = rootgrp.createGroup("fo
2021-10-29 11:30:49
594
多时间尺度、多分辨率、多PET计算方式的 日/周/月干旱指标SPEI计算代码及测试文件
2024-09-15
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人