空間插值分析課件

實(shí)習(xí)七空間插值分析1.用AvenueScript程序腳本進(jìn)行趨勢面分析2.核密度估算3.用IDW進(jìn)行空間插值4.兩種薄板樣條函數(shù)法的比較5.用AvenueScript程序腳本進(jìn)行普通克里金插值分析1實(shí)習(xí)七空間插值分析1.用AvenueScript程序腳本今天請大家完成007和008的內(nèi)容(全部內(nèi)容)，另外007的內(nèi)容進(jìn)行了修改完善，請大家更新數(shù)據(jù)，如果有全部完成了的同學(xué)請告知我；下周我們進(jìn)行上機(jī)實(shí)踐考試（共兩道題）；請大家運(yùn)行“飛鴿傳書”軟件。2今天請大家完成007和008的內(nèi)容(全部內(nèi)容)，另外007的1.用AvenueScript程序腳本進(jìn)行趨勢面分析所需數(shù)據(jù)：stations.shp，包含愛達(dá)荷州105個氣象站數(shù)據(jù)的SHAPE文件；idoutl.shp，一個顯示愛達(dá)荷州外輪廓的SHAPE文件；trend.ave，進(jìn)行趨勢面分析用的AvenueScript程序腳本。該習(xí)作要求調(diào)用AvenueScript程序腳本，使用stations.shp文件中名為ann_prec的屬性表，運(yùn)行趨勢面分析。ann_prec是自1961－1990年的年平均降水量。31.用AvenueScript程序腳本進(jìn)行趨勢面分析所需數(shù)啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。4啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。4打開一個新視圖，將stations.shp和idoutl.shp加到視圖中，并從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。5打開一個新視圖，將stations.shp和idoutl.s6677在Project窗口中點(diǎn)擊Script和new打開Script1；8在Project窗口中點(diǎn)擊Script和new打開Scrip在Script菜單中點(diǎn)擊按鈕，瀏覽到trend.ave的路徑并雙擊它。9在Script菜單中點(diǎn)擊按鈕，瀏覽到trend.ave的路徑在按鈕欄中點(diǎn)擊編譯按鈕對該腳本程序進(jìn)行編譯。10在按鈕欄中點(diǎn)擊編譯按鈕對該腳本程序進(jìn)行編譯。10回到視圖界面，雙擊工具或按鈕欄空白處，出現(xiàn)對話框，進(jìn)行相關(guān)設(shè)置，并返回到視圖。11回到視圖界面，雙擊工具或按鈕欄空白處，出現(xiàn)對話框，進(jìn)行相關(guān)設(shè)按照以下步驟進(jìn)行：12按照以下步驟進(jìn)行：12Grid1是基于三階多項(xiàng)式方程輸出的趨勢面分析網(wǎng)格13Grid1是基于三階多項(xiàng)式方程輸出的趨勢面分析網(wǎng)格13在趨勢面分析網(wǎng)格的基礎(chǔ)上繪制等值線。14在趨勢面分析網(wǎng)格的基礎(chǔ)上繪制等值線。14按照設(shè)置可自動生成等值線專題ContoursofGrid1，將其激活并顯示出來，并可進(jìn)行自動標(biāo)注。結(jié)果圖如下：15按照設(shè)置可自動生成等值線專題ContoursofGrid2.核密度估算所需數(shù)據(jù)：deer.shp，一個顯示鹿點(diǎn)的點(diǎn)狀SHAPE文件。該習(xí)作用核估算法依據(jù)deer.shp計算每公頃看到的鹿的平均數(shù)目。鹿點(diǎn)的數(shù)據(jù)具有50m的最小分辨距離，因此，有些點(diǎn)多次看到。162.核密度估算所需數(shù)據(jù)：deer.shp，一個顯示鹿點(diǎn)的點(diǎn)狀啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。17啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。17打開新視圖并把deer.shp添加到視圖中，并設(shè)置地圖單位為米。18打開新視圖并把deer.shp添加到視圖中，并設(shè)置地圖單位為1919按照圖解進(jìn)行相關(guān)操作和設(shè)置：20按照圖解進(jìn)行相關(guān)操作和設(shè)置：202121從分析菜單中選取計算密度菜單項(xiàng)，并進(jìn)行相關(guān)設(shè)置：22從分析菜單中選取計算密度菜單項(xiàng)，并進(jìn)行相關(guān)設(shè)置：222323生成的成果圖如下，為了使該圖面更簡潔，可進(jìn)行后頁的設(shè)置。24生成的成果圖如下，為了使該圖面更簡潔，可進(jìn)行后頁的設(shè)置。24252526263.用IDW進(jìn)行空間插值所需數(shù)據(jù)：stations.shp，包含愛達(dá)荷州105個氣象站數(shù)據(jù)的SHAPE文件；idoutl.shp，一個顯示愛達(dá)荷州外輪廓的SHAPE文件；本習(xí)作包括三個步驟：首先將idoutl.shp轉(zhuǎn)化成格網(wǎng)idoutlgd；其次用反距離權(quán)重法（IDW）并以idoutlgd作為屏蔽掩膜格網(wǎng)創(chuàng)建降水量格網(wǎng)；第三從降水量格網(wǎng)生成等值線圖

273.用IDW進(jìn)行空間插值所需數(shù)據(jù)：stations.shp啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。28啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。28打開一個新視圖，將stations.shp和idoutl.shp加到視圖中，并從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。29打開一個新視圖，將stations.shp和idoutl.s3030激活專題并打開其屬性表，屬性表中的ann_prec字段是用于空間插值的Z值。31激活專題并打開其屬性表，屬性表中的ann_prec字段是用于3232333334343535363637373838393940404141424243434.兩種薄板樣條函數(shù)法的比較所需數(shù)據(jù)：習(xí)作3的stations.shp和格網(wǎng)idoutlgd。該習(xí)作對兩種薄板樣條函數(shù)產(chǎn)生的結(jié)果作比較。本習(xí)作有三個部分：①用規(guī)則樣條法創(chuàng)建插值格網(wǎng)；②用薄板張力樣條法創(chuàng)建插值格網(wǎng)；③用局部運(yùn)算比較這兩種格網(wǎng)，從中可以看出兩種插值方法產(chǎn)生的不同結(jié)果。444.兩種薄板樣條函數(shù)法的比較所需數(shù)據(jù)：習(xí)作3的station啟動系統(tǒng)，加載空間分析擴(kuò)展模塊。45啟動系統(tǒng)，加載空間分析擴(kuò)展模塊。45打開一個新視圖，并把stations.shp和idoutlgd加到視圖中。從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。46打開一個新視圖，并把stations.shp和idoutlg4747484849495050運(yùn)行規(guī)則薄板樣條插值。51運(yùn)行規(guī)則薄板樣條插值。51系統(tǒng)生成格網(wǎng)，然后修改生成的專題名為Regularized。52系統(tǒng)生成格網(wǎng)，然后修改生成的專題名為Regularized。5353用張力薄板樣條法創(chuàng)建插值格網(wǎng)。并按照前述方法查看格網(wǎng)并重命名為Tension。54用張力薄板樣條法創(chuàng)建插值格網(wǎng)。并按照前述方法查看格網(wǎng)并重命名比較兩種方式創(chuàng)建的格網(wǎng)。55比較兩種方式創(chuàng)建的格網(wǎng)。555656575758585.用AvenueScript程序腳本進(jìn)行普通克里金插值分析所需數(shù)據(jù)：習(xí)作3的stations.shp和格網(wǎng)idoutlgd；kriging.ave，運(yùn)行克里金法的AvenueScript程序腳本。該習(xí)作要求用線型模型的kriging.ave運(yùn)行普通克里金法。595.用AvenueScript程序腳本進(jìn)行普通克里金插值分啟動系統(tǒng)，加載空間分析擴(kuò)展模塊。60啟動系統(tǒng)，加載空間分析擴(kuò)展模塊。60打開一個新視圖，并把stations.shp和idoutlgd加到視圖中。從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。61打開一個新視圖，并把stations.shp和idoutlg6262在Project窗口中點(diǎn)擊Script和new打開Script1；63在Project窗口中點(diǎn)擊Script和new打開Scrip在Script菜單中點(diǎn)擊按鈕，瀏覽到kriging.ave的路徑并雙擊它。64在Script菜單中點(diǎn)擊按鈕，瀏覽到kriging.ave的在按鈕欄中點(diǎn)擊編譯按鈕對該腳本程序進(jìn)行編譯。65在按鈕欄中點(diǎn)擊編譯按鈕對該腳本程序進(jìn)行編譯。65回到視圖界面，雙擊工具或按鈕欄空白處，出現(xiàn)對話框，進(jìn)行相關(guān)設(shè)置，并返回到視圖。66回到視圖界面，雙擊工具或按鈕欄空白處，出現(xiàn)對話框，進(jìn)行相關(guān)設(shè)按照以下步驟進(jìn)行：67按照以下步驟進(jìn)行：67系統(tǒng)自動生成兩個格網(wǎng)：Vargrid和grid2（如下圖），其中Vargrid為估算的離差格網(wǎng)，以下介紹如何將其改變?yōu)楣浪愕臉?biāo)準(zhǔn)差格網(wǎng)。注意：如果在這一步?jīng)]有自動生成的話，請檢查D:\PAGIS\DATA\07路徑中是否有vargrid和info文件夾，將這兩個文件夾刪除即可。68系統(tǒng)自動生成兩個格網(wǎng)：Vargrid和grid2（如下圖），696970707171TheEnd72TheEnd72實(shí)習(xí)七空間插值分析1.用AvenueScript程序腳本進(jìn)行趨勢面分析2.核密度估算3.用IDW進(jìn)行空間插值4.兩種薄板樣條函數(shù)法的比較5.用AvenueScript程序腳本進(jìn)行普通克里金插值分析73實(shí)習(xí)七空間插值分析1.用AvenueScript程序腳本今天請大家完成007和008的內(nèi)容(全部內(nèi)容)，另外007的內(nèi)容進(jìn)行了修改完善，請大家更新數(shù)據(jù)，如果有全部完成了的同學(xué)請告知我；下周我們進(jìn)行上機(jī)實(shí)踐考試（共兩道題）；請大家運(yùn)行“飛鴿傳書”軟件。74今天請大家完成007和008的內(nèi)容(全部內(nèi)容)，另外007的1.用AvenueScript程序腳本進(jìn)行趨勢面分析所需數(shù)據(jù)：stations.shp，包含愛達(dá)荷州105個氣象站數(shù)據(jù)的SHAPE文件；idoutl.shp，一個顯示愛達(dá)荷州外輪廓的SHAPE文件；trend.ave，進(jìn)行趨勢面分析用的AvenueScript程序腳本。該習(xí)作要求調(diào)用AvenueScript程序腳本，使用stations.shp文件中名為ann_prec的屬性表，運(yùn)行趨勢面分析。ann_prec是自1961－1990年的年平均降水量。751.用AvenueScript程序腳本進(jìn)行趨勢面分析所需數(shù)啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。76啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。4打開一個新視圖，將stations.shp和idoutl.shp加到視圖中，并從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。77打開一個新視圖，將stations.shp和idoutl.s786797在Project窗口中點(diǎn)擊Script和new打開Script1；80在Project窗口中點(diǎn)擊Script和new打開Scrip在Script菜單中點(diǎn)擊按鈕，瀏覽到trend.ave的路徑并雙擊它。81在Script菜單中點(diǎn)擊按鈕，瀏覽到trend.ave的路徑在按鈕欄中點(diǎn)擊編譯按鈕對該腳本程序進(jìn)行編譯。82在按鈕欄中點(diǎn)擊編譯按鈕對該腳本程序進(jìn)行編譯。10回到視圖界面，雙擊工具或按鈕欄空白處，出現(xiàn)對話框，進(jìn)行相關(guān)設(shè)置，并返回到視圖。83回到視圖界面，雙擊工具或按鈕欄空白處，出現(xiàn)對話框，進(jìn)行相關(guān)設(shè)按照以下步驟進(jìn)行：84按照以下步驟進(jìn)行：12Grid1是基于三階多項(xiàng)式方程輸出的趨勢面分析網(wǎng)格85Grid1是基于三階多項(xiàng)式方程輸出的趨勢面分析網(wǎng)格13在趨勢面分析網(wǎng)格的基礎(chǔ)上繪制等值線。86在趨勢面分析網(wǎng)格的基礎(chǔ)上繪制等值線。14按照設(shè)置可自動生成等值線專題ContoursofGrid1，將其激活并顯示出來，并可進(jìn)行自動標(biāo)注。結(jié)果圖如下：87按照設(shè)置可自動生成等值線專題ContoursofGrid2.核密度估算所需數(shù)據(jù)：deer.shp，一個顯示鹿點(diǎn)的點(diǎn)狀SHAPE文件。該習(xí)作用核估算法依據(jù)deer.shp計算每公頃看到的鹿的平均數(shù)目。鹿點(diǎn)的數(shù)據(jù)具有50m的最小分辨距離，因此，有些點(diǎn)多次看到。882.核密度估算所需數(shù)據(jù)：deer.shp，一個顯示鹿點(diǎn)的點(diǎn)狀啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。89啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。17打開新視圖并把deer.shp添加到視圖中，并設(shè)置地圖單位為米。90打開新視圖并把deer.shp添加到視圖中，并設(shè)置地圖單位為9119按照圖解進(jìn)行相關(guān)操作和設(shè)置：92按照圖解進(jìn)行相關(guān)操作和設(shè)置：209321從分析菜單中選取計算密度菜單項(xiàng)，并進(jìn)行相關(guān)設(shè)置：94從分析菜單中選取計算密度菜單項(xiàng)，并進(jìn)行相關(guān)設(shè)置：229523生成的成果圖如下，為了使該圖面更簡潔，可進(jìn)行后頁的設(shè)置。96生成的成果圖如下，為了使該圖面更簡潔，可進(jìn)行后頁的設(shè)置。24972598263.用IDW進(jìn)行空間插值所需數(shù)據(jù)：stations.shp，包含愛達(dá)荷州105個氣象站數(shù)據(jù)的SHAPE文件；idoutl.shp，一個顯示愛達(dá)荷州外輪廓的SHAPE文件；本習(xí)作包括三個步驟：首先將idoutl.shp轉(zhuǎn)化成格網(wǎng)idoutlgd；其次用反距離權(quán)重法（IDW）并以idoutlgd作為屏蔽掩膜格網(wǎng)創(chuàng)建降水量格網(wǎng)；第三從降水量格網(wǎng)生成等值線圖

993.用IDW進(jìn)行空間插值所需數(shù)據(jù)：stations.shp啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。100啟動系統(tǒng)并加載空間分析擴(kuò)展模塊。28打開一個新視圖，將stations.shp和idoutl.shp加到視圖中，并從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。101打開一個新視圖，將stations.shp和idoutl.s10230激活專題并打開其屬性表，屬性表中的ann_prec字段是用于空間插值的Z值。103激活專題并打開其屬性表，屬性表中的ann_prec字段是用于1043210533106341073510836109371103811139112401134111442115434.兩種薄板樣條函數(shù)法的比較所需數(shù)據(jù)：習(xí)作3的stations.shp和格網(wǎng)idoutlgd。該習(xí)作對兩種薄板樣條函數(shù)產(chǎn)生的結(jié)果作比較。本習(xí)作有三個部分：①用規(guī)則樣條法創(chuàng)建插值格網(wǎng)；②用薄板張力樣條法創(chuàng)建插值格網(wǎng)；③用局部運(yùn)算比較這兩種格網(wǎng)，從中可以看出兩種插值方法產(chǎn)生的不同結(jié)果。1164.兩種薄板樣條函數(shù)法的比較所需數(shù)據(jù)：習(xí)作3的station啟動系統(tǒng)，加載空間分析擴(kuò)展模塊。117啟動系統(tǒng)，加載空間分析擴(kuò)展模塊。45打開一個新視圖，并把stations.shp和idoutlgd加到視圖中。從視圖下拉菜單中選擇屬性，將地圖單位設(shè)置為米。118打開一個新視圖，并把stations.shp和idoutlg11947120481214912250運(yùn)行規(guī)則薄板樣條插值。123運(yùn)行規(guī)則薄板樣條插值。51系統(tǒng)生成格網(wǎng)，然后修改生成的專題名為Regularized。124系統(tǒng)生成格網(wǎng)，然后修改生成的專題名為Regularized。12553用張力薄板樣條法創(chuàng)建插值格網(wǎng)。并按照前述方法查看格網(wǎng)并重命名為Tension。126用張力薄板樣條法創(chuàng)建插值格網(wǎng)。并按照前述方法查看格網(wǎng)并重命名比較兩種方式創(chuàng)建的格網(wǎng)。127比較兩種方式創(chuàng)建的格網(wǎng)。551285612957130585.用AvenueScript程序腳本進(jìn)行普通克里金插值分析所需數(shù)據(jù)：習(xí)作3的stations.shp和格網(wǎng)idoutlgd；kriging.ave，運(yùn)行克里金法的AvenueScript程序腳本。該習(xí)作要求用線型模型的kriging.ave運(yùn)行普通