基于DEM的黃土高原地貌分形美學(xué)圖案研究

基于DEM的黃土高原地貌分形美學(xué)圖案研究【摘要】為進(jìn)一步探索黃土高原千溝萬壑、支離破碎的自然景觀特征中蘊(yùn)含的分形美學(xué)價(jià)值，本文基于黃土高原6個(gè)不同地貌類型區(qū)的數(shù)字高程模型，在分形圖案設(shè)計(jì)軟件UltraFracta1中，利用多種標(biāo)準(zhǔn)化分形變換和Newton分形公式，設(shè)計(jì)并對(duì)比使用原始DEM數(shù)據(jù)提取出的地貌因子數(shù)據(jù)與不同分形變換集合得到的分形美學(xué)圖案。結(jié)果表明，對(duì)蘊(yùn)含地貌特征的數(shù)據(jù)進(jìn)行分形設(shè)計(jì)，可以將自然景觀特征與公式化的數(shù)學(xué)思維相結(jié)合，得到具有多樣性、美觀性、裝飾性的圖案。研究結(jié)果為自然景觀的美學(xué)研究提供了一個(gè)新的方向，也可為服飾設(shè)計(jì)、家居用品、文創(chuàng)產(chǎn)品等提供靈感與創(chuàng)意?！娟P(guān)鍵詞】黃土高原；分形；美學(xué)圖案設(shè)計(jì)；DEMAbstract：InordertoexplorethefractalaestheticvaluecontainedintheLoessPlateau'sthousandsofgullyandfragmentednaturallandscapefeatures,basedontheDigitalElevationModelof6differentgeomorphologicalareasoftheLoessPlateau，bythefractalpatterndesignsoftwareUltraFractal,avarietyofstandardizedfractaltransformationsandNewtonfractalformulasareusedtodesignandcomparethefractalaestheticpatternsobtainedbyusingtheoriginalDEMdata,theextractedgeomorphicfactordataanddifferentfractaltransformationsets.Theresultsshowthatthefractaldesignofdatacontaininggeomorphicfeaturescancombinenaturallandscapefeatureswithformulaicmathematicalthinkingtoobtainpatternsthatarediverse,beautifulanddecorative.Theresearchresultsprovideanewdirectionfortheaestheticresearchofnaturallandscapes,aswellasinspirationandcreativityforappareldesign,householdgoods,creativeproducts,etc.Keywords：loessplateau;fractal;aestheticpatterndesign;DEM0引言黃土高原是世界上最大的黃土沉積區(qū)，包含日月山以東、太行山以西，長(zhǎng)城以南、秦嶺以北的廣闊地區(qū)，平均海拔為1500?2000m，面積約為62.14萬km2'1］。這一極具特色的地形地貌區(qū)，經(jīng)流水長(zhǎng)期強(qiáng)烈侵蝕，逐漸形成千溝萬壑、支離破碎的特殊自然景觀，豐富多樣的自然景觀不僅具有極高的地學(xué)研究意義，也蘊(yùn)含著相當(dāng)豐富的美學(xué)價(jià)值。分形常常被定義為一個(gè)粗糙或零碎的幾何形狀，可以分成數(shù)個(gè)部分，且每一部分都（至少近似地）是整體縮小后的形狀，具有自相似性、自仿射性和標(biāo)度不變性。分形幾何作為一門新興的橫斷學(xué)科，自建立以來，很快受到各個(gè)學(xué)科領(lǐng)域的廣泛關(guān)注。孫靜等利用陰香木超微觀構(gòu)造圖案設(shè)計(jì)分形美學(xué)圖案，將豐富的構(gòu)圖關(guān)系融入陰香木超微觀美學(xué)圖案中［2］，韋曉丹等將香樟樹不同微觀構(gòu)造得到的分形圖案設(shè)計(jì)為多種風(fēng)格和用途的裝飾圖案［3］，王旭婷等利用二值化及邊界提取榆木微觀構(gòu)造圖像后，獲得榆木分形圖像的紋理特征［4］;El-DARWISHII介紹了分形在建筑中的演變以及當(dāng)今建筑中分形設(shè)計(jì)面臨的挑戰(zhàn)［5］;蔡燕燕等通過改變參數(shù)探索Mandelbrot集的絲綢圖案設(shè)計(jì)［6］。分形在各個(gè)領(lǐng)域賦予了圖形創(chuàng)作新的形式，但在地形地貌方面尚少被人們發(fā)掘利用。為此，將分形與黃土高原自然景觀所表現(xiàn)出的自相似性結(jié)

合，進(jìn)而探索開發(fā)黃土高原的自然美學(xué)價(jià)值，將地學(xué)美融入生產(chǎn)生活中，對(duì)于豐富地貌美學(xué)研究、拓展美學(xué)設(shè)計(jì)具有十分重要的意義。1材料與方法1.1實(shí)驗(yàn)材料根據(jù)黃土高原地區(qū)綜合治理規(guī)劃大綱，黃土高原可劃分為6個(gè)大的地貌類型區(qū)：黃土高塬溝壑區(qū)、黃土丘陵溝壑區(qū)、土石山區(qū)、河谷平原區(qū)、沙地和沙漠區(qū)、農(nóng)灌區(qū)［7］o本文分別從這6個(gè)地貌類型區(qū)中選取甘肅省永靖縣、山西省臨縣、河南省嵩縣、陜西省藍(lán)田縣、內(nèi)蒙古自治區(qū)杭錦旗、寧夏回族自治區(qū)鹽池縣為研究對(duì)數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺(tái)SRTMDEM90m原始高程數(shù)據(jù)(圖2)。象(圖1),并收集了6個(gè)縣的數(shù)字高程模型(DEM)，圖例 ,01rv~innmmLiJr^~i審圖號(hào)：GS(2020)4631號(hào)圖1象(圖1),并收集了6個(gè)縣的數(shù)字高程模型(DEM)，圖例 ,01rv~innmmLiJr^~i審圖號(hào)：GS(2020)4631號(hào)圖1黃土高原綜合治理區(qū)及研究區(qū)域分布情況Fig.1DistributionofcomprehensivemanagementareasandresearchareasontheLoessPlateau圖26個(gè)研究區(qū)的DEM數(shù)據(jù)及其窗口裁剪情況Fig.2DEMdataof6studyareasandtheirwindowclipping山地，該區(qū)域中部高，東西兩側(cè)較低，高程分布具有一定對(duì)稱性和規(guī)律性，窗口大小為47.541X48.664km2；圖2d為沙地和沙漠區(qū)中的內(nèi)蒙古杭錦旗，風(fēng)沙地貌顯著，圖像紋理較少，高程變化規(guī)律性較強(qiáng)且地面較為平緩，窗口大小為84.825X84.601km2；圖2e為黃土丘陵溝壑區(qū)中的山西臨縣，該區(qū)域內(nèi)地貌類型主要為黃土梁和黃土溝，地表支離破碎，溝谷縱橫，細(xì)節(jié)豐富，分形美感較高，窗口大小為38.148X34.137km2；圖2f為土石山區(qū)中的河南嵩縣，地形破碎復(fù)雜［9］，圖像紋理無顯著規(guī)律，區(qū)域范圍內(nèi)相對(duì)高差較大，窗口大小為34.621X36.113km2。1.2實(shí)驗(yàn)方法將實(shí)驗(yàn)區(qū)jpg格式的原始素材導(dǎo)入U(xiǎn)ltraFractal軟件，選擇Kaleidoscope、Twist、Ripples等標(biāo)準(zhǔn)化分形變換并修改所對(duì)應(yīng)的參數(shù)，然后加載Newton分形公式進(jìn)行非線性迭代，通過軟件預(yù)覽窗口觀察分形圖案并導(dǎo)出最終的結(jié)果圖案。標(biāo)準(zhǔn)化分形變換是通過分形圖案整體的變形和扭曲來獲得不同的分形形狀。本文使用的Kaleidoscope變換將原圖案的徑向切片進(jìn)行旋轉(zhuǎn)對(duì)稱，Twist變換為原圖案增加一個(gè)小范圍的螺旋形變形，Ripples變換為原圖案增加波紋效果。還可以通過調(diào)整變換中心坐標(biāo)的實(shí)部(Re)和虛部(Im)獲得不同形狀的分形圖案。標(biāo)準(zhǔn)化分形變換可以疊加使用，多重變換疊加的效果是在前一個(gè)變換的基礎(chǔ)上進(jìn)行的新變換，而不是對(duì)原圖案分別進(jìn)行單獨(dú)變換后再疊加。分形公式用于創(chuàng)建分形的基本形狀和樣式。本文采用Newton分形公式，該公式是通過牛頓法求解多項(xiàng)式方程而生成的。設(shè)是關(guān)于復(fù)數(shù)z的非線性函數(shù)，且f'(z) 0。牛頓迭代法是將f(z)=0逐次線性化的一種近似方法，迭代公式為：Zn+1=Zn-常,（n=0,1,2，…） （1）對(duì)于給定的初值點(diǎn)％,產(chǎn)生的序列｛%｝將出現(xiàn)收斂和發(fā)散兩種情況，可以通過改變方程f（z）=0來繪制分形圖形M。2結(jié)果與分析2.1基于DEM數(shù)據(jù)的分形圖案設(shè)計(jì)在Newton公式集合下，使用Kaleidoscope變換對(duì)6個(gè)來自不同地貌區(qū)的原始DEM圖像進(jìn)行兩次變換設(shè)計(jì)，得到黃土高塬溝壑區(qū)（圖3A）、黃土丘陵溝壑區(qū)（圖3B）、土石山區(qū)（圖3C）、河谷平原區(qū)（圖3D）、沙地和沙漠區(qū)（圖3E）、農(nóng)灌區(qū)（圖3F）6個(gè)不同的分形圖案，設(shè)計(jì)中使用的參數(shù)見表1。圖36種不同地貌設(shè)計(jì)出的分形圖案Fig.3Fractalpatternsdesignedfrom6differentlandforms表16種不同地貌分形圖案的設(shè)計(jì)參數(shù)Tab.1Designparametersof6differenttopographicfractalpatterns編號(hào)地貌區(qū)所在行政區(qū)-參數(shù)Center(Re)Center(Im)A黃土高塬甘肅省永靖縣-1.484086,-0.998371,溝壑區(qū)0.3933410.136198B黃土丘陵山西省臨縣-2.988848,-0.807895,溝壑區(qū)0.19334-0.120945C土石山區(qū)河南省嵩縣-0.703133,-1.826943,-0.520944-0.035231D河谷平原區(qū)陜西省藍(lán)田縣-0.440884,-0.202789,-0.1517650.162521E沙地和沙漠區(qū)內(nèi)蒙古自治區(qū)-0.24599,-1.884086,杭錦旗-0.482849-0.111422F農(nóng)灌區(qū)寧夏回族自治-0.722181,-1.884086,區(qū)鹽池縣-0.406659-0.111422分形圖案的自相似性與中國傳統(tǒng)對(duì)稱美學(xué)中以點(diǎn)面為核心、色彩和諧統(tǒng)一、兩側(cè)圖案對(duì)稱一致等特征有著極高的相似性，能夠被廣泛地應(yīng)用到各類作品當(dāng)中。通過在不同地貌區(qū)疊加相同分形函數(shù)集合來實(shí)現(xiàn)不同樣式的對(duì)稱模擬，不同地貌類型區(qū)的DEM圖像的紋理、灰度均不同，因此即使是同一種分形函數(shù)集合操作也使得設(shè)計(jì)出的圖案千差萬別，并且每種圖案都有其獨(dú)特的美感。由圖3可以看出，以上6種不同地貌類型區(qū)的DEM數(shù)據(jù)分形圖案在構(gòu)成、紋理、灰度等方面均呈現(xiàn)出不同的風(fēng)格，且具有極強(qiáng)的對(duì)稱性。在DEM數(shù)據(jù)中，高程以灰度的形式反映，灰度越大，高程越大，因此高程偏小的地區(qū)在圖像中明度較低，如甘肅永靖縣的分形圖案（圖3A），亮度整體偏暗但紋理清晰，傳遞出有條不紊、簡(jiǎn)潔干練之感；高程分布較為均勻地區(qū)生成的圖案具有較低的對(duì)比度，圖像的灰度變化較小，如內(nèi)蒙古自治區(qū)杭錦旗生成的圖案（圖3E），其整體均勻和諧，宛如荷葉的葉脈，極具朦朧端莊的傳統(tǒng)美；高差較大的地區(qū)，即坡度變化大的地區(qū)具有更加豐富的紋理特征，以位于土石山區(qū)的河南嵩縣為例，其坡度從0°?46.9°均有分布，并且集中分布于10°—25°之間，因此由它生成的圖案（圖3C）紋理豐富且對(duì)比明顯，具有較強(qiáng)的立體感，加之其明暗交替出現(xiàn)形成了極富動(dòng)感的形態(tài)。在山脊山谷發(fā)育明顯的地區(qū)，其圖案輪廓細(xì)致分明，位于黃土丘陵溝壑區(qū)的山西臨縣其境內(nèi)山脈連綿延長(zhǎng)，因此生成的圖案（圖3B）輪廓明晰，層次分明，具有更高的美感度。2.2基于地貌因子數(shù)據(jù)的分形圖案設(shè)計(jì)為進(jìn)一步說明不同地貌類型區(qū)地形特征的美學(xué)價(jià)值，以藍(lán)田縣為例對(duì)其等值線、水系、坡度、山脊線、山谷線等不同地貌因子在Newton公式集合下，使用Kaleidoscope變換進(jìn)行了分形圖案設(shè)計(jì)。所得分形結(jié)果如圖4所示。a數(shù)字高程模型 b100m間距等值線 c200m間距等值線adigitalelevationmodel b100mspacingcontourc200mspacingcontourd坡度 e山谷線 f0m閾值水系dslope evalleyline f0mthresholdhydrographicneth山脊線hridgelineh山脊線hridgelineg100m閾值水系

g100mthreshold

hydrographicnet圖4不同地貌因子設(shè)計(jì)出的分形圖案Fig.4DesignedfractalpatternsunderdifferentterrainfactorsDEM制圖與地貌因子制圖之間的差異主要體現(xiàn)在圖像灰度變化上。由于地面的高低起伏是連續(xù)的，因此圖像中用于表示高程的灰度值的變化也是連續(xù)的。而等值線、山脊線、山谷線等地貌因子經(jīng)過軟件計(jì)算，它們之間的變化是間斷的，那么相應(yīng)圖像灰度值的變化也不連續(xù)，在圖像上更多表現(xiàn)的是線條樣式圖案。另外，雖然地貌因子數(shù)據(jù)設(shè)計(jì)時(shí)使用的分形公式和變換都相同，最終得到的圖案也有較相似的樣式，但是在灰度、紋理、輪廓等方面仍有較大差別，且各具特色。其中，將不同間距等值線設(shè)計(jì)出的圖案對(duì)比可以發(fā)現(xiàn)：等值線間距越小，圖案的線條越密、灰度值越小、局部紋理越豐富多樣。不同閾值水系數(shù)據(jù)設(shè)計(jì)出的圖案也有相似的規(guī)律。2.3不同變換集合的疊加圖案設(shè)計(jì)黃土高塬溝壑區(qū)作為黃土高原生態(tài)分區(qū)中分布面積最廣的區(qū)域，塬面廣闊平坦、溝壑縱橫深切，極具黃土高原的自然景觀特色，且地形支離破碎具有良好的天然分形特征。本文以黃土高塬溝壑區(qū)甘肅永靖縣DEM數(shù)據(jù)為例，在Newton公式集合下，采用不同分形變換進(jìn)行疊加來構(gòu)造分形圖案，說明通過不同分形公式疊加變換得到的圖案在結(jié)構(gòu)和內(nèi)容上的差別（圖5、表2）。圖案A1有規(guī)律的重復(fù)和排列，細(xì)膩的紋路和精致的花紋，凸顯出古典氣息；圖案A2由四周向中央層層遞進(jìn)，原素材中的地形起伏轉(zhuǎn)變?yōu)橛擅靼到诲e(cuò)的線條填充的花紋，宛如鏡像，遠(yuǎn)看有立體高低起伏之感；圖案A3如同樹木的年輪層層展開，有著嚴(yán)格的中心對(duì)稱，飽含莊重嚴(yán)肅感；圖案A4輪廓簡(jiǎn)潔大方，宛如天空祥云，中心圖案波紋漣漪，宛如湖水碧波，極具中國古代花紋蜿蜒曲折之美；圖案A5猶如花瓣的有序排列，其中的明暗紋理皆是錯(cuò)落有序，雜而不亂，蘊(yùn)含著平衡、穩(wěn)定之美；圖案A6是旋轉(zhuǎn)對(duì)稱圖形，其簡(jiǎn)練單純的形狀充滿了運(yùn)動(dòng)的姿態(tài)，具有很好韻律感以及齊

圖5不同公式變換下黃土高塬溝壑區(qū)甘肅省永靖縣DEM數(shù)據(jù)的分形圖案Fig.5FractalpatternsofDEMdatainYongjingCounty,GansuProvinceunderthetransformationofdifferentformulas表2黃土高塬溝壑區(qū)甘肅永靖縣DEM數(shù)據(jù)的分形圖案設(shè)計(jì)相關(guān)參數(shù)

Tab.2RelatedparametersoffractalpatterndesignforDEMdatain

YongjingCounty,GansuProvince編號(hào)變換類型變換參數(shù)Center(Re)Center(Im)Twist-1.090476-0.166667A1Kaleidoscope-1.036466-0.435229Kaleidoscope-2.1221810.06953Twist-0.509524-0.947619A2Kaleidoscope-2.265037-0.320943Ripples-2.30307981750.11644357Kaleidoscope-2.5602760.421911Inverse00A3Kaleidoscope-1.998370.18382Kaleidoscope-1.541829-0.269481Inverse0.0380960.257142A4Kale