博文
Handbook of Broadleaf.Korean.pine.LAI——基于分层空间拓扑的叶面积指数空间分布模拟
|
Handbook of Broadleaf.Korean.pine.LAI
软件包简介
叶面积指数（LAI）指单位水平地面面积上总叶面积的一半，为植物冠层结构的最常用参数之一。我团队提出了一种基于分层泰森多边形估算复层林叶面积指数的方法，相较于传统基于公式的叶面积指数估计方法，该法可通过泰森多边形估算林分中各林木的树冠下叶面积指数，并在此基础上较为准确的计算林分中任意位点处的若干与叶面积指数相关的冠层结构参数（局部叶面积指数、局部叶面积指数标注差、各类树种对局部叶面积指数的贡献率，叶面积指数垂直结构特征等）。
我团队将基于分层泰森多边形估算复层林叶面积指数的方法用于阔叶红松林叶面积指数的估算，开发了R包Broadleaf.Korean.pine.LAI。
依赖关系
本包依赖于 ape；sp ；sf；raster；grDevices；rgeos；gstat ；tcltk ；ggplot2；plot3D ；deldir ；grDevices；stats 运行前请先安装
软件包安装方式
library(devtools);install_github(“DuXinChina/Broadleaf.Korean.pine.LAI/Broadleaf.Korean.pine.LAI”)
1 示例数据展示
示例数据为虚拟数据，数据为50m×50m样地，样地中包含色木、紫椴、冷杉、红松四种树种的空间坐标、胸径、树高信息。依据树高，将阔叶红松林分为不同林层。其中，下木层树高小于10m；亚林层树高为10m—16m；主林层树高为16m—30m；超冠层树高大于30m。虚拟数据中，超冠层包含两株红松，为随机分布；主林层中包含15株红松、10株冷杉、10株紫椴、5株色木，均为随机分布；亚林层中包含50株色木，为以（10,40）和（40,10）为中心的各25株的聚集分布；林下层包含四类树种各20株植株，各树种在小尺度为聚集分布，在大尺度为随机分布。超冠层中，林木胸径大于60cm；主林层中，林木胸径为30cm—60 cm；亚林层中林木胸径为10 cm-30cm；下木层胸径小于10 cm。
数据中x与y列分别为林木横纵坐标、DBH为胸径、H为树高、Species列为树种种名。
Species列可包含的常见树种种名为HS 红松、YS 云杉、LS 冷杉、ZD 紫椴、KD 糠椴、MGL 蒙古栎、SQL 水曲柳、HTQ 胡桃楸、HBL 黄檗、SMQ 色木槭、QKQ 青楷槭。HKQ 花楷槭、JSQ 假色槭、NJQ 拧筋槭、BNQ 白牛槭、CY 春榆、BH 白桦、HH 怀槐、FH 枫桦、LYY 裂叶榆、QT 其他树种。
data = Broadleaf.Korean.pine.LAI::b data
## x y DBH H Species ## 1 33.4288480 18.9819806 65.580755 35 HS ## 2 6.4346078 42.0391756 66.180769 35 HS ## 3 18.1460637 19.2660497 52.264410 25 HS ## 4 27.5928278 13.8968076 32.366493 25 HS ## 5 6.1081642 25.6396705 35.555799 25 HS ## 6 39.5111778 5.1855691 53.037166 25 HS ## 7 18.2991147 36.4641006 31.796712 25 HS ## 8 25.6670756 4.6233300 31.566986 25 HS ## 9 25.2746910 5.0032010 34.373929 25 HS ## 10 11.1674790 8.0279473 58.859233 25 HS ## 11 10.9212414 11.5923407 35.492684 25 HS ## 12 22.6771733 12.8692340 50.351396 25 HS ## 13 3.1528245 40.5892981 40.749832 25 HS ## 14 48.9035375 10.3960484 43.531746 25 HS ## 15 44.4565091 36.6393705 30.292931 25 HS ## 16 18.4087182 19.8480956 34.429790 25 HS ## 17 37.7821442 18.2801837 38.516582 25 HS ## 18 10.2223582 0.3512422 35.396181 25 LS ## 19 23.3121148 6.1817106 40.818709 25 LS ## 20 3.7540585 32.3250375 33.630106 25 LS ## 21 28.3364101 6.5744882 39.985535 25 LS ## 22 6.8510692 9.6053103 54.610498 25 LS ## 23 33.3752221 36.7657443 34.227696 25 LS ## 24 44.9338280 0.6791399 56.573134 25 LS ## 25 1.5318491 24.5214906 34.453612 25 LS ## 26 40.2367382 11.2009961 58.903647 25 LS ## 27 16.3113690 25.0485350 48.073396 25 LS ## 28 31.4198530 6.2785452 52.088891 25 ZD ## 29 46.0771567 27.8802443 51.481795 25 ZD ## 30 0.2149013 46.6323809 59.699024 25 ZD ## 31 48.3191063 30.0955467 52.972707 25 ZD ## 32 25.6461456 8.1469459 30.973053 25 ZD ## 33 37.0948580 34.3532822 45.913191 25 ZD ## 34 15.7041204 49.4823055 57.685524 25 ZD ## 35 9.0926424 24.8328917 45.382739 25 ZD ## 36 10.1880918 36.0375070 58.164637 25 ZD ## 37 6.6190212 46.7338556 44.880348 25 ZD ## 38 27.8663572 22.0768596 51.736290 25 SMQ ## 39 15.3366110 16.2149379 42.419761 25 SMQ ## 41 27.9685288 10.7013771 31.988869 25 SMQ ## 42 30.4606082 4.4862904 31.223996 25 SMQ ## 43 27.4182794 5.0037950 20.837063 15 SMQ ## 44 44.4396113 28.8874433 10.625898 15 SMQ ## 45 34.8701260 10.9891497 26.209086 15 SMQ ## 46 29.3663570 24.8532922 21.750956 15 SMQ ## 47 39.0940332 21.3669462 14.558203 15 SMQ ## 48 37.8461923 0.8798485 25.238661 15 SMQ ## 49 40.8974968 15.1046335 15.410212 15 SMQ ## 50 43.7170631 12.4452323 24.453914 15 SMQ ## 52 34.7943280 15.7893710 13.230423 15 SMQ ## 53 44.7095683 16.8244097 11.873488 15 SMQ ## 54 32.1441950 10.4155232 16.800018 15 SMQ ## 55 40.1892333 20.6559082 24.111454 15 SMQ ## 56 45.1569390 11.7849286 29.909211 15 SMQ ## 57 31.1358896 4.4016615 25.619205 15 SMQ ## 58 41.9316033 12.6385274 19.933430 15 SMQ ## 59 38.1089669 6.2206536 27.308530 15 SMQ ## 60 44.1061557 15.5097960 13.451536 15 SMQ ## 61 43.9375359 4.8815030 15.789073 15 SMQ ## 62 43.9779725 1.6313955 16.654748 15 SMQ ## 63 42.8673806 10.6359013 23.262358 15 SMQ ## 64 40.8195774 8.6401999 29.290362 15 SMQ ## 65 35.0804810 24.0241466 22.261274 15 SMQ ## 66 23.9579239 14.2034792 12.860445 15 SMQ ## 67 45.5492528 13.8757384 28.784073 15 SMQ ## 68 27.7498276 42.5004315 27.557451 15 SMQ ## 69 16.9513238 39.2511895 11.311346 15 SMQ ## 70 7.8428609 35.1662446 18.100418 15 SMQ ## 71 21.5821916 46.1810341 15.722701 15 SMQ ## 72 9.8932887 38.4724617 18.591866 15 SMQ ## 73 21.4302064 47.5721941 24.531588 15 SMQ ## 74 16.0871787 38.7437988 17.157289 15 SMQ ## 75 27.3022366 30.3321379 21.895134 15 SMQ ## 76 11.6537314 37.9788635 17.060371 15 SMQ ## 77 8.9436178 40.0632726 17.962532 15 SMQ ## 78 6.2872774 39.4914437 28.778782 15 SMQ ## 79 5.3213029 47.4720506 26.740215 15 SMQ ## 80 4.4828161 40.1074313 23.338306 15 SMQ ## 81 17.4019105 37.3673022 22.306889 15 SMQ ## 82 14.4704528 41.0565769 20.476022 15 SMQ ## 83 15.4274825 48.2122797 19.459245 15 SMQ ## 84 10.7972558 43.9528950 24.105329 15 SMQ ## 85 1.3316003 48.4359719 27.001667 15 SMQ ## 86 19.4026872 38.7502684 24.929608 15 SMQ ## 87 20.8525298 43.6054503 12.541657 15 SMQ ## 88 10.3695183 42.6703659 20.995959 15 SMQ ## 89 11.9001001 23.5515179 25.591832 15 SMQ ## 90 7.8739939 39.0204771 16.820479 15 SMQ ## 91 12.9090867 40.7360448 15.755752 15 SMQ ## 92 14.8874215 28.9533239 19.382401 15 SMQ ## 93 39.6436424 40.9811422 9.934820 8 LS ## 94 49.6122418 34.0210738 8.551500 8 LS ## 95 46.5867969 37.3322312 5.193316 8 LS ## 96 49.0000000 48.3085792 8.632041 8 LS ## 97 29.0852543 41.1179223 6.165639 8 LS ## 98 29.0206565 34.8073513 9.694307 8 LS ## 99 28.7770563 37.6178602 9.306080 8 LS ## 100 38.6189406 33.9430026 8.822112 8 LS ## 101 45.9365553 22.4422532 6.603459 8 LS ## 102 36.6124557 23.8282281 6.482322 8 LS ## 103 42.8625806 7.1329634 8.337840 8 LS ## 104 32.6357700 21.9731506 7.542603 8 LS ## 105 11.7305549 2.2269401 6.781676 8 LS ## 106 5.1875622 1.0000000 8.702628 8 LS ## 107 8.4335608 0.1246017 9.224902 8 LS ## 108 6.2689242 3.7764784 5.028092 8 LS ## 109 49.7399440 13.3581418 5.351386 8 LS ## 110 44.1681187 13.2980574 9.975786 8 LS ## 111 49.0000000 21.0523839 7.792711 8 LS ## 112 47.4105839 1.0000000 8.252535 8 LS ## 113 49.0000000 1.0000000 9.106646 8 HS ## 114 48.3429683 5.6185863 7.605254 8 HS ## 115 40.5465709 8.1345473 8.666798 8 HS ## 116 49.0000000 11.0004526 9.061282 8 HS ## 117 6.2794533 4.2384272 5.919919 8 HS ## 118 1.0000000 5.9588289 7.419849 8 HS ## 119 3.8418102 12.5395119 6.657954 8 HS ## 120 1.7257937 7.9206185 8.538493 8 HS ## 121 41.8513062 16.1607322 7.772832 8 HS ## 122 40.2428301 19.8903321 5.406923 8 HS ## 123 36.5047511 15.0316257 9.628256 8 HS ## 124 36.7982127 11.3218356 9.553490 8 HS ## 125 31.1297617 16.2011382 8.440440 8 HS ## 126 36.6997820 25.2196392 7.903824 8 HS ## 127 32.7186422 11.4235729 9.530176 8 HS ## 128 28.2869781 20.5159761 6.201615 8 HS ## 129 17.7794649 14.8297948 8.094559 8 HS ## 130 13.0119315 15.5301695 8.087604 8 HS ## 131 12.8823991 15.8329528 8.933832 8 HS ## 132 10.1283477 10.6649023 5.769270 8 HS ## 133 3.0805107 29.2901996 8.609461 8 ZD ## 134 0.9071648 25.2446763 8.010450 8 ZD ## 135 1.0000000 29.8599316 5.312565 8 ZD ## 136 1.0000000 31.6702470 8.310919 8 ZD ## 137 18.3311714 49.0000000 7.540457 8 ZD ## 138 21.4111116 46.7843710 5.504480 8 ZD ## 139 18.8694693 43.5405025 8.373520 8 ZD ## 140 19.7880615 46.1762995 5.829197 8 ZD ## 141 7.4560338 44.4200915 8.199286 8 ZD ## 142 16.6670340 47.2940320 6.212923 8 ZD ## 143 12.5297550 49.0000000 8.027718 8 ZD ## 144 15.8983970 49.0000000 7.926028 8 ZD ## 145 38.9231748 41.7665160 6.291640 8 ZD ## 146 37.4657599 47.2459319 8.682012 8 ZD ## 147 49.6668905 39.7330850 8.203196 8 ZD ## 148 31.8817342 39.6527677 5.795613 8 ZD ## 149 5.6448613 28.8535785 5.703050 8 ZD ## 150 4.2739190 17.6424662 6.565414 8 ZD ## 151 10.9522913 19.9384604 9.095775 8 ZD ## 152 9.6877684 14.2428889 6.697895 8 ZD ## 153 7.8820318 49.0000000 9.916860 8 SMQ ## 154 1.0000000 44.4448761 6.899585 8 SMQ ## 155 2.3695873 48.7241619 6.060310 8 SMQ ## 156 5.6845312 42.1728318 7.417518 8 SMQ ## 157 20.6701393 33.8747094 8.585438 8 SMQ ## 158 17.3413093 28.6130006 7.200311 8 SMQ ## 159 24.7963305 26.5485412 6.546525 8 SMQ ## 160 25.8814922 23.0839561 7.570419 8 SMQ ## 161 38.1404212 18.1701050 7.064143 8 SMQ ## 162 44.4997075 17.7732229 8.575503 8 SMQ ## 163 41.8010491 30.4660159 8.362883 8 SMQ ## 164 43.8052430 18.7144681 7.256104 8 SMQ ## 165 27.9742202 48.2190508 7.243758 8 SMQ ## 166 27.1037613 49.0000000 8.991668 8 SMQ ## 167 33.9898945 48.8409792 8.534445 8 SMQ ## 168 37.3406029 49.0000000 6.711967 8 SMQ ## 169 10.4386867 34.3300143 5.309605 8 SMQ ## 170 8.0524403 27.1832789 9.161863 8 SMQ ## 171 1.0000000 25.2146483 9.680147 8 SMQ ## 172 12.1497107 29.6801073 6.772546 8 SMQ
2 示例
2.1 LAI.single(a,b,r)
2.1.1功能介绍： 基于传统方法，计算阔叶红松林内以任意位点为中心，一定半径的圆形区域内的叶面积指数，并给出不同树种在叶面积指数中的占比，及针叶树种与阔叶树种叶面积指数的占比。
2.1.2参数意义： a—需要计算叶面积指数的位点坐标 b—样地中林木的坐标、树种、胸径、树高信息 r—以位点中心圆的半径
2.1.3例：
计算以（25,25）为中心5m为半径的圆内的叶面积
a=matrix(c(25,25),1,2)
colnames(a)=c("x","y")
a## x y ## [1,] 25 25
b=Broadleaf.Korean.pine.LAI::b[,-4] head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
Broadleaf.Korean.pine.LAI::LAI.single(a=a,b=b,r=5)
## $LAI ## [1] 12.23283 ## ## $Species_LAI ## HS_LAI LS_LAI YS_LAI ZD_LAI KD_LAI MGL_LAI SQL_LAI HTQ_LAI HBL_LAI ## [1,] 0 0 0 0 0 0 0 0 0 ## SMQ_LAI QKQ_LAI HKQ_LAI JSQ_LAI NJQ_LAI BNQ_LAI CY_LAI BH_LAI HH_LAI ## [1,] 12.23283 0 0 0 0 0 0 0 0 ## FH_LAI LYY_LAI QT_LAI ## [1,] 0 0 0 ## ## $N_B_LAI ## Needles_LAI Broadleaf_LAI N_L_percent B_L_percent ## [1,] 0 12.23283 0 100
总叶面积指数为12.23283、树种色木槭叶面积为12.23283、阔叶树种叶面积12.23283、阔叶占比100%
2.2 LAI.mult(a,b,r)
2.2.1功能介绍： 计算阔叶红松林内的多个位点为中心，一定半径的圆形区域内的叶面积指数，针叶树种与阔叶树种叶面积指数以及针阔叶面积分别的占比。
2.2.2参数意义：
a—需要计算叶面积指数的多个位点坐标
b—样地中林木的坐标、树种、胸径、树高信息
r—以位点中心圆的半径
2.2.3例：
计算多个样点的叶面积指数，针叶叶面积指数，阔叶叶面积指数以及针阔叶树种叶面积指数的占比
a=matrix(runif(10,5,45),5,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 22.620917 23.239634 ## 2 10.191897 9.485000 ## 3 11.337259 6.751937 ## 4 9.478818 34.610415 ## 5 40.675446 40.486808
b=Broadleaf.Korean.pine.LAI::b[,-4] head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
Broadleaf.Korean.pine.LAI::LAI.mult(a=a,b=b,r=5)
## $LAI ## x y LAI ## 1 22.620917 23.239634 0.6615122 ## 2 10.191897 9.485000 18.1454956 ## 3 11.337259 6.751937 12.4831165 ## 4 9.478818 34.610415 13.6670219 ## 5 40.675446 40.486808 0.3555538 ## ## $B_N_LAI ## x y Needles_LAI Broadleaf_LAI N_L_percent B_L_percent ## 1 22.620917 23.239634 0.0000000 0.6615122 0.00000 100.000000 ## 2 10.191897 9.485000 17.9309321 0.2145634 98.81754 1.182461 ## 3 11.337259 6.751937 12.4831165 0.0000000 100.00000 0.000000 ## 4 9.478818 34.610415 0.0000000 13.6670219 0.00000 100.000000 ## 5 40.675446 40.486808 0.1618193 0.1937345 45.51191 54.488093
2.3 Plot.LAI.Krig (minx,maxx,miny,maxy,b,seq,r)
2.3.1功能介绍： 绘制阔叶红松林叶面积指数的克里格插值图。计算阔叶红松林叶面积指数半变异函数
2.3.2参数意义：
minx—绘制克里格插值图的最小横坐标
maxx—绘制克里格插值图的最大横坐标
miny—绘制克里格插值图的最小纵坐标
maxy—绘制克里格插值图的最大纵坐标
b—样地中的林木坐标、树种、胸径
seq—绘制克里格图的空间分辨率
r—以位点中心圆的半径
2.3.3例：
####
b=Broadleaf.Korean.pine.LAI::b[,-4] head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
p=Broadleaf.Korean.pine.LAI::Plot.LAI.Krig(5,45,5,45,b=b,seq=20,r=3)
## model psill range ## 1 Nug 1.77716 0.000000 ## 2 Sph 118.41951 6.636184
## [1] "Coefficient_of_Determination= 0.976014299370735" ## [using ordinary kriging]
p1=p+geom_vline(xintercept = c(5,45),linetype=2)+geom_hline(yintercept =c(5,45),linetype=2) p1=p+scale_x_continuous(expand= c(0, 5))+scale_y_continuous(expand= c(0, 5))
## Scale for 'x' is already present. Adding another scale for 'x', which will ## replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which will ## replace the existing scale.
p2=p1+geom_point(data=b,aes(x=x,y=y),size=b$DBH/8,col="grey4",alpha=0.2)+ scale_x_continuous(expand= c(0, 0))+scale_y_continuous(expand= c(0, 0))
## Scale for 'x' is already present. Adding another scale for 'x', which will ## replace the existing scale. ## Scale for 'y' is already present. Adding another scale for 'y', which will ## replace the existing scale.
p2
2.4 Standing_Forest_LAI(minx,maxx,miny,maxy,b)
2.4.1功能介绍：
计算阔叶红松林全林分的叶面积指数，并给出不同树种叶面积指数的量与占比，以及针、阔叶树种叶面积指数的量与占比。
2.4.2 参数意义：
minx—计算叶面积指数样地横坐标最小边界
maxx—计算叶面积指数样地横坐标最大边界
miny—计算叶面积指数样地纵坐标最小边界
maxy—计算叶面积指数样地纵坐标最大边界
b—林木分布信息
2.4.3例：
b=Broadleaf.Korean.pine.LAI::b[,-4] head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
Broadleaf.Korean.pine.LAI::Standing_Forest_LAI(minx=5,maxx=45,miny=5,maxy=45,b)
## $LAI ## [1] 11.25676 ## ## $Species_LAI ## HS_LAI LS_LAI YS_LAI ZD_LAI KD_LAI MGL_LAI SQL_LAI HTQ_LAI HBL_LAI ## [1,] 3.948259 1.242444 0 1.352644 0 0 0 0 0 ## SMQ_LAI QKQ_LAI HKQ_LAI JSQ_LAI NJQ_LAI BNQ_LAI CY_LAI BH_LAI HH_LAI ## [1,] 4.713409 0 0 0 0 0 0 0 0 ## FH_LAI LYY_LAI QT_LAI ## [1,] 0 0 0 ## ## $N_B_LAI ## Needles_LAI Broadleaf_LAI N_L_percent B_L_percent ## [1,] 5.190703 6.066053 46.11189 53.88811
林分叶面积指数11.896。红松叶面积指数3.948；冷杉叶面积指数1.242；紫椴叶面积指数1.353；色木叶面积指数5.352。针叶树叶面积指5.191；阔叶树叶面积指数6.705。针叶树种叶面积指数占比43.636%；阔叶树种叶面积指数占比56.364%。
2.5 LSD_LAI(a,b,r,Lr)
2.5.1功能介绍：
在以一定半径圆的面积为分辨率计算阔叶红松林叶面指数后。计算以某一位点为中心，一定半径圆内叶面积指数的局部标准差。通常情况下，林缘处及大树小树的冠层交错处的叶面积指数具有较大的局部标准差。
2.5.2参数意义：
a—需要计算叶面积指数局部标准差的位点坐标
b—林木分布信息及胸径、树高
r—计算林木内叶面积指数时所用圆的半径
Lr—计算位点叶面积指数局部标准差时所用圆的半径
2.5.3例：
a=data.frame(28,25)
colnames(a)=c("x","y")
a## x y ## 1 28 25
b=Broadleaf.Korean.pine.LAI::b[,-4] head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
result=Broadleaf.Korean.pine.LAI::LSD_LAI(a=a,b=b,r=3,Lr=1.5) result
## x y LSD_LAI ## 1 28 25 13.82343
2.6 LSD_LAI_mult(a,b,r,Lr)
2.6.1功能介绍：
计算多个位点叶面积指数的局部标准差
2.6.2参数意义：
a—需要计算叶面积指数局部标准差的位点坐标
b—林木分布信息及胸径、树高
r—计算林木内叶面积指数时所用圆的半径
Lr—计算位点叶面积指数局部标准差时所用圆的半径
2.6.3例：
a=matrix(runif(10,5,45),5,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 11.910299 33.28216 ## 2 38.291663 14.58496 ## 3 17.615739 33.31745 ## 4 7.594751 24.31607 ## 5 14.801852 19.40255
b=Broadleaf.Korean.pine.LAI::b[,-4] head(b)
## x y DBH Species ## 1 33.428848 18.981981 65.58075 HS ## 2 6.434608 42.039176 66.18077 HS ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS
result=Broadleaf.Korean.pine.LAI::LSD_LAI_mult(a=a,b=b,r=3,Lr=1.5) formattable::formattable(result)
| x | y | LSD_LAI |
|---|---|---|
| 11.910299 | 33.28216 | 10.155311 |
| 38.291663 | 14.58496 | 7.173425 |
| 17.615739 | 33.31745 | 4.786455 |
| 7.594751 | 24.31607 | 5.277411 |
| 14.801852 | 19.40255 | 14.379131 |
2.7 Plot.LSD_LAI.Krig(minx,maxx,miny,maxy,b,seq,r,Lr)
2.7.1功能介绍：
绘制阔叶红松林叶面积指数局部标准差的克里格插值图。计算阔叶红松林叶面积指数局部标准差的半变异函数。
2.7.2参数意义：
minx—绘制克里格插值图的最小横坐标
maxx—绘制克里格插值图的最大横坐标
miny—绘制克里格插值图的最小纵坐标
maxy—绘制克里格插值图的最大纵坐标
b—样地中的林木坐标、树种、胸径
seq—绘制克里格图的空间分辨率
r—计算林分叶面积指数时以位点中心圆的半径
Lr—计算林分叶面积指数局部标准差时取样圆的半径
2.7.3例：
p=Broadleaf.Korean.pine.LAI::Plot.LSD_LAI.Krig(minx=5,maxx=45,miny=5,maxy=45,b=b, seq=20,r=3,Lr=1)
## model psill range ## 1 Nug 4.822950 0.000000 ## 2 Sph 9.399724 3.202621
## [1] "Coefficient_of_Determination= 0.620507794377322" ## [using ordinary kriging]
p1=p+geom_vline(xintercept = c(5,45),linetype=2)+geom_hline(yintercept =c(5,45),linetype=2) p1=p+scale_x_continuous(expand= c(0, 5))+scale_y_continuous(expand= c(0, 5))
## Scale for 'x' is already present. Adding another scale for 'x', which will ## replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which will ## replace the existing scale.
p2=p1+geom_point(data=b,aes(x=x,y=y),size=b$DBH/8,col="grey4",alpha=0.2)+ scale_x_continuous(expand= c(0, 0))+scale_y_continuous(expand= c(0, 0))
## Scale for 'x' is already present. Adding another scale for 'x', which will ## replace the existing scale. ## Scale for 'y' is already present. Adding another scale for 'y', which will ## replace the existing scale.
p2
2.8 Voronoi.LAI(minx,maxx,miny,maxy,boundary,b,r)
2.8.1功能介绍：
基于泰森多边形法计算林地内（单一林层）各株树木下的叶面积指数。
2.8.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
2.8.3例：
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H<10) b=b[,-4] head(b)
## x y DBH Species ## 93 39.64364 40.98114 9.934820 LS ## 94 49.61224 34.02107 8.551500 LS ## 95 46.58680 37.33223 5.193316 LS ## 96 49.00000 48.30858 8.632041 LS ## 97 29.08525 41.11792 6.165639 LS ## 98 29.02066 34.80735 9.694307 LS
result=Broadleaf.Korean.pine.LAI::Voronoi.LAI(minx=0,maxx=50,miny=0,maxy=50,boundary=5,b=b,r=2)
## deldir 0.1-25
##
## PLEASE NOTE: The components "delsgs" and "summary" of the
## object returned by deldir() are now DATA FRAMES rather than
## matrices (as they were prior to release 0.0-18).
## See help("deldir").
##
## PLEASE NOTE: The process that deldir() uses for determining
## duplicated points has changed from that used in version
## 0.0-9 of this package (and previously). See help("deldir").formattable::formattable(result)
| x | y | DBH | Species | deldir_area | LAI | |
|---|---|---|---|---|---|---|
| 1131 | 49.0000000 | 1.0000000 | 9.106646 | HS | 6.351193 | 3.5706003 |
| 1141 | 48.3429683 | 5.6185863 | 7.605254 | HS | 10.526040 | 1.5490501 |
| 1151 | 40.5465709 | 8.1345473 | 8.666798 | HS | 8.673114 | 2.3881160 |
| 1161 | 49.0000000 | 11.0004526 | 9.061282 | HS | 9.023746 | 2.4902250 |
| 1171 | 6.2794533 | 4.2384272 | 5.919919 | HS | 6.923454 | 1.4886705 |
| 1181 | 1.0000000 | 5.9588289 | 7.419849 | HS | 11.101522 | 1.4038545 |
| 1191 | 3.8418102 | 12.5395119 | 6.657954 | HS | 9.804799 | 1.3034956 |
| 1201 | 1.7257937 | 7.9206185 | 8.538493 | HS | 9.751063 | 2.0668947 |
| 1211 | 41.8513062 | 16.1607322 | 7.772832 | HS | 10.233926 | 1.6581341 |
| 1221 | 40.2428301 | 19.8903321 | 5.406923 | HS | 10.203823 | 0.8556186 |
| 1231 | 36.5047511 | 15.0316257 | 9.628256 | HS | 9.955278 | 2.5225137 |
| 1241 | 36.7982127 | 11.3218356 | 9.553490 | HS | 11.079585 | 2.2344182 |
| 1251 | 31.1297617 | 16.2011382 | 8.440440 | HS | 9.774625 | 2.0187643 |
| 1261 | 36.6997820 | 25.2196392 | 7.903824 | HS | 8.563030 | 2.0432608 |
| 1271 | 32.7186422 | 11.4235729 | 9.530176 | HS | 11.262486 | 2.1883197 |
| 1281 | 28.2869781 | 20.5159761 | 6.201615 | HS | 9.819162 | 1.1429130 |
| 1291 | 17.7794649 | 14.8297948 | 8.094559 | HS | 10.289700 | 1.7762754 |
| 1301 | 13.0119315 | 15.5301695 | 8.087604 | HS | 8.745888 | 2.0865342 |
| 1311 | 12.8823991 | 15.8329528 | 8.933832 | HS | 6.865428 | 3.1892889 |
| 1321 | 10.1283477 | 10.6649023 | 5.769270 | HS | 10.686881 | 0.9199679 |
| 931 | 39.6436424 | 40.9811422 | 9.934820 | LS | 6.717356 | 1.8920034 |
| 941 | 49.6122418 | 34.0210738 | 8.551500 | LS | 10.466771 | 0.8900805 |
| 951 | 46.5867969 | 37.3322312 | 5.193316 | LS | 9.450425 | 0.3508748 |
| 961 | 49.0000000 | 48.3085792 | 8.632041 | LS | 9.759273 | 0.9733234 |
| 971 | 29.0852543 | 41.1179223 | 6.165639 | LS | 10.572632 | 0.4475085 |
| 981 | 29.0206565 | 34.8073513 | 9.694307 | LS | 9.591146 | 1.2595182 |
| 991 | 28.7770563 | 37.6178602 | 9.306080 | LS | 10.123716 | 1.0964011 |
| 1001 | 38.6189406 | 33.9430026 | 8.822112 | LS | 11.587956 | 0.8575519 |
| 1011 | 45.9365553 | 22.4422532 | 6.603459 | LS | 12.710238 | 0.4290834 |
| 1021 | 36.6124557 | 23.8282281 | 6.482322 | LS | 8.639266 | 0.6075241 |
| 1031 | 42.8625806 | 7.1329634 | 8.337840 | LS | 10.265588 | 0.8611864 |
| 1041 | 32.6357700 | 21.9731506 | 7.542603 | LS | 10.991263 | 0.6535290 |
| 1051 | 11.7305549 | 2.2269401 | 6.781676 | LS | 9.651604 | 0.5971063 |
| 1061 | 5.1875622 | 1.0000000 | 8.702628 | LS | 7.232443 | 1.3357205 |
| 1071 | 8.4335608 | 0.1246017 | 9.224902 | LS | 10.791388 | 1.0100682 |
| 1081 | 6.2689242 | 3.7764784 | 5.028092 | LS | 6.016010 | 0.5154787 |
| 1091 | 49.7399440 | 13.3581418 | 5.351386 | LS | 10.095259 | 0.3495082 |
| 1101 | 44.1681187 | 13.2980574 | 9.975786 | LS | 10.572373 | 1.2124097 |
| 1111 | 49.0000000 | 21.0523839 | 7.792711 | LS | 10.145719 | 0.7574857 |
| 1121 | 47.4105839 | 1.0000000 | 8.252535 | LS | 7.036175 | 1.2299673 |
| 1331 | 3.0805107 | 29.2901996 | 8.609461 | ZD | 7.397749 | 3.4315899 |
| 1341 | 0.9071648 | 25.2446763 | 8.010450 | ZD | 7.861971 | 2.8704561 |
| 1351 | 1.0000000 | 29.8599316 | 5.312565 | ZD | 7.944869 | 1.4531956 |
| 1361 | 1.0000000 | 31.6702470 | 8.310919 | ZD | 8.341543 | 2.8729964 |
| 1371 | 18.3311714 | 49.0000000 | 7.540457 | ZD | 7.238526 | 2.8247306 |
| 1381 | 21.4111116 | 46.7843710 | 5.504480 | ZD | 9.344123 | 1.3092564 |
| 1391 | 18.8694693 | 43.5405025 | 8.373520 | ZD | 8.301285 | 2.9225026 |
| 1401 | 19.7880615 | 46.1762995 | 5.829197 | ZD | 7.549090 | 1.7794798 |
| 1411 | 7.4560338 | 44.4200915 | 8.199286 | ZD | 8.359500 | 2.8042480 |
| 1421 | 16.6670340 | 47.2940320 | 6.212923 | ZD | 8.655479 | 1.7221941 |
| 1431 | 12.5297550 | 49.0000000 | 8.027718 | ZD | 9.927875 | 2.2811416 |
| 1441 | 15.8983970 | 49.0000000 | 7.926028 | ZD | 9.152429 | 2.4234644 |
| 1451 | 38.9231748 | 41.7665160 | 6.291640 | ZD | 5.676560 | 2.6804741 |
| 1461 | 37.4657599 | 47.2459319 | 8.682012 | ZD | 7.404217 | 3.4758700 |
| 1471 | 49.6668905 | 39.7330850 | 8.203196 | ZD | 11.848520 | 1.9800244 |
| 1481 | 31.8817342 | 39.6527677 | 5.795613 | ZD | 9.619317 | 1.3834013 |
| 1491 | 5.6448613 | 28.8535785 | 5.703050 | ZD | 9.278930 | 1.3969576 |
| 1501 | 4.2739190 | 17.6424662 | 6.565414 | ZD | 11.075167 | 1.4727726 |
| 1511 | 10.9522913 | 19.9384604 | 9.095775 | ZD | 11.712855 | 2.3707069 |
| 1521 | 9.6877684 | 14.2428889 | 6.697895 | ZD | 11.577869 | 1.4555159 |
| 1531 | 7.8820318 | 49.0000000 | 9.916860 | SMQ | 9.756934 | 4.7011145 |
| 1541 | 1.0000000 | 44.4448761 | 6.899585 | SMQ | 9.666726 | 2.5786428 |
| 1551 | 2.3695873 | 48.7241619 | 6.060310 | SMQ | 10.536778 | 1.9022829 |
| 1561 | 5.6845312 | 42.1728318 | 7.417518 | SMQ | 11.126609 | 2.5301780 |
| 1571 | 20.6701393 | 33.8747094 | 8.585438 | SMQ | 9.558668 | 3.7658803 |
| 1581 | 17.3413093 | 28.6130006 | 7.200311 | SMQ | 9.952959 | 2.6906948 |
| 1591 | 24.7963305 | 26.5485412 | 6.546525 | SMQ | 10.111355 | 2.2569088 |
| 1601 | 25.8814922 | 23.0839561 | 7.570419 | SMQ | 10.975876 | 2.6544254 |
| 1611 | 38.1404212 | 18.1701050 | 7.064143 | SMQ | 8.354546 | 3.1042401 |
| 1621 | 44.4997075 | 17.7732229 | 8.575503 | SMQ | 11.831227 | 3.0366083 |
| 1631 | 41.8010491 | 30.4660159 | 8.362883 | SMQ | 9.328978 | 3.6919476 |
| 1641 | 43.8052430 | 18.7144681 | 7.256104 | SMQ | 9.298452 | 2.9177030 |
| 1651 | 27.9742202 | 48.2190508 | 7.243758 | SMQ | 8.766556 | 3.0858836 |
| 1661 | 27.1037613 | 49.0000000 | 8.991668 | SMQ | 7.925497 | 4.9089445 |
| 1671 | 33.9898945 | 48.8409792 | 8.534445 | SMQ | 11.451778 | 3.1120164 |
| 1681 | 37.3406029 | 49.0000000 | 6.711967 | SMQ | 9.818730 | 2.4237531 |
| 1691 | 10.4386867 | 34.3300143 | 5.309605 | SMQ | 9.887581 | 1.6231128 |
| 1701 | 8.0524403 | 27.1832789 | 9.161863 | SMQ | 8.734993 | 4.5966489 |
| 1711 | 1.0000000 | 25.2146483 | 9.680147 | SMQ | 2.613125 | 16.8545260 |
| 1721 | 12.1497107 | 29.6801073 | 6.772546 | SMQ | 10.718650 | 2.2540481 |
2.9 Plot.Voronoi.LAI(minx,maxx,miny,maxy,boundary,b,r)
2.9.1功能介绍：
基于泰森多边形法计算单层林或复层林内单一林层的林地内各位置的叶面积指数，并绘图。
2.9.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
2.9.3例：
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
p=Broadleaf.Korean.pine.LAI::Plot.Voronoi.LAI(minx=0,maxx=50,miny=0,maxy=50,boundary=5, b=b,r=3.5) p1=p+geom_point(data=b,aes(x=x,y=y),size=b$DBH/8,col="grey4",alpha=0.2) p1
2.10 Single.point.Voronoi.LAI(minx,maxx,miny,maxy,boundary,a,b,r)
2.10.1功能介绍：
基于泰森多边形法计算单层林或复层林内单一林层的林地内任意单一位点的叶面积指数。
2.10.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要计算叶面积指数的位点的坐标
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
2.10.3例：
a=matrix(runif(2,5,45),1,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 9.806964 40.51719
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
result=Broadleaf.Korean.pine.LAI::Single.point.Voronoi.LAI(minx=0,maxx=50,miny=0,maxy=50, boundary=5,a=a,b=b,r=3.5) formattable::formattable(result)
| pointx | treex | treey | Species | LAI | |
|---|---|---|---|---|---|
| 29 | 9.806964 | 8.155445 | 40.25 | LX | 0 |
2.11 mult.point.Voronoi.LAI(minx,maxx,miny,maxy,boundary,a,b,r)
2.11.1功能介绍：
基于泰森多边形法计算单层林或复层林内单一林层的林地内任意多个位点的叶面积指数。
2.11.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要计算叶面积指数的位点坐标
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
2.11.3例：
set.seed(100)
a=matrix(runif(40,5,45),20,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
head(a)## x y ## 1 17.310644 26.43245 ## 2 15.306900 33.43215 ## 3 27.092897 26.53395 ## 4 7.255326 34.95889 ## 5 23.741971 21.80406 ## 6 24.350829 11.85681
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
result=Broadleaf.Korean.pine.LAI::mult.point.Voronoi.LAI(minx=0,maxx=50,miny=0,maxy=50, boundary=5,a=a,b=b,r=3.5) formattable::formattable(result)
| pointx | pointy | treex | treey | Species | LAI |
|---|---|---|---|---|---|
| 17.310644 | 26.43245 | 16.31137 | 25.048535 | LS | 10.05055 |
| 15.306900 | 33.43215 | 14.21762 | 29.750000 | LX | 0.00000 |
| 27.092897 | 26.53395 | 26.34198 | 29.750000 | LX | 0.00000 |
| 7.255326 | 34.95889 | 10.18809 | 36.037507 | ZD | 15.39914 |
| 23.741971 | 21.80406 | 23.31089 | 24.500000 | LX | 0.00000 |
| 24.350829 | 11.85681 | 22.67717 | 12.869234 | HS | 15.18164 |
| 37.496105 | 35.81206 | 37.09486 | 34.353282 | ZD | 12.69569 |
| 19.812821 | 40.27814 | 20.27980 | 40.250000 | LX | 0.00000 |
| 26.862344 | 26.96387 | 26.34198 | 29.750000 | LX | 0.00000 |
| 11.810482 | 16.10895 | 15.33661 | 16.214938 | SMQ | 17.00122 |
| 29.999859 | 24.53224 | 27.86636 | 22.076860 | SMQ | 25.00448 |
| 40.286621 | 42.14020 | 38.46633 | 40.250000 | LX | 0.00000 |
| 16.214154 | 18.94768 | 18.14606 | 19.266050 | HS | 35.88500 |
| 20.939516 | 43.16631 | 20.27980 | 40.250000 | LX | 0.00000 |
| 35.502043 | 32.81097 | 37.09486 | 34.353282 | ZD | 12.69569 |
| 31.760868 | 40.57814 | 32.40416 | 40.250000 | LX | 0.00000 |
| 13.184486 | 12.21629 | 10.92124 | 11.592341 | HS | 10.21953 |
| 19.300994 | 30.17563 | 20.27980 | 29.750000 | LX | 0.00000 |
| 19.379005 | 44.58257 | 17.24871 | 45.500000 | LX | 0.00000 |
| 32.611621 | 10.21155 | 31.41985 | 6.278545 | ZD | 20.20589 |
2.12 Local.single.point.Voronoi.LAI(minx,maxx,miny,maxy,boundary,a,b,r,Lr)
2.12.1功能介绍：
基于泰森多边形法计算单层林或复层林内单一林层的林地内任意单一位点的与叶面积有关的局部冠层结构参数。
2.12.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要模拟局部冠层结构参数的位置坐标
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
Lr—提取局部冠层结构参数时的圆的半径
2.12.3例：
set.seed(100)
a=matrix(runif(2,5,45),1,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.31064 15.3069
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
result=Broadleaf.Korean.pine.LAI::Local.single.point.Voronoi.LAI(minx=0,maxx=50,miny=0, maxy=50,boundary=5,a=a,b=b,r=3.5,Lr=1.5) formattable::formattable(result)
| x | y | LAI | Local_min_LAI | Local_max_LAI | Local_mean_LAI | Local_sd_LAI | Gap_percent | Canopy_percent | N_L_Percent | B_L_Percent |
|---|---|---|---|---|---|---|---|---|---|---|
| 17.31064 | 15.3069 | 17.00122 | 17.00122 | 35.885 | 17.20648 | 1.961598 | 0 | 1 | 0.02266904 | 0.977331 |
2.13 Local.mult.point.Voronoi.LAI(minx,maxx,miny,maxy,boundary,a,b,r,Lr)
2.13.1功能介绍：
基于泰森多边形法计算单层林或复层林内单一林层的林地内任意多个位点的与叶面积有关的局部冠层结构参数。
2.13.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要模拟局部冠层结构参数的位置坐标
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
Lr—提取局部冠层结构参数时的圆的半径
2.13.3例：
set.seed(100)
a=matrix(runif(40,5,45),20,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
head(a)## x y ## 1 17.310644 26.43245 ## 2 15.306900 33.43215 ## 3 27.092897 26.53395 ## 4 7.255326 34.95889 ## 5 23.741971 21.80406 ## 6 24.350829 11.85681
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
result=Broadleaf.Korean.pine.LAI::Local.mult.point.Voronoi.LAI(minx=0,maxx=50,miny=0,maxy=50, boundary=5,a=a,b=b,r=3.5,Lr=1.5) formattable::formattable(result)
2.14
Semivariogram.Voronoi.LAI.Single(minx,maxx,miny,maxy,boundary,b,r,seq)
2.14.1功能介绍：
基于泰森多边形法模拟单层林或复层林内单一林层的叶面积指数的半变异函数
2.14.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
seq—模拟半变异函数时的空间分辨率
2.14.3例：
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
Broadleaf.Korean.pine.LAI::Semivariogram.Voronoi.LAI.Single(minx=0,maxx=50,miny=0,maxy=50, boundary=5,b=b,r=3.5,seq=20)
## model psill range ## 1 Nug 0.00000 0.000000 ## 2 Sph 69.79591 6.601658
## [1] "Coefficient_of_Determination= 0.981974004649492"
2.15 Plot.Voronoi.LAI.Sum(minx,maxx,miny,maxy,boundary,b,strata,r)
2.15.1功能介绍：
基于分层泰森多边形法绘制复层林林分叶面积指数图
2.15.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
2.15.3例：
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
p=Broadleaf.Korean.pine.LAI::Plot.Voronoi.LAI.Sum(minx=0,maxx=50,miny=0,maxy=50,boundary=5, b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4)) p1=p+geom_vline(xintercept = c(5,45),linetype=2)+geom_hline(yintercept =c(5,45),linetype=2) p1=p+scale_x_continuous(expand= c(0, 5))+scale_y_continuous(expand= c(0, 5))
## Scale for 'x' is already present. Adding another scale for 'x', which will ## replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which will ## replace the existing scale.
p2=p1+geom_point(data=b,aes(x=x,y=y),size=b$DBH/8,col="grey4",alpha=0.2)+ scale_x_continuous(expand= c(0, 0))+scale_y_continuous(expand= c(0, 0))
## Scale for 'x' is already present. Adding another scale for 'x', which will ## replace the existing scale. ## Scale for 'y' is already present. Adding another scale for 'y', which will ## replace the existing scale.
p2
2.16 Single.point.Voronoi.LAI.sum(minx,maxx,miny,maxy,boundary,a,b,strata,r)
2.16.1功能介绍：
基于分层泰森多边形法计算复层林内单一任意位点的叶面积指数
2.16.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要计算叶面积指数的位点的坐标
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
2.16.3例：
set.seed(100)
a=matrix(runif(2,5,45),1,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.31064 15.3069
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
result=Broadleaf.Korean.pine.LAI::Single.point.Voronoi.LAI.sum(minx=0,maxx=50,miny=0, maxy=50,boundary=5,a=a,b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4)) result
## x y Species LAI Species LAI Species LAI Species LAI ## 1291 17.31064 15.3069 HS 1.776275 LX 0 SMQ 17.00122 LX 0 ## SumLAI Strata_cont ## 1291 18.77749 1.124323
2.17 mult.point.Voronoi.LAI.sum(minx,maxx,miny,maxy,boundary,a,b,strata,r)
2.17.1功能介绍：
基于分层泰森多边形法计算复层林内多个任意位点的叶面积指数
2.17.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要计算叶面积指数的位点的坐标
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
2.17.3 例：
set.seed(100)
a=matrix(runif(40,5,45),20,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.310644 26.43245 ## 2 15.306900 33.43215 ## 3 27.092897 26.53395 ## 4 7.255326 34.95889 ## 5 23.741971 21.80406 ## 6 24.350829 11.85681 ## 7 37.496105 35.81206 ## 8 19.812821 40.27814 ## 9 26.862344 26.96387 ## 10 11.810482 16.10895 ## 11 29.999859 24.53224 ## 12 40.286621 42.14020 ## 13 16.214154 18.94768 ## 14 20.939516 43.16631 ## 15 35.502043 32.81097 ## 16 31.760868 40.57814 ## 17 13.184486 12.21629 ## 18 19.300994 30.17563 ## 19 19.379005 44.58257 ## 20 32.611621 10.21155
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
result=Broadleaf.Korean.pine.LAI::mult.point.Voronoi.LAI.sum(minx=0,maxx=50,miny=0, maxy=50,boundary=5,a=a,b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4)) head(result)
## x y Species LAI Species LAI Species LAI Species LAI ## 1 17.310644 26.43245 LX 0 LX 0.000000 LS 10.05055 LX 0 ## 2 15.306900 33.43215 LX 0 LX 0.000000 LX 0.00000 LX 0 ## 3 27.092897 26.53395 LX 0 LX 0.000000 LX 0.00000 LX 0 ## 4 7.255326 34.95889 LX 0 SMQ 7.813829 ZD 15.39914 LX 0 ## 5 23.741971 21.80406 LX 0 LX 0.000000 LX 0.00000 LX 0 ## 6 24.350829 11.85681 LX 0 LX 0.000000 HS 15.18164 LX 0 ## SumLAI Strata_cont ## 1 10.05055 1.200000 ## 2 0.00000 NaN ## 3 0.00000 NaN ## 4 23.21297 1.065354 ## 5 0.00000 NaN ## 6 15.18164 1.200000
2.18 Local.single.point.Voronoi.LAI.sum(minx,maxx,miny,maxy,boundary,a,b,
strata,r,Lr)
2.18.1功能介绍：
基于分层泰森多边形法计算复层林内单一任意位点的与叶面积指数相关的局部冠层结构特征
2.18.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要模拟局部冠层结构参数的位点的坐标
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
Lr—模拟局部冠层结构参数的圆的半径
2.18.3 例：
set.seed(100)
a=matrix(runif(2,5,45),1,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.31064 15.3069
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
result=Broadleaf.Korean.pine.LAI::Local.single.point.Voronoi.LAI.sum(minx=0,maxx=50, miny=0,maxy=50,boundary=5,a=a,b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4),Lr=1.5) formattable::formattable(result)
2.19
Local.mult.point.Voronoi.LAI.sum(minx,maxx,miny,maxy,boundary,a,b,strata,r,Lr)
2.19.1功能介绍：
基于分层泰森多边形法计算复层林内多个任意位点的与叶面积指数相关的局部冠层结构特征
2.19.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—样地边界缓冲区宽度
a—需要模拟局部冠层结构参数的位点的坐标
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
Lr—模拟局部冠层结构参数的圆的半径
2.19.3 例：
set.seed(100)
a=matrix(runif(40,5,45),20,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.310644 26.43245 ## 2 15.306900 33.43215 ## 3 27.092897 26.53395 ## 4 7.255326 34.95889 ## 5 23.741971 21.80406 ## 6 24.350829 11.85681 ## 7 37.496105 35.81206 ## 8 19.812821 40.27814 ## 9 26.862344 26.96387 ## 10 11.810482 16.10895 ## 11 29.999859 24.53224 ## 12 40.286621 42.14020 ## 13 16.214154 18.94768 ## 14 20.939516 43.16631 ## 15 35.502043 32.81097 ## 16 31.760868 40.57814 ## 17 13.184486 12.21629 ## 18 19.300994 30.17563 ## 19 19.379005 44.58257 ## 20 32.611621 10.21155
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
result=Broadleaf.Korean.pine.LAI::Local.mult.point.Voronoi.LAI.sum(minx=0,maxx=50,miny=0, maxy=50,boundary=5,a=a,b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4),Lr=1.5) formattable::formattable(result)
2.20 Semivariogram.Voronoi.LAI(minx,maxx,miny,maxy,boundary,b,strata,r,seq )
2.20.1功能介绍：
基于分层泰森多边形法模拟复层林的叶面积指数的半变异函数
2.20.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—缓冲区宽度
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
seq—模拟半变异函数时的空间分辨率
2.20.3例：
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
Broadleaf.Korean.pine.LAI::Semivariogram.Voronoi.LAI(minx=0,maxx=50,miny=0,maxy=50, boundary =5,b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4),seq = 20)
## model psill range ## 1 Nug 13.27126 0.000000 ## 2 Sph 99.49735 6.691065
## [1] "Coefficient_of_Determination= 0.977115571818784"
2.21 plot.Local.single.point.Voronoi.LAI(minx,maxx,miny,maxy,a,b,r,Lr)
2.21.1功能介绍：
基于泰森多边形法计算单层林或复层林内单一林层的林地内任意单一位点的与叶面积有关的局部冠层结构参数,并绘图验证。
2.21.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
a—需要模拟局部冠层结构参数的位置坐标
b—样地中的林木坐标、树种、胸径
r—将在林分中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为添加样点构成泰森多边形外接圆的半径
Lr—提取局部冠层结构参数时的圆的半径
2.21.3例：
set.seed(100)
a=matrix(runif(2,5,45),1,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.31064 15.3069
b=Broadleaf.Korean.pine.LAI::b b=subset(b,b$H>20 & b$H<30) b=b[,-4] head(b)
## x y DBH Species ## 3 18.146064 19.266050 52.26441 HS ## 4 27.592828 13.896808 32.36649 HS ## 5 6.108164 25.639671 35.55580 HS ## 6 39.511178 5.185569 53.03717 HS ## 7 18.299115 36.464101 31.79671 HS ## 8 25.667076 4.623330 31.56699 HS
Broadleaf.Korean.pine.LAI::plot.Local.single.point.Voronoi.LAI(minx=0,maxx=50,miny=0, maxy=50,a=a,b=b,r=3.5,Lr=1.5)
## x y LAI Local_min_LAI Local_max_LAI Local_mean_LAI ## 1 17.31064 15.3069 17.00122 17.00122 35.885 17.20648 ## Local_sd_LAI Gap_percent Canopy_percent N_L_Percent B_L_Percent ## 1 1.961598 0 1 0.02266904 0.977331
2.22 plot.Local.single.point.Voronoi.LAI.sum(minx,maxx,miny,maxy,a,b,strata,r,Lr)
2.22.1功能介绍：
基于分层泰森多边形法计算复层林内单一任意位点的与叶面积指数相关的局部冠层结构特征，并绘图
2.22.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
a—需要模拟局部冠层结构参数的位点的坐标
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
Lr—模拟局部冠层结构参数的圆的半径
2.22.3 例：
set.seed(100)
a=matrix(runif(2,5,45),1,2)
colnames(a)=c("x","y")
a=as.data.frame(a)
a## x y ## 1 17.31064 15.3069
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
gc()
## used (Mb) gc trigger (Mb) max used (Mb) ## Ncells 1709832 91.4 2998418 160.2 2998418 160.2 ## Vcells 3033143 23.2 8388608 64.0 8388608 64.0
Broadleaf.Korean.pine.LAI::plot.Local.single.point.Voronoi.LAI.sum(minx=0,maxx=50, miny=0,maxy=50,a=a,b=b,strata=c(10,16,30),r=c(2,2.5,3.5,4),Lr=1.5)
## x y LAI Local_min_LAI Local_max_LAI Local_mean_LAI ## 1 17.31064 15.3069 18.77749 17.00122 35.885 18.58303 ## Local_sd_LAI Gap_percent Canopy_percent N_L_Percent B_L_Percent Strata_cont ## 1 1.961076 0 1 0.09510175 0.9048982 1.140772
2.23 Voronoi.LAI.sum.ISAA(minx, maxx, miny, maxy, boundary, b, strata, r, indis, lag)
2.23.1功能介绍：
基于分层泰森多边形法进行复层林内叶面积指数的增量空间自相关分析
2.23.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—缓冲区宽度
b—样地中的林木坐标、树种、胸径
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
indis—初始距离
lag—滞后距增量
2.23.3 例：
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
gc()
## used (Mb) gc trigger (Mb) max used (Mb) ## Ncells 2322732 124.1 4342478 232 3373384 180.2 ## Vcells 14021044 107.0 50324907 384 50324907 384.0
Broadleaf.Korean.pine.LAI::Voronoi.LAI.sum.ISAA(minx=0, maxx=50, miny=0, maxy=50, boundary=5, b=b, strata=c(10,16,30),r=c(2,2.5,3.5,4), indis=1, lag=1)
2.24 Voronoi.pointcloud(minx, maxx, miny, maxy, boundary, b, seq, strata, r,
S, theta, phi)
2.24.1功能介绍：
基于分层泰森多边形反演叶面积的三维空间分布，并以点云绘图
2.24.2参数意义：
minx—样地横坐标最小范围
maxx—样地横坐标最大范围
miny—样地纵坐标最小范围
maxy—样地纵坐标最大范围
boundary—缓冲区宽度
b—样地中的林木坐标、树种、胸径
seq—绘制点云时的空间分辨率
strata—不同林层间的分割高度
r—将在林分各林层中没有林木分布的位置添加样点，以避免以林木为中心的泰森多边形过大。r为不同林层添加样点构成泰森多边形外接圆的半径
S—林分内平均冠形率，冠长与冠幅的比例
theta—绘图水平翻转角度
phi—绘图垂直翻转角度
2.24.3 例：
b=Broadleaf.Korean.pine.LAI::b head(b)
## x y DBH H Species ## 1 33.428848 18.981981 65.58075 35 HS ## 2 6.434608 42.039176 66.18077 35 HS ## 3 18.146064 19.266050 52.26441 25 HS ## 4 27.592828 13.896808 32.36649 25 HS ## 5 6.108164 25.639671 35.55580 25 HS ## 6 39.511178 5.185569 53.03717 25 HS
Broadleaf.Korean.pine.LAI::Voronoi.pointcloud(minx=0, maxx=50, miny=0, maxy=50, boundary=5, b=b, seq=100, strata=c(10,16,30),r=c(2,2.5,3.5,4), S=1.5, theta=120, phi=20)
Broadleaf.Korean.pine.LAI::Voronoi.pointcloud(minx=0, maxx=50, miny=0, maxy=50, boundary=5, b=b, seq=100, strata=c(10,16,30),r=c(2,2.5,3.5,4), S=1.5, theta=0, phi=90)
https://wap.sciencenet.cn/blog-3511023-1378221.html
下一篇：一种能够消除尺度累积效应的空间权重矩阵，及其基于莫兰指数（Moran‘S I）应用的初步想法
扫一扫，分享此博文