Đề tài Nghiên cứu một số bài toán tính chuyển toạ độ ứng dụng trong trắc địa công trình

Ón lùa chän cÇn ph¶i ®¶m b¶o c¶ hai yªu cÇu ®· nªu trªn. MÆt kh¸c, ®Ó tÝnh chuyÓn to¹ ®é vµ ®é cao c¸c ®iÓm khèng chÕ gi÷a hai hÖ quy chiÕu, cÇn ph¶i cã mét sè ®iÓm song trïng, tøc lµ c¸c ®iÓm ®· biÕt to¹ ®é vµ ®é cao trong c¶ hai hÖ to¹ ®é. KÕt qu¶ ®o GPS sÏ cho to¹ ®é c¸c ®iÓm khèng chÕ trong hÖ to¹ ®é vu«ng gãc kh«ng gian ®Þa t©m WGS – 84 (X, Y, Z). Tõ gi¸ trÞ (X, Y, Z) trong hÖ WGS – 84 sÏ tÝnh chuyÓn vÒ hÖ to¹ ®é tr¾c ®Þa trong cïng hÖ quy chiÕu (B, L, H), sÏ cã to¹ ®é tr¾c ®Þa (B, L) vµ ®é cao (H) cña ®iÓm ®ã. Täa ®é c¸c ®iÓm khèng chÕ thi c«ng trong tr¾c ®Þa c«ng tr×nh th­êng ®­îc tÝnh trong hÖ to¹ ®é gi¶ ®Þnh (x, y), ®é cao tÝnh trong hÖ ®é cao th­êng H . Nh­ vËy ®èi víi c¸c ®iÓm khèng chÕ thi c«ng, to¹ ®é mÆt b»ng vµ ®é cao cña c¸c ®iÓm khèng chÕ th­êng kh«ng cïng n»m trong mét hÖ quy chiÕu. §Ó tÝnh chuyÓn tõ ®é cao th­êng vÒ ®é cao tr¾c ®Þa, th­êng sö dông c«ng thøc tÝnh: H =  H (3.17) Trong ®ã: H lµ ®é cao tr¾c ®Þa cña ®iÓm xÐt H lµ ®é cao th­êng lÊy tõ kÕt qu¶ ®o thuû chuÈn  lµ dÞ th­êng ®é cao t¹i ®iÓm quan s¸t. Tõ c«ng thøc (3.17) cho thÊy, x¸c ®Þnh gi¸ trÞ dÞ th­êng ®é cao lµ mét vÊn ®Ò hÕt søc khã kh¨n. Nh­ vËy trong c¸c gi¸ trÞ to¹ ®é vµ ®é cao cña c¸c ®iÓm song trïng cã mét yÕu tè to¹ ®é vµ ®é cao H kh«ng thÓ x¸c ®Þnh chÝnh x¸c ë trong mét hÖ quy chiÕu (hÖ to¹ ®é vµ ®é cao thi c«ng), do ®ã sÏ lµm gi¶m ®é chÝnh x¸c tÝnh chuyÓn c¸c yÕu tè to¹ ®é mÆt b»ng (x, y), kh«ng ®¸p Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4848 øng ®­îc c¸c yªu cÇu cÇn thiÕt vÒ ®é chÝnh x¸c cÇn thiÕt cña l­íi khèng chÕ thi c«ng c«ng tr×nh. HÖ to¹ ®é th­êng dïng trong tr¾c ®Þa c«ng tr×nh lµ hÖ to¹ ®é vu«ng gãc ph¼ng Gauss – Kruger ®­îc thµnh lËp trªn mÆt ph¼ng dùa vµo phÐp chiÕu h×nh trô ngang ®ång gãc Gauss hoÆc UTM víi mói chiÕu 60 (hoÆc 30). Nh­ vËy nÕu sö dông ph­¬ng ph¸p tÝnh chuyÓn tõ to¹ ®é tr¾c ®Þa (B, L) sang hÖ to¹ ®é vu«ng gãc ph¼ng ph¶i thùc hiÖn phÐp chiÕu Gauss g©y ra sù biÕn d¹ng chiÒu dµi c¸c c¹nh l­íi khèng chÕ thi c«ng. §Ó tÝnh chuyÓn c¸c trÞ ®o GPS vÒ hÖ to¹ ®é thi c«ng th× cÇn ph¶i lùa chän ph­¬ng ph¸p tÝnh chuyÓn sao cho t¸ch biÖt ®­îc gi¸ trÞ to¹ ®é (X, Y) vµ gi¸ trÞ ®é cao H, khi ®ã ®é chÝnh x¸c cña phÐp tÝnh chuyÓn to¹ ®é sÏ ®¸p øng ®­îc yªu cÇu cÇn thiÕt cña l­íi khèng chÕ thi c«ng c«ng tr×nh ®­îc thµnh lËp b»ng c«ng nghÖ GPS. Bµi to¸n tÝnh chuyÓn to¹ ®é tõ hÖ to¹ ®é vu«ng gãc kh«ng gian ®Þa t©m vÒ hÖ to¹ ®é ®Þa diÖn ch©n trêi t¹i ®iÓm quan s¸t kÕt hîp víi phÐp tÝnh chuyÓn gi÷a hai hÖ to¹ ®é vu«ng gãc ph¼ng (phÐp tÝnh chuyÓn Helmert) sÏ ®¸p øng ®­îc yªu cÇu ®é chÝnh x¸c cÇn thiÕt cña l­íi khèng chÕ thi c«ng c«ng tr×nh [6]. 3.2.3 ThuËt to¸n tÝnh chuyÓn B­íc 1: TÝnh chuyÓn to¹ ®é c¸c ®iÓm ®o GPS tõ hÖ to¹ ®é vu«ng gãc kh«ng gian ®Þa t©m WGS – 84 (X, Y, Z) sang hÖ to¹ ®é tr¾c ®Þa cïng hÖ quy chiÕu (B, L, H) theo c¸c c«ng thøc tõ (2.32) ®Õn (2.35) B­íc 2: Chän ®iÓm gèc cña hÖ to¹ ®é ®Þa diÖn lµ ®iÓm träng t©m cña l­íi khèng chÕ thi c«ng c«ng tr×nh (®iÓm cã to¹ ®é vµ ®é cao lµ gi¸ trÞ trung b×nh cña c¸c ®iÓm khèng chÕ thi c«ng). §iÓm gèc cña hÖ to¹ ®é ®Þa diÖn còng cã thÓ chän lµ mét trong c¸c ®iÓm song trïng cña hai hÖ to¹ ®é vµ ®é cao lµ ®é cao cña mÆt chiÕu ®· sö dông ®Ó thiÕt kÕ c«ng tr×nh. n B B n i o   1 (3.18) n H H n i o   1 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4849 B­íc 3: Sö dông c¸c c«ng thøc tõ (2.51) ®Õn (2.56) ®Ó tÝnh chuyÓn c¸c ®iÓm khèng chÕ trong l­íi khèng chÕ vÒ hÖ to¹ ®é ®Þa diÖn ®· chän. Nh­ vËy víi mçi mét ®iÓm tÝnh chuyÓn sÏ cã thµnh phÇn to¹ ®é z trong hÖ to¹ ®é ®Þa diÖn. Trong tr¾c ®Þa c«ng tr×nh, ph¹m vi thi c«ng cña c¸c c«ng tr×nh th­êng nhá vµ ®Þa h×nh t­¬ng ®èi b»ng ph¼ng, nªn cã thÓ bá qua sè hiÖu chØnh do ®é lÖch d©y däi. Khi ®ã gi¸ trÞ z cña c¸c ®iÓm tÝnh chuyÓn sÏ lµ ®é cao cña c¸c ®iÓm ®ã trong hÖ to¹ ®é ®Þa diÖn, tøc lµ ®é cao c¸c ®iÓm so víi ®é cao trung b×nh cña khu vùc( hoÆc so víi mÆt ph¼ng chiÕu ®· lùa chän khi thiÕt kÕ c«ng tr×nh). Trong hÖ to¹ ®é ®Þa diÖn chªnh cao cña c¸c ®iÓm khèng chÕ so víi mÆt chiÕu ®· ®­îc gi¶m ®i rÊt nhiÒu. MÆt kh¸c, khi ®é cao c¸c ®iÓm khèng chÕ so víi mÆt chiÕu kh«ng v­ît qu¸ 32m, trong tr­êng hîp nµy sai sè do phÐp chiÕu cña c¸c c¹nh ®o vÒ mÆt chiÕu sÏ kh«ng v­ît qu¸ 1/200000 vµ cã thÓ bá qua. Khi ®ã chóng ta cã thÓ c¸c ®iÓm cña l­íi khèng chÕ thi c«ng cïng n»m trªn mÆt ph¼ng chiÕu. B­íc 4: Dùa vµo c¸c ®iÓm song trïng trong hÖ to¹ ®é thi c«ng, tÝnh c¸c tham sè tÝnh chuyÓn theo phÐp tÝnh chuyÓn gi÷a hai hÖ to¹ ®é vu«ng gãc ph¼ng ( phÐp tÝnh chuyÓn Helmert) theo c«ng thøc tõ (2.15) ®Õn (2.22). B­íc 5: TÝnh to¹ ®é cho c¸c ®iÓm ®o GPS cßn l¹i trong hÖ to¹ ®é thi c«ng theo c¸c tham sè tÝnh chuyÓn ë môc (4). ThuËt to¸n vµ tr×nh tù tÝnh chuyÓn theo quy tr×nh trªn ®· ®­îc lËp tr×nh thµnh ch­¬ng tr×nh tÝnh to¸n. Ch­¬ng tr×nh ®­îc viÕt b»ng ng«n ng÷ FORTRAN 77 ( phô lôc) 3.2.4 TÝnh to¸n thùc nghiÖm §Ó ®¸nh gi¸ møc ®é tin cËy cña ph­¬ng ph¸p tÝnh chuyÓn c¸c ®iÓm ®o GPS vÒ hÖ to¹ ®é thi c«ng t«i tiÕn hµnh tÝnh chuyÓn ®èi víi l­íi khèng chÕ thi c«ng khu c«ng nghiÖp Yªn Phong (B¾c Ninh). S¬ ®å l­íi nh­ h×nh vÏ (h×nh 3.2). Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4850 B¶n ®å ®Þa h×nh dïng ®Ó thiÕt kÕ khu c«ng nghiÖp Yªn Phong ®­îc thµnh lËp dùa vµo l­íi khèng chÕ thµnh lËp trong giai ®o¹n kh¶o s¸t ( hÖ to¹ ®é VN2000, mói chiÕu 30, kinh tuyÕn trôc 105030’) víi 2 ®iÓm ®Þnh vÞ lµ IV- 15 vµ IV- 16. Trong giai ®o¹n thi c«ng, ng­êi ta ®· sö dông c«ng nghÖ GPS ®Ó thµnh lËp l­íi khèng chÕ thi c«ng bao gåm 18 ®iÓm trong hÖ to¹ ®é VN2000, mói 30, kinh tuyÕn trôc 105008’. §Ó phï hîp víi hÖ to¹ ®é thiÕt kÕ c«ng tr×nh, ®¬n vÞ thi c«ng ®· tiÕn hµnh tÝnh chuyÓn to¹ ®é c¸c ®iÓm ®o GPS tõ khu vùc cã kinh tuyÕn trôc 105008’ hÖ to¹ ®é VN2000 Ellipxoid WGS 84, l­íi chiÕu UTM vÒ khu vùc cã kinh tuyÕn trôc 105030’, mói chiÕu 30 b»ng phÇn mÒm H×nh 3.3 L­íi GPS khu c«ng nghiÖp Yªn Phong – B¾c Ninh YP2 YP1 104511 IV-14 YP7 IV-YP4 YP3 YP4 YP5 YP6 IV-16 IV-15 IV-YP3 YP8 YP9 YP10 YP11 YP12 IV-YP2 IV-YP1 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4851 Geotool cña Bé Tµi Nguyªn M«i tr­êng vµ ®o l¹i tÊt c¶ c¸c c¹nh cña l­íi b»ng m¸y Toµn ®¹c ®iÖn tö TCR - 303. B¶ng 3: To¹ ®é c¸c ®iÓm ®o GPS trªn kinh tuyÕn 105030’ (tr­íc khi tÝnh chuyÓn) Sè TT Tªn ®iÓm X(m) Y(m) 1 104511 2345742.250 550095.739 2 105486 2343045.007 552883.025 3 IV14 2345912.888 550834.734 4 IV15 2344375.319 551780.890 5 IV16 2344662.770 551279.180 6 IVYP1 2343442.276 552726.462 7 IVIP2 2343792.718 552356.510 8 IVYP3 2345352.691 551102.448 9 YP-01 2345558.569 549920.372 10 YP-02 2345308.288 549634.047 11 YP-03 2344947.767 549737.296 12 YP-04 2344905.533 550106.462 13 YP-05 2344886.251 550272.613 14 YP-06 2344813.265 550797.272 15 YP-07 2345459.307 550399.422 16 YP-08 2345267.399 551381.897 17 YP-09 2345171.892 551701.792 18 YP-10 2344821.579 551992.610 19 YP-11 2344705.435 552338.966 20 YP-12 2344279.283 552342.122 Tuy nhiªn kh«ng thÓ sö dông to¹ ®é c¸c ®iÓm ®o GPS ®Ó ph¸t triÓn l­íi cÊp thÊp vµ bè trÝ c«ng tr×nh ®­îc do sù biÕn d¹ng vÒ chiÒu dµi c¹nh cña l­íi Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4852 GPS so víi chiÒu dµi c¹nh ®o trùc tiÕp trªn mÆt ®Êt b»ng m¸y Toµn ®¹c ®iÖn tö. So s¸nh kÕt qu¶ gi÷a 2 c¸ch (b¶ng 4). B¶ng 4: So s¸nh chiÒu dµi c¹nh theo ph­¬ng ¸n ch­a tÝnh chuyÓn C¹nh ®o S ®o b»ng m¸y TC- 303 (m) S tÝnh trong kinh tuyÕn 105030’(m) Sai kh¸c (mm)§iÓm ®Çu §iÓm cuèi YP5 YP6 529.745 529.711 34.0 YP5 YP7 586.956 586.918 38.0 YP7 IV-YP4 711.104 711.064 40.0 YP7 IV-14 628.722 628.675 47.0 YP7 105511 415.088 415.066 22.0 YP1 YP2 380.337 380.292 45.0 YP2 YP3 375.042 375.014 28.0 YP3 105511 871.653 871.598 55.0 YP6 YP8 740.345 740.286 59.0 YP8 IV-14 846.265 846.193 72.0 YP9 YP10 455.311 455.295 16.0 YP11 YP12 426.188 426.163 25.0 YP6 IV-16 504.881 504.860 21.0 YP9 IV-16 661.701 661.669 32.0 Tõ kÕt qu¶ so s¸nh ë b¶ng 1 cho thÊy sù sai kh¸c vÒ chiÒu dµi gi÷a c¸c c¹nh ®­îc tÝnh tõ to¹ ®é ®o GPS vµ chiÒu dµi ®o trùc tiÕp trªn mÆt ®Êt lµ kh¸ lín. Trong l­íi ®Æc biÖt cã c¹nh YP8 - IV14 sai kh¸c gi÷a c¹nh ®o b»ng m¸y TCR- 303 vµ c¹nh tÝnh trong mói 105030’ sai kh¸c tíi 72 mm. Sù sai kh¸c nµy sÏ lµm cho l­íi khèng chÕ thi c«ng kh«ng ®¸p øng ®­îc ®é chÝnh x¸c cÇn thiÕt khi bè trÝ c«ng tr×nh. V× vËy t«i ®· sö dông quy tr×nh tÝnh chuyÓn ®· nªu trªn vµ theo ch­¬ng tr×nh m¸y tÝnh ®­îc thµnh lËp t«i ®· tiÕn hµnh tÝnh chuyÓn c¸c ®iÓm ®o GPS Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4853 cña l­íi khèng chÕ thi c«ng vÒ hÖ to¹ ®é thi c«ng c«ng tr×nh. Sau khi tÝnh chuyÓn to¹ ®é c¸c ®iÓm ®o GPS vÒ hÖ to¹ ®é thi c«ng ®­îc ®Þnh vÞ b»ng 2 ®iÓm IV- 15 vµ IV-16 ta cã to¹ ®é c¸c ®iÓm l­íi khèng chÕ thi c«ng trong hÖ to¹ ®é thi c«ng c«ng tr×nh vµ chiÒu dµi c¸c c¹nh cña l­íi so víi chiÒu dµi ®o trùc tiÕp trªn mÆt ®Êt b»ng m¸y toµn ®¹c ®iÖn tö TCR-303. B¶ng 5: To¹ ®é GPS sau khi tÝnh chuyÓn Sè TT Tªn ®iÓm X(m) Y(m) 1 104511 2345742.169 550095.529 2 105486 2343044.856 552882.995 3 IV14 2345912.936 550834.589 4 IV15 2344375.319 551780.890 5 IV16 2344662.784 551279.156 6 IVYP1 2343442.17 552726.457 7 IVIP2 2343792.672 552356.516 8 IVYP3 2345352.721 551102.368 9 YP-01 2345558.452 549920.172 10 YP-02 2345308.089 549633.868 11 YP-03 2344947.570 549737.162 12 YP-04 2344905.421 550106.365 13 YP-05 2344886.16 550272.519 14 YP-06 2344813.241 550797.214 15 YP-07 2345459.266 550399.274 16 YP-08 2345267.448 551381.852 17 YP-09 2345171.933 551701.769 18 YP-10 2344821.602 551992.614 19 YP-11 2344705.427 552338.983 20 YP-12 2344279.260 552342.137 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4854 B¶ng 6: So s¸nh chiÒu dµi c¹nh sau khi tÝnh chuyÓn C¹nh ®o S ®o b»ng m¸y TC- 303 (m) S tÝnh trong kinh tuyÕn 105030’(m) Sai kh¸c (mm)§iÓm ®Çu §iÓm cuèi YP5 YP6 529.745 529.737 8.0 YP5 YP7 586.956 586.955 1.0 YP7 IV-YP4 711.104 711.122 -18.0 YP7 IV-14 628.722 628.741 -19.0 YP7 105511 415.088 415.084 4.0 YP1 YP2 380.337 380.331 6.0 YP2 YP3 375.042 375.024 18.0 YP3 105511 871.653 871.673 -20.0 YP6 YP8 740.345 740.341 4.0 YP8 IV-14 846.265 846.258 7.0 YP9 YP10 455.311 455.326 -15.0 YP11 YP12 426.188 426.178 10.0 YP6 IV-16 504.881 504.890 -9.0 YP9 IV-16 661.701 661.690 11.0 KÕt qu¶ sai sè chiÒu dµi c¹nh cña l­íi GPS sau khi tÝnh chuyÓn cho thÊy: sù biÕn d¹ng vÒ chiÒu dµi ®· gi¶m ®i ®¸ng kÓ (sai kh¸c lín nhÊt lµ 20mm vµ nhá nhÊt lµ 1mm). Sù sai kh¸c nµy sÏ ®¶m b¶o ®­îc tÝnh chÊt quan träng lµ biÕn d¹ng chiÒu dµi c¸c c¹nh l­íi khèng chÕ thi c«ng so víi chiÒu dµi thùc trªn mÆt ®Êt n»m trong h¹n sai cho phÐp. Nh­ vËy b»ng gi¶i ph¸p tÝnh chuyÓn to¹ ®é ®· ®¸p øng ®­îc 2 yªu cÇu cÇn thiÕt cña l­íi khèng chÕ thi c«ng: ®ã lµ sù ®ång nhÊt gi÷a hÖ to¹ ®é thiÕt kÕ vµ hÖ to¹ ®é thi c«ng c«ng tr×nh mµ vÉn ®¶m b¶o ®­îc sù biÕn d¹ng cña chiÒu dµi c¹nh n»m trong giíi h¹n cho phÐp nh»m ®¶m b¶o ®é chÝnh x¸c bè trÝ c«ng tr×nh. §©y còng chÝnh lµ hiÖu qu¶ cña gi¶i ph¸p tÝnh chuyÓn to¹ ®é c¸c ®iÓm ®o GPS vÒ hÖ to¹ ®é thi c«ng c«ng tr×nh. Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4855 3.4 Bµi to¸n tÝnh chuyÓn vÒ ®é cao khu vùc 3.4.1.Lựa chọn phương ¸n tÝnh chuyÓn B¶n chÊt cña ph­¬ng ph¸p lµ tÝnh chuyÓn to¹ ®é c¸c ®iÓm cña l­íi khèng chÕ thi c«ng b»ng c¸ch tÝnh chuyÓn chiÒu dµi c¸c c¹nh cña l­íi theo ®é cao mÆt chiÕu. Mặt chiếu được chọn trong x©y dựng c«ng tr×nh là mặt cã độ cao trung b×nh của khu vực x©y dựng c«ng tr×nh. Khi hiệu số độ cao mặt đất và mặt chiếu < 32m th× cã thể bỏ qua số hiệu chỉnh SH , nếu lớn hơn th× phải tÝnh số hiệu chỉnh do độ cao. Trªn thùc tÕ, khi x©y dùng c¸c c«ng tr×nh cã ph¹m vi kh«ng lín vµ ë khu vùc ®Þa h×nh t­¬ng ®èi b»ng ph¼ng (cÇu v­ît s«ng, tuyÕn ®­êng giao th«ng ë ®ång b»ng…), hiệu số độ cao mặt đất và mặt chiếu < 32m th× cã thể bỏ qua số hiệu chỉnh SH. §èi víi nh÷ng c«ng tr×nh x©y dùng ë vïng nói nh­: thuû ®iÖn, c¸c tuyÕn ®­êng giao th«ng, ®­êng hÇm…®é cao cña khu vùc x©y dùng so víi mÆt n­íc biÓn trung b×nh rÊt lín, sù biÕn d¹ng cña c¸c trÞ ®o trªn mÆt ®Êt lín lµm gi¶m ®é chÝnh x¸c cña c¸c c«ng tr×nh, v× vËy ®Ó ®¶m b¶o ®é chÝnh x¸c bè trÝ c¸c c«ng tr×nh th× ta ph¶i tÝnh to¹ ®é c¸c ®iÓm cña l­íi khèng chÕ thi c«ng vÒ ®é cao trung b×nh cña khu vùc x©y dùng. 3.4.2 ThuËt to¸n tÝnh chuyÓn - B×nh sai l­íi khèng chÕ thi c«ng ë mÆt chiÕu H0 = 0.0 m - TÝnh to¹ ®é träng t©m cña l­íi khèng chÕ thi c«ng ë mÆt chiÕu H0 = 0.0 m theo c«ng thøc   n x x i0   n y y i0 (3.19) n: tæng sè ®iÓm cÇn tÝnh chuyÓn. Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4856 - TÝnh gia sè to¹ ®é vµ chiÒu dµi gi÷a tÊt c¶ c¸c c¹nh cÇn tÝnh chuyÓn so víi ®iÓm träng t©m cña l­íi khèng chÕ thi c«ng. 2 0 2 0 )()( yyxxS iii  (3.20) - TÝnh chuyÓn chiÒu dµi c¹nh cña l­íi khèng chÕ thi c«ng vÒ ®é cao trung b×nh HTB cña khu vùc x©y dùng theo c«ng thøc: KSS ii ' (3.21) Víi R HR K TB  : hÖ sè biÕn d¹ng chiÒu dµi R: B¸n kÝnh tr¸i ®Êt (= 6370 km) - TÝnh to¹ ®é cña c¸c ®iÓm l­íi khèng chÕ thi c«ng trong hÖ to¹ ®é míi. iii SxX cos ' 0  iii SyY sin ' 0  (3.22) 0 0 xx yy arctg i i i   (3.23) ThuËt to¸n vµ tr×nh tù tÝnh chuyÓn theo quy tr×nh trªn ®· ®­îc lËp tr×nh thµnh quy tr×nh tÝnh to¸n. Ch­¬ng tr×nh ®­îc viÕt b»ng ng«n ng÷ FORTRAN 77 (phô lôc). 3.4.3 TÝnh to¸n thùc nghiÖm C«ng tr×nh thuû lîi - thuû ®iÖn Qu¶ng TrÞ ®­îc x©y dùng trªn s«ng Rµo Qu¸n thuéc ®Þa phËn x· T©n Hîp huyÖn H­íng Ho¸ tØnh Qu¶ng TrÞ. Do ®Æc ®iÓm c«ng tr×nh c¸c h¹ng môc bè trÝ trªn ph¹m vi réng tr¶i dµi trªn 6 km chªnh cao ®Þa h×nh lín, tuyÕn ®Ëp vµ nhµ m¸y c¸ch xa nhau, ®­êng hÇm dµi vµ nhá thi c«ng tõ 3 cöa hÇm phô, ®iÒu kiÖn thi c«ng khã kh¨n tiÕn ®é gÊp rót do vËy l­íi tam gi¸c thuû c«ng vµ thuû chuÈn thuû c«ng lµ hÖ thèng l­íi khèng chÕ mÆt b»ng vµ ®é cao chÝnh x¸c nèi liÒn tÊt c¶ c¸c h¹ng môc cña c«ng tr×nh phôc vô thi c«ng. §Ó ®¸p øng ®­îc c¸c yªu cÇu kÜ thuËt cña c«ng tr×nh, l­íi tam gi¸c thuû c«ng ®­îc x©y dùng gåm 19 ®iÓm cã kÝ hiÖu tõ TC01…TC19 t¹o thµnh ®å h×nh chuçi tam gi¸c. L­íi ®­îc ®o nèi víi c¸c ®iÓm (RQ-2, RQ-3, RQ-6, RQ- Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4857 7, RQ-12, NM-2, NM-3, IV23) ®· ®­îc thµnh lËp trong giai ®o¹n tr­íc thuéc hÖ to¹ ®é thi c«ng c«ng tr×nh. S¬ ®å l­íi nh­ h×nh vÏ (3.4). NM03 H×nh 3.4 L­íi khèng chÕ thuû c«ng thuû ®iÖn Qu¶ng TrÞ TC14 TC18 TC19 TC16 TC17 TC15 TC13 TC11 TC09 TC07 TC05 RQ02 TC02 RQ03 TC01 TC03 TC04 TC06 TC08 TC10 TC12 RQ07 IV23 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4858 L­íi ®­îc ®o b»ng m¸y Toµn ®¹c ®iÖn tö TC – 1700 (®é chÝnh x¸c ®o gãc Mβ = 2’’, ®é chÝnh x¸c ®o c¹nh MS = 3mm  1 ppm) kÕt hîp víi m¸y ®Þnh vÞ GPS 4600LS mét tÇn sè (®é chÝnh x¸c ®o c¹nh MS =  5mm +1.10-6D). §Ó kÕt nèi m¹ng l­íi thuû c«ng vÒ hÖ to¹ ®é thiÕt kÕ c«ng tr×nh ng­êi ta tiÕn hµnh ®o GPS tÊt c¶ c¸c ®iÓm cña l­íi thuû c«ng bao gåm c¶ c¸c ®iÓm ®Þnh vÞ c«ng tr×nh theo hÖ quy chiÕu ®­îc chän lµ Ellipsoid Krasowki, phÐp chiÕu Gauss, kinh tuyÕn trung ­¬ng 106015’, hÖ to¹ ®é HN - 72. KÕt qu¶ ®o GPS nh­ b¶ng 7. B¶ng 7: B¶ng thµnh qu¶ to¹ ®é ph¼ng vµ ®é cao b×nh sai l­íi GPS Sè TT KÝ hiÖu ®iÓm To¹ ®é, ®é cao X(m) Y(m) h(m) 1 IV23 1841048.047 552276.066 439.619 2 NM03 1841951.005 552828.500 133.609 3 RQ02 1847011.622 548278.756 498.579 4 RQ03 1846356.134 547862.587 470.848 5 RQ07 1843489.247 548652.200 525.217 6 TC01 1846589.592 548334.242 516.959 7 TC02 1846339.288 547838.343 472.629 8 TC03 1844547.526 547706.836 481.909 9 TC04 1845490.866 548315.949 564.498 10 TC05 1846056.375 549029.681 490.704 11 TC06 1844977.195 548547.866 520.154 12 TC07 1845122.788 549320.350 506.534 13 TC08 1844063.599 549058.368 489.770 14 TC09 1844532.300 549949.881 487.674 15 TC10 1843631.528 549908.075 504.494 16 TC11 1844231.644 550726.325 449.672 17 TC12 1842747.715 550387.493 505.700 18 TC13 1843253.460 551117.908 495.661 19 TC14 1842123.255 550981.437 486.234 20 TC15 1842564.458 551982.967 441.984 21 TC16 1841871.529 552210.565 425.618 22 TC17 1842343.014 552662.639 200.351 23 TC18 1841685.178 553078.315 160.613 24 TC19 1841533.436 552349.219 422.170 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4859 Sau khi cã sè liÖu ®o ®¹c vµ tÝnh to¸n cña l­íi tam gi¸c thuû c«ng theo hai ph­¬ng ph¸p: Toµn ®¹c ®iÖn tö vµ c«ng nghÖ GPS, tiÕn hµnh b×nh sai tæng hîp cho m¹ng l­íi trªn bÒ mÆt H0 = 0.0 m. KÕt qu¶ b×nh sai nh­ b¶ng 8. B¶ng 8: B¶ng thµnh qu¶ to¹ ®é b×nh sai l­íi kÕt hîp GPS – tam gi¸c ( trªn mÆt ph¼ng H0 = 0.0 m) Sè TT KÝ hiÖu ®iÓm Toa ®é X (m) Y(m) 1 TC01 1846589.595 548334.243 2 TC02 1846339.291 547838.341 3 TC03 1845487.523 547706.836 4 TC04 1845490.864 548315.948 5 TC05 1846056.376 549029.683 6 TC06 1844977.194 548547.865 7 TC07 1845122.788 549320.350 8 TC08 1844063.599 549058.368 9 TC09 1844532.300 549949.882 10 TC10 1843631.528 549908.075 11 TC11 1844231.643 550726.325 12 TC12 1842747.715 550387.494 13 TC13 1843253.461 551117.908 14 TC14 1842123.255 550981.438 15 TC15 1842564.461 551982.969 16 TC16 1841871.529 552210.565 17 TC17 1842343.012 552662.635 18 TC18 1841685.178 553078.314 19 TC19 1841533.434 552349.221 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4860 Tuy nhiªn kh«ng thÓ sö dông to¹ ®é c¸c ®iÓm sau b×nh sai ®Ó bè trÝ c«ng tr×nh ®­îc do sù biÕn d¹ng vÒ chiÒu dµi c¸c c¹nh cña l­íi so víi chiÒu dµi c¹nh ®o trùc tiÕp trªn mÆt ®Êt b»ng m¸y Toµn ®¹c ®iÖn tö. So s¸nh chiÒu dµi c¹nh tr­íc khi tÝnh chuyÓn víi chiÒu dµi c¹nh ®­îc ®o trªn mÆt ®Êt b»ng m¸y TC – 700. B¶ng 9: So s¸nh chiÒu dµi c¹nh theo ph­¬ng ¸n ch­a tÝnh chuyÓn (trªn mÆt ph¼ng H0 = 0.0 m) C¹nh ®o S ®o b»ng m¸y TC-700 S tr­íc tÝnh chuyÓn Sai kh¸c (mm) TC01-TC02 555.525 555.488 37.0 TC01-TC04 1098.943 1098.878 65.0 TC01-TC05 876.381 876.331 50.0 TC02-TC03 861.905 861.854 51.0 TC02-TC04 973.665 973.615 50.0 TC04-TC05 910.673 910.612 61.0 TC07-TC04 1069.776 1069.721 55.0 TC05-TC06 1181.915 1181.852 63.0 TC07-TC09 863.173 863.125 48.0 TC06-TC07 786.132 786.085 47.0 TC09-TC11 832.663 832.621 42.0 TC11-TC13 1053.689 1053.649 40.0 TC13-TC15 1105.962 1105.916 46.0 TC15-TC18 1404.650 1404.606 44.0 TC15-TC19 1094.192 1094.147 45.0 TC16-TC17 653.221 653.195 26.0 TC16-TC19 365.435 365.423 12.0 TC18-TC19 744.746 744.716 30.0 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4861 Tõ kÕt qu¶ so s¸nh ë b¶ng 9 cho thÊy sù sai kh¸c vÒ chiÒu dµi gi÷a c¸c c¹nh ®­îc tÝnh tõ to¹ ®é sau b×nh sai hçn hîp GPS – tam gi¸c vµ chiÒu dµi ®o trùc tiÕp trªn mÆt ®Êt lµ kh¸ lín (lín nhÊt lµ 65.0 mm vµ nhá nhÊt lµ 12.0 mm). Sù sai kh¸c nµy lµm cho l­íi khèng chÕ thuû c«ng kh«ng ®¸p øng ®­îc ®é chÝnh x¸c khi bè trÝ c«ng tr×nh. §é chªnh cao gi÷a c¸c h¹ng môc cña nhµ m¸y thuû ®iÖn Qu¶ng TrÞ lµ rÊt lín nh»m môc ®Ých lîi dông c¸c cét n­íc cao ®Ó t¹o n¨ng l­îng. §é cao trung b×nh cña thuû ®iÖn Qu¶ng TrÞ lµ 300m, v× vËy to¹ ®é c¸c ®iÓm sau khi b×nh sai hçn hîp ®­îc tÝnh chuyÓn lªn bÒ mÆt cã ®é cao trung b×nh cña khu vùc x©y dùng theo thuËt to¸n ®· nªu trªn. Sau khi tiÕn hµnh tÝnh chuyÓn to¹ ®é c¸c ®iÓm cña l­íi khèng chÕ thuû c«ng vÒ hÖ to¹ ®é c«ng tr×nh ë ®é cao H0 = 300 m ta thu ®­îc kÕt qu¶ nh­ b¶ng 10. B¶ng 10: B¶ng thµnh qu¶ tÝnh chuyÓn to¹ ®é (trªn mÆt ph¼ng H0 =300 m) STT Tªn ®iÓm To¹ ®é tr­íc tÝnh chuyÓn To¹ ®é c«ng tr×nh X(m) Y(m) X(m) Y(m) 1 TC01 1846589.595 548334.243 1846589.720 548334.156 2 TC02 1846339.291 547838.341 1846339.405 547838.231 3 TC03 1845487.523 547706.836 1845487.596 547706.719 4 TC04 1845490.864 548315.948 1845490.938 548315.860 5 TC05 1846056.376 549029.683 1846056.476 549029.628 6 TC06 1844977.194 548547.865 1844977.243 548547.788 7 TC07 1845122.788 549320.350 1845122.844 549320.309 8 TC08 1844063.599 549058.368 1844063.605 549058.315 9 TC09 1844532.300 549949.882 1844532.328 549949.871 10 TC10 1843631.528 549908.075 1843631.514 550461.039 11 TC11 1844231.643 550726.325 1844231.657 550726.351 12 TC12 1842747.715 550387.494 1842747.659 550387.504 13 TC13 1843253.461 551117.908 1843253.429 551117.952 14 TC14 1842123.255 550981.438 1842123.17 550981.475 15 TC15 1842564.461 551982.969 1842564.397 551983.054 16 TC16 1841871.529 552210.565 1841871.432 552210.660 17 TC17 1842343.012 552662.635 1842342.937 552662.752 18 TC18 1841685.178 553078.314 1841685.072 553078.450 19 TC19 1841533.434 552349.221 1841533.321 552349.323 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4862 Tõ to¹ ®é c¸c ®iÓm sau tÝnh chuyÓn, so s¸nh chiÒu dµi mét sè c¹nh ®o trªn mÆt ®Êt b»ng m¸y Toµn ®¹c ®iÖn tö TC- 700 víi chiÒu dµi c¸c c¹nh sau khi tÝnh chuyÓn (b¶ng 11). B¶ng 11: B¶ng so s¸nh chiÒu dµi c¹nh sau tÝnh chuyÓn (trªn mÆt ph¼ng H0 = 300m) C¹nh ®o S ®o b»ng m¸y TC-700 (m) S sau tÝnh chuyÓn (m) Sai kh¸c (mm) TC01-TC02 555.525 555.517 7.9 TC01-TC04 1098.943 1098.935 8.3 TC01-TC05 876.381 876.374 7.2 TC02-TC03 861.905 861.900 4.8 TC02-TC04 973.665 973.666 1.1 TC04-TC05 910.673 910.659 13.8 TC07-TC04 1069.776 1069.772 4.5 TC05-TC06 1181.915 1181.911 4.1 TC07-TC09 863.173 863.165 7.4 TC06-TC07 786.132 786.122 9.3 TC09-TC11 832.663 832.660 2.4 TC11-TC13 1053.689 1053.699 - 10.0 TC13-TC15 1105.962 1105.968 - 6.7 TC15-TC18 1404.650 1404.672 - 22.6 TC15-TC19 1094.192 1094.198 - 6.3 TC16-TC17 653.221 653.225 - 4.6 TC16-TC19 365.435 365.439 - 4.8 TC18-TC19 744.746 744.751 - 5.7 KÕt qu¶ sai sè chiÒu dµi c¹nh cña l­íi tam gi¸c thuû c«ng sau khi tÝnh chuyÓn cho thÊy: sù biÕn d¹ng vÒ chiÒu dµi ®· gi¶m ®i ®¸ng kÓ (lín nhÊt lµ 22.6 mm vµ nhá nhÊt lµ 1.1 mm) . Sù sai kh¸c nµy sÏ ®¶m b¶o ®­îc tÝnh chÊt quan träng lµ biÕn d¹ng chiÒu dµi c¸c c¹nh l­íi tam gi¸c thuû c«ng so víi chiÒu dµi thùc trªn mÆt ®Êt n»m trong h¹n sai cho phÐp. Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4863 KÕt luËn Sau mét thêi gian nghiªn cøu, häc hái vµ ®Æc biÖt ®­îc sù h­íng dÉn tËn t×nh cña thÇy gi¸o TrÇn ViÕt TuÊn cïng c¸c thÇy c« trong bé m«n, ®Õn nay b¶n ®å ¸n tèt nghiÖp víi ®Ò tµi:“ Nghiªn cøu mét sè bµi to¸n tÝnh chuyÓn to¹ ®é øng dông trong tr¾c ®Þa c«ng tr×nh” ®· hoµn thµnh. Tõ b¶n ®å ¸n nµy, t«i rót ra ®­îc mét sè kÕt luËn sau: 1./ Tõ kÕt qu¶ nghiªn cøu vµ tÝnh to¸n thùc nghiÖm cho thÊy cÇn ph¶i tÝnh chuyÓn to¹ ®é c¸c ®iÓm cña l­íi khèng chÕ thi c«ng (theo c«ng nghÖ truyÒn thèng hoÆc c«ng nghÖ GPS) vÒ hÖ to¹ ®é thi c«ng c«ng tr×nh. C¸c bµi to¸n tÝnh chuyÓn to¹ ®é ®ãng vai trß quan träng trong tr¾c ®Þa c«ng tr×nh. 2./ Khi sö dông c«ng nghÖ GPS ®Ó thµnh lËp l­íi khèng chÕ thi c«ng trong tr¾c ®Þa c«ng tr×nh, th× cÇn ph¶i tÝnh chuyÓn to¹ ®é c¸c ®iÓm ®o GPS vÒ hÖ to¹ ®é thi c«ng c«ng tr×nh. KÕt qu¶ tÝnh chuyÓn sÏ ®¶m b¶o tÝnh ®ång nhÊt gi÷a hÖ to¹ ®é thiÕt kÕ vµ hÖ to¹ ®é thi c«ng c«ng tr×nh còng nh­ ®¶m b¶o ®é chÝnh x¸c cÇn thiÕt khi bè trÝ c«ng tr×nh. 3./ Tuú thuéc vµo tõng h¹ng môc c«ng tr×nh mµ lùa chän c¸c bµi to¸n tÝnh chuyÓn to¹ ®é phï hîp. Bµi to¸n tÝnh chuyÓn to¹ ®é phï hîp ph¶i ®¶m b¶o 2 ®iÒu kiÖn: - Sù ®ång nhÊt gi÷a hÖ to¹ ®é thiÕt kÕ vµ hÖ to¹ ®é thi c«ng c«ng tr×nh. - Sù biÕn d¹ng chiÒu dµi c¹nh lµ nhá nhÊt. Cuèi cïng t«i xin ch©n thµnh c¶m ¬n thÇy gi¸o TrÇn ViÕt TuÊn cïng toµn bé thÇy c« trong bé m«n vµ c¸c b¹n ®· gióp ®ì t«i hoµn thµnh ®å ¸n. Hµ Néi, th¸ng 6 n¨m 2008 Sinh viªn thùc hiÖn: Vò ThÞ Hµ Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4864 Tµi liÖu tham kh¶o [1]. Phan V¨n HiÕn (chñ biªn)(2004), Tr¾c ®Þa c«ng tr×nh. Nhµ xuÊt b¶n giao th«ng vËn t¶i.Hµ Néi [2]. §ç Ngäc §­êng - §Æng Nam Chinh (2007), Bµi gi¶ng c«ng nghÖ GPS. [3]. Ph¹m Hoµng L©n (chñ biªn) (1999), Gi¸o tr×nh tr¾c ®Þa cao cÊp (häc phÇn 3: tr¾c ®Þa mÆt cÇu). [4]. Hoµng Ngäc Hµ (2006), B×nh sai tÝnh to¸n l­íi Tr¾c ®Þa vµ GPS. Nhµ xuÊt b¶n khoa häc kü thuËt. Hµ Néi [5]. NguyÔn Quang Th¾ng – TrÇn ViÕt TuÊn (2007), Tr¾c ®Þa c«ng tr×nh c«ng nghiÖp – thµnh phè. [6]. TrÇn ViÕt TuÊn (2007), Nghiªn cøu øng dông c«ng nghÖ GPS trong Tr¾c ®Þa c«ng tr×nh ë ViÖt Nam, LuËn ¸n tiÕn sÜ kü thuËt, Th­ viÖn tr­êng §¹i häc Má - §Þa chÊt, Hµ Néi. [7]. §Æng Nam Chinh (2005), B¸o c¸o kü thuËt ®o ®¹c tÝnh to¸n m¹ng l­íi GPS kiÓm tra tuyÕn n¨ng l­îng thuû ®iÖn A-V­¬ng – Qu¶ng Nam, Tr­êng ®¹i häc Má - §Þa chÊt, Hµ Néi. [8]. C«ng ty t­ vÊn x©y dùng ®iÖn 1 (2003), B¸o c¸o x©y dùng l­íi tam gi¸c thuû c«ng vµ l­íi thuû chuÈn thuû c«ng c«ng tr×nh thuû lîi – thuû ®iÖn Qu¶ng TrÞ. Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4865 1./ Ch­¬ng tr×nh tÝnh chuyÓn gi÷a c¸c hÖ to¹ ®é ph¼ng (L­íi GPS tuyÕn n¨ng l­îng thuû ®iÖn A - V­¬ng). $Debug $Large C Program: BAI TOAN TINH CHUYEN GIUA CAC HE TOA DO PHANG C Ho va ten: Vu Thi Ha _ Trac dia A_K48 Implicit Real *8(a-h,o-z) Character ff*50 Character tenCT *50,Tr*8 Character Khd(1000) *7,Ten(1000) *7 Common /V5/HV(1000),PTC(5000) Dimension Xc(1000),Yc(1000),Xm(1000),Ym(1000),Yp1(1000) Dimension Xp(1000),Yp(1000),X(1000),Y(1000),Xp1(1000) Write (*,'(a\)') ' Nhap ten FILE du lieu : ' Read (*, '(a)') ff Open (1, file= ff, status= 'old') Write (*,'(a\)') ' Nhap ten FILE ket qua : ' Read (*, '(a)') ff Open (4, file= ff, status= 'new') b1= 1.d0 Pi= 4.d0*dAtan(b1) Pi2= 2.d0*Pi Ro= 3600.d0*180.d0/Pi Read(1,*) Nst,Ntc Read(1,50) (Khd(I),I=1,Ntc) 50 Format(10A7) Do 1 I=1,Nst Read(1,*) I,Xc(I),Yc(I),Xm(I),Ym(I) 1 Continue Read(1,51) (Ten(I),I=1,Nst) 51 Format(10A7) Write(4,40) 40 Format(///22x,'KET QUA TINH CHUYEN TOA DO PHANG' *' THEO PHUONG PHAP HELMERT'// *30x,' BANG 1 - TOA DO CAC DIEM SONG TRUNG'/ *30x,' -----------------------------------'/ *10x,82('=')/ *10x,'| | | TOA DO TRONG HE CU |' *' TOA DO TRONG HE MOI |'/ *10x,'| STT | TEN DIEM |-------------------------------|' *'-------------------------------|'/ *10x,'| | | X(m) | Y(m) |' *' X(m) | Y(m) |'/ *10x,'|-----|----------|---------------|---------------|' *'---------------|---------------|') Do 41 I=1,Nst Write(4,42) I,Ten(I),Xc(I),Yc(I),Xm(I),Ym(I) 42 Format(10x,'|',I3,' | ',A7,' | ',F12.4,' | ',F12.4,' | ', *F12.4,' | ',F12.4,' |') 41 Continue Write(4,43) 43 Format(10x,82('=')/// *28x,' BANG 2 - KET QUA CAC THAM SO TINH CHUYEN'/ *28x,' ----------------------------------------'/ *28x,43('=')/ *28x,'| STT | THAM SO | GIA TRI |'/ *28x,'|-------|-----------|---------------------|') X1=0.d0 Y1=0.d0 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4866 Do 3 I=1,Nst X1=X1+Xc(I) Y1=Y1+Yc(I) 3 Continue Do 2 I=1,Ntc Read(1,*) I,X(I),Y(I) X1=X1+X(I) Y1=Y1+Y(I) 2 Continue Xo=X1/(Nst+Ntc) Yo=Y1/(Nst+Ntc) Do 4 I=1,Nst Xp(I)=Xc(I)-Xo Yp(I)=Yc(I)-Yo 4 Continue Do 29 I=1,Ntc Xp1(I)=X(I)-Xo Yp1(I)=Y(I)-Yo 29 Continue Ian=4 Do 33 I=1,(Ian+1) Do 34 J=I,(Ian+1) K=J*(J-1)/2+I PTC(K)=0.d0 34 Continue 33 Continue Do 5 I=1,Nst Do 6 L=1,(Ian+1) HV(L)=0.d0 6 Continue HV(1)=1.d0 HV(2)=0.D0 HV(3)= Xp(I) HV(4)=-Yp(I) HV(5)=-Xm(I) P=1.d0 Do 11 I1= 1,(Ian+1) Do 12 J1= I1,(Ian+1) K= J1*(J1-1)/2 + I1 PTC(K)= PTC(K)+HV(I1)*HV(J1)*P 12 Continue 11 Continue 5 Continue Do 8 I=1,Nst Do 7 L1=1,(Ian+1) HV(L1)=0.d0 7 Continue HV(1)=0.d0 HV(2)=1.D0 HV(3)= Yp(I) HV(4)= Xp(I) HV(5)=-Ym(I) P=1.d0 Do 23 I1= 1,(Ian+1) Do 24 J1= I1,(Ian+1) K= J1*(J1-1)/2 + I1 PTC(K)= PTC(K)+HV(I1)*HV(J1)*P 24 Continue 23 Continue 8 Continue Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4867 Call dGPTC(Ian) call dNDAO(Ian) KK =Ian*(Ian-1)/2+Ian X0=PTC(KK+1) Y0=PTC(KK+2) P1=PTC(KK+3) Q1=PTC(KK+4) afa=dMod(dAtan2(Q1,P1)+Pi2,Pi2) If(afa.Lt.0.d0) afa=afa+Pi2 af1=afa Call Drado(af1,Id,Iph,g,Ro) DM=dsqrt(P1*P1+Q1*Q1) Write(4,20) X0,Y0,Id,Iph,g,DM 20 Format(28x,'| 1 | Xo | ',F15.7,' |'/ * 28x,'| 2 | Yo | ',F15.7,' |'/ * 28x,'| 3 | GOC PHI | ',I3,I5,F7.2,' |'/ * 28x,'| 4 | M | ',F15.7,' |'/ * 28x,43('=')) Vx=0.d0 Vy=0.d0 Do 60 I=1,Nst Xt=X0+DM*Xp(I)*dcos(afa)-DM*Yp(I)*dsin(afa) Vx1=Xm(I)-Xt Vx=Vx+Vx1*Vx1 Yt=Y0+DM*Yp(I)*dcos(afa)+DM*Xp(I)*dsin(afa) Vy1=Ym(I)-Yt Vy=Vy+Vy1*Vy1 60 Continue VV=Vx+Vy Nk=2*Nst If(Nk.Eq.Ian) Goto 62 Dmo=dsqrt(VV/(2*Nst-4)) Goto 63 62 Write(4,61) 61 Format(28x,'SAI SO TRUNG PHUONG TRONG SO DON VI: 0.0000') Goto 64 63 Write(4,65) Dmo 65 Format(28x,'SAI SO TRUNG PHUONG TRONG SO DON VI:',f6.4) 64 Write(4,70) 70 Format(//30x,' BANG 3 - THANH QUA TOA DO TINH CHUYEN'/ * 30x,' -------------------------------------'/ *10x,82('=')/ *10x,'| | | TOA DO TRONG HE CU |' *' TOA DO TRONG HE MOI |'/ *10x,'| STT | TEN DIEM |-------------------------------|' *'-------------------------------|'/ *10x,'| | | X(m) | Y(m) |' *' X(m) | Y(m) |'/ *10x,'|-----|----------|---------------|---------------|' *'---------------|---------------|') Do 71 I=1,Ntc Xtc=X0+DM*Xp1(I)*dcos(afa)-DM*Yp1(I)*dsin(afa) Ytc=Y0+DM*Yp1(I)*dcos(afa)+DM*Xp1(I)*dsin(afa) Write(4,42) I,Khd(I),X(I),Y(I),Xtc,Ytc 71 Continue Write(4,75) 75 Format(10x,82('=')) Stop End Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4868 Subroutine Drado (Beta,Id,Iph,g,Ro) Implicit Real *8(a-h,o-z) Beta= Beta*Ro Id =Beta/3600.d0 Iph=(Beta-Id*3600d0)/60.d0 g =(Beta-Id*3600.d0-Iph*60.d0) Return End c CHUONG TRINH CON LAP PHUONG TRINH CHUAN Subroutine dPTC(Ian, P) Implicit Real *8(a-h, o-z) Common /V5/HV(1000),PTC(5000) Do 1 I= 1,(Ian+1) Do 2 J= I,(Ian+1) K= J*(J-1)/2 + I PTC(K)= PTC(K)+HV(I)*HV(J)*P 2 Continue 1 Continue Return End c Giai he phuong trinh chuan Subroutine dGPTC(iv) Implicit Real *8(a-h, o-z) Common /V5/HV(1000),PTC(5000) in= iv+1 Do 1 j= 2, in m = (j-1)*j/2 jn=j-1 Do 2 i=1, jn L=((i-1)*i)/2 kn= m+i Ln=L+i r=ptc(kn)/ptc(ln) Do 3 k= j, in n1=((k-1)*k)/2 nn=n1+i mn=n1+j ptc(mn)=ptc(mn)-ptc(nn)*r 3 Continue 2 Continue 1 Continue Do 4 i= iv, 1, -1 j=((i-1)*i)/2 k=iv L=((k+1)*k)/2 im=i+L D=-ptc(im) ii=iv*in/2+iv 5 If(k-i) 7, 7, 6 6 L=(k-1)*k/2 L1=(k+1)*k/2-k jn=L+i If(jn.le.L1) Goto 9 D=D-ptc(jn)*ptc(ii) 9 k=k-1 ii=ii-1 Goto 5 7 kn= i+j Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4869 ptc(im)=D/ptc(kn) 4 Continue 8 Return End c Chuong trinh con tinh nghich dao ma tran SUBROUTINE DNDAO(IV) IMPLICIT REAL *8(A-H,O-Z) Common /V5/HV(1000),PTC(5000) DIMENSION SS1(100) IN=IV+1 I=IN 504 I=I-1 J=(I-1)*I/2 K=IV IN=(K-1)*K/2+I D=1.D0 505 IF(K-I) 507,507,506 506 D=0.d0 JN=(K-1)*K/2+I 508 IM=IN+K-I SS1(K)=PTC(JN) 519 D=D-SS1(K)*PTC(IM) L=(K-2)*(K-1)/2 JN=I+L KN=I+J K=K-1 IF(KN-JN)508,505,505 507 KN=I+J PTC(IN)=D/PTC(KN) IF(KN-IN)511,512,511 511 K=IV 513 IN=IN-K+1 IF(KN-IN)530,531,530 531 D=1.D0 GOTO 533 530 D=0.d0 533 KK=(K-1)*K/2 II=K-1 K1=IV 558 L=(K1-1)*K1/2 IM=II+L IF(KK-IM)560,563,560 560 D=D-SS1(K1)*PTC(IM) K1=K1-1 GOTO 558 561 D=D-SS1(K1)*PTC(IM) K1=K1-1 IM=IM-1 563 IF(IM-IN)561,562,561 562 PTC(IN)=D/PTC(KN) K=K-1 IF(K-I)512,512,513 512 IF(I-1)504,509,504 509 RETURN END Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4870 KET QUA TINH CHUYEN TOA DO PHANG THEO PHUONG PHAP HELMERT BANG 1 - TOA DO CAC DIEM SONG TRUNG ----------------------------------- =========================================================================== | | TOA DO TRONG HE CU | TOA DO TRONG HE MOI | | TEN DIEM |-------------------------------|-------------------------------| | | X(m) | Y(m) | X(m) | Y(m) | |----------|---------------|---------------|---------------|---------------| | CNN-II | 1748638.5660 | 498539.6610 | 1750429.4300 | 784327.9410 | | CNN-IV | 1748249.8760 | 498340.2030 | 1750038.2110 | 784133.4530 | | CS-1 | 1745367.0010 | 501871.8930 | 1747200.5600 | 787701.7820 | | CS-III | 1745299.2740 | 501405.8830 | 1747126.9090 | 787236.6500 | | NM-II | 1744043.0570 | 501995.6750 | 1745878.2630 | 787842.4330 | ============================================================================ BANG 2 - KET QUA CAC THAM SO TINH CHUYEN ---------------------------------------- =========================================== | STT | THAM SO | GIA TRI | |-------|-----------|---------------------| | 1 | Xo | 1747963.4817323 | | 2 | Yo | 786381.3557567 | | 3 | GOC PHI | 359 16 7.01 | | 4 | M | 1.0000315 | =========================================== SAI SO TRUNG PHUONG TRONG SO DON VI: .010 BANG 3 - THANH QUA TOA DO TINH CHUYEN ------------------------------------- ============================================================================ | | TOA DO TRONG HE CU | TOA DO TRONG HE MOI | | TEN DIEM |-------------------------------|-------------------------------| | | X(m) | Y(m) | X(m) | Y(m) | |----------|---------------|---------------|---------------|---------------| | CNN-II | 1748638.5660 | 498539.6610 | 1750429.4311 | 784327.9417 | | CNN-IV | 1748249.8760 | 498340.2030 | 1750038.2144 | 784133.4553 | | CS-1 | 1745367.0010 | 501871.8930 | 1747200.5659 | 787701.7693 | | CS-III | 1745299.2740 | 501405.8830 | 1747126.8936 | 787236.6472 | | NM-II | 1744043.0570 | 501995.6750 | 1745878.2681 | 787842.4455 | | N1C-0 | 1748375.2820 | 498672.9350 | 1750167.8615 | 784464.5699 | | NM-P1 | 1744164.1700 | 502338.5190 | 1746003.7515 | 788183.7264 | | NM-P2 | 1744171.9030 | 501979.6840 | 1746006.9035 | 787824.8106 | =========================================================================== Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4871 2./ Ch­¬ng tr×nh tÝnh chuyÓn to¹ ®é c¸c ®iÓm ®o GPS vÒ hÖ to¹ ®é c«ng tr×nh (l­íi khèng chÕ thi c«ng khu c«ng nghiÖp Yªn Phong – B¾c Ninh). $debug $large c CHUONG TRINH TINH CHUYEN CAC DIEM DO GPS VE HE TOA DO THI CONG c Ho va ten : Vu Thi Ha c Lop : Trac dia A - K48 Implicit Real *8(a-h,o-z) Character ff*30 Character Tenct*100, td(1000)*10 Common /V5/HV(1000),PTC(5000) Dimension X(1000),Y(1000),Z(1000) Dimension X1(1000),Y1(1000),Z1(1000) Dimension GB(1000),GL(1000),H1(1000),Xtc(1000),Ytc(1000) Dimension Xm(1000),Ym(1000),Xp(1000),Yp(1000) Write (*,'(a\)') ' Nhap ten FILE so lieu : ' Read (*, '(a)') ff Open (1, file= ff, status= 'old') Write (*,'(a\)') ' Nhap ten FILE ket qua : ' Read (*, '(a)') ff Open (4, file= ff, status= 'new') b1= 1.d0 Pi= 4.d0*dAtan(b1) Pi2= 2.d0*Pi Ro= 3600.d0*180.d0/Pi Read (1,'(a)') Tenct Read(1,*) Nd,Nst Write(4,15) Tenct,Nd 15 Format(///40x,'KET QUA TINH TOAN'/ * 22x,A80/ * 40x,'=======o0o======='/// * 30x,'Tong so diem can tinh chuyen: ',I3/) Read(1,50)(Td(i),i=1,Nd) 50 Format(10A7) Do 1 I=1,Nd Read(1,*) X(I),Y(I),Z(I) 1 Continue Write(4,300) 300 Format(//30x,'I. HE TOA DO VUONG GOC KHONG GIAN'/ * 30x,'================================='/ *15x,64('=')/ *15x,'| | | TOA DO VUONG' *' GOC KHONG GIAN |'/ *15x,'| STT | TEN |---------------------' *'--------------------------|'/ *15x,'| | DIEM | X (m) | Y (m) |' *' Z (m) |'/ *15x,'|-----|--------|---------------|---------------|' *'---------------|') Do 301 I=1,Nd Write(4,302) I,Td(I),X(I),Y(I),Z(I) 302 Format(15x,'|',I3,' | ',A7,'|',F14.4,' |',F14.4,' |', *F14.4,' |') 301 Continue Write(4,303) 303 Format(15x,64('=')/) Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4872 Write(4,360) 360 Format(//36x,'II. HE TOA DO TRAC DIA'/ * 36x,'======================'/ *16x,61('=')/ *16x,'| | | TOA DO TRAC ' *' DIA |'/ *16x,'| STT | TEN |---------------------' *'-----------------------|'/ *16x,'| | DIEM | B | L |' *' H |'/ *16x,'| | | (o , ") | (o , ") |' *' (m) |'/ *16x,'|-----|--------|----------------|---------------|' *'-----------|') e2=0.006693421623 a=6378137 gh=0.00000000000001 Do 2 I=1,Nd P=Dsqrt(X(I)*X(I)+Y(I)*Y(I)) Bo=DAtan(Z(I)/P/(1-e2)) 6 T1=Dcos(Bo) T2=Dsin(Bo) Can=Dsqrt(1-e2*T2*T2) DNo=a/Can H=P/Dcos(Bo)-DNo T3=1-e2*DNo/(DNo+H) PT=P*T3 B=dMod(dAtan2(Z(I),PT)+Pi2,Pi2) If(B.Lt.0.d0) B=B+Pi2 e=ABS(B-Bo) If(e.LE.(gh/Ro)) Goto 5 Bo=B Goto 6 5 DL=DAtan(Y(I)/X(I)) DL=DL+Pi GL(I)=DL GB(I)=B HH=P/Dcos(B)-a/Dsqrt(1-e2*(dsin(B)**2)) H1(I)=HH Call Drado(B,Id1,Iph1,g1,Ro) Call Drado(DL,Id2,Iph2,g2,Ro) Write(4,80) I,Td(I),Id1,Iph1,g1,Id2,Iph2,g2,HH 80 Format(16x,'|',I3,' | ',A7,'| ',I4,I3,F6.2,' | ', *I4,I3,F6.2,' |',F10.3' |') 2 Continue Write(4,56) 56 Format(16x,61('=')///) TB=0.d0 TL=0.d0 TH=0.d0 Do 30 I=1,Nd TB=TB+GB(I) TL=TL+GL(I) TH=TH+H1(I) 30 Continue T1=TB/Nd T2=TL/Nd T3=TH/Nd FN=a/Dsqrt(1-e2*(dsin(T1)**2)) Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4873 a1=(FN+T3)*Dcos(T1)*Dcos(T2) a2=(FN+T3)*Dcos(T1)*Dsin(T2) a3=(FN*(1-e2)+T3)*Dsin(T1) Do 31 I=1,Nd q11=dcos(GB(I)) q12=dsin(GB(I)) q21=dcos(GL(I)) q22=dsin(GL(I)) x1(I)=-((X(I)-a1)*q12*q21+(Y(I)-a2)*q12*q22-(Z(I)-a3)*q11) y1(I)=-((X(I)-a1)*q22-(Y(I)-a2)*q21) z1(I)= (X(I)-a1)*q11*q21+(Y(I)-a2)*q11*q22+(Z(I)-a3)*q12 31 Continue M1=Nd-Nst Do 60 I=M1+1,Nd Read(1,*) I,Xm(I),Ym(I) 60 Continue X11=0.d0 Y11=0.d0 Do 104 I=1,Nd X11=X11+x1(I) Y11=Y11+y1(I) 104 Continue Xo=X11/Nd Yo=Y11/Nd Do 106 I=1,Nd Xp(I)=x1(I)-Xo Yp(I)=y1(I)-Yo 106 Continue Ian=4 Do 108 I=1,(Ian+1) Do 109 J=I,(Ian+1) K=J*(J-1)/2+I PTC(K)=0.d0 109 Continue 108 Continue Do 110 I=M1+1,Nd Do 111 L=1,(Ian+1) HV(L)=0.d0 111 Continue HV(1)=1.d0 HV(2)=0.D0 HV(3)= Xp(I) HV(4)=-Yp(I) HV(5)=-Xm(I) P=1.d0 Do 112 I1= 1,(Ian+1) Do 113 J1= I1,(Ian+1) K= J1*(J1-1)/2 + I1 PTC(K)= PTC(K)+HV(I1)*HV(J1)*P 113 Continue 112 Continue 110 Continue Do 114 I=M1+1,Nd Do 115 L1=1,(Ian+1) HV(L1)=0.d0 115 Continue HV(1)=0.d0 HV(2)=1.D0 HV(3)= Yp(I) Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4874 HV(4)= Xp(I) HV(5)=-Ym(I) P=1.d0 Do 116 I1= 1,(Ian+1) Do 117 J1= I1,(Ian+1) K= J1*(J1-1)/2 + I1 PTC(K)= PTC(K)+HV(I1)*HV(J1)*P 117 Continue 116 Continue 114 Continue Call dGPTC(Ian) call dNDAO(Ian) KK =Ian*(Ian-1)/2+Ian X0=PTC(KK+1) Y0=PTC(KK+2) P1=PTC(KK+3) Q1=PTC(KK+4) afa=dMod(dAtan2(Q1,P1)+Pi2,Pi2) If(afa.Lt.0.d0) afa=afa+Pi2 af1=afa Call Drado(af1,Id,Iph,g,Ro) DM=dsqrt(P1*P1+Q1*Q1) Vx=0.d0 Vy=0.d0 Do 69 I=M1+1,Nd Xt=X0+DM*Xp(I)*dcos(afa)-DM*Yp(I)*dsin(afa) Vx1=Xm(I)-Xt Vx=Vx+Vx1*Vx1 Yt=Y0+DM*Yp(I)*dcos(afa)+DM*Xp(I)*dsin(afa) Vy1=Ym(I)-Yt Vy=Vy+Vy1*Vy1 69 Continue Do 150 I=1,Nd Xtc(I)=X0+DM*Xp(I)*dcos(afa)-DM*Yp(I)*dsin(afa) Ytc(I)=Y0+DM*Yp(I)*dcos(afa)+DM*Xp(I)*dsin(afa) 150 Continue Write(4,330) 330 Format(36x,'III. HE TOA DO DIA DIEN'/ * 36x,'======================='/ *16x,64('=')/ *16x,'| | | TOA DO DIA DIEN ' *' |'/ *16x,'| STT | TEN |--------------------------------------' *'---------|'/ *16x,'| | DIEM | x (m) | y (m) |' *' z (m) |'/ *16x,'|-----|--------|---------------|---------------|' *'---------------|') Do 331 I=1,Nd Write(4,332) I,Td(I),X1(I),Y1(I),Z1(I) 332 Format(16x,'|',I3,' | ',A7,'| ',F10.4,' | ',F10.4, *' | ',F9.4,' |') 331 Continue Write(4,335) 335 Format(16x,64('=')///) Write(4,430) 430 Format(35x,'IV. HE TOA DO THI CONG'/ * 35x,'======================'/ *22x,48('=')/ Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4875 *22x,'| | | TOA DO THI CONG |'/ *22x,'| STT | TEN |-------------------------------|'/ *22x,'| | DIEM | X (m) | Y (m) |'/ *22x,'|-----|--------|---------------|---------------|') Do 531 I=1,Nd Write(4,532) I,Td(I),Xtc(I),Ytc(I) 532 Format(22x,'|',I3,' | ',A7,'| ',F12.4,' | ', *F12.4,' | ') 531 Continue Write(4,535) 535 Format(22x,48('=')/) Stop End Subroutine Drado (Beta,Id,Iph,g,Ro) Implicit Real *8(a-h,o-z) Beta= Beta*Ro Id =Beta/3600.d0 Iph=(Beta-Id*3600d0)/60.d0 g =(Beta-Id*3600.d0-Iph*60.d0) Return End c CHUONG TRINH CON LAP PHUONG TRINH CHUAN Subroutine dPTC(Ian, P) Implicit Real *8(a-h, o-z) Common /V5/HV(1000),PTC(5000) Do 1 I= 1,(Ian+1) Do 2 J= I,(Ian+1) K= J*(J-1)/2 + I PTC(K)= PTC(K)+HV(I)*HV(J)*P 2 Continue 1 Continue Return End c Giai he phuong trinh chuan Subroutine dGPTC(iv) Implicit Real *8(a-h, o-z) Common /V5/HV(1000),PTC(5000) in= iv+1 Do 1 j= 2, in m = (j-1)*j/2 jn=j-1 Do 2 i=1, jn L=((i-1)*i)/2 kn= m+i Ln=L+i r=ptc(kn)/ptc(ln) Do 3 k= j, in n1=((k-1)*k)/2 nn=n1+i mn=n1+j ptc(mn)=ptc(mn)-ptc(nn)*r 3 Continue 2 Continue 1 Continue Do 4 i= iv, 1, -1 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4876 j=((i-1)*i)/2 k=iv L=((k+1)*k)/2 im=i+L D=-ptc(im) ii=iv*in/2+iv 5 If(k-i) 7, 7, 6 6 L=(k-1)*k/2 L1=(k+1)*k/2-k jn=L+i If(jn.le.L1) Goto 9 D=D-ptc(jn)*ptc(ii) 9 k=k-1 ii=ii-1 Goto 5 7 kn= i+j ptc(im)=D/ptc(kn) 4 Continue 8 Return End c Chuong trinh con tinh nghich dao ma tran SUBROUTINE DNDAO(IV) IMPLICIT REAL *8(A-H,O-Z) Common /V5/HV(1000),PTC(5000) DIMENSION SS1(100) IN=IV+1 I=IN 504 I=I-1 J=(I-1)*I/2 K=IV IN=(K-1)*K/2+I D=1.D0 505 IF(K-I) 507,507,506 506 D=0.d0 JN=(K-1)*K/2+I 508 IM=IN+K-I SS1(K)=PTC(JN) 519 D=D-SS1(K)*PTC(IM) L=(K-2)*(K-1)/2 JN=I+L KN=I+J K=K-1 IF(KN-JN)508,505,505 507 KN=I+J PTC(IN)=D/PTC(KN) IF(KN-IN)511,512,511 511 K=IV 513 IN=IN-K+1 IF(KN-IN)530,531,530 531 D=1.D0 GOTO 533 530 D=0.d0 533 KK=(K-1)*K/2 II=K-1 K1=IV 558 L=(K1-1)*K1/2 IM=II+L Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4877 IF(KK-IM)560,563,560 560 D=D-SS1(K1)*PTC(IM) K1=K1-1 GOTO 558 561 D=D-SS1(K1)*PTC(IM) K1=K1-1 IM=IM-1 563 IF(IM-IN)561,562,561 562 PTC(IN)=D/PTC(KN) K=K-1 IF(K-I)512,512,513 512 IF(I-1)504,509,504 509 RETURN END Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4878 KET QUA TINH TOAN TINH CHUYEN CAC DIEM DO GPS VE HE TOA DO THI CONG =======o0o======= Tong so diem can tinh chuyen: 20 I. HE TOA DO VUONG GOC KHONG GIAN ================================= ================================================================ | | | TOA DO VUONG GOC KHONG GIAN | | STT | TEN |-----------------------------------------------| | | DIEM | X (m) | Y (m) | Z (m) | |-----|--------|---------------|---------------|---------------| | 1 | 104511 | -1637975.5132 | 5718877.5746 | 2292679.3997 | | 2 | 105486 | -1640916.3957 | 5719051.7850 | 2290156.1370 | | 3 | IV14 | -1638669.8897 | 5718616.0022 | 2292836.6566 | | 4 | IVYP1 | -1640728.4742 | 5718959.9023 | 2290528.4510 | | 5 | IVYP2 | -1640339.6048 | 5718942.0091 | 2290857.3628 | | 6 | IVYP4 | -1638981.2098 | 5718736.8111 | 2292313.1677 | | 7 | YP1 | -1637824.5796 | 5718989.5548 | 2292508.5914 | | 8 | YP2 | -1637573.4493 | 5719155.4256 | 2292276.0888 | | 9 | YP3 | -1637708.5128 | 5719256.0365 | 2291940.9737 | | 10 | YP4 | -1638067.8775 | 5719170.3412 | 2291900.8996 | | 11 | YP5 | -1638229.5098 | 5719131.4450 | 2291882.4235 | | 12 | YP6 | -1638740.6065 | 5719011.1266 | 2291812.1393 | | 13 | YP7 | -1638294.8952 | 5718892.8216 | 2292414.7043 | | 14 | YP8 | -1639258.4543 | 5718690.7380 | 2292233.2048 | | 15 | YP9 | -1639575.0585 | 5718635.3487 | 2292142.8653 | | 16 | YP10 | -1639888.6480 | 5718677.8311 | 2291815.4892 | | 17 | YP11 | -1640232.9191 | 5718623.1514 | 2291706.1491 | | 18 | YP12 | -1640277.1302 | 5718770.7519 | 2291308.7924 | | 19 | IV15 | -1639729.3281 | 5718895.4959 | 2291401.6398 | | 20 | IV16 | -1639219.1804 | 5718933.3868 | 2291671.2684 | ================================================================ Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4879 III. HE TOA DO DIA DIEN ======================= ================================================================ | | | TOA DO DIA DIEN | | STT | TEN |-----------------------------------------------| | | DIEM | x (m) | y (m) | z (m) | |-----|--------|---------------|---------------|---------------| | 1 | 104511 | 923.1233 | -1085.6856 | -1.5191 | | 2 | 105486 | -1783.0525 | 1693.2512 | -1.6876 | | 3 | IV14 | 1091.5478 | -346.0780 | -.8344 | | 4 | IVYP1 | -1385.2392 | 1537.9716 | 2.0367 | | 5 | IVYP2 | -1033.5606 | 1169.1388 | 4.7483 | | 6 | IVYP4 | 530.4787 | -80.0723 | -2.0369 | | 7 | YP1 | 739.9602 | -1261.6260 | -1.6737 | | 8 | YP2 | 490.5034 | -1548.7266 | -1.5136 | | 9 | YP3 | 129.6540 | -1446.5753 | 2.0969 | | 10 | YP4 | 86.3335 | -1077.5031 | 2.9672 | | 11 | YP5 | 66.5451 | -911.4082 | 2.8931 | | 12 | YP6 | -8.0381 | -386.9397 | .8254 | | 13 | YP7 | 639.2539 | -782.8354 | -1.6813 | | 14 | YP8 | 444.3184 | 199.1436 | -1.0338 | | 15 | YP9 | 347.7889 | 518.7604 | -1.9899 | | 16 | YP10 | -3.4676 | 808.4967 | -1.6951 | | 17 | YP11 | -120.7425 | 1154.5001 | -1.6938 | | 18 | YP12 | -546.9230 | 1156.3032 | -1.6838 | | 19 | IV15 | -449.0833 | 595.3557 | 2.8034 | | 20 | IV16 | -160.0247 | 94.5289 | 3.1305 | ================================================================ IV. HE TOA DO THI CONG Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4880 ====================== ================================================ | | | TOA DO THI CONG | | STT | TEN |-------------------------------| | | DIEM | X (m) | Y (m) | |-----|--------|---------------|---------------| | 1 | 104511 | 2345742.1692 | 550095.5290 | | 2 | 105486 | 2343044.8562 | 552882.9950 | | 3 | IV14 | 2345912.9359 | 550834.5886 | | 4 | IVYP1 | 2343442.1697 | 552726.4568 | | 5 | IVYP2 | 2343792.6721 | 552356.5157 | | 6 | IVYP4 | 2345352.7208 | 551102.3685 | | 7 | YP1 | 2345558.4517 | 549920.1727 | | 8 | YP2 | 2345308.0891 | 549633.8685 | | 9 | YP3 | 2344947.5704 | 549737.1621 | | 10 | YP4 | 2344905.4211 | 550106.3648 | | 11 | YP5 | 2344886.1598 | 550272.5194 | | 12 | YP6 | 2344813.2410 | 550797.2145 | | 13 | YP7 | 2345459.2656 | 550399.2737 | | 14 | YP8 | 2345267.4475 | 551381.8524 | | 15 | YP9 | 2345171.9333 | 551701.7692 | | 16 | YP10 | 2344821.6022 | 551992.6140 | | 17 | YP11 | 2344705.4267 | 552338.9827 | | 18 | YP12 | 2344279.2599 | 552342.1372 | | 19 | IV15 | 2344375.3190 | 551780.8900 | | 20 | IV16 | 2344662.7840 | 551279.1560 | ================================================ Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4881 3./ Ch­¬ng tr×nh tÝnh chuyÓn to¹ ®é vÒ ®é cao khu vùc (L­íi tam gi¸c thuû c«ng thuû ®iÖn Qu¶ng TrÞ). $debug $large c CHUONG TRINH TINH CHUYEN TOA DO VE DO CAO KHU VUC c Ho va ten : Vu Thi Ha c Lop : Trac dia A - K48 Implicit Real *8(a-h, o-z) Character Tenct*100, td(1000)*10,ff*100 Dimension X(1000),Y(1000),X1(1000),Y1(1000) Write (*,'(a\)') ' Nhap ten FILE du lieu : ' Read (*, '(a)') ff Open (1, file= ff, status= 'old') Write (*,'(a\)') ' Nhap ten FILE ket qua : ' Read (*, '(a)') ff Open (4, file= ff, status= 'new') b1= 1.d0 Pi= 4.d0*dAtan(b1) Pi2= 2.d0*Pi Ro= 3600.d0*180.d0/Pi Read(1,*) Nd,Ho,H Read(1,50) (td(I),I=1,Nd) 50 Format(10A7) IHo=INT(Ho) IH=INT(H) Write(4,20) IHo,IH 20 Format(///32x,' KET QUA TINH CHUYEN TOA DO'/ * 32x,' GIUA CAC MAT CHIEU '/ * 32x,'============oOo============'/ * 32x,' ========O======== '/// * 34x,' THANH QUA TINH CHUYEN '/ * 34x,' =========***========= '/ *5x,84('=')/ *5x,'| | | DO CAO BAN DAU Ho = ',I4,' |' *' DO CAO TB KHU VUC H = ',I4,' |'/ *5x,'| STT | TEN DIEM |-------------------------------|' *'---------------------------------|'/ *5x,'| | | X (m) | Y (m) |' *' X (m) | Y (m) |'/ *5x,'|-----|----------|---------------|---------------|' *'----------------|----------------|') DH=H-Ho xo=0.d0 yo=0.d0 Do 1 i=1,Nd Read(1,*) I,X(I),Y(I) xo=xo+X(I) yo=yo+Y(I) 1 Continue xo=xo/Nd yo=yo/Nd Rm=6370000 dK=(Rm+DH)/Rm Do 2 I=1,Nd a2=(X(I)-xo)*dK b2=(Y(I)-yo)*dK x1(I)=xo+a2 y1(I)=yo+b2 Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4882 Write(4,30) I,td(I),X(I),Y(I),x1(I),y1(I) 30 Format(5x,'|',I3,' | ',A7,' |',F14.4,' |',f14.4,' |', *F14.4,' |',f14.4,' |') 2 Continue Write(4,31) 31 Format(5x,84('=')) Stop End c DOI GOC TU RADIAN RA DO Subroutine Drado(Beta,Id,Iph,g,Ro) Implicit Real *8(a-h,o-z) Beta= Beta*Ro Id =Beta/3600d0 Iph =(Beta-Id*3600d0)/60d0 g =(Beta-Id*3600d0-Iph*60d0) Return End Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4883 KET QUA TINH CHUYEN TOA DO GIUA CAC MAT CHIEU ============oOo============ ========O======== THANH QUA TINH CHUYEN =========***========= ============================================================================= | | DO CAO BAN DAU Ho = 0 | DO CAO TB KHU VUC H = 300 | | TEN DIEM |-------------------------------|---------------------------------| | | X (m) | Y (m) | X (m) | Y (m) | |----------|---------------|---------------|----------------|----------------| | TC-01 | 1846589.5950 | 548334.2430 | 1846589.7203 | 548334.1559 | | TC-02 | 1846339.2910 | 547838.3410 | 1846339.4045 | 547838.2305 | | TC-03 | 1845487.5230 | 547706.8360 | 1845487.5964 | 547706.7193 | | TC-04 | 1845490.8640 | 548315.9480 | 1845490.9376 | 548315.8600 | | TC-05 | 1846056.3760 | 549029.6830 | 1846056.4762 | 549029.6286 | | TC-06 | 1844977.1940 | 548547.8650 | 1844977.2434 | 548547.7879 | | TC-07 | 1845122.7880 | 549320.3500 | 1845122.8442 | 549320.3093 | | TC-08 | 1844063.5990 | 549058.3680 | 1844063.6054 | 549058.3150 | | TC-09 | 1844532.3000 | 549949.8820 | 1844532.3284 | 549949.8709 | | TC-10 | 1843631.5280 | 549908.0750 | 1843631.5140 | 549908.0620 | | TC-11 | 1844231.6430 | 550726.3250 | 1844231.6573 | 550726.3505 | | TC-12 | 1842747.7150 | 550387.4940 | 1842747.6594 | 550387.5036 | | TC-13 | 1843253.4610 | 551117.9080 | 1843253.4292 | 551117.9520 | | TC-14 | 1842123.2550 | 550981.4380 | 1842123.1700 | 550981.4755 | | TC-15 | 1842564.4610 | 551982.6960 | 1842564.3968 | 551982.7807 | | TC-16 | 1841871.5290 | 552210.5650 | 1841871.4321 | 552210.6604 | | TC-17 | 1842343.0120 | 552662.6350 | 1842342.9373 | 552662.7517 | | TC-18 | 1841685.1780 | 553078.3140 | 1841685.0723 | 553078.4503 | | TC-19 | 1841533.4340 | 552349.2210 | 1841533.3212 | 552349.3229 | ============================================================================== Khoa Tr¾c ®Þa §å ¸n tèt nghiÖp Sinh viªn: Vò ThÞ Hµ Líp: Tr¾c ®Þa A – K4884 PHô LôC