x=[114 104 108 113 119 116 123 105 88 89 102 124 135 136 127 145 141 115 90 77 79 96 122 147 155 143 112 81 67 65 75 97 111 152 163 153 123 141 79 60 55 67 88 287 290 306 324 339 341 333 293 272 293 293 265 240 221 230 251 275 84 131 160 167 173 164 140 102 74 54 44 45 59 58 101 150 172 194 222 253 304 322 320 295 263 231 192 178 177 182 182 172 154 132 102 63 39 29 33 308 345 357 353 331 303 278 264 284 45 27 19 12 9 11 17 25 40 52 78 118 141 164 186 207 228 260 281 311 329 347 344 335 326 313 285 253 202 191 194 200 175 152 110 68 231 255 285 314 347 371 375 361 324 278 253 232];
y=[266 259 243 241 251 263 294 299 275 239 214 215 237 267 313 269 221 192 205 234 280 311 334 299 244 203 179 196 228 275 306 332 354 324 257 195 165 161 178 218 278 315 341 202 226 234 229 205 180 167 169 194 421 390 373 376 392 425 437 435 360 375 342 303 241 196 161 143 163 194 237 280 330 358 394 403 409 434 458 468 451 423 387 357 351 352 360 335 293 251 222 188 158 135 124 146 191 237 304 256 216 176 147 143 162 187 222 247 132 165 185 210 251 294 320 342 375 392 414 435 444 459 474 487 495 500 499 486 472 440 408 379 365 347 332 330 328 307 261 220 173 136 98 108 223 253 372 272 232 172 129 112 126 167 190 219];
z=[35 35 35 35 35 35 30 30 30 30 30 30 30 30 25 25 25 25 25 25 25 25 20 20 20 20 20 20 20 20 20 20 15 15 15 15 15 15 15 15 15 15 15 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
xi=[0:0.1:512];
yi=[0:0.1:512];
ZI=griddata(x,y,z,xi,yi,'cubic');
mesh(xi,yi,ZI)
[XI,YI]=meshgrid(xi,yi);
[Fx,Fy]=gradient(ZI,0.001,0.001);
S=sqrt(1+Fx.^2+Fy.^2)*0.000001.*( ~isnan(ZI) ) ;
sum(S(~isnan(S)))
z =
Columns 1 through 17
35 35 35 35 35 35 30 30 30 30 30 30 30 30 25 25 25
Columns 18 through 34
25 25 25 25 25 20 20 20 20 20 20 20 20 20 20 15 15
Columns 35 through 51
15 15 15 15 15 15 15 15 15 10 10 10 10 10 10 10 10
Columns 52 through 68
10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
Columns 69 through 85
10 10 10 10 10 5 5 5 5 5 5 5 5 5 5 5 5
Columns 86 through 102
5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
Columns 103 through 119
5 5 5 5 5 5 0 0 0 0 0 0 0 0 0 0 0
Columns 120 through 136
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Columns 137 through 153
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Columns 154 through 156
0 0 0
??? Error using ==> set
View [azimuth elevation] must be a two element vector.
Error in ==> mesh at 128
set(cax,'View',y);
Warning: Z must not be a scalar or vector, not rendering surface.
Warning: Z must not be a scalar or vector, not rendering surface.
Warning: Z must not be a scalar or vector, not rendering surface.