Top

個別メモ

Latest update on 2012年3月5日 (月) at 10:54:46.

【第1617回】 Just do what to do（2010年6月22日）

• 6:00に起きたがやるべき家事が多くて家を出るのが遅くなった。予報では雨になるらしいが自転車で長野駅に出て，往路あさま512号。群馬側は新前橋からバスにした。
• ともかくいろいろやらねばならないことが山積みなので，ただやるべきことをやるのみ。
• 来週の月曜日にChoropleth mapの作り方をインドネシアのデータでやろうと思い，準備を始めた。インドネシアの地図データシェイプファイルはDIVA-GISから入手できるが（大抵の国のデータはある。ソロモン諸島はないのだが……），いくつか問題がある。インドネシアの地図データzipファイルを展開すると，3種類のシェイプファイルができ，それぞれ区分の細かさが異なる。IDN_ADM0は島ごとの区分だけ，IDN_ADM1はProvinceごと，IDN_ADM2はRegency/Cityごとの区分となっている。しかし，最近の統計データが示されているProvince区分と，IDN_ADM1のProvince区分が合わないのだ。例えば，2010年3月のAIDS統計には33のProvinceが載っているが，IDN_ADM1は26しかなく，PapuaとWest PapuaをひっくるめてIrian Jayaとして扱われている。
• IDN_ADM2は非常に細かい区分なので，全部表示すると何が何やらわけがわからなくなる（ちょうど日本の市町村区分で日本全国を塗り分けるような感じになる）。そこで，その中から1つのProvince分だけを抜き出そうと思ったのだが，これがうまくいかない。maptoolsライブラリがS3クラスを扱えた頃はsubset関数で簡単に部分の取り出しもできたのだが，read.ShapePolyで読みこまれるS4クラス地図オブジェクトの部分集合をどうやって取り出すのかがわからない。EPIINFOのEpiMapを使えば当然できるので，それで部分集合を作って別ファイルにしておき，それをRから読むことにすれば問題は解決するだろうけれども，どうにも悔しいところ。
• バスで前橋に出て，復路あさま549号。

▼前【1616】（見学実習など（2010年6月21日） ）　▲次【1618】（SPARC Japanセミナー（2010年6月23日） ）　●Top

△Read/Write COMMENTS

Notice to cite or link here | [TOP PAGE]