请教各位,如下代码报错

singtaoo

a=shaperead('bou2_4p.shp');%读取shp文件的内容
bou2_4lx=[a(:).X];%提取经度信息
bou2_4ly=[a(:).Y];%提取纬度信息
m_proj('Lambert Conformal Conic','lon',[70,140],'lat',[0,60])%选择投影方式
m_plot(bou2_4lx,bou2_4ly)%绘图
m_grid%添加坐标
hold on
lat=[116.41667 117.2 114.48333 112.53333 111.68 ...
    123.38333 125.35 126.63333 121.43333 118.78333 ...
    120.2 117.28 119.3 115.9 117 ...
    113.65 114.31667 113 113.23333 108.32 ...
    110.35 106.45 104.06667 106.71667 102.73333 ...
    91  108.95 103.73333 101.75 106.26667 87.68333];
lon=[39.91667 39.13333 38.03333 37.86667 40.8 ...
    41.8 43.88333 45.75 34.5 32.05 ...
    30.26667 31.87 26.08333 28.68333 36.67 ...
    34.76667  30.51667 28.21667 23.16667 22.8 ...
    20.01667 29.56667 30.66667 26.56667 25.05 ...
    29.6  34.26667  36.03333 36.56667 38.46667 43.76667];
T=[35.2 22.5 193.7 63.8 181.9...
    143.7 120.0 325.0 126.3 552.2 ...
    203.0 293.1 202.5 239.7 222.5 ...
    224.6 288.0 325.2 469.0 301.6 ...
    44.4 86.4 230.3 101.4 147.5 ...
    35.2 83.4 121.8 30.8 72.4 535.1];
surfm(lat,lon,T);
请教各位,为何上述代码运行时报错  谢谢各位
Error using gcm (line 26)
Not a map axes.

Error in nextmap (line 25)
    gcm; % Will error if gca is not a valid map axes

Error in nextmap (line 44)
        nextmap

Error in surfm (line 77)
nextmap(varargin);

Aires

你前面用的是M_Map的函数画图，后面又用matlab自带的函数surfm画，你没有建立新的坐标轴

又是那隻貓

用m_map包中自带的画等值线的命令试试