氣象數(shù)據(jù)處理流程

編號(hào)：時(shí)間：2021年x月x日書山有路勤為徑，學(xué)海無涯苦作舟頁碼：第頁氣象數(shù)據(jù)處理流程數(shù)據(jù)下載登錄中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)注冊(cè)用戶選擇地面氣象資料選擇中國地面國際交換站日值數(shù)據(jù)選擇所需數(shù)據(jù)點(diǎn)擊預(yù)覽（本次氣象數(shù)據(jù)為：降水量、日最高氣溫、日最低氣溫、平均濕度、輻射度、積雪厚度等；地區(qū)為：黑龍江省、吉林省、遼寧省、內(nèi)蒙古）下載數(shù)據(jù)并同時(shí)下載文檔說明網(wǎng)站數(shù)據(jù)粘貼并保存為TXT文檔建立屬性庫存儲(chǔ)后的TXT文檔用Excel打開并將第一列按逗號(hào)分列站點(diǎn)數(shù)據(jù)處理由于站點(diǎn)數(shù)據(jù)為經(jīng)緯度數(shù)據(jù)為方便插值數(shù)據(jù)設(shè)置分辨率（1公里）減少投影變換次數(shù)，先將站點(diǎn)坐標(biāo)轉(zhuǎn)為大地坐標(biāo)并添加X、Y列存儲(chǔ)大地坐標(biāo)值后將各項(xiàng)數(shù)據(jù)按照站點(diǎn)字段年月日合成總數(shù)據(jù)庫（注意：數(shù)據(jù)庫存儲(chǔ)為DBF3格式，個(gè)字段均為數(shù)值型坐標(biāo)需設(shè)置小數(shù)位數(shù)）為填補(bǔ)插值后北部和東部數(shù)據(jù)的空缺采用最鄰近法將漠河北部、富錦東部補(bǔ)齊2點(diǎn)數(shù)據(jù)。利用VBA程序Subwe()i=6Forj=1To30Windows("chengle.dbf").ActivateRows("1:1").SelectSelection.AutoFilterSelection.AutoFilterField:=5,Criteria1:=iSelection.AutoFilterField:=6,Criteria1:=jCells.SelectSelection.CopyWorkbooks.AddActiveSheet.PasteWindows("chengle.dbf").ActivateRows("1:1").SelectSelection.AutoFilterWindows("book"+CStr(j)).ActivateRange("A1:n100").SelectRange("I14").ActivateCells.EntireColumn.AutoFitChDir"C:\DocumentsandSettings\王\桌面"ActiveWorkbook.SaveAsFilename:="C:\DocumentsandSettings\王\桌面\6\"&InputBox("輸入保存名",Title="保存名字","20070"+CStr(i)+"0"+CStr(j)),_FileFormat:=xlDBF4,CreateBackup:=FalseActiveWorkbook.CloseSaveChanges:=TrueNextjEndSub將數(shù)據(jù)庫按照日期分為365個(gè)文件建立回歸模型增加點(diǎn)密度由于現(xiàn)有的日輻射值數(shù)據(jù)不能覆蓋東三省（如圖），需要對(duì)現(xiàn)有數(shù)據(jù)建模分析，以增加氣象數(shù)據(jù)各點(diǎn)密度。已有數(shù)據(jù)10個(gè)太陽輻射站點(diǎn)，為了實(shí)現(xiàn)回歸模型更好擬合效果，將10個(gè)樣本全部作為回歸參數(shù)。利用SPSS軟件建模步驟：打開一月中旬?dāng)?shù)據(jù)選擇分析→回歸→非線性回歸將輻射值設(shè)為因變量將經(jīng)度（X）和緯度（Y）作為自變量，采用二次趨勢(shì)面模型（f=b0+b1*x+b2*y+b3*x2+b4*x*y+b5*y2）進(jìn)行回歸，回歸方法采用強(qiáng)迫引入法。如圖，在模型表達(dá)式中輸入模型方程。在參數(shù)中設(shè)置參數(shù)初始值在右上方”保存”中選擇殘差在選項(xiàng)中選擇序列二次編程，迭代次數(shù)為50。確定，在文本輸出對(duì)話框中可以顯示迭代次數(shù)和誤差值，以下圖為例確定參數(shù)為。精度R2為0.844二次趨勢(shì)面模型的擬合精度在0.5以上就認(rèn)為可以進(jìn)行建模模擬，本次模型精度為0.844，認(rèn)為可以對(duì)一月份其他站點(diǎn)進(jìn)行模擬，得到所有站點(diǎn)的日輻射值。利用同樣的方法對(duì)每個(gè)月進(jìn)行一次建模估計(jì)，得到東三省加內(nèi)蒙古地區(qū)101個(gè)點(diǎn)的輻射值，再ArcGIS中進(jìn)行插值分析。插值（最高低溫、最低氣溫、降水、濕度、建模后的輻射亮度）將每日的數(shù)據(jù)按照X、Y字段生成點(diǎn)文件使用ArcToolbox樹中的3DAnalystTools工具→RasterInterpolation→Kriging設(shè)置插值數(shù)據(jù)輸出路徑，選擇插值字段，設(shè)置分辨率（1公里與遙感數(shù)據(jù)結(jié)合）批量處理可使用Kriging中的BatchGrid生成插值數(shù)據(jù)為方便ENVI打開令存為TIFF格式裁切使用ENVI軟件同時(shí)打開要裁剪的TIFF圖像和東北三省邊界將矢量邊界設(shè)為ROI感興趣區(qū)域：File→ExportLayerstoROI裁剪T