478 Views
April 03, 24
スライド概要
[第6回大阪sas勉強会]筒井杏奈
SAS言語を中心として,解析業務担当者・プログラマなのコミュニティを活性化したいです
2020年1月31日 第6回大阪SAS勉強会 SASの 日本地図データ MAPSGFK.JAPAN に関して 大阪大学大学院医学系研究科保健学専攻 総合ヘルスプロモーション科学講座 数理保健学教室 筒井杏奈
本日の内容 2 一般的に、SASで日本地図を描く場合、デフォルトの地図データを利用する MAPSGFK.JAPANデータ MAPS.JAPANデータ ただし、あの地域がデータにないため描画できない 例: SAS Institute Japan (株) 東京本社の場所* 本日はその対応方法を紹介します *データステップ100万回 SAS新手一生 (2017) より転載 sgplotで日本地図を書いて新入社員の気をひこうって話 http://sas-tumesas.blogspot.com/2017/04/sgplot.html
3 結論 デフォルトのJAPANデータでなく、 他の境界情報を持つデータセットを使用する シェープファイルをMAPIMPORTプロシジャによりインポートして作成可能 参考サイト: SAS/GRAPH® and Base SAS® 9.4: Mapping Referenceより Map Data Sets, Map Preparation Procedures, and Tools Provided by SAS https://go.documentation.sas.com/?docsetId=grmapref&docsetVersion=9.4_01&docsetTarget=p1jwzmv15qfcg2n1luotzf7qkfud.htm&locale =en#n1sys ryvnwg1m4n1qx4ehe6fx8mu What Does the MAPIMPORT Procedure Do? https://go.documentation.sas.com/?docsetId=grmapref&docsetVersion=9.4_01&docsetTarget=p1rvpewhocha4yn1gipf920vchn3.htm&locale =en ※ SASのテクニカルサポートにてご教示いただきました
4 実際に試してみる (1) 1. シェープファイルをダウンロード 国土数値情報ダウンロードサービス「国土数値情報 行政区域データ」 http://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-N03-v2_3.html ファイルサイズが大きいため今回は北海道データのみ使用 2. 一つのフォルダに格納 一シェープファイルは複数のファイルを保有
実際に試してみる (2) 3. 読み込み PROC MAPIMPORT OUT=myworld DATAFILE=“../XXX.shp"; ID N03_007 ; run; 5
6 実際に試してみる (3) MAPSGFK.JAPAN (10,645 obs. for Hokkaido data) シェープファイルから作成 (1,422,498 obs.) 0 1 • より太い境界線 • 横長の形 北海道の地図をGMAPプロシジャで描画 • データは適当な乱数を発生させて作成したもの
7 補足 MAPIMPORTプロシジャのIDステートメントを未指定の場合 不要な直線が発生 0 1 GMAPプロシジャで指定した 地図データのIDキーが、 各ポリゴンに対してUniqueで ないために生じた *SAS Support Communities (2011) proc mapimport and gmap in SAS 9.1 https://communities.sas.com/t5/Graphics-Programming/proc-mapimport-and-gmap-in-SAS-9-1/m-p/8322/highlight/true#M13
Thank you 今回の発表に政治的な意図はありません 8