Subject: [ts-mag:00022] OpenGIS Mail Magazine for T.S.U 【第 227号】
Date: Fri, 25 Nov 2005 23:28:06 +0900
From: WAZA Toshihiko
OpenGIS Mail Magazine for Technical Support Users 2005/11/25
Keywords:OpenGIS,TNTmips,MailMagazine,RasterProperty,Extract,SML
オープンGIS
テクニカル・サポート専用メールマガジン
第227号「 ラ ス タ の 集 計 機 能 チェックしました！ 」
#################################################################
株式会社 オープンGIS
先週ご報告しました「ラスタの集計」機能のチェックを行いましたので、
今週は予定を変更して、その結果についてご報告します。
チェックしたメニューは「解析処理」>「ベクタ」>「属性」>「ラスタの
集計...」です。
テストデータとして、１から７の値(ＤＮ値)を持つ縦横 22 x 22 のラス
タと、円形のポリゴンベクタを用意しました。
ラスタのセルサイズは、100 x 100(m)です。
・・・
「境界のセル」のオプションが少し変わりました。
１．部分的に内側なら含む
２．include if center inside（セルの中心が内側であれば含む）
３．比例配分計算
４．部分的に飛び出していたら除外
２番目は以前「面積の大きい方を含む」でした。今回定義がより正確に
なりました。今回は、各オプションについてもう少し補足説明をしてお
きます。
●１の「部分的に内側なら含む」とは、このようなカウントのことです。
はみ出たセルもカウントしますので、面積は大きめに出ます。
●２の「セルの中心が内側であれば含む」とは、このようなカウントで
す。
セル中心の点と集計する図形(この場合円)の包含関係によって、カウン
トしたりされなかったりします。面積的に大きくなるか小さくなるかは、
図形の境界の形状に依存します。
●３の「比例配分計算」は、
処理が正しく行われれば、セルの面積は常に集計する図形の面積と同じ
にならなくてはいけません。
●４の「部分的に飛び出していたら除外」は、こうです。
セルが集計図形から少しでもはみ出ていたらカウントされません。面積
的には常に小さめに出ます。用途としては、境界付近のデータを採用し
たくない場合です。
これらのオプションを１枚の絵にまとめました。
それぞれの手法で計算されるセルの面積をS1, S2, S3, S4, 集計するポ
リゴン(今の例では円）の面積をSc とすると、一般的に次のような大小
関係が成り立ちます。
S1 > S3 = Sc > S4
S1 > S2 > S4
S2 と S3 の大小は、境界の形状などによって変わるでしょう。
またラスタのセルサイズを横軸にとって、各手法ごとのセルの合計面積
の挙動を予想すると、これも大体次のようになるでしょう。
３の「比例配分」以外の手法において、ラスタのセルサイズが相対的に
大きくなればなるほど、集計する図形との面積との乖離(かいり)が大き
くなります。この辺に気をつけてこれらのオプションを使うようにしま
しょう。
・・・
さて解説はこの位にして、結果を示します。
１）まずは手法ごとに結果を整理しました（ケースＡ)。
集計は、大きい方の円(10と描いたほう)について行いました。
１の「部分的に内側なら含む」、４の「部分的に飛び出していたら除外」
の集計結果(カウント数)は、円の左側の接し方がきわどいですが、目視
確認の結果正しいことが分かります。
２）３の「比例配分計算」については他の方法を使ってさらに計算して
みました。
【第222号】【第223号】で、インターセクトやポリゴンの集計機能を使
って、土地利用図ベクタの集計を行いました。ラスタを一度ベクタに変
換してこの方法で面積を求めてみました。
▼ まずは"ラスベク"変換します。
▼ "ラスベク"変換のパラメータは簡単です。
▼ チェックには「解析処理」>「ベクタ」>「属性」>「ポリゴンの集計
...」を使用しました。
▼ データなど処理に必要なパラメータをセットします。
▼ 表示ウィンドウです。
▽ 実行後の集計結果です。円10の集計結果です。
(黄色の小計の部分のみ抽出しました)
これは、ケースＡの集計結果の３と同じです。
当たり前といったら何ですが、ベクタに変換しても同じ結果が得られま
した。
３）次に２の「セルの中心が内側であれば含む」も他の方法で検証して
みました。セルの中心を求めるために、ラスタをポイントに変換しまし
た。
▼ ラスタをポイントに変換するには次のＳＭＬを使用します。
このＳＭＬは「解析処理」>「ＳＭＬ」>「実行...」で開きます。
▽ 変換後のポイントはこのような状況です。ラスタのＤＮ値はInternal
テーブルのＺ座標に入っています。
▼ 切り出します。「解析処理」>「ベクタ」>「組み合わせ(演算)...」
の「インターセクト(AND)」を使いました。
▽ 処理後のポイントはこのような感じです。円の内側に含まれる点が切
り出されました。
▼ １０番の円に対するポイントを全部ボックスで選択し、Internalテー
ブルを表示して集計します。並べ替えて、点の個数をカウントしました。
(黄色のカウントの部分のみ抽出しました)
ケースＡの２の場合のセルの個数と同じになっています。
４）セルサイズの効果
最後に、ラスタをさらに細かく分割(リサンプル)して、セルサイズを５
０ｍ(ケースＢ), ２０ｍ(ケースＣ)の場合について集計しました。
リサンプルには、セルの値を変えない"最近隣法 Nearest Neighbour"を
使用しました。
セルサイズと合計面積の関係をプロットしました。
予想どうりですが、手法２の変動が小さかったのは、円のようなあまり
極端な凸凹のない図形を集計単位に選んだからだと思います。
・・・
以上、結果は良好のようです。
これまでご不便、ご迷惑をおかけいたしました。
□■Ｔｉｐｓ□■
===============================================================
ラスタの座標系によらずに、機械的にセルを再分割(リサンプル)するに
は、「解析処理」>「ラスタ」>「抜き出し...」メニューの[拡大/方向]
タブを使います。
(メルマガ【第198号】で既出)
================================================================
エラー・バグ・トラブル
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
我々がジオリファレンスなどでラスタのセル座標を扱う時、ラスタの原
点は、左上隅のセルの左上隅で、(0, 0)です。これは今でも正しいので
すが、最近になってＳＭＬの中では(1, 1)である、に変わったことが分
かりました。マイクロイメージ社は、V7.0からこのような変更をしたよ
うです。
V7.0や7.1を使って、地図画像などを取り込んだ場合、１ピクセルほどコ
ントロールポイントの位置がずれている可能性があります。
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
◆◇ お 知 ら せ ◆◇
..................................................................
１年以上前のメールマガジンを公開しました。
http://www.opengis.co.jp/htm/kako_mail/index_kako.htm
..................................................................
今週の話題 / メーリングリスト
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
特にありませんでした。
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
メーリングリストに参加されたい方は、
info@opengis.co.jp までご連絡ください。
弊社からの対応は バグ対応と、新製品告知のみですが、
ユーザーさまの間での情報交換の場として
利用されています。
■■■■■■■ 現 在 の T N T m i p s の 価 格 ■■■■■■
平成17年11月21日付けでＴＮＴ製品の価格を改定致しました。
現在の Ver 7.1 の価格は
・TNTmips シングルライセンス 890,400円(送料・税込)
・TNTedit シングルライセンス 537,600円(送料・税込)
・TNTview シングルライセンス 81,900円(送料・税込)
・TNTmips 年間バージョンアップ(2回分) 147,000円(税込)
・年間テクニカルサポート 94,500円(税込)
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
困ったときの問題解決フローチャート
▼STEP1「ウェブサイトで検索してみる」
▼STEP2「ウェブサイトにある情報を眺めてみる」
▼STEP3「ドキュメントを読んでみる」
▼STEP4「メーリングリストへ聞いてみる」
▼STEP5「テクニカルサポート」
STEP1 「ウェブ検索テクニック」
▽オープンＧＩＳ検索システム
http://www.opengis.co.jp/htm/namazu.cgi
▽メールマガジン検索システム
http://www.opengis.co.jp/mail_mag/namazu.cgi
▽マイクロイメージ検索サイト
http://www.microimages.com/search/
▽Google 検索
http://www.google.co.jp/
とくに、and検索(キーワードの間に空白)、not検索(キーワードの前に
- や not を入れる)、とは検索(キーワードの後ろに "とは" をつける)
を使いこなすと、便利です！
STEP2 「ウェブサイトの有用な情報源」
▽オープンGIS ウェブサイト「基本操作」
http://www.opengis.co.jp/htm/basic/title.htm [日本語]
▽オープンGIS ウェブサイト「メニュー解説(Flash)」
http://www.opengis.co.jp/flash/index.htm [日本語]
▽オープンGIS ウェブサイト「SML スクリプトの解説」
http://www.opengis.co.jp/sml/sml_index.htm [日本語]
▽オープンGIS ウェブサイト「作業で必要となる情報を集約したポスター類」
http://www.opengis.co.jp/htm/gallery/index.htm [日本語]
▽オープンGIS ウェブサイト「メルマガバックナンバー」
http://www.opengis.co.jp/mail_mag/index.htm [日本語]
STEP3 「ドキュメントを読んでみる」
▽オープンGIS ウェブサイト「その他ドキュメント一覧ページ」
http://www.opengis.co.jp/htm/documents.htm [日本語]
▽オープンGIS ウェブサイト「チュートリアルテキスト＆リモセンテキスト」
http://www.opengis.co.jp/htm/info/book.htm [日本語]
▽マイクロイメージウェブサイト「ワンポイントテクニック集」
http://www.microimages.com/didyouknow/ [英文]
▽マイクロイメージウェブサイト「TNT入門(Getting Started)」
http://www.opengis.co.jp/htm/getstart/getstart.htm [日本語]
http://www.microimages.com/getstart/ [英文]
▽マイクロイメージウェブサイト「リファレンスマニュアル」
http://www.microimages.com/refman/
STEP4「メーリングリストへ聞いてみる」
メーリングリストに参加されたい方は、
info@opengis.co.jp までご連絡ください。
弊社からの対応は バグ対応と、新製品告知のみになります。
STEP5「テクニカルサポート」
テクニカル・サポート・ユーザー専用 メールアドレス
############# このメールマガジンの配信について ################
このメールマガジンは、
年間テクニカル・サポートに加入されている、
ユーザーさまのみを対象としております。
基本的には、登録アカウント１つに対して
１ユーザーを配信対象としておりますが、
例えば、研究室の予算で加入された場合などは、
ご購入いただいた指導教官と、実際に TNTmips で作業を行う学生1名の、
計２ユーザーを対象として配信させていただきます。
１アカウントに対して合計３名以上の配信は行いません。
学生さん２名が同じくらい使用する場合は、
申し訳ありませんが、じゃんけん等で１人に決めてください。
また、テクニカル・サポートの期限を３ヶ月以上過ぎてしまった場合は、
メールマガジンの配信を停止させていただきます。
メールマガジンが配信されなくなった時は、
サポート期限が３ヶ月過ぎたとお考えください。
転勤、ドメイン変更などでメールアドレスが変わった場合は、
お手数ですが、弊社までご連絡ください。
テクニカル・サポート活用案内
>サービスその１：緊急なバグ情報はすぐにメール配信します。
>サービスその２：TNTmips に関するどんな質問にも TEL/e-mail/FAX で対応。
>サービスその３：対応は最優先で行います。
>サービスその４：専用FAX質問用紙をお送りします。
>サービスその５：週に１回メールマガジンを配信します。
>サービスその６：展示ブースへご来場の方には特製オリジナルグッズを進呈。
>サービスその７：サポート対応で FTP サイトを 20MB まで利用できます。
>サービスその８：メルマガ・バックナンバー・サイトも利用できます。
>サービスその９：FTP サイトを 500MB に拡張しました...
>サービスその１０：現在計画中です...
＝＝＝＝＝＝＝＝＝＝＝＝ おことわり ＝＝＝＝＝＝＝＝＝＝＝＝＝
※OpenGIS Mail Magazine for T.S.Uでは、
TNTmipsに関わる新しいニュースを、
毎週金曜日に、皆さまへご提供させていただきます。
※このメールマガジンの転載・転送はご遠慮ください。
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝