1. 内容
——————————————————————————————————————————
——————————————————————————————————————————
这里主要没有用到什么新的知识点,就是先获取数据,然后找到经纬度的最大值最小值,以此得到行列数去构建装数据的数组,然后数据按顺序放入到数组中。另外,如果你的数据源明确说明数据的分辨率是准确的而不是近似,那么你可以不进行均值插值,当然为了确保你可以先出图看看缺失值是否存在再来判断是否需要均值插值。做好均值插值后最后输出即可。
——————————————————————————————————————————————————————————————————————————————————————
2. 编程
pro geo_out_tiff, lon, lat, geo_data, resolution, geo_data_out, geo_info
; 本程序用于将经纬度、xx数据集转化为可输出为tiff的数据
; 确定最大最小的经纬度
lon_min = min(lon)
lon_max = max(lon)
lat_min = min(lat)
lat_max = max(lat)
; 获取转化为tiff数据所需要的行列数
column = ceil((lon_max - lon_min) / resolution) ; ceil()函数表示向上取整
row = ceil((lat_max - lat_min) / resolution)
; 得到每一个点的行列位置
column_array = floor((lon - lon_min) / resolution) ; 向下取整
row_array = floor((lat_max - lat) / resolution)
; 创建装数据的box
data_box