ameng2009 发表于 2009-3-22 14:47

DEM数据怎么截取一小块数据另存出来啊

DEM数据怎么截取一小块数据另存出来啊源文件太大

clear,clc
%读DEM文件
=uigetfile('*.dem','选择文件');
if isequal(path,0) | isequal(file,0)
   errordlg('请选择要打开的数据文件','打开文件错误');
   return
else
   =fopen(,'r');
end
for i=1:12
   mam{i}=fgetl(fid);
end
dem_unit=mam{3};%单位
start_x=str2num(mam{4});%列向坐标起点
start_y=str2num(mam{5});%行向坐标起点
totle_x=str2num(mam{6});%列向总长
totle_y=str2num(mam{7});%行向总长
dx=str2num(mam{8});      %列向采样间隔
dy=str2num(mam{9});      %行向采样间隔
volumn_num=str2num(mam{10}); %数据列数
row_num=str2num(mam{11});    %数据行数
mydem=zeros(row_num,volumn_num);
for i=1:row_num
    myline=str2num(fgetl(fid));
    ll=length(myline);
    mydem(i,1:ll)=myline;
end
mydem1=mydem(row_num-650:row_num-350,volumn_num-650:volumn_num-350);显示的这一部分怎么单独存储出来啊
fclose(fid)

ChaChing 发表于 2009-3-22 22:44

首先DEM数据是什麼? 不清楚
還有個人水平有限, 建議樓主看下本版規則並加強發問題方式!

ameng2009 发表于 2009-3-24 13:51

DEM是数字高程模型数据

DEM是数字高程模型数据,在下水平洼,楼上见谅:@$

:@D 实际就是很大的数字矩阵有901*1023这么大,我知道要取得是:例如650行-----550行,750列到650列,怎么把这一部分另存为一个文件,以后直接用matlab直接读这个文件即可出地形图等等。

[ 本帖最后由 ChaChing 于 2009-3-24 15:23 编辑 ]

ydlcsu 发表于 2009-3-24 14:00

回复 楼主 ameng2009 的帖子

save filename mydem1

ameng2009 发表于 2009-3-24 14:05

谢谢你 楼上的

谢谢你 楼上的回去试试:lol
页: [1]
查看完整版本: DEM数据怎么截取一小块数据另存出来啊