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

ameng2009

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

clear,clc
%读DEM文件
[file,path]=uigetfile('*.dem','选择文件');
if isequal(path,0) | isequal(file,0)
   errordlg('请选择要打开的数据文件','打开文件错误');
   return
else
   [fid,message]=fopen([path,file],'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

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

ameng2009

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

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

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

ydlcsu

save filename mydem1

ameng2009

谢谢你 楼上的  回去试试:lol