Tìm hiểu cơ sở lý thuyết về trọng lực và các phương pháp hiệu chỉnh trọng lực.Xây dựng quy trình tính toán hiệu chỉnh trọng lực bao gồm: hiệu chỉnh khoảng không (Freeair reduction), hiệu chỉnh Bouguer (Bouguer reduction), hiệu chỉnh địa hình (Terrain correction), hiệu chỉnh đẳng tĩnh địa hình (TopographicIsostatic reduction) theo mô hình của AiryHeiskanen và hiệu chỉnh ảnh hưởng gián tiếp (Indirect effect).Xây dựng chương trình tính toán xử lý số liệu đo trọng lực viết bằng ngôn ngữ Matlab nhằm hiệu chỉnh trọng lực theo quy trình đã đưa ra.
ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA -oOo -* BK ĐẶNG NGỌC DANH TÊN ĐỀ TÀI: XÂY DỰNG CHƯƠNG TRÌNH TÍNH HIỆU CHỈNH TRỌNG LỰC Chuyên ngành: Kỹ thuật Trắc địa - Bản đồ Mã số: 60520503 LUẬN VẪN THẠC SĨ TP HỒ CHÍ MINH, tháng năm 2018 m Cơng trình hồn thành tại: Trường Đại học Bách Khoa - ĐHQG-HCM Cán hướng dẫn khoa học: TS Lưong Bảo Bình Cán chấm nhận xét 1: PGS.TS Đào Xuân Lộc Cán chấm nhận xét 2: TS Đỗ Minh Tuấn Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp.HCM ngày 16 tháng 08 năm 2018 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: PGS.TS Nguyễn Ngọc Lâu PGS.TS Đào Xuân Lộc TS Đỗ Minh Tuấn PGS.TS Lê Trung Chon TS Phan Thị Anh Thư Xác nhận Chủ tịch Hội đồng đánh giá LV Trưởng Khoa quản lý chuyên ngành sau luận văn sửa chữa CHỦ TỊCH HỘI ĐÒNG TRƯỞNG KHOA ĐẠI HỌC QUỐC GIA TP.HỒ CHÍ MINH CỘNG HỊA XÃ HỘI CHỦ NGHĨA VIỆT NAM TRƯỜNG ĐẠI HỌC BÁCH KHOA Độc Lập - Tự Do - Hạnh Phúc LỜI CẢM ƠN NHIỆM VỤ LUẬN VĂN CAO HỌC Họ tên học : Đặng Ngọc Danh Tôi viên xin gửi lời cảm ơn chân thành đến: MSHV: 1570197 - TS.năm Lương Bình trực tiếp hướng dẫn,Nơi gópsinh: ý vàBình giúpThuận đỡ tơi Ngày, tháng, sinhBảo : 3/11/1987 suốt quángành trình thực :luận văn Trắc địa - Bản đồ Chuyên Kỹ thuật - Các Thầy cô giảng dạy lớp Cao học Trắc địa - Bản đồ cung cấp cho 60520503 Mã sồ kiến thức quý báu trình học tập Trường Đại học Bách Khoa I TÊN ĐỀ TÀI: Xây dựng chương trình tính hiệu chỉnh trọng lực Tp.HCM II MỤC TIÊU NGHIÊN CỨU: - Người thân gia đĩnh, bạn bè, đồng nghiệp giúp đỡ, động viên sởtơi lý hồn thuyếtthành trọng lực phương pháp hiệu chỉnh trọng lực tạoTìm điềuhiểu kiệncơ cho tốt luận văncác Xây dựng quy trình tính toán hiệu chỉnh trọng lực bao gồm: hiệu chỉnh khoảng Tp.HCM, ngày reduction), 01 tháng năm không (Free-air reduction), hiệu chỉnh Bouguer (Bouguer hiệu 2018 chỉnh địa hình (Terrain correction), hiệu chỉnh đẳng tĩnh địa hình (Topographic-Isostatic reduction) theo mơ hình Airy-Heiskanen hiệu chỉnh ảnh hưởng gián tiếp (Indirect effect) Đặng Ngọc Danh Xây dựng chương trình tính tốn xử lý số liệu đo trọng lực viết ngôn ngữ Matlab nhằm hiệu chỉnh trọng lực theo quy trình đưa III NGÀY GIAO NHIỆM VỤ: 26/02/2018 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 01/08/2018 V CÁN BỘ HƯỚNG DẪN: TS Lương Bảo Bình Tp.HCM, ngày thảng .năm 2018 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO TRƯỞNG KHOA TĨM TẮT Luận văn tìm hiểu lý thuyết trọng lực phương pháp hiệu chỉnh trọng lực bao gồm: hiệu chỉnh khoảng không, hiệu chỉnh Bouguer, hiệu chỉnh địa hình, hiệu chỉnh Prey, hiệu chỉnh đẳng tĩnh theo hai giả thuyết Pratt-Hayford Airy-Heiskanen, hiệu chỉnh ảnh hưởng gián tiếp, hiệu chỉnh nghịch đảo Rudzki, hiệu chỉnh cô đặc Helmert vấn đề đặt luận văn ứng dụng lý thuyết hiệu chỉnh trọng lực để tính toán hiệu chỉnh cho giá trị trọng lực đo bề mặt đất gồm hiệu chỉnh: hiệu chỉnh khoảng khơng, hiệu chỉnh Bouguer, hiệu chỉnh địa hình, hiệu chỉnh đẳng tĩnh địa hình theo mơ hình Airy-Heiskanen hiệu chỉnh ảnh hưởng gián tiếp Cuối tác giả xây dựng chương trĩnh tính tốn phần mềm Matlab chạy thực nghiệm cho chương trĩnh xây dựng SUMMARY This thesis is to find out about the theory of gravity and gravity reduction methods includes: Free-air reduction, Bouguer reduction, Terrain correction, Prey reduction, Isostatic reduction addcording to two hypotheses of Pratt-Hayford and Airy-Heiskanen, the indirect effect, the inversion reduction of Rudzki, the condensation reduction of Helmert The question of this thesis is using the theory of gravity reduction to correct for the gavity values measured on the surface of the earth includes: Free-air reduction, Bouguer reduction, Terrain correction, Isostatic reduction addcording to hypothesis of Airy-Heiskanen, the indirect effect Finally, the author build up a calculation program by Matlab software and run experiments for program built up LỜI CAM ĐOAN Tôi xin cam đoan đề tài “ X â y d ự n g c h n g t r ì n h t ỉ n h H i ệ u c h ỉ n h t r ọ n g l ự c ’ ' ’ TS Lương Bảo Bình hướng dẫn cơng trình nghiên cứu riêng tơi Các số liệu, kết luận văn trung thực chưa công bố công trình khác Tác giả luận văn Đặng Ngọc Danh MUC LUC •• CHƯƠNG 1: MỞ ĐẦU 1.1 Lý chọn đề tài 1.2 Các nghiên cứu liên quan 1.2.1 Tình hình nghiên cứu giới 1.2.2 Tình hình nghiên cứu nước 1.3 Mục tiêu nghiên cứu 1.4 Nội dung nghiên cứu 1.5 Phương pháp nghiên cứu 1.6 Đối tượng phạm vi nghiên cứu 1.7 Ý nghĩa khoa học thực tiễn „1 „1 „2 „2 5 6 CHƯƠNG 2: Cơ SỞ LÝ THUYẾT VỀ TRỌNG LỰC VÀ HIỆU CHỈNH TRỌNG LỰC 2.1 Lực hấp dẫn hấp dẫn 2.2 Trọng lực dị thường trọng lực 2.2.1 Trọng lực 2.2.2 Dị thường trọng lực 2.3 Hiệu chỉnh trọng lực (Gravity Reduction) 2.3.1 Các công thức bổ trợ .7 10 10 14 16 17 17 2.3.1.1 p nằm khối trụ 2.3.1.2 p nằm mặt khối trụ 19 2.3.1.3 p nằm khối trụ 19 2.3.1.4 Đĩa tròn 20 2.3.1.5 Một phần vành khuyên 20 2.3.2 Hiệu chỉnh khoảng không (Free-air reduction) .21 2.3.3 Hiệu chỉnh Bouguer (Bouguer reduction) 22 2.3.3.1 Đĩa Bouguer 22 2.3.3.2 Hiệu chỉnh địa hình (Terrain correction) 23 2.3.3.3 Phương pháp họp 24 2.3.4 Các giả thuyết đẳng tĩnh 25 2.3.4.1 Mơ hình đẳng tĩnh theo giả thuyết Pratt - Hayford 25 2.3.4.2 thống Aữy - Heiskanen 27 2.3.5 Hiệu chỉnh đẳng tĩnh địa hình (Topographic-isostatic reductions) 29 2.3.6 Hiệu chỉnh ảnh hưởng gián tiếp (The indữect effect) .30 2.3.7 Một số hiệu chỉnh tham khảo 32 2.3.7.1 Hiệu chỉnh Prey (Prey reduction) 32 23.1.2 Hiệu chỉnh nghịch đảo Rudzki (The inversion reduction of Rudzki) .34 2.3.7.3 Helmert) Hiệu chỉnh cô đặc Helmert (The condensation reduction of 36 CHƯƠNG 3: XÂY DỤNG CHƯƠNG TRÌNH HIỆU CHỈNH TRỌNG Lực 39 3.1 Khái quát chương trình hiệu chỉnh trọng lực (Gravity Reduction) 39 Hệ 3.2.3 3.3 Đối với liệu đầu vào dạng tọa độ trắc địa 45 3.4 Xây dựng chương trình tính 45 3.3.1 Thiết kế chương trình 45 3.3.2 Tổ chức file liệu 50 3.3.2.1 Dữ liệu file DTM 50 3.3.2.2 Dữ liệu đo trọng lực 51 Tính tốn thực nghiệm 52 3.4.1 Kết tính tốn chương trình 52 3.4.2 So sánh kết chương trình với phần mềm Excel 58 CHƯƠNG 4: KẾT LUẬN 62 4.1 vấn đề đạt .62 4.2 4.3 Khó khăn hạn chế 62 4.2.1 Khó khăn 62 4.2.2 Hạn chế .63 Hướng phát triển đề tài 63 CÁC TÀI LIỆU THAM KHẢO 64 Các LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % Lap ma tran cao gom 4*n*n hang va k cot cho moi diem luc for i=l:k HHt (:,i)=ht ( (i-1)*4*n*n+l:i*4*n*n, :); end % % Hieu chinh dia hình At for i=l:k for j=l:4*n*n HH_gach(j,i)=abs(HHt(j,i)-Hp(i)); DD1(j,i)=sqrt(((XXt (j,i)-XP(i))A2+(YYt(j,i)-YP(i))A2)+((HHt(j,i) Hp(i))/2)A2); % Tinh hieu chinh dia hinh cho moi don vi khoi luong Delta m dAAt(j,i)=(G*roO*(10A5)*((p*10A3)A2)*(HH_gach(j , i))A2)/(2*(DD1 (j,i)A3) ) ; end end % Tinh hieu chinh dia hinh AAt=sum(dAAt); % % Hieu chinh lop bu tru Ac % Tinh khoang each tu diem luc den tarn cac ROOT T=30000; t=4.45*HHt; % % Hieu chinh lop bu tru Ac for i=l:k for j=l:4*n*n % Tinh khoang each tu diem luc den tarn cac ROOT DD2 (j,i)=sqrt( ( (XXt(j,i)-XP(i))A2+(YYt(j,i)- YP (i))A2) + (HHt(j,i)+T+ (t(j,i))/2)A2); % Tinh anh huong cua phan bu tru cho moi don vi KL Delta m dAAc(j,i)=((G*dro*((p*10A3)A2)*t(j,i)*(Hp(i)+T+(t(j,i))/2))*10A5)/(DD2(j, i)A3); end end AAc=sum(dAAc); % Ban kinh tinh toan bk=n*p*10A3; % (m) % % Anh huong gian tiep for i=l:k for j =1:4*n*n % Tinh UT (m2/s2) dung uo duut(j,i)=pi*G*roO*((HHt(j,i)A2)+HHt(j,i)* sqrt(bkA2+HHt(j,i)A2) + (bkA2)*log( (HHt(j,i)+sqrt(bkA 2+HHt (j,i)A2) )/bk)) ; % Tinh Uc dung Ue dUUc(j,i)=pi*G*dro*(TA2- ( (t (j,i)+T)A2)- T*sqrt(bkA2+TA2)+(t(j,i)+T)*sqrt(bkA2+ (t(j,i)+T)A2)(bkA2)*log(+T+sqrt(bkA2+(+T)A2))+(bkA2)*log(t(j,i)+T+sqrt(bkA2+(t(j,i)+T) 2) ) ) ; end end uut=sum(duut); Trang 72 LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH UUc=sum(dUUc); % % Cac hieu chinh F=0.308 6*Hp; Ab=0.1119*Hp; At=transpose(AAt); Ac=transpose(AAc); DELTA_N=(UUt'-UUc')./(ỵ*10A-5); DELTA=+0.3Ũ86*DELTA_N; % (mgal) % % Cac di thuong % Di thuong khoang khong gO=g+F; delta_gF=gO-ỵ; % Di thuong Bouguer don gian gB_dg=g-Ab+F; de1ta_gB_dg=gB_dg-ỵ; % Di thuong Bouguer hoan chinh gB_hc=g-Ab+At+F; delta_gB_hc=gB_hc-ỵ; % Di thuong Dang tinh dia hinh AT=Ab-At; gti=g-AT+Ac+F; delta_gti=gti-ỵ; delta_gc=g-AT+Ac+F+DELTA-ỵ; kq=l; if AH1==0 for i=l: k F (i)=0;delta_gF(i)=0; end end if AH2==0 for i=l:k Ab(i)=0;delta_gB_dg(i)=0; end end if AH3==0 for i=l:k At(i)=0;delta_gB_hc(i)=0; end end if AH4==0 for i=l:k Ac (i)=0;delta_gti(i)=0; end end if AH5==0 for i=l:k DELTA(i)=0;delta_gc(i)=0; end end end Trang 73 LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % Goi ham xuat ket qua anh huong if TD1==1 && kq==l ketquatheo_XY(k,Xp,Yp,F,delta_gF,Ab,delta_gB_dg,At,delta_gB_hc,Ac,delta_g ti,DELTA,delta_gc) end if TD2==1 && kq==l ketquatheo_BL(k,Xp,Yp,F,delta_gF,Ab,delta_gB_dg,At,delta_gB_hc,Ac,delta_g ti,DELTA,delta_gc) end % % ket thuc tinh toan va mo file msgbox('File Results.txt was saved in the current f o l d e r I n f o r m a t i o n ' ) GravityReduction.m function varargout = GravityReduction(varargin) % GravitỵReduction MATLAB code for GravityReduction.fig % GravityReduction, by itself, creates a new GravityReduction or raises the existing % singleton* % % H = GravityReduction returns the handle to a new GravityReduction or the handle to % the existing singleton* % % GravityReduction('CALLBACK',hobject,eventData, handles, ) calls the local % function named CALLBACK in GravitỵReduction.M with the given input arguments % % G r a v i t y R e d u c t i o n ( ' P r o p e r t y ' , ' V a l u e ) creates a new GravityReduction or raises the % existing singleton* starting from the left, property value pairs are % applied to the GUI before GravityReduction_OpeningFcn gets called An % unrecognized property name or invalid value makes property application % stop All inputs are passed to GravityReduction_OpeningFcn via varargin o, o % one % o, o *See GUI Options on GUIDE'S Tools menu Choose "GUI allows only instance to run (singleton)" % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help GravityReduction % Last Modified by GUIDE v2.5 26-Jul-2018 09:51:41 Trang 74 LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % Begin initialization code - DO NOT gui_Singleton = 1; gui_state = struct('gui_Name', 'gui_Singleton', 'gui_OpeningFcn', 'gui_OutputFcn', 'gui_LayoutFcn', 'gui Callback', EDIT mfilename, gui_Singleton, @GravityReduction_OpeningFcn @GravityReduction_OutputFcn, [ ] , [] ) ; if nargin && ischar(varargin{1}) gui_state.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_state, varargin{:}); else gui_mainfcn(gui_state, varargin{:}); end % End initialization code - DO NOT EDIT % - - Executes just before GravityReduction is made visible function GravityReduction_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % varargin command line arguments to GravityReduction (see VARARGIN) % Choose default handles.output = % Update handles guidata(hObject, command line output for GravityReduction hObject; structure handles); % UIWAIT makes GravityReduction wait for user response (see UIRESUME) % uiwait(handles.figurel); % - - Outputs from this function are returned to the command line function varargout = GravityReduction_OutputFcn(hobject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Get default command line output from handles structure varargoutfl} = handles.output; % - - Executes on selection change in popupmenul Trang 75 LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH function popupmenul_Callback(hObject, eventdata, handles) % hObject handle to popupmenul (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: contents = cellstr(get(hObject,'string')) returns popupmenul contents as cell array % contents{get(hObject,'Value')} returns selected item from popupmenul % - - Executes during object creation, after setting all properties function popupmenul_CreateFcn(hObject, eventdata, handles) % hObject handle to popupmenul (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: popupmenu controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isegual(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function duongdanl_Callback(hObject, eventdata, handles) % hObject handle to duongdanl (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'string') returns contents of duongdanl as text % str2double(get(hObject,'String')) returns contents of duongdanl as a double % - - Executes during object creation, after setting all properties function duongdanl_CreateFcn(hObject, eventdata, handles) % hObject handle to duongdanl (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isegual(get (hObject, 'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end Trang 76 LUẬN VĂN TỐT NGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % - - Executes on button press in Chon_l function Chon_l_Callback(hObject, eventdata, handles) % hObject handle to Chon_l (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [ten DD]=uigetfile('*.txt','Select File'); if isequal (ten,0) II isequal(DD,0) dispt'User pressed cancel') else full=strcat (DD,ten); disp(['User selected ',full]) set(handles.duonqdanl,'strinq',full); qlobal A A=textread(full); end % - - Executes on button press in Chon_2 function Chon_2_Callback(hObject, eventdata, handles) % hObject handle to Chon_2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) [ten DD]=uigetfile('*.txt','Select File'); if isequal (ten,0) II isequal(DD,0) dispt'User pressed cancel') else full=strcat (DD,ten); disp(['User selected ',full]) set(handles.duonqdan2,'strinq',full); qlobal B B=textread(full); end % - - Executes on button press in Tinhtoan function Tinhtoan_Callback(hobject, eventdata, handles) % hObject handle to Tinhtoan (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % ssss=str2num(qet(handles.duonqdanl,'strinq')); qlobal A B % Goi ham tinh cacloai_hieuchinh_va_dithuonq % - - Executes on button press in chk_AHl function chk_AHl_Callback(hObject, eventdata, handles) % hObject handle to chk_AHl (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) Trang 77 LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % Hint: get(hObject,'Value') returns toggle state of chk_AHl % - - Executes on button press in chk_AH2 function chk_AH2_Callback(hObject, eventdata, handles) % hObject handle to chk_AH2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of chk_AH2 % - - Executes on button press in chk_AH3 function chk_AH3_Callback(hObject, eventdata, handles) % hObject handle to chk_AH3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject,'Value') returns toggle state of chk_AH3 % - - Executes on button press in chk_AH4 function chk_AH4_Callback(hObject, eventdata, handles) % hObject handle to chk_AH4 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of chk_AH4 % - - Executes on button press in chk_AH5 function chk_AH5_Callback(hObject, eventdata, handles) % hObject handle to chk_AH5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject, 'Value1) returns toggle state of chk_AH5 function duongdan2_Callback(hObject, eventdata, handles) % hObject handle to duongdan2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject, String') returns contents of duongdan2 as text % str2double(get(hObject,'string')) returns contents of duongdan2 as a double % - - Executes during object creation, after setting all properties function duongdan2_CreateFcn(hObject, eventdata, handles) % hObject handle to duongdan2 (see GCBO) Trang 78 LUẬN VĂN TỐT NGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal (get(hObject, 'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - - Executes during object deletion, before destroying properties function duongdan2_DeleteFcn(hObject, eventdata, handles) % hObject handle to duongdan2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % - - Executes during object creation, after setting all properties function Chon_l_CreateFcn(hObject, eventdata, handles) % hObject handle to Chon_l (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called function canh_Callback(hObject, eventdata, handles) % hObject handle to canh (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get (hObject,'string') returns contents of canh as text % str2double(get(hObject,'string')) returns contents of canh as a double % - - Executes during object creation, after setting all properties function canh_CreateFcn(hObject, eventdata, handles) % hObject handle to canh (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isequal(get(hObject,1BackgroundColor1), get(0,'defaultUicontrolBackgroundColor1)) set(hObject,1BackgroundColor1,1 white 1); end Trang 79 LUẬN VĂN TỐT NGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH function bankinh_Callback(hObject, eventdata, handles) % hObject handle to bankinh (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'string') returns contents of bankinh as text % str2double(get(hObject,'string')) returns contents of bankinh as a double % - - Executes during object creation, after setting all properties function bankinh_CreateFcn(hObject, eventdata, handles) % hObject handle to bankinh (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows % See ISPC and COMPUTER if ispc && isegual(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % - - Executes on button press in Close function Close_Callback(hObject, eventdata, handles) % hObject handle to Close (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) ar=guestdlg('Do you want exit the p r o g r a m ? E x i t Y e s N o Y e s ' ) ; if strcmp(ar,'Yes') close; end % - - Executes on button press in graphic function graphic_Callback(hObject, eventdata, handles) % hObject handle to graphic (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)z global LC1 LC2 LCl=get(handles.XY, 'value'); LC2=get(handles.BL, 'value'); run Graphic.m Graphic.m function varargout = Graphic(varargin) % GRAPHIC MATLAB code for Graphic.fig % GRAPHIC, by itself, creates a new GRAPHIC or raises the existing % singleton* Trang 80 LUẬN VĂN TÔTNGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % % % H = GRAPHIC returns the handle to a new GRAPHIC or the handle to the existing singleton* % % GRAPHIC('CALLBACK',hObject,eventData,handles, ) calls the local % function named CALLBACK in GRAPHIC.M with the given input arguments % % GRAPHIC('Property','Value', ) creates a new GRAPHIC or raises the % existing singleton* starting from the left, property value pairs are % applied to the GUI before Graphic_OpeningFcn gets called An % unrecognized property name or invalid value makes property application % stop All inputs are passed to Graphic_OpeningFcn via varargin % % *See GUI Options rdt_gra GUIDE'S Tools menu Choose "GUI allows only one % instance to run (singleton)" % % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help Graphic % Last Modified by GUIDE v2.5 30-Jul-2018 08:24:41 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_state = struct('gui_Name', mfilename, 'gui_Singleton', gui_Singleton, 'gui_OpeningFcn', @Graphic_OpeningFcn, 'gui_OutputFcn', @Graphic_OutputFcn, 'gui_LayoutFcn', [] , 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_state.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_state, varargin{:}); else gui_mainfcn(gui_state, varargin{:}); end % End initialization code - DO NOT EDIT % - - Executes during object creation, after setting all properties function axes2_CreateFcn(hObject, eventdata, handles) Trang 81 LUẬN VĂN TÔTNGHIỆP THẠC sĩ % hObject % eventdata % handles called HVTH: ĐẶNG NGỌC DANH handle to axes2 (see GCBO) reserved - to be defined in a future version of MATLAB empty - handles not created until after all CreateFcns % Hint: place code in OpeningFcn to populate axes2 % - - Executes just before Graphic is made visible function Graphic_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % varargin command line arguments to Graphic (see VARARGIN) % Choose default command line output for Graphic handles.output = hObject; global A B LC1 LC2 % Update handles structure guidata(hObject, handles); % UIWAIT makes Graphic wait for user response (see UIRESUME) % uiwait(handles.figurel); % - - Outputs from this function are returned to the command line function varargout = Graphic_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Get default command line output from handles structure varargoutfl} = handles.output; % - - Executes rdt_gra button press in rdt_l function rdt_l_Callback(hObject, eventdata, handles) % hObject handle to rdt_l (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject, 'Value1) returns toggle state of rdt_l % - - Executes rdt_gra button press in rdt_2 function rdt_2_Callback(hObject, eventdata, handles) % hObject handle to rdt_2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB Trang 82 LUẬN VĂN TỐT NGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject,'Value') returns toggle state of rdt_2 % - - Executes rdt_gra button press in rdt_3 function rdt_3_Callback(hObject, eventdata, handles) % hObject handle to rdt_3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject,'Value') returns toggle state of rdt_3 %- - -Executes rdt_gra button press in checkbox'! function checkbox4_Callback(hObject, eventdata, handles) % hObject handle to checkbox'! (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject, 'Value') returns toggle state of checkbox'! % - - Executes on button press in chk_GRID function chk_GRID_Callback(hobject, eventdata, handles) % hObject handle to chk_GRID (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) %Lay gia tri check box % Hint: get (hObject,'Value') returns toggle state of chk_GRID % - - Executes on button press in chk_TEN function chk_TEN_Callback(hObject, eventdata, handles) % hObject handle to chk_TEN (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get (hObject, 'Value1) returns toggle state of chk_TEN % - - Executes rdt_gra button press in close function close_Callback(hObject, eventdata, handles) % hObject handle to close (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) close % - - Executes rdt_gra button press in Done function Done_Callback(hObject, eventdata, handles) % hObject handle to Done (see GCBO) Trang 83 LUẬN VĂN TỐT NGHIỆP THẠC sĩ HVTH: ĐẶNG NGỌC DANH % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global A B LC1 LC2 [ k, 1, Xp, Yp, Hp, g]=docfile_trongluc(A); [r, s, X, Y, H] =docfile_DTM (B) ; RDTl=get(handles.rdt_l, 'value'); RDT2=get(handles.rdt_2, 'value'); RDT3=get(handles.rdt_3, 'value'); GRAD3D=get(handles.GRA3D, 'value'); DTM3D=get(handles.DTM3D, 'value'); GRID=get(handles.chk_GRID, 'value'); TEN=get(handles.chk_TEN, 'value'); TEN_GRA=[l:k]'; TEN_DTM=[l:r]'; % Chon he toa if LC1==1 XP=Xp;YP=Yp;XD=X;YD=Y; end if LC2==1 [Xp_X,Yp_Y]=BL2XY(A); [X_X,Y_Y]=BL2XY(B); XP^lamtronso(Xp_X,4); YP=lamtronso(Yp_Y,4); XD^lamtronso(X_X,4); YD^lamtronso(Y_Y,4); end % Xac dinh tong so diem luc nam ngoai khu vuc tinh Dn= [ ]; for i=l:k if XP(i)>max(XD) II YP(i)>max(YD) II XP(i)