fields Loading required package spam Package spam is

Σημείωση: Όταν καλούμε τη βιβλιοθήκη fields εμφανίζονται οι συνοδοί βιβλιοθήκες Loading required package: spam

(Διευρυμένη) Περιοχή μελέτης Γεωγραφικό μήκος: -30ο (δυτικό), +45ο (ανατολικό) Γεωγραφικό πλάτος: +20ο (βόρειο), +80ο

Εντολές σχεδίασης και τίτλος (με βάση την περίληψη): > as. image(max, loc, na. rm=T)->look

Αποτελέσματα παρεμβολής και διασταυρωτικής επικύρωσης : > summary(out) CALL: Krig(x = loc, Y =

Αποτελέσματα παρεμβολής και διασταυρωτικής επικύρωσης: Residual Summary: min 1 st Q median 3 rd

Χωρική κατανομή από παρεμβολή: >par(mfrow=c(1, 1)) > surface. Krig(out, xlab="lon", ylab="lat", xlim=c(30, 45), ylim=c(20,

Εντολές σχεδίασης χωρικής κατανομής από παρεμβολή χωρίς ισοπληθείς: > predict. surface(nx=272, ny=272, out)->look 1

Παρεμβολή Kriging με Arcmap Βαλκάνια Ελλάδα 2003 -2006 mean August NO 2 tropospheric VCD-GR

Slides: 28

Download presentation

Σημείωση: Όταν καλούμε τη βιβλιοθήκη fields εμφανίζονται οι συνοδοί βιβλιοθήκες Loading required package: spam Package 'spam' is loaded. Spam version 0. 23 -0 (2010 -09 -01). Type demo( spam) for some demos, help( Spam) for an overview of this package. Help for individual functions is optained by adding the suffix '. spam' to the function name, e. g. 'help(chol. spam)'. Attaching package: 'spam' The following object(s) are masked from 'package: base': backsolve, forwardsolve, norm Use help(fields) for an overview of this library( fields, keep. source=TRUE) retains comments in the source code.

(Διευρυμένη) Περιοχή μελέτης Γεωγραφικό μήκος: -30ο (δυτικό), +45ο (ανατολικό) Γεωγραφικό πλάτος: +20ο (βόρειο), +80ο (βόρειο) Μελετώμενο μέγεθος: Μέγιστο μέσου εποχικού κύκλου ύψους υετού (mm) 1920 -2016 > summary(Rain) Lat Lon Alt December Min. : 29. 55 Min. : -20. 283 Min. : -3. 00 Min. : 0. 0 1 st Qu. : 44. 76 1 st Qu. : 5. 963 1 st Qu. : 17. 75 1 st Qu. : 562. 6 Median : 51. 17 Median : 13. 183 Median : 90. 00 Median : 617. 8 Mean : 51. 41 Mean : 13. 520 Mean : 188. 00 Mean : 630. 6 3 rd Qu. : 59. 35 3 rd Qu. : 23. 771 3 rd Qu. : 207. 50 3 rd Qu. : 716. 3 Max. : 78. 25 Max. : 39. 250 Max. : 3106. 00 Max. : 1529. 0 January February March April Min. : 167. 3 Min. : 189. 7 Min. : 147. 0 Min. : 279. 3 1 st Qu. : 419. 2 1 st Qu. : 574. 2 1 st Qu. : 383. 6 1 st Qu. : 437. 1 Median : 518. 7 Median : 618. 5 Median : 484. 1 Median : 501. 3 Mean : 622. 2 Mean : 675. 4 Mean : 588. 8 Mean : 577. 8 3 rd Qu. : 767. 6 3 rd Qu. : 754. 8 3 rd Qu. : 702. 5 3 rd Qu. : 693. 0 Max. : 1898. 8 Max. : 1627. 2 Max. : 1994. 5 Max. : 1428. 8 May June July August Min. : 0. 0 Min. : 344. 6 Min. : 173. 8 Min. : 222. 1 1 st Qu. : 610. 8 1 st Qu. : 505. 5 1 st Qu. : 369. 6 1 st Qu. : 350. 7 Median : 663. 6 Median : 609. 9 Median : 466. 6 Median : 440. 9 Mean : 676. 7 Mean : 711. 0 Mean : 566. 1 Mean : 535. 8 3 rd Qu. : 755. 5 3 rd Qu. : 860. 4 3 rd Qu. : 656. 0 3 rd Qu. : 636. 2 Max. : 1706. 9 Max. : 1827. 2 Max. : 2013. 6 Max. : 1833. 4 September October November max Min. : 224. 7 Min. : 160. 6 Min. : 71. 69 Min. : 383. 0 1 st Qu. : 377. 3 1 st Qu. : 420. 6 1 st Qu. : 482. 06 1 st Qu. : 638. 6 Median : 468. 8 Median : 468. 4 Median : 548. 29 Median : 689. 8 Mean : 533. 2 Mean : 538. 0 Mean : 588. 58 Mean : 809. 7 3 rd Qu. : 616. 2 3 rd Qu. : 645. 8 3 rd Qu. : 678. 14 3 rd Qu. : 931. 2 Max. : 1369. 4 Max. : 1368. 5 Max. : 1495. 26 Max. : 2013. 6 min Amp Monthofmax Monthofmin Min. : 0. 0 Min. : 112. 7 Min. : 1. 000 Min. : 3. 000 1 st Qu. : 322. 2 1 st Qu. : 240. 9 1 st Qu. : 5. 000 1 st Qu. : 8. 000 Median : 410. 2 Median : 326. 8 Median : 5. 000 Median : 8. 000 Mean : 434. 8 Mean : 374. 9 Mean : 5. 743 Mean : 8. 562 3 rd Qu. : 512. 6 3 rd Qu. : 414. 3 3 rd Qu. : 6. 000 3 rd Qu. : 10. 000 Max. : 1280. 2 Max. : 1389. 0 Max. : 12. 000

Εντολές σχεδίασης και τίτλος (με βάση την περίληψη): > as. image(max, loc, na. rm=T)->look > image. plot(look, xlab="lon", ylab="lat", xlim=c(-30, 45), ylim=c(20, 80), zlim=c(350, 2100), + main="1920 -2016 rainfall (mm) + mean seasonal cycle max")

Αποτελέσματα παρεμβολής και διασταυρωτικής επικύρωσης : > summary(out) CALL: Krig(x = loc, Y = max) Number of Observations: 272 Number of unique points: 267 Number of parameters in the null space 3 Parameters for fixed spatial drift 3 Effective degrees of freedom: 76. 8 Residual degrees of freedom: 195. 2 MLE sigma 191 GCV sigma 188. 3 Pure error sigma 20. 04 MLE rho 17010 Scale passed for covariance (rho) <NA> Scale passed for nugget (sigma^2) <NA> Smoothing parameter lambda 2. 145

Αποτελέσματα παρεμβολής και διασταυρωτικής επικύρωσης: Residual Summary: min 1 st Q median 3 rd Q max -442. 30 -77. 96 -23. 81 32. 55 695. 50 Covariance Model: stationary. cov Covariance function is Names of non-default covariance arguments: DETAILS ON SMOOTHING PARAMETER: Method used: GCV Cost: 1 lambda tr. A GCV. one GCV. model shat 2. 145 76. 844 51858. 293 49405. 764 51468. 206 188. 276 Summary of all estimates found for lambda tr. A GCV shat REML GCV 2. 145 76. 84 51858 188. 3 1810 GCV. model 1. 868 83. 61 51440 184. 8 NA GCV. one 1. 410 98. 25 49186 177. 3 NA RMSE NA NA NA pure error NA NA NA REML 1. 365 99. 99 52225 176. 4 1810

Χωρική κατανομή από παρεμβολή: >par(mfrow=c(1, 1)) > surface. Krig(out, xlab="lon", ylab="lat", xlim=c(30, 45), ylim=c(20, 80), zlim=c(350, 2100), + main="1920 -2016 rainfall (mm) + mean seasonal cycle max by kriging")

Εντολές σχεδίασης χωρικής κατανομής από παρεμβολή χωρίς ισοπληθείς: > predict. surface(nx=272, ny=272, out)->look 1 > image. plot(look 1, xlab="lon", ylab="lat", xlim=c(-30, 45), ylim=c(20, 80), zlim=c(350, 2100), + main="1920 -2016 rainfall (mm) + mean seasonal cycle max by kriging") > addland() > grid() # εμφάνιση θέσης σταθμών μετρήσεων >points(loc)

Παρεμβολή Kriging με Arcmap Βαλκάνια Ελλάδα 2003 -2006 mean August NO 2 tropospheric VCD-GR 1997 -2003 mean August NO 2 tropospheric VCD-GR 2003 -2006 mean December NO 2 tropospheric VCD-GR 1997 -2003 mean December NO 2 tropospheric VCD-GR