Page 21 - 中国仿真学会通讯2020第1期
P. 21
ｉｎ ｓｅｃｔｉｏｎ ４． Ｔｈｅ ｌａｓｔ ｓｅｃｔｉｏｎ ｄｒａｗｓ ｃｏｎｃｌｕｓｉｏｎｓ ｘ － ｘ１ｎ ＝ ｙ － ｙ ｎ ＝ ｚ － ｚ１ｎ ＝ ｋｎｇ ，ｋ ｎ ∈ ［０，１］ ．
ａｎｄ ｐｕｔｓ ｆｏｒｗａｒｄ ｅｘｐｅｃｔａｔｉｏｎｓ ｆｏｒ ｆｕｒｔｈｅｒ ｓｔｕｄｙ． ｘ２ｎ － ｘ１ｎ ｙ２ｎ － １ ｇ
ｙ１ｎ ｚ ｎ － ｚ１ｎ
２
２． Ｉｎｔｅｒｖｉｓｉｂｉｌｉｔｙ Ａｌｇｏｒｉｔｈｍ Ｆｏｒ ｔｈｅ ｔｅｒｒａｉｎ ｄａｔａ ｓｔｏｒｉｎｇ ｉｎ ＲＳＧ ｍｏｄｅｌ
ｗｈｏｓｅ ｉｎｔｅｒｖａｌ （ ｐｒｅｃｉｓｉｏｎ ） ｉｓ Δｘ ＝ Δｙ ＝ ａ，
Ｔｈｅ ｉｎｔｅｒｖｉｓｉｂｉｌｉｔｙ ｉｓ ａｃｑｕｉｒｅｄ ｂｙ ｃｏｍｐｕｔｉｎｇ ｉｎｔｅｒｖｉｓｉｂｉｌｉｔｙ ｃｏｍｐｕｔｉｎｇ ｉｎ ｏｎｅ ＧＣＳ ｔｉｌｅ ｉｓ ｂｒｏｋｅｎ
ｔｈｅ ｉｎｔｅｒｓｅｃｔｉｏｎｓ ｂｅｔｗｅｅｎ ａｌｌ ｏｆ ｔｈｅ ｓｅｇｍｅｎｔｓ ｉｎ ｕｐ ｉｎｔｏ ｔｈａｔ ｉｎ ｓｅｖｅｒａｌ ＲＳＧｓ． Ｔｈｅ ｉｎｔｅｒｓｅｃｔｉｏｎ
ｄｉｆｆｅｒｅｎｔ ＧＣＳ ｔｉｌｅｓ ａｎｄ ｔｈｅ ｔｅｒｒａｉｎ ｃｕｒｅｄ ｓｕｒｆａｃｅｓ ｐｏｉｎｔ ｂｅｔｗｅｅｎ ｔｈｅ ＬｏＳ ｓｅｇｍｅｎｔ ａｎｄ ｔｈｅ ｃｕｒｒｅｎｔ
ｏｒｄｅｒｌｙ． Ａｓｓｕｍｉｎｇ ｔｈｅ ｉｎｔｅｒｓｅｃｔｉｏｎ ｐｏｉｎｔ ＲＳＧ （ ｐｘ，ｐｙ） ｉｓ Ｍ（ ｘ，ｙ） ａｎｄ ｔｈｅ ｌｉｎｅ ｓｌｏｐｅ ｉｓ ｋｎ
ｃｏｏｒｄｉｎａｔｅｓ ｏｆ ｔｈｅ ｎｔｈ ｓｅｇｍｅｎｔ ａｎｄ ｔｈｅ ｎｔｈ ＧＣＳ ｔｉｌｅ ＝ （ ｙ２ｎ － ｙ１ｎ ） ／ （ ｘ ｎ － ｘ１ｎ ） ． Ｔｈｅｎ ｔｈｅ ｎｅｘｔ ｉｎｔｅｒｓｅｃｔｉｏｎ
２
ａｒｅ Ｐ ｎ （ ｘ１ｎ ， ｙ１ｎ ， ｚ１ｎ ） ａｎｄ Ｐ２ｎ （ ｘ２ｎ ， ｙ ｎ ， ｚ２ｎ ） ｐｏｉｎｔ Ｎ （ ｘ∗， ｙ∗ ） ｏｆ ｔｈｅ ｓｅｇｍｅｎｔ ａｎｄ ｔｈｅ ｎｅｘｔ
１ ２
ｒｅｓｐｅｃｔｉｖｅｌｙ， ｔｈｅ ｅｑｕａｔｉｏｎ ｏｆ ｔｈｉｓ ｓｅｇｍｅｎｔ ｋ ｎ ｃａｎ ＲＳＧ （ ｐ∗ｘ ， ｐ ∗ ） ｃａｎ ｂｅ ｓｈｏｗｅｄ ｇｒａｐｈｉｃａｌｌｙ ｉｎ
ｇ ｙ
ｂｅ ｅｘｐｒｅｓｓｅｄ ａｓ Ｆｉｇｕｒｅ １．
Ｆｉｇ．１．Ａｌｌ ｓｉｔｕａｔｉｏｎｓ ｏｆ ｔｈｅ ｎｅｘｔ ｉｎｔｅｒｓｅｃｔｉｏｎ ｐｏｉｎｔ ａｎｄ ｔｈｅ ｎｅｘｔ ＲＳＧ．
Ｔｈｅ ｆｏｕｒ ｖｅｒｔｅｘｅｓ ｏｆ ｏｎｅ ＲＳＧ ａｒｅ （ ｐｘ， ｐｙ， ｃｏｎｄｉｔｉｏｎ ｔｈｒｏｕｇｈ ｂｉｌｉｎｅａｒ ｉｎｔｅｒｐｏｌａｔｉｏｎ． Ｔｈｅ
ｐｚ１ ） ，（ ｐｘ ＋ａ，ｐｙ ，ｐｚ２ ） ，（ ｐｘ ＋ａ，ｐｙ ＋ａ，ｐｚ３ ） ａｎｄ （ ｐｘ ， ｅｑｕａｔｉｏｎ ｏｆ ｔｈｅ ｔｅｒｒａｉｎ ｓｕｒｆａｃｅ ｉｓ ｅｘｐｒｅｓｓｅｄ ａｓ
ｐｙ ＋ ａ， ｐｚ４ ） ｒｅｓｐｅｃｔｉｖｅｌｙ． Ｂａｓｅｄ ｏｎ ｔｈｅｓｅ
ｃｏｏｒｄｉｎａｔｅｓ， ｔｈｅ ｔｅｒｒａｉｎ ｓｕｒｆａｃｅ ｉｓ ｃｏｎｓｔｒｕｃｔｅｄ ｔｏ Ａｘｙ ＋ Ｂｘ ＋ Ｃｙ ＋ Ｄ － ｚ ＝ ０，
Ｗｈｅｒｅ
ｒｅｓｔｏｒｅ ｔｈｅ ｒｅｌａｔｉｖｅｌｙ ａｃｔｕａｌ ｔｅｒｒａｉｎ ｅｌｅｖａｔｉｏｎ
Ａ ＝ ｐｚ１ － ｐｚ２ ＋ ｐｚ３ － ｐｚ４ ，
ａ２
Ｂ ＝ （ｐｙ ＋ ａ） （ ｐｚ２ － ｐｚ１ ） ＋ ｐｘ（ ｐｚ４ － ｐｚ３ ） ，
ａ２
１８
