Luận văn Công nghệ đo ảnh trên cơ sở ảnh chụp từ máy chụp ảnh số phổ thông

Ðo dµi trong mét kho¶ng thêi gian nhÊt ®Þnh nªn c¸c sè liÖu kh«ng ®−îc thèng nhÊt trong cïng mét thêi ®iÓm, g©y nªn khã kh¨n cho c«ng t¸c tÝnh to¸n khèi l−îng. §Ó cã c¬ së hç trî ng−êi tr¾c ®Þa má gi¶m bít c«ng viÖc ngo¹i nghiÖp, ®ång thêi t¨ng ®é chÝnh x¸c tÝnh to¸n khèi l−îng khai th¸c, chóng t«i tiÕn hµnh nhiÒu c«ng t¸c thùc nghiÖm ®o vÏ ®Þa h×nh dïng ¶nh chôp tõ m¸y ¶nh sè phæ th«ng. Chóng t«i ®· tiÕn hµnh chôp ¶nh trong ba ®ît vµo c¸c ngµy 11-10- 2003, ngµy 02-12-2003 vµ ngµy 20-02-2004. Khu vùc chän ®Ó chôp lµ khu §«ng Nam (khu1) vµ khu §«ng B¾c (khu 2) cña má Cäc S¸u. 4.2.2. C«ng t¸c chôp ¶nh vµ c«ng t¸c tr¾c ®Þa ngo¹i nghiÖp Chóng t«i chän tr¹m chôp ë bèn ®iÓm ®−îc ch«n mèc trªn bê tÜnh. Trªn c¬ së m¹ng l−íi khèng chÕ má, bèn ®iÓm nµy ®−îc truyÒn täa ®é tr¾c ®Þa ®Ó hßa chung hÖ täa ®é khu má vµ lµ c¬ së ®Ó x¸c ®Þnh täa ®é c¸c b¶ng ng¾m t¹i c¸c ®iÓm kiÓm tra trong m« h×nh lËp thÓ cña ¶nh chôp. C¸c ®iÓm tr¹m chôp nµy ®−îc thiÕt kÕ ®Ó cã thÓ ®o chôp ¶nh tíi kho¶ng c¸ch 200m-700m tõ ®−êng ®¸y chôp. Ngµy 11-10-2003, chóng t«i tiÕn hµnh chôp ¶nh ë hai khu vùc: - Khu §«ng Nam má chôp 4 ¶nh, mçi tr¹m chôp 2 tÊm cã d¹ng chôp gÇn ®óng lµ th¼ng vµ th¼ng nghiªng, t¹o thµnh hai cÆp ¶nh lËp thÓ lµ 637-643 vµ 638-644 cña cïng mét khu vùc ®Ó cã ®iÒu kiÖn kiÓm tra lÉn nhau. - Khu §«ng B¾c má chôp 4 tÊm, mçi tr¹m chôp 2 tÊm ë d¹ng chôp gÇn ®óng lµ th¼ng vµ giao nhau, x©y dùng ®−îc hai m« h×nh lËp thÓ ®éc lËp nhau lµ 647-645 vµ 648-646. Trªn mçi khu ®o bè trÝ 10 ®iÓm kiÓm tra thµnh 3 hµng dïng ®Þnh h−íng vµ võa ®Ó kiÓm tra m« h×nh. Trong ®ît nµy chóng t«i còng tiÕn hµnh ®o vÏ trùc tiÕp ngoµi thùc ®Þa b»ng m¸y toµn ®¹c ®iÖn tö ë c¸c ®−êng ®Æc tr−ng cña má 83 nh− ch©n tÇng, mÐp taly. C¸c ®iÓm nµy phôc vô viÖc ®¸nh gi¸ ®é chÝnh x¸c kÕt qu¶ ®o ¶nh ®−îc thùc hiÖn ®éc lËp víi qu¸ tr×nh chôp ¶nh. Ngµy 02-12-2003, chóng t«i thùc hiÖn t−¬ng tù ®ît ®Çu nh−ng do nhÇm lÉn ë mét sè kh©u ®o tr¾c ®Þa vµo c¸c b¶ng ng¾m nªn kÕt qu¶ ®ît nµy kh«ng ®−îc sö dông. §©y còng lµ mét bµi häc cÇn rót kinh nghiÖm cña chóng t«i. Ngµy 20-02-2004, chóng t«i vÉn tiÕn hµnh chôp ¶nh vµ ®o c¸c b¶ng ng¾m ë cïng mét khu vùc nh− ®ît ®Çu, ®Ó cã ®iÒu kiÖn tÝnh to¸n sù thay ®æi ®Þa h×nh trªn c¬ së b×nh ®å vµ khèi l−îng khai th¸c: - Khu §«ng Nam má t¹i 2 tr¹m chôp, mçi tr¹m chôp 2 tÊm cã d¹ng chôp th¼ng vµ xiªn ph¶i, t¹o thµnh hai cÆp ¶nh lËp thÓ lµ 823-819 vµ 826-820 kiÓm tra lÉn nhau. - Khu §«ng B¾c má chôp 4 tÊm, mçi tr¹m chôp 2 tÊm ë d¹ng chôp gÇn ®óng lµ th¼ng vµ xiªn ph¶i, x©y dùng ®−îc hai m« h×nh lËp thÓ ®éc lËp nhau lµ 837-839 vµ 838-840 kiÓm tra lÉn nhau. Trong ®ît nµy sè l−îng ®iÓm kiÓm tra ®−îc chóng t«i bè trÝ dµy ®Æc h¬n, ®Ó cã sè l−îng ®iÓm kh«ng tham gia x©y dùng m« h×nh mµ tham gia ®¸nh gi¸ ®é chÝnh x¸c m« h×nh Ýt nhÊt lµ 7 ®iÓm. C¸c tÊm ¶nh nh×n toµn c¶nh trong hai ®ît cña khu thùc nghiÖm §«ng Nam (khu1) ë phô lôc 9, khu §«ng B¾c (khu 2) ë phô lôc 10. 4.2.3. C«ng t¸c ®o vÏ xö lý ¶nh trong phßng C«ng t¸c ®o vÏ xö lý ¶nh ®−îc tiÕn hµnh trªn tr¹m ®o vÏ ¶nh sè ImageStation ë Côc §o ®¹c vµ B¶n ®å. T¸c gi¶ lµ ng−êi trùc tiÕp x©y dùng vµ kiÓm tra m« h×nh lËp thÓ, c«ng t¸c ®o vÏ chi tiÕt ®−îc thùc hiÖn bëi ng−êi t¸c nghiÖp cã tr×nh ®é ®o c¾t lËp thÓ cao. Cô thÓ c¸c b−íc tiÕn hµnh theo quy tr×nh c«ng nghÖ ®· tr×nh bµy ë h×nh 3.9. - X©y dùng vµ kiÓm tra m« h×nh lËp thÓ ®−îc thùc hiÖn b»ng modul ISDM. Khi x©y dùng m« h×nh lËp thÓ chóng t«i sö dông c¸c ®iÓm kiÓm tra cã täa ®é tr¾c ®Þa chÝnh x¸c ®Ó ®Þnh h−íng m« h×nh (control point), mét sè ®iÓm kh«ng tham gia x©y dùng m« h×nh th× dïng lµm ®iÓm kiÓm tra ®é chÝnh x¸c 84 (check point). Trong c«ng ®o¹n nµy chóng t«i còng ®¸nh gi¸ sù hiÖu qu¶ cña c¸c hÖ sè mÐo h×nh kiÓm ®Þnh tham gia x©y dùng m« h×nh b»ng c¸ch so s¸nh ®é chÝnh x¸c cña c¸c ®iÓm kiÓm tra trong 1 m« h×nh víi 2 tr−êng hîp cã vµ kh«ng cã hÖ sè mÐo h×nh. Cïng víi c¸ch kiÓm tra nµy, chóng t«i cßn sö dông mét c¸ch kiÓm tra ®é chÝnh x¸c m« h×nh n÷a khi sö dông c¸c ®iÓm chi tiÕt ®−îc ®o b»ng m¸y toµn ®¹c ®iÖn tö trùc tiÕp t¹i thùc ®Þa sau ®ã ®−a lªn m« h×nh sè cña ¶nh trªn c¬ së täa ®é XY råi so s¸nh ®é chªnh gi÷a Z ®o ¶nh vµ Z ®o ngo¹i nghiÖp, c¸c kÕt qu¶ nµy thÓ hiÖn ë phô lôc 5. C¸c sè liÖu thèng kª vÒ x©y dùng vµ kiÓm tra c¸c m« h×nh lËp thÓ ®−îc tr×nh bµy v¾n t¾t ë b¶ng 3.6. Cô thÓ c¸c sè liÖu nµy ®−îc xuÊt ra d−íi d¹ng v¨n b¶n cña phÇn mÒm Integraph ®−îc giíi thiÖu chi tiÕt ë c¸c phô lôc 1, 2, 3, 4. - C«ng t¸c ®o vÏ chi tiÕt ®−îc thùc hiÖn b»ng modul ISSD. Khi ®o vÏ trªn m« h×nh lËp thÓ cã tÝnh bao qu¸t vµ c¾t lËp thÓ s¸t c¸c ch©n tÇng, ®Ønh tÇng, bæ sung thªm c¸c ®−êng chuyÓn biÕn ®é cao ë g−¬ng tÇng, taly, mÆt tÇng... víi sè l−îng ®iÓm ®o kh«ng h¹n chÕ, nh÷ng chç ng−êi ®i mia, ®i m¸y kh«ng tíi ®−îc, b×nh ®å ®Þa h×nh khu má lé thiªn ®−îc x©y dùng mét c¸ch chÝnh x¸c cho tÝnh to¸n khèi l−îng khai th¸c. Khi tÝnh to¸n khèi l−îng khai th¸c trªn c¬ së c¸c ®iÓm ®o cã mËt ®é dµy ®Æc, chóng t«i dïng phÇn mÒm MGE Terrain Analys x©y dùng DTM theo l−íi TIN cña khu vùc tÝnh gi÷a c¸c lÇn ®o vÏ kh¸c nhau, chång phñ lªn nhau, khai b¸o biªn tÝnh khèi l−îng, ch−¬ng tr×nh sÏ cho ngay khèi l−îng khai th¸c theo ®¬n vÞ mÐt khèi (m3). C¸c b×nh ®å ®Þa h×nh khu má ®−îc biªn tËp s¬ bé vµ in Ên ë phô lôc 11, 12, 13, 14. C¸c DTM t−¬ng øng ®−îc in ë phô lôc 15, 16, 17, 18. TÝnh to¸n khèi l−îng khai th¸c ®−îc tr×nh bµy ë phô lôc 6. Ngoµi ra trong c«ng ®o¹n nµy, v× ®· cã thªm m« h×nh lËp thÓ ®éc lËp kh¸c nh−ng cïng khu vùc nªn chóng t«i tiÕp tôc kiÓm tra ®é chÝnh x¸c cña m« h×nh lËp thÓ còng nh− kÕt qu¶ ®o vÏ b»ng c¸ch: sau khi ®o vÏ xong m« h×nh thø nhÊt th× chång xÕp kÕt qu¶ ®o vÏ nµy lªn m« h×nh cßn l¹i xem xÐt sù b¸m s¸t bÒ mÆt m« h×nh. NhËn thÊy r»ng, c¸c kÕt qu¶ nµy rÊt b¸m s¸t m« h×nh ®éc lËp cßn l¹i. §©y còng lµ mét c¸ch kiÓm tra ®¬n gi¶n mµ rÊt hiÖu qu¶. 85 4.2.4. NhËn xÐt vÒ c¸c kÕt qu¶ ®o vÏ ®Þa h×nh - KÕt qu¶ x©y dùng vµ kiÓm tra m« h×nh lËp thÓ ®¹t ®−îc ®é chÝnh x¸c ®¹t yªu cÇu. §é chÝnh x¸c nµy n»m trong giíi h¹n thµnh lËp b×nh ®å ®Þa h×nh tû lÖ 1:2000 víi kho¶ng c¸ch 500m. + §ît 10-11-2003: T¹i c¸c m« h×nh 637-643, 638-644 cã 4 ®iÓm kiÓm tra, m« h×nh 647-645, 648-646 cã 3 ®iÓm kiÓm tra. Sai sè t¹i c¸c ®iÓm kiÓm tra phÇn lín lµ MX<0,1m, MY <0,35m, MZ <0,15m. T¹i m« h×nh 648-646 c¸ biÖt chØ cã ®iÓm kiÓm tra 2008 do n»m qu¸ s¸t mÐp ¶nh nªn cã sai sè lín ®Õn 1,1m. TÊt c¶ c¸c sai sè cßn l¹i cña c¸c ®iÓm tham gia x©y dùng m« h×nh ®Òu nhá h¬n 0.1m, ®a phÇn lµ rÊt nhá 0,01m. + §ît 20-02-2004: C¸c m« h×nh 823-819, 826-820 cã 6 ®iÓm kiÓm tra, m« h×nh 837-839, 838-840 cã 7 ®iÓm kiÓm tra. Sè l−îng ®iÓm kiÓm tra ®−îc t¨ng c−êng nhiÒu h¬n ®ît ®Çu. PhÇn lín c¸c ®iÓm kiÓm tra n»m trong ®é chÝnh x¸c MX<0,1m, MY <0,3m, MZ <0,25m. - Sù tham gia cña c¸c hÖ sè mÐo h×nh kÝnh vËt ®· lµm t¨ng ®é chÝnh x¸c lªn ®¸ng kÓ. Sè liÖu so s¸nh ®é chÝnh x¸c t¹i c¸c ®iÓm kiÓm tra kh«ng tham gia x©y dùng m« h×nh khi kh«ng cã vµ cã hÖ sè mÐo h×nh sù kh¸c biÖt rÊt lín. §a phÇn sai sè MP =4m-8m, c¸ biÖt cã ®iÓm lªn tíi 16m khi kh«ng cã hÖ sè mÐo h×nh. Khi cho hÖ sè mÐo h×nh th× c¸c sai sè nµy gi¶m ®i râ rÖt xuèng cßn 0,4m chiÕm ®a sè lµ nhá h¬n 0,3m. - KÕt qu¶ so s¸nh gi÷a ®iÓm ®o trªn ¶nh vµ ®iÓm ®o trùc tiÕp t¹i thùc ®Þa còng rÊt kh¶ quan. Täa ®é XY cña c¸c ®iÓm ®o ngo¹i nghiÖp b»ng m¸y toµn ®¹c ®−îc ®−a lªn m« h×nh lËp thÓ vµ so s¸nh ®é cao gi÷a Z®o ¶nh vµ Zt®. Víi sè l−îng ®iÓm ®o lµ 19, sù chªnh lÖch ®−îc thèng kª trong b¶ng 4.1 B¶ng 4.1. Thèng kª sù so s¸nh kÕt qu¶ ®o ngo¹i nghiÖp vµ kÕt qu¶ ®o ¶nh [ Sai sè MZ (m) 0,0 ÷ 0,1 0,1 ÷ 0,2 0,2 ÷ 0,3 0,3 ÷ 0,4 0,4 ÷ 0,5 Sè l−îng 1 3 5 9 1 NÕu h¹n sai ®é cao ®èi víi b¶n ®å ®Þa h×nh cã kho¶ng cao ®Òu lµ 1m lµ 0,33m th× sù so s¸nh gi÷a hai ph−¬ng ph¸p ®o ®éc lËp cã giíi h¹n sai sè lµ 86 m46,0233,0 =± th× c¸c kÕt qu¶ nµy hoµn toµn n»m trong ®é chÝnh x¸c cho phÐp. - KÕt qu¶ tÝnh to¸n khèi l−îng gi÷a hai DTM cña cïng mét khu vùc ®−îc x©y dùng trªn rÊt nhiÒu ®iÓm ®o cã ®é chÝnh x¸c tèt v× kÕt qu¶ cña DTM ®· ®−îc kiÓm tra ®é chÝnh x¸c chÆt chÏ ë c¸c phÇn x©y dùng m« h×nh vµ kiÓm tra ngo¹i nghiÖp. C«ng t¸c tÝnh to¸n ®¬n gi¶n vµ nhanh chãng. 4.3. thùc nghiÖm ®o vÏ phi ®Þa h×nh c«ng tr×nh kiÕn tróc C«ng t¸c thùc nghiÖm ®o vÏ c«ng tr×nh kiÕn tróc thuéc lÜnh vùc ®o ¶nh phi ®Þa h×nh ®−îc t¸c gi¶ tiÕn hµnh t¹i Nhµ h¸t Lín Hµ Néi. §©y lµ mét c«ng tr×nh kiÕn tróc næi tiÕng cña thñ ®« cã vÎ ®Ñp cæ kÝnh mang ®Ëm nÐt kiÕn tróc Ph¸p cña thÕ kØ 19. C«ng t¸c chôp ¶nh ®−îc t¸c gi¶ thiÕt kÕ c¸ch c«ng tr×nh kho¶ng 20m, ®−êng ®¸y chôp kho¶ng 4,5m, d¹ng chôp tù do vµ h¬i nghiªng lªn cao. TiÕn hµnh chôp 6 tÊm ¶nh b»ng m¸y chôp ¶nh sè Sony F707 hîp thµnh 3 cÆp ¶nh lËp thÓ ®éc lËp ®Ó kiÓm tra lÉn nhau, ®é ph©n gi¶i ®−îc lùa chän lµ 2560x1920pixel (5MP), chÕ ®é thu phãng v« cùc. C«ng t¸c ®o c¸c kho¶ng c¸ch kh«ng gian dïng ®Ó quy tû lÖ vµ kiÓm tra ®é chÝnh x¸c kÕt qu¶ ®o ®−îc tiÕn hµnh t¹i bËc thang, c¸c khung cöa, ®é réng vµ chiÒu cao cña c¸c cét b»ng th−íc 30m chÝnh x¸c cña NhËt B¶n. C«ng t¸c ngo¹i nghiÖp nµy lµ rÊt nhÑ nhµng vµ ®¬n gi¶n. C«ng t¸c ®o vÏ trong phßng tiÕn hµnh trªn phÇn mÒm Photomodeler Pro5. C«ng t¸c khai b¸o c¸c th«ng sè cña m¸y chôp ¶nh cã thÓ nhËp b»ng bµn phÝm tõ kÕt qu¶ kiÓm ®Þnh vµ còng cã thÓ gi¶i nhanh c¸c th«ng sè nµy d−íi sù trî gióp cña phÇn mÒm. ViÖc nµy ®−îc tiÕn hµnh ®¬n gi¶n b»ng c¸ch chôp mét tê giÊy cã kÝch cì x¸c ®Þnh (A4, A3, A2...) víi ®é ph©n gi¶i vµ chÕ ®é thu phãng gièng nh− lóc chôp ¶nh, sau ®ã ®Þnh vÞ mét « vu«ng cho khíp víi c¸c c¹nh cña h×nh ¶nh tê giÊy, ch−¬ng tr×nh sÏ tÝnh ngay ®−îc c¸c th«ng sè cña m¸y ¶nh (tiªu cù, täa ®é ®iÓm chÝnh ¶nh...). ViÖc ®Þnh h−íng cÆp ¶nh hoµn toµn kh«ng dïng ®iÓm khèng chÕ cã täa ®é mµ chØ quy tû lÖ b»ng kho¶ng c¸ch g¸n sè liÖu ®o mét c¹nh ®o th−íc vµo c¹nh ®o t−¬ng øng trªn m« h×nh. 87 §o vÏ trong hÖ täa ®é c«ng tr×nh lµ trôc X däc theo mÆt ngang c«ng tr×nh, trôc Z däc theo chiÒu ®øng c«ng tr×nh, trôc Y vu«ng gãc víi hai trôc XZ, ®iÓm gèc täa ®é lÊy tù do. C«ng t¸c quy tû lÖ vµ kiÓm tra ®−îc tr×nh bµy trong b¶ng 4.2, vµ nhË thÊy r»ng c¸c sù chªnh lÖch nµy lµ rÊt nhá. B¶ng 4.2. Quy tû lÖ vµ kiÓm tra m« h×nh c«ng tr×nh Nhµ h¸t Lín Hµ Néi STT S®o th−íc (m) S®o ¶nh (m) DS=S®o th−íc-S®o ¶nh Ghi chó 1 18,25 18,25 0,00 Quy tû lÖ 2 18,74 18,78 -0,05 BËc thang 3 20,91 20,92 -0,01 BËc thang 4 0,08 0,08 0,00 BiÓn b¸o 5 0,06 0,055 0,01 BiÓn b¸o 6 1,38 1,34 0,04 Cét 7 1,38 1,35 0,03 Cét C«ng t¸c ®o vÏ chi tiÕt c«ng tr×nh lµ sù ®o c¸c ®iÓm cïng tªn trªn hai tÊm ¶nh ë c¸c ®iÓm ®Æc tr−ng cña c«ng tr×nh nh− c¸c cét, c¸c cöa, c¸c ®−êng c¹nh t−êng, c¸c hoa v¨n häa tiÕt. V× trªn c«ng tr×nh cã rÊt nhiÒu yÕu tè gièng nhau nh− ®é réng chiÒu cao cña c¸c cöa sæ, c¸c cét, c¸c « kiÕn tróc nªn cã thÓ ®o mét « vµ sao chÐp cho c¸c vÞ trÝ cßn l¹i. NÕu cã ®o thªm th× sÏ lÊy mét sè liÖu trung b×nh cho h×nh d¸ng c¸c bé phËn gièng nhau nµy. Sau khi ®o vÏ xong th× xuÊt tÖp nµy sang ch−¬ng tr×nh AutoCad ®Ó biªn tËp chi tiÕt cho ®óng quy c¸ch cña b¶n vÏ kiÕn tróc. KÕt qu¶ vÒ ¶nh chôp vµ b¶n vÏ kiÕn tróc cña c«ng tr×nh Nhµ h¸t Lín Hµ Néi ®−îc tr×nh bµy trong phô lôc 19. LÜnh vùc ®o ¶nh phi ®Þa h×nh cßn cã thÓ tiÕn hµnh mét c¸ch ®¬n gi¶n khi ®o vÏ c¸c ®å vËt nh− bµn ghÕ, tñ g−êng, lä hoa... 88 kÕt luËn vμ kiÕn nghÞ Tæng hîp c¸c néi dung ®· tr×nh bµy trong luËn v¨n kÕt hîp qu¸ tr×nh thùc nghiÖm chøng minh, trªn c¬ së nh÷ng thµnh c«ng ®· ®¹t ®−îc, nh÷ng thuËn lîi, khã kh¨n, nh÷ng tån t¹i kh¸ch quan, t¸c gi¶ rót ra mét sè kÕt luËn vµ kiÕn nghÞ cña ®Ò tµi “C«ng nghÖ ®o ¶nh trªn c¬ së ¶nh chôp tõ m¸y chôp ¶nh sè phæ th«ng”. kÕt luËn 1. Sù ph¸t triÓn cña c«ng nghÖ ®o ¶nh mµ ®Ønh cao lµ c«ng nghÖ ¶nh sè cã ¶nh h−ëng rÊt lín ®Õn kh¶ n¨ng ®o vÏ hiÖn ®¹i nãi chung cña tr¾c ®Þa ¶nh. 2. LÜnh vùc øng dông cña ®o ¶nh sè ngµy cµng ®a d¹ng vµ t− liÖu ¶nh sè phæ th«ng lµ mét sù lùa chän hîp lý trong thêi ®iÓm hiÖn t¹i vµ t−¬ng lai nh»m ®o vÏ ®Þa h×nh tû lÖ lín, ®o ®¹c cËp nhËt vµ ®o vÏ phi ®Þa h×nh. 3. M¸y chôp ¶nh sè phæ th«ng rÎ tiÒn mµ tiÖn lîi khi chôp ¶nh. ViÖc kiÓm ®Þnh c¸c th«ng sè m¸y chôp ¶nh lo¹i nµy gåm tiªu cù, täa ®é ®iÓm chÝnh ¶nh vµ c¸c hÖ sè mÐo h×nh kÝnh vËt lµ c«ng viÖc tÊt yÕu. ViÖc nµy lµ biÖn ph¸p kh«ng thÓ thiÕu nh»m n©ng cao ®é chÝnh x¸c ®o vÏ ¶nh vµ tõ ®ã kh¼ng ®Þnh tÝnh kh¶ thi cña ph−¬ng ph¸p. 4. TÝnh kh¶ thi cña t− liÖu ¶nh sè phæ th«ng lµ rÊt râ rµng v× vËt t− thiÕt bÞ rÎ, dÔ mua, dÔ chôp vµ ®¹t ®é chÝnh x¸c cÇn thiÕt. ¶nh sè phæ th«ng gi¶i quyÕt ®−îc khã kh¨n vÒ th«ng tin ®Çu vµo cho ®o ¶nh vµ ®· thùc sù lo¹i bá m¸y quÐt ¶nh trong quy tr×nh c«ng nghÖ ®o vÏ ¶nh sè. 5. C«ng nghÖ ®o ¶nh sè sö dông ¶nh sè chôp tõ m¸y chôp ¶nh sè phæ th«ng ®· thµnh c«ng trong thµnh lËp b×nh ®å ®Þa h×nh tû lÖ lín, ®o ®¹c cËp nhËt vµ ®o vÏ phi ®Þa h×nh, tiÕt kiÖm thêi gian, n©ng cao ®é chÝnh x¸c. Kh¼ng ®Þnh ®−îc −u thÕ vÒ ®o vÏ vµ tÝnh to¸n cña ph−¬ng ph¸p ®o ¶nh sè so víi c¸c ph−¬ng ph¸p ®o vÏ truyÒn thèng kh¸c. 6. Sù ph¸t triÓn nhanh chãng cña c«ng nghÖ s¶n xuÊt m¸y chôp ¶nh sè vÒ sù n©ng cao ®é ph©n gi¶i, n©ng cao chÊt l−îng hÖ thèng kÝnh vËt vµ sù c¹nh tranh vÒ gi¸ c¶ sÏ lµ c¬ së vËt chÊt kü thuËt kh¼ng ®Þnh tÝnh hîp lý cña h−íng 89 nghiªn cøu. Cïng víi c¸c c¬ së to¸n häc chÆt chÏ, ch−¬ng tr×nh xö lý sè liÖu kiÓm ®Þnh hîp lý khiÕn cho tÝnh kh¶ thi cao vÒ kinh tÕ kü thuËt cña ®o vÏ ¶nh sè phæ th«ng cµng ®−îc thuyÕt phôc râ rÖt. 7. Môc ®Ých vµ ph¹m vi nghiªn cøu cña ®Ò tµi ®Æt ra ban ®Çu ®· ®−îc gi¶i quyÕt trän vÑn. KiÕn nghÞ 1. Trªn c¬ së c¸c thµnh qu¶ ®· ®¹t ®−îc vµ kÕt luËn cña ®Ò tµi, v× tÝnh kü thuËt ®−îc ®¶m b¶o, tÝnh −u thÕ cña ®o vÏ tÝnh to¸n, tÝnh kinh tÕ khi ®Çu t− tr¹m ¶nh sè vµ m¸y chôp ¶nh t¹i thêi ®iÓm hiÖn t¹i, nªn øng dông réng r·i c«ng t¸c ®o vÏ sö dông ¶nh sè phæ th«ng trong s¶n xuÊt thùc tiÔn, nhÊt lµ c«ng t¸c ®o vÏ b×nh ®å tÝnh to¸n khèi l−îng khai th¸c cËp nhËt cña má than lé thiªn. Phèi hîp víi c¸c c¬ quan chøc n¨ng kh¸c nh− ngµnh c«ng an, ngµnh v¨n hãa, b¶o tån b¶o tµng, kiÕn tróc x©y dùng, n«ng l©m nghiÖp, y tÕ... nghiªn cøu ®Çu t− triÓn khai øng dông ®o ¶nh sè phæ th«ng nh»m ®o vÏ chÝnh x¸c, kh¸ch quan hiÖn tr−êng tai n¹n, phôc chÕ kiÕn tróc v¨n hãa, tÝnh to¸n sù ph¸t triÓn gia sóc, c©y trång... 2. CÇn ph¶i x©y dùng mét b·i kiÓm ®Þnh chuÈn víi sè l−îng ®iÓm lín, tr¶i ®Òu vµ l−u tr÷ ®−îc l©u dµi ®Ó cã ®−îc kÕt qu¶ kiÓm ®Þnh kh¸ch quan h¬n vµ ®Þnh k× trong qu¸ tr×nh sö dông m¸y chôp ¶nh sè phæ th«ng. 3. CÇn nghiªn cøu vµ tiÕp tôc kh¶o s¸t c¸c ®iÒu kiÖn ®Ó ®¶m b¶o kh¶ n¨ng t¨ng kho¶ng c¸ch ®Õn ®èi t−îng ®o vÏ h¬n n÷a nhÊt lµ trong lÜnh vùc ®Þa h×nh. 4. §Ó cã h−íng ph¸t triÓn ®a d¹ng vµ hoµn thiÖn trong ®o ¶nh, cÇn ph¶i øng dông ®−îc ph−¬ng ph¸p t¨ng dµy tam gi¸c ¶nh kh«ng gian ®èi víi ¶nh sè phæ th«ng nh»m gi¶m thiÓu sè l−îng ®iÓm ®o ngo¹i nghiÖp t¹i thùc ®Þa. §©y lµ mét sè kiÕn nghÞ c¬ b¶n trªn c¬ së tÝnh x©y dùng cho sù øng dông thùc tiÔn, t¸c gi¶ cã nguyÖn väng r»ng sÏ ®−îc tiÕp tôc nghiªn cøu theo h−íng nµy ®Ó lµm cho tÝnh ®a d¹ng cña ®o ¶nh ®−îc ph¸t triÓn thùc chÊt trong ®iÒu kiÖn ViÖt Nam. 90 danh môc c«ng tr×nh cña t¸c gi¶ 1. Phan V¨n Léc, TrÇn Trung Anh vµ nnk (2004), øng dông chôp ¶nh sè phæ th«ng vµo viÖc ®o vÏ b¶n ®å ®Þa h×nh vµ x¸c ®Þnh khèi l−îng khai th¸c ë má lé thiªn, §Ò tµi hîp t¸c khoa häc ViÖt Nam-Ba Lan, m· sè HTNC-02, Hµ Néi. 2. Phan V¨n Léc, TrÇn Trung Anh (2004), “C¸c c¬ së lý thuyÕt ®¶m b¶o viÖc ¸p dông ¶nh sè cù ly gÇn vµo c«ng t¸c ®o vÏ ë má lé thiªn”, T¹p chÝ C«ng NghiÖp Má sè 6-2004, Héi Khoa häc vµ C«ng nghÖ Má ViÖt Nam. 3. TrÇn Trung Anh, Phan V¨n Léc (2004), “KiÓm ®Þnh m¸y chôp ¶nh sè phæ th«ng vµ kh¶ n¨ng øng dông trong thµnh lËp b×nh ®å tû lÖ lín”, T¹p chÝ Khoa häc Kü thuËt Má §Þa chÊt sè (8), tr 33-37. 4. TrÇn Trung Anh (2004), Nghiªn cøu mét sè ®iÒu kiÖn c¬ b¶n ®¶m b¶o ®o chôp ¶nh sè phæ th«ng cù ly gÇn, §Ò tµi Nghiªn cøu Khoa häc cÊp tr−êng, §¹i häc Má §Þa chÊt, Hµ Néi. 91 Tμi liÖu tham kh¶o 1. Lª V¨n H−êng (2000), §o ¶nh mÆt ®Êt, Nxb Giao Th«ng VËn T¶i, Hµ Néi. 2. Tr−¬ng Anh KiÖt (2000), C¬ së ®o ¶nh, Nxb Giao Th«ng VËn T¶i, Hµ Néi. 3. Tr−¬ng Anh KiÖt (2000), Ph−¬ng ph¸p ®o ¶nh gi¶i tÝch vµ ®o ¶nh sè, Bµi gi¶ng cao häc ngµnh Tr¾c ®Þa, Tr−êng §¹i häc Má §Þa chÊt, Hµ Néi. 4. Phan V¨n Léc (2000), §o ¶nh lËp thÓ, Nxb Giao Th«ng VËn T¶i, Hµ Néi. 5. Phan V¨n Léc (2000), Tù ®éng hãa ®o ¶nh, Bµi gi¶ng cao häc ngµnh Tr¾c §Þa, Tr−êng §¹i häc Má §Þa chÊt, Hµ Néi. 6. Ph¹m Väng Thµnh (2000), “Nghiªn cøu kh¶ n¨ng øng dông m¸y ¶nh kü thuËt sè trong c«ng t¸c kiÕn tróc x©y dùng”, TuyÓn tËp c¸c c«ng tr×nh Khoa häc sè(11-2000), Tr−êng §¹i häc Má §Þa chÊt, Hµ Néi. 7. TrÇn §×nh TrÝ (2004), X©y dùng thuËt to¸n vµ qui tr×nh c«ng nghÖ ®o vÏ gi¶i tÝch ¶nh ®−îc chôp tõ kho¶ng c¸ch gÇn b»ng m¸y chôp ¶nh kh«ng chuyªn trong lÜnh vùc phi ®Þa h×nh, B¸o c¸o §Ò tµi Nghiªn cøu Khoa häc cÊp Bé B2002-36-31, Hµ Néi. 8. Tiªu chuÈn ngµnh 96 TCN 43-90 (1990), Quy ph¹m ®o vÏ b¶n ®å ®Þa h×nh tû lÖ 1:500, 1:1000, 1:2000, 1:5000 (phÇn ngoµi trêi), Côc §o ®¹c vµ B¶n ®å Nhµ n−íc. 9. A.H. Lobanov (1984), §o ¶nh, Nxb Nedra, M¸txc¬va. 10. Surveing Engineering Department (2005), History of Photogrammetry, Ferris State University, USA. 11. 3DMAnalyst (2005), Mining Analysis Software, ADAM technology, Australia. 12. EosSystems Inc (2003), Photomodeler Pro 5 User’s guide, Canada. 13. Photometrix Pty Ltd (2004), Report on Accuracy Test of iWitness 3D Measurement System, Australia. 14. Christopher Dean (2003 ), “Calculation of Wood Volume and Stem Taper Using Terrestrial Single Image Close Range Photogrammetry and Contemporary Software Tools”, Silva Fennica 37(3), pp 359–380. 15. Vassilios T (2000), “An Automatic process for the extraction of the 3D model of the human back surface for scoliosis treatment”, IAPRS, Vol. XXXIII, Amsterdam. 92 16. Heinrich S, Jurgen H, Lothar G (2004), LIMEZ-Photogrammetric Measurement of Railroad Clearance Obstacles, Germany. 17. Dimitar Jechev (2004), Close Range Photogrammetry with amateur camera, GIS Sofia Ltd., Bulgaria. 18. Rainer S, Bernhard B, Hans D, Andreas E, Stefan H, Joe H, Werner K, Erik L, Ralf R, Stefan W (2000), “Design principle of the LH systems ADS40 airborne digital sensor”, IAPRS, Vol. XXXIII, Amsterdam. 19. Jeff Lower, Jeremy Conner (2004), “Photogrammetric Mapping With the Z/I Digital Mapping Camera (DMC)”, GeoSpatial World. 20. Eric Meisenzahl (2004), Charge-Coupled Device Image Sensors, Kodak, USA. 21. Jim Janesick (2004), Scientific CMOS pixels, Sarnoff Corporation CN5300, 201 Washington Road Princeton, NJ 08543-5300. 22. Dipl.-Ing. Robert Godding (2004), Geometric Calibration and Orientation of Digital Imaging Systems, AICON 3D Systems GmbH, Braunschweig. 23. Ayman F.H, Michel M, Young R.L (2002), “Bundle adjustment with self- calibration using straight lines”, Photogrammetric Record 17(100), pp 635-650. 24. Moumen T, Aly A (2003), “Nonmetric Lens Distortion Calibration: Closed- form Solutions, Robust Estimation and Model Selection”, Proceedings of the Ninth IEEE International Conference on Computer Vision Volume 2, Computer Society. 25. Michael Cramer (2004), “Digital Camera Calibration”, XXth ISPRS congress, EuroSDR network, Istalbul, Turkey. 26. Z/I imaging (2001), ISPM & ISDM User’s guide, Zeiss&Integraph, USA. 27. Naser El Sheimy (1999), Adjustment of Observations, Department of Geomatics Engineering, The University of Calragy. 28. Digital-Photogrammetry.org (2005), 20 Megapixel Fujifilm CCD sensor and 22 Megapixel Kodak-Sinar CCD sensor, FLAARnetwork. 29. Canon Consumer Products (2004), The highest resolution sensor ever to be used in an EOS 16,7 Megapixel, Full Frame CMOS Sensor with DIGIC II Image Processor, USA. phô lôc Phô lôc 1 KÕt qu¶ ®Þnh h−íng vµ kiÓm tra m« h×nh: ¶nh chôp 11-10-2003 khu 1 Photo-T Parameters and Results for Project 637-643 (KHU 1) Object Statistics Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 1001 Control -0.010 -0.001 0.015 0.018 429406.158 26238.660 -37.110 429406.168 26238.661 -37.125 1002 Check -0.046 -0.146 0.037 0.158 429452.445 26276.071 -41.331 429452.491 26276.217 -41.368 1003 Control 0.005 -0.003 -0.012 0.013 429512.877 26360.249 -48.750 429512.872 26360.252 -48.738 1004 Control 0.009 0.003 -0.009 0.013 429378.759 26308.209 -25.351 429378.750 26308.206 -25.342 1005 Check -0.087 -0.062 0.002 0.107 429427.314 26372.302 -25.828 429427.401 26372.364 -25.830 1006 Check -0.040 -0.075 -0.012 0.086 429479.728 26427.614 -27.486 429479.768 26427.689 -27.474 1007 Control -0.001 0.000 0.003 0.003 429549.447 26510.799 -23.542 429549.448 26510.799 -23.545 1008 Control -0.003 0.000 0.003 0.004 429544.311 26670.974 32.713 429544.314 26670.974 32.710 1009 Check -0.023 0.396 0.022 0.397 429436.760 26532.036 17.938 429436.783 26531.640 17.916 Exterior Orientation Strip Id Photo Id X Y Z Omega Phi Kappa Status 1 643 429467.933 26126.443 -5.626 1.518 0.092 -0.026 Used 1 637 429431.729 26130.934 -6.112 1.561 0.036 0.006 Used Photo-T Parameters and Results for Project 638-644 (KHU 1) Object Statistics Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 1001 Control -0.010 0.001 0.015 0.018 429406.158 26238.662 -37.110 429406.168 26238.661 -37.125 1002 Check -0.031 -0.151 0.028 0.156 429452.460 26276.066 -41.340 429452.491 26276.217 -41.368 1003 Control 0.005 -0.003 -0.012 0.014 429512.877 26360.249 -48.750 429512.872 26360.252 -48.738 1004 Control 0.009 0.002 -0.009 0.013 429378.759 26308.208 -25.351 429378.750 26308.206 -25.342 1005 Check -0.089 -0.044 0.006 0.099 429427.312 26372.320 -25.824 429427.401 26372.364 -25.830 1006 Check -0.030 -0.064 -0.006 0.071 429479.738 26427.625 -27.480 429479.768 26427.689 -27.474 1007 Control -0.001 0.000 0.003 0.004 429549.447 26510.799 -23.542 429549.448 26510.799 -23.545 1008 Control -0.003 0.000 0.003 0.004 429544.311 26670.974 32.713 429544.314 26670.974 32.710 1009 Check -0.031 0.377 0.031 0.379 429436.752 26532.017 17.947 429436.783 26531.640 17.916 Exterior Orientation Strip Id Photo Id X Y Z Omega Phi Kappa Status 1 644 429467.946 26126.473 -5.637 1.528 0.125 -0.022 Used 1 638 429431.746 26130.888 -6.084 1.561 0.036 0.006 Used Phô lôc 2 KÕt qu¶ ®Þnh h−íng vµ kiÓm tra m« h×nh: ¶nh chôp 11-10-2003 khu 2 Photo-T Parameters and Results for Project 647-645 (KHU 2) Object Statistics Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 2001 Control -0.010 -0.001 0.004 0.011 430421.110 26536.449 -28.726 430421.120 26536.450 -28.730 2002 Control 0.013 -0.003 0.003 0.014 430515.386 26508.379 -20.978 430515.373 26508.382 -20.981 2004 Control 0.002 0.000 -0.009 0.009 430372.423 26653.879 -7.465 430372.421 26653.879 -7.456 2005 Check 0.007 0.084 0.045 0.096 430450.538 26657.019 -7.704 430450.531 26656.935 -7.749 2006 Check 0.040 -0.250 0.031 0.255 430548.519 26648.380 -7.645 430548.479 26648.630 -7.676 2007 Control -0.011 0.004 -0.005 0.013 430589.603 26603.067 -9.730 430589.614 26603.063 -9.725 2008 Check -0.188 -0.603 -0.076 0.636 430670.104 26714.711 33.618 430670.292 26715.314 33.694 2009 Control 0.005 -0.001 0.007 0.009 430477.287 26764.158 32.440 430477.282 26764.159 32.433 Exterior Orientation Strip Id Photo Id X Y Z Omega Phi Kappa Status 2 647 430423.622 26234.720 0.227 1.599 -0.056 -0.005 Used 2 645 430493.030 26224.566 2.886 1.522 0.037 0.029 Used Photo-T Parameters and Results for Project 648-646 (KHU 2) Object Statistics Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 2001 Control -0.001 0.000 0.004 0.004 430421.119 26536.450 -28.726 430421.120 26536.450 -28.730 2002 Control -0.002 0.001 -0.003 0.004 430515.371 26508.383 -20.984 430515.373 26508.382 -20.981 2004 Control 0.002 0.001 -0.002 0.003 430372.423 26653.880 -7.458 430372.421 26653.879 -7.456 2005 Check 0.014 0.333 0.033 0.335 430450.545 26657.268 -7.716 430450.531 26656.935 -7.749 2006 Control 0.004 -0.001 0.001 0.005 430548.483 26648.629 -7.675 430548.479 26648.630 -7.676 2008 Check -0.525 -1.100 -0.112 1.224 430669.767 26714.214 33.582 430670.292 26715.314 33.694 2009 Control -0.004 0.000 0.000 0.004 430477.278 26764.159 32.433 430477.282 26764.159 32.433 Exterior Orientation Strip Id Photo Id X Y Z Omega Phi Kappa Status 2 646 430493.030 26224.676 2.817 1.522 0.051 0.027 Used 2 648 430423.592 26234.840 0.246 1.593 -0.135 -0.004 Used Phô lôc 3 KÕt qu¶ ®Þnh h−íng vµ kiÓm tra m« h×nh: ¶nh chôp 20-02-2004 khu 1 Photo-T Parameters and Results for Project 823-819 (KHU 1) Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 1002 Control -0.003 -0.002 0.001 0.004 429393.407 26326.788 -29.599 429393.410 26326.790 -29.600 1003 Check -0.003 -0.140 0.054 0.150 429453.569 26376.259 -30.918 429453.572 26376.399 -30.972 1004 Check -0.026 -0.213 0.050 0.221 429482.466 26416.869 -33.159 429482.492 26417.082 -33.209 1005 Control -0.007 0.001 0.000 0.008 429503.930 26437.846 -36.408 429503.937 26437.845 -36.408 1006 Check 0.032 0.034 0.024 0.052 429422.787 26352.685 -30.494 429422.755 26352.651 -30.518 1007 Control 0.007 -0.002 -0.001 0.007 429548.248 26466.742 -40.646 429548.241 26466.744 -40.645 1008 Control 0.004 0.004 -0.001 0.006 429420.273 26249.177 -55.923 429420.269 26249.173 -55.922 1009 Check 0.011 -0.102 0.085 0.133 429439.296 26266.427 -54.774 429439.285 26266.529 -54.859 1010 Control 0.000 -0.142 0.082 0.164 429466.719 26295.215 -55.957 429466.719 26295.357 -56.039 1011 Check -0.001 -0.236 0.144 0.276 429470.084 26338.295 -52.883 429470.085 26338.531 -53.027 1012 Check 0.068 0.006 0.054 0.087 429494.254 26371.277 -49.803 429494.186 26371.271 -49.857 Strip Id Photo Id X Y Z Omega Phi Kappa Status 1 823 429431.799 26130.705 -5.879 1.244 0.020 -0.035 Used 1 819 429465.904 26140.087 -9.370 1.418 0.069 -0.034 Used Photo-T Parameters and Results for Project 826-820 (KHU 1) Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 1002 Control 0.001 0.000 0.004 0.004 429393.411 26326.790 -29.596 429393.410 26326.790 -29.600 1003 Check 0.009 0.068 0.033 0.077 429453.581 26376.467 -30.939 429453.572 26376.399 -30.972 1004 Check 0.033 0.290 0.012 0.292 429482.525 26417.372 -33.197 429482.492 26417.082 -33.209 1005 Control -0.006 0.000 -0.010 0.012 429503.931 26437.845 -36.418 429503.937 26437.845 -36.408 1006 Check 0.027 0.142 0.065 0.159 429422.782 26352.793 -30.453 429422.755 26352.651 -30.518 1007 Control 0.004 0.000 0.008 0.009 429548.245 26466.744 -40.637 429548.241 26466.744 -40.645 1008 Control -0.001 0.000 -0.001 0.002 429420.268 26249.173 -55.923 429420.269 26249.173 -55.922 1009 Check 0.062 0.083 0.006 0.104 429439.347 26266.612 -54.853 429439.285 26266.529 -54.859 1010 Control 0.092 0.111 0.002 0.144 429466.811 26295.468 -56.037 429466.719 26295.357 -56.039 1011 Check 0.053 0.004 0.108 0.121 429470.138 26338.535 -52.919 429470.085 26338.531 -53.027 1012 Check 0.062 0.160 0.032 0.175 429494.248 26371.431 -49.825 429494.186 26371.271 -49.857 Strip Id Photo Id X Y Z Omega Phi Kappa Status 1 826 429431.921 26130.588 -5.790 1.347 -0.114 -0.094 Used 1 820 429466.113 26139.941 -9.316 1.417 0.076 -0.034 Used Phô lôc 4 KÕt qu¶ ®Þnh h−íng vµ kiÓm tra m« h×nh: ¶nh chôp 20-02-2004 khu 2 Photo-T Parameters and Results for Project 837-839 (KHU 2) Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 2001 Check 0.043 0.046 -0.268 0.276 430461.165 26658.101 -8.378 430461.122 26658.055 -8.110 2002 Control 0.000 -0.002 -0.055 0.055 430504.443 26661.901 -8.285 430504.443 26661.903 -8.230 2003 Check 0.012 0.094 -0.146 0.174 430546.805 26657.095 -6.823 430546.793 26657.001 -6.677 2005 Check 0.080 0.144 0.085 0.185 430609.464 26590.033 -12.524 430609.384 26589.889 -12.609 2006 Check 0.049 0.216 0.176 0.283 430635.912 26552.055 -14.007 430635.863 26551.839 -14.183 2007 Control -0.006 0.008 0.077 0.078 430646.801 26508.682 -14.333 430646.807 26508.674 -14.410 2012 Control 0.002 -0.019 -0.129 0.130 430594.145 26400.403 -21.228 430594.143 26400.422 -21.099 2104 Check -0.027 -0.040 0.009 0.049 430592.876 26627.385 -10.004 430592.903 26627.425 -10.013 2108 Control 0.005 0.012 0.106 0.107 430509.193 26456.347 -23.987 430509.188 26456.335 -24.093 2109 Check 0.054 0.170 0.224 0.286 430523.991 26450.770 -22.602 430523.937 26450.600 -22.826 2110 Check 0.130 0.239 0.235 0.359 430551.073 26438.124 -21.861 430550.943 26437.885 -22.096 Strip Id Photo Id X Y Z Omega Phi Kappa Status 2 839 430493.011 26224.530 3.859 1.474 -0.158 0.002 Used 2 837 430423.627 26234.832 0.647 1.534 -0.333 -0.022 Used Photo-T Parameters and Results for Project 838-840 (KHU 2) Point Id Type VX VY VZ V(XYZ) ComputX ComputY ComputZ Given X Given Y Given Z 2001 Check 0.123 0.709 -0.320 0.788 430461.245 26658.764 -8.430 430461.122 26658.055 -8.110 2002 Control -0.007 0.000 -0.052 0.053 430504.436 26661.903 -8.282 430504.443 26661.903 -8.230 2003 Check -0.031 0.091 -0.054 0.110 430546.762 26657.092 -6.731 430546.793 26657.001 -6.677 2005 Check 0.072 0.189 0.242 0.315 430609.456 26590.078 -12.367 430609.384 26589.889 -12.609 2006 Check -0.187 -0.047 0.369 0.416 430635.676 26551.792 -13.814 430635.863 26551.839 -14.183 2007 Control 0.005 0.000 0.073 0.073 430646.812 26508.674 -14.337 430646.807 26508.674 -14.410 2012 Control -0.008 -0.008 -0.123 0.123 430594.135 26400.414 -21.222 430594.143 26400.422 -21.099 2104 Check -0.100 -0.139 0.122 0.210 430592.803 26627.286 -9.891 430592.903 26627.425 -10.013 2108 Control 0.010 0.008 0.102 0.102 430509.198 26456.343 -23.991 430509.188 26456.335 -24.093 2109 Check 0.039 0.069 0.222 0.236 430523.976 26450.669 -22.604 430523.937 26450.600 -22.826 2110 Check 0.087 0.113 0.291 0.324 430551.030 26437.998 -21.805 430550.943 26437.885 -22.096 Strip Id Photo Id X Y Z Omega Phi Kappa Status 2 840 430492.983 26224.513 3.624 1.481 -0.348 -0.016 Used 2 838 430423.616 26234.637 0.658 1.528 -0.465 -0.022 Used Phô lôc 5 Täa ®é mét sè ®iÓm chi tiÕt ®o trªn m« h×nh vµ ®o trùc tiÕp ngoµi thùc ®Þa STT X (M) Y(M) Z(M) ®o ¶nh Z(M) ®o trùc tiÕp DZ 1 429396.881 26328.413 -31.10 -30.84 -0.26 2 429408.617 26290.336 -54.19 -54.53 0.34 3 429425.094 26310.630 -54.17 -54.55 0.38 4 429442.223 26303.348 -54.80 -54.93 0.13 5 429441.912 26270.887 -55.26 -55.66 0.4 6 429464.217 26291.637 -56.27 -56.41 0.14 7 429473.129 26331.171 -58.33 -58.58 0.25 8 429466.360 26333.089 -53.56 -53.94 0.38 9 429515.356 26353.528 -60.62 -60.97 0.35 10 429500.369 26380.393 -49.71 -50.08 0.37 11 429477.449 26369.914 -50.73 -51.18 0.45 12 429447.667 26367.474 -32.39 -32.69 0.3 13 429460.643 26393.150 -33.78 -33.63 -0.15 14 429495.301 26407.874 -47.86 -48.25 0.39 15 429494.828 26430.820 -36.38 -36.59 0.21 16 429511.883 26444.674 -37.85 -38.13 0.28 17 429525.234 26447.936 -43.61 -43.94 0.33 18 429528.111 26461.262 -39.88 -39.89 0.01 19 429442.628 26330.964 -53.91 -54.26 0.35 Phô lôc 6 KÕt qu¶ tÝnh to¸n khèi l−îng khai th¸c t¹i má lé thiªn Cäc S¸u (tõ ngµy 11-10-2003 ®Õn ngµy 20-02-2004) TÝnh khèi l−îng gi÷a hai DTM b»ng phÇn mÒm MGE Terrain analys Giíi h¹n tÝnh khu §«ng Nam (khu 1): #1 429408.438 26327.745 #2 429479.149 26398.456 #3 429514.504 26363.101 #4 429443.793 26292.390 Area: 5000 m2 DTM 637-643 (9h30 ngay 11-10-2003) DTM 826-840 (10h05 ngay 20-02-2004) Volume Result by MGE Terrain analys V = 46067 m3 Giíi h¹n tÝnh khu §«ng B¾c (khu 2): #1 430424.050 26624.250 #2 430424.050 26674.250 #3 430574.050 26624.250 #4 429443.793 26674.250 Area: 7500 m2 DTM 648-646 (11h50 ngay 11-10-2003) DTM 837-839 (11h57 ngay 20-02-2004) Volume Result by MGE Terrain analys V = 5254 m3 Phô lôc 7 M· nguån mét sè hµm vµ thñ tôc gi¶i bµi to¸n kiÓm ®Þnh ¶nh unit Calculation_Calibration; (*APP: CALIBRATION CAMERA FOR NON-METRIC DIGITAL CAMERA AND AERIAL CAMERA PROGRAMER : TRAN TRUNG ANH OFF : DEPARMENT OF PHOTOGRAMMETRY. HA NOI UNIVERSITY OF MINING&GEOLOGY ADD : DONG NGAC - TU LIEM - HA NOI TEL : (84)4.8387987 - (84)912090304 EMAIL: TRUNGANH7885@YAHOO.COM PROGRAM LANGUAGE : BORLAND DELPHI. ENTERPRISE VERSION 7.0 *) interface uses Classes,SysUtils,Dialogs ; type TMyThread = class(TThread) private Procedure Calculation; protected procedure Execute; override; // Main thread execution published constructor CreateIt; destructor Destroy; override; end; type TIC=Record PN:String; x,y:Real; PE:String; end; TCtr=Record PN:String; X,Y,Z:Real; PE:String; end; Const KTra='KT';Loai='Loai';Nghin=1000;Nghin2=1000000;Nghin4=1000000000000; var IC:Array of TIC; KC:Array of TCtr; AE,ER:Array[0..17]of Real; AX:Array[0..10]of Real; NE:Array[0..160]of Real; implementation uses windows,Camera,ExCamera, MAIN, Hansai, KQtinh, VCF1; Procedure RotateMatrix(Matdat:Boolean;phi,omega,kapa:Real;Var AX:Array of Real); Begin if not Matdat then begin AX[1]:=Cos(phi)*Cos(kapa)-Sin(phi)*Sin(omega)*Sin(kapa); AX[2]:=-Cos(phi)*Sin(kapa)-Sin(phi)*Sin(omega)*Cos(kapa);; AX[3]:=-Sin(phi)*Cos(omega); AX[4]:=Cos(omega)*Sin(kapa); AX[5]:=Cos(omega)*Cos(kapa); AX[6]:=-Sin(omega); AX[7]:=Sin(phi)*Cos(kapa)+Cos(phi)*Sin(omega)*Sin(kapa); AX[8]:=-Sin(phi)*Sin(kapa)+Cos(phi)*Sin(omega)*Cos(kapa);; AX[9]:=Cos(phi)*Cos(omega); end else begin //According to LOBANOV for Close range AX[1]:=Cos(phi)*Cos(kapa)-Sin(phi)*Sin(omega)*Sin(kapa); AX[2]:=Sin(phi)*Cos(omega); AX[3]:=-Cos(phi)*Sin(kapa)-Sin(phi)*Sin(omega)*Cos(kapa); AX[4]:=-Sin(phi)*Cos(kapa)-Cos(phi)*Sin(omega)*Sin(kapa); AX[5]:=Cos(phi)*Cos(omega); AX[6]:=Sin(phi)*Sin(kapa)-Cos(phi)*Sin(omega)*Cos(kapa); AX[7]:=Cos(omega)*Sin(kapa); AX[8]:=Sin(omega); AX[9]:=Cos(omega)*Cos(kapa); end; End; Procedure CorrectDistortion(xa,ya,x0,y0,K1,K2,K3,P1,P2:Real;Var xn,yn:Real); Var rd,dr,dx,dy:Real; Begin //According to Brown (Intergraph) rd:=Sqrt(Sqr(xa-x0)+Sqr(ya-y0)); dr:=K1+K2*Sqr(rd)+K3*Sqr(Sqr(rd)); dx:=P1*(Sqr(rd)+2*Sqr(xa-x0))+2*P2*(xa-x0)*(ya-y0); dy:=2*P1*(xa-x0)*(ya-y0)+P2*(Sqr(rd)+2*Sqr(ya-y0)); xn:=xa-(xa-x0)*dr+dx; yn:=ya-(ya-y0)*dr+dy; End; //Image coordinate point on photogrammetric spatial coordinate system Procedure XYZtoXNYNZN(Matdat:Boolean;X,Y,Z,f,x0,y0,XS,YS,ZS,hk:Real;AX:Array of Real;Var XN,YN,ZN,xt,yt:Real); begin XN:=AX[1]*(X-XS)+AX[4]*(Y-YS)+AX[7]*(Z-ZS); YN:=AX[2]*(X-XS)+AX[5]*(Y-YS)+AX[8]*(Z-ZS); ZN:=AX[3]*(X-XS)+AX[6]*(Y-YS)+AX[9]*(Z-ZS); if not Matdat then begin xt:=-f*XN/ZN+x0; yt:=-hk*f*YN/ZN+y0; end else begin xt:=f*XN/YN+x0; yt:=hk*f*ZN/YN+y0; end; end; //Setup image coordinate adjustment equations: one point has two equations Procedure SetupAEx(Matdat:Boolean;xa,ya,X,Y,Z,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk:Re al;Var AE:Array of Real); var XN,YN,ZN,xt,yt:Real; begin RotateMatrix(Matdat,phi,omega,kapa,AX); XYZtoXNYNZN(Matdat,X,Y,Z,f,x0,y0,XS,YS,ZS,hk,AX,XN,YN,ZN,xt,yt); if not Matdat then begin AE[1]:=(AX[3]*(xa-x0)+AX[1]*f)/ZN;// dx/dXS AE[2]:=(AX[6]*(xa-x0)+AX[4]*f)/ZN;// dx/dYS AE[3]:=(AX[9]*(xa-x0)+AX[7]*f)/ZN;// dx/dZS AE[4]:=(AX[7]*(X-XS)-AX[1]*(Z-ZS))*f/ZN+(AX[9]*(X-XS)-AX[3]*(Z- ZS))*(xa-x0)/ZN;// dx/dphi AE[5]:=(AX[6]*(X-XS)*Sin(phi)+(Y-YS)*Cos(omega)-AX[6]*(Z- ZS)*Cos(phi))*(xa-x0)/ZN-f*Sin(kapa); // dx/domega end else begin AE[1]:=(AX[2]*(xa-x0)-AX[1]*f)/YN;// dx/dXS AE[2]:=(AX[5]*(xa-x0)-AX[4]*f)/YN;// dx/dYS AE[3]:=(AX[8]*(xa-x0)-AX[7]*f)/YN;// dx/dZS AE[4]:=(AX[4]*(X-XS)-AX[1]*(Y-YS))*f/YN-(AX[5]*(X-XS)-AX[2]*(Y- YS))*(xa-x0)/YN;// dx/dphi AE[5]:=(AX[8]*(X-XS)*Sin(phi)+AX[8]*(Y-YS)*Cos(phi)-(Z- ZS)*Cos(omega))*(xa-x0)/YN-f*Sin(kapa); // dx/domega end; AE[6]:=(ya-y0)/hk;// dx/dkapa AE[7]:=(xa-x0)/f;// dx/df AE[8]:=1;// dx/dx0 AE[9]:=0;// dx/dz0 AE[10]:=0;// dx/dhk AE[11]:=xt-xa;// Lx end; Procedure SetupAEy(Matdat:Boolean;xa,ya,X,Y,Z,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk:Re al;Var AE:Array of Real); var XN,YN,ZN,xt,yt:Real; begin RotateMatrix(Matdat,phi,omega,kapa,AX); XYZtoXNYNZN(Matdat,X,Y,Z,f,x0,y0,XS,YS,ZS,hk,AX,XN,YN,ZN,xt,yt); if not Matdat then begin AE[1]:=(AX[3]*(ya-y0)+hk*AX[2]*f)/ZN;// dy/dXS AE[2]:=(AX[6]*(ya-y0)+hk*AX[5]*f)/ZN;// dy/dYS AE[3]:=(AX[9]*(ya-y0)+hk*AX[6]*f)/ZN;// dy/dZS AE[4]:=(AX[8]*(X-XS)-AX[2]*(Z-ZS))*hk*f/ZN+(AX[9]*(X-XS)-AX[3]*(Z- ZS))*(ya-y0)/ZN;// dy/dphi AE[5]:=(AX[6]*(X-XS)*Sin(phi)+(Y-YS)*Cos(omega)-AX[6]*(Z- ZS)*Cos(phi))*(ya-y0)/ZN-hk*f*Cos(kapa); // dy/domega end else begin AE[1]:=(AX[2]*(ya-y0)-hk*AX[3]*f)/YN;// dz/dXS AE[2]:=(AX[5]*(ya-y0)-hk*AX[6]*f)/YN;// dz/dYS AE[3]:=(AX[8]*(ya-y0)-hk*AX[9]*f)/YN;// dz/dZS AE[4]:=(AX[6]*(X-XS)-AX[3]*(Y-YS))*hk*f/YN-(AX[5]*(X-XS)-AX[2]*(Y- YS))*(ya-y0)/YN;// dz/dphi AE[5]:=(AX[8]*(X-XS)*Sin(phi)+AX[8]*(Y-YS)*Cos(phi)-(Z- ZS)*Cos(omega))*(ya-y0)/YN-hk*f*Cos(kapa); // dz/domega end; AE[6]:=-hk*(xa-x0);// dy/dkapa AE[7]:=(ya-y0)/f;// dy/df AE[8]:=0;// dy/dx0 AE[9]:=1;// dy/dz0 AE[10]:=(ya-y0)/hk;// dx/dhk AE[11]:=yt-ya;// Lz end; // Setup adjustment equations radial and tangential distortions Procedure SetupRTDx(bg:byte;x0,y0,xa,ya,xt,yt:Real;Var AE:Array of Real); var rd:Real; begin rd:=Sqrt(Sqr(xa-x0)+Sqr(ya-y0)); AE[bg+1]:=xa-x0; //dK1 AE[bg+2]:=(xa-x0)*Sqr(rd);//dK2 AE[bg+3]:=(xa-x0)*Sqr(Sqr(rd)); AE[bg+4]:=-(Sqr(rd)+2*Sqr(xa-x0)); AE[bg+5]:=-2*(xa-x0)*(ya-y0); AE[bg+6]:=xt-xa; end; Procedure SetupRTDz(bg:byte;x0,y0,xa,ya,xt,yt:Real;Var AE:Array of Real); var rd:Real; Begin rd:=Sqrt(Sqr(xa-x0)+Sqr(ya-y0)); AE[bg+1]:=ya-y0; AE[bg+2]:=(ya-y0)*Sqr(rd); AE[bg+3]:=(ya-y0)*Sqr(Sqr(rd)); AE[bg+4]:=-2*(xa-x0)*(ya-y0); AE[bg+5]:=-(Sqr(rd)+2*Sqr(ya-y0)); AE[bg+6]:=yt-ya; End; Procedure GiaiDH3(Matdat:Boolean;MaxNR:Word;Var LT,ndh:Word;Var f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk:Real;K1,K2,K3,P1,P2,df,dx0,dy0,dhk1,dh k2,dhk3,dhp1,dhp2:Real;Var OK:Boolean); Const IV=10; Var i,j,i1,m:Word; giaiOK,done:Boolean; xa,ya,X,Y,Z,PVV,Mo:Real; Begin //MessageBox(FrmKQtinh.Handle,PAnsiChar(IntToStr(IV)),PAnsiChar('DH '+IntToStr(ndh)),MB_OK+MB_ICONINFORMATION); LT:=0; OK:=False;done:=False; Repeat //Processing one Image until to satisfy sigma limit Inc(LT); InitialMT(NE); RotateMatrix(Matdat,phi,omega,kapa,AX); i:=0;ndh:=0; While i<=High(IC) do begin For j:=0 to High(KC) do if (IC[i].PN=KC[j].PN)and(IC[i].PELoai)and(KC[j].PEKTra) then //Start process begin xa:=IC[i].x;ya:=IC[i].y; CorrectDistortion(xa,ya,x0,y0,K1,K2,K3,P1,P2,xa,ya); X:=KC[j].X;Y:=KC[j].Y;Z:=KC[j].Z; SetupAEx(Matdat,xa,ya,X,Y,Z,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,AE); SetupNE(IV,1,AE,NE); SetupAEy(Matdat,xa,ya,X,Y,Z,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,AE); SetupNE(IV,1,AE,NE); Inc(ndh); end; Inc(i); end;//end while. i>High(IC) If (ndh*2)<IV then Exit;// Khong du diem de giai SolveNE(IV,True,NE,giaiOK,PVV); Mo:=0; if ndh*2>IV then Mo:=Sqrt(PVV/(ndh*2-IV)); GetSolution(IV,NE,AE,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P2); ER[0]:=Mo;//SSTP nhom 1 if giaiOK then // Tinh SS cua an begin InverseNE(IV,NE);// khong suy bien For i1:=1 to IV do begin m:=Trunc(i1*(i1-1)/2+i1); if NE[m]<0 then Exit;// Co SS tho ER[i1]:=Mo*SQRT(NE[m]); end; end; if LT<=MaxNR then done:=CompareSigmaLimit(IV,NE,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2) else Exit;// Co SS tho until done;//satisfy sigma limit for image. orientation image. OK:=True;// Neu giai thanh cong End; ////////////// Procedure GiaiMH5(Matdat:Boolean;MaxNR:Word;Var LT,nmh:Word;f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk:Real;Var K1,K2,K3,P1,P2:Real;df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2:Real;Var OK:Boolean); Const IV=5; Var i,j,i1,m:Word; giaiOK,done:Boolean; xa,ya,X,Y,Z,PVV,Mo,XN,YN,ZN,xt,yt:Real; Begin LT:=0; OK:=False;done:=False; RotateMatrix(Matdat,phi,omega,kapa,AX); Repeat //Processing one Image until to satisfy sigma limit Inc(LT); InitialMT(NE); i:=0;nmh:=0; While i<=High(IC) do begin For j:=0 to High(KC) do if (IC[i].PN=KC[j].PN)and(IC[i].PELoai)and(KC[j].PEKTra) then //Start process begin xa:=IC[i].x;ya:=IC[i].y; CorrectDistortion(xa,ya,x0,y0,K1,K2,K3,P1,P2,xa,ya); X:=KC[j].X;Y:=KC[j].Y;Z:=KC[j].Z; XYZtoXNYNZN(Matdat,X,Y,Z,f,x0,y0,XS,YS,ZS,hk,AX,XN,YN,ZN,xt,yt); SetupRTDx(0,x0,y0,xa,ya,xt,yt,AE); SetupNE(IV,1,AE,NE); SetupRTDz(0,x0,y0,xa,ya,xt,yt,AE); SetupNE(IV,1,AE,NE); Inc(nmh); end; Inc(i); end;//end while. i>High(IC) If nmh*2<IV then Exit;// Khong du diem de giai SolveNE(IV,True,NE,giaiOK,PVV); Mo:=0; if nmh*2>IV then Mo:=Sqrt(PVV/(nmh*2-IV)); // GetSolution(IV,NE,AE,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P2); ER[16]:=Mo;//SSTP nhom 2 if giaiOK then // Tinh SS cua an begin InverseNE(IV,NE);// khong suy bien For i1:=1 to IV do begin m:=Trunc(i1*(i1-1)/2+i1); if NE[m]<0 then Exit;// Co SS tho ER[10+i1]:=Mo*SQRT(NE[m]); end; end; if LT<=MaxNR then done:=CompareSigmaLimit(IV,NE,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2) else Exit;// Co SS tho until done;//satisfy sigma limit for image. orientation image. OK:=True;// Neu giai thanh cong End; //////////////////////// Procedure GiaiDT8(Matdat:Boolean;MaxNR:Word;Var LT,ndt:Word;Var f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P2:Real;df,dx0,dy0,dhk1,dh k2,dhk3,dhp1,dhp2:Real;Var OK:Boolean); Const IV=15; Var i,j,i1,m,LT3,LT5,ndh,nmh:Word; giaiOK,done,OK3,OK5:Boolean; xa,ya,X,Y,Z,PVV,Mo,XN,YN,ZN,xt,yt,dK1,dK2,dK3,dP1,dP2:Real; Begin LT:=0; OK:=False;done:=False; Repeat //Processing one Image until to satisfy sigma limit Inc(LT); dK1:=K1;dK2:=K2;dK3:=K3;dP1:=P1;dP2:=P2; GiaiDH3(Matdat,MaxNR,LT3,ndh,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3, P1,P2,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,OK3); GiaiMH5(Matdat,MaxNR,LT5,nmh,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3, P1,P2,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,OK5); dK1:=dK1-K1;dK2:=dK2-K2;dK3:=dK3-K3;dP1:=dP1-P1;dP2:=dP2-P2; InitialMT(NE); RotateMatrix(Matdat,phi,omega,kapa,AX); i:=0;ndt:=0; While i<=High(IC) do begin For j:=0 to High(KC) do if (IC[i].PN=KC[j].PN)and(IC[i].PELoai)and(KC[j].PEKTra) then //Start process begin xa:=IC[i].x;ya:=IC[i].y; CorrectDistortion(xa,ya,x0,y0,K1,K2,K3,P1,P2,xa,ya); X:=KC[j].X;Y:=KC[j].Y;Z:=KC[j].Z; XYZtoXNYNZN(Matdat,X,Y,Z,f,x0,y0,XS,YS,ZS,hk,AX,XN,YN,ZN,xt,yt); SetupAEx(Matdat,xa,ya,X,Y,Z,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,AE); SetupRTDx(10,x0,y0,xa,ya,xt,yt,AE); SetupNE(IV,1,AE,NE); SetupAEy(Matdat,xa,ya,X,Y,Z,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,AE); SetupRTDz(10,x0,y0,xa,ya,xt,yt,AE); SetupNE(IV,1,AE,NE); Inc(ndt); end; Inc(i); end;//end while. i>High(IC) SetupMRMH(dK1,dK2,dK3,dP1,dP2,NE); If (ndt*2+5)<IV then Exit;// Khong du diem de giai SolveNE(IV,True,NE,giaiOK,PVV); Mo:=0; if (ndt*2+5)>IV then Mo:=Sqrt(PVV/(ndt*2+5-IV)); GetSolution(IV,NE,AE,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P2); ER[16]:=Mo;//SSTP ER[0]:=Mo; if giaiOK then // Tinh SS cua an begin InverseNE(IV,NE);// khong suy bien For i1:=1 to IV do begin m:=Trunc(i1*(i1-1)/2+i1); if NE[m]<0 then Exit;// Co SS tho ER[i1]:=Mo*SQRT(NE[m]); end; end; if LT<=MaxNR then done:=CompareSigmaLimit(IV,NE,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2) else Exit;// Co SS tho until done;//satisfy sigma limit for image. orientation image. OK:=True;// Neu giai thanh cong End; //////////////////////// Function LoaiSSnho(Matdat:Boolean;f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P 2,DxzErr:Real):Word; Var i,j,L:Word; xa,ya,X,Y,Z,XN,YN,ZN,xt,yt:Real; Begin Result:=0;L:=0; RotateMatrix(Matdat,phi,omega,kapa,AX); i:=0; While i<=High(IC) do begin For j:=0 to High(KC) do if (IC[i].PN=KC[j].PN)and(KC[j].PEKTra) then //Start process begin xa:=IC[i].x;ya:=IC[i].y;IC[i].PE:=''; CorrectDistortion(xa,ya,x0,y0,K1,K2,K3,P1,P2,xa,ya); X:=KC[j].X;Y:=KC[j].Y;Z:=KC[j].Z; XYZtoXNYNZN(Matdat,X,Y,Z,f,x0,y0,XS,YS,ZS,hk,AX,XN,YN,ZN,xt,yt); if (ABS(xt-xa)>dxzErr)or(ABS(yt-ya)>dxzErr) then begin IC[i].PE:=Loai;Inc(L);end; end; Inc(i); end;//end Result:=L; End; Procedure TMyThread.Calculation; const Khongdu='Kh«ng ®ñ ®iÓm tÝnh'; giaiPTC='Gi¶i phu¬ng tr×nh chuÈn...'; NgichdaoPTC='NghÞch ®¶o phu¬ng tr×nh chuÈn...'; TinhSS='TÝnh sai sè cña Èn...'; DHanh='§Þnh huíng ¶nh ...#'; TinhHSMH='TÝnh hÖ sè mÐo h×nh...#'; Dongthoi='§H ¶nh ®ång thêi TÝnh hÖ sè mÐo h×nh...'; SSTP='SSTP mo ...'; Label 1,2,3; var beginLine,stopLine,MaxNR,ndh,nmh,ndt,LL,LT3,LT5,LT8,nLoai,nLoaiP:Word; TypeFile,TypeCalMethod,IV,Nimage:Byte; ImageName:string; EndSLD,LoaiErr,KiemDinh,OK3,OK5,OK8:Boolean; df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,dxzErr, f,x0,y0,XS,YS,ZS,hk,phi,omega,kapa,K1,K2,K3,P1,P2:Real; BEGIN TinhXong:=False; //Lay kieu tinh, han sai GetSigmaLimit(TypeCalMethod,MaxNR,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,dxz Err,LoaiErr); if FrmCamera.RGkdmh.ItemIndex=0 then KiemDinh:=True else KiemDinh:=False; if (TypeCalMethod=0) then IV:=10 else IV:=15; // MessageBox(FrmKQtinh.Handle,PAnsiChar(IntToStr(IV)),'gan',MB_OK+MB_ICONIN FORMATION); ER[16]:=-99; /// //Read image control points file With FrmKQtinh.REtientrinh.Lines do Begin if Matdat then Add('D¹ng chôp ¶nh mÆt ®Êt.......') else Add('D¹ng chôp ¶nh hµng kh«ng.......'); Add(TimeToStr(now)+' -> KiÓm tra ®Þnh d¹ng sè liÖu KC kiÓm ®Þnh...'); if FrmMain.CheckF1BookKC > 0 then begin Add(TimeToStr(now)+' -> Lçi ®Þnh d¹ng ®iÓm khèng chÕ kiÓm ®Þnh');Exit;end; Add(TimeToStr(now)+' -> KiÓm tra ®Þnh d¹ng sè liÖu ®o täa ®é ¶nh...'); if FrmMain.CheckF1BookSLD > 0 then begin Add(TimeToStr(now)+' -> Lçi ®Þnh d¹ng sè liÖu ®o täa ®é ¶nh');Exit;end; Add(TimeToStr(now)+' -> §äc sè liÖu khèng chÕ kiÓm ®Þnh...');ReadKC; if not CheckErrKC(IV)then Exit; //Read multi images coordinate data file Nimage:=0;beginLine:=1; Repeat //multi image Add(TimeToStr(now)+' -> §äc täa ®é ¶nh...');ReadICfile(beginLine,ImageName,stopLine,EndSLD); Add(TimeToStr(now)+' -> ¶nh "'+ImageName+'" Sè ®iÓm ®o '+IntToStr(High(IC)+1)); Add(TimeToStr(now)+' -> Khëi t¹o Èn sè...');InitialVar(f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P2); Inc(Nimage); //orientation image LL:=0;nLoai:=0;nLoaiP:=0; Repeat 2:Inc(LL);if IV=10 then Add(TimeToStr(now)+' -> TÝnh riªng lÎ : §Þnh huíng ¶nh...'+IntToStr(LL)) else Add(TimeToStr(now)+' -> TÝnh s¬ bé...'); GiaiDH3(Matdat,MaxNR,LT3,ndh,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3, P1,P2,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,OK3); if (OK3 and LoaiErr) then nLoai:=LoaiSSnho(Matdat,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3,P1,P2 ,dxzErr); //Loai SS nho //Add(TimeToStr(now)+' -> lo¹i...'+IntToStr(nLoai)); if not OK3 then begin Add(TimeToStr(now)+' -> "'+ImageName+'"...Kh«ng thÓ xö lý. Bá qua'); Nimage:=Nimage-1;Goto 1;end; if ((LoaiErr)and((nLoai-nLoaiP)0)and(LL<MaxNR)) then begin nLoaiP:=nLoai; if ndh*2 Sè ®iÓm bÞ lo¹i...'+IntToStr(nLoai)+'...Kh«ng ®ñ ®iÓm ®Ó tÝnh');Nimage:=Nimage- 1;Goto 1;end; Goto 2; end; // Loai SS nho if (not KiemDinh)and(OK3) then begin DislaySolution(Matdat,beginLine,Nimage,ImageName,AE,ER); Goto 1;end; if IV=10 then Add(TimeToStr(now)+' -> TÝnh riªng lÎ : HÖ sè mÐo h×nh ...'+IntToStr(LL)); GiaiMH5(Matdat,MaxNR,LT5,nmh,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3, P1,P2,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,OK5); if IV=15 then begin if (LoaiErr and((ndh*2+5) Sè ®iÓm bÞ lo¹i...'+IntToStr(nLoai)+'...Kh«ng ®ñ ®iÓm ®Ó tÝnh');Nimage:=Nimage-1;Goto 1;end; Goto 3;//Dong thoi end; until ((LL>=MaxNR)or(LT3<=2)); // Rieng le if OK3 and OK5 then begin DislaySolution(Matdat,beginLine,Nimage,ImageName,AE,ER);Goto 1;end else begin Nimage:=Nimage-1;Goto 1;end; //Dong thoi 3:Add(TimeToStr(now)+' -> TÝnh ®ång thêi: §H ¶nh vµ hÖ sè mÐo h×nh..."'+ImageName+'"'); GiaiDT8(Matdat,MaxNR,LT8,ndt,f,x0,y0,XS,YS,ZS,phi,omega,kapa,hk,K1,K2,K3, P1,P2,df,dx0,dy0,dhk1,dhk2,dhk3,dhp1,dhp2,OK8); if OK8 then begin DislaySolution(Matdat,beginLine,Nimage,ImageName,AE,ER);Goto 1;end else begin Nimage:=Nimage-1;Goto 1;end; 1:beginLine:=StopLine+1;//start new image Until EndSLD;//End of Data Add(TimeToStr(now)+' -> TÝnh to¸n xong!'); End;//end with TinhXong:=True; END; end. Phô lôc 8 ¶nh dïng kiÓm ®Þnh chôp b»ng m¸y Sony F707 Phô lôc 9 ¶nh chôp khu §«ng Nam (khu 1) má Cäc S¸u ngµy 11 - 10 - 2003 ¶nh chôp khu §«ng Nam (khu 1) má Cäc S¸u ngµy 20 - 02 - 2004 Phô lôc 10 ¶nh chôp khu §«ng B¾c (khu 2) má Cäc S¸u ngµy 11 - 10 - 2003 ¶nh chôp khu §«ng B¾c (khu 2) má Cäc S¸u ngµy 20 - 02 - 2004 Phô lôc 19 CÆp ¶nh lËp thÓ chôp Nhµ H¸t Lín Hµ Néi ngµy 20-7-2005 (Chôp b»ng m¸y ¶nh sè Sony F707, ®é ph©n gi¶i 2560x1920 pixel) B¶n vÏ kiÕn tróc mÆt tiÒn Nhµ H¸t Lín Hµ Néi (§o vÏ b»ng phÇn mÒm Photomodeler Pro5, biªn tËp trªn AutoCad 2004)