「2007-11-29:色分けと統計情報の表示」の編集履歴（バックアップ）一覧はこちら
追加された行は緑色になります。
削除された行は赤色になります。
*0.イントロ
前回の授業で行ったアスペクトの再分類したデータの色分けを行う。
また再分類した結果の統計情報を表示の仕方を説明する。
*1.色分け
・マップセットを開く
>位置 ：mie
>Mapset：ユーザ名
・aspect_newsを開く
add a grass layerからaspect_newsを選択し画像データを開く(白黒になっているはず)
・色分けするデータを取り込む
授業のwikiにaspect_news.txtがアップロードされています。
[[このページ>2007-11-29]]の下のほうにaspect_news.txtがあります。
このテキストデータをC:/GRASS/Mie/(mapset)/colrへコピーします。
テキストファイルになっているので、名前をaspect_newsに変更します。
そしてGISの地図画面をクリックすると色が付きます。
※aspect_news.txtの中身と説明
% 4 categories →4つに分類
0:0:255:255 →No aspectはシアンブルー
1:255:255:0 →北は黄色
2:0:0:255 →西は濃い青色
3:0:255:0 →南は緑色
4:255:0:0 →東は赤色
*2.統計情報の表示
**統計情報
再分類した結果、あるカテゴリー(例えばある地域の風向きや植生、降水量の違い)
がどのくらいの割合を占めたかを表示する。
今回は東西南北の傾斜方向の割合を表示します。
**統計情報の表示(r.report)の使い方
プラグイン→GRASS→openGRASStoolsを開く
一番上にあるGRASS SHELLを開く
r.reportの使い方をみる
>r.report help
色々描いてありますが、units(単位)のカテゴリーで
mi(マイル),me(メートル),k,a(エーカー),h(ヘクタール),pと書いてありますが、今回用いるのはk(キロメートル)とp(パーセント)です。
なので4つの傾斜方向の面積の割合を、平方キロメートルとパーセンテージで表示します。
実際に統計情報を表示してみましょう
>r.report map=aspect_news u=k,p
u=k,pが単位を決めるコマンドです。
**ヘッダーを消す
実行日時などの情報を表示せず、求めたい割合だけを表示したい場合は-hをつけます。
>r.report map=aspect_news u=k,p -h
するとヘッダーがなくなります。
**陸の部分だけを考える(MASKをかける)
前回の授業と同じようにマスクをコピーします
>g.copy rast=maskland,MASK
地図画面をクリックすると海の部分が白くなります。
もう一度統計情報を表示します
>r.report map=aspect_news u=k,p -h
するとNo aspectの割合は減少し、no data(海の部分)ができています。
**マスクのかかっている部分を統計対象から外す
今までの情報はNo aspect(海)の割合が存在していたが、マスクのかかっている部分を
統計対象から外すことで陸域だけを考えます。
マスクのかかっている部分を統計対象から外したい場合は-nをつけます。
>r.report map=aspect_news u=k,p -hn
するとno dataが消えます。
**統計データの保存方法
統計データを保存したい場合、まずどこにデータを保存したいかを考えます。
今回はC:/GRASSにAnalysisというフォルダを作ります。
・GRASS SHELLでのフォルダの作り方
C:/GRASSに場所を移動します
>cd C:GRASS
ここ(C:/GRASS)でAnalysisのフォルダを作ります
>mkdir Analysis
これでフォルダを作ることができました。
・統計データを保存する
さっき作ったフォルダの中(Analysis)に統計データをaspect_newsという名前で保存します。
Analysisのフォルダの中に入ります。
>cd Analysis
統計データ(r.report map=aspect_news u=k,p -hn)をaspect_newsの名前で
Analysisに保存します。
>r.report map=aspect_news u=k,p -hn > aspect_news
lsで確認してみましょう
>ls
aspect_newsができています。
保存したデータの表を表示してみましょう
>less aspect_news
表が表示されます。qで終わります。
本日の内容は以上です
*0.イントロ
前回の授業で行ったアスペクトの再分類したデータの色分けを行う。
また再分類した結果の統計情報を表示の仕方を説明する。
*1.色分け
・マップセットを開く
>位置 ：mie
>Mapset：ユーザ名
・aspect_newsを開く
add a grass layerからaspect_newsを選択し画像データを開く(白黒になっているはず)
・色分けするデータを取り込む
授業のwikiにaspect_news.txtがアップロードされています。
[[このページ>2007-11-29]]の下のほうにaspect_news.txtがあります。
このテキストデータをC:/GRASS/Mie/(mapset)/colrへコピーします。
テキストファイルになっているので、名前をaspect_newsに変更します。
そしてGISの地図画面をクリックすると色が付きます。
※aspect_news.txtの中身と説明
% 4 categories →4つに分類
0:0:255:255 →No aspectはシアンブルー
1:255:255:0 →北は黄色
2:0:0:255 →西は濃い青色
3:0:255:0 →南は緑色
4:255:0:0 →東は赤色
*2.統計情報の表示
**統計情報
再分類した結果、あるカテゴリー(例えばある地域の風向きや植生、降水量の違い)
がどのくらいの割合を占めたかを表示する。
今回は東西南北の傾斜方向の割合を表示します。
**統計情報の表示(r.report)の使い方
プラグイン→GRASS→openGRASStoolsを開く
一番上にあるGRASS SHELLを開く
r.reportの使い方をみる
>r.report help
色々描いてありますが、units(単位)のカテゴリーで
mi(マイル),me(メートル),k,a(エーカー),h(ヘクタール),pと書いてありますが、今回用いるのはk(キロメートル)とp(パーセント)です。
なので4つの傾斜方向の面積の割合を、平方キロメートルとパーセンテージで表示します。
実際に統計情報を表示してみましょう
>r.report map=aspect_news u=k,p
u=k,pが単位を決めるコマンドです。
**ヘッダーを消す
実行日時などの情報を表示せず、求めたい割合だけを表示したい場合は-hをつけます。
>r.report map=aspect_news u=k,p -h
するとヘッダーがなくなります。
**陸の部分だけを考える(MASKをかける)
前回の授業と同じようにマスクをコピーします
>g.copy rast=maskland,MASK
地図画面をクリックすると海の部分が白くなります。
もう一度統計情報を表示します
>r.report map=aspect_news u=k,p -h
するとNo aspectの割合は減少し、no data(海の部分)ができています。
**マスクのかかっている部分を統計対象から外す
今までの情報はNo aspect(海)の割合が存在していたが、マスクのかかっている部分を
統計対象から外すことで陸域だけを考えます。
マスクのかかっている部分を統計対象から外したい場合は-nをつけます。
>r.report map=aspect_news u=k,p -hn
するとno dataが消えます。
**統計データの保存方法
統計データを保存したい場合、まずどこにデータを保存したいかを考えます。
今回はC:/GRASSにAnalysisというフォルダを作ります。
・GRASS SHELLでのフォルダの作り方
C:/GRASSに場所を移動します
>cd C:GRASS
ここ(C:/GRASS)でAnalysisのフォルダを作ります
>mkdir Analysis
これでフォルダを作ることができました。
・統計データを保存する
さっき作ったフォルダの中(Analysis)に統計データをaspect_newsという名前で保存します。
Analysisのフォルダの中に入ります。
>cd Analysis
統計データ(r.report map=aspect_news u=k,p -hn)をaspect_newsの名前で
Analysisに保存します。
>r.report map=aspect_news u=k,p -hn > aspect_news
lsで確認してみましょう
>ls
aspect_newsができています。
保存したデータの表を表示してみましょう
>less aspect_news
表が表示されます。qで終わります。
本日の内容は以上です
