YouTube | Facebook | X(Twitter) | RSS

ラスターのジオリファレンス(多項式変換)残差・RMS エラーの取得

2016/9/1 (木)

//元座標値
IPoint[] pPointSource = new PointClass[4];
pPointSource[0] = new PointClass();
pPointSource[0].PutCoords(13411566.749205, 392094.84965);
 
//補正座標値
IPoint[] pPointTarget = new PointClass[4];
pPointSource[0] = new PointClass();
pPointTarget[0].PutCoords(13411793.050983, 392030.191999);
 
//元座標と補正座標でPointCollectionを作成
IPointCollection pPointCollectionSource = new MultipointClass();
pPointCollectionSource.AddPoint(pPointSource[0]);
IPointCollection pPointCollectionTarget = new MultipointClass();
pPointCollectionTarget.AddPoint(pPointTarget[0]);
 
object pResidual;
double dblRMS;
IPolynomialXform pPolynomialXform = new PolynomialXformClass();
pPolynomialXform.DefineFromControlPoints(pPointCollectionSource, pPointCollectio    //第3引数で何次多項式変換かを指定(1-3)
pPolynomialXform.CheckResidualRMS(pPointCollectionSource, pPointCollectionTargetormDirection.esriTransformForward,out pResidual,out dblRMS);
 
double[,] dblArray = (double[,])pResidual;  //二次元配列にキャスト
 
//配列から残差を取得(1次元目の要素がレコード、2次元目の配列が順に 残差_x、残差_y, 残差)
for(long i= 0;i == dblArray.GetLength(0);++i)
{
    Console.WriteLine(dblArray[i, 0] +"," + dblArray[i, 1] + "," + dblArray[i, 2]);
}
  • この記事を書いた人

羽田 康祐

伊達と酔狂のGISエンジニア。GIS上級技術者、Esri認定インストラクター、CompTIA CTT+ Classroom Trainer、潜水士、PADIダイブマスター、四アマ。WordPress は 2.1 からのユーザーで歴だけは長い。 代表著書『"地図リテラシー入門―地図の正しい読み方・描き方がわかる』 GIS を使った自己紹介はこちら。ESRIジャパン(株)所属、元青山学院大学非常勤講師を兼務。日本地図学会第31期常任委員。発言は個人の見解です。

-プログラミング, ArcGIS
-,