登山社出隊地圖

2001~2009我在登山社出隊地圖

在較大的地圖上查看Mountains
註：綠色樹木為勘查或郊山；藍色波浪為溯溪隊；黃色太陽為高山傳統隊伍；紅色山崩屬於技術攀岩隊伍。

我發現爬山主要集中在宜蘭和烏來一帶，這也反應出操、鳥、賤的登山風格。而中、南部的隊伍數不到十隊，但幾乎都是長程的隊伍。

雖然我去過的隊伍已勘查居多，但一個不會游泳的人居然去過那麼多溯溪隊!!。

這幾年，登山社爬山的記憶在似乎退化的特別快。今天怪叔叔和子容從美國回台，從學術研究聊到美國生活，再聊到爬山。問到我過去爬過哪些台灣中部的山，我居然完全想不來。如果不常整理，可能再過幾年就會忘光了。

ArcGIS將Raster網格資料轉成xyz文字檔

以前在把Raster網格資料轉成xyz的表格，通常經過許多繁複的程序： Raster to point > add xy > open attribute table >export table，這幾天從ESRI網站發現原來只要一個指令就可以做到以上的步驟：

資料來源： http://support.esri.com/en/knowledgebase/techarticles/detail/26568
The following steps require the Spatial Analyst extension to be active and the Spatial Analyst toolbar to be turned on.

1. Add the grid to ArcMap.
2. Select the Raster Calculator option from the Spatial Analyst menu.
3. Use the following expression in the Raster Calculator where [ingrid] is the name of the grid layer:
   
   

   outfile.txt = SAMPLE([ingrid],[ingrid])

4. The text file called outfile.txt is written to the working directory. The working directory can be determined or set from the Spatial Analyst > Options > General tab. 
   
   
    The format of the output text file may have to be modified to be compatible with specific needs. The first and last columns in the text file are the z-values. 
   
   

    ESRI does not provide the VBA code for the above procedure. The necessary code must be developed by the user.

   
   

   翻譯

   1. 將你的grid加入ArcMap.
   2. 從 Spatial Analyst 開啟   Raster Calculator .
   3. 輸入以下的code，其中ingrid為你要轉出的檔案，例如 abc.txt = SAMPLE(C:\a1,  C:\a1 )。
      
      

      outfile.txt = SAMPLE([ingrid],[ingrid])

   4. 如果你要改輸出路徑，可以從 Spatial Analyst > Options > General tab去修改，但是其實只要在輸出檔名增加路徑即可，例如: D:\abc.txt = SAMPLE(C:\a1,  C:\a1 )
   5. 
      

如何批次製作河道橫斷面

 
在模擬水理模式時，常需要由DEM資料收產橫斷面資料，如何自動生產垂直縱斷面的polyline常需要花費許多人力，以下介紹當有一個河道縱斷面的polyline，並且已經切割為數個河段(feature)，如何生批次產出河道橫斷面。
1. 到ESRI網站下載Create Perpendicular Lines to Each Segment of a Shapefile

http://arcscripts.esri.com/details.asp?dbid=15756

ArcGIS ArcMap建立徐昇氏多邊形

如何利用ArcGIS產生徐昇氏多邊形?

Google Earth 6.0版推出了

這次新版的Google Earth整合Google Map的街景服務，在導覽工具上出現街景人形圖示，操作方法跟Google Map一樣把圖案拖進即可。

Google Earth 6.0 新功能

與Google Map不同的是，在Google Earth瀏覽街景更為順暢、也更直覺。

新版本個另一項賣點就是3D樹景，但是目前似乎僅於在東京(Tokyo)等大都市郊區有。

3D樹景

當然新版本還加強其他原有的功能，詳情見 Google Earth http://www.google.com/intl/zh-TW/earth/index.html
Google Earth 6.0 下載位置  http://www.google.com/intl/zh-TW/earth/download/ge/agree.html

利用ArcGIS Model Builder建立差異化植生指數(NDVI)模組

差異化植生指數(Normalized Difference Vegetation Index；NDVI)是一個常見的環境因子，利用多光譜之衛星影像用來評估地表上的植生覆蓋情形，用來計算植物生長與生產力，另外也常被用來應用在崩塌地判釋或土壤沖蝕量估算或是土地利用變化之監測，原理是利用植生在不同波段之反射率的差異，將此差異標準化，如圖A所示，植生(綠線)在近紅外光段反射率約50%；而在紅光段的反射率卻只有10%左右。

圖 A 土壤、水和植生的典型光譜反射曲線
利用上述之特性，NDVI方程式表為:  NDVI = (NIR – RED) / (NIR + RED) ，其中NIR為光譜近紅外光段；RED為光譜紅光段，計算的值範圍會在 -1 至 1 之間。

統計面資料內polygon的網格

我們在使用 ArcGIS時，常用到需要向量式的shapefile和網格式資料要交互分析，例如:一個集水區的平均坡度或平均高程，由於集水區通常為shapefile，坡度圖通常圍網格資料，因此一個集水區通常會包含許多網格資料，要如何去計算最大值、最小值、加總或平均值，這時就需要用到zonal stastics，以下就計算坡單元平均坡度為例，來說明如何使用zonal stastics


 
1. 開DEM算slope

2.把坡單元也加進來

 
3.check坡單元欄位中，必須有一個欄位是可以分別每個feature（FID不算），如果沒有，新增一個long integer的欄位，再用field calculator，欄位等於FID即可。

 
4.zonal

 
按OK


 
5.EXPORT

完成


 
開啟 Table欄位中的mean則代表，平均坡度。


 
坡向的算法也是類似，但是坡向比較麻煩，因為1度和359應該相似，但是計算出的結果=(1+359)/2 = 180，導致北向坡與坡向變化大的坡有的錯誤，如果懶的時候也是可以這樣算。正確做法要先算出坡向圖後，把角度換成向量ex:(0,-1)，去做zonal，再用sin與cos去分辨像現，則可以計算出平均坡向180。


 

如何在ArcGIS中設定TWD97與TWD67座標系統(2)

繼上一篇  如何在ArcMAP設定TWD97與TWD67座標 ，因為這個方法實在有點麻煩本身在設定座標，而且實在沒把握到底參數設定的正確性，建議可以採用別的高手已經用設好的座標檔。推薦中研院計算中心的一個網站，裡面有許多很好用的GIS小工具。

1. 開啟 研究支援工具(http://www.ascc.sinica.edu.tw/gis/ISTIS/tools.html)

[caption id="attachment_142" align="alignnone" width="503" caption="研究支援工具"][/caption]

1. 開啟壓縮檔，將檔案解壓縮到 C:\Program File\ArcGIS\Coordinate Systems

[caption id="attachment_144" align="alignnone" width="608" caption="解壓縮至Coordinate Systems"][/caption]

1. 開啟 ArcMAP >> Layers  >> Properties >>  Coordinate Systems >> Predefined就可以看到剛剛課訂好的坐標系系統。
   
   [caption id="attachment_145" align="alignnone" width="638" caption="開啟Properties "][/caption]
   
   [caption id="attachment_148" align="alignnone" width="486" caption="看到剛剛加入的座標檔"][/caption]

如何在ArcMAP設定TWD97與TWD67座標 1

TWD97和TWD67是台灣常用座標系統，在ArcGIS上面如何設定？

如何讓GridView刪除資料時，出現確認對話框

如何讓GridView刪除資料時，出現確認對話框：

網路上有一推範例，找到一種最簡單的方式。

<asp:LinkButton ID="EditButton" runat="server" CausesValidation="False"

CommandName="Edit" Text="編輯" />

&nbsp;<asp:LinkButton ID="DeleteButton" runat="server" CausesValidation="False"

CommandName="Delete" Text="刪除" OnClientClick="if(!window.confirm('確定要刪除此項目嗎？'))return false;" />

&nbsp;<asp:LinkButton ID="NewButton" runat="server" CausesValidation="False"

CommandName="New" Text="新增" />

MODIES看台灣

MODIES是NASA科學衛星，免費提供250m、500m、1000m解析度之影像

CAD to shapefile

試過許多ArtoCAD轉shapefile的轉檔方式，發現一個好用又強力的工具Extract CAD Layers to shapefiles，可以轉point, polyline, polygon，包括Annotation也可以轉換（之前遇到問題就是 CAD Annotation讀不出來），並且保留完整的attribute，下載網址：http://arcscripts.esri.com/details.asp?dbid=14679

ArcSWAT安裝好了

今天終於把ArcSWAT裝好了！

準備把他拿來跑降雪分布看看

不過還要做好多初始設定

而且一直有ERROR產生

不知道為什麼GIS的軟體好像一直都不太穩定

像是Desktop的Model Builder不可解的問題一堆

用Model Builder來做真的會把人倒到瘋掉