clc,clear,close all
warning off
A=[99,110,0;433,111,0;254,196,0;330,256,0;408,69,0;501,257,0;412,345,0;112,354,0;223,232,0;276,306,0;206,359,0;122,53,0;171,25,0;204,14,0;319,16,0;354,29,0;393,54,0;435,114,0;463,171,0;489,209,0;500,241,0;492,304,0;458,341,0;405,343,0;350,317,0;334,228,0;331,204,0;286,191,0;250,197,0;219,201,0;191,186,0;163,168,0;126,144,0;103,118,0;130,321,0;166,281,0;194,249,0;246,246,0;271,283,0;264,326,0;233,347,0;179,370,0;131,375,0;104,76,0;115,60,0;126,50,0;143,38,0;152,33,0;167,26,0;178,22,0;199,15,0;216,12,0;229,10,0;234,10,0;239,10,0;243,10,0;245,9,0;250,9,0;257,9,0;263,9,0;269,10,0;277,10,0;285,10,0;300,13,0;314,15,0;328,19,0;343,25,0;358,32,0;375,41,0;390,51,0;398,58,0;404,65,0;411,77,0;418,81,0;425,91,0;432,103,0;435,115,0;439,126,0;445,139,0;449,150,0;456,161,0;464,172,0;475,187,0;483,198,0;489,210,0;493,216,0;496,228,0;500,239,0;501,247,0;501,253,0;501,260,0;501,270,0;500,282,0;496,292,0;494,299,0;489,308,0;484,319,0;475,328,0;466,336,0;458,341,0;447,345,0;436,347,0;419,346,0;405,343,0;394,340,0;383,336,0;374,332,0;351,318,0;344,311,0;336,299,0;333,290,0;331,284,0;330,279,0;330,273,0;330,263,0;330,253,0;331,247,0;331,243,0;333,232,0;334,227,0;334,220,0;334,214,0;333,210,0;331,206,0;329,201,0;326,198,0;321,195,0;317,193,0;313,192,0;306,191,0;298,191,0;290,191,0;282,191,0;277,192,0;273,193,0;266,193,0;259,195,0;252,196,0;248,198,0;240,199,0;234,201,0;227,201,0;219,201,0;214,198,0;206,194,0;195,189,0;181,180,0;174,175,0;168,173,0;158,166,0;151,161,0;142,156,0;135,151,0;126,144,0;120,138,0;108,126,0;102,117,0;99,109,0;98,102,0;98,95,0;99,91,0;101,82,0;107,71,0;114,62,0;122,54,0;131,46,0;234,236,0;241,242,0;249,249,0;254,256,0;262,268,0;269,278,0;274,290,0;275,294,0;276,297,0;277,301,0;276,306,0;275,312,0;271,317,0;265,325,0;257,332,0;248,338,0;237,345,0;228,350,0;219,355,0;209,358,0;200,362,0;173,371,0;164,374,0;157,375,0;150,375,0;146,376,0;140,376,0;135,376,0;133,375,0;125,374,0;119,370,0;114,365,0;113,360,0;112,351,0;116,340,0;125,326,0;140,308,0;153,294,0;167,280,0;182,263,0;189,255,0;201,244,0;206,239,0;210,237,0;214,234,0;218,233,0;220,232,0;224,232,0;227,232,0;230,234,0;
126,111,50;254,28,50;400,111,50;254,182,50;352,256,50;408,166,50;468,258,50;411,323,50;143,339,50;218,263,50;256,307,50;206,350,50;133,80,50;149,63,50;171,47,50;188,41,50;203,36,50;216,32,50;240,29,50;287,31,50;310,36,50;328,42,50;344,48,50;354,54,50;362,60,50;371,68,50;382,76,50;390,90,50;397,105,50;401,118,50;403,129,50;403,148,50;405,163,50;416,179,50;426,187,50;438,196,50;449,209,50;456,220,50;463,234,50;468,249,50;467,275,50;461,290,50;457,300,50;451,307,50;439,317,50;452,212,50;462,228,50;468,249,50;467,275,50;458,297,50;446,311,50;432,319,50;420,323,50;398,322,50;387,319,50;379,314,50;370,308,50;361,299,50;356,289,50;354,277,50;352,268,50;351,251,50;353,237,50;356,223,50;361,210,50;362,197,50;362,193,50;361,189,50;357,186,50;353,183,50;348,181,50;340,179,50;333,177,50;324,177,50;314,177,50;300,177,50;293,178,50;286,178,50;280,179,50;275,180,50;265,181,50;253,182,50;242,183,50;229,183,50;219,181,50;208,178,50;198,174,50;191,172,50;184,167,50;175,162,50;166,157,50;160,154,50;156,151,50;149,144,50;142,138,50;134,130,50;132,125,50;128,118,50;127,113,50;126,109,50;126,99,50;128,89,50;132,82,50;135,77,50;222,263,50;227,264,50;235,268,50;241,274,50;246,281,50;251,287,50;254,293,50;255,299,50;256,307,50;255,311,50;254,315,50;253,318,50;248,323,50;244,328,50;240,332,50;223,343,50;213,347,50;206,350,50;195,354,50;188,356,50;182,357,50;174,357,50;171,358,50;165,358,50;161,357,50;156,357,50;152,355,50;146,350,50;143,345,50;143,340,50;143,334,50;147,323,50;150,319,50;153,314,50;158,308,50;164,302,50;170,295,50;174,290,50;181,283,50;187,278,50;194,272,50;200,268,50;204,266,50;207,265,50;211,264,50;213,263,50;218,263,50;
145,111,100;254,43,100;375,111,100;254,173,100;373,256,100;402,219,100;439,258,100;411,297,100;168,330,100;212,286,100;236,308,100;205,338,100;146,94,100;152,84,100;158,77,100;166,69,100;178,62,100;191,55,100;205,50,100;214,48,100;227,46,100;235,44,100;243,43,100;250,43,100;250,43,100;258,43,100;269,43,100;282,45,100;295,47,100;313,51,100;325,57,100;339,65,100;350,72,100;362,83,100;370,97,100;374,105,100;375,110,100;376,116,100;376,124,100;376,130,100;373,139,100;368,145,100;363,151,100;355,155,100;347,159,100;338,161,100;329,163,100;322,164,100;318,165,100;313,165,100;309,167,100;298,167,100;290,168,100;284,169,100;277,170,100;274,171,100;268,171,100;262,172,100;255,172,100;251,173,100;244,173,100;237,173,100;231,172,100;224,171,100;216,170,100;206,167,100;198,164,100;188,158,100;179,154,100;168,145,100;162,139,100;153,130,100;147,121,100;145,114,100;143,106,100;143,102,100;147,93,100;373,256,100;373,250,100;377,238,100;382,230,100;388,224,100;395,220,100;399,219,100;406,219,100;414,221,100;425,228,100;431,234,100;436,245,100;439,251,100;439,258,100;439,265,100;436,272,100;432,283,100;424,290,100;414,296,100;409,297,100;402,297,100;392,294,100;384,288,100;378,279,100;375,272,100;374,265,100;373,262,100;219,287,100;227,291,100;230,294,100;232,297,100;235,300,100;236,304,100;236,308,100;236,314,100;234,317,100;231,322,100;228,326,100;224,329,100;221,329,100;215,335,100;207,338,100;195,341,100;187,342,100;185,341,100;177,340,100;173,339,100;170,337,100;169,334,100;167,331,100;167,327,100;167,324,100;170,319,100;172,314,100;177,308,100;186,299,100;193,293,100;200,289,100;204,287,100;207,286,100;211,286,100;
162,111,150;254,54,150;356,111,150;254,163,150;162,107,150;164,98,150;168,91,150;174,82,150;183,75,150;193,68,150;203,64,150;213,60,150;227,57,150;239,55,150;243,54,150;251,54,150;258,54,150;269,54,150;283,56,150;295,59,150;306,63,150;318,67,150;330,75,150;338,82,150;348,91,150;353,100,150;356,110,150;356,119,150;354,127,150;349,137,150;340,144,150;331,149,150;316,154,150;308,156,150;299,157,150;292,159,150;283,160,150;277,162,150;270,162,150;263,163,150;253,163,150;244,163,150;238,164,150;230,163,150;223,162,150;214,160,150;200,154,150;194,151,150;187,147,150;182,143,150;178,139,150;171,131,150;165,115,150;162,111,150;161,109,150;
178,112,200;254,65,200;337,111,200;254,155,200;177,105,200;181,95,200;185,90,200;191,84,200;195,83,200;200,79,200;208,74,200;219,70,200;227,67,200;233,66,200;240,65,200;250,65,200;260,65,200;273,65,200;280,66,200;290,69,200;305,75,200;314,79,200;322,84,200;330,94,200;334,100,200;337,109,200;338,114,200;338,117,200;336,124,200;329,133,200;318,140,200;306,146,200;293,150,200;288,151,200;272,154,200;266,155,200;246,155,200;237,155,200;230,154,200;218,150,200;208,147,200;198,141,200;190,136,200;184,129,200;179,118,200;178,111,200;177,109,200;
193,112,250;254,75,250;320,111,250;254,146,250;193,108,250;194,104,250;197,99,250;203,92,250;209,88,250;214,84,250;225,80,250;234,77,250;240,76,250;245,76,250;250,75,250;256,75,250;263,75,250;269,76,250;277,77,250;292,82,250;304,88,250;310,94,250;319,104,250;320,108,250;320,117,250;319,120,250;315,126,250;307,133,250;297,138,250;287,141,250;275,144,250;267,145,250;259,146,250;251,146,250;236,145,250;231,144,250;219,141,250;208,135,250;198,126,250;196,121,250;194,117,250;193,115,250;193,111,250;
211,112,300;254,86,300;300,111,300;254,137,300;211,111,300;214,104,300;219,98,300;227,94,300;236,89,300;244,87,300;250,86,300;257,86,300;263,86,300;268,87,300;281,91,300;291,98,300;297,102,300;299,108,300;300,112,300;298,117,300;293,125,300;269,135,300;262,136,300;256,137,300;216,122,300;212,117,300;211,113,300;294,124,300;289,128,300;285,130,300;280,132,300;268,136,300;264,136,300;249,137,300;245,136,300;240,136,300;235,134,300;232,134,300;229,132,300;225,130,300;221,127,300;217,124,300;215,121,300;212,118,300;211,114,300;
242,113,350;254,104,350;268,112,350;254,119,350;242,112,350;245,107,350;248,106,350;252,104,350;256,104,350;262,105,350;265,107,350;267,109,350;268,111,350;268,114,350;264,117,350;261,119,350;257,120,350;254,119,350;250,119,350;246,118,350;243,105,350;242,113,350;248,105,350];

x=A(:,1);y=A(:,2);z=A(:,3);
scatter(x,y,5,z)%散点图

%%
figure(‘color’,[1,1,1]),
[X,Y,Z]=griddata(x,y,z,linspace(1,512)’,linspace(1,512),‘v4’);%插值
n1=size(Y);
for i=1:n1(1,1)
for j=1:n1(1,2)
Y(i,j)=450-Y(i,j);
end
end
n=size(Z);
for i=1:n(1,1)
for j=1:n(1,2)
if Z(i,j)<=0
Z(i,j)=0;
end
end
end
pcolor(X,Y,Z);
shading interp%伪彩色图
%%
figure,
contourf(X,Y,Z) %等高线图
figure(‘color’,[1,1,1]),
surf(X,Y,Z)%三维曲面