kmlラボ
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
&color(black,#CCFFCC){LINK:[[[GoogleEarthとR]]][[[ShapeFileライブラリ]]][[[kmlラボ]]][[[空間的なデータの分析]]][[[Rでジオコーディング]]][[[RでGPS]]]};~
----
目次
#contents
----
*目的 [#f6347ee0]
rgdalパッケージのwriteOGR関数を使うと、ShapeFileからKMLファイルを作成することができますが、単にKMLファイルを出力するだけで、細かい設定ができず、あまり利用価値がありません。~
そこでXMLパッケージを用いて、作成されたKMLファイルを操作することで、KMLファイルの利用価値を高めることができると考えました。（ついでにXMLパッケージをつかいこなすぞ！）~
注）ここは、「KMLの公開」を目的としたサイトではありません。公開されているKMLファイルは、あくまでもRで加工するための基礎データにすぎません。~
*対応表 [#d03e0725]
RとXMLパッケージのバージョンによってはうまく動かないので、注意してください。~
推奨環境：Windows: R2.9.2 - XML1.96~
-Windows: R2.9.2 - XML1.9.6 O
-Windows: R2.7.2 - XML1.94 O (CDATA内部の日本語文字化けなし）~
install.packages("XML", repos = "http://www.omegahat.org/R")~
-Windows: R2.7.2 - XML1.96 X （CDATA内部の日本語文字化け）~
-Windows: R2.6.2 - XML1.94 O~
-MacOSX10.5.5: R2.7.2 O~
~
*関連リンク [#w425a615]
すばらしいサイトです。ここで勉強させてもらってます。~
[[TME:http://thematicmapping.org/]]~
[[EarthAtlas:http://earthatlas.info/?kml=http://www.nature.com/nature/googleearth/avianflu1.kml]]~
*kmlライブラリ [#qca3a38d]
公開データの著作権等について:~
＊「地域境界KMLFileデータの作成に当たっては、ESRIジャパン株式会社の全国市区町村界データを使用しました。本データの著作権はESRIジャパン株式会社に帰属します。なお、地域境界KMLFileデータは営利目的には利用できません。」~
＊「地域メッシュKMLFileデータの作成にあたっては、財団法人 地域地盤環境研究所 村上貴志氏のメッシュデータを使用しました。本データの著作権は村上貴志氏に帰属します。また、地域メッシュKMLFileデータの利用に関してですが，著作者の表記（書式自由）をして頂ければ自由にご利用ください。公開されたデータが非営利・営利いずれの使われ方をしても構いません。」~
|50|150|250|100|300|100|c
|県コード|CENTER:地域名称|CENTER:地域境界KMLFile|CENTER:地域メッシュ名称|CENTER:地域メッシュKMLFile|CENTER:備考|h
|#|日本全図|#ref(kmlライブラリ/japank.kml)|第１次地域区画|#ref(kmlライブラリ/StdMesh01.kml)||
|#|日本全図||第２次地域区画|#ref(kmlライブラリ/StdMesh02.zip)||
|０１|北海道|#ref(kmlライブラリ/hokkaido.kml)|1kmメッシュ|||
|０２|青森県|#ref(kmlライブラリ/aomori.kml)|1kmメッシュ|#ref(kmlライブラリ/StdMesh04_JGD_02.zip)||
|０３|岩手県|#ref(kmlライブラリ/iwate.kml)|1kmメッシュ|||
|０４|宮城県|#ref(kmlライブラリ/miyagi.kml)|1kmメッシュ|||
|０５|秋田県|#ref(kmlライブラリ/akita.kml)|1kmメッシュ|||
|１０|群馬県|#ref(kmlライブラリ/gunma.kml)||||
**国勢調査のデータのダウンロード [#o3ee385b]
http://www.e-stat.go.jp/SG2/toukeichiri/TopFrame.do?fromPage=init&toPage=download
*コメント欄 [#w89db461]
#comment(below)
- XML対応表更新 -- [[okinawa]] &new{2009-11-30 (月) 18:21:41};
- maptoolsにkmlLine,kmlPolygon関数が追加されました。これらの関数にて生成されるkmlはwriteOGRで生成されるkmlと構造が違うためkmlFunctions.rはそのままでは使えません。 -- [[okinawa]] &new{2009-02-20 (金) 11:26:41};
- さて、そろそろ２次メッシュにいこか・・・。 -- [[okinawa]] &new{2009-02-17 (火) 17:31:12};
- メッシュkmlデータの作成:1次メッシュ生成関数作成終了。 -- [[okinawa]] &new{2009-02-13 (金) 16:52:40};
- メッシュkmlデータの作成更新。メッシュを連続生成できるようにしました。まだ作りかけですが。 -- [[okinawa]] &new{2009-02-13 (金) 09:07:23};
- メッシュkmlデータの作成を追加。現在、初期段階のコードです。 -- [[okinawa]] &new{2009-02-12 (木) 14:50:36};
- GE5向けにkmlFunctions.rを変更しました。今後、ここの内容はGE5向けに変更していきますのでご了承ください。 -- [[okinawa]] &new{2009-02-05 (木) 09:33:38};
- データを読み込んでkmlを再構成する（２） 那覇市国勢調査編のコードを若干変更（ポリゴンホールを回避します） -- [[okinawa]] &new{2009-02-03 (火) 21:59:56};
- 今後、TMEの技術を取り入れていこうと思います。・・・まだ自動でメッシュデータ作るのが残ってますが・・・。 -- [[okinawa]] &new{2008-12-19 (金) 09:22:06};
- ここの名称を「kmlラボ」に変更しました。 -- [[okinawa]] &new{2008-12-08 (月) 08:41:24};
- kmlFunctions.rバグ修正。 -- [[okinawa]] &new{2008-11-26 (水) 18:25:28};
- useInternalNodes = TRUEを使わんといかんのか・・・。utf-8が標準じゃないwindowsの方がめんどくさいですね。 -- [[okinawa]] &new{2008-11-18 (火) 09:37:57};
- 那覇市の周辺市町村（浦添市、豊見城市、南風原町）を同時に表示してみた。 -- [[okinawa]] &new{2008-11-17 (月) 17:15:05};
- getCDATA関数復活。これでPlaceMark名にCDATA内の指定データをnameタグに挿入できます。これをsetAltitude()と一緒に使うことで、placemarkに日本語が表示できるようになると思います。 （Macでは成功、Windowsではうまくいかないなあ・・・）-- [[okinawa]] &new{2008-11-16 (日) 13:50:16};
- MacOSXでrgdalを使えるようにするを追加。これでMacでもKML自由自在ですね。 -- [[okinawa]] &new{2008-11-15 (土) 15:26:01};
- 明日、那覇市長選がある。那覇市国勢調査編として分析したが、３Dグラフから那覇市はすでに人口が空洞化していると考えられる。このような地域は企業が集中している反面、取り残された住民に高齢者の割合が高くなる。いわゆる都市構造的な問題は、一自治体の政策だけでは解決できないと考えられるが、いかがだろうか？ -- [[okinawa]] &new{2008-11-15 (土) 12:41:58};
- kmlFunctions.rのsetAltitude()にラベルのパラメータ指定を追加。ラベルを指定することで、KML内のPlaceMark名をセットする。ただし、現状では日本語をセットするとエラーになる。XMLパッケージのロケール問題っぽい。 -- [[okinawa]] &new{2008-11-14 (金) 14:53:02};
- kmlFunctions.rバグ取り＆ちょっと修正。 -- [[okinawa]] &new{2008-11-12 (水) 15:30:15};
- データを読み込んでkmlを再構成する（２） 那覇市国勢調査編追加 -- [[okinawa]] &new{2008-11-07 (金) 14:49:09};
- kmlAltitudeNode関数修正、setAltitudeMesh関数をsetAltitude関数に統合。MGパラメータでメッシュKMLと切り替えます。 -- [[okinawa]] &new{2008-11-07 (金) 11:56:58};
- 次はメッシュデータを自動でつくってみます。 -- [[okinawa]] &new{2008-11-06 (木) 15:20:27};
- データを読み込んで、kmlを再構成するを追加。shapefileから直接つくるのにくらべ簡素に？なりました。 -- [[okinawa]] &new{2008-11-06 (木) 14:50:28};
- kmlFunctions.r追加。 -- [[okinawa]] &new{2008-11-06 (木) 11:06:04};
- colors.kml関数追加しました。 -- [[okinawa]] &new{2008-11-05 (水) 16:19:15};
- getKmlNode関数追加。ノードを指定してKMLの切り抜きができます。 -- [[okinawa]] &new{2008-10-24 (金) 14:47:32};
- 関数追加。あとは色の変換関数がいるな。 -- [[okinawa]] &new{2008-10-23 (木) 21:50:06};
- 関数を若干修正 -- [[okinawa]] &new{2008-10-23 (木) 08:56:03};
- 地域メッシュKMLデータの応用追加 -- [[okinawa]] &new{2008-10-22 (水) 18:22:05};
- kmlファイルをRから直接いじってみました。今後は、ライブラリを充実させ、色々な機能追加をしてみようと思います。 -- [[okinawa]] &new{2008-10-21 (火) 18:03:29};
*kmlデータをXMLパッケージで再構成する関数 [#p7c92b92]
writeOGRで出力されたkmlデータをXMLパッケージを用いて再構成し、高さ、ポリゴンの色、ポリゴンの線の色を指定できる関数を作成しました。~
(構造決めうちなので、[[ShapeFileライブラリ]]のデータをwriteOGRでkml出力したデータしか使えないと思います）~
***kmlFunctions.r [#rb911ac6]
kmlを操作する関数をまとめました。下記のように、適当な場所に保存して、souceで呼び出してください。（GE5向けに新しくしました）~
#ref(kmlラボ/kmlFunctions.r)~
古いものは↓
#ref(kmlライブラリ/kmlFunctions.r)~
library(XML) #宣言
source("c:/kmlFunctions.r")#ソースの読み込み
doc <- xmlTreeParse("c:/gunma.kml")#土台となるkmlファイル読み込み
data<-sample(1000:10000,39) #高度データを自動生成
fillcolor<-colors.kml(heat.colors(39),"EE")#色の割り当て（ポリゴン）
linecolor<-NA#色の割り当て（線）
x<-kmlAltitudeNode(doc,data,fillcolor,linecolor)#kmlファイルにデータを書き込み
saveXML(x,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')#ファイルに出力
~
#ref(kmlライブラリ/color2.JPG,30%)~
***kmlAltitudeNode関数(GE5) [#n3382d3b]
kmlAltitudeNode<-function(doc,data,fillcolor,linecolor,extrude="1",tessellate="1",altitudeMode="absolute",MG=T){
kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
kmlver5<-'http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2'
kmlnode<-xmlNode("kml",node0,attrs=c(xmlns=kmlver5))
len1<-length(node0[[1]])
for(i in 2:len1){
if(MG==T){
node1<-node0[[1]][[i]][[2]][[1]]
}else if(MG==F){
node1<-node0[[1]][[i]][[2]]
}
node1<-addChildren(node1,xmlNode("extrude",extrude),xmlNode("tessellate",tessellate),xmlNode("gx:altitudeMode",altitudeMode))
node2<-node0[[1]][[i]][[3]][[2]]
node2<-xmlNode("PolyStyle",xmlNode("color",fillcolor[i-1]))
#linecolor
node3<-node0[[1]][[i]][[3]][[1]]
node3<-xmlNode("LineStyle",xmlNode("color",linecolor))
#polygon
len2<-length(node0[[1]][[i]][[2]][[1]])
for(j in 1:len2){
x1<-xmlValue(node1[[j]][[1]][[1]])
x2<-rev(unlist(strsplit(x1," ")))
x3<-paste(x2,data[[i-1]],sep=",",collapse=" ")
node1[[j]][[1]][[1]]<-xmlNode("coordinates",x3)
}
#replace node
if(MG==T){
node0[[1]][[i]][[2]][[1]]<-node1
}else if(MG==F){
node0[[1]][[i]][[2]]<-node1
}
node0[[1]][[i]][[3]][[2]]<-node2
node0[[1]][[i]][[3]][[1]]<-node3
kmlnode[[1]]<-node0
}
return(kmlnode)
}
#sample
library(XML)
doc <- xmlTreeParse("c:/gunma.kml")
data<-sample(1000:10000,39)
fillcolor<-c("ff000000","ffff0000","ff00ff00","ff0000ff","ffffff00","ffff00ff","ff00ffff","ffffffff","00000000")
linecolor<-"NA"
x<-kmlAltitudeNode(doc,data,fillcolor,linecolor)
saveXML(x,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/xml1.JPG,30%)~
#ref(kmlラボ/gunma2.kml)~
***setAltitude関数（GE5) [#nb393eab]
ノード番号を指定してピンポイントでAltitudeをセットする。~
setAltitude<-function(kmlnode,nodeno,data,fillcolor,linecolor,extrude="1",tessellate="1",altitudeMode="absolute",MG=T,label=""){
#kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
kmlver5<-'http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2'
kmlnode<-xmlNode("kml",node0,attrs=c(xmlns=kmlver5))
if(MG==T){
#altitude
node1<-node0[[1]][[nodeno]][[2]][[1]]
node1<-addChildren(node1,xmlNode("extrude",extrude),xmlNode("tessellate",tessellate),xmlNode("gx:altitudeMode",altitudeMode))
#fillcolor
node2<-node0[[1]][[nodeno]][[3]][[2]]
node2<-xmlNode("PolyStyle",xmlNode("color",fillcolor))
#linecolor
node3<-node0[[1]][[nodeno]][[3]][[1]]
node3<-xmlNode("LineStyle",xmlNode("color",linecolor))
#label
node4<-node0[[1]][[nodeno]]
node4<-addChildren(node4,xmlNode("name",label))
#polygon
len2<-length(node0[[1]][[nodeno]][[2]][[1]])
j<-1:len2
text1<-paste('x1<-xmlValue(node1[[',j,']][[1]][[1]]);x2<-rev(unlist(strsplit(x1," ")));',sep='')
text2<-paste('x3<-paste(x2,data,sep=",",collapse=" ");node1[[',j,']][[1]][[1]]<-xmlNode("coordinates",x3)',sep='')
text<-paste(text1,text2,sep="")
eval(parse(text=text))
node0[[1]][[nodeno]]<-node4
node0[[1]][[nodeno]][[2]][[1]]<-node1
node0[[1]][[nodeno]][[3]][[2]]<-node2
node0[[1]][[nodeno]][[3]][[1]]<-node3
kmlnode[[1]]<-node0
}else if(MG==F){
#altitude
node1<-node0[[1]][[nodeno]][[2]]
node1<-addChildren(node1,xmlNode("extrude",extrude),xmlNode("tessellate",tessellate),xmlNode("gx:altitudeMode",altitudeMode))
#fillcolor
node2<-node0[[1]][[nodeno]][[3]][[2]]
node2<-xmlNode("PolyStyle",xmlNode("color",fillcolor))
#linecolor
node3<-node0[[1]][[nodeno]][[3]][[1]]
node3<-xmlNode("LineStyle",xmlNode("color",linecolor))
#label
node4<-node0[[1]][[nodeno]]
node4<-addChildren(node4,xmlNode("name",label))
#polygon
len2<-length(node0[[1]][[nodeno]][[2]][[1]])
j<-1:len2
text1<-paste('x1<-xmlValue(node1[[',j,']][[1]][[1]]);x2<-rev(unlist(strsplit(x1," ")));',sep='')
text2<-paste('x3<-paste(x2,data,sep=",",collapse=" ");node1[[',j,']][[1]][[1]]<-xmlNode("coordinates",x3)',sep='')
text<-paste(text1,text2,sep="")
eval(parse(text=text))
node0[[1]][[nodeno]]<-node4
node0[[1]][[nodeno]][[2]]<-node1
node0[[1]][[nodeno]][[3]][[2]]<-node2
node0[[1]][[nodeno]][[3]][[1]]<-node3
kmlnode[[1]]<-node0
}
return(kmlnode)
}
getCDATANo関数と組み合わせて必要な場所のみ設定~
#sample
doc <- xmlTreeParse("c:/gunma.kml")
x<-xmlRoot(doc)
#1つめの地域
name<-"JCODE"
value<-"10366"
nodeno<-getCDATANo(x,name,value)
data<-1000
fillcolor<-c("ffff0000")
linecolor<-"ff000000"
x<-setAltitude(x,nodeno,data,fillcolor,linecolor)
#２つめの地域
name<-"JCODE"
value<-"10201"
nodeno<-getCDATANo(x,name,value)
data<-1000
fillcolor<-c("ffff00ff")
linecolor<-"ff000000"
x<-setAltitude(x,nodeno,data,fillcolor,linecolor)
#保存
saveXML(x,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/xml4.JPG,30%)~
メッシュのKMLはMultiGeometoryノードが無いので、MGパラメータでF指定する~
#sample
library(XML)
doc <- xmlTreeParse("c:/StdMesh01.kml")
kmlnode<-xmlRoot(doc)
#1
name<-"CODE"
value<-"3622"
nodeno<-getCDATANo(kmlnode,name,value)
data<-1000000
fillcolor<-c("ffff0000")
linecolor<-"ff000000"
kmlnode<-setAltitude(kmlnode,nodeno,data,fillcolor,linecolor,MG=F)
#2
name<-"CODE"
value<-"3623"
nodeno<-getCDATANo(kmlnode,name,value)
data<-800000
fillcolor<-c("ffff00ff")
linecolor<-"ff000000"
kmlnode<-setAltitude(kmlnode,nodeno,data,fillcolor,linecolor,MG=F)
#save
saveXML(kmlnode,"c:/StdMesh012.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/color4.JPG,30%)~
***kmlNodeLength関数 [#y7fb430a]
kmlAltitudeNode関数を使うときに、kmlのノード数を知りたい時に使う。~
#kmlNodeLength() ver 0.1.2
kmlNodeLength<-function(doc){
library(XML)
kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
len1<-length(node0[[1]])
return(len1)
}
#sample
library(XML)
doc <- xmlTreeParse("c:/gunma.kml")
kmlNodeLength(doc)
> #sample
> library(XML)
> doc <- xmlTreeParse("c:/gunma.kml")
> kmlNodeLength(doc)
[1] 40
注）先頭のノードに<name>gunma</name>が入っていますので、市町村数はノード数−１となります。~
***getCDATAList関数 [#ic0a4023]
CDATA(HTML)の中身を取り出す。nameに対応する項目名を指定する。（必須）~
getCDATAList<-function(doc,name){
getCDATA<-function(node,name){
x<-gsub(" ","",iconv(xmlValue(node),"utf-8","cp932"))
x<-unlist(strsplit(x,"<br/>\n"))
x<-x[[grep(name,x)]]
x<-gsub("<b>","",gsub("</b>","",gsub("<i>","",gsub("</i>","",gsub("<br/>","",gsub(paste(name,":",sep=""),"",x))))))
#x<-iconv(x,"","utf-8")
return(x)
}
kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
len1<-length(node0[[1]])
x<-NULL
for(i in 2:len1){
x[i-1]<-getCDATA(node0[[1]][[i]][[1]][[1]],name)
}
return(x)
}
#sample
doc <- xmlTreeParse("c:/gunma.kml")
name<-"CITY1"
getCDATAList(doc,name)
> #sample
> doc <- xmlTreeParse("c:/gunma.kml")
> name<-"CITY1"
> getCDATAList(doc,name)
[1] "前橋市" "高崎市" "桐生市" "伊勢崎市" "太田市" "沼田市" "館林市" "渋川市" "藤岡市" "富岡市"
[11] "安中市" "みどり市" "富士見村" "榛名町" "榛東村" "吉岡町" "吉井町" "上野村" "神流町" "下仁田町"
[21] "南牧村" "甘楽町" "中之条町" "長野原町" "嬬恋村" "草津町" "六合村" "高山村" "東吾妻町" "片品村"
[31] "川場村" "昭和村" "みなかみ町" "玉村町" "板倉町" "明和町" "千代田町" "大泉町" "邑楽町"
> #sample
> doc <- xmlTreeParse("c:/gunma.kml")
> name<-"JCODE"
> getCDATAList(doc,name)
[1] "10201" "10202" "10203" "10204" "10205" "10206" "10207" "10208" "10209" "10210" "10211" "10212" "10303" "10321" "10344" "10345" "10363"
[18] "10366" "10367" "10382" "10383" "10384" "10421" "10424" "10425" "10426" "10427" "10428" "10429" "10443" "10444" "10448" "10449" "10464"
[35] "10521" "10522" "10523" "10524" "10525"
***getCDATA関数 [#vf724883]
CDATAのデータを取り出す。cdatanoは、grepで検索した場合、複数の位置が取り出された場合、初期位置からの相対位置を指定するパラメータ。~
getCDATA<-function(kmlnode,name,nodeno=1,sep=":",cdatano=1){
node0<-kmlnode[[1]]
node1<-node0[[1]][[nodeno]][[1]][[1]]
#x<-gsub(" ","",xmlValue(node1))#mac
x<-gsub(" ","",iconv(xmlValue(node1),"utf-8","cp932"))
#print(x)
x<-unlist(strsplit(x,"<br/>\n"))
xno<-grep(name,x)
len<-length(xno)
if(len>1){
x<-x[[xno[cdatano]]]
}else if(len==1){
x<-x[[xno]]
}
#print(x)
x<-gsub("<b>","",gsub("</b>","",gsub("<i>","",gsub("</i>","",gsub("<br/>","",x)))))
x<-gsub(paste(name,sep,sep=""),"",x)
return(x)
}
> #sample
> doc <- xmlTreeParse("/gunma.kml")
> kmlnode<-xmlRoot(doc)
> name<-"JCODE"
> nodeno<-2
> getCDATA(kmlnode,name,nodeno)
[1] "10201"
***getCDATANo関数 [#x4aa4dac]
nameとvalueを指定することで、CDATA(HTML)内のデータを検索し、一致したノード番号を返す。~
#getCDATANo
getCDATANo<-function(kmlnode,name,value){
getCDATA2<-function(node,name){
x<-gsub(" ","",iconv(xmlValue(node),"utf-8","cp932"))
x<-unlist(strsplit(x,"<br/>\n"))
x<-x[[grep(name,x)]]
x<-gsub("<b>","",gsub("</b>","",gsub("<i>","",gsub("</i>","",gsub("<br/>","",x)))))
#x<-iconv(x,"","utf-8")
return(x)
}
#kmlnode <- xmlRoot(doc)
namevalue<-paste(name,":",value,sep="")
node0<-kmlnode[[1]]
len1<-length(node0[[1]])
for(i in 2:len1){
if(namevalue==getCDATA2(node0[[1]][[i]][[1]][[1]],name)){
x<-i
break
}else{
x<-NA
}
}
return(x)
}
#sample
doc <- xmlTreeParse("c:/gunma.kml")
kmlnode<-xmlRoot(doc)
name<-"JCODE"
value<-"10366"
getCDATANo(kmlnode,name,value)
> #sample
> doc <- xmlTreeParse("c:/gunma.kml")
> name<-"JCODE"
> value<-"10366"
> getCDATANo(kmlnode,name,value)
[1] 19
***getKmlNode関数(GE5) [#udad8667]
ノードをnodenoで指定し、指定したノードの部分だけ抽出してKMLを再構成~
#getKmlNode
getKmlNode<-function(kmlnode,nodeno){
len<-length(nodeno)
node0<-kmlnode[[1]]
kmlver5<-'http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2'
kmlnode<-xmlNode("kml",node0,attrs=c(xmlns=kmlver5))
node2<-xmlNode("Folder",xmlNode("name",as.character(nodeno)))
for(i in 1:len){
node1<-node0[[1]][[nodeno[[i]]]]
node2<-addChildren(node2,node1)
}
kmlnode[[1]][[1]]<-node2
return(kmlnode)
}
#sample
library(XML)
doc <- xmlTreeParse("c:/gunma.kml")
kmlnode<-xmlRoot(doc)
nodeno<-c(2,4,6)
expdata<-getKmlNode(kmlnode,nodeno)
saveXML(expdata,"c:/gunma3.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
***colors.kml関数 [#r236dd3e]
colors.kml("色名","アルファチャンネル")
colors.kml<-function(rgb=colors(),alpha="FF"){
colors.hex <- function( x=colors() ) {
color.hex <- function(x) do.call( "rgb", as.list(col2rgb(x)/255) )
sapply( x, color.hex )
}
col.hex<-colors.hex(rgb)
red<-substr(col.hex,2,3)
green<-substr(col.hex,4,5)
blue<-substr(col.hex,6,7)
alpha<-alpha
col.kml<-paste(alpha,blue,green,red,sep="")
return(col.kml)
}
> colors.kml("red")
[1] "FF0000FF"
> colors.kml(rainbow(10))
[1] "FF0000FF" "FF0099FF" "FF00FFCC" "FF00FF33" "FF66FF00" "FFFFFF00" "FFFF6600" "FFFF0033" "FFFF00CC"
[10] "FF9900FF"
参考：[[色見本]]~
*メッシュkmlデータの作成 [#t6141a60]
メッシュkmlデータを自動作成します。生成されたkmlはもちろん他の関数で操作できます。
**makeMesh1関数：第１次地域区画（１次メッシュ）生成関数 [#f324b924]
指定された地域メッシュコードを基に１次メッシュをブロック生成します。~
下記のDefaultMesh1.kmlをタネにして作成しますので、ダウンロードして利用してください。~
#ref(kmlラボ/DefaultMesh1.kml)~
makeMesh1<-function(kmlnode,codea,codeb){
transMeshCode<-function(code){
lat0<-as.integer(substr(code,1,2))/1.5
lon0<-as.integer(substr(code,3,4))+100
lat1<-lat0;lon1<-lon0+1;lat2<-lat0+(2/3);lon2<-lon1
lat3<-lat2;lon3<-lon0
x<-paste(lon0,",",lat0," ",lon1,",",lat1," ",lon2,",",lat2," ",lon3,",",lat3," ",lon0,",",lat0,sep="")
return(x)
}
code1<-as.integer(substr(codea,1,2))
code3<-as.integer(substr(codea,3,4))
code2<-as.integer(substr(codeb,1,2))
code4<-as.integer(substr(codeb,3,4))
print(code1)
print(code2)
print(code3)
print(code4)
node0<-kmlnode[[1]]
node1<-node0[[1]]
node2<-node0[[1]][[2]]
x<<-node0[[1]]
nodex<-removeChildren(x,kids=c("Placemark"))
for(j in code1:code2){
code01<-paste(as.character(j),code3,sep="")
code02<-paste(as.character(j),code4,sep="")
print(j)
for(i in code01:code02){
print(i)
node3<-node0[[1]][[2]][[1]][[1]]
node4<-node0[[1]][[2]][[2]][[1]][[1]][[1]]
CDATA<-paste("<b>CODE:</b><i>",i,"</i><br />\n",sep="")
node3<-xmlCDataNode(CDATA)
node4<-xmlNode("coordinates",transMeshCode(i))
node0[[1]][[2]][[2]][[1]][[1]][[1]]<-node4
node0[[1]][[2]][[1]][[1]]<-node3
node2<-node0[[1]][[2]]
nodex<-addChildren(nodex,node2)
}
node0[[1]]<-nodex
}
kmlnode[[1]]<-node0
return(kmlnode)
}
#Sample
library(XML)
doc <- xmlTreeParse("c:/DefaultMesh1.kml")
kmlnode <- xmlRoot(doc)
kmlnode2<-makeMesh1(kmlnode,3622,3643)#2つの区画を指定することで、２区画間の区画をブロック生成する
saveXML(kmlnode2,"c:/meshdata.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlラボ/mesh2.JPG,30%)~
１次メッシュ生成~
#Sample2:少々時間がかかります
library(XML)
doc <- xmlTreeParse("c:/DefaultMesh1.kml")
kmlnode <- xmlRoot(doc)
kmlnode2<-makeMesh1(kmlnode,3622,6853)
saveXML(kmlnode2,"c:/meshtest.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlラボ/mesh4.JPG,30%)~
*応用編 [#me6a710e]
***データを読み込んでkmlを再構成する（１） [#cce8c1cb]
csvデータを読み込んで、kmlを再構成しています。JCODEでマッチングさせていますので、事前にデータをノード順と一致させる必要はありません。~
#ref(kmlラボ/agedprop.csv)~
library(XML)
source("c:/kmlFunctions.r")
#kmlファイルの読み込み
doc <- xmlTreeParse("c:/gunma.kml")
kmlnode<-xmlRoot(doc)
#csvデータの取得
dat<-read.csv("c:/agedprop.csv")
data<-dat$AP2006
Altdata<-data*1000#高度
name<-"JCODE"#マッチング用指定
#色階層の割り当て
cutnum<-16
classes <- cut(data,seq(min(data),max(data),length=cutnum+1),include.lowest=T)
cols <- colors.kml(rev(heat.colors(cutnum)))
cols2<-cols[classes]
#kmlへの書き込み
len<-kmlNodeLength(doc)-1
for(i in 1:len){
value<-dat$JCODE[[i]]
print(value)
nodeno<-getCDATANo(kmlnode,name,value)#JCODEからノードを取得
print(nodeno)
fillcolor<-cols2[[i]]
linecolor<-NA
kmlnode<-setAltitude(kmlnode,nodeno,Altdata[[i]],fillcolor,linecolor)
}
#kml保存
saveXML(kmlnode,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/color3.JPG,30%)~
***データを読み込んでkmlを再構成する（２） 那覇市国勢調査編 [#e7808495]
まず、那覇市国勢調査のデータ（平成１７年）を地図で見る統計（統計ＧＩＳ）のデータダウンロードから落とします。~
http://www.e-stat.go.jp/SG2/toukeichiri/TopFrame.do?fromPage=init&toPage=download~
（１）ダウンロードページ
#ref(kmlライブラリ/naha1.JPG,30%)
（２）世界測地系緯度経度Shape形式を落とします。
#ref(kmlライブラリ/naha2.JPG,30%)
（３）前準備~
落としたShapeFileから、writeOGRでKMLデータ（h17ka47201.kml)を作成します。~
library(maptools)
library(rgdal)
x<-readShapePoly("c:/nahashi/h17ka47201.shp")
x@data$KEN_NAME<-iconv(x@data$KEN_NAME,"","utf-8")
x@data$GST_NAME<-iconv(x@data$GST_NAME,"","utf-8")
x@data$MOJI<-iconv(x@data$MOJI,"","utf-8")
tf <- "c:/nahashi/"
writeOGR(x, paste(tf, "h17ka47201.kml", sep=""), "h17ka47201", "KML")
（４）3Dグラフの作成~
作成したkmlとShapeFileに付属しているdbfを使って3Dグラフを作成します。~
国勢調査データから作成したkmlはMultiGeometryタグが無いので、setAltitude関数をMG=Fで取り込みます。~
ここで、表示させるのは人口のデータです。~
library(XML)
library(foreign)
#path<-"/"#Mac
path<-"c:/"
source(paste(path,"kmlFunctions.r",sep=""))
doc <- xmlTreeParse(paste(path,"nahashi/h17ka47201.kml",sep=""))
kmlnode<-xmlRoot(doc)
#データの取得
dat<-read.dbf(paste(path,"nahashi/h17ka47201.dbf",sep=""))
data<-dat$JINKO#人口データ指定(ここを変えることで他の項目も作れます）
Altdata<-data+10#高度:ポリゴンホールが出来ないように全体を１０m上げる
name<-"SEQ_NO2"#マッチング用指定
#色階層の割り当て
cutnum<-16
classes <- cut(data,seq(min(data),max(data),length=cutnum+1),include.lowest=T)
cols <- colors.kml(rev(heat.colors(cutnum)))
cols2<-cols[classes]
#kmlへの書き込み
len<-kmlNodeLength(doc)-1
for(i in 1:len){
value<-dat$SEQ_NO2[[i]]#マッチング用指定
print(value)
nodeno<-getCDATANo(kmlnode,name,value)#SEQ_NO2からノードを取得
#nodeno<-i+1 #こちらを使うとかなり速いです（ノードとＤＢＦのデータの並びが同じ場合のみ。SEQ_NO2でマッチングさせてないので速いです。）
print(nodeno)
fillcolor<-cols2[[i]]
linecolor<-NA
#Macの場合以下のコードで実行可能
#label<-getCDATA(kmlnode,"MOJI",nodeno,cdatano=2)#PlaceMark名にMOJIを指定。NMOJIも検索してくるので、相対位置パラメータcdatano=2を指定。
#Windowsの場合、いまのところlabelに日本語指定不可
label<-getCDATA(kmlnode,"KEYCODE1",nodeno)
kmlnode<-setAltitude(kmlnode,nodeno,Altdata[[i]],fillcolor,linecolor,MG=F,label=label)
}
#kml保存
saveXML(kmlnode,paste(path,"nahashi/nahashi3.kml",sep=""),prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
作成には、かなりの時間がかかりますので、気長にお待ちください。（私のＰＣで１０分）~
ただし、nodeno<-i+1の方を使うとすぐに作成できます。（１分くらい）ただし、ノードとＤＢＦのデータの並びが同じ場合のみです。~
データは、nahashi3.kmlで落ちてきます。~
赤くて高い場所ほど人口が多い地域です。那覇港は手作業で着色しました。~
#ref(kmlラボ/naha5.JPG,30%)
那覇市の周辺市町村（浦添市、豊見城市、南風原町）を同時に表示してみた。~
明らかに周辺市町村に人口が流れているのがわかる。~
#ref(kmlラボ/city2.JPG,30%)
なお、国勢調査データには利用の制限がありますので、ご注意ください。~
~
利用の制限~
~
１． 利用者は、本システムでダウンロードしたデータ及び画像データをそのまま複製（ファイル形式を変換しての複製を含む。）して第三者に譲渡することを禁じます。
なお、商用目的で複製する場合は、e-Statの「お問い合わせ」よりご連絡ください。
~
２．法律、政令、規則、省令その他すべての法令及び条例等の法規に違反する目的、手段又は方法で、本システムでダウンロードしたデータ及び画像データを利用することを一切禁じます。また、他人の権利を侵害する目的、手段又は方法での利用、公序良俗に反する利用についても一切禁じます。
*MacOSXでrgdalを使えるようにする [#l1ecb7c3]
rgdalパッケージはMacOSXではバイナリで提供されていませんが、下記のサイトを参考にすると、rgdalをMacOSXに導入することができます。（#GDALのインストールと#rgdalのインストールだけやってください）~
//↓の記事書いたものです．BLOGを移転しましたので，新しいアドレスに変更させていただきました．移転前のは近いうちに消すと思います．
//何か問題があればお手数をおかけしてしまい申し訳ありませんが，修正いただければ幸いです．何卒宜しくお願い申し上げます．
//http://www.0093-style.net/0093/2008/02/mac_os_x_leopar_31f8.html~
http://blog.0093.tv/2008/02/mac-os-x-leopardgrassgdalrgdalspgrass6_4571.html~
writeOGRもきちんと動きました。(^_^)/~
** MacOS X でrgdalを使えるようにするもうひとつの方法 [#l9a7d706]
+[[http://www.kyngchaos.com/software:frameworks]]より、GDAL completeをインストール
+ /Library/Frameworks/GDAL.framework/Programs/gdal-config を /usr/local/bin などにシンボリックリンク
+ rgdalのソースパッケージをCRANから取得
+ $ R32 CMD INSTALL rgdal_0.6-20.tar.gz --configure-args='--with-proj-include=/Library/Frameworks/PROJ.framework/Headers --with-proj-lib=/Library/Frameworks/PROJ.Framework/unix/lib'
*アクセス数： [#tdc929b1]
&counter;
人
終了行:
&color(black,#CCFFCC){LINK:[[[GoogleEarthとR]]][[[ShapeFileライブラリ]]][[[kmlラボ]]][[[空間的なデータの分析]]][[[Rでジオコーディング]]][[[RでGPS]]]};~
----
目次
#contents
----
*目的 [#f6347ee0]
rgdalパッケージのwriteOGR関数を使うと、ShapeFileからKMLファイルを作成することができますが、単にKMLファイルを出力するだけで、細かい設定ができず、あまり利用価値がありません。~
そこでXMLパッケージを用いて、作成されたKMLファイルを操作することで、KMLファイルの利用価値を高めることができると考えました。（ついでにXMLパッケージをつかいこなすぞ！）~
注）ここは、「KMLの公開」を目的としたサイトではありません。公開されているKMLファイルは、あくまでもRで加工するための基礎データにすぎません。~
*対応表 [#d03e0725]
RとXMLパッケージのバージョンによってはうまく動かないので、注意してください。~
推奨環境：Windows: R2.9.2 - XML1.96~
-Windows: R2.9.2 - XML1.9.6 O
-Windows: R2.7.2 - XML1.94 O (CDATA内部の日本語文字化けなし）~
install.packages("XML", repos = "http://www.omegahat.org/R")~
-Windows: R2.7.2 - XML1.96 X （CDATA内部の日本語文字化け）~
-Windows: R2.6.2 - XML1.94 O~
-MacOSX10.5.5: R2.7.2 O~
~
*関連リンク [#w425a615]
すばらしいサイトです。ここで勉強させてもらってます。~
[[TME:http://thematicmapping.org/]]~
[[EarthAtlas:http://earthatlas.info/?kml=http://www.nature.com/nature/googleearth/avianflu1.kml]]~
*kmlライブラリ [#qca3a38d]
公開データの著作権等について:~
＊「地域境界KMLFileデータの作成に当たっては、ESRIジャパン株式会社の全国市区町村界データを使用しました。本データの著作権はESRIジャパン株式会社に帰属します。なお、地域境界KMLFileデータは営利目的には利用できません。」~
＊「地域メッシュKMLFileデータの作成にあたっては、財団法人 地域地盤環境研究所 村上貴志氏のメッシュデータを使用しました。本データの著作権は村上貴志氏に帰属します。また、地域メッシュKMLFileデータの利用に関してですが，著作者の表記（書式自由）をして頂ければ自由にご利用ください。公開されたデータが非営利・営利いずれの使われ方をしても構いません。」~
|50|150|250|100|300|100|c
|県コード|CENTER:地域名称|CENTER:地域境界KMLFile|CENTER:地域メッシュ名称|CENTER:地域メッシュKMLFile|CENTER:備考|h
|#|日本全図|#ref(kmlライブラリ/japank.kml)|第１次地域区画|#ref(kmlライブラリ/StdMesh01.kml)||
|#|日本全図||第２次地域区画|#ref(kmlライブラリ/StdMesh02.zip)||
|０１|北海道|#ref(kmlライブラリ/hokkaido.kml)|1kmメッシュ|||
|０２|青森県|#ref(kmlライブラリ/aomori.kml)|1kmメッシュ|#ref(kmlライブラリ/StdMesh04_JGD_02.zip)||
|０３|岩手県|#ref(kmlライブラリ/iwate.kml)|1kmメッシュ|||
|０４|宮城県|#ref(kmlライブラリ/miyagi.kml)|1kmメッシュ|||
|０５|秋田県|#ref(kmlライブラリ/akita.kml)|1kmメッシュ|||
|１０|群馬県|#ref(kmlライブラリ/gunma.kml)||||
**国勢調査のデータのダウンロード [#o3ee385b]
http://www.e-stat.go.jp/SG2/toukeichiri/TopFrame.do?fromPage=init&toPage=download
*コメント欄 [#w89db461]
#comment(below)
- XML対応表更新 -- [[okinawa]] &new{2009-11-30 (月) 18:21:41};
- maptoolsにkmlLine,kmlPolygon関数が追加されました。これらの関数にて生成されるkmlはwriteOGRで生成されるkmlと構造が違うためkmlFunctions.rはそのままでは使えません。 -- [[okinawa]] &new{2009-02-20 (金) 11:26:41};
- さて、そろそろ２次メッシュにいこか・・・。 -- [[okinawa]] &new{2009-02-17 (火) 17:31:12};
- メッシュkmlデータの作成:1次メッシュ生成関数作成終了。 -- [[okinawa]] &new{2009-02-13 (金) 16:52:40};
- メッシュkmlデータの作成更新。メッシュを連続生成できるようにしました。まだ作りかけですが。 -- [[okinawa]] &new{2009-02-13 (金) 09:07:23};
- メッシュkmlデータの作成を追加。現在、初期段階のコードです。 -- [[okinawa]] &new{2009-02-12 (木) 14:50:36};
- GE5向けにkmlFunctions.rを変更しました。今後、ここの内容はGE5向けに変更していきますのでご了承ください。 -- [[okinawa]] &new{2009-02-05 (木) 09:33:38};
- データを読み込んでkmlを再構成する（２） 那覇市国勢調査編のコードを若干変更（ポリゴンホールを回避します） -- [[okinawa]] &new{2009-02-03 (火) 21:59:56};
- 今後、TMEの技術を取り入れていこうと思います。・・・まだ自動でメッシュデータ作るのが残ってますが・・・。 -- [[okinawa]] &new{2008-12-19 (金) 09:22:06};
- ここの名称を「kmlラボ」に変更しました。 -- [[okinawa]] &new{2008-12-08 (月) 08:41:24};
- kmlFunctions.rバグ修正。 -- [[okinawa]] &new{2008-11-26 (水) 18:25:28};
- useInternalNodes = TRUEを使わんといかんのか・・・。utf-8が標準じゃないwindowsの方がめんどくさいですね。 -- [[okinawa]] &new{2008-11-18 (火) 09:37:57};
- 那覇市の周辺市町村（浦添市、豊見城市、南風原町）を同時に表示してみた。 -- [[okinawa]] &new{2008-11-17 (月) 17:15:05};
- getCDATA関数復活。これでPlaceMark名にCDATA内の指定データをnameタグに挿入できます。これをsetAltitude()と一緒に使うことで、placemarkに日本語が表示できるようになると思います。 （Macでは成功、Windowsではうまくいかないなあ・・・）-- [[okinawa]] &new{2008-11-16 (日) 13:50:16};
- MacOSXでrgdalを使えるようにするを追加。これでMacでもKML自由自在ですね。 -- [[okinawa]] &new{2008-11-15 (土) 15:26:01};
- 明日、那覇市長選がある。那覇市国勢調査編として分析したが、３Dグラフから那覇市はすでに人口が空洞化していると考えられる。このような地域は企業が集中している反面、取り残された住民に高齢者の割合が高くなる。いわゆる都市構造的な問題は、一自治体の政策だけでは解決できないと考えられるが、いかがだろうか？ -- [[okinawa]] &new{2008-11-15 (土) 12:41:58};
- kmlFunctions.rのsetAltitude()にラベルのパラメータ指定を追加。ラベルを指定することで、KML内のPlaceMark名をセットする。ただし、現状では日本語をセットするとエラーになる。XMLパッケージのロケール問題っぽい。 -- [[okinawa]] &new{2008-11-14 (金) 14:53:02};
- kmlFunctions.rバグ取り＆ちょっと修正。 -- [[okinawa]] &new{2008-11-12 (水) 15:30:15};
- データを読み込んでkmlを再構成する（２） 那覇市国勢調査編追加 -- [[okinawa]] &new{2008-11-07 (金) 14:49:09};
- kmlAltitudeNode関数修正、setAltitudeMesh関数をsetAltitude関数に統合。MGパラメータでメッシュKMLと切り替えます。 -- [[okinawa]] &new{2008-11-07 (金) 11:56:58};
- 次はメッシュデータを自動でつくってみます。 -- [[okinawa]] &new{2008-11-06 (木) 15:20:27};
- データを読み込んで、kmlを再構成するを追加。shapefileから直接つくるのにくらべ簡素に？なりました。 -- [[okinawa]] &new{2008-11-06 (木) 14:50:28};
- kmlFunctions.r追加。 -- [[okinawa]] &new{2008-11-06 (木) 11:06:04};
- colors.kml関数追加しました。 -- [[okinawa]] &new{2008-11-05 (水) 16:19:15};
- getKmlNode関数追加。ノードを指定してKMLの切り抜きができます。 -- [[okinawa]] &new{2008-10-24 (金) 14:47:32};
- 関数追加。あとは色の変換関数がいるな。 -- [[okinawa]] &new{2008-10-23 (木) 21:50:06};
- 関数を若干修正 -- [[okinawa]] &new{2008-10-23 (木) 08:56:03};
- 地域メッシュKMLデータの応用追加 -- [[okinawa]] &new{2008-10-22 (水) 18:22:05};
- kmlファイルをRから直接いじってみました。今後は、ライブラリを充実させ、色々な機能追加をしてみようと思います。 -- [[okinawa]] &new{2008-10-21 (火) 18:03:29};
*kmlデータをXMLパッケージで再構成する関数 [#p7c92b92]
writeOGRで出力されたkmlデータをXMLパッケージを用いて再構成し、高さ、ポリゴンの色、ポリゴンの線の色を指定できる関数を作成しました。~
(構造決めうちなので、[[ShapeFileライブラリ]]のデータをwriteOGRでkml出力したデータしか使えないと思います）~
***kmlFunctions.r [#rb911ac6]
kmlを操作する関数をまとめました。下記のように、適当な場所に保存して、souceで呼び出してください。（GE5向けに新しくしました）~
#ref(kmlラボ/kmlFunctions.r)~
古いものは↓
#ref(kmlライブラリ/kmlFunctions.r)~
library(XML) #宣言
source("c:/kmlFunctions.r")#ソースの読み込み
doc <- xmlTreeParse("c:/gunma.kml")#土台となるkmlファイル読み込み
data<-sample(1000:10000,39) #高度データを自動生成
fillcolor<-colors.kml(heat.colors(39),"EE")#色の割り当て（ポリゴン）
linecolor<-NA#色の割り当て（線）
x<-kmlAltitudeNode(doc,data,fillcolor,linecolor)#kmlファイルにデータを書き込み
saveXML(x,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')#ファイルに出力
~
#ref(kmlライブラリ/color2.JPG,30%)~
***kmlAltitudeNode関数(GE5) [#n3382d3b]
kmlAltitudeNode<-function(doc,data,fillcolor,linecolor,extrude="1",tessellate="1",altitudeMode="absolute",MG=T){
kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
kmlver5<-'http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2'
kmlnode<-xmlNode("kml",node0,attrs=c(xmlns=kmlver5))
len1<-length(node0[[1]])
for(i in 2:len1){
if(MG==T){
node1<-node0[[1]][[i]][[2]][[1]]
}else if(MG==F){
node1<-node0[[1]][[i]][[2]]
}
node1<-addChildren(node1,xmlNode("extrude",extrude),xmlNode("tessellate",tessellate),xmlNode("gx:altitudeMode",altitudeMode))
node2<-node0[[1]][[i]][[3]][[2]]
node2<-xmlNode("PolyStyle",xmlNode("color",fillcolor[i-1]))
#linecolor
node3<-node0[[1]][[i]][[3]][[1]]
node3<-xmlNode("LineStyle",xmlNode("color",linecolor))
#polygon
len2<-length(node0[[1]][[i]][[2]][[1]])
for(j in 1:len2){
x1<-xmlValue(node1[[j]][[1]][[1]])
x2<-rev(unlist(strsplit(x1," ")))
x3<-paste(x2,data[[i-1]],sep=",",collapse=" ")
node1[[j]][[1]][[1]]<-xmlNode("coordinates",x3)
}
#replace node
if(MG==T){
node0[[1]][[i]][[2]][[1]]<-node1
}else if(MG==F){
node0[[1]][[i]][[2]]<-node1
}
node0[[1]][[i]][[3]][[2]]<-node2
node0[[1]][[i]][[3]][[1]]<-node3
kmlnode[[1]]<-node0
}
return(kmlnode)
}
#sample
library(XML)
doc <- xmlTreeParse("c:/gunma.kml")
data<-sample(1000:10000,39)
fillcolor<-c("ff000000","ffff0000","ff00ff00","ff0000ff","ffffff00","ffff00ff","ff00ffff","ffffffff","00000000")
linecolor<-"NA"
x<-kmlAltitudeNode(doc,data,fillcolor,linecolor)
saveXML(x,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/xml1.JPG,30%)~
#ref(kmlラボ/gunma2.kml)~
***setAltitude関数（GE5) [#nb393eab]
ノード番号を指定してピンポイントでAltitudeをセットする。~
setAltitude<-function(kmlnode,nodeno,data,fillcolor,linecolor,extrude="1",tessellate="1",altitudeMode="absolute",MG=T,label=""){
#kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
kmlver5<-'http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2'
kmlnode<-xmlNode("kml",node0,attrs=c(xmlns=kmlver5))
if(MG==T){
#altitude
node1<-node0[[1]][[nodeno]][[2]][[1]]
node1<-addChildren(node1,xmlNode("extrude",extrude),xmlNode("tessellate",tessellate),xmlNode("gx:altitudeMode",altitudeMode))
#fillcolor
node2<-node0[[1]][[nodeno]][[3]][[2]]
node2<-xmlNode("PolyStyle",xmlNode("color",fillcolor))
#linecolor
node3<-node0[[1]][[nodeno]][[3]][[1]]
node3<-xmlNode("LineStyle",xmlNode("color",linecolor))
#label
node4<-node0[[1]][[nodeno]]
node4<-addChildren(node4,xmlNode("name",label))
#polygon
len2<-length(node0[[1]][[nodeno]][[2]][[1]])
j<-1:len2
text1<-paste('x1<-xmlValue(node1[[',j,']][[1]][[1]]);x2<-rev(unlist(strsplit(x1," ")));',sep='')
text2<-paste('x3<-paste(x2,data,sep=",",collapse=" ");node1[[',j,']][[1]][[1]]<-xmlNode("coordinates",x3)',sep='')
text<-paste(text1,text2,sep="")
eval(parse(text=text))
node0[[1]][[nodeno]]<-node4
node0[[1]][[nodeno]][[2]][[1]]<-node1
node0[[1]][[nodeno]][[3]][[2]]<-node2
node0[[1]][[nodeno]][[3]][[1]]<-node3
kmlnode[[1]]<-node0
}else if(MG==F){
#altitude
node1<-node0[[1]][[nodeno]][[2]]
node1<-addChildren(node1,xmlNode("extrude",extrude),xmlNode("tessellate",tessellate),xmlNode("gx:altitudeMode",altitudeMode))
#fillcolor
node2<-node0[[1]][[nodeno]][[3]][[2]]
node2<-xmlNode("PolyStyle",xmlNode("color",fillcolor))
#linecolor
node3<-node0[[1]][[nodeno]][[3]][[1]]
node3<-xmlNode("LineStyle",xmlNode("color",linecolor))
#label
node4<-node0[[1]][[nodeno]]
node4<-addChildren(node4,xmlNode("name",label))
#polygon
len2<-length(node0[[1]][[nodeno]][[2]][[1]])
j<-1:len2
text1<-paste('x1<-xmlValue(node1[[',j,']][[1]][[1]]);x2<-rev(unlist(strsplit(x1," ")));',sep='')
text2<-paste('x3<-paste(x2,data,sep=",",collapse=" ");node1[[',j,']][[1]][[1]]<-xmlNode("coordinates",x3)',sep='')
text<-paste(text1,text2,sep="")
eval(parse(text=text))
node0[[1]][[nodeno]]<-node4
node0[[1]][[nodeno]][[2]]<-node1
node0[[1]][[nodeno]][[3]][[2]]<-node2
node0[[1]][[nodeno]][[3]][[1]]<-node3
kmlnode[[1]]<-node0
}
return(kmlnode)
}
getCDATANo関数と組み合わせて必要な場所のみ設定~
#sample
doc <- xmlTreeParse("c:/gunma.kml")
x<-xmlRoot(doc)
#1つめの地域
name<-"JCODE"
value<-"10366"
nodeno<-getCDATANo(x,name,value)
data<-1000
fillcolor<-c("ffff0000")
linecolor<-"ff000000"
x<-setAltitude(x,nodeno,data,fillcolor,linecolor)
#２つめの地域
name<-"JCODE"
value<-"10201"
nodeno<-getCDATANo(x,name,value)
data<-1000
fillcolor<-c("ffff00ff")
linecolor<-"ff000000"
x<-setAltitude(x,nodeno,data,fillcolor,linecolor)
#保存
saveXML(x,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/xml4.JPG,30%)~
メッシュのKMLはMultiGeometoryノードが無いので、MGパラメータでF指定する~
#sample
library(XML)
doc <- xmlTreeParse("c:/StdMesh01.kml")
kmlnode<-xmlRoot(doc)
#1
name<-"CODE"
value<-"3622"
nodeno<-getCDATANo(kmlnode,name,value)
data<-1000000
fillcolor<-c("ffff0000")
linecolor<-"ff000000"
kmlnode<-setAltitude(kmlnode,nodeno,data,fillcolor,linecolor,MG=F)
#2
name<-"CODE"
value<-"3623"
nodeno<-getCDATANo(kmlnode,name,value)
data<-800000
fillcolor<-c("ffff00ff")
linecolor<-"ff000000"
kmlnode<-setAltitude(kmlnode,nodeno,data,fillcolor,linecolor,MG=F)
#save
saveXML(kmlnode,"c:/StdMesh012.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/color4.JPG,30%)~
***kmlNodeLength関数 [#y7fb430a]
kmlAltitudeNode関数を使うときに、kmlのノード数を知りたい時に使う。~
#kmlNodeLength() ver 0.1.2
kmlNodeLength<-function(doc){
library(XML)
kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
len1<-length(node0[[1]])
return(len1)
}
#sample
library(XML)
doc <- xmlTreeParse("c:/gunma.kml")
kmlNodeLength(doc)
> #sample
> library(XML)
> doc <- xmlTreeParse("c:/gunma.kml")
> kmlNodeLength(doc)
[1] 40
注）先頭のノードに<name>gunma</name>が入っていますので、市町村数はノード数−１となります。~
***getCDATAList関数 [#ic0a4023]
CDATA(HTML)の中身を取り出す。nameに対応する項目名を指定する。（必須）~
getCDATAList<-function(doc,name){
getCDATA<-function(node,name){
x<-gsub(" ","",iconv(xmlValue(node),"utf-8","cp932"))
x<-unlist(strsplit(x,"<br/>\n"))
x<-x[[grep(name,x)]]
x<-gsub("<b>","",gsub("</b>","",gsub("<i>","",gsub("</i>","",gsub("<br/>","",gsub(paste(name,":",sep=""),"",x))))))
#x<-iconv(x,"","utf-8")
return(x)
}
kmlnode <- xmlRoot(doc)
node0<-kmlnode[[1]]
len1<-length(node0[[1]])
x<-NULL
for(i in 2:len1){
x[i-1]<-getCDATA(node0[[1]][[i]][[1]][[1]],name)
}
return(x)
}
#sample
doc <- xmlTreeParse("c:/gunma.kml")
name<-"CITY1"
getCDATAList(doc,name)
> #sample
> doc <- xmlTreeParse("c:/gunma.kml")
> name<-"CITY1"
> getCDATAList(doc,name)
[1] "前橋市" "高崎市" "桐生市" "伊勢崎市" "太田市" "沼田市" "館林市" "渋川市" "藤岡市" "富岡市"
[11] "安中市" "みどり市" "富士見村" "榛名町" "榛東村" "吉岡町" "吉井町" "上野村" "神流町" "下仁田町"
[21] "南牧村" "甘楽町" "中之条町" "長野原町" "嬬恋村" "草津町" "六合村" "高山村" "東吾妻町" "片品村"
[31] "川場村" "昭和村" "みなかみ町" "玉村町" "板倉町" "明和町" "千代田町" "大泉町" "邑楽町"
> #sample
> doc <- xmlTreeParse("c:/gunma.kml")
> name<-"JCODE"
> getCDATAList(doc,name)
[1] "10201" "10202" "10203" "10204" "10205" "10206" "10207" "10208" "10209" "10210" "10211" "10212" "10303" "10321" "10344" "10345" "10363"
[18] "10366" "10367" "10382" "10383" "10384" "10421" "10424" "10425" "10426" "10427" "10428" "10429" "10443" "10444" "10448" "10449" "10464"
[35] "10521" "10522" "10523" "10524" "10525"
***getCDATA関数 [#vf724883]
CDATAのデータを取り出す。cdatanoは、grepで検索した場合、複数の位置が取り出された場合、初期位置からの相対位置を指定するパラメータ。~
getCDATA<-function(kmlnode,name,nodeno=1,sep=":",cdatano=1){
node0<-kmlnode[[1]]
node1<-node0[[1]][[nodeno]][[1]][[1]]
#x<-gsub(" ","",xmlValue(node1))#mac
x<-gsub(" ","",iconv(xmlValue(node1),"utf-8","cp932"))
#print(x)
x<-unlist(strsplit(x,"<br/>\n"))
xno<-grep(name,x)
len<-length(xno)
if(len>1){
x<-x[[xno[cdatano]]]
}else if(len==1){
x<-x[[xno]]
}
#print(x)
x<-gsub("<b>","",gsub("</b>","",gsub("<i>","",gsub("</i>","",gsub("<br/>","",x)))))
x<-gsub(paste(name,sep,sep=""),"",x)
return(x)
}
> #sample
> doc <- xmlTreeParse("/gunma.kml")
> kmlnode<-xmlRoot(doc)
> name<-"JCODE"
> nodeno<-2
> getCDATA(kmlnode,name,nodeno)
[1] "10201"
***getCDATANo関数 [#x4aa4dac]
nameとvalueを指定することで、CDATA(HTML)内のデータを検索し、一致したノード番号を返す。~
#getCDATANo
getCDATANo<-function(kmlnode,name,value){
getCDATA2<-function(node,name){
x<-gsub(" ","",iconv(xmlValue(node),"utf-8","cp932"))
x<-unlist(strsplit(x,"<br/>\n"))
x<-x[[grep(name,x)]]
x<-gsub("<b>","",gsub("</b>","",gsub("<i>","",gsub("</i>","",gsub("<br/>","",x)))))
#x<-iconv(x,"","utf-8")
return(x)
}
#kmlnode <- xmlRoot(doc)
namevalue<-paste(name,":",value,sep="")
node0<-kmlnode[[1]]
len1<-length(node0[[1]])
for(i in 2:len1){
if(namevalue==getCDATA2(node0[[1]][[i]][[1]][[1]],name)){
x<-i
break
}else{
x<-NA
}
}
return(x)
}
#sample
doc <- xmlTreeParse("c:/gunma.kml")
kmlnode<-xmlRoot(doc)
name<-"JCODE"
value<-"10366"
getCDATANo(kmlnode,name,value)
> #sample
> doc <- xmlTreeParse("c:/gunma.kml")
> name<-"JCODE"
> value<-"10366"
> getCDATANo(kmlnode,name,value)
[1] 19
***getKmlNode関数(GE5) [#udad8667]
ノードをnodenoで指定し、指定したノードの部分だけ抽出してKMLを再構成~
#getKmlNode
getKmlNode<-function(kmlnode,nodeno){
len<-length(nodeno)
node0<-kmlnode[[1]]
kmlver5<-'http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2'
kmlnode<-xmlNode("kml",node0,attrs=c(xmlns=kmlver5))
node2<-xmlNode("Folder",xmlNode("name",as.character(nodeno)))
for(i in 1:len){
node1<-node0[[1]][[nodeno[[i]]]]
node2<-addChildren(node2,node1)
}
kmlnode[[1]][[1]]<-node2
return(kmlnode)
}
#sample
library(XML)
doc <- xmlTreeParse("c:/gunma.kml")
kmlnode<-xmlRoot(doc)
nodeno<-c(2,4,6)
expdata<-getKmlNode(kmlnode,nodeno)
saveXML(expdata,"c:/gunma3.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
***colors.kml関数 [#r236dd3e]
colors.kml("色名","アルファチャンネル")
colors.kml<-function(rgb=colors(),alpha="FF"){
colors.hex <- function( x=colors() ) {
color.hex <- function(x) do.call( "rgb", as.list(col2rgb(x)/255) )
sapply( x, color.hex )
}
col.hex<-colors.hex(rgb)
red<-substr(col.hex,2,3)
green<-substr(col.hex,4,5)
blue<-substr(col.hex,6,7)
alpha<-alpha
col.kml<-paste(alpha,blue,green,red,sep="")
return(col.kml)
}
> colors.kml("red")
[1] "FF0000FF"
> colors.kml(rainbow(10))
[1] "FF0000FF" "FF0099FF" "FF00FFCC" "FF00FF33" "FF66FF00" "FFFFFF00" "FFFF6600" "FFFF0033" "FFFF00CC"
[10] "FF9900FF"
参考：[[色見本]]~
*メッシュkmlデータの作成 [#t6141a60]
メッシュkmlデータを自動作成します。生成されたkmlはもちろん他の関数で操作できます。
**makeMesh1関数：第１次地域区画（１次メッシュ）生成関数 [#f324b924]
指定された地域メッシュコードを基に１次メッシュをブロック生成します。~
下記のDefaultMesh1.kmlをタネにして作成しますので、ダウンロードして利用してください。~
#ref(kmlラボ/DefaultMesh1.kml)~
makeMesh1<-function(kmlnode,codea,codeb){
transMeshCode<-function(code){
lat0<-as.integer(substr(code,1,2))/1.5
lon0<-as.integer(substr(code,3,4))+100
lat1<-lat0;lon1<-lon0+1;lat2<-lat0+(2/3);lon2<-lon1
lat3<-lat2;lon3<-lon0
x<-paste(lon0,",",lat0," ",lon1,",",lat1," ",lon2,",",lat2," ",lon3,",",lat3," ",lon0,",",lat0,sep="")
return(x)
}
code1<-as.integer(substr(codea,1,2))
code3<-as.integer(substr(codea,3,4))
code2<-as.integer(substr(codeb,1,2))
code4<-as.integer(substr(codeb,3,4))
print(code1)
print(code2)
print(code3)
print(code4)
node0<-kmlnode[[1]]
node1<-node0[[1]]
node2<-node0[[1]][[2]]
x<<-node0[[1]]
nodex<-removeChildren(x,kids=c("Placemark"))
for(j in code1:code2){
code01<-paste(as.character(j),code3,sep="")
code02<-paste(as.character(j),code4,sep="")
print(j)
for(i in code01:code02){
print(i)
node3<-node0[[1]][[2]][[1]][[1]]
node4<-node0[[1]][[2]][[2]][[1]][[1]][[1]]
CDATA<-paste("<b>CODE:</b><i>",i,"</i><br />\n",sep="")
node3<-xmlCDataNode(CDATA)
node4<-xmlNode("coordinates",transMeshCode(i))
node0[[1]][[2]][[2]][[1]][[1]][[1]]<-node4
node0[[1]][[2]][[1]][[1]]<-node3
node2<-node0[[1]][[2]]
nodex<-addChildren(nodex,node2)
}
node0[[1]]<-nodex
}
kmlnode[[1]]<-node0
return(kmlnode)
}
#Sample
library(XML)
doc <- xmlTreeParse("c:/DefaultMesh1.kml")
kmlnode <- xmlRoot(doc)
kmlnode2<-makeMesh1(kmlnode,3622,3643)#2つの区画を指定することで、２区画間の区画をブロック生成する
saveXML(kmlnode2,"c:/meshdata.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlラボ/mesh2.JPG,30%)~
１次メッシュ生成~
#Sample2:少々時間がかかります
library(XML)
doc <- xmlTreeParse("c:/DefaultMesh1.kml")
kmlnode <- xmlRoot(doc)
kmlnode2<-makeMesh1(kmlnode,3622,6853)
saveXML(kmlnode2,"c:/meshtest.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlラボ/mesh4.JPG,30%)~
*応用編 [#me6a710e]
***データを読み込んでkmlを再構成する（１） [#cce8c1cb]
csvデータを読み込んで、kmlを再構成しています。JCODEでマッチングさせていますので、事前にデータをノード順と一致させる必要はありません。~
#ref(kmlラボ/agedprop.csv)~
library(XML)
source("c:/kmlFunctions.r")
#kmlファイルの読み込み
doc <- xmlTreeParse("c:/gunma.kml")
kmlnode<-xmlRoot(doc)
#csvデータの取得
dat<-read.csv("c:/agedprop.csv")
data<-dat$AP2006
Altdata<-data*1000#高度
name<-"JCODE"#マッチング用指定
#色階層の割り当て
cutnum<-16
classes <- cut(data,seq(min(data),max(data),length=cutnum+1),include.lowest=T)
cols <- colors.kml(rev(heat.colors(cutnum)))
cols2<-cols[classes]
#kmlへの書き込み
len<-kmlNodeLength(doc)-1
for(i in 1:len){
value<-dat$JCODE[[i]]
print(value)
nodeno<-getCDATANo(kmlnode,name,value)#JCODEからノードを取得
print(nodeno)
fillcolor<-cols2[[i]]
linecolor<-NA
kmlnode<-setAltitude(kmlnode,nodeno,Altdata[[i]],fillcolor,linecolor)
}
#kml保存
saveXML(kmlnode,"c:/gunma2.kml",prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
#ref(kmlライブラリ/color3.JPG,30%)~
***データを読み込んでkmlを再構成する（２） 那覇市国勢調査編 [#e7808495]
まず、那覇市国勢調査のデータ（平成１７年）を地図で見る統計（統計ＧＩＳ）のデータダウンロードから落とします。~
http://www.e-stat.go.jp/SG2/toukeichiri/TopFrame.do?fromPage=init&toPage=download~
（１）ダウンロードページ
#ref(kmlライブラリ/naha1.JPG,30%)
（２）世界測地系緯度経度Shape形式を落とします。
#ref(kmlライブラリ/naha2.JPG,30%)
（３）前準備~
落としたShapeFileから、writeOGRでKMLデータ（h17ka47201.kml)を作成します。~
library(maptools)
library(rgdal)
x<-readShapePoly("c:/nahashi/h17ka47201.shp")
x@data$KEN_NAME<-iconv(x@data$KEN_NAME,"","utf-8")
x@data$GST_NAME<-iconv(x@data$GST_NAME,"","utf-8")
x@data$MOJI<-iconv(x@data$MOJI,"","utf-8")
tf <- "c:/nahashi/"
writeOGR(x, paste(tf, "h17ka47201.kml", sep=""), "h17ka47201", "KML")
（４）3Dグラフの作成~
作成したkmlとShapeFileに付属しているdbfを使って3Dグラフを作成します。~
国勢調査データから作成したkmlはMultiGeometryタグが無いので、setAltitude関数をMG=Fで取り込みます。~
ここで、表示させるのは人口のデータです。~
library(XML)
library(foreign)
#path<-"/"#Mac
path<-"c:/"
source(paste(path,"kmlFunctions.r",sep=""))
doc <- xmlTreeParse(paste(path,"nahashi/h17ka47201.kml",sep=""))
kmlnode<-xmlRoot(doc)
#データの取得
dat<-read.dbf(paste(path,"nahashi/h17ka47201.dbf",sep=""))
data<-dat$JINKO#人口データ指定(ここを変えることで他の項目も作れます）
Altdata<-data+10#高度:ポリゴンホールが出来ないように全体を１０m上げる
name<-"SEQ_NO2"#マッチング用指定
#色階層の割り当て
cutnum<-16
classes <- cut(data,seq(min(data),max(data),length=cutnum+1),include.lowest=T)
cols <- colors.kml(rev(heat.colors(cutnum)))
cols2<-cols[classes]
#kmlへの書き込み
len<-kmlNodeLength(doc)-1
for(i in 1:len){
value<-dat$SEQ_NO2[[i]]#マッチング用指定
print(value)
nodeno<-getCDATANo(kmlnode,name,value)#SEQ_NO2からノードを取得
#nodeno<-i+1 #こちらを使うとかなり速いです（ノードとＤＢＦのデータの並びが同じ場合のみ。SEQ_NO2でマッチングさせてないので速いです。）
print(nodeno)
fillcolor<-cols2[[i]]
linecolor<-NA
#Macの場合以下のコードで実行可能
#label<-getCDATA(kmlnode,"MOJI",nodeno,cdatano=2)#PlaceMark名にMOJIを指定。NMOJIも検索してくるので、相対位置パラメータcdatano=2を指定。
#Windowsの場合、いまのところlabelに日本語指定不可
label<-getCDATA(kmlnode,"KEYCODE1",nodeno)
kmlnode<-setAltitude(kmlnode,nodeno,Altdata[[i]],fillcolor,linecolor,MG=F,label=label)
}
#kml保存
saveXML(kmlnode,paste(path,"nahashi/nahashi3.kml",sep=""),prefix = '<?xml version="1.0" encoding="UTF-8"?>\n')
作成には、かなりの時間がかかりますので、気長にお待ちください。（私のＰＣで１０分）~
ただし、nodeno<-i+1の方を使うとすぐに作成できます。（１分くらい）ただし、ノードとＤＢＦのデータの並びが同じ場合のみです。~
データは、nahashi3.kmlで落ちてきます。~
赤くて高い場所ほど人口が多い地域です。那覇港は手作業で着色しました。~
#ref(kmlラボ/naha5.JPG,30%)
那覇市の周辺市町村（浦添市、豊見城市、南風原町）を同時に表示してみた。~
明らかに周辺市町村に人口が流れているのがわかる。~
#ref(kmlラボ/city2.JPG,30%)
なお、国勢調査データには利用の制限がありますので、ご注意ください。~
~
利用の制限~
~
１． 利用者は、本システムでダウンロードしたデータ及び画像データをそのまま複製（ファイル形式を変換しての複製を含む。）して第三者に譲渡することを禁じます。
なお、商用目的で複製する場合は、e-Statの「お問い合わせ」よりご連絡ください。
~
２．法律、政令、規則、省令その他すべての法令及び条例等の法規に違反する目的、手段又は方法で、本システムでダウンロードしたデータ及び画像データを利用することを一切禁じます。また、他人の権利を侵害する目的、手段又は方法での利用、公序良俗に反する利用についても一切禁じます。
*MacOSXでrgdalを使えるようにする [#l1ecb7c3]
rgdalパッケージはMacOSXではバイナリで提供されていませんが、下記のサイトを参考にすると、rgdalをMacOSXに導入することができます。（#GDALのインストールと#rgdalのインストールだけやってください）~
//↓の記事書いたものです．BLOGを移転しましたので，新しいアドレスに変更させていただきました．移転前のは近いうちに消すと思います．
//何か問題があればお手数をおかけしてしまい申し訳ありませんが，修正いただければ幸いです．何卒宜しくお願い申し上げます．
//http://www.0093-style.net/0093/2008/02/mac_os_x_leopar_31f8.html~
http://blog.0093.tv/2008/02/mac-os-x-leopardgrassgdalrgdalspgrass6_4571.html~
writeOGRもきちんと動きました。(^_^)/~
** MacOS X でrgdalを使えるようにするもうひとつの方法 [#l9a7d706]
+[[http://www.kyngchaos.com/software:frameworks]]より、GDAL completeをインストール
+ /Library/Frameworks/GDAL.framework/Programs/gdal-config を /usr/local/bin などにシンボリックリンク
+ rgdalのソースパッケージをCRANから取得
+ $ R32 CMD INSTALL rgdal_0.6-20.tar.gz --configure-args='--with-proj-include=/Library/Frameworks/PROJ.framework/Headers --with-proj-lib=/Library/Frameworks/PROJ.Framework/unix/lib'
*アクセス数： [#tdc929b1]
&counter;
人
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
