ArcGIS で R

ArcGIS で R

VBA での利用方法†

R(D)COM の ArcGIS 内でのセットアップ†

１．R(D)COMのインストール。

２．ArcMapの起動

３．ArcMapのVisual Basic Editor の起動

４．「Tools」- 「Refrences」を開く

５．「References」内の StatConnector を選ぶ

サンプルプログラム１†

　１．VBA フォームでボタンを作成し、そこに以下のスクリプトを埋め込む。

Private Sub CommandButton1_Click()
 Dim RLink As StatConnector
 Dim mean As Double
 Dim n As Long
 
 n = 150
 Set RLink = New StatConnector
 RLink.Init ("R")'R の初期化
 RLink.SetSymbol "n", n
 RLink.EvaluateNoReturn ("x <- rnorm(n,5,2)")
 ' mean = RLink.GetSymbol("mean(x)") Incorrect: mean(x) is not a symbol
 RLink.EvaluateNoReturn ("m <- mean(x)")
 mean = RLink.GetSymbol("m")
 RLink.Close
 MsgBox "平均値： " & mean, vbOKOnly, "Powered by R"
End Sub

　２．ボタンを実行する。

サンプルプログラム２†

　画像の利用例。

サンプルプログラム３†

　コロプレス図

Python からの利用方法†

　Python からArcGISとRのそれぞれの機能を呼び出すことができ、両者を組み合わせて、空間回帰などができる。

• RPyGeo パッケージ:Python 経由で、R から、ArcGIS の処理を呼び出す
• Rpy
• RScript

リンク†

　R / Arcgis Repository R(D)COM のStatConnector とVBAを使った例がある。仕掛け人はあの Bivand先生（６月頃にも来日予定。ユーザ会にも参加？）

利用例†

• 佐藤俊明・岡部篤行(2006):ネットワークボロノイクロスＫ関数法の提案とそのツール開発, CSIS Discussion paper，82.

参考†

• ArcGIS 9.3 beta の　Spatial Analysis Extenstion の空間統計分析機能 これで、ArcGIS ユーザ（特に日本の）は R 言語の空間統計機能は使わなくなるかな？
• ArcGIS 9.3 の help の地理的加重回帰の説明
• RPyGeo(Python 経由による、R での ArcGIS のジオプロセッシング)パッケージ中のオブジェクト中の一覧

最新の20件