FeatureClassToNumPyArray


NumPyからラスタ
NumPy⇒Raster @ArcPy

テーブルからNumPy
http://resources.arcgis.com/ja/help/main/10.1/index.html#//018w00000015000000

まあGDALでもできるが。

SciPy標準で入れてくれれば、主成分とかクラスタとか楽なのだが。。

取り敢えず順列とって相関係数確認コード
つかえるカラム全てぶちこめば相関行列がとれるので無駄コードだが。

# -*- coding: cp932 -*-

import arcpy
import itertools
import numpy

fclass = r"C:\xxx.gdb\pref"

fieldList = arcpy.ListFields(fclass)

lstFname = []
for field in fieldList:
    if field.type != "Double" and field.type != "Single" and field.type != "Integer" and field.type != "SmallInteger":
        continue
    lstFname.append(field.name)
    
for comb in itertools.combinations(lstFname, 2):
    print(comb)
    arr = arcpy.da.FeatureClassToNumPyArray(fclass, comb)
    print(numpy.corrcoef((arr[comb[0]], arr[comb[1]])))
カテゴリー: 開発 タグ: パーマリンク