网站的开发费用吗,做一静态网站 多少钱,淘宝做网站的店,wordpress优酷视频插件【WRF数据准备】MODIS静态地理数据下载及制备 MODIS数据介绍数据下载数据拼接MRT工具介绍基于MRT软件完成数据拼接 格式转换#xff1a;tif文件转二进制格式编写INDEX修改GEOGRID.TBL以及namelist.wps修改GEOGRID.TBL修改namelist.wps 参考 MODIS数据介绍
MODIS-MCD12Q1 v061… 【WRF数据准备】MODIS静态地理数据下载及制备 MODIS数据介绍数据下载数据拼接MRT工具介绍基于MRT软件完成数据拼接 格式转换tif文件转二进制格式编写INDEX修改GEOGRID.TBL以及namelist.wps修改GEOGRID.TBL修改namelist.wps 参考 MODIS数据介绍
MODIS-MCD12Q1 v061
时间范围2001-2022年尺度数据空间分辨率500m
数据下载
1、点击【Access Data】如下 2、筛选数据 在其中选择对应时间、矩形选择方式并框选所需要的范围以下示例中国大部分区域 稍等片刻加载时间选择下载全部 第一步选择时间范围 第二步选择空间范围 网格范围经度106 120.58 纬度16.626 29.748 SW16.626,106 NE29.748,120.58 确定后界面如下 登录账号后下载数据。 需下载Earthdata Download软件进行数据下载界面如下 Earthdata Download中数据下载界面如下
数据拼接
MRT工具介绍
采用MRT(MODIS Reprojection Tool) 进行数据拼接MRT软件的下载及安装参见另一博客-【MODIS处理工具】MRT(MODIS Reprojection Tool) 下载及安装。
安装完成后打开MRT的安装路径文件夹打开bin文件夹双击ModisTool.bat或者ModisTool.jar即可进 入MRT的GUI界面 GUI界面如下
基于MRT软件完成数据拼接
1、导入并查看数据 点击【Open Input File】导入下载的hdf文件如下 可以使用【View Selected Tile】来查看所下载的数据的大致位置如果有出现在另一半球可以在此时删掉最终效果如下很奇怪15个文件就是有一个一直没下载好这会导致后续处理失败数据拼接要求导入的数据必须得是可以拼接起来的相邻影像 重新下载了所需的数据数据范围如下 2、 设置输出路径、文件格式、重采样方法 1Specify Output File设置输出文件路径。 2Output File Type设置输出文件的格式格式有 GEOTIFF、HDFEOS、RAW BINARY三种常用 GEOTIFF 格式。 3Resampling Type设置重采样方法有最邻近Nearest Neighbor、双线性插值Bilinear、三次卷积Cubic Convolution三种方法 3、设置投影、空间分辨率 1Output Projection Type中可以选择多种投影方式大家根据自己的需要来即可。此处笔者选择通用横墨卡托格网系统UTM。在Edit Projection Parameters中可以选择投影分带如果研究区正好在该分带中可以在UTM Zone中填入对应带号若研究区跨度很大可以不必填写。在Datum中选择WGS84坐标。 2Output Pixel Size中设置输出影像的空间分辨率单位是m在填写分辨率时最好升尺度不要降尺度比如说原影像分辨率为1000m那可以填1000、5000等等不要填500、250等小于原分辨率的数值。
其中波段使用LC_TYPE1重采样方式为Bilinear投影类型为Geographic像素分辨率为0.004991576420597609投影参数点进去选择WGS84 补充 在WGS84地理坐标系中分辨率通常以经纬度表示。MCD12Q1产品的原始分辨率是500米而0.004991576420597609度在赤道附近大约相当于500米的地面距离。 4、选择导出的数据集 本文利用的是MCD12Q1 v061数据该数据中包含多种数据集。但我若想只提取其中某一数据集该如何做 1首先将在 Selected Bands中的所有数据集通过 “ ” 符号全部移入 Avaliable Bands 栏中。 2将所需要的数据集通过 “ ” 符号再次导入到 Selected Bands栏中。 5、保存批处理的prm文件 点击Save Parameters file得到保存的prm文件。prm文件可以应用于后续的批处理。
6、Run运行 点击Run运行在弹出的窗口中看到Finished processing便表明运行完成啦。可以在所设置的路径中找到处理完成的影像和prm文件。到此需要处理单次MODIS数据的步骤就完成了。
格式转换tif文件转二进制格式
生成的TIF使用以下python程序处理成二进制格式Python代码如下
import rasterio as ra
from osgeo import gdaltifPath rC:/Users/NING\Desktop/modis2017/拼接/pinjie.LC_Type1.tif
rasterDataset ra.open(tifPath)xstart 00001
ystart 00001
xend str(rasterDataset.width).zfill(5)
yend str(rasterDataset.height).zfill(5)outputFile xstart - xend . ystart - yenddata rasterDataset.read(1)[::-1]
data.tofile(rC:/Users/NING\Desktop/modis2017/拼接/ outputFile)#INDEX数据提取
def get_tif_info(tif):dataset gdal.Open(tif)band dataset.GetRasterBand(1)xsize band.XSizeysize band.YSizemaxmin band.ComputeRasterMinMax()geo dataset.GetGeoTransform()print(地理变换6元素, geo)print(栅格矩阵的列数:, xsize, 栅格矩阵的行数:, ysize)print(最小最大值, maxmin)if __name__ __main__:get_tif_info(tifPath)编写INDEX
以上程序一同生成数据要素使用GetGeoTransform()输出tiff文件的地理信息六要素。
typecategorical
category_min1
category_max21
projectionregular_ll
dx0.004491576420597609
dy0.004491576420597609
known_x1.0
known_y11130.0
known_lat60.00009587599426
known_lon30.462848532138274
wordsize1
tile_x33291
tile_y11130
tile_z1
unitscategory
descriptionMODIS2017 landuse
mminluMODIFIED_IGBP_MODIS_NOAH
iswater17
islake21
isice15
isurban13各项参数含义
type为文件描述类型category_min分类代码的最小值category_max分类代码的最大值projection投影类型dx横向格点间的间隔即栅格影像的横向分辨率dy纵向格点间的间隔即栅格影像的纵向分辨率known_x指定一个标记点的横轴坐标known_y指定一个标记点的纵轴坐标known_lat标记点横向坐标的纬度known_lon标记点纵向坐标的经度tile_x横向格点数tile_y纵向格点数units格点值的单位description文件描述iswater水体类别的编号islake湖泊类别的编号isice冰川类别的编号isurban城市类别的编号
修改GEOGRID.TBL以及namelist.wps
修改GEOGRID.TBL
修改namelist.wps
参考
1、CSDN博客-制作MODIS静态地理数据并在WRF模式中使用 2、CSDN博客-MRT工具下载、安装、批处理MODIS数据
