การพลอต Shuttle Radar Topography Mission (SRTM) map by using MATLab
SRTM data เป็นข้อมูล digital terrain models (DEM) รูปแบบหนึ่ง ซึ่งสามารถ download ได้ฟรี
อ้างอิงตาม
clear; close all ; clc;
latmin = 13-1;
latmax = 16;
lonmin = 98;
lonmax = 99;
fid = fopen('N14E098.hgt','r');
SRTM1 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N14E099.hgt','r');
SRTM2 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N13E098.hgt','r');
SRTM3 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N13E099.hgt','r');
SRTM4 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N15E098.hgt','r');
SRTM5 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N15E099.hgt','r');
SRTM6 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N16E098.hgt','r');
SRTM7 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
fid = fopen('N16E099.hgt','r');
SRTM8 = fread(fid,[1201,inf],'int16','b');
fclose(fid);
SRTM1 = SRTM1'; SRTM1 = flipud(SRTM1);
SRTM2 = SRTM2'; SRTM2 = flipud(SRTM2);
SRTM3 = SRTM3'; SRTM3 = flipud(SRTM3);
SRTM4 = SRTM4'; SRTM4 = flipud(SRTM4);
SRTM5 = SRTM5'; SRTM5 = flipud(SRTM5);
SRTM6 = SRTM6'; SRTM6 = flipud(SRTM6);
SRTM7 = SRTM7'; SRTM7 = flipud(SRTM7);
SRTM8 = SRTM8'; SRTM8 = flipud(SRTM8);
SRTM1(SRTM1 == -32768) = NaN;
SRTM2(SRTM2 == -32768) = NaN;
SRTM3(SRTM3 == -32768) = NaN;
SRTM4(SRTM4 == -32768) = NaN;
SRTM5(SRTM5 == -32768) = NaN;
SRTM6(SRTM6 == -32768) = NaN;
SRTM7(SRTM7 == -32768) = NaN;
SRTM8(SRTM8 == -32768) = NaN;
SRTM1(:,1201) = [];
SRTM3(:,1201) = [];
SRTM5(:,1201) = [];
SRTM7(:,1201) = [];
SRTM1(1201,:) = [];
SRTM2(1201,:) = [];
SRTM7(1201,:) = [];
SRTM8(1201,:) = [];
SRTM5(1201,:) = [];
SRTM6(1201,:) = [];
% SRTM = horzcat(SRTM1,SRTM2; SRTM3,SRTM4);
SRTM = [ SRTM3,SRTM4; SRTM1,SRTM2; SRTM5,SRTM6; SRTM7,SRTM8;];
clear SRTM1 SRTM2 SRTM3 SRTM4 SRTM5 SRTM6 SRTM7 SRTM8
for i = 2 : 1200*4
for j = 2 : 2400
if isnan(SRTM(i,j))
SRTM(i,j) = nanmean(nanmean(...
SRTM(i-1:i+1,j-1:j+1)));
end
end
end
clear i j
B = 1/25 * ones(5,5);
SRTM = filter2(B,SRTM);
[LON,LAT] = meshgrid(lonmin:1/1200:lonmax+1,latmin:1/1200:latmax);
v = -100:500:2000;
% v = [700 800 900 1000 1100 1200 ...
% 1300 1500 2000 2500 3000];
%% plot
%
figure1 = figure('Color',[1 1 1],...
'Position',[50 50 1200 600]);
axes1 = axes('Visible','off',...
'Units','centimeters',...
'FontSize',8);
hold(axes1,'all');
surf1 = surf(LON,LAT,SRTM,...
'SpecularExponent',20,...
'FaceLighting','phong',...
'FaceColor','interp',...
'EdgeColor','none');
light1 = light('Parent',axes1,...
'Style','local',...
'Position',[145 70 900000]);
set(gca,'View',[25 20])
daspect([1 1 20000])
% colormap(flipud(hsv))
% zlimits = [-100 2000];
% demcmap(zlimits);
demcmap(SRTM)
colorbar;
% contour3(LON,LAT,SRTM,v,...
% 'Color','w')
โฆษณา
- ดาวน์โหลดแอปพลิเคชัน
- © 2025 Blockdit
