-
数据转换 HDF --> TIFF
- MOD09GA、MCD43A1
- 工具：MRT
- 拼接：H25V4, H26V4, H26V5, H27V5四个区域
- 投影转换：geographic
- 分辨率：0.01 decimal degrees
- MOD11_L2
- 工具：
MRTSwathIDL - ditto
- 对每个hdf文件，运行IDL转成dat，然后用ENVI打开dat转成tiff
- 工具：
- MOD09GA、MCD43A1
-
研究区域的确定
- 计算出研究区域对应的经纬度BBOX
- minLon = 105.78432759485935
- maxLon = 111.88621521850175
- minLat = 36.95825861049653
- maxLat = 41.45847714860637
- 对每一组数据
- L2 Swath数据转换后的图像：根据经纬度BBOX计算对应图像中的BBOX坐标
- L3 网格化的数据转换后的图像：直接使用0812数据计算的BBOX坐标
- 图像中任意（x, y）点的经纬度计算
from osgeo import gdal ds = gdal.Open('path/to/file') width = ds.RasterXSize height = ds.RasterYSize gt = ds.GetGeoTransform() # bbox minLon = gt[0] minLat = gt[3] + width*gt[4] + height*gt[5] maxLon = gt[0] + width*gt[1] + height*gt[2] maxLat = gt[3] # any point (xy以左上角为原点) lat = gt[3] + x*gt[4] + y*gt[5] lon = gt[0] + x*gt[1] + y*gt[2]
- 由经纬度坐标计算图像坐标，即以上公式的反过程
- 计算出研究区域对应的经纬度BBOX
2020-07-18 数据预处理
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
