孫洪兵　蔣舜媛　孫輝　王紅蘭　周毅　朱文濤

[摘要]為解決藥材生長最適宜生長區(qū)域與高品質(zhì)區(qū)域空間分離的問題，建立具有針對性生產(chǎn)指導(dǎo)的區(qū)劃，該文以羌活藥材為例提出綜合考慮藥材生長適宜性和品質(zhì)適宜性的中藥材功能型生產(chǎn)區(qū)劃方法?；谇捌谇蓟钌L適宜性和品質(zhì)適宜性區(qū)劃的研究結(jié)果，將羌活生長適宜性空間矢量數(shù)據(jù)與藥材品質(zhì)適宜性空間矢量數(shù)據(jù)進行空間加權(quán)總和處理，獲得羌活生產(chǎn)適宜性空間矢量數(shù)據(jù)，發(fā)現(xiàn)優(yōu)質(zhì)高產(chǎn)區(qū)域主要分布在道地產(chǎn)區(qū)四川省境內(nèi)；利用研究區(qū)域的土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)進一步對羌活生產(chǎn)適宜區(qū)進行功能劃分，發(fā)現(xiàn)適宜羌活規(guī)模化種植和野生林地撫育區(qū)主要分布在四川阿壩州和甘孜州境內(nèi)，與實地調(diào)研結(jié)果一致。該研究結(jié)果將為羌活的道地性研究、野生資源保護、產(chǎn)業(yè)化種植和野生林地撫育提供科學(xué)的理論依據(jù)與技術(shù)支持。

[關(guān)鍵詞]羌活； 3S； 生長適宜性； 品質(zhì)適宜性； 生產(chǎn)適宜性； 功能型區(qū)劃

[Abstract]Based on the research results of suitability evaluation of growth and production quality， a functional production cultivation regionalization with high feasibility and operability to solve the problem of the separation of high quality and high yield were formatted for protection， wild monitoring， and cultivation of this plant by weighted sum of spatial suitability data of ecology and quality， as well as integrated with land use and cover data The results of the study revealed that good quality and high yield area were mainly distributed in the original production areas in Sichuan province and where could carry out monitoring and commercialization cultivation This method is expected to overcome shortage of traditional regionalization methods difficult to distinguish the contribution of ecological factors and quality factors， which provide an innovative theory and methodology for regionalization， and is helpful to practical application of wild resource protection， nature research， monitoring， and commercialization cultivation for Chinese material medica

[Key words]Notopterygium incisum； 3S； growth suitability； quality suitability； production quality； functional regionalization

中藥區(qū)劃主要是研究中藥材品質(zhì)及產(chǎn)量空間屬性及其地域系統(tǒng)的空間分異規(guī)律，最終綜合社會、經(jīng)濟及人為因素進行區(qū)域劃分，其研究成果可為中藥材的生產(chǎn)提供基礎(chǔ)數(shù)據(jù)和技術(shù)支撐[1]。由于區(qū)域之間中藥材所在地的自然生態(tài)和人為因素等“立地條件”不同，中藥材生長及品質(zhì)適宜性的空間分異規(guī)律或許不同[2]。即植物積累次生代謝產(chǎn)物所需的適宜生境（藥材品質(zhì)適宜性）與其生長發(fā)育的適宜生境（植物生長適宜性）可能并不一致，甚至相反，即藥材的高產(chǎn)中心有可能與高質(zhì)中心分離或者偏離[3]。因此，在利用中藥材生產(chǎn)區(qū)劃技術(shù)對所謂的“立地條件”進行研究的過程中，不但要體現(xiàn)原植物的生長適宜性，即完成生長適宜性區(qū)劃，還應(yīng)從藥材品質(zhì)（次生代謝產(chǎn)物的積累）與區(qū)域相關(guān)的角度完成藥材品質(zhì)適宜性區(qū)劃，然后選取適宜的空間疊加分析模型，對二者進行空間疊加，才能獲得較為合理的藥材生產(chǎn)區(qū)劃結(jié)果。這不是單純的生態(tài)（生境）適宜性區(qū)劃所能完成的任務(wù)。

為解決上述問題，本文以羌活藥材為例，在前期羌活生長適宜性[4]和品質(zhì)適宜性區(qū)劃研究結(jié)果[5]的基礎(chǔ)上，通過分析兩者空間分異規(guī)律，建立更加具有產(chǎn)業(yè)化實踐指導(dǎo)意義的中藥材功能型生產(chǎn)區(qū)劃的評價體系和方法，以期為羌活的資源保護、野生撫育和產(chǎn)業(yè)化種植實踐指導(dǎo)，也為其他中藥材產(chǎn)業(yè)規(guī)劃提供科學(xué)理論和參考依據(jù)。

1材料

11樣點數(shù)據(jù)于2013年8月下旬至9月中旬和2014年8月下旬至9月中旬，在四川、甘肅、青海、西藏等省采集野生羌活樣品，共設(shè)定38個采樣點（采樣點定位信息經(jīng)由GPS實地獲得，并經(jīng)ArcGIS100軟件進行圖層篩選得到32個采樣點），每個采樣點不少于5株，收集了217份樣本[4]。經(jīng)四川省中醫(yī)藥科學(xué)院蔣舜媛研究員和周毅副研究員鑒定為傘形科植物羌活（Notopterygium incisum，憑證標本保存于四川省中醫(yī)藥科學(xué)院）。樣品按單株陰干、粉碎，過30目篩后測定其有效成分含量。利用超高效液相色譜法（UPLCPDA）測定各樣品中羌活醇（notopterol）和異歐前胡素（isoimperatorin）的含量[6]。

12地理及環(huán)境數(shù)據(jù)生態(tài)因子數(shù)據(jù)來源于“中藥資源空間信息網(wǎng)格數(shù)據(jù)庫”（http：//wwwtcmresourcescom/）；縣級以上行政邊界數(shù)據(jù)下載于國家基礎(chǔ)地理信息系統(tǒng)網(wǎng)（http：//nfgisnsdigovcn/），比例尺為1∶400萬；2010年中國土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：//wwwresdccn）。

2方法

21羌活生產(chǎn)適宜性劃分及分析基于前期羌活生長適宜性和品質(zhì)適宜性區(qū)劃的研究結(jié)果[4]，將羌活生長適宜性和品質(zhì)適宜性空間分布矢量數(shù)據(jù)進行空間關(guān)聯(lián)性分析，分析兩者空間分異規(guī)律，在此基礎(chǔ)上選取適宜的疊加模型，獲得綜合藥材品質(zhì)與生長適宜性的羌活生產(chǎn)適宜性空間分布數(shù)據(jù)，并進行區(qū)劃分析。

22羌活功能型區(qū)劃構(gòu)建及分析借助ArcGIS強大數(shù)據(jù)處理平臺，從2010年中國土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)中提取研究區(qū)域農(nóng)用耕地遙感監(jiān)測數(shù)據(jù)和林地遙感監(jiān)測數(shù)據(jù)，經(jīng)重分類后，以此為掩膜分別對羌活生產(chǎn)適宜性空間分布矢量數(shù)據(jù)進行空間數(shù)據(jù)提取，建立羌活規(guī)?；斯しN植區(qū)劃和羌活野生林地撫育區(qū)域。最后利用ArcGIS中重分類處理方法對各適宜分布區(qū)面積進行分省統(tǒng)計。

3結(jié)果與分析

31羌活品質(zhì)及生長適宜性空間差異性分析借助ArcGIS數(shù)據(jù)處理平臺的Spatial Analyst多元分析（multivariate）中的波段集統(tǒng)計（band collection statistics）模塊，選用計算協(xié)方差和相關(guān)矩陣選項，利用式（1）公式求取羌活生長適宜性圖層和品質(zhì)適宜性圖層之間的協(xié)方差。

Covij=Nk=1（Zij-μi）（Zij-μj）N-1（1）

式中，Z為像元值；μ為圖層平均值；N為像元的數(shù)量；k為特定像元；i代指羌活生長適宜性圖層；j代指羌活品質(zhì)適宜性圖層。

此基礎(chǔ)上利用式（2）將羌活生長適宜性圖層和品質(zhì)適宜性圖層間的協(xié)方差與標準差乘積進行空間運算，空間相關(guān)系數(shù)Corrij為0449 5，呈正相關(guān)。表明羌活生長適宜性與品質(zhì)適宜性空間分布數(shù)據(jù)從整體上來講在空間變化趨勢上具有一定的相似性。

Corrij=Covijδiδj（2）

按羌活生長適宜性與品質(zhì)適宜性空間分布數(shù)據(jù)的屬性值（value）分別提取（利用ArcGIS Spatial Analyst中提取分析模塊的按屬性提取功能）各自的“高、中、低”區(qū)（采用自然間斷點分級法劃分）空間矢量數(shù)據(jù)，并兩兩做空間相關(guān)分析，分析結(jié)果表明各區(qū)段的空間相關(guān)性極低（幾乎為零），也就是說羌活植物積累次生代謝產(chǎn)物所需的適宜生境（藥材品質(zhì)適宜性）與其生長發(fā)育的適宜生境（植物生長適宜性）存在空間差異性，即羌活藥材的高產(chǎn)中心與高質(zhì)中心分離或者偏離，因此進行羌活生產(chǎn)區(qū)劃研究不能單純以生長適宜性或者簡單疊加品質(zhì)適宜性后的數(shù)據(jù)進行分析和評價，需要綜合考慮羌活的生長適宜性與品質(zhì)適宜性的各自比重，建立綜合藥材產(chǎn)量和品質(zhì)最大化的羌活生產(chǎn)區(qū)劃。

32羌活生產(chǎn)適宜性區(qū)劃模型構(gòu)建及適宜性分析目前，常見的空間疊加分析方法主要有加權(quán)疊加（weighted overlay）、加權(quán)總和（weighted sum）以及模糊疊加（fuzzy overlay）這3種，每種方法都有不同的基本前提和假設(shè)。模糊疊加以集合論為基礎(chǔ)，分析時將數(shù)據(jù)值重分類或變換到相同等級，變換后的值可用于定義隸屬于指定集合的可能性。由于模糊疊加法中變換后的值表示隸屬可能性，因此輸入柵格不需要進行加權(quán)，將轉(zhuǎn)換后的數(shù)據(jù)進行簡單的柵格運算處理（常用的算子有取平均值、最大值、最小值、加減、乘除、交集以及并集等）。疊加后的數(shù)值表示隸屬于羌活生長適宜性或品質(zhì)適宜性的可能性[7]，因此模糊疊加模型的疊加值不能用于羌活生產(chǎn)適宜性分析。

相對而言，加權(quán)疊加法與加權(quán)總和法隨著分析假設(shè)有利因素的增多，其最終輸出柵格中生成的值越高，因此這些值較高位置將確定為最佳位置，較適用于進行生產(chǎn)適宜性模型構(gòu)建。加權(quán)疊加（weighted overlay）使用常用測量比例來疊加多個柵格數(shù)據(jù)，并根據(jù)各柵格數(shù)據(jù)的重要性分配權(quán)重，分配給各輸入柵格的權(quán)重之和必須等于100%，將圖層乘以相應(yīng)的倍數(shù)，然后將每個像元生成的值相加。而加權(quán)總和（weighted sum）是通過將柵格各自乘以指定的權(quán)重并合計在一起來疊加多個柵格。與“加權(quán)疊加”分析有所不同，其使用的權(quán)重不同，而且這些值不會按照定義的等級重新安排。因此，模型中所輸入值的屬性分辨率將保持不變。這2種工具主要區(qū)別如下：加權(quán)總和工具不能將重分類值重設(shè)為評估等級，因此可保持其分辨率。加權(quán)總和工具允許使用浮點型和整型值作為輸入，而加權(quán)疊加工具只接受整型柵格作為輸入。加權(quán)疊加工具最常用于適宜性建模，并且可用于確保遵循正確的方法。加權(quán)總和工具在需要保持模型分辨率，或需要浮點型輸出或小數(shù)權(quán)重時很有用。加權(quán)總和可將每個輸入柵格的指定字段值與指定權(quán)重相乘。然后將所有輸入柵格相加來創(chuàng)建輸出柵格。因為羌活生長適宜性和品質(zhì)適宜性空間矢量數(shù)據(jù)為浮點型數(shù)據(jù)，而且疊加結(jié)果需要保持模型的分辨率，以便于進行產(chǎn)地適宜性分析，不宜在數(shù)據(jù)運算中將變量進行重分類作為評估等級。因此，最終選取加權(quán)總和模型（weighted sum）進行羌活生產(chǎn)適宜性區(qū)劃研究。

由加權(quán)總和的疊加結(jié)果可知，羌活生產(chǎn)適宜性高度區(qū)段主要分布在四川省阿壩州境內(nèi)，西藏的昌都地區(qū)及甘肅省的甘南地區(qū)有零星分布；中度適宜區(qū)段主要分布在四川省阿壩州以及甘孜州境內(nèi)，以甘孜州為主；低度適宜區(qū)主要分布在青海省東部地區(qū)、甘肅隴南地區(qū)、西藏昌都等地區(qū)以及三?。ùā⑶嗪筒兀┙唤缣帲▓D1）。

33羌活功能型區(qū)劃研究中藥材生產(chǎn)區(qū)劃研究的目的是為中藥材的生產(chǎn)做科學(xué)性指導(dǎo)，目前中藥材生產(chǎn)方式主要有人工栽培和野生撫育（仿野生栽培）2種形式。為建立更具指導(dǎo)性、科學(xué)性和實用性的中藥材生產(chǎn)區(qū)劃，需要按照其所需的土地利用類型進行功能性劃分。以羌活為例，其野生種群的生境主要是亞高山針葉林上沿、林線、林窗、樹線交錯帶及高山灌叢草地，生境特征主要是林地（前期羌活生理生態(tài)的相關(guān)研究也表明羌活是森林成分而非草地成分），因此野生撫育區(qū)域應(yīng)該以林地為主，灌叢草地因為及生境脆弱和難以到達性而難以成為野生撫育的區(qū)域；羌活人工栽培主要是在高寒山區(qū)耕地上進行，高寒區(qū)域因生態(tài)保護而嚴禁墾殖林地和草地，因此人工栽培區(qū)域的土地利用方式主要是耕地。以此為依據(jù)，羌活的生產(chǎn)區(qū)劃按功能性可劃分為適宜羌活野生撫育的林地生產(chǎn)區(qū)劃和適宜人工栽培的耕地生產(chǎn)區(qū)劃。

羌活功能型區(qū)劃構(gòu)建中首先借助ArcGIS平臺按屬性值從2010年土地利用遙感檢測矢量數(shù)據(jù)中分別提取耕地和林地分布數(shù)據(jù)，并進行重分類，導(dǎo)出SHP格式矢量數(shù)據(jù)，然后按此掩膜提取生產(chǎn)適宜性分布數(shù)據(jù)；采用自然間斷點分級法（Jenks）進行等級劃分，分別建立適合羌活規(guī)?；斯しN植的高寒區(qū)域的農(nóng)用耕地（圖2a）和野生撫育的高寒林地區(qū)劃分布圖（圖2b）。

由圖2a中可知，適宜羌活生產(chǎn)的農(nóng)耕用地分為3個區(qū)段：高適宜區(qū)段、中適宜區(qū)段和低適宜區(qū)段。高適宜區(qū)主要分布在四川阿壩州和甘孜州，該區(qū)域因為環(huán)境的空間異質(zhì)性高而且高海拔適宜區(qū)域的耕地面積相對較大，因而有一定面積的適宜耕地；中適宜區(qū)分布在甘肅省隴南地區(qū)、青海省東部地區(qū)以及兩省交界處，該區(qū)域盡管海拔普遍較甘孜阿壩低但因為緯度偏北而且耕地面積較大，因而有較大面積的適宜耕地（甘肅省和青海?。?；低適宜區(qū)分布在西藏昌都地區(qū)以及與青海省交界處，這部分面積較小，一是因為海拔高、寒旱問題嚴重、土壤質(zhì)量不高，二是因為該區(qū)域以草地為主而耕地面積非常有限。

適合羌活野生撫育林地（圖2b）也大致劃分為3個區(qū)段：高適宜區(qū)段、中適宜區(qū)段和低適宜區(qū)段。高適宜區(qū)段主要分布在四川阿壩州境內(nèi)；中適宜區(qū)段分布在青海以及四川甘孜州境內(nèi)；低適宜區(qū)段主要分布在甘肅省隴南地區(qū)和西藏省的昌都地區(qū)。

利用ArcGIS空間統(tǒng)計功能對其分布區(qū)域進行縣域統(tǒng)計（表1），結(jié)合圖2可知，四川省、西藏省和云南省適合羌活野生撫育種植區(qū)劃（林地）面積高于規(guī)模化人工種植區(qū)劃（耕地）面積。據(jù)考查[8]，該區(qū)域多為羌活N incisum的歷史產(chǎn)區(qū)，種質(zhì)資源較佳，因此可在該區(qū)域大力扶持羌活野生撫育種植，以保障藥材質(zhì)量，并提高當?shù)剞r(nóng)民收入。甘肅省和青海省適宜羌活種植的耕地面積較為寬廣，交通相對便利，可合理扶持羌活規(guī)?；斯しN植。

4討論

本文通過生長適宜性與品質(zhì)適宜性的空間疊加，構(gòu)建中藥材功能型生產(chǎn)區(qū)劃的評價體系及方法，具有實踐的針對性、結(jié)果的定量化、與總體規(guī)劃耦合等優(yōu)點，特別是對于那些道地性明顯與產(chǎn)量適宜區(qū)域分離的中藥材品種。如地黃、白前、中華栝樓、羅漢果和川貝母[9]等，很多中藥材生產(chǎn)與農(nóng)作物的生產(chǎn)不同，中藥材質(zhì)量主要體現(xiàn)在藥效，其表現(xiàn)形式是某些活性次生代謝產(chǎn)物的含量高低，適宜生長并不表明能夠同時促進次生產(chǎn)物的合成代謝旺盛，很多情況下恰恰相反，某些不利的環(huán)境條件才是一些次生產(chǎn)物（如一些萜類、揮發(fā)油、單寧類物質(zhì)等）

合成的必要條件[3]。同時，作為栽培，也不能只考慮藥材質(zhì)量而忽視經(jīng)濟可行性，即在保證質(zhì)量的前提下需要一定的經(jīng)濟產(chǎn)量，具有比較經(jīng)濟效益，才能實現(xiàn)產(chǎn)業(yè)化栽培，因此需要在生長適宜性方面進行平衡。一般而言，在區(qū)劃實際中應(yīng)充分考慮藥材生長的次適宜區(qū)，這些區(qū)域藥材有一定的生物產(chǎn)量和藥材產(chǎn)量，如果同時具有促進特定活性成分次生產(chǎn)物形成的脅迫因子，則是產(chǎn)業(yè)化栽培的潛在區(qū)域。

應(yīng)當看到，所有基于統(tǒng)計學(xué)的區(qū)劃方法都是經(jīng)驗方法，而不是從植物代謝等機制方面去研究藥材的區(qū)劃問題，在一定空間范圍內(nèi)，取樣調(diào)查的樣點數(shù)據(jù)和樣本數(shù)據(jù)越多，區(qū)劃的準確性也就越高。因此，區(qū)劃成果只是源于已有數(shù)據(jù)的統(tǒng)計模型，而不是源于機制模型，再完美的統(tǒng)計模型也只是對野生撫育、栽培產(chǎn)業(yè)化的一個方向性的指引，它的功能在于明確區(qū)劃范圍內(nèi)開展工作的可靠性比較高，而不是表明實踐中必然如此，這也是中藥材區(qū)劃與農(nóng)作物區(qū)劃不一樣的地方。在實際應(yīng)用中，對于中藥材區(qū)劃結(jié)果，首先是對區(qū)劃范圍不同等級適宜區(qū)域進行驗證和試驗示范，然后才能大規(guī)模推廣應(yīng)用，以避免對中藥材次生產(chǎn)物代謝機制、中藥材藥理藥效機制、道地性形成機制研究不足，僅僅依靠統(tǒng)計模型建立的區(qū)劃可能帶來的產(chǎn)業(yè)化栽培中的中藥材質(zhì)量和經(jīng)濟風險。

將藥材生長適宜性因子、藥材道地性因子、藥材質(zhì)量控制特征因子、土地利用方式因子等諸多因素結(jié)合，從地統(tǒng)計學(xué)角度可明確界定適宜的野生資源分布區(qū)域、野生種群保護區(qū)域、藥材生產(chǎn)質(zhì)量保證區(qū)域、野生撫育區(qū)域等，在此基礎(chǔ)上形成一個結(jié)合實際和功能需求的完整區(qū)劃，傳統(tǒng)區(qū)劃方法難以將這些具體的空間實際和現(xiàn)實需求相互結(jié)合并定量實現(xiàn)。這種綜合考慮多因素的功能性生產(chǎn)區(qū)劃理論和方法，可為中藥材戰(zhàn)略資源保護和可持續(xù)利用的理論構(gòu)建提供新思路，也對中藥材的生產(chǎn)實踐具有重要的指導(dǎo)意義。

通過以羌活為例建立的功能型區(qū)劃分析方法，有效的解決因藥材植物積累次生代謝產(chǎn)物所需的適宜生境（藥材品質(zhì)適宜性）與其生長發(fā)育的適宜生境（植物生長適宜性）生長適宜性與品質(zhì)適宜性空間分離問題，而導(dǎo)致的生產(chǎn)區(qū)劃難以兼顧藥材產(chǎn)量及品質(zhì)的問題。通過對羌活的功能型生產(chǎn)區(qū)劃研究表明，羌活的生長適宜性和品質(zhì)適宜性整體上具有一點的吻合性，但是就其高品質(zhì)區(qū)段和生長高密度區(qū)段來說，在空間上存在較明顯偏離。綜合來看，羌活的藥材品質(zhì)及生長適宜性最優(yōu)區(qū)段主要分布在道地產(chǎn)區(qū)的四川省阿壩州及甘孜州境內(nèi)，與文獻描述的傳統(tǒng)產(chǎn)區(qū)較為一致，因此四川省較適宜開展羌活開展規(guī)?；斯しN植和野生撫育的區(qū)段。

[注]本文地圖插圖審圖號：GS（2017）2128號。

[參考文獻]

[1]黃璐琦 中藥區(qū)劃專題編者按[J] 中國中藥雜志，2016，41（17）： 3113

[2]銀玲， 彭月， 劉榮， 等 產(chǎn)地生態(tài)環(huán)境要素與中藥品質(zhì)相關(guān)性研究[J] 中藥與臨床， 2012， 3（6）： 9

[3]黃璐琦， 郭蘭萍 環(huán)境脅迫下次生代謝產(chǎn)物的積累及道地藥材的形成[J] 中國中藥雜志， 2007， 33（4）： 277

[4]孫洪兵， 孫輝， 蔣舜媛， 等 基于3G技術(shù)的羌活區(qū)劃研究I 基于MaxEnt和ArcGIS的羌活生長適宜性分析及評價[J]. 中國中藥雜志， 2015， 40（5）： 853

[5]孫洪兵，蔣舜媛，孫輝， 等基于3S技術(shù)的羌活區(qū)劃研究Ⅱ.基于協(xié)同克里金插值的羌活品質(zhì)適宜性分析及評價價[J] 中國中藥雜志， 2017， 42（13）：2633

[6]蔣舜媛， 孫洪兵， 孫輝， 等 羌活藥材商品規(guī)格與質(zhì)量等級的對應(yīng)分析[J] 中國中藥雜志， 2016， 41（5）： 793

[7]黃波， 徐冠華， 閻守邕 GIS中空間模糊疊加模型的設(shè)計[J]. 測繪學(xué)報，1996， 25 （1）： 53

[8]單鋒， 袁媛， 郝近大， 等 獨活、羌活的本草源流考[J] 中國中藥雜志， 2014， 39（17）： 3399

[9]蔣舜媛， 孫洪兵， 秦紀洪， 等 基于生長適宜性和品質(zhì)適宜性的川貝母功能型生產(chǎn)區(qū)劃研究[J] 中國中藥雜志， 2016， 41（17）： 3194

[責任編輯呂冬梅]