YouTube | Facebook | Twitter  Feed

ラスターのジオリファレンス(多項式変換)残差・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]);
}

関連記事

  • この記事を書いた人

羽田 康祐

Esri認定インストラクター、GIS上級技術者、測量士補、潜水士。GISy / GISc とその関連分野である地理学・地図学について日々の出来事で学んだ記憶を記録するためにブログを書いています。行動原理は伊達と酔狂。好きな地形は圏谷。好きな地図投影法はパースクインカンシャル図法。呉市生まれ広島市出身。GIS を使った自己紹介はこちら

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

© 2020 WINGFIELD since1981