Q3_1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "616e35e8-0c98-401f-a1c0-3823be4d93fa",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import warnings \n",
    "warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "7081be33-abb8-4720-88eb-c2afcc4bcb4f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>year</th>\n",
       "      <th>放牧小区（plot）</th>\n",
       "      <th>放牧强度（intensity）</th>\n",
       "      <th>SOC土壤有机碳</th>\n",
       "      <th>SIC土壤无机碳</th>\n",
       "      <th>STC土壤全碳</th>\n",
       "      <th>全氮N</th>\n",
       "      <th>土壤C/N比</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2012</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>17.051564</td>\n",
       "      <td>11.300000</td>\n",
       "      <td>28.351564</td>\n",
       "      <td>2.192376</td>\n",
       "      <td>12.931887</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2012</td>\n",
       "      <td>G19</td>\n",
       "      <td>NG</td>\n",
       "      <td>19.529439</td>\n",
       "      <td>11.333333</td>\n",
       "      <td>30.862772</td>\n",
       "      <td>1.883051</td>\n",
       "      <td>16.389768</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2012</td>\n",
       "      <td>G21</td>\n",
       "      <td>NG</td>\n",
       "      <td>14.549908</td>\n",
       "      <td>9.600000</td>\n",
       "      <td>24.149908</td>\n",
       "      <td>1.871001</td>\n",
       "      <td>13.818022</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2012</td>\n",
       "      <td>G6</td>\n",
       "      <td>LGI</td>\n",
       "      <td>18.672931</td>\n",
       "      <td>12.666667</td>\n",
       "      <td>31.339598</td>\n",
       "      <td>2.407612</td>\n",
       "      <td>12.160230</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2012</td>\n",
       "      <td>G12</td>\n",
       "      <td>LGI</td>\n",
       "      <td>13.904905</td>\n",
       "      <td>8.133333</td>\n",
       "      <td>22.038239</td>\n",
       "      <td>1.933916</td>\n",
       "      <td>11.395652</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   year 放牧小区（plot） 放牧强度（intensity）   SOC土壤有机碳   SIC土壤无机碳    STC土壤全碳       全氮N  \\\n",
       "0  2012        G17              NG  17.051564  11.300000  28.351564  2.192376   \n",
       "1  2012        G19              NG  19.529439  11.333333  30.862772  1.883051   \n",
       "2  2012        G21              NG  14.549908   9.600000  24.149908  1.871001   \n",
       "3  2012         G6             LGI  18.672931  12.666667  31.339598  2.407612   \n",
       "4  2012        G12             LGI  13.904905   8.133333  22.038239  1.933916   \n",
       "\n",
       "      土壤C/N比  \n",
       "0  12.931887  \n",
       "1  16.389768  \n",
       "2  13.818022  \n",
       "3  12.160230  \n",
       "4  11.395652  "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data=pd.read_excel('附件14不同放牧强度土壤碳氮监测数据.xlsx')\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "3e207394-55f1-4e15-abc2-753a841ecf87",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(132, 8)"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "0f66a5a8-ad66-42f9-b7bb-bc3ee271fc67",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>年份</th>\n",
       "      <th>轮次</th>\n",
       "      <th>处理</th>\n",
       "      <th>日期</th>\n",
       "      <th>植物种名</th>\n",
       "      <th>植物群落功能群</th>\n",
       "      <th>放牧小区Block</th>\n",
       "      <th>重复</th>\n",
       "      <th>营养苗</th>\n",
       "      <th>生殖苗</th>\n",
       "      <th>株/丛数</th>\n",
       "      <th>丛幅1</th>\n",
       "      <th>丛幅2</th>\n",
       "      <th>鲜重(g)</th>\n",
       "      <th>干重(g)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2019</td>\n",
       "      <td>牧前</td>\n",
       "      <td>轻牧（3天）</td>\n",
       "      <td>2019.5.10</td>\n",
       "      <td>大针茅</td>\n",
       "      <td>PB</td>\n",
       "      <td>G6</td>\n",
       "      <td>1</td>\n",
       "      <td>24.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>15.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>59.43</td>\n",
       "      <td>25.41</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2019</td>\n",
       "      <td>牧前</td>\n",
       "      <td>轻牧（3天）</td>\n",
       "      <td>2019.5.10</td>\n",
       "      <td>羊草</td>\n",
       "      <td>PR</td>\n",
       "      <td>G6</td>\n",
       "      <td>1</td>\n",
       "      <td>12.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>17.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.70</td>\n",
       "      <td>0.64</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2019</td>\n",
       "      <td>牧前</td>\n",
       "      <td>轻牧（3天）</td>\n",
       "      <td>2019.5.10</td>\n",
       "      <td>苔草</td>\n",
       "      <td>PR</td>\n",
       "      <td>G6</td>\n",
       "      <td>1</td>\n",
       "      <td>11.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>13.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.27</td>\n",
       "      <td>0.99</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2019</td>\n",
       "      <td>牧前</td>\n",
       "      <td>轻牧（3天）</td>\n",
       "      <td>2019.5.10</td>\n",
       "      <td>糙隐子草</td>\n",
       "      <td>PB</td>\n",
       "      <td>G6</td>\n",
       "      <td>1</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>16.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>2.60</td>\n",
       "      <td>1.25</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2019</td>\n",
       "      <td>牧前</td>\n",
       "      <td>轻牧（3天）</td>\n",
       "      <td>2019.5.10</td>\n",
       "      <td>知母</td>\n",
       "      <td>PF</td>\n",
       "      <td>G6</td>\n",
       "      <td>1</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>9.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.62</td>\n",
       "      <td>0.30</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     年份  轮次      处理         日期  植物种名 植物群落功能群 放牧小区Block 重复   营养苗  生殖苗  株/丛数  \\\n",
       "0  2019  牧前  轻牧（3天）  2019.5.10   大针茅      PB        G6  1  24.0  NaN  15.0   \n",
       "1  2019  牧前  轻牧（3天）  2019.5.10    羊草      PR        G6  1  12.0  NaN  17.0   \n",
       "2  2019  牧前  轻牧（3天）  2019.5.10    苔草      PR        G6  1  11.0  NaN  13.0   \n",
       "3  2019  牧前  轻牧（3天）  2019.5.10  糙隐子草      PB        G6  1   4.0  NaN  16.0   \n",
       "4  2019  牧前  轻牧（3天）  2019.5.10    知母      PF        G6  1   4.0  NaN   9.0   \n",
       "\n",
       "    丛幅1   丛幅2  鲜重(g)  干重(g)  \n",
       "0  11.0  11.0  59.43  25.41  \n",
       "1   NaN   NaN   1.70   0.64  \n",
       "2   NaN   NaN   1.27   0.99  \n",
       "3   4.0   4.0   2.60   1.25  \n",
       "4   NaN   NaN   0.62   0.30  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "f=open('附件15群落结构监测数据集.csv')\n",
    "temp1=pd.read_csv(f,encoding='utf-8')\n",
    "temp1.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "6b8c58b1-b917-4356-9eff-709dce45e366",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['G6-2019', 'G8-2019', 'G9-2019', 'G11-2019', 'G12-2019',\n",
       "       'G13-2019', 'G16-2019', 'G17-2019', 'G18-2019', 'G19-2019',\n",
       "       'G20-2019', 'G21-2019', 'G6-2020', 'G8-2020', 'G9-2020',\n",
       "       'G11-2020', 'G12-2020', 'G13-2020', 'G16-2020', 'G17-2020',\n",
       "       'G18-2020', 'G19-2020', 'G20-2020', 'G21-2020', 'G6-2018',\n",
       "       'G12-2018', 'G18-2018', 'G8-2018', 'G11-2018', 'G16-2018',\n",
       "       'G9-2018', 'G13-2018', 'G20-2018', 'G17-2018', 'G19-2018',\n",
       "       'G21-2018', 'G6-2017', 'G12-2017', 'G18-2017', 'G8-2017',\n",
       "       'G11-2017', 'G16-2017', 'G9-2017', 'G13-2017', 'G20-2017',\n",
       "       'G17-2017', 'G19-2017', 'G21-2017', 'G6-2016', 'G12-2016',\n",
       "       'G18-2016', 'G8-2016', 'G11-2016', 'G16-2016', 'G9-2016',\n",
       "       'G13-2016', 'G17-2016', 'G19-2016', 'G21-2016', 'G20-2016'],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "temp1.apply(lambda x:x['放牧小区Block']+'-'+str(x['年份']),axis=1).unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "3e221ca4-7a05-422e-9ac0-3ef70ef63797",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['G17-2012', 'G19-2012', 'G21-2012', 'G6-2012', 'G12-2012',\n",
       "       'G18-2012', 'G8-2012', 'G11-2012', 'G16-2012', 'G9-2012',\n",
       "       'G13-2012', 'G20-2012', 'G17-2014', 'G19-2014', 'G21-2014',\n",
       "       'G6-2014', 'G12-2014', 'G18-2014', 'G8-2014', 'G11-2014',\n",
       "       'G16-2014', 'G9-2014', 'G13-2014', 'G20-2014', 'G17-2016',\n",
       "       'G19-2016', 'G21-2016', 'G6-2016', 'G12-2016', 'G18-2016',\n",
       "       'G8-2016', 'G11-2016', 'G16-2016', 'G9-2016', 'G13-2016',\n",
       "       'G20-2016', 'G17-2018', 'G19-2018', 'G21-2018', 'G6-2018',\n",
       "       'G12-2018', 'G18-2018', 'G8-2018', 'G11-2018', 'G16-2018',\n",
       "       'G9-2018', 'G13-2018', 'G20-2018', 'G17-2020', 'G19-2020',\n",
       "       'G21-2020', 'G6-2020', 'G12-2020', 'G18-2020', 'G8-2020',\n",
       "       'G11-2020', 'G16-2020', 'G9-2020', 'G13-2020', 'G20-2020'],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.apply(lambda x:x['放牧小区（plot）']+'-'+str(x['year']),axis=1).unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "306c208d-56b5-45d3-862f-fb8f8568f543",
   "metadata": {},
   "outputs": [],
   "source": [
    "temp1['inner']=temp1.apply(lambda x:x['放牧小区Block']+'-'+str(x['年份']),axis=1)\n",
    "data['inner']=data.apply(lambda x:x['放牧小区（plot）']+'-'+str(x['year']),axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "04fde587-21dd-4566-ac56-549c67966257",
   "metadata": {},
   "outputs": [],
   "source": [
    "data2=pd.merge(data,temp1,on='inner')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "bd55a76a-3313-4d27-90c2-5ca9e9c90d8b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>year</th>\n",
       "      <th>放牧小区（plot）</th>\n",
       "      <th>放牧强度（intensity）</th>\n",
       "      <th>SOC土壤有机碳</th>\n",
       "      <th>SIC土壤无机碳</th>\n",
       "      <th>STC土壤全碳</th>\n",
       "      <th>全氮N</th>\n",
       "      <th>土壤C/N比</th>\n",
       "      <th>inner</th>\n",
       "      <th>年份</th>\n",
       "      <th>...</th>\n",
       "      <th>植物群落功能群</th>\n",
       "      <th>放牧小区Block</th>\n",
       "      <th>重复</th>\n",
       "      <th>营养苗</th>\n",
       "      <th>生殖苗</th>\n",
       "      <th>株/丛数</th>\n",
       "      <th>丛幅1</th>\n",
       "      <th>丛幅2</th>\n",
       "      <th>鲜重(g)</th>\n",
       "      <th>干重(g)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>G17-2016</td>\n",
       "      <td>2016</td>\n",
       "      <td>...</td>\n",
       "      <td>PB</td>\n",
       "      <td>G17</td>\n",
       "      <td>1</td>\n",
       "      <td>40.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>10.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>54.90</td>\n",
       "      <td>27.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>G17-2016</td>\n",
       "      <td>2016</td>\n",
       "      <td>...</td>\n",
       "      <td>PR</td>\n",
       "      <td>G17</td>\n",
       "      <td>1</td>\n",
       "      <td>25.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>74.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>33.26</td>\n",
       "      <td>16.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>G17-2016</td>\n",
       "      <td>2016</td>\n",
       "      <td>...</td>\n",
       "      <td>PF</td>\n",
       "      <td>G17</td>\n",
       "      <td>1</td>\n",
       "      <td>16.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>22.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>17.18</td>\n",
       "      <td>4.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>G17-2016</td>\n",
       "      <td>2016</td>\n",
       "      <td>...</td>\n",
       "      <td>PB</td>\n",
       "      <td>G17</td>\n",
       "      <td>1</td>\n",
       "      <td>7.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>9.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>4.70</td>\n",
       "      <td>2.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>G17-2016</td>\n",
       "      <td>2016</td>\n",
       "      <td>...</td>\n",
       "      <td>PF</td>\n",
       "      <td>G17</td>\n",
       "      <td>1</td>\n",
       "      <td>30.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3.26</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14292</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>G20-2020</td>\n",
       "      <td>2020</td>\n",
       "      <td>...</td>\n",
       "      <td>PB</td>\n",
       "      <td>G20</td>\n",
       "      <td>笼2</td>\n",
       "      <td>12.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>51.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>14.79</td>\n",
       "      <td>8.12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14293</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>G20-2020</td>\n",
       "      <td>2020</td>\n",
       "      <td>...</td>\n",
       "      <td>PR</td>\n",
       "      <td>G20</td>\n",
       "      <td>笼2</td>\n",
       "      <td>24.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>121.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>37.58</td>\n",
       "      <td>19.07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14294</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>G20-2020</td>\n",
       "      <td>2020</td>\n",
       "      <td>...</td>\n",
       "      <td>PF</td>\n",
       "      <td>G20</td>\n",
       "      <td>笼2</td>\n",
       "      <td>11.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>8.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>4.69</td>\n",
       "      <td>0.85</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14295</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>G20-2020</td>\n",
       "      <td>2020</td>\n",
       "      <td>...</td>\n",
       "      <td>AB</td>\n",
       "      <td>G20</td>\n",
       "      <td>笼2</td>\n",
       "      <td>13.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.09</td>\n",
       "      <td>0.05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14296</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>G20-2020</td>\n",
       "      <td>2020</td>\n",
       "      <td>...</td>\n",
       "      <td>AB</td>\n",
       "      <td>G20</td>\n",
       "      <td>笼2</td>\n",
       "      <td>13.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>66.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>13.48</td>\n",
       "      <td>9.56</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>14297 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       year 放牧小区（plot） 放牧强度（intensity）   SOC土壤有机碳  SIC土壤无机碳    STC土壤全碳  \\\n",
       "0      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "1      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "2      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "3      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "4      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "...     ...        ...             ...        ...       ...        ...   \n",
       "14292  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14293  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14294  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14295  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14296  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "\n",
       "            全氮N     土壤C/N比     inner    年份  ... 植物群落功能群 放牧小区Block  重复   营养苗  \\\n",
       "0      2.332079   9.040030  G17-2016  2016  ...      PB       G17   1  40.0   \n",
       "1      2.332079   9.040030  G17-2016  2016  ...      PR       G17   1  25.0   \n",
       "2      2.332079   9.040030  G17-2016  2016  ...      PF       G17   1  16.0   \n",
       "3      2.332079   9.040030  G17-2016  2016  ...      PB       G17   1   7.0   \n",
       "4      2.332079   9.040030  G17-2016  2016  ...      PF       G17   1  30.0   \n",
       "...         ...        ...       ...   ...  ...     ...       ...  ..   ...   \n",
       "14292  2.054600  10.039253  G20-2020  2020  ...      PB       G20  笼2  12.0   \n",
       "14293  2.054600  10.039253  G20-2020  2020  ...      PR       G20  笼2  24.0   \n",
       "14294  2.054600  10.039253  G20-2020  2020  ...      PF       G20  笼2  11.0   \n",
       "14295  2.054600  10.039253  G20-2020  2020  ...      AB       G20  笼2  13.0   \n",
       "14296  2.054600  10.039253  G20-2020  2020  ...      AB       G20  笼2  13.0   \n",
       "\n",
       "      生殖苗   株/丛数  丛幅1  丛幅2  鲜重(g)  干重(g)  \n",
       "0     NaN   10.0  5.0  5.0  54.90  27.00  \n",
       "1     NaN   74.0  NaN  NaN  33.26  16.00  \n",
       "2     NaN   22.0  NaN  NaN  17.18   4.00  \n",
       "3     NaN    9.0  3.0  3.0   4.70   2.00  \n",
       "4     NaN    1.0  NaN  NaN   3.26   1.00  \n",
       "...    ..    ...  ...  ...    ...    ...  \n",
       "14292 NaN   51.0  7.0  7.0  14.79   8.12  \n",
       "14293 NaN  121.0  NaN  NaN  37.58  19.07  \n",
       "14294 NaN    8.0  NaN  NaN   4.69   0.85  \n",
       "14295 NaN    2.0  NaN  NaN   0.09   0.05  \n",
       "14296 NaN   66.0  NaN  NaN  13.48   9.56  \n",
       "\n",
       "[14297 rows x 24 columns]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "95e5446d-b24c-4ad1-bcb5-a1f9ff0d2cc4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "year                   0\n",
       "放牧小区（plot）             0\n",
       "放牧强度（intensity）        0\n",
       "SOC土壤有机碳               0\n",
       "SIC土壤无机碳               0\n",
       "STC土壤全碳                0\n",
       "全氮N                    0\n",
       "土壤C/N比                 0\n",
       "inner                  0\n",
       "年份                     0\n",
       "轮次                     0\n",
       "处理                     0\n",
       "日期                     0\n",
       "植物种名                   0\n",
       "植物群落功能群                0\n",
       "放牧小区Block              0\n",
       "重复                     0\n",
       "营养苗                   24\n",
       "生殖苗                12581\n",
       "株/丛数                 357\n",
       "丛幅1                 9438\n",
       "丛幅2                10054\n",
       "鲜重(g)                705\n",
       "干重(g)                  0\n",
       "dtype: int64"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2.isnull().sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "d6b41475-faf5-434e-accb-13c870475a31",
   "metadata": {},
   "outputs": [],
   "source": [
    "data2=data2[[ '年份', '放牧小区（plot）', '放牧强度（intensity）', 'SOC土壤有机碳', 'SIC土壤无机碳',\n",
    "       'STC土壤全碳', '全氮N', '土壤C/N比',  '轮次', '处理',  '植物种名',\n",
    "       '植物群落功能群',  '重复', '营养苗',  '株/丛数',\n",
    "       '鲜重(g)', '干重(g)']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "304243f4-7605-4c75-b1a6-1bc04eb673d0",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>年份</th>\n",
       "      <th>放牧小区（plot）</th>\n",
       "      <th>放牧强度（intensity）</th>\n",
       "      <th>SOC土壤有机碳</th>\n",
       "      <th>SIC土壤无机碳</th>\n",
       "      <th>STC土壤全碳</th>\n",
       "      <th>全氮N</th>\n",
       "      <th>土壤C/N比</th>\n",
       "      <th>轮次</th>\n",
       "      <th>处理</th>\n",
       "      <th>植物种名</th>\n",
       "      <th>植物群落功能群</th>\n",
       "      <th>重复</th>\n",
       "      <th>营养苗</th>\n",
       "      <th>株/丛数</th>\n",
       "      <th>鲜重(g)</th>\n",
       "      <th>干重(g)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>大针茅</td>\n",
       "      <td>PB</td>\n",
       "      <td>1</td>\n",
       "      <td>40.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>54.90</td>\n",
       "      <td>27.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>羊草</td>\n",
       "      <td>PR</td>\n",
       "      <td>1</td>\n",
       "      <td>25.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>33.26</td>\n",
       "      <td>16.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>知母</td>\n",
       "      <td>PF</td>\n",
       "      <td>1</td>\n",
       "      <td>16.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>17.18</td>\n",
       "      <td>4.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>糙隐子草</td>\n",
       "      <td>PB</td>\n",
       "      <td>1</td>\n",
       "      <td>7.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>4.70</td>\n",
       "      <td>2.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.40000</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.040030</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>黄花葱</td>\n",
       "      <td>PF</td>\n",
       "      <td>1</td>\n",
       "      <td>30.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>3.26</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14292</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>第四轮牧后</td>\n",
       "      <td>重牧（12天）</td>\n",
       "      <td>糙隐子草</td>\n",
       "      <td>PB</td>\n",
       "      <td>笼2</td>\n",
       "      <td>12.0</td>\n",
       "      <td>51.0</td>\n",
       "      <td>14.79</td>\n",
       "      <td>8.12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14293</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>第四轮牧后</td>\n",
       "      <td>重牧（12天）</td>\n",
       "      <td>羊草</td>\n",
       "      <td>PR</td>\n",
       "      <td>笼2</td>\n",
       "      <td>24.0</td>\n",
       "      <td>121.0</td>\n",
       "      <td>37.58</td>\n",
       "      <td>19.07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14294</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>第四轮牧后</td>\n",
       "      <td>重牧（12天）</td>\n",
       "      <td>矮葱</td>\n",
       "      <td>PF</td>\n",
       "      <td>笼2</td>\n",
       "      <td>11.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>4.69</td>\n",
       "      <td>0.85</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14295</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>第四轮牧后</td>\n",
       "      <td>重牧（12天）</td>\n",
       "      <td>画眉草</td>\n",
       "      <td>AB</td>\n",
       "      <td>笼2</td>\n",
       "      <td>13.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>0.09</td>\n",
       "      <td>0.05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14296</th>\n",
       "      <td>2020</td>\n",
       "      <td>G20</td>\n",
       "      <td>HGI</td>\n",
       "      <td>15.969600</td>\n",
       "      <td>4.65705</td>\n",
       "      <td>20.626650</td>\n",
       "      <td>2.054600</td>\n",
       "      <td>10.039253</td>\n",
       "      <td>第四轮牧后</td>\n",
       "      <td>重牧（12天）</td>\n",
       "      <td>刺藜</td>\n",
       "      <td>AB</td>\n",
       "      <td>笼2</td>\n",
       "      <td>13.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>13.48</td>\n",
       "      <td>9.56</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>14297 rows × 17 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         年份 放牧小区（plot） 放牧强度（intensity）   SOC土壤有机碳  SIC土壤无机碳    STC土壤全碳  \\\n",
       "0      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "1      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "2      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "3      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "4      2016        G17              NG  18.682066   2.40000  21.082066   \n",
       "...     ...        ...             ...        ...       ...        ...   \n",
       "14292  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14293  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14294  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14295  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "14296  2020        G20             HGI  15.969600   4.65705  20.626650   \n",
       "\n",
       "            全氮N     土壤C/N比     轮次       处理  植物种名 植物群落功能群  重复   营养苗   株/丛数  \\\n",
       "0      2.332079   9.040030     牧前   无牧（0天）   大针茅      PB   1  40.0   10.0   \n",
       "1      2.332079   9.040030     牧前   无牧（0天）    羊草      PR   1  25.0   74.0   \n",
       "2      2.332079   9.040030     牧前   无牧（0天）    知母      PF   1  16.0   22.0   \n",
       "3      2.332079   9.040030     牧前   无牧（0天）  糙隐子草      PB   1   7.0    9.0   \n",
       "4      2.332079   9.040030     牧前   无牧（0天）   黄花葱      PF   1  30.0    1.0   \n",
       "...         ...        ...    ...      ...   ...     ...  ..   ...    ...   \n",
       "14292  2.054600  10.039253  第四轮牧后  重牧（12天）  糙隐子草      PB  笼2  12.0   51.0   \n",
       "14293  2.054600  10.039253  第四轮牧后  重牧（12天）    羊草      PR  笼2  24.0  121.0   \n",
       "14294  2.054600  10.039253  第四轮牧后  重牧（12天）    矮葱      PF  笼2  11.0    8.0   \n",
       "14295  2.054600  10.039253  第四轮牧后  重牧（12天）   画眉草      AB  笼2  13.0    2.0   \n",
       "14296  2.054600  10.039253  第四轮牧后  重牧（12天）    刺藜      AB  笼2  13.0   66.0   \n",
       "\n",
       "       鲜重(g)  干重(g)  \n",
       "0      54.90  27.00  \n",
       "1      33.26  16.00  \n",
       "2      17.18   4.00  \n",
       "3       4.70   2.00  \n",
       "4       3.26   1.00  \n",
       "...      ...    ...  \n",
       "14292  14.79   8.12  \n",
       "14293  37.58  19.07  \n",
       "14294   4.69   0.85  \n",
       "14295   0.09   0.05  \n",
       "14296  13.48   9.56  \n",
       "\n",
       "[14297 rows x 17 columns]"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e6e56370-3dbf-45d8-a5b8-01446be09e0c",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "d416b9b4-fa70-4cd6-ae3b-badd6b6f586d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(14297, 17)"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "3fc9b721-bab9-45f3-81d1-0ff5258ba56b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>10cm湿度(kg/m2)</th>\n",
       "      <th>40cm湿度(kg/m2)</th>\n",
       "      <th>100cm湿度(kg/m2)</th>\n",
       "      <th>200cm湿度(kg/m2)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.45</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.10</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.96</td>\n",
       "      <td>52.13</td>\n",
       "      <td>93.44</td>\n",
       "      <td>164.48</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.17</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>13.93</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   月份    年份  经度(lon)  纬度(lat)  10cm湿度(kg/m2)  40cm湿度(kg/m2)  100cm湿度(kg/m2)  \\\n",
       "0   1  2022  115.375   44.125          12.45          52.14           93.45   \n",
       "1   2  2022  115.375   44.125          12.10          52.14           93.45   \n",
       "2   3  2022  115.375   44.125          14.96          52.13           93.44   \n",
       "3   1  2021  115.375   44.125          15.17          50.40           82.44   \n",
       "4   2  2021  115.375   44.125          13.93          50.40           82.44   \n",
       "\n",
       "   200cm湿度(kg/m2)  \n",
       "0          164.48  \n",
       "1          164.48  \n",
       "2          164.48  \n",
       "3          165.92  \n",
       "4          165.92  "
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "temp2=pd.read_excel('附件3、土壤湿度2022—2012年.xls')\n",
    "temp2.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "79f68212-b03d-4c74-829c-b0290e3a815e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>土壤蒸发量(W/m2)</th>\n",
       "      <th>土壤蒸发量(mm)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.41</td>\n",
       "      <td>0.44</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.28</td>\n",
       "      <td>0.27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>5.58</td>\n",
       "      <td>5.98</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>1.05</td>\n",
       "      <td>1.13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>1.58</td>\n",
       "      <td>1.53</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   月份    年份  经度(lon)  纬度(lat)  土壤蒸发量(W/m2)  土壤蒸发量(mm)\n",
       "0   1  2022  115.375   44.125         0.41       0.44\n",
       "1   2  2022  115.375   44.125         0.28       0.27\n",
       "2   3  2022  115.375   44.125         5.58       5.98\n",
       "3   1  2021  115.375   44.125         1.05       1.13\n",
       "4   2  2021  115.375   44.125         1.58       1.53"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "temp3=pd.read_excel('附件4、土壤蒸发量2012—2022年.xls')\n",
    "temp3.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "8cedf56f-9fe7-4f39-9a0b-912ec34a84b8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>植被指数(NDVI)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.010</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.010</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.187</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.191</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>-0.002</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   月份    年份  经度(lon)  纬度(lat)  植被指数(NDVI)\n",
       "0   1  2022  115.375   44.125       0.010\n",
       "1   2  2022  115.375   44.125       0.010\n",
       "2   3  2022  115.375   44.125       0.187\n",
       "3   4  2022  115.375   44.125       0.191\n",
       "4   1  2021  115.375   44.125      -0.002"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "temp4=pd.read_excel('附件6、植被指数-NDVI2012-2022年.xls')\n",
    "temp4.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "10c9e810-24ec-426b-8140-c44d439dd1de",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>径流量(m3/s)</th>\n",
       "      <th>径流量(m3)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.160786</td>\n",
       "      <td>430650.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.116071</td>\n",
       "      <td>280800.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.825605</td>\n",
       "      <td>2211300.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.066532</td>\n",
       "      <td>178200.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>0.033482</td>\n",
       "      <td>81000.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   月份    年份  经度(lon)  纬度(lat)  径流量(m3/s)    径流量(m3)\n",
       "0   1  2022  115.375   44.125   0.160786   430650.0\n",
       "1   2  2022  115.375   44.125   0.116071   280800.0\n",
       "2   3  2022  115.375   44.125   0.825605  2211300.0\n",
       "3   1  2021  115.375   44.125   0.066532   178200.0\n",
       "4   2  2021  115.375   44.125   0.033482    81000.0"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "temp5=pd.read_excel('附件9、径流量2012-2022年.xlsx')\n",
    "temp5.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "3ba79597-0200-4f55-92f7-f224f0353644",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>10cm湿度(kg/m2)</th>\n",
       "      <th>40cm湿度(kg/m2)</th>\n",
       "      <th>100cm湿度(kg/m2)</th>\n",
       "      <th>200cm湿度(kg/m2)</th>\n",
       "      <th>土壤蒸发量(W/m2)</th>\n",
       "      <th>土壤蒸发量(mm)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.45</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "      <td>0.41</td>\n",
       "      <td>0.44</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.10</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "      <td>0.28</td>\n",
       "      <td>0.27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.96</td>\n",
       "      <td>52.13</td>\n",
       "      <td>93.44</td>\n",
       "      <td>164.48</td>\n",
       "      <td>5.58</td>\n",
       "      <td>5.98</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.17</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "      <td>1.05</td>\n",
       "      <td>1.13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>13.93</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "      <td>1.58</td>\n",
       "      <td>1.53</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   月份    年份  经度(lon)  纬度(lat)  10cm湿度(kg/m2)  40cm湿度(kg/m2)  100cm湿度(kg/m2)  \\\n",
       "0   1  2022  115.375   44.125          12.45          52.14           93.45   \n",
       "1   2  2022  115.375   44.125          12.10          52.14           93.45   \n",
       "2   3  2022  115.375   44.125          14.96          52.13           93.44   \n",
       "3   1  2021  115.375   44.125          15.17          50.40           82.44   \n",
       "4   2  2021  115.375   44.125          13.93          50.40           82.44   \n",
       "\n",
       "   200cm湿度(kg/m2)  土壤蒸发量(W/m2)  土壤蒸发量(mm)  \n",
       "0          164.48         0.41       0.44  \n",
       "1          164.48         0.28       0.27  \n",
       "2          164.48         5.58       5.98  \n",
       "3          165.92         1.05       1.13  \n",
       "4          165.92         1.58       1.53  "
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tempdf=pd.merge(temp2,temp3)\n",
    "tempdf.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "66eb4e76-3737-4a04-a7c6-81e8d18813d4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>10cm湿度(kg/m2)</th>\n",
       "      <th>40cm湿度(kg/m2)</th>\n",
       "      <th>100cm湿度(kg/m2)</th>\n",
       "      <th>200cm湿度(kg/m2)</th>\n",
       "      <th>土壤蒸发量(W/m2)</th>\n",
       "      <th>土壤蒸发量(mm)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.45</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "      <td>0.41</td>\n",
       "      <td>0.44</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.10</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "      <td>0.28</td>\n",
       "      <td>0.27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.96</td>\n",
       "      <td>52.13</td>\n",
       "      <td>93.44</td>\n",
       "      <td>164.48</td>\n",
       "      <td>5.58</td>\n",
       "      <td>5.98</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.17</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "      <td>1.05</td>\n",
       "      <td>1.13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>13.93</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "      <td>1.58</td>\n",
       "      <td>1.53</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>118</th>\n",
       "      <td>8</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.76</td>\n",
       "      <td>44.94</td>\n",
       "      <td>61.25</td>\n",
       "      <td>168.44</td>\n",
       "      <td>13.32</td>\n",
       "      <td>14.27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>119</th>\n",
       "      <td>9</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>17.12</td>\n",
       "      <td>45.67</td>\n",
       "      <td>55.29</td>\n",
       "      <td>168.38</td>\n",
       "      <td>12.88</td>\n",
       "      <td>13.36</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>120</th>\n",
       "      <td>10</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.04</td>\n",
       "      <td>42.49</td>\n",
       "      <td>54.69</td>\n",
       "      <td>168.33</td>\n",
       "      <td>7.59</td>\n",
       "      <td>8.13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>121</th>\n",
       "      <td>11</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.79</td>\n",
       "      <td>42.71</td>\n",
       "      <td>55.02</td>\n",
       "      <td>168.32</td>\n",
       "      <td>0.88</td>\n",
       "      <td>0.92</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122</th>\n",
       "      <td>12</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.42</td>\n",
       "      <td>42.71</td>\n",
       "      <td>55.02</td>\n",
       "      <td>168.32</td>\n",
       "      <td>0.32</td>\n",
       "      <td>0.34</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>123 rows × 10 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     月份    年份  经度(lon)  纬度(lat)  10cm湿度(kg/m2)  40cm湿度(kg/m2)  100cm湿度(kg/m2)  \\\n",
       "0     1  2022  115.375   44.125          12.45          52.14           93.45   \n",
       "1     2  2022  115.375   44.125          12.10          52.14           93.45   \n",
       "2     3  2022  115.375   44.125          14.96          52.13           93.44   \n",
       "3     1  2021  115.375   44.125          15.17          50.40           82.44   \n",
       "4     2  2021  115.375   44.125          13.93          50.40           82.44   \n",
       "..   ..   ...      ...      ...            ...            ...             ...   \n",
       "118   8  2012  115.375   44.125          15.76          44.94           61.25   \n",
       "119   9  2012  115.375   44.125          17.12          45.67           55.29   \n",
       "120  10  2012  115.375   44.125          15.04          42.49           54.69   \n",
       "121  11  2012  115.375   44.125          14.79          42.71           55.02   \n",
       "122  12  2012  115.375   44.125          14.42          42.71           55.02   \n",
       "\n",
       "     200cm湿度(kg/m2)  土壤蒸发量(W/m2)  土壤蒸发量(mm)  \n",
       "0            164.48         0.41       0.44  \n",
       "1            164.48         0.28       0.27  \n",
       "2            164.48         5.58       5.98  \n",
       "3            165.92         1.05       1.13  \n",
       "4            165.92         1.58       1.53  \n",
       "..              ...          ...        ...  \n",
       "118          168.44        13.32      14.27  \n",
       "119          168.38        12.88      13.36  \n",
       "120          168.33         7.59       8.13  \n",
       "121          168.32         0.88       0.92  \n",
       "122          168.32         0.32       0.34  \n",
       "\n",
       "[123 rows x 10 columns]"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tempdf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "733ef691-6a9b-4f74-8e7b-f9a365e432be",
   "metadata": {},
   "outputs": [],
   "source": [
    "tempdf=pd.merge(tempdf,temp4)\n",
    "tempdf=pd.merge(tempdf,temp5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "465f5684-1c4d-4c42-a0a9-af3b243d06fe",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>10cm湿度(kg/m2)</th>\n",
       "      <th>40cm湿度(kg/m2)</th>\n",
       "      <th>100cm湿度(kg/m2)</th>\n",
       "      <th>200cm湿度(kg/m2)</th>\n",
       "      <th>土壤蒸发量(W/m2)</th>\n",
       "      <th>土壤蒸发量(mm)</th>\n",
       "      <th>植被指数(NDVI)</th>\n",
       "      <th>径流量(m3/s)</th>\n",
       "      <th>径流量(m3)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.45</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "      <td>0.41</td>\n",
       "      <td>0.44</td>\n",
       "      <td>0.010</td>\n",
       "      <td>0.160786</td>\n",
       "      <td>430650.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.10</td>\n",
       "      <td>52.14</td>\n",
       "      <td>93.45</td>\n",
       "      <td>164.48</td>\n",
       "      <td>0.28</td>\n",
       "      <td>0.27</td>\n",
       "      <td>0.010</td>\n",
       "      <td>0.116071</td>\n",
       "      <td>280800.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2022</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.96</td>\n",
       "      <td>52.13</td>\n",
       "      <td>93.44</td>\n",
       "      <td>164.48</td>\n",
       "      <td>5.58</td>\n",
       "      <td>5.98</td>\n",
       "      <td>0.187</td>\n",
       "      <td>0.825605</td>\n",
       "      <td>2211300.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.17</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "      <td>1.05</td>\n",
       "      <td>1.13</td>\n",
       "      <td>-0.002</td>\n",
       "      <td>0.066532</td>\n",
       "      <td>178200.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2</td>\n",
       "      <td>2021</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>13.93</td>\n",
       "      <td>50.40</td>\n",
       "      <td>82.44</td>\n",
       "      <td>165.92</td>\n",
       "      <td>1.58</td>\n",
       "      <td>1.53</td>\n",
       "      <td>0.002</td>\n",
       "      <td>0.033482</td>\n",
       "      <td>81000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>118</th>\n",
       "      <td>8</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.76</td>\n",
       "      <td>44.94</td>\n",
       "      <td>61.25</td>\n",
       "      <td>168.44</td>\n",
       "      <td>13.32</td>\n",
       "      <td>14.27</td>\n",
       "      <td>0.597</td>\n",
       "      <td>0.621472</td>\n",
       "      <td>1664550.10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>119</th>\n",
       "      <td>9</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>17.12</td>\n",
       "      <td>45.67</td>\n",
       "      <td>55.29</td>\n",
       "      <td>168.38</td>\n",
       "      <td>12.88</td>\n",
       "      <td>13.36</td>\n",
       "      <td>0.475</td>\n",
       "      <td>0.444010</td>\n",
       "      <td>1150875.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>120</th>\n",
       "      <td>10</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>15.04</td>\n",
       "      <td>42.49</td>\n",
       "      <td>54.69</td>\n",
       "      <td>168.33</td>\n",
       "      <td>7.59</td>\n",
       "      <td>8.13</td>\n",
       "      <td>0.282</td>\n",
       "      <td>0.361643</td>\n",
       "      <td>968624.94</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>121</th>\n",
       "      <td>11</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.79</td>\n",
       "      <td>42.71</td>\n",
       "      <td>55.02</td>\n",
       "      <td>168.32</td>\n",
       "      <td>0.88</td>\n",
       "      <td>0.92</td>\n",
       "      <td>0.176</td>\n",
       "      <td>0.290104</td>\n",
       "      <td>751950.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122</th>\n",
       "      <td>12</td>\n",
       "      <td>2012</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>14.42</td>\n",
       "      <td>42.71</td>\n",
       "      <td>55.02</td>\n",
       "      <td>168.32</td>\n",
       "      <td>0.32</td>\n",
       "      <td>0.34</td>\n",
       "      <td>0.018</td>\n",
       "      <td>0.213710</td>\n",
       "      <td>572400.00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>123 rows × 13 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     月份    年份  经度(lon)  纬度(lat)  10cm湿度(kg/m2)  40cm湿度(kg/m2)  100cm湿度(kg/m2)  \\\n",
       "0     1  2022  115.375   44.125          12.45          52.14           93.45   \n",
       "1     2  2022  115.375   44.125          12.10          52.14           93.45   \n",
       "2     3  2022  115.375   44.125          14.96          52.13           93.44   \n",
       "3     1  2021  115.375   44.125          15.17          50.40           82.44   \n",
       "4     2  2021  115.375   44.125          13.93          50.40           82.44   \n",
       "..   ..   ...      ...      ...            ...            ...             ...   \n",
       "118   8  2012  115.375   44.125          15.76          44.94           61.25   \n",
       "119   9  2012  115.375   44.125          17.12          45.67           55.29   \n",
       "120  10  2012  115.375   44.125          15.04          42.49           54.69   \n",
       "121  11  2012  115.375   44.125          14.79          42.71           55.02   \n",
       "122  12  2012  115.375   44.125          14.42          42.71           55.02   \n",
       "\n",
       "     200cm湿度(kg/m2)  土壤蒸发量(W/m2)  土壤蒸发量(mm)  植被指数(NDVI)  径流量(m3/s)     径流量(m3)  \n",
       "0            164.48         0.41       0.44       0.010   0.160786   430650.00  \n",
       "1            164.48         0.28       0.27       0.010   0.116071   280800.00  \n",
       "2            164.48         5.58       5.98       0.187   0.825605  2211300.00  \n",
       "3            165.92         1.05       1.13      -0.002   0.066532   178200.00  \n",
       "4            165.92         1.58       1.53       0.002   0.033482    81000.00  \n",
       "..              ...          ...        ...         ...        ...         ...  \n",
       "118          168.44        13.32      14.27       0.597   0.621472  1664550.10  \n",
       "119          168.38        12.88      13.36       0.475   0.444010  1150875.00  \n",
       "120          168.33         7.59       8.13       0.282   0.361643   968624.94  \n",
       "121          168.32         0.88       0.92       0.176   0.290104   751950.00  \n",
       "122          168.32         0.32       0.34       0.018   0.213710   572400.00  \n",
       "\n",
       "[123 rows x 13 columns]"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tempdf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "fba41ae5-992b-4abe-b62c-30192cc90fec",
   "metadata": {},
   "outputs": [],
   "source": [
    "data2=data2[['年份', '放牧小区（plot）', '放牧强度（intensity）', 'SOC土壤有机碳', 'SIC土壤无机碳',\n",
    "       'STC土壤全碳', '全氮N', '土壤C/N比', '轮次', '处理', '植物种名', '植物群落功能群', '重复', '营养苗',\n",
    "       '株/丛数', '鲜重(g)', '干重(g)']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "e80984a2-5662-42be-9cdf-874a74e1ac6a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>年份</th>\n",
       "      <th>放牧小区（plot）</th>\n",
       "      <th>放牧强度（intensity）</th>\n",
       "      <th>SOC土壤有机碳</th>\n",
       "      <th>SIC土壤无机碳</th>\n",
       "      <th>STC土壤全碳</th>\n",
       "      <th>全氮N</th>\n",
       "      <th>土壤C/N比</th>\n",
       "      <th>轮次</th>\n",
       "      <th>处理</th>\n",
       "      <th>植物种名</th>\n",
       "      <th>植物群落功能群</th>\n",
       "      <th>重复</th>\n",
       "      <th>营养苗</th>\n",
       "      <th>株/丛数</th>\n",
       "      <th>鲜重(g)</th>\n",
       "      <th>干重(g)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.4</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.04003</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>大针茅</td>\n",
       "      <td>PB</td>\n",
       "      <td>1</td>\n",
       "      <td>40.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>54.90</td>\n",
       "      <td>27.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.4</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.04003</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>羊草</td>\n",
       "      <td>PR</td>\n",
       "      <td>1</td>\n",
       "      <td>25.0</td>\n",
       "      <td>74.0</td>\n",
       "      <td>33.26</td>\n",
       "      <td>16.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.4</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.04003</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>知母</td>\n",
       "      <td>PF</td>\n",
       "      <td>1</td>\n",
       "      <td>16.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>17.18</td>\n",
       "      <td>4.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.4</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.04003</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>糙隐子草</td>\n",
       "      <td>PB</td>\n",
       "      <td>1</td>\n",
       "      <td>7.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>4.70</td>\n",
       "      <td>2.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2016</td>\n",
       "      <td>G17</td>\n",
       "      <td>NG</td>\n",
       "      <td>18.682066</td>\n",
       "      <td>2.4</td>\n",
       "      <td>21.082066</td>\n",
       "      <td>2.332079</td>\n",
       "      <td>9.04003</td>\n",
       "      <td>牧前</td>\n",
       "      <td>无牧（0天）</td>\n",
       "      <td>黄花葱</td>\n",
       "      <td>PF</td>\n",
       "      <td>1</td>\n",
       "      <td>30.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>3.26</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     年份 放牧小区（plot） 放牧强度（intensity）   SOC土壤有机碳  SIC土壤无机碳    STC土壤全碳       全氮N  \\\n",
       "0  2016        G17              NG  18.682066       2.4  21.082066  2.332079   \n",
       "1  2016        G17              NG  18.682066       2.4  21.082066  2.332079   \n",
       "2  2016        G17              NG  18.682066       2.4  21.082066  2.332079   \n",
       "3  2016        G17              NG  18.682066       2.4  21.082066  2.332079   \n",
       "4  2016        G17              NG  18.682066       2.4  21.082066  2.332079   \n",
       "\n",
       "    土壤C/N比  轮次      处理  植物种名 植物群落功能群 重复   营养苗  株/丛数  鲜重(g)  干重(g)  \n",
       "0  9.04003  牧前  无牧（0天）   大针茅      PB  1  40.0  10.0  54.90   27.0  \n",
       "1  9.04003  牧前  无牧（0天）    羊草      PR  1  25.0  74.0  33.26   16.0  \n",
       "2  9.04003  牧前  无牧（0天）    知母      PF  1  16.0  22.0  17.18    4.0  \n",
       "3  9.04003  牧前  无牧（0天）  糙隐子草      PB  1   7.0   9.0   4.70    2.0  \n",
       "4  9.04003  牧前  无牧（0天）   黄花葱      PF  1  30.0   1.0   3.26    1.0  "
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "952e634d-3156-4bff-ad63-44c2364bedd8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>plot__G11</th>\n",
       "      <th>plot__G12</th>\n",
       "      <th>plot__G13</th>\n",
       "      <th>plot__G16</th>\n",
       "      <th>plot__G17</th>\n",
       "      <th>plot__G18</th>\n",
       "      <th>plot__G19</th>\n",
       "      <th>plot__G20</th>\n",
       "      <th>plot__G21</th>\n",
       "      <th>plot__G6</th>\n",
       "      <th>plot__G8</th>\n",
       "      <th>plot__G9</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14292</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14293</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14294</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14295</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14296</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>14297 rows × 12 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       plot__G11  plot__G12  plot__G13  plot__G16  plot__G17  plot__G18  \\\n",
       "0              0          0          0          0          1          0   \n",
       "1              0          0          0          0          1          0   \n",
       "2              0          0          0          0          1          0   \n",
       "3              0          0          0          0          1          0   \n",
       "4              0          0          0          0          1          0   \n",
       "...          ...        ...        ...        ...        ...        ...   \n",
       "14292          0          0          0          0          0          0   \n",
       "14293          0          0          0          0          0          0   \n",
       "14294          0          0          0          0          0          0   \n",
       "14295          0          0          0          0          0          0   \n",
       "14296          0          0          0          0          0          0   \n",
       "\n",
       "       plot__G19  plot__G20  plot__G21  plot__G6  plot__G8  plot__G9  \n",
       "0              0          0          0         0         0         0  \n",
       "1              0          0          0         0         0         0  \n",
       "2              0          0          0         0         0         0  \n",
       "3              0          0          0         0         0         0  \n",
       "4              0          0          0         0         0         0  \n",
       "...          ...        ...        ...       ...       ...       ...  \n",
       "14292          0          1          0         0         0         0  \n",
       "14293          0          1          0         0         0         0  \n",
       "14294          0          1          0         0         0         0  \n",
       "14295          0          1          0         0         0         0  \n",
       "14296          0          1          0         0         0         0  \n",
       "\n",
       "[14297 rows x 12 columns]"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.get_dummies(data2['放牧小区（plot）'],prefix='plot_')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "a8c2d22d-3927-41fe-80b1-0d36a1f35098",
   "metadata": {},
   "outputs": [],
   "source": [
    "data2=pd.concat([data2,pd.get_dummies(data2['放牧小区（plot）'],prefix='plot_')],axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "4b28d7af-49e3-45e5-bbd2-0c960e33153d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['NG', 'LGI', 'MGI', 'HGI'], dtype=object)"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2['放牧强度（intensity）'].unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "c3487aa7-b680-4c1f-815e-f73b9e188ca0",
   "metadata": {},
   "outputs": [],
   "source": [
    "#对照（NG， 0羊/天/公顷 ）、轻度放牧强度（LGI， 2羊/天/公顷 ）、中度放牧强度（MGI，4羊/天/公顷 ）和重度放牧强度（HGI，8羊/天/公顷 ）\n",
    "data2['放牧强度']=data2['放牧强度（intensity）'].map(\n",
    "{\n",
    "    'NG':0,\n",
    "    'LGI':1, \n",
    "    'MGI':2, \n",
    "    'HGI':3\n",
    "})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "b33e6091-c7d3-45cb-897a-997515ef3fa5",
   "metadata": {},
   "outputs": [],
   "source": [
    "data2=pd.concat([data2,pd.get_dummies(data2['轮次'],prefix='轮次_')],axis=1)\n",
    "\n",
    "data2=pd.concat([data2,pd.get_dummies(data2['处理'],prefix='处理_')],axis=1)\n",
    "\n",
    "data2=pd.concat([data2,pd.get_dummies(data2['植物群落功能群'],prefix='功能群_')],axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "f85c89d8-9c03-4da0-8008-acebb163d679",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['年份', '放牧小区（plot）', '放牧强度（intensity）', 'SOC土壤有机碳', 'SIC土壤无机碳',\n",
       "       'STC土壤全碳', '全氮N', '土壤C/N比', '轮次', '处理', '植物种名', '植物群落功能群', '重复', '营养苗',\n",
       "       '株/丛数', '鲜重(g)', '干重(g)', 'plot__G11', 'plot__G12', 'plot__G13',\n",
       "       'plot__G16', 'plot__G17', 'plot__G18', 'plot__G19', 'plot__G20',\n",
       "       'plot__G21', 'plot__G6', 'plot__G8', 'plot__G9', '放牧强度', '轮次__牧前',\n",
       "       '轮次__第一轮牧后', '轮次__第三轮牧后', '轮次__第二轮牧后', '轮次__第四轮牧后', '处理__中牧（6天）',\n",
       "       '处理__无牧（0天）', '处理__轻牧（3天）', '处理__重牧（12天）', '功能群__AB', '功能群__PB',\n",
       "       '功能群__PF', '功能群__PR'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data2.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "e0c4e617-471a-46a7-aa58-914c1d0935d3",
   "metadata": {},
   "outputs": [],
   "source": [
    "data2_process=data2[['年份', '放牧强度', 'SOC土壤有机碳', 'SIC土壤无机碳',\n",
    "       'STC土壤全碳', '全氮N', '土壤C/N比',  '营养苗',\n",
    "       '株/丛数', '鲜重(g)', '干重(g)', 'plot__G11', 'plot__G12', 'plot__G13',\n",
    "       'plot__G16', 'plot__G17', 'plot__G18', 'plot__G19', 'plot__G20',\n",
    "       'plot__G21', 'plot__G6', 'plot__G8', 'plot__G9', '轮次__牧前', '轮次__第一轮牧后',\n",
    "       '轮次__第三轮牧后', '轮次__第二轮牧后', '轮次__第四轮牧后', '处理__中牧（6天）', '处理__无牧（0天）',\n",
    "       '处理__轻牧（3天）', '处理__重牧（12天）', '功能群__AB', '功能群__PB', '功能群__PF',\n",
    "       '功能群__PR']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "8e870d5d-6b33-4d11-abcf-5be3582299a4",
   "metadata": {},
   "outputs": [],
   "source": [
    "df=pd.merge(tempdf,data2_process,on=['年份'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "a1ae55be-f700-4fb6-8582-b58736072285",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>月份</th>\n",
       "      <th>年份</th>\n",
       "      <th>经度(lon)</th>\n",
       "      <th>纬度(lat)</th>\n",
       "      <th>10cm湿度(kg/m2)</th>\n",
       "      <th>40cm湿度(kg/m2)</th>\n",
       "      <th>100cm湿度(kg/m2)</th>\n",
       "      <th>200cm湿度(kg/m2)</th>\n",
       "      <th>土壤蒸发量(W/m2)</th>\n",
       "      <th>土壤蒸发量(mm)</th>\n",
       "      <th>...</th>\n",
       "      <th>轮次__第二轮牧后</th>\n",
       "      <th>轮次__第四轮牧后</th>\n",
       "      <th>处理__中牧（6天）</th>\n",
       "      <th>处理__无牧（0天）</th>\n",
       "      <th>处理__轻牧（3天）</th>\n",
       "      <th>处理__重牧（12天）</th>\n",
       "      <th>功能群__AB</th>\n",
       "      <th>功能群__PB</th>\n",
       "      <th>功能群__PF</th>\n",
       "      <th>功能群__PR</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2020</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.21</td>\n",
       "      <td>33.85</td>\n",
       "      <td>47.07</td>\n",
       "      <td>166.7</td>\n",
       "      <td>0.33</td>\n",
       "      <td>0.36</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>2020</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.21</td>\n",
       "      <td>33.85</td>\n",
       "      <td>47.07</td>\n",
       "      <td>166.7</td>\n",
       "      <td>0.33</td>\n",
       "      <td>0.36</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1</td>\n",
       "      <td>2020</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.21</td>\n",
       "      <td>33.85</td>\n",
       "      <td>47.07</td>\n",
       "      <td>166.7</td>\n",
       "      <td>0.33</td>\n",
       "      <td>0.36</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>2020</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.21</td>\n",
       "      <td>33.85</td>\n",
       "      <td>47.07</td>\n",
       "      <td>166.7</td>\n",
       "      <td>0.33</td>\n",
       "      <td>0.36</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1</td>\n",
       "      <td>2020</td>\n",
       "      <td>115.375</td>\n",
       "      <td>44.125</td>\n",
       "      <td>12.21</td>\n",
       "      <td>33.85</td>\n",
       "      <td>47.07</td>\n",
       "      <td>166.7</td>\n",
       "      <td>0.33</td>\n",
       "      <td>0.36</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 48 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "   月份    年份  经度(lon)  纬度(lat)  10cm湿度(kg/m2)  40cm湿度(kg/m2)  100cm湿度(kg/m2)  \\\n",
       "0   1  2020  115.375   44.125          12.21          33.85           47.07   \n",
       "1   1  2020  115.375   44.125          12.21          33.85           47.07   \n",
       "2   1  2020  115.375   44.125          12.21          33.85           47.07   \n",
       "3   1  2020  115.375   44.125          12.21          33.85           47.07   \n",
       "4   1  2020  115.375   44.125          12.21          33.85           47.07   \n",
       "\n",
       "   200cm湿度(kg/m2)  土壤蒸发量(W/m2)  土壤蒸发量(mm)  ...  轮次__第二轮牧后  轮次__第四轮牧后  \\\n",
       "0           166.7         0.33       0.36  ...          0          0   \n",
       "1           166.7         0.33       0.36  ...          0          0   \n",
       "2           166.7         0.33       0.36  ...          0          0   \n",
       "3           166.7         0.33       0.36  ...          0          0   \n",
       "4           166.7         0.33       0.36  ...          0          0   \n",
       "\n",
       "   处理__中牧（6天）  处理__无牧（0天）  处理__轻牧（3天）  处理__重牧（12天）  功能群__AB  功能群__PB  功能群__PF  \\\n",
       "0           0           1           0            0        0        1        0   \n",
       "1           0           1           0            0        0        1        0   \n",
       "2           0           1           0            0        0        0        0   \n",
       "3           0           1           0            0        0        0        0   \n",
       "4           0           1           0            0        0        0        1   \n",
       "\n",
       "   功能群__PR  \n",
       "0        0  \n",
       "1        0  \n",
       "2        1  \n",
       "3        1  \n",
       "4        0  \n",
       "\n",
       "[5 rows x 48 columns]"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "id": "91069e90-f66c-4f62-a622-9d78847ef06f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(171564, 48)"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "b881db73-c817-492b-8355-7b1dda2baf57",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "月份                0.000000\n",
       "年份                0.000000\n",
       "经度(lon)           0.000000\n",
       "纬度(lat)           0.000000\n",
       "10cm湿度(kg/m2)     0.000000\n",
       "40cm湿度(kg/m2)     0.000000\n",
       "100cm湿度(kg/m2)    0.000000\n",
       "200cm湿度(kg/m2)    0.000000\n",
       "土壤蒸发量(W/m2)       0.000000\n",
       "土壤蒸发量(mm)         0.000000\n",
       "植被指数(NDVI)        0.000000\n",
       "径流量(m3/s)         0.000000\n",
       "径流量(m3)           0.000000\n",
       "放牧强度              0.000000\n",
       "SOC土壤有机碳          0.000000\n",
       "SIC土壤无机碳          0.000000\n",
       "STC土壤全碳           0.000000\n",
       "全氮N               0.000000\n",
       "土壤C/N比            0.000000\n",
       "营养苗               0.001679\n",
       "株/丛数              0.024970\n",
       "鲜重(g)             0.049311\n",
       "干重(g)             0.000000\n",
       "plot__G11         0.000000\n",
       "plot__G12         0.000000\n",
       "plot__G13         0.000000\n",
       "plot__G16         0.000000\n",
       "plot__G17         0.000000\n",
       "plot__G18         0.000000\n",
       "plot__G19         0.000000\n",
       "plot__G20         0.000000\n",
       "plot__G21         0.000000\n",
       "plot__G6          0.000000\n",
       "plot__G8          0.000000\n",
       "plot__G9          0.000000\n",
       "轮次__牧前            0.000000\n",
       "轮次__第一轮牧后         0.000000\n",
       "轮次__第三轮牧后         0.000000\n",
       "轮次__第二轮牧后         0.000000\n",
       "轮次__第四轮牧后         0.000000\n",
       "处理__中牧（6天）        0.000000\n",
       "处理__无牧（0天）        0.000000\n",
       "处理__轻牧（3天）        0.000000\n",
       "处理__重牧（12天）       0.000000\n",
       "功能群__AB           0.000000\n",
       "功能群__PB           0.000000\n",
       "功能群__PF           0.000000\n",
       "功能群__PR           0.000000\n",
       "dtype: float64"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.isnull().sum()/df.shape[0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "2b74c99f-2554-493a-b5c4-ba85b87a9846",
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in ['营养苗','株/丛数', '鲜重(g)']:\n",
    "    df[i]=df[i].fillna(df[i].mean())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "id": "9f00aa79-bc8d-4cc2-8abf-7fd2cc5cdf8f",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "def filter_extreme_MAD(dataframe,n=3):\n",
    "    # MAD:绝对中位差去极值\n",
    "    for i in dataframe.columns:\n",
    "        median = dataframe[i].quantile(0.5)\n",
    "        new_median = ((dataframe[i] - median).abs()).quantile(0.50)\n",
    "        max_range = median + n * new_median\n",
    "        min_range = median - n * new_median\n",
    "        dataframe[i]=pd.DataFrame(np.clip(dataframe[i].values, min_range, max_range),columns=None)\n",
    "\n",
    "    return dataframe\n",
    "def filter_extreme_3sigma(dataframe,n=3):\n",
    "    # 3sigma 去极值\n",
    "    for i in dataframe.columns:\n",
    "        mean=dataframe[i].mean()\n",
    "        std=dataframe[i].std()\n",
    "        max_range=mean+n*std\n",
    "        min_range=mean-n*std\n",
    "        dataframe[i] = pd.DataFrame(np.clip(dataframe[i].values, min_range, max_range), columns=None)\n",
    "    return dataframe\n",
    "\n",
    "def filter_extreme_percentile(dataframe,min=0.025,max=0.975):\n",
    "    # 百分位法 去极值\n",
    "    for i in dataframe.columns:\n",
    "        Temp=dataframe[i].sort_values()\n",
    "        q=Temp.quantile([min,max])\n",
    "        dataframe[i] = pd.DataFrame(np.clip(dataframe[i].values, q.iloc[0], q.iloc[1]), columns=None)\n",
    "    return"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "id": "5bc855b1-7704-4e8e-accd-1dc60e374ba2",
   "metadata": {},
   "outputs": [],
   "source": [
    "df=filter_extreme_3sigma(df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "068c99f2-1b2e-4701-9a6d-901c47dd2d3e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "月份                0\n",
       "年份                0\n",
       "经度(lon)           0\n",
       "纬度(lat)           0\n",
       "10cm湿度(kg/m2)     0\n",
       "40cm湿度(kg/m2)     0\n",
       "100cm湿度(kg/m2)    0\n",
       "200cm湿度(kg/m2)    0\n",
       "土壤蒸发量(W/m2)       0\n",
       "土壤蒸发量(mm)         0\n",
       "植被指数(NDVI)        0\n",
       "径流量(m3/s)         0\n",
       "径流量(m3)           0\n",
       "放牧强度              0\n",
       "SOC土壤有机碳          0\n",
       "SIC土壤无机碳          0\n",
       "STC土壤全碳           0\n",
       "全氮N               0\n",
       "土壤C/N比            0\n",
       "营养苗               0\n",
       "株/丛数              0\n",
       "鲜重(g)             0\n",
       "干重(g)             0\n",
       "plot__G11         0\n",
       "plot__G12         0\n",
       "plot__G13         0\n",
       "plot__G16         0\n",
       "plot__G17         0\n",
       "plot__G18         0\n",
       "plot__G19         0\n",
       "plot__G20         0\n",
       "plot__G21         0\n",
       "plot__G6          0\n",
       "plot__G8          0\n",
       "plot__G9          0\n",
       "轮次__牧前            0\n",
       "轮次__第一轮牧后         0\n",
       "轮次__第三轮牧后         0\n",
       "轮次__第二轮牧后         0\n",
       "轮次__第四轮牧后         0\n",
       "处理__中牧（6天）        0\n",
       "处理__无牧（0天）        0\n",
       "处理__轻牧（3天）        0\n",
       "处理__重牧（12天）       0\n",
       "功能群__AB           0\n",
       "功能群__PB           0\n",
       "功能群__PF           0\n",
       "功能群__PR           0\n",
       "dtype: int64"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.isnull().sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f94f9e46-2525-49ff-8246-b05877f62f49",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "id": "32ab8d1b-0383-45a5-b8d2-e0febb94b6f9",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Looking in indexes: https://pypi.tuna.tsinghua.edu.cn/simple\n",
      "Requirement already satisfied: xgboost in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (1.5.2)\n",
      "Requirement already satisfied: numpy in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from xgboost) (1.22.3)\n",
      "Requirement already satisfied: scipy in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from xgboost) (1.8.1)\n",
      "Looking in indexes: https://pypi.tuna.tsinghua.edu.cn/simple\n",
      "Requirement already satisfied: lightgbm in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (3.3.2)\n",
      "Requirement already satisfied: scipy in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from lightgbm) (1.8.1)\n",
      "Requirement already satisfied: scikit-learn!=0.22.0 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from lightgbm) (0.24.2)\n",
      "Requirement already satisfied: wheel in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from lightgbm) (0.37.1)\n",
      "Requirement already satisfied: numpy in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from lightgbm) (1.22.3)\n",
      "Requirement already satisfied: threadpoolctl>=2.0.0 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from scikit-learn!=0.22.0->lightgbm) (3.1.0)\n",
      "Requirement already satisfied: joblib>=0.11 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from scikit-learn!=0.22.0->lightgbm) (1.1.0)\n",
      "Looking in indexes: https://pypi.tuna.tsinghua.edu.cn/simple\n",
      "Requirement already satisfied: matplotlib in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (3.5.2)\n",
      "Requirement already satisfied: python-dateutil>=2.7 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (2.8.2)\n",
      "Requirement already satisfied: packaging>=20.0 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (21.3)\n",
      "Requirement already satisfied: kiwisolver>=1.0.1 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (1.4.2)\n",
      "Requirement already satisfied: pillow>=6.2.0 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (9.1.0)\n",
      "Requirement already satisfied: cycler>=0.10 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (0.11.0)\n",
      "Requirement already satisfied: pyparsing>=2.2.1 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (3.0.6)\n",
      "Requirement already satisfied: numpy>=1.17 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (1.22.3)\n",
      "Requirement already satisfied: fonttools>=4.22.0 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from matplotlib) (4.33.3)\n",
      "Requirement already satisfied: six>=1.5 in c:\\users\\administrator\\appdata\\roaming\\spsspro\\spsspro\\lib\\python\\lib\\site-packages (from python-dateutil>=2.7->matplotlib) (1.16.0)\n"
     ]
    }
   ],
   "source": [
    "import sys\n",
    "!{sys.executable} -m pip install xgboost -i https://pypi.tuna.tsinghua.edu.cn/simple\n",
    "!{sys.executable} -m pip install lightgbm -i https://pypi.tuna.tsinghua.edu.cn/simple\n",
    "!{sys.executable} -m pip install matplotlib -i https://pypi.tuna.tsinghua.edu.cn/simple"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "id": "8835dec6-b239-41e2-815e-9a4a65c540c1",
   "metadata": {},
   "outputs": [],
   "source": [
    "import xgboost as xgb\n",
    "import lightgbm as lgb\n",
    "from sklearn.ensemble import RandomForestClassifier\n",
    "from sklearn.linear_model import LogisticRegression\n",
    "from sklearn.metrics import classification_report,f1_score\n",
    "from sklearn.model_selection import train_test_split"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "id": "020ebdf9-131c-4fc3-ac14-247c74590799",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "逻辑 回归\n",
      "              precision    recall  f1-score   support\n",
      "\n",
      "         0.0       0.29      0.50      0.37      6498\n",
      "         1.0       0.24      0.24      0.24     14986\n",
      "         2.0       0.49      0.45      0.47     13947\n",
      "         3.0       0.47      0.38      0.42     16039\n",
      "\n",
      "    accuracy                           0.37     51470\n",
      "   macro avg       0.37      0.39      0.37     51470\n",
      "weighted avg       0.39      0.37      0.37     51470\n",
      "\n",
      "Randomforest分类\n",
      "              precision    recall  f1-score   support\n",
      "\n",
      "         0.0       1.00      1.00      1.00     11085\n",
      "         1.0       1.00      1.00      1.00     14430\n",
      "         2.0       1.00      1.00      1.00     12893\n",
      "         3.0       1.00      1.00      1.00     13062\n",
      "\n",
      "    accuracy                           1.00     51470\n",
      "   macro avg       1.00      1.00      1.00     51470\n",
      "weighted avg       1.00      1.00      1.00     51470\n",
      "\n",
      "lgbm分类\n",
      "              precision    recall  f1-score   support\n",
      "\n",
      "         0.0       1.00      1.00      1.00     11085\n",
      "         1.0       1.00      1.00      1.00     14430\n",
      "         2.0       1.00      1.00      1.00     12893\n",
      "         3.0       1.00      1.00      1.00     13062\n",
      "\n",
      "    accuracy                           1.00     51470\n",
      "   macro avg       1.00      1.00      1.00     51470\n",
      "weighted avg       1.00      1.00      1.00     51470\n",
      "\n",
      "[18:47:32] WARNING: C:/Users/Administrator/workspace/xgboost-win64_release_1.5.1/src/learner.cc:1115: Starting in XGBoost 1.3.0, the default evaluation metric used with the objective 'multi:softprob' was changed from 'merror' to 'mlogloss'. Explicitly set eval_metric if you'd like to restore the old behavior.\n",
      "XGboost分类\n",
      "              precision    recall  f1-score   support\n",
      "\n",
      "         0.0       1.00      1.00      1.00     11085\n",
      "         1.0       1.00      1.00      1.00     14430\n",
      "         2.0       1.00      1.00      1.00     12893\n",
      "         3.0       1.00      1.00      1.00     13062\n",
      "\n",
      "    accuracy                           1.00     51470\n",
      "   macro avg       1.00      1.00      1.00     51470\n",
      "weighted avg       1.00      1.00      1.00     51470\n",
      "\n"
     ]
    }
   ],
   "source": [
    "\n",
    "# X=df[['SOC土壤有机碳', 'SIC土壤无机碳',\n",
    "#        'STC土壤全碳', '全氮N', '土壤C/N比',  '营养苗',\n",
    "#        '株/丛数', '鲜重(g)', '干重(g)', 'plot__G11', 'plot__G12', 'plot__G13',\n",
    "#        'plot__G16', 'plot__G17', 'plot__G18', 'plot__G19', 'plot__G20',\n",
    "#        'plot__G21', 'plot__G6', 'plot__G8', 'plot__G9', '轮次__牧前', '轮次__第一轮牧后',\n",
    "#        '轮次__第三轮牧后', '轮次__第二轮牧后', '轮次__第四轮牧后', '处理__中牧（6天）', '处理__无牧（0天）',\n",
    "#        '处理__轻牧（3天）', '处理__重牧（12天）', '功能群__AB', '功能群__PB', '功能群__PF',\n",
    "#        '功能群__PR']]\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "X=df[['SOC土壤有机碳', 'SIC土壤无机碳',\n",
    "       'STC土壤全碳', '全氮N', '土壤C/N比']]\n",
    "\n",
    "y=df['放牧强度']\n",
    "# 拆分训练集和测试集\n",
    "x_train, x_test, y_train, y_test = train_test_split(X, y, test_size=0.3, random_state=12)\n",
    "    \n",
    "#逻辑 回归\n",
    "model = LogisticRegression()\n",
    "\n",
    "model.fit(x_train, y_train)\n",
    "print('逻辑 回归')\n",
    "print(classification_report(model.predict(x_test),y_test))\n",
    "\n",
    "\n",
    "# Randomforest分类\n",
    "\n",
    "model = RandomForestClassifier()\n",
    "\n",
    "model.fit(x_train, y_train)\n",
    "print('Randomforest分类')\n",
    "print(classification_report(model.predict(x_test),y_test))\n",
    "\n",
    "# lgbm分类\n",
    "model = lgb.LGBMClassifier()\n",
    "model.fit(x_train, y_train)\n",
    "print('lgbm分类')\n",
    "print(classification_report(model.predict(x_test),y_test))\n",
    "\n",
    "\n",
    "# XGboost分类\n",
    "model = xgb.XGBClassifier()\n",
    "model.fit(x_train, y_train)\n",
    "print('XGboost分类')\n",
    "print(classification_report(model.predict(x_test),y_test))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "id": "03352e36-e1b1-4d48-8508-e285d3f24c81",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "逻辑 回归\n",
      "              precision    recall  f1-score   support\n",
      "\n",
      "         0.0       0.24      0.35      0.29     25476\n",
      "         1.0       0.26      0.25      0.25     48816\n",
      "         2.0       0.60      0.49      0.54     53460\n",
      "         3.0       0.41      0.41      0.41     43812\n",
      "\n",
      "    accuracy                           0.38    171564\n",
      "   macro avg       0.38      0.37      0.37    171564\n",
      "weighted avg       0.40      0.38      0.39    171564\n",
      "\n"
     ]
    }
   ],
   "source": [
    "#逻辑 回归\n",
    "model = LogisticRegression()\n",
    "\n",
    "model.fit(X, y)\n",
    "print('逻辑 回归')\n",
    "print(classification_report(model.predict(X),y))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "id": "c6c12cb6-bed4-4120-a02a-ed0e87a9d8ee",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ 0.18478838, -0.10679748,  0.1906822 , -1.88945995,  0.02358119],\n",
       "       [ 0.00475829,  0.02593296,  0.11509955, -0.58182518, -0.10627045],\n",
       "       [-0.03134059,  0.02497467, -0.2447123 ,  0.908844  ,  0.12710039],\n",
       "       [-0.15820608,  0.05588984, -0.06106944,  1.56244114, -0.04441114]])"
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.coef_"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "id": "c796fbb7-82c1-44d2-b11d-f965f86b9b79",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-3.16704975, -0.19630133,  2.32242855,  1.04092253])"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.intercept_"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5e4f8234-b911-4e92-ba3c-e453318fd9f3",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4717a7ee-1ed0-47a2-9be1-097b27da5851",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (Bundle)",
   "language": "python",
   "name": "pyb"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}