JavaScript8

１３．３　開放トラバース 　　　　　　　　　[もくじへ] [トップページへ] [前ページへ]

===================================================

プログラム例題２４

開放トラバースのプログラムを作成せよ。

=========================================================

考え方

初点方向角α_A、初点座標 A（X_A,Y_A）、点Aでの測角 β₁、測距S₁ とすると、

点Aより点１への方向角α₁は、

　　　α₁ = α_A + β₁ 　　　　（ α₁ が３６０°を超えたら３６０°を減ずる　）

点１の座標は、

　　　X₁ = X_A + S₁(α₁)

　　　Y₁ = Y_A + S₁(α₁)

となります。

また、点1より点2への方向角α₂は、

　　　α₂ = α₁ + β₂ + 180°（ α₂ が３６０°を超えたら３６０°を減ずる　）

方向角α₂が求められれば、点2の座標は点1の座標と同様に、

　　　X₂ = X₁ + S₂(α₂)

　　　Y₂ = Y₁ + S₂(α₂)

となります。

以下同様にして各点の座標を求めることができます。したがって、一般式は

次のようになります。

　　　α_n = α_n-1 + β_n + 180°（ α_n が３６０°を超えたら３６０°を減ずる　）

　　　X_n = X_n-1 + S_n(α_n)

　　　Y_n = Y_n-1 + S_n(α_n)

プログラム

<HTML>
<HEAD>
　　　<TITLE>例題２４</TITLE>

<SCRIPT LANGUAGE="JavaScript">

</SCRIPT>

</HEAD>

<BODY>

<FORM NAME="Form1">
　＊＊＊　トラバース計算　＊＊＊<BR>
<BR>
初点名= <INPUT TYPE="text" NAME="AN" size=15>
　　 XA= <INPUT TYPE="text" NAME="XA" size=15>
　　 YA= <INPUT TYPE="text" NAME="YA" size=15>
　　αA= <INPUT TYPE="text" NAME="AA" size=15>
<BR>
終点名= <INPUT TYPE="text" NAME="BN" size=15>
　　 XB= <INPUT TYPE="text" NAME="XB" size=15>
　　 YB= <INPUT TYPE="text" NAME="YB" size=15>
　　αB= <INPUT TYPE="text" NAME="AB" size=15>
<BR>
　番号= <INPUT TYPE="text" NAME="NO" size=15>
　点名= <INPUT TYPE="text" NAME="NA" size=15>
　測角= <INPUT TYPE="text" NAME="BT" size=15>
　距離= <INPUT TYPE="text" NAME="DS" size=15>
<BR>
<BR>
　　　　　 <INPUT TYPE="button" NAME="BCK" VALUE="　＜　" onClick="onBCK( )">
<INPUT TYPE="button" NAME="FRD" VALUE="　＞　" onClick="onFRD( )">
　　　　　 <INPUT TYPE="button" NAME="INI" VALUE="クリア" onClick="onINI( )">
<INPUT TYPE="button" NAME="CAL" VALUE="計　算" onClick="onCAL( )">
<HR>

</FORM>

</BODY>
</HTML>

入力データと実行結果

　　　　　　　　X_A=-6132.14, Y_A=8462.15

閉合トラバースや、結合トラバースにも使える表にしました。空白になっている

箇所があるのは、そのためです。

プログラムの説明

        ai = new Array(50);                                // 方向角
        x = new Array(50);                                 // X座標
        y = new Array(50);                                // Y座標
        dx = new Array(50);                                 // X座標差
        dy = new Array(50);                                // Y座標差
        var doo, fun, byo;                                    // 度分秒

入力画面のプログラムに、追加したグローバル変数です。

//　　角度の変換(度を秒に）
//
        function Dobyo( x ) {
                var d, f, b, r;
                d = Math.floor(x);
                f = Math.floor((x - d)*100);
                b = Math.floor(((x - d)*100 - f)*100 + 0.5);
                r = d*3600 + f*60 + b;
                return r;
        }
//
//
//　　角度の変換(秒を度分秒に）
//
        function Byodfb( x ) {
                doo = Math.floor(x/3600);
                fun = Math.floor((x - doo*3600)/60);
                byo = x - doo*3600 - fun*60;
                if (byo<0) byo = 0;
        }

測量計算では、角度は1秒も無視できません。ラジアンや度のままだと、場合

によっては数秒の誤差が出る場合があります。したがって、秒に変換して計

算誤差が出ないようにしています。これらは、そのための関数プロシージャ

です。

//　　四捨五入
//
        function Rund(x,y) {
                var sg, r;
                if (x>=0) {sg = 1;}
                else {sg =-1;}
                r = sg*Math.floor(Math.abs(x)*Math.pow(10,y) + 0.5)/Math.pow(10,y);
                return r;
        }

変数 sg は、四捨五入するデータの符号（±１）を代入します。abs( ) は、絶

対値を求める関数です。したがって、このプロシージャでは、データの絶対値

を小数第y+1位を四捨五入した後、符号を元に戻します。JavaScript には、

四捨五入の関数があります。

r = Math.round(x*Math.pow(10,y))/Math.pow(10,y) でこのプログラムと

同じ結果が得られます。

//　　　＊＊＊　計算部　＊＊＊
// 　　角度を秒に変換
                alpa = Dobyo(alpa);　　①
                for (i = 1; i<=n; i++) {
                        bt[i] = Dobyo(bt[i]);　　②
                }
//
                ai[1] = alpa + bt[1];　　③
                while (ai[1]>=1296000) {　　④
                           ai[1] = ai[1] - 1296000;
                }
                for (i = 2; i<=n; i++) {
                        ai[i] = ai[i-1] + bt[i] + 648000;　　⑤
                        while (ai[i]>=1296000) {
                                   ai[i] = ai[i] - 1296000;
                        }
                }
//
                dx[1] = Rund(ds[1]*Math.cos(ai[1]/648000*Math.PI),3);　　⑥
                dy[1] = Rund(ds[1]*Math.sin(ai[1]/648000*Math.PI),3);　　⑦
                x[1] = xa + dx[1];　　⑧
                y[1] = ya + dy[1];　　⑨
                for (i = 2; i<=n; i++) {
                        dx[i] = Rund(ds[i]*Math.cos(ai[i]/648000*Math.PI),3);　　⑩
                        dy[i] = Rund(ds[i]*Math.sin(ai[i]/648000*Math.PI),3);　　⑪
                        x[i] = x[i-1] + dx[i];　　⑫
                        y[i] = y[i-1] + dy[i];　　⑬
                }
                Kaiho( );　　⑭

開放トラバースの計算部です。

①②　入力された角度データをすべて秒に変換します。

③　　　初点から測点1への方向角の計算です。

④　　　③の方向角が３６０°（129600秒）を超えたら、３６０°を減じます。

while 文

ｆor 文はくり返し回数の決まっている場合に使いますが、 while 文はくり返し

回数の決まっていない場合に用います。

　　　　while ( 条件式 ) {

　　　　　　くりかえす文

　　　　}

のように記述し、条件式を満たしている間くり返します。

⑤　　　③以外の方向角の計算です。測点2以降の方向角計算は同じなので、

ループに入れます。

⑥⑦　初点からの座標差を求めます。座標差は小数第3位に丸めます。

⑧⑨　測点1の座標を求めます。

⑩～⑬　測点2以降の座標差と座標を求めます。　

⑭　開放トラバース計算結果出力のプロシージャを実行します。

// 開放トラバース計算結果出力
//
function Kaiho( ) {
　　・
　　・
　　・
｝

開放トラバース計算結果出力のプロシージャです。閉合トラバースと結合トラ

バースにも使えるように設計してあります。

JavaScript には、表示桁をそろえるフォーマット文がないので、HTML の表作

成のタグを使って表を作成します。しかし、表作成はなかなか面倒でなもので

す。このプログラムを作る時も、何回か試行錯誤をくりかえして、なんとか表ら

しいものにすることができました。

したがって、ここでは表作成のタグの説明をしますので、それを理解した上で

プログラムを解析し、すっきりした表を作ってみてください。

表作成のタグについて

表の始まりと終わりのタグは、<TABLE> </TABLE> です。その中に行を作

るには <TR> </TR> 、列を作るには <TD> </TD> のタグを使います。

また、HTML では行の中に列を作るという考え方をしますので、たとえ行の中

に列が１列しかなくても、<TR> <TD> </TD> </TR> のように使います。

表に枠を付ける時は <TABLE BORDER="1"> のようにします。BORDER

の数値は枠の太さを表します。

枠の中に入れるデータの位置を調整するには、<TD ALIGN="right"> のよ

うに、ALIGN を使います。左揃えは left、中央揃えは center、右揃えは

right を指定します。

たとえば、次のように HTML を書くと、

次のような、2行3列の表ができます。

BORDER などに値を指定する場合、"1" のようにダブルクォーテーション

で囲みますが、これを JavaScript で記述する場合、次のように内側はシン

グルクォーテーションにしなければなりません。

document.write("<CENTER><TABLE BORDER='1'>");

[もくじへ] [トップページへ]

次回は、『閉合トラバース』の予定です。（１０月１０日ごろ）

この講座に対する、ご意見やご要望はメールでお願いします。

『JavaScriptによる測量計算プログラミング入門』　

講師　大愚良寛

[もくじへ] [トップページへ] [次ページへ]