GB∕T 43605-2023 大氣腐蝕圖繪制方法（正式版）

大氣腐蝕圖繪制方法2023-12-28發(fā)布國家市場監(jiān)督管理總局國家標(biāo)準(zhǔn)化管理委員會I前言 l2規(guī)范性引用文件 3術(shù)語和定義 4一般規(guī)定 5繪制流程 26繪圖數(shù)據(jù)準(zhǔn)備 37地理空間插值模型構(gòu)建 48數(shù)據(jù)源空白區(qū)域賦值 59圖面繪制 510圖面輸出 711編制繪制說明 7附錄A(資料性)地理空間插值模型 8附錄B(資料性)留一交叉驗(yàn)證方法 附錄C(資料性)大氣腐蝕圖色標(biāo) 附錄D(資料性)大氣腐蝕圖圖面示例 附錄E(資料性)《山東省大氣腐蝕圖繪制說明》示例 參考文獻(xiàn) ⅢGB/T43605—2023本文件按照GB/T1.1—2020《標(biāo)準(zhǔn)化工作導(dǎo)則第1部分：標(biāo)準(zhǔn)化文件的結(jié)構(gòu)和起草規(guī)則》的規(guī)定起草。請注意本文件的某些內(nèi)容可能涉及專利。本文件的發(fā)布機(jī)構(gòu)不承擔(dān)識別專利的責(zé)任。本文件由中國機(jī)械工業(yè)聯(lián)合會提出。本文件由全國金屬與非金屬覆蓋層標(biāo)準(zhǔn)化技術(shù)委員會(SAC/TC57)歸口。本文件起草單位：國網(wǎng)山東省電力公司電力科學(xué)研究院、中國機(jī)械總院集團(tuán)武漢材料保護(hù)研究所有限公司、國網(wǎng)湖北省電力有限公司電力科學(xué)研究院、國網(wǎng)浙江省電力有限公司電力科學(xué)研究院、冶金工業(yè)信息標(biāo)準(zhǔn)研究院、朗松珂利(上海)儀器儀表有限公司。1大氣腐蝕圖繪制方法本文件描述了基于大氣環(huán)境腐蝕試驗(yàn)及數(shù)據(jù)繪制區(qū)域性大氣腐蝕圖的方法，包括一般規(guī)定、繪制流本文件適用于以金屬材料腐蝕速率評價(jià)大氣腐蝕性進(jìn)而繪制的大氣腐蝕圖。注：本文件所述金屬材料是指金屬、合金和金屬類覆蓋層。2規(guī)范性引用文件下列文件中的內(nèi)容通過文中的規(guī)范性引用而構(gòu)成本文件必不可少的條款。其中，注日期的引用文件，僅該日期對應(yīng)的版本適用于本文件；不注日期的引用文件，其最新版本(包括所有的修改單)適用于本文件。GB/T13989國家基本比例尺地形圖分幅和編號GB/T14165金屬和合金大氣腐蝕試驗(yàn)現(xiàn)場試驗(yàn)的一般要求GB/T14511地圖印刷規(guī)范GB/T19292.1—2018金屬和合金的腐蝕大氣腐蝕性第1部分：分類、測定和評估GB/T19292.3—2018金屬和合金的腐蝕大氣腐蝕性第3部分：影響大氣腐蝕性環(huán)境參數(shù)的測量GB/T19292.4—2018金屬和合金的腐蝕大氣腐蝕性第4部分：用于評估腐蝕性的標(biāo)準(zhǔn)試樣的腐蝕速率的測定GB22021國家大地測量基本技術(shù)規(guī)定GB35650國家基本比例尺地圖測繪基本技術(shù)規(guī)定3術(shù)語和定義GB/T19292.1—2018界定的以及下列術(shù)語和定義適用于本文件。大氣腐蝕圖atmosphericcorrosionmap描述某區(qū)域在某時(shí)間段內(nèi)大氣環(huán)境對金屬材料腐蝕性強(qiáng)弱的地理分布圖。注：大氣腐蝕圖包括大氣腐蝕性等級圖和金屬材料大氣腐蝕速率圖。描述大氣環(huán)境腐蝕性等級的地理分布圖。金屬材料大氣腐蝕速率圖atmosphericcorrosionratemapofmetallicmaterials描述某種金屬材料在大氣環(huán)境中腐蝕速率大小的地理分布圖。繪制大氣腐蝕圖所使用的地理地圖。2在大氣腐蝕圖中，相鄰兩個(gè)區(qū)域的大氣腐蝕性等級或金屬材料的腐蝕速率分段出現(xiàn)2個(gè)及以上的數(shù)據(jù)源點(diǎn)datasourcesite用于繪制大氣腐蝕圖的實(shí)測繪圖數(shù)據(jù)的地理坐標(biāo)。在大氣腐蝕圖的繪制區(qū)域中，數(shù)據(jù)源點(diǎn)以外的其他區(qū)域。地理空間插值模型geospatialinterpolationmodel在大氣腐蝕圖的繪制區(qū)域范圍內(nèi)，采用數(shù)據(jù)源點(diǎn)的實(shí)測數(shù)據(jù)估算數(shù)據(jù)源空白區(qū)域數(shù)據(jù)的計(jì)算方法。用相等數(shù)值點(diǎn)的連線表示大氣腐蝕數(shù)據(jù)空間連續(xù)分布且逐漸變化的地理分布圖。色斑圖spottedmap選取不同色域填充表示大氣腐蝕數(shù)據(jù)空間連續(xù)分布且逐漸變化的地理分布圖。4一般規(guī)定繪制大氣腐蝕圖應(yīng)滿足以下要求：a)繪制大氣腐蝕圖應(yīng)符合GB22021、GB35650等國家測繪標(biāo)準(zhǔn)，中國全圖應(yīng)詳列南海諸島和釣魚島；b)大氣腐蝕圖分為大氣腐蝕性等級圖和金屬材料大氣腐蝕速率圖；c)大氣腐蝕性等級圖中腐蝕性分級應(yīng)符合GB/T19292.1—2018的規(guī)定；d)金屬材料大氣腐蝕速率圖應(yīng)使用金屬材料的年腐蝕速率繪制；e)大氣腐蝕圖應(yīng)包含圖面和繪制說明；f)大氣腐蝕圖應(yīng)覆蓋繪制的全部區(qū)域；g)毗鄰地區(qū)交界區(qū)域的大氣腐蝕性數(shù)據(jù)不應(yīng)出現(xiàn)躍變5繪制流程大氣腐蝕圖的繪制主要包括繪圖數(shù)據(jù)準(zhǔn)備、地理空間插值模型構(gòu)建、數(shù)據(jù)源空白區(qū)域賦值、圖面繪制、圖面輸出、編制繪制說明等步驟，如圖1所示。3數(shù)據(jù)獲取數(shù)據(jù)獲取賦值數(shù)據(jù)的誤差評價(jià)繪圖軟件大氣腐蝕性等級圖圖面繪制圖面輸山據(jù)備數(shù)準(zhǔn)繪圖值圖1大氣腐蝕圖繪制流程6繪圖數(shù)據(jù)準(zhǔn)備6.1數(shù)據(jù)獲取6.1.1用于繪制大氣腐蝕圖的數(shù)據(jù)應(yīng)至少包括數(shù)據(jù)源點(diǎn)的地理坐標(biāo)、金屬材料大氣腐蝕的年腐蝕速率和獲取數(shù)據(jù)的時(shí)間。數(shù)據(jù)的格式見表1。6.1.2在大氣腐蝕性等級圖繪制中，用于評定大氣腐蝕性強(qiáng)弱的數(shù)據(jù)應(yīng)按照GB/T14165和GB/T19292.4—2018中規(guī)定的方法獲取，或在繪制區(qū)域內(nèi)檢測規(guī)定金屬材料在戶外暴露第1年的腐蝕速率。也可根據(jù)GB/T19292.1—2018中的劑量-響應(yīng)函數(shù)，通過計(jì)算獲得該地區(qū)規(guī)定金屬材料的年腐蝕速率；或采用經(jīng)驗(yàn)證的其他劑量-響應(yīng)函數(shù)計(jì)算獲取。6.1.3在金屬材料大氣腐蝕速率圖繪制中，繪圖數(shù)據(jù)應(yīng)根據(jù)GB/T14165和GB/T19292.4—2018中規(guī)定的方法獲取金屬材料在某一時(shí)間段的腐蝕速率，或采用經(jīng)驗(yàn)證的劑量-響應(yīng)函數(shù)計(jì)算獲取，其計(jì)量單位均應(yīng)為年腐蝕速率。4表1繪圖數(shù)據(jù)格式數(shù)據(jù)名稱單位/符號備注獲取數(shù)據(jù)的時(shí)間年/月/日—年/月/日自然暴露試驗(yàn)的時(shí)間周期；劑量-響應(yīng)函數(shù)中環(huán)境因子獲取周期；或者由其他經(jīng)驗(yàn)證的方法獲取數(shù)據(jù)的周期地理坐標(biāo)××°××'××”(S/N),數(shù)據(jù)源點(diǎn)的經(jīng)緯度坐標(biāo)金屬材料的年腐蝕速率g/(m2·a)或μm/a采用自然暴露試驗(yàn)、劑量-響應(yīng)函數(shù)或者其他經(jīng)驗(yàn)證方法計(jì)算獲取的年化腐蝕速率數(shù)據(jù)6.1.4采用的劑量-響應(yīng)函數(shù)和大氣腐蝕監(jiān)/檢測儀器等應(yīng)經(jīng)過該地區(qū)至少1年的自然暴露試驗(yàn)驗(yàn)證，影響大氣腐蝕性的環(huán)境參數(shù)應(yīng)按照GB/T19292.3—2018等方法獲取。6.1.5繪圖數(shù)據(jù)亦可采用經(jīng)驗(yàn)證的大氣腐蝕監(jiān)/檢測儀器等試驗(yàn)方法獲取。當(dāng)檢測周期不足1年時(shí)，可采用插值計(jì)算方法補(bǔ)充為年腐蝕速率，但原始數(shù)據(jù)的檢測周期應(yīng)不少于3個(gè)月，且每日不少于1組數(shù)據(jù)。6.1.6采用大氣腐蝕監(jiān)/檢測儀器等試驗(yàn)方法獲取繪圖數(shù)據(jù)時(shí)，繪制的大氣腐蝕圖應(yīng)隨數(shù)據(jù)變化而動(dòng)態(tài)更新。6.1.7不同方法獲取的同一時(shí)間段的繪圖數(shù)據(jù)，可在同一張大氣腐蝕圖的繪制中同時(shí)使用。6.2數(shù)據(jù)驗(yàn)證與確認(rèn)6.2.1獲取繪圖數(shù)據(jù)后，應(yīng)對數(shù)據(jù)的準(zhǔn)確性進(jìn)行驗(yàn)證與確認(rèn)，包括獲取方法、測試過程和數(shù)據(jù)誤差。6.2.2若采用GB/T14165和GB/T19292.4—2018中規(guī)定之外的方法獲取繪圖數(shù)據(jù)，應(yīng)在繪制說明中詳述繪圖數(shù)據(jù)的獲取方法和數(shù)據(jù)誤差驗(yàn)證結(jié)果。6.3其他要求6.3.1繪制大氣腐蝕性等級圖時(shí)，應(yīng)優(yōu)先采用GB/T19292.1—2018中規(guī)定的標(biāo)準(zhǔn)金屬鋅、碳鋼、銅和鋁經(jīng)自然暴露試驗(yàn)獲取的年腐蝕速率數(shù)據(jù)。6.3.2地理相鄰的繪圖數(shù)據(jù)所對應(yīng)的腐蝕性等級或腐蝕速率分段之間不應(yīng)出現(xiàn)躍變。否則，應(yīng)在該相鄰數(shù)據(jù)的地理坐標(biāo)之間的位置增加數(shù)據(jù)源點(diǎn)。6.3.3繪圖數(shù)據(jù)宜覆蓋繪制區(qū)域內(nèi)的全部氣候帶以及城市大氣、鄉(xiāng)村大氣、工業(yè)大氣、海洋大氣等不同自然環(huán)境。數(shù)據(jù)源點(diǎn)宜均勻分布。6.3.4應(yīng)在圖面上的合適位置注明繪圖數(shù)據(jù)來源。7地理空間插值模型構(gòu)建7.1模型準(zhǔn)備7.1.1應(yīng)根據(jù)已有繪圖數(shù)據(jù)選擇或建立地理空間插值模型，模型應(yīng)能適用于繪制區(qū)域內(nèi)全部數(shù)據(jù)源空白區(qū)域的賦值。7.1.2用于數(shù)據(jù)源空白區(qū)域賦值的地理空間插值模型包括反距離加權(quán)法、克里金法、趨勢面法和泰森多邊形法等，參見附錄A。在可提高插值精度的前提下，繪圖者可自行建立插值模型。7.2模型的評價(jià)與優(yōu)化7.2.1應(yīng)對地理空間插值模型進(jìn)行插值誤差的評價(jià)。57.2.2宜采用留一交叉驗(yàn)證、自然暴露試驗(yàn)驗(yàn)證等方法定量評價(jià)模型的插值誤差，并根據(jù)評價(jià)結(jié)果進(jìn)行模型優(yōu)化。留一交叉驗(yàn)證方法參見附錄B,自然暴露試驗(yàn)驗(yàn)證方法應(yīng)按照GB/T14165執(zhí)行。7.2.3模型的平均插值誤差應(yīng)不大于繪制的圖面中金屬材料腐蝕速率的分段跨度值。否則，應(yīng)重新優(yōu)化或選擇。7.2.4宜通過增加數(shù)據(jù)源點(diǎn)密度、優(yōu)化地理空間插值模型等方法提升數(shù)據(jù)源空白區(qū)域的賦值精度。8數(shù)據(jù)源空白區(qū)域賦值8.1賦值數(shù)據(jù)的計(jì)算8.1.1應(yīng)采用經(jīng)優(yōu)化的地理空間插值模型進(jìn)行數(shù)據(jù)源空白區(qū)域賦值，并保證覆蓋全部繪制區(qū)域。8.1.2確定空間插值的網(wǎng)格形式和網(wǎng)格密度時(shí)，插值的網(wǎng)格形式宜優(yōu)先選用離散型三角網(wǎng)格，在賦值曲率較大的區(qū)域宜做加密插值，在曲率較小的區(qū)域可不進(jìn)行加密處理。8.1.3賦值的網(wǎng)格密度宜在1km×1km至100km×100km范圍內(nèi)選取。8.1.4應(yīng)根據(jù)賦值計(jì)算結(jié)果和8.2中的誤差評價(jià)結(jié)果對插值模型、網(wǎng)格形式和網(wǎng)格密度進(jìn)行優(yōu)化，以保證網(wǎng)格相鄰節(jié)點(diǎn)的賦值結(jié)果之間不出現(xiàn)躍變。8.2賦值數(shù)據(jù)的誤差評價(jià)8.2.1宜采用留一交叉驗(yàn)證的方法，對繪圖區(qū)域內(nèi)的全部數(shù)據(jù)源點(diǎn)的賦值數(shù)據(jù)進(jìn)行驗(yàn)證，計(jì)算全部驗(yàn)證結(jié)果并進(jìn)行誤差的統(tǒng)計(jì)計(jì)算。8.2.2在數(shù)據(jù)源空白區(qū)域宜采用自然暴露試驗(yàn)的方法，對賦值數(shù)據(jù)進(jìn)行誤差評價(jià)。8.2.3宜采用平均誤差、平均絕對誤差、均方根誤差描述數(shù)據(jù)源空白區(qū)域賦值數(shù)據(jù)誤差的大小，并應(yīng)滿足7.2.3的要求。否則，應(yīng)依照8.1中的要求進(jìn)行優(yōu)化處理。8.2.4應(yīng)在大氣腐蝕圖繪制說明中詳述數(shù)據(jù)源空白區(qū)域賦值的誤差評價(jià)結(jié)果。9圖面繪制9.1繪圖軟件9.1.1繪圖軟件宜選用能夠滿足繪圖數(shù)據(jù)處理、空白區(qū)域賦值、圖面渲染等要求的專業(yè)性計(jì)算機(jī)輔助繪圖軟件。9.1.2繪制者亦可自行開發(fā)繪圖軟件。9.2底圖的選擇確認(rèn)繪圖區(qū)域并選擇底圖。所選底圖應(yīng)符合GB22021、GB35650等國家測繪標(biāo)準(zhǔn)要求。9.3大氣腐蝕性等級圖腐蝕性分級與色標(biāo)的選擇9.3.1繪制大氣腐蝕性等級圖時(shí)，對應(yīng)各個(gè)大氣腐蝕性級別應(yīng)采用不同色標(biāo)在圖中標(biāo)識。9.3.2當(dāng)大氣腐蝕性按照GB/T19292.1—2018分級時(shí)，色標(biāo)宜根據(jù)表2執(zhí)行。9.3.3根據(jù)不同的使用目的，大氣腐蝕圖繪制者在GB/T19292.1—2018劃分的C1～CX基礎(chǔ)上，可進(jìn)一步細(xì)化大氣腐蝕性分級。細(xì)化分級色標(biāo)的選擇可參照附錄C執(zhí)行。6表2大氣腐蝕性分級及色標(biāo)大氣腐蝕性分級色標(biāo)名稱填充顏色(RGB)色標(biāo)藍(lán)綠黃橙紅黑9.4金屬材料大氣腐蝕速率圖的腐蝕速率分段與色標(biāo)的選擇9.4.1繪制金屬材料大氣腐蝕速率圖時(shí)，應(yīng)對所繪金屬材料的大氣腐蝕速率值加以分段，并采用不同的色標(biāo)進(jìn)行標(biāo)識。大氣腐蝕性等級所對應(yīng)的腐蝕速率分段相同。色標(biāo)的選擇可參照表2和附錄C執(zhí)行。9.4.3金屬材料為GB/T19292.1—2018中規(guī)定之外的金屬時(shí)，應(yīng)預(yù)先估計(jì)在繪圖區(qū)域內(nèi)可能出現(xiàn)的最小和最大腐蝕速率值，宜根據(jù)需要合理分段，并根據(jù)最終的空白區(qū)域賦值結(jié)果加以調(diào)整。色標(biāo)的選擇可參照附錄C執(zhí)行。9.5圖面渲染9.5.1繪制的大氣腐蝕圖圖面包括等值線圖和色斑圖兩種形式，繪制者可依據(jù)需要加以選擇。具體圖面形式參見附錄D。9.5.2依據(jù)選擇的圖面形式，按數(shù)據(jù)源點(diǎn)和像元填值，將對應(yīng)的色標(biāo)渲染在底圖上。9.6圖面內(nèi)容9.6.1大氣腐蝕圖圖面中除展示大氣腐蝕性分級或金屬材料大氣腐蝕速率結(jié)果外，還應(yīng)包含標(biāo)題、時(shí)段、圖例、發(fā)布機(jī)構(gòu)、發(fā)布時(shí)間、底圖中的主要地理信息以及其他必要的信息。9.6.2圖面中應(yīng)具有顯示圖中任意位置的經(jīng)緯度坐標(biāo)以及大氣腐蝕性等級或金屬材料年腐蝕速率的功能。9.6.3圖面標(biāo)題應(yīng)簡明扼要，中文標(biāo)注“××區(qū)域大氣腐蝕性等級圖”“××區(qū)域××大氣腐蝕速率圖”,宜置于圖的上方位置。9.6.5可根據(jù)圖面形式，選擇腐蝕分級色塊、等值線、數(shù)據(jù)源點(diǎn)和行政區(qū)劃等元素標(biāo)注圖例。比例尺等信息可位于圖例下方。色斑圖圖例中的色塊可橫向分布，對應(yīng)的文字說明在色塊下方；也可縱向分布，對應(yīng)的文字信息在色塊右側(cè)。9.6.6在圖面的下方應(yīng)標(biāo)注發(fā)布機(jī)構(gòu)全稱和發(fā)布時(shí)間，并使用“yyyy年mm月dd日”格式標(biāo)注發(fā)布時(shí)間。字的大小應(yīng)比標(biāo)題小一號或者小兩號。9.6.7除上述內(nèi)容以外，圖中還可標(biāo)注其他有必要特別說明的信息。710圖面輸出大氣腐蝕圖的輸出形式按數(shù)據(jù)結(jié)構(gòu)分類，分為矢量圖、柵格圖和矢柵混合圖。大氣腐蝕圖的分幅輸出應(yīng)按照GB/T13989的要求執(zhí)行，或根據(jù)實(shí)際需要分幅。大氣腐蝕圖的印刷出版發(fā)行應(yīng)符合GB/T14511的要求。11編制繪制說明11.1大氣腐蝕圖的繪制說明是大氣腐蝕圖的必要組成部分，其文件名稱應(yīng)為《××區(qū)域大氣腐蝕圖繪11.2大氣腐蝕圖繪制說明中應(yīng)詳述以下主要內(nèi)容：a)繪圖數(shù)據(jù)的來源，包括數(shù)據(jù)的獲取方式、數(shù)據(jù)源點(diǎn)的分布、試驗(yàn)材料的種類、數(shù)據(jù)獲取的周期和數(shù)據(jù)準(zhǔn)確性的驗(yàn)證情況；b)采用的地理空間插值模型及其誤差評價(jià)方法和評價(jià)結(jié)果；c)大氣腐蝕性分級依據(jù)與色標(biāo)選擇、金屬材料大氣腐蝕速率的分段依據(jù)與色標(biāo)選擇；d)繪圖過程中采用相關(guān)標(biāo)準(zhǔn)的情況；e)繪圖數(shù)據(jù)與大氣腐蝕圖的更新時(shí)間周期；f)大氣腐蝕圖有效期的建議；g)對大氣腐蝕圖的使用可能產(chǎn)生影響的其他信息；h)數(shù)據(jù)提供者和發(fā)布者以及版權(quán)歸屬的有關(guān)信息。具體編制可參照附錄E。11.3大氣腐蝕圖繪制說明應(yīng)與大氣腐蝕圖同步發(fā)布。8(資料性)地理空間插值模型A.1反距離加權(quán)法反距離加權(quán)法的理論基礎(chǔ)是“地理學(xué)第一定律”,即鄰近的區(qū)域比距離遠(yuǎn)的區(qū)域更相似。它是以插值點(diǎn)與樣本點(diǎn)間的距離為權(quán)重進(jìn)行加權(quán)平均，離插值點(diǎn)越近的樣本點(diǎn)賦予的權(quán)重越大。具體見公式(A.1)。Z?=Z;-Z;a;……………(A.1)式中：Z?——插值點(diǎn)的估計(jì)值；n——估算中用到控制點(diǎn)的數(shù)目；Z;——樣本點(diǎn)i的實(shí)際值；λ;——控制點(diǎn)的權(quán)重。控制點(diǎn)的權(quán)重見公式(A.2)。式中：d;——控制點(diǎn)i與插值點(diǎn)間的距離，當(dāng)插值點(diǎn)與某一樣本點(diǎn)距離為0時(shí)，插值點(diǎn)的值等于該樣本點(diǎn)實(shí)測值，其他情況下，插值點(diǎn)估值為控制點(diǎn)數(shù)值加權(quán)求和；p——加權(quán)指數(shù)。由于距離較近的站點(diǎn)對預(yù)測值的影響權(quán)重更大，且站點(diǎn)之間的距離越大，它們之間的關(guān)聯(lián)性越小。為縮短計(jì)算時(shí)間，通常采用指定搜索鄰域的方法來優(yōu)化模型，即在插值過程中，只考慮距離預(yù)測點(diǎn)一定距離范圍內(nèi)的站點(diǎn)影響，將距離較遠(yuǎn)的對預(yù)測影響微弱的站點(diǎn)排除在外。A.2克里金法克里金法(Kriging)是以變異函數(shù)理論和結(jié)構(gòu)分析為基礎(chǔ)，衡量所選已知點(diǎn)之間空間相關(guān)程度的測度是半變異函數(shù)，在有限區(qū)域內(nèi)對區(qū)域化變量進(jìn)行無偏最優(yōu)估計(jì)的一種方法。衡量各點(diǎn)之間空間相關(guān)程度的測度是半方差，見公式(A.3)。式中：γ(h)——已知點(diǎn)間的半方差；h——已知點(diǎn)之間的距離，常用于作為滯后系數(shù)；n——被h分開的成對樣本點(diǎn)的數(shù)量；z(x;)——區(qū)域化變量g(x)在空間位置x;處的真實(shí)值；z(x;+h)——區(qū)域化變量z(x)在空間位置x;+h處的真實(shí)值。普通克里金插值法估算某點(diǎn)的Z值的通用方程見公式(A.4)和公式(A.5)。Z?=E?=?Z?W?9GB/T43605—2023式中：Z?！烙?jì)值；n用于估算的已知點(diǎn)的數(shù)目；Z,——已知點(diǎn)的值；W——與每個(gè)已知點(diǎn)關(guān)聯(lián)的權(quán)重。A.3趨勢面法趨勢面法采用多項(xiàng)式表示線或面，并按最小二乘法原理對樣本點(diǎn)屬性值進(jìn)行擬合，得到空間坐標(biāo)與樣本點(diǎn)屬性值的關(guān)系，并用于估算插值點(diǎn)的屬性值。具體見公式(A.6)。Z(x,y)=22n2^an.kx?-y1+e…………(A.6)式中：Z(x,y)——插值點(diǎn)的估計(jì)值；x,y——平面直接坐標(biāo)系中的橫、縱坐標(biāo)；n?——多項(xiàng)式階數(shù)；k,i——多項(xiàng)式的次數(shù)；a.;——樣本點(diǎn)值確定的系數(shù)?——趨勢面與真實(shí)面存在的隨機(jī)誤差。A.4泰森多邊形法泰森多邊形法(ThiessenPolygon)又叫作最近鄰點(diǎn)插值法，其基本原理是插值點(diǎn)的最佳值由最鄰近的樣本點(diǎn)產(chǎn)生。具體見公式(A.7)。式中：Z?——插值點(diǎn)的估計(jì)值；Z;——樣本點(diǎn)i的實(shí)際值。i點(diǎn)必須滿足公式(A.8)和公式(A.9)條件：do=min(do,do?,…,do)…………(A.8)d;;=√(x;-x;)2+(y;-y;)2…………(A.9)式中：do;——插值點(diǎn)與樣本點(diǎn)i間的距離；d;;——樣本點(diǎn)i與樣本點(diǎn)j間的距離；x;,y;——樣本點(diǎn)i平面直角坐標(biāo)系中的橫、縱坐標(biāo)；x;,y;——樣本點(diǎn)j平面直角坐標(biāo)系中的橫、縱坐標(biāo)。(資料性)留一交叉驗(yàn)證方法B.1交叉驗(yàn)證的一般定義交叉驗(yàn)證是把得到的樣本數(shù)據(jù)進(jìn)行切分并組合為不同的訓(xùn)練集和測試集，用訓(xùn)練集來訓(xùn)練模型，用測試集來評估預(yù)測模型的好壞。重復(fù)地使用數(shù)據(jù)，可以得到多組不同的訓(xùn)練集和測試集。交叉驗(yàn)證是在機(jī)器學(xué)習(xí)建立預(yù)測模型和驗(yàn)證模型參數(shù)時(shí)常用的方法。根據(jù)數(shù)據(jù)切分方式的不同，交叉驗(yàn)證主要分為下面三種方法。a)簡單交叉驗(yàn)證(CrossValidation)。隨機(jī)地將樣本數(shù)據(jù)分為兩部分，然后用訓(xùn)練集來訓(xùn)練模型，在測試集上驗(yàn)證模型及參數(shù)。一次驗(yàn)證之后再把樣本打亂，重新選擇訓(xùn)練集和測試集，繼續(xù)訓(xùn)練數(shù)據(jù)和檢驗(yàn)?zāi)Ｐ?，最后選擇損失函數(shù)評估最優(yōu)的模型和參數(shù)。b)S折交叉驗(yàn)證(S-FolderCrossValidation)。與“簡單交叉驗(yàn)證”不同，S折交叉驗(yàn)證會把樣本數(shù)據(jù)隨機(jī)地分成S份，每次隨機(jī)地選擇S-1份作為訓(xùn)練集，剩下的1份做測試集。當(dāng)這一輪完成后，重新隨機(jī)選擇S-1份來訓(xùn)練數(shù)據(jù)。若干輪(小于S)之后，選擇損失函數(shù)評估最優(yōu)的模型和參數(shù)。c)留一交叉驗(yàn)證(Leave-one-outCrossValidation),它是“S折交叉驗(yàn)證”的特例，此時(shí)S等于樣本數(shù)N,這樣對于N個(gè)樣本，每次選擇N-1個(gè)樣本來訓(xùn)練數(shù)據(jù)，留一個(gè)樣本來驗(yàn)證模型預(yù)測的好壞。此方法主要用于樣本量相對不多的情況，該方法也是使用最普遍的方法。B.2驗(yàn)證方法本文件所采用的留一交叉驗(yàn)證法對于空間插值模型進(jìn)行誤差估計(jì)，是每次移除1個(gè)數(shù)據(jù)位置，然后使用其他位置的數(shù)據(jù)來預(yù)測與其相關(guān)聯(lián)的數(shù)據(jù)，將預(yù)測值與實(shí)測值相比較獲得平均誤差、平均絕對誤差和均方根誤差。以平均誤差、平均絕對誤差和均方根誤差作為常用驗(yàn)證指標(biāo)，驗(yàn)證插值結(jié)果的準(zhǔn)確性。GB/T43605—2023(資料性)大氣腐蝕圖色標(biāo)C.1色標(biāo)色域及選擇原則大氣腐蝕圖的色標(biāo)宜在圖C.1中選取。當(dāng)大氣腐蝕性分級和金屬材料大氣腐蝕速率在GB/T19292.1劃分的基礎(chǔ)上需要更細(xì)劃分時(shí)，選擇的色標(biāo)宜連續(xù)漸變，依據(jù)大氣腐蝕性由弱至強(qiáng)，從左至右依次在橙色黑色綠色黃色橙色黑色綠色黃色藍(lán)色圖C.1色標(biāo)色域C.2色標(biāo)選擇示例在繪制大氣腐蝕圖時(shí)，在滿足C.1要求的前提下，可根據(jù)需要從表C.1中選擇色標(biāo)。表C.1色標(biāo)示例序號色標(biāo)填充顏色序號色標(biāo)填充顏色1255,199,020,86,239255.181.030,155,204255,165,040,205,152255,155,050,239,84255.133,06255.100.0784,255,0255,55,089204,255,3239,255,0255,255,0255,236,0255.218,0(資料性)大氣腐蝕圖圖面示例D.1大氣腐蝕圖色斑圖示例圖D.1和圖D.2分別是山東省大氣腐蝕性等級圖和金屬鋅大氣腐蝕速率圖繪制成色斑圖的示例。兩圖均是將大氣腐蝕性等級和金屬鋅的大氣腐蝕速率細(xì)分為13個(gè)等級繪制而成。詳細(xì)的繪制說明參見附錄E。山東省地圖來自山東省自然資源廳官方網(wǎng)站。山東省大氣腐蝕性等級圖圖D.1山東省大氣腐蝕性等級圖山東省金屬鋅大氣腐蝕速率圖陳度；37.186888省級行政中心數(shù)據(jù)源點(diǎn)界線圖D.2山東省金屬鋅大氣腐蝕速率圖上面兩圖中出現(xiàn)的小窗口，分別是讀者通過位置坐標(biāo)選擇之后顯示的該位置的大氣腐蝕性等級和金屬鋅的年腐蝕速率。D.2大氣腐蝕圖等值線圖實(shí)例圖D.3和圖D.4分別是對應(yīng)圖D.1和圖D.2的等值線圖。山東省大氣腐蝕性等級圖C3-3C4-1-C4-2—C4-3圖D.3山東省大氣腐蝕性等級圖山東省金屬鋅大氣腐蝕速率圖0訇@—21g/(m2·a)—0.7g/(m2a)—24g/(m3-a)—7g/(m2-a)—60g/(m2-a)-9g/(m2a)-11g(m2a)—180g/(m2-a)圖D.4山東省金屬鋅大氣腐蝕速率圖(資料性)E.1繪圖數(shù)據(jù)的說明附錄D中山東省大氣腐蝕圖的繪圖數(shù)據(jù)，來源于山東省100個(gè)站點(diǎn)的儀器檢測，其傳感器陽極材料為金屬鋅。用于山東省大氣腐蝕圖繪制的100個(gè)數(shù)據(jù)源點(diǎn)分布和腐蝕性等級見圖E.1。圖E.1用于山東省大氣腐蝕圖繪制的100個(gè)數(shù)據(jù)源點(diǎn)分布圖D.1繪圖數(shù)據(jù)源自2018年6月至2019年6月儀器檢測的第1年腐蝕速率數(shù)據(jù)，與同時(shí)期金屬鋅試樣實(shí)際暴露檢測結(jié)果對比，平均絕對誤差為0.2597g/(m2·a)。圖D.2的繪圖數(shù)據(jù)源自2018年6月至2022年6月儀器檢測的年化腐蝕速率。E.2大氣腐蝕性分級與顯示色標(biāo)選擇的說明基于山東省的繪圖數(shù)據(jù)，將大氣腐蝕性和金屬鋅大氣腐蝕速率都細(xì)分為13個(gè)等級，并依據(jù)附錄C表E.1大氣腐蝕性與金屬鋅大氣腐蝕速率分級及色標(biāo)序號大氣腐蝕性分級金屬鋅大氣腐蝕速率/[g/(m2·a)]填充顏色(RGB)色標(biāo)1C1-1reon≤0.70,0,2552C2-10.7<rcr≤50,155,2043C3-15<rn≤70,186,694C3-27<ro≤9153,255,05C3-39<ror≤11204,255,06C3-411<r-or≤13255,255,07C3-513<ror≤15255,230,0GB/T43605—2023表E.1大氣腐蝕性與金屬鋅大氣腐蝕速率分級及色標(biāo)(續(xù))序號大氣腐蝕性分級金屬鋅大氣腐蝕速率/[g/(m2·a)]填充顏色(RGB)色標(biāo)8C4-115<rmr≤18255.165,09C4-218<r-or≤21255,128,0C4-321<rcor≤24255,0,0C4-424<rcor≤30C5-130<rcor≤60CX-160<r-or≤1800,0,0E.3選取的地理空間插值模型在山東省大氣腐蝕圖繪制過程中，采用了反距離加權(quán)插值模型對數(shù)據(jù)源空白區(qū)域進(jìn)行插值。模型計(jì)算見公式(E.1)。Z?=Z′=?Z;λ;…式中：Z?——插值點(diǎn)的估計(jì)值；n——在估算中用到的控制點(diǎn)的數(shù)目；Z;——樣本點(diǎn)i的實(shí)際值；λ;——控制點(diǎn)的權(quán)重?？刂泣c(diǎn)的權(quán)重計(jì)算見公式(E.2)。式中：……………(E.2