无标题文档
呃…您访问的页面居然不存在！
您浏览的页面可能被删除，转移或暂时不可用
遇到此情况的概率和中彩票差不多哦！(据我们所知)基于TIMESAT的3种时序NDVI拟合方法比较研究_以藏北草地为例30
上亿文档资料，等你来发现
基于TIMESAT的3种时序NDVI拟合方法比较研究_以藏北草地为例30
第26卷第2期2011年4月；遥感技术与应用；Vol.26No.2Apr.2011；基于TIMESAT的3种时序NDVI拟合方法比较；)))以藏北草地为例；宋春桥1,2,柯灵红1,2,游松财3,刘高焕1,；(1.中国科学院地理科学与资源研究所,资源与环境；摘要:以藏北地区年MODIS1；1引言；卫星传感器所获取的NDVI(归一化差值植被指数
第26卷 第2期2011年4月遥 感 技 术 与 应 用 Vol.26 No.2Apr.2011基于TIMESAT的3种时序NDVI拟合方法比较研究)))以藏北草地为例宋春桥1,2,柯灵红1,2,游松财3,刘高焕1,钟新科1,2(1.中国科学院地理科学与资源研究所,资源与环境信息系统国家重点实验室,北京 .中国科学院研究生院,北京 .中国农业科学院农业环境与可持续发展研究所,北京 100081)摘要:以藏北地区年MODIS16d合成的NDVI时间序列为例,介绍了基于TIMESAT2.3软件的3种主要拟合算法)))非对称高斯函数(AG)拟合、双Logistic曲线(D-L)拟合和Savitzky-Golay(S-G)滤波法的基本原理和实现流程;重点从拟合重建NDVI时间序列对原始ND-VI值上包络线的拟合效果及保持原始高质量NDVI点值真实值的程度两个方面,分析比较3种算法的特点。结果表明:13种拟合算法均能不同程度提高整个区域的NDVI平均值,AG与D-L拟合法处理后的NDVI时间序列与原始NDVI曲线的整体特征较S-G滤波方法更加吻合;oAG与D-L拟合重建的NDVI时间曲线在生长季峰期高于上包络线,S-G滤波法处理结果低于上包络线,3种方法中AG拟合结果与上包络线最为接近;?在保持原始高质量NDVI值真实性方面,AG与D-L拟合法处理结果相似,除生长季曲线的峰期外,均优于Savitzky-Golay滤波法。该研究结论为基于NDVI时间序列进行陆地系统生态环境各方面研究中数据去噪预处理的方法选择提供参考。关 键 词:MODISNDVI;TIMESAT;拟合;噪声去除;时间序列数据;藏北中图分类号:TP79
文献标志码:A
文章编号:11)02-0147-091 引 言卫星传感器所获取的NDVI(归一化差值植被指数)时间序列数据能够精确地反映陆地生态系统植被的生长状态和季相、年际变化特征,已广泛应用于全球与区域生态环境变化监测与模拟、植被覆盖动态变化研究、植被物候特征识别与信息提取等许多方面[1-2]。目前常用的NDVI时间序列主要来自NOAAAVHRR、SPOTVEGETATION和TER-RA/AQUAMODIS等传感器。NDVI反映植被生长过程的时间序列曲线通常是连续平滑的,但由于传感器本身性能、数据传输过程失误、太阳光照角度、观测视角、地物双向性反射以及云、大气气溶胶等观测条件因时间而异,此外受地表水、冰雪等随机干扰因素,得到的观测值包含很多不可预测的噪音,得到的NDVI时间序列呈锯齿状的不规则波动,反映季节变化趋势不明显,使陆地生态系统监测、植被覆盖变化分析和信息提取等研究结论引进了诸多干扰[3-4]。尽管常用的NDVI数据集广泛基于最大值合成法(MVC)来降低这些噪声,但是一些假高值、合成期内连续的NDVI低值等噪声残差还是存在,为此研究发展了一系列ND-VI时间序列重建的方法。目前NDVI时间序列曲线重建的主要方法大致归纳为3类:阈值法、基于滤波函数的平滑法和曲线拟合法[5]。3种方法既可以基于原始的NDVI时间序列曲线重建,也可以对MVC处理后的数据集进行再处理。阈值法包括Viovy等数斜率提取法(BISE)和Lovell等[6]提出的最佳指在此基础上提[7]出的改进算法;基于滤波的平滑方法主要包括Savitzky-Golay(S-G)滤波法[8]、均值迭代滤波法(MVI)[9]、傅立叶变换[10]和时间序列谐波分析法(HANTS)[11];曲线拟合法以双Logistic函数拟合收稿日期:;修订日期:基金项目:藏北高原地区土壤水分与土壤温度时空变化模拟分析()。作者简介:宋春桥(1986-),男,湖南衡阳人,硕士研究生,主要从事遥感与GIS应用研究。E-mail:chunqiao_。法(D-LCF)[12]和非对称性高斯函数拟合法(AG-FF)[13]为代表。以上基本涵盖了当前遥感NDVI时间序列曲线重建的主要方法。目前,对于NDVI数据集预处理的方法除MVC得到大家基本认可,广泛应用于NOAA、MODIS、SPOT等遥感数据生产合成时间序列数据初级产品,其他方法分别用于国内外学者在全球或不同区域的研究,尚未达成共识。研究者在不断改进或者提出新的算法,希望能够克服现有方法本身的适应性或缺陷,有效滤除噪声、反映植被真实的生长状况和物候规律等。同时,部分国内外学者选择多种或两种不同的重建算法进行了比较研究,分析各自的优劣[14]。比如:Julien等[15]对其提出的IDR(IterativeInterpolationforDataReconstruction)方法与HANTS、D-LCF进行实例样点对比研究,IDR在多数样例中具有最好的效果,但是在强噪声中对NDVI低值过度拟合;Sakamoto等[16]沧江等主要江河的发源地,图1为藏北地理位置及草地类型空间分布。区域幅员辽阔,纬度位置为30b27c25d~35b39c13dN,83b41c14d~95b10c46dE。由于疆域跨度大,加之多变的垂直地带性影响,气候及水热条件差异较大,所以从东南向西北,植被类型分布与生物气候带相适应,依次出现亚高山疏林)灌丛草甸、高寒草甸、高寒草原、高寒荒漠草原四大基带及相邻的过渡亚带[21]。以日本为研究区,对小波变换和傅立叶变换方法进行了比较分析,结果表明,小波变换总体上优于傅立叶变换;Jennifer[17]选择加拿大阿尔伯塔中西部的洛矶山脉中一个典型区域为示范区,比较两种曲线拟合方法(D-L拟合法和AG拟合法)与4种滤波方法(S-G滤波法、MVI、ARMD3-ARMA5滤波及4253H两次滤波法)对NDVI时间序列重建效果,认为两种曲线拟合方法总体上优于4种滤波方法。在国内,仅有极少数相关的研究,且多针对上述各种方法的原理、优缺点进行定性分析。如顾娟[18]、李儒[19]对目前主要的几种NDVI时间序列方法进行概述;吴文斌等研究发现,S-G滤波法与非对称高斯拟合法对NDVI时间序列数据集拟合重建效果具有区域差异性。由于不同区域大气及地面条件、植被覆盖类型及作物制度的差异,以及缺乏定量比较手段,本研究以中国西部地区高寒和干旱环境具有代表性的藏北草地为例,利用非对称性高斯函数拟合、双Logistic曲线拟合和Savitzky-Golay滤波3种最常用的拟合重建算法对年MODIS16d合成的NDVI时间序列进行处理,对比分析拟合结果曲线与NDVI序列上包络曲线拟合效果以及保留高质量原始真值程度两方面的效果,为今后NDVI应用研究中选择时间序列拟合的方法提供参考。[5]图1 研究实例区域Fig.1 Thegeograhpicpositionofstudyinstance3 数据来源研究使用的实例数据选择EOS/Terra卫星的MODIS产品之一MOD13Q1v005,包括基于MVC方法16d合成的250m分辨率NDVI及其质量控制数据,该产品经过几何校正和大气校正。该时间序列数据集从美国NASALPDAAC工作组(LANDPROCESSESDISTRIBUTEDACTIVEARCHIVECENTER)网站下载(https://lpdaac.usgs.gov/lpdaac/products/modis_products_table),时间覆盖2007年1月~2009年12月,3a时间共69期影像。NDVI值标准范围为-1.0~1.0,MOD13Q1NDVI产品是-的DN值,-3000为填充值,从DN值转化成NDVI值的关系式为:NDVI=0.0001@DN
(1)此外,还需从MOD13Q1v005产品中提取datapixelreliability数据集,其空间、时间分辨率与ND-VI数据集匹配。该数据是对区域每个像元NDVI数据质量进行说明(如表1),分为5个不同等级,为该NDVI产品不同应用级别提供参考。本研究中,2 研究区域藏北地区位于青藏高原腹地,素有/江河源0、/生[20],SAT2.3平台对原始NDVI时间序列进行平滑处理,其关键步骤为设置参数,包括NDVI有效值域(Range,-0.05~1.00)、原始NDVI值权重(Weightvalue,详见表1)、噪声去除阈值(Spike,2)、滑动窗口大小(W-size,5-6-7)、拟合峰值参数(Alptitude,5)、迭代次数(Fittingsteps,3)等,基于S-G滤波法、非对称高斯函数拟合法及双Logistic曲线拟合法产生拟合后的NDVI时间序列数据集;?基于图1中的草地类型斑块,各种类型草地随机选取60个左右(共计300个)的年时间序列样点,得到各样点年原始及3种算法拟合后的NDVI时间序列数据集;?利用Matlab程序获取原始NDVI时间序列曲线的上包络线,并计算3种方法拟合后的NDVI序列对原始高质量NDVI值和上包络线的拟合效果,评价3种方法的优劣。拟合重建处理中的任何时刻像元NDVI值依据datapixelreliability数据设置相应的权重,图2中像元在3、4、6、10、16、19、20、21期为噪声点(圆圈标识)。图中3种拟合算法的加权与未加权处理的NDVI时间序列曲线显示,前者拟合后曲线均有不同程度的提升,因为受云掩(点10、16)及地面冰雪覆盖(其他圈出的点)的干扰造成NDVI值降低,使重建曲线对上包络线拟合不够,拟合后的NDVI最大值偏低。另外,由于春季返青之前和进入枯黄期阶段的连续噪声,生长季开始阶段和结束阶段的拟合曲线都发生右移,如果未进行不同权重处理,对于植被物候信息(如返青期、生长峰值、枯黄期等)的提取都将引入新的噪声。本文对3种拟合算法的比较研究都将引入数据质量权重系数进行处理。[22]NDVI时间序列拟合算法中填充NDVI点不参与计算,DN为0的NDVI点权值为1,DN为1的点权值设为0.8,而DN为2、3的点值权重为0.2。表1 MOD13Q1v005产品datapixelreliability数据集DN值说明Table1 DescriptionsofpixelDNvaluetypesof/datapixelreliabilitydataset0inMOD13Q1v005product像素值-10123含义填充值质量好质量较好冰雪覆盖云掩描述说明像元未被处理具有极大可信度质量比较可靠,具体可参考QA信息受冰雪覆被噪声影响受大气云量影响,目标被遮盖本研究所取权重01.00.80.20.24 研究方法4.1 基于TIMESAT的NDVI时间序列重建流程TIMESAT软件包是Jênsson与Eklundh共同开发的用于植被指数时间序列数据集重建及提取植被生长物候信息的程序包,以便于获取季节性植被生长规律信息等。该软件平台以S-G滤波法、AG拟合法及D-L拟合法为核心算法,能对NO-AA/NOAAPAL、TERRA或AQUA/MODIS等NDVI时间序列数据进行有效地处理提取信息。基于上述3种方法对年共23期NDVI时间序列重建效果对比分析处理主要包括以下4个步骤:1通过MRT软件对从NASA网站上获取的HDF格式MOD13Q1产品进行转换,并进行区域拼接、投影转换等预处理;o基于TIME-图2 单个像元NDVI年时间序列数据质量加权)未加权3种算法拟合效果对比Fig.2 Comparisonofunweightedandweighted4.2 3种拟合重建算法基本思想法(对称高斯函数拟合方法是一个从局部拟合到整体拟合的过程,使用分段高斯函数来模拟植被生长过程,最后通过平滑连接各高斯拟合曲线实现时间序列重建[13]。其过程大致分为区间提取、局部拟合和整体(2)双Logistic函数拟合法(D-L)。双Logistic期拟合处理前后的NDVI值,NDVIo、NDVIp分别为处理前后的年内NDVI时间序列的均值。通过计算每个样本像元的拟合重建前后的整年NDVI时间序列相关系数,统计所有像元相关系数的平均值,反映3种拟合算法保持原始NDVI时间曲线整体特征的能力。基于这3种主流算法的重建前后NDVI时间序列的拟合效果采用两种方式表示:一是拟合前后NDVI偏差分布的箱须图(Box-whiskerPlot)表示法,利用极小值、第一四分位数、中位数、第四分位数和极大值5个表示两组NDVI序列残差分布的离散程度及对称性等;二是回归估计标准差(如式3),描述重建后的NDVI年内时间序列与原始值之间的平均差异程度,表明重建前后NDVI值的代表性强弱,其值越小,拟合值的代表性越强。2(NDVIpi-NDVIoi)连接三步骤[22]。曲线拟合法是由Beck等[22]于2006年提出的一种新算法,与非对称高斯函数拟合法类似,双逻辑函数滤波也是一种半局部拟合方法。首先将整个时间序列中时间点对应的值按极大或极小值分成多个区间,分别对该区间进行DoubleLogistic函数局部拟合,其处理过程与AG方法类似。(3)Savitzky-Golay滤波法(S-G)。Savitzky-Golay滤波法最初由Savitzky和Golay于1960年提出,是一种移动窗口的加权平均算法,但其加权系数不是简单的常数窗口,而是通过在滑动窗口内对给定高阶多项式的最小二乘拟合得出。4.3 不同算法拟合效果比较方法在基于植被NDVI的各种陆地生态系统应用研究中,通常选取生长季峰期的NDVI值反映全年植被生长状况。因此有必要比较分析各种算法拟合重建后的NDVI值在7、8月份相对于初始NDVI时间曲线的上包络线的拟合效果。此外,利用ND-VI时间序列提取物候信息及应用全年或生长季期的NDVI平均状况的研究中,选择一种在去除原始NDVI数据集中噪声的同时,能够有效地保持质量好的NDVI点值的真实性的算法也十分关键。为了比较AG拟合、D-L拟合与S-G滤波法3种方法对高寒干旱环境的藏北草地MODIS16d合成的NDVI时间序列曲线拟合重建的效果,将着重从3个方面分析:1分析3种算法对整个区域NDVI均值在时间序列上改变的程度,以及拟合的时间曲线与原始NDVI曲线的相关系数;o比较各像元拟合后的NDVI时间序列曲线在成熟期对上包络线的拟合效果;?分析全年、生长季内、生长季开始和结束时段以及生长季峰期重建后的NDVI与高质量NDVI原始值的偏离程度。其中,相关系数用于反映拟合后的NDVI时间序列曲线保持原始植被生长季特征的能力,如式(2)所示:23(3)N其中:NDVIoi、NDVIpi分别是年内时间序列中第i期拟合处理前后的NDVI值,N表示质量控制文件中标识高质量(DN=0,1)的NDVI像元样本数量。RMSE=15 结果与分析5.1 全局特征拟合效果分析选择3种算法对年全区300个样点像元的NDVI时间序列曲线拟合重建,均将不同程度地去除各类噪声,通常情况一定区域内的像元受地面水分或冰雪覆盖以及大气云层、气溶胶影响造成的NDVI值降低的情况远多于/假0高值的情形。因此,各种重建方法将使区域内的像元NDVI整体提高,如图3所示,5条曲线分别为900个(300@3年)像元NDVI重建前、3种算法拟合后以及像元点原始NDVI序列的上包络线取均值的整年时间序列。其中,3种算法分别拟合重建后的NDVI值在全年23期较原始NDVI时间序列曲线均不同幅度提升。特别在植被返青(Number8~10)及进入枯黄期(Number18~21)时段,NDVI值重建后的整体水平相对提高11.9%~30.2%。这与实例区域的本身特征有关,这两个时段分别是草地覆被NDVI序列发生突变的/拐点0,同时也是藏北地区自然条件发生转变的时刻,土壤分别处于开始解冻和发生冻融作用的阶段,受冰雪覆盖变化影响较大;藏北高原的大部分地区气温、降水以及大气状况也3r=E(NDVIp-ii=1NDVIp)(NDVIoi-NDVIo)21(NDVIpi-NDVIp)@E(2)(NDVIoi-NDVIo)2i=1p建前后的相关系数。其中非对称高斯函数拟合方法重建的时间序列与植被原始NDVI时间曲线相关性最好,全年23期的NDVI点值及生长季内NDVI序列的拟合相关系数分别为0.791、0.803。D-L拟合法比AG拟合略低,相关程度分别为0.788、0.798;而S-G滤波法受噪声影响对整个时间序列曲线整体特征的保持性不及前两种方法。相对于全年的拟合效果,生长季内的时间序列采用任何一种算法重建前后相关程度都更好,因为实例区域在生长季以外的时间内地表覆被NDVI受噪声影响更严重,重建结果可靠性相对稍低。5.2 成熟期上包络线拟合效果分析生长季NDVI最大值(藏北草地时间为7~8月份)通常反映出一个地区植被生长最好的状况,而NDVI受噪声污染引起的NDVI值严重偏低将给研究结论带来偏差。图5反映了基于AG拟合、D-L拟合及S-G滤波3种算法重建的NDVI时间序列在13~15期(7月12日~8月28日)相对于原始NDVI序列上包络线的拟合效果比较。图3中显示,AG拟合与D-L拟合的结果超出上包络线曲线,S-G滤波法拟合重建结果低于上包络曲线。图5中不难发现,AG函数拟合与D-L曲线拟合结果相近,效果均优于S-G滤波法。种拟合重建方法对原始NDVI值的改变程度不一致。D-L曲线拟合与非对称高斯函数拟合方法拟合后的结果比较吻合,仅在生长季峰期D-L拟合法较AG拟合法的NDVI均值略高0.002左右。S-G滤波法的/过度0拟合特征对噪声比较敏感,其处理的ND-VI均值与前两种方法有较大的差异,在生长季开始和结束时间,拟合的曲线相对于原始及前两种算法重建的NDVI时间序列分别发生较大程度的/左偏0和/右移0,这对基于NDVI时间序列提取植被物候信息(如返青期或枯黄期)将引入新的噪声。在生长季峰期(Number13~15),3种算法拟合后的NDVI平均值较初始值均有提高,AG函数与D-L曲线拟合处理结果甚至超出原始NDVI曲线的上包络线,而S-G滤波法重建结果低于上包络线值。图3 3种算法重建前后全区样点NDVI均值时间序列曲线及原始上包络线比较Fig.3 Theenvelopecurve,meanNDVItime-seriesbeforeandafterthethreereconstructionalgorithms图4所示为900个全年及生长季内(5月9日~10月15日)NDVI时间序列基于上述3种算法拟合重图5 基于3种算法重建结果相对原始NDVI序列上包络线在成熟期的回归标准差Fig.5 TheS.E.ofregressionofNDVIsamplesintheseasonpeaksbetweentheenvelopecurveandreconstructedtime-seriesbasedondoublelogisticfunction,asymmetricGuassianfunction,andSavitzky-Golayfilteringmethods图4 基于3种算法的全年及生长季内NDVI序列的重建前后相关系数比较Fig.4 ThecomparisonofcorrelationcoefficientsofNDVItime-seriesbeforeandafterreconstructioninthewholeyearorgrowingseasonperiodbasedonthethreealgorithmabove5.3 原始曲线拟合效果分析图6所示为900个(年300个均匀分布于全区域的样点像元)全年NDVI时间序列样本点,基于上述3种算法拟合后的NDVI时间序列,包含各类专业文献、各类资格考试、专业论文、生活休闲娱乐、高等教育、外语学习资料、幼儿教育、小学教育、基于TIMESAT的3种时序NDVI拟合方法比较研究_以藏北草地为例30等内容。　
您可在本站搜索以下内容：
别人正在看什么？【焦作手机号码交易|焦作手机靓号转让求购】- 焦作百姓网
&&关注百姓微信公众号，尽享快捷刷新、信息推送、抽奖活动等众多微信功能您只需要：1 &用微信“扫一扫”右边的二维码2 &在微信中“关注”我们&|&|||||地区全焦作
QQ号 运营商不限
自定义区间：～元供求转让
只看有图只看今天发布的信息浏览方式：
1月7日解放10000元&/&移动1月5日武陟600元&/&其他1月2日温县 - 武德镇1600元&/&移动14年12月28日山阳30-50元14年12月19日解放400元&/&其他14年11月30日温县600元14年11月30日博爱400元14年10月22日博爱700元14年10月11日博爱 - 博爱中山路30000元&/&其他14年9月30日沁阳5000元&/&其他14年9月25日山阳1元14年8月22日山阳1500元14年8月21日沁阳200元&/&其他14年5月4日沁阳60元13年7月5日山阳50元周边城市或地区：1月7日郑州 - 中原 - 华山路99元1月7日郑州 - 金水 - 花园路18元1月7日郑州 - 管城2750元&/&移动1月7日郑州3000元&/&移动1月7日郑州 - 管城2750元&/&移动1月7日郑州 - 管城1100元&/&移动1月7日郑州 - 管城1100元&/&移动1月7日郑州5500元&/&移动1月7日郑州4950元&/&移动1月7日郑州4180元&/&移动1月7日郑州14300元&/&移动1月7日郑州19800元&/&移动1月7日郑州25300元&/&移动1月7日郑州82500元&/&移动1月7日郑州68200元&/&移动1月7日郑州33000元&/&移动1月6日洛阳 - 洛龙 - 伊川县2200元&/&其他1月6日新乡 - 封丘1月6日郑州3300元&/&移动1月6日郑州3300元&/&移动1月6日郑州3080元&/&移动1月6日郑州93500元&/&移动1月6日郑州30800元&/&移动1月6日郑州93500元&/&移动1月6日郑州3080元&/&移动1月6日郑州4070元&/&移动1月6日郑州4400元&/&移动1月6日郑州4400元&/&移动1月6日郑州4180元&/&移动1月6日郑州 - 二七4180元&/&移动1月5日洛阳 - 伊川 - 伊川明皋10000元1月5日郑州4950元&/&移动1月5日郑州4950元&/&移动1月5日郑州 - 二七4950元&/&移动1月5日郑州 - 二七4950元&/&移动1月4日郑州 - 金水 - 东风路500元&/&其他1月4日新乡 - 封丘1000元&/&其他1月1日郑州 - 二七 - 大学路5000元1月1日郑州25300元&/&移动1月1日郑州30800元&/&移动1月1日郑州2750元&/&移动1月1日郑州2420元&/&移动1月1日郑州3300元&/&移动1月1日郑州3520元&/&移动1月1日郑州6600元&/&移动1月1日郑州2200元&/&移动1月1日郑州2200元&/&移动1月1日郑州1100元&/&移动1月1日郑州2200元&/&移动1月1日郑州1760元&/&移动1月1日郑州5500元&/&移动1月1日郑州1980元&/&移动1月1日郑州7700元&/&移动1月1日郑州2200元&/&移动1月1日郑州1980元&/&移动1月1日郑州1980元&/&移动1月1日郑州2860元&/&移动1月1日郑州990元&/&移动1月1日郑州990元&/&移动1月1日郑州2860元&/&移动14年12月31日新乡 - 新乡1100元14年12月31日新乡 - 卫辉 - 后河元&/&其他14年12月31日郑州41800元&/&移动14年12月31日郑州1980元&/&移动14年12月31日郑州6050元&/&移动14年12月31日郑州1980元&/&移动14年12月31日郑州 - 管城132000元&/&移动14年12月31日郑州2200元&/&移动14年12月31日郑州2200元&/&移动14年12月31日郑州30800元&/&移动14年12月31日郑州60500元&/&移动14年12月31日郑州1980元&/&移动14年12月31日郑州24200元&/&移动14年12月31日郑州17600元&/&移动14年12月31日郑州74800元&/&移动14年12月31日郑州18700元&/&移动14年12月31日郑州990元&/&移动14年12月31日郑州14300元&/&移动14年12月31日郑州24200元&/&移动14年12月31日郑州990元&/&移动14年12月31日郑州11000元&/&移动14年12月31日郑州1980元&/&移动14年12月31日郑州 - 二七13200元&/&移动14年12月31日郑州17050元&/&移动14年12月31日郑州4400元&/&移动找不到合适信息？试试免费订阅吧，“焦作 手机号码”有新信息时会优先通知你哦！
赞助商链接
××手机App订阅扫码下载百姓网APP，随时订阅信息，还能私信沟通！微信订阅打开微信扫一扫二维码，一步轻松订阅！订阅后，“焦作 手机号码”有新信息时将通过微信优先通知您。没有微信？您可以使用QQ订阅&或&邮件订阅 &×手机App订阅扫码下载百姓网APP，随时订阅信息，还能私信沟通！邮件订阅请输入QQ号或电子邮箱，最新"焦作 手机号码"信息将发送到您的邮箱。点此通过微信订阅，消息接受更及时
反馈建议描述：请填写描述手机号：请填写手机号请填写手机号 只看有图40.}