2013年10月28日 星期一

散佈圖、相關分析、迴歸分析

地理學計量分析中的科學化推論,建立因果關係的模式。

三個重要主題:散佈圖、相關分析、迴歸分析
自耕農數 耕地面積(公頃)
屏東市 3661 2777
潮州鎮 1531 2570
東港鎮 1375 949
萬丹鄉 3012 2754
長治鄉 2965 3653
麟洛鄉 1281 1020
九如鄉 1352 2657
里港鄉 1528 2925
鹽埔鄉 1707 3881
高樹鄉 4671 7506
萬巒鄉 2343 4577
內埔鄉 3712 5390
竹田鄉 2320 2352
新埤鄉 1306 4724
枋寮鄉 1832 3697
新園鄉 3249 2001
崁頂鄉 1476 1932
林邊鄉 1312 593
南州鄉 854 1577
佳冬鄉 1553 1449


以上的相關係數是0.64
相關係數的公式是這樣的

現在用電腦計算非常方便
只用兩個程序就可以了



沒有電腦那個年代,就要用一張表格一個一個數字計算,如果不怕驚嚇,可以參考一下.






2013年10月22日 星期二

社會科科展的題目參考


X年X班學生的通勤圈分析
台電電箱的編號解碼
花蓮縣的外配分布
7-11 VS 全家,以花蓮市區分布為例
找尋花蓮市區舊東線鐵路的足跡
由花蓮2010年風花圖分析風力發電的可行性
利用DEM分析吉安鄉的潛勢土石流研究
哈!我家住在斷層上
花蓮縣原住民的分布
花蓮縣農產分布
花蓮縣日照時數分析,太陽光電的可能性
花蓮縣等雨線的分布
國小班級及學區人口的關係

班級學生分布


按此連結填寫表單

https://docs.google.com/forms/d/1UBsjS8UPQar9HXVq53pslB1i3gdO1cM1lKggd1wTvMU/viewform

2013年10月19日 星期六

GIS練習場--24--等值線的繪製

等值線的繪製,以香港的等雨量線繪製為例

這次以香港天文台及另一個提供各國SHP圖檔的網站資料來製作
連上香港天文台的網站,把資料下載
http://www.weather.gov.hk/cis/data/drf_summary_c.htm
網站上有各個雨量測站的經緯度坐標


點進各個測站,有2000年到2013年的雨量資料,就選一個2010年


把每一個測站的資料都下載,分4個欄位,x,y,name,
整理成EXCEL,要存CSV的格式。HONGKONG_STATION



http://www.gadm.org/country
到這個網站把SHAPEFILE下載回來


把圖讀進來


要用讀CSV的功能選單,【圖層】,【新增】


把剛才整理好的CSV讀進來。


讀進來了。

要另存為shapefile,目的是要把轉成向量圖檔,才可以計算。


轉成向量圖檔時會要求坐標,一定要跟原本下載的圖相同,原本是WGS84


新的向量圖檔,就用hongkongrain。


經過一最後的確認,存檔。


完成匯出向量圖,存檔成功了,


貼上標籤,看一下各個測站的站名。


看一下屬性表,4個欄位,坐標,站名及雨量紀錄。


以hongkongrain為工作圖層,點選【影像工具】,【分析】。【grid(interpolation)】,interpolation是【內插】。


把【interpolation】選單填妥, z-field指定為mm,。輸出檔案取名為rain_interpolation。


完成後,把raster圖也讀進來。

新產生的raster讀進來,顏色深淺,表示雨量值的高低。


再利用raster圖層,輸出一個等值線的shp檔,點選【影像工具】,【extraction】,【等高線】。


把等高線選單一一填妥。等高線以100為單位,輸出檔名為rain100。


等值線完成了


把各個圖層調整一下,完成了一張等值線圖了。



2013年10月18日 星期五

GIS練習場--23--Qgis結合GRASS外掛程式以輸出3D效能

Qgis結合GRASS外掛程式以輸出3D效能

把日月潭所在的魚池鄉圖層讀入,分別為向量圖層及網格式的DEM圖層。


選取,【附加元件】【GRASS】【新增地圖集】


系統要求產生一個新的地圖集,所謂地圖集,就是一個新的資料夾,本例子把地圖集放在現存資料夾中,【yuchih】是現存資料夾。原本的魚池鄉的向量圖層及DEM圖層也在裡面。

GRASS位置:表示在現存資料夾再開一個資料夾,就取名為yuchihnew。


到【yuchih】資料夾中,可以看到增加了一個叫yuchihnew的資料夾,而GRASS把它叫地圖集。


投影的坐標系統,要設定相同的坐標,否則會位置有誤差。


出現圖層的坐標,再確認。

地圖集,表示由GRASS產生的新圖層會放進地圖集中,就取名【yuchihnew】


表示新建立的地圖集是放在下列的地方。


打開資料夾,看到有三層的資料夾,剛才所做的設定,就是要建立這些層狀目錄。





最後,終於完成了。



自動關閉視窗。




選取【附加元件】【GRASS】【開啟GRASS工具】


選【樹狀模組】【檔案管理】


【檔案管理】以下有三個選項,【IMPORT INTO GRASS】


再選【Import  raster into GRASS 】


選【Import raster into GRASS from QGIS view】


出現【SELECT  A  LAYER】


要輸出一個圖層,意思是說,把yuchihdem圖層再輸出成GRASS專用的圖層。所以說明文字是:NAME FOR OUTPUT RASTER MAP

選取【yuchihdem】


點選[運行]後輸出完成,SUCCESSFULLY  FINISHED。



 CLOSE

 重複上一個動作,現在要INPUT VECTOR 圖層




有數個選項第一個選項是[由QGIS中讀入向量圖層]

選一個圖層來INPUT



就選YUCHIHA


輸出圖層的名字為VECTORNEW



點選[運行],輸出完成了,關掉視窗。



現在需要把剛才產生的兩個圖層讀進來



選取,【附加元件】【GRASS】【新增GRASS向量圖層集】
向量圖就是剛才輸出完成的VECTORNEW


到YUCHIH的YUCHIHNEW中去找




找到了VECTORNEW

讀進來了,在圖層視窗中出現了VECTORNEW


 再讀取網格資料,選取,【附加元件】【GRASS】【新增GRASS影像圖層】

 

選DEMNEW



兩個圖層都讀進來了,圖層視窗中出現了VECTORNEW及DEMNEW


【樹狀模組】選單中有一個【3D視覺化】工具。





 在【3D視覺化】工具中,出現取RASTER圖層的選單


RASTER圖層選取DEMNEW資料。






 再選向量資料圖層,選取VECTORNEW




出現3D的控制視窗【NVIZ】


利用左側的控制按鍵,可以旋轉圖層,行政區的界線是向量圖層.


也可以輸出一張TIF格式的圖