Fragstats4.2景观参数计算与景观生态学意义
“ Fragstats4.2是目前最新的生态学参数计算软件，可以很方便的计算景观生态学相关参数，本文以Fragstats4.2介绍生态参数计算及其意义，还有R语言对计算结果进行制图可视化”
01
景观指数计算
—
打开Analysis paraments选项卡，在Sampling strategy中选择No sampling, 复选框选择Class metrics、Landscape metrics和Generate patchID file
然后从右侧选择不同的景观参数进行计算。
点击开始计算按钮开始计算，然后再点一下Proceed按钮，这样fragstats就开始计算了
切换到Result面板查看计算结果，然后可以对计算结果进行保存。
02
常用的景观参数介绍
—
斑块类型面积(Total Area, CA/TA)。是度量景观组分和计算其他指标的基础, 其值的大小制约着以某种景观类型斑块为聚集地的物种丰度、数量、食物链及其次生种的繁殖等。
式中，ai为某景观斑块面积，n为斑块数。
斑块个数(Number of Patches, NP)。反映景观的空间格局, 经常被用来描述整个景观的异质性, 其值的大小与景观的破碎度有很好的正相关性, 一般NP 越大, 破碎度越高; NP 越低, 破碎度越低。
NP=n
式中，n为某种景观类型总斑块数。
最大斑块指数（LargestPatch Index, LPI），表征了某一类型的最大斑块占整个景观区域面积都比例。
03
运算结果ArcGIS后处理
—
SetNull("shei.tif" == - 999,"shei.tif")需要注意的是，代码函数和符号应该使用鼠标点击，尤其是-999的负号，如果自己使用键盘输入会报错ERROR 000539
由于生成的TIFF栅格数据很多，可以使用模型构建器批量处理，模型搭建情况如下，有关模型构建器循环处理栅格更多信息请参阅ArcGIS批量转投影
04
景观指数R语言可视化
—
如图，绘制一个斑块类型面积的旋风图，需要将数据转化为如下形式：
| TYPE | CA | COUNTRY |
| 其它 | 157.95 | CN |
| 草地 | 45633.78 | CN |
| 其它 | 160360.56 | RU |
| 建设用地 | 165136.32 | CN |
| 建设用地 | 227573.55 | RU |
| 水体 | 271418.04 | CN |
| 湿地 | 295790.13 | CN |
| 草地 | 575991.81 | RU |
| 水体 | 589802.31 | RU |
| 耕地 | 702245.7 | RU |
| 湿地 | 1016431.74 | RU |
| 林地 | 1832953.86 | CN |
| 耕地 | 3250981.17 | CN |
| 林地 | 10242924.66 | RU |
使用R语言ggplot2进行绘图：
library(ggplot2)library(readxl)CAdata <- read_excel(path = "CA.xlsx", sheet = 2, col_names = T)p1 <- ggplot(CAdata, aes(x = TYPE, y = ifelse(COUNTRY== 'RU', CA, -CA), fill=COUNTRY))+ geom_bar(stat = "identity") + coord_flip() p1 + scale_y_continuous(limits = c(-11000000, 11000000), breaks = seq(-11000000, 11000000, 10000000), labels = abs)+ labs(y = 'CA(ha)')p1是未进行坐标轴对称的形式：
完整的代码是修整坐标轴后。