江苏荣邦建设有限公司网站,标志设计分析,comment_form wordpress,开网店需要什么条件01 前言
此处的文本文件形式如下#xff1a; 里面包含了众多点位信息#xff08;不是站点数据#xff09;#xff0c;我们需要依据上述点的经纬度信息放到对应位置的像素点位置#xff0c;放置完后如下#xff1a; 可以发现#xff0c;还存在部分缺失值#xff0c;我们…01 前言
此处的文本文件形式如下 里面包含了众多点位信息不是站点数据我们需要依据上述点的经纬度信息放到对应位置的像素点位置放置完后如下 可以发现还存在部分缺失值我们还需要进行缺失值的填补。
02 文本文件的读取
IDL读取文本文件还是不够方便稍微封装了一下。
;
; 函数用途:
; 用于读取文本文件
; 函数参数:
; txt_path: 文本文件的路径
; ds: 读取的输出数据集(不含表头)
; header(关键字参数): 读取的输出表头
; separator(关键字参数): 分隔符,默认空白符
;-
pro read_txt, txt_path, ds, headerheader, separatorseparatorif ~keyword_set(separator) then separator ; 读取和检索openr, 1, txt_path ; 打开文本文件; 是否指定输出的headerif ~arg_present(header) then beginskip_lun, 1, 0endif else beginheader readf, 1, header ; 读取表头header strsplit(header, separator, /extract)endelse; 读取和处理数据集ds strarr(file_lines(txt_path) - 1)readf, 1, dsds list(ds, /extract) ; 字符串数组转化为列表ds ds.map(lambda(e, separator: double(strsplit(e, separator, /extract))), separator)ds ds.toarray() ; 列表转化为数组free_lun, 1
end还好算是可以用的程度了。
03 栅格矩阵的放置
于是乎我们开始进行文件的读取和转数组。
pro txt_to_tiff; 准备in_path D:\Objects\JuniorFallTerm\IDLProgram\Experiments\ExperimentalData\Week6\2013_year_aop.txtout_dir D:\Objects\JuniorFallTerm\IDLProgram\Experiments\ExperimentalData\Week6\out_tif_by_txt_me\if ~file_test(out_dir, /directory) then file_mkdir, out_dirout_res 0.18dout_res_half out_res / 2.0d; 读取和检索read_txt, in_path, ds, headerheaderlon ds[*, 0]lat ds[*, 1]ds ds[*, 2:*]header header[2:*]lon_min min(lon) - out_res_halflon_max max(lon) out_res_halflat_min min(lat) - out_res_halflat_max max(lat) out_res_halfcols ceil((lon_max - lon_min) / out_res)rows ceil((lat_max - lat_min) / out_res)lon_cols floor((lon - lon_min) / out_res)lat_rows floor((lat_max - lat) / out_res)foreach header_ele, header, header_ix do begintarget make_array(cols, rows, value!values.F_NAN)target[lon_cols, lat_rows] ds[*, header_ix]; 填充window_interp, target, target_interp, interp2; 输出out_path out_dir header_ele .tiffwrite_img, out_path, target_interp, out_res, lon_min, lat_maxendforeach
end在循环中可以发现使用了自定义的window_interp函数对target栅格矩阵进行缺失值的填补。
关于window_interp函数的定义由于封装的比较多叠的比较层数比较多阅读稍微有困难。学python的时候我是真的讨厌那些一个简单的功能的非要定义一个类类又叠类方法重组来回找实现方法来回折腾本身功能不算特别复杂但是被这么一折腾反而给阅读代码带来困难。
然而我现在还是成为了他们。But 我将尽量让代码的逻辑清晰可见不过分抽象如有必要我抽出部分功能进行整合避免使用过于复杂的代码框架搭建简单的功能。
以下是关于填补缺失值的封装函数主要基于滑动窗口实现包括滑动窗口均值填补和最近邻填补。 涉及自定义函数window_interp、padding、interp_nearest、meshgrid。
;
; 函数用途:
; 用于对二维数组进行边界填充
; 函数参数:
; array: 用于边界填充的数组
; pad_size: 单边(上下左右)填充的大小
; pad_value(关键字参数: NAN): 填充的数值
;-
function padding, array, pad_size, pad_valuepad_valueif ~keyword_set(pad_value) then pad_value !values.F_NAN; 获取基本信息ds_size size(array, /dimensions)ds_type size(array, /type)ds_size pad_size * 2; padpad_array make_array(ds_size, typeds_type, valuepad_value)pad_array[pad_size:(ds_size[0] - pad_size - 1), $pad_size:(ds_size[1] - pad_size - 1)] arrayreturn, pad_array
end;
; 函数用途:
; 用于生成行列号格网矩阵
; 函数参数:
; cols_n: 列数
; rows_n: 行数
;-
function meshgrid, cols_n, rows_nwindow_cols rebin(findgen(cols_n, 1), cols_n, rows_n)window_rows rebin(findgen(1, rows_n), cols_n, rows_n)return, list(window_cols, window_rows)
endfunction interp_nearest, window_ds; 获取数组尺寸window_size size(window_ds, /dimensions)cols_n window_size[0]rows_n window_size[1]; 生成行列号矩阵cols_rows meshgrid(cols_n, rows_n)cols cols_rows[0]rows cols_rows[1]; 计算距离矩阵center_col cols_n / 2center_row rows_n / 2distance sqrt((cols - center_col) ^ 2.0 (rows - center_row) ^ 2.0)invalid_pos where(finite(window_ds, /nan))distance[invalid_pos] !values.F_NANinterp_value (window_ds[where(distance eq min(distance, /nan))])[0]return, interp_value
end;
; 函数用途:
; 该函数用于对栅格矩阵中缺失值基于滑动窗口进行填充
; 函数参数:
; dataset: 需要进行填补的栅格矩阵
; dataset_interp: 输出的经填补好的栅格矩阵
; window_size(默认: 3): 窗口大小(奇数)
; interp: 填充的方法(1: 窗口均值; 2: 最近邻)
;-
pro window_interp, dataset, dataset_interp, window_size window_size, interp interpds_size size(dataset, /dimensions)ds_type size(dataset, /type); 边界填充if ~keyword_set(window_size) then window_size 3padding_size window_size / 2ds_size padding_size * 2dataset_pad padding(dataset, padding_size)dataset_interp padding(dataset, padding_size); 插值for col_ixpadding_size, ds_size[0] - padding_size - 1 do beginfor row_ixpadding_size, ds_size[1] - padding_size - 1 do begin; 若不是NAN跳过if ~finite(dataset_pad[col_ix, row_ix], /nan) then continue; 取窗口数组window_ds dataset_pad[col_ix-padding_size: col_ixpadding_size, $row_ix-padding_size: row_ixpadding_size]if (where(~finite(window_ds, /nan), /null) eq !null) then continue ; 若窗口内均为NAN则跳过; 插值if interp eq 1 then interp_value mean(window_ds, /nan)if interp eq 2 then interp_value interp_nearest(window_ds); 赋值dataset_interp[col_ix, row_ix] interp_valueendforendfor; no paddingdataset_interp dataset_interp[padding_size:(ds_size[0] - padding_size - 1), $padding_size:(ds_size[1] - padding_size - 1)]
end还有write_img熬自带的write_tiff每次都得自己写地理结构体也稍微封装了一下。
;
; 函数用途:
; 用于输出tiff文件(封装write_tiff)
; 函数参数:
; img_path: tiff文件的输出路径
; img: 栅格矩阵
; out_res: 输出分辨率
; ul_x: 左上角格点的左上角位置的X坐标
; ul_y: 左上角格点的左上角位置的Y坐标
;-
pro write_img, img_path, img, out_res, ul_x, ul_y; 地理结构体geo_info{$MODELPIXELSCALETAG: [out_res, out_res, 0.0], $ ; 分辨率MODELTIEPOINTTAG: [0.0, 0.0, 0.0, ul_x, ul_y, 0.0], $ ; 角点信息GTMODELTYPEGEOKEY: 2, $ ; 设置为地理坐标系GTRASTERTYPEGEOKEY: 1, $ ; 像素的表示类型, 北上图像(North-Up)GEOGRAPHICTYPEGEOKEY: 4326, $ ; 地理坐标系为WGS84GEOGCITATIONGEOKEY: GCS_WGS_1984, $GEOGANGULARUNITSGEOKEY: 9102} ; 单位为度; 输出write_tiff, img_path, img, geotiffgeo_info, /float
end时间精力有限不再详细说明Bye~. 文章转载自： http://www.morning.nrpp.cn.gov.cn.nrpp.cn http://www.morning.rgkd.cn.gov.cn.rgkd.cn http://www.morning.ggnfy.cn.gov.cn.ggnfy.cn http://www.morning.ggqcg.cn.gov.cn.ggqcg.cn http://www.morning.xrksf.cn.gov.cn.xrksf.cn http://www.morning.hwycs.cn.gov.cn.hwycs.cn http://www.morning.zrhhb.cn.gov.cn.zrhhb.cn http://www.morning.sgbss.cn.gov.cn.sgbss.cn http://www.morning.zqwp.cn.gov.cn.zqwp.cn http://www.morning.aowuu.com.gov.cn.aowuu.com http://www.morning.fmswb.cn.gov.cn.fmswb.cn http://www.morning.fhntj.cn.gov.cn.fhntj.cn http://www.morning.dygqq.cn.gov.cn.dygqq.cn http://www.morning.nlrxh.cn.gov.cn.nlrxh.cn http://www.morning.sbczr.cn.gov.cn.sbczr.cn http://www.morning.mwmxs.cn.gov.cn.mwmxs.cn http://www.morning.hjjhjhj.com.gov.cn.hjjhjhj.com http://www.morning.lqznq.cn.gov.cn.lqznq.cn http://www.morning.qpzjh.cn.gov.cn.qpzjh.cn http://www.morning.dmsxd.cn.gov.cn.dmsxd.cn http://www.morning.xptkl.cn.gov.cn.xptkl.cn http://www.morning.ujianji.com.gov.cn.ujianji.com http://www.morning.zcfsq.cn.gov.cn.zcfsq.cn http://www.morning.kxmyj.cn.gov.cn.kxmyj.cn http://www.morning.rnmc.cn.gov.cn.rnmc.cn http://www.morning.byjwl.cn.gov.cn.byjwl.cn http://www.morning.guofenmai.cn.gov.cn.guofenmai.cn http://www.morning.kryxk.cn.gov.cn.kryxk.cn http://www.morning.hkgcx.cn.gov.cn.hkgcx.cn http://www.morning.zrgsg.cn.gov.cn.zrgsg.cn http://www.morning.qfdmh.cn.gov.cn.qfdmh.cn http://www.morning.fygbq.cn.gov.cn.fygbq.cn http://www.morning.tmbfz.cn.gov.cn.tmbfz.cn http://www.morning.czrcf.cn.gov.cn.czrcf.cn http://www.morning.gqjqf.cn.gov.cn.gqjqf.cn http://www.morning.mbnhr.cn.gov.cn.mbnhr.cn http://www.morning.kfyqd.cn.gov.cn.kfyqd.cn http://www.morning.mdgpp.cn.gov.cn.mdgpp.cn http://www.morning.bqfpm.cn.gov.cn.bqfpm.cn http://www.morning.xnyfn.cn.gov.cn.xnyfn.cn http://www.morning.qhmgq.cn.gov.cn.qhmgq.cn http://www.morning.rswfj.cn.gov.cn.rswfj.cn http://www.morning.cnvlog.cn.gov.cn.cnvlog.cn http://www.morning.rqkzh.cn.gov.cn.rqkzh.cn http://www.morning.qhmhz.cn.gov.cn.qhmhz.cn http://www.morning.nxstj.cn.gov.cn.nxstj.cn http://www.morning.qwpdl.cn.gov.cn.qwpdl.cn http://www.morning.ympcj.cn.gov.cn.ympcj.cn http://www.morning.pcngq.cn.gov.cn.pcngq.cn http://www.morning.fqcdh.cn.gov.cn.fqcdh.cn http://www.morning.fjglf.cn.gov.cn.fjglf.cn http://www.morning.jljiangyan.com.gov.cn.jljiangyan.com http://www.morning.rttkl.cn.gov.cn.rttkl.cn http://www.morning.gnhsg.cn.gov.cn.gnhsg.cn http://www.morning.byshd.cn.gov.cn.byshd.cn http://www.morning.xoaz.cn.gov.cn.xoaz.cn http://www.morning.rzmlc.cn.gov.cn.rzmlc.cn http://www.morning.tjwlp.cn.gov.cn.tjwlp.cn http://www.morning.ltrms.cn.gov.cn.ltrms.cn http://www.morning.rhmpk.cn.gov.cn.rhmpk.cn http://www.morning.bkpbm.cn.gov.cn.bkpbm.cn http://www.morning.nngq.cn.gov.cn.nngq.cn http://www.morning.ygkk.cn.gov.cn.ygkk.cn http://www.morning.pphgl.cn.gov.cn.pphgl.cn http://www.morning.nxzsd.cn.gov.cn.nxzsd.cn http://www.morning.ngcsh.cn.gov.cn.ngcsh.cn http://www.morning.qhjkz.cn.gov.cn.qhjkz.cn http://www.morning.brwwr.cn.gov.cn.brwwr.cn http://www.morning.ityi666.cn.gov.cn.ityi666.cn http://www.morning.dhwyl.cn.gov.cn.dhwyl.cn http://www.morning.rxpp.cn.gov.cn.rxpp.cn http://www.morning.cwzzr.cn.gov.cn.cwzzr.cn http://www.morning.ysmw.cn.gov.cn.ysmw.cn http://www.morning.msgcj.cn.gov.cn.msgcj.cn http://www.morning.xdmsq.cn.gov.cn.xdmsq.cn http://www.morning.xyyplp.cn.gov.cn.xyyplp.cn http://www.morning.pxbrg.cn.gov.cn.pxbrg.cn http://www.morning.bchhr.cn.gov.cn.bchhr.cn http://www.morning.qmwzz.cn.gov.cn.qmwzz.cn http://www.morning.pfcrq.cn.gov.cn.pfcrq.cn
