map工具箱显示问题

lhm901112

在把计算数据贴到地图上时遇到的奇怪问题，循环显示多幅图像，但奇怪的是第一幅数据不能显示，后面的都可以正常显示，求大神指点迷津，程序附后（小菜鸟写的请多了解哈

1. 
   
2. listing = dir('I:\new_data\Insar\Ascending\*.dat');
   
3. path = cd;
   
4. for p = 1:2
   
5.    filename = listing(p).name;
   
6.    U_x = dlmread(filename);
   
7.    file = strcat(filename(1:end-3),'N1');
   
8.    file = strcat(path,'\',file);
   
9.    fid = fopen(file,'r','b');
   
10.    
    
11.    fseek(fid,5314,'bof');
    
12.    a = fscanf(fid,'%c',10);
    
13.    number = str2num(a);
    
14. 
    
15.     %length
    
16.    fseek(fid,5335,'bof');
    
17.    a=fscanf(fid,'%c',10);
    
18.    length=str2num(a);
    
19.    length=(length-17)/4;
    
20. 
    
21.    fseek(fid,4546+133,'bof');
    
22.    a=fscanf(fid,'%c',21);
    
23.    offset=str2num(a);
    
24.    
    
25.    sample_num=zeros(12,11);
    
26.    lats=zeros(12,11);
    
27.    lons=zeros(12,11);
    
28.    line_num=zeros(12,1);
    
29.    
    
30.    for j=1:11
    
31.        fseek(fid,offset+(j-1)*521+13,'bof');
    
32.        line_num(j)=fread(fid,1,'uint32');
    
33.        fseek(fid,8,'cof');
    
34.        sample_num(j,:)=fread(fid,11,'uint32');
    
35.        fseek(fid,88,'cof');
    
36.        lats(j,:)=fread(fid,11,'int32')/1e6;
    
37.        lons(j,:)=fread(fid,11,'int32')/1e6;
    
38.    end
    
39.    
    
40.    fseek(fid,521+34,'cof');
    
41.    line_num(12)=number;
    
42.    sample_num(12,:)=fread(fid,11,'uint32');
    
43.    fseek(fid,88,'cof');
    
44.    lats(12,:)=fread(fid,11,'int32')/1e6;
    
45.    lons(12,:)=fread(fid,11,'int32')/1e6;
    
46.    
    
47.    x=sample_num;y=repmat(line_num,1,11);
    
48.    xi=repmat(1:length,number,1);
    
49.    yi=repmat((1:number)',1,length);
    
50.    lat=interp2(x,y,lats,xi,yi);
    
51.    lon=interp2(x,y,lons,xi,yi);
    
52.    clear xi & yi;
    
53.    
    
54.    lat=downsample(lat,floor(number/length)^2);lat=downsample(lat',floor(number/length));lat=lat';
    
55.    lon=downsample(lon,floor(number/length)^2);lon=downsample(lon',floor(number/length));lon=lon';
    
56.   
    
57.    
    
58.    % pass title
    
59.    fseek(fid,1946,'bof');
    
60.    a=fscanf(fid,'%c',10);
    
61.    ind=strcmp(a,'DESCENDING');
    
62.    if ind
    
63.        lat=fliplr(lat);
    
64.        lon=fliplr(lon);
    
65.        %U_x=fliplr(U_x);
    
66.    else
    
67.        lat=flipud(lat);
    
68.        lon=flipud(lon);
    
69.        %U_x=flipud(U_x);
    
70.    end
    
71.    
    
72.    figure;
    
73.    clf;
    
74.    hold on
    
75.    m_proj('UTM','lon',[min(min(lon))-0.2 max(max(lon))+0.2],'lat',[min(min(lat))-0.2 max(max(lat))+0.2]);         
    
76.    m_grid('xlabeldir','middle','fontsize',10,'tickdir','in','box','fancy','ylabeldir','middle');   
    
77.    h=m_pcolor(lon,lat,U_x);
    
78.    m_gshhs_f('patch',[.7 .7 .7]);
    
79.    set(h,'edgecolor','none')
    
80.    caxis([-1.0 1.0])
    
81.    
    
82.    h=colorbar('v');
    
83.    set(get(h,'title'),'string','m/s','fontsize',10);
    
84.    set(h,'pos',[0.8570 0.1560 0.02 0.7440])
    
85. 
    
86.    
    
87. %    clearvars -except p listing
    
88.    
    
89. end

复制代码

lhm901112

C:\Users\lihuimin\Desktop\QQ截图20140227193531.jpg
图没上去