DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 摘 要 乡村振兴战略五项“总要求”的关键环节是“产业兴旺”，没有产业就没有发展。有了产业的振兴，乡村的建设才有基础。如何破解产业振兴难题，明确回答“产业兴旺”到底依靠谁、为了谁，方法与路径怎么解决，显得格外重要又十分迫切。“公司 + 基地 + 农户”模式作为最常用的农业产业化经营模式能否响应时代召唤、担起“产业兴旺”大任，正面临大考。而首当其冲的是，模式核心竞争力何在、如何建构，是能否胜任的关键所在。本文试图以安徽霍山白芨中药材产业基地为例子，探索总结乡村振兴战略背景下“公司 + 基地 + 农户”模式的核心竞争力到底是什么，怎么建构，可选择路径在哪里，为实施乡村振兴战略提供路径选择、可操作方案参考。 关键词 :乡村振兴战略，“公司 + 基地 + 农户”模式，核心竞争力，建构路径 Copyright © 2019 by author and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 习总书记在党的十九大报告提出实施乡村振兴战略，并明确“产业兴旺、生态宜居、乡风文明、治理有效、生活富裕”的五项总要求，这项强国富民的重大战略能否实现，关乎复兴中华、实现伟大中国梦的成败 [ 1 ] 。 要想真正解决 “三农”问题，主要依靠农业产业化的推进，以实现农民致富目标。同时，农业产业化也是国际上促进农业发展、农民致富目标的主要路径。要加快农村发展就必须要建立一种科学的发展模式，结合各地区的实际情况，因地制宜地加快推进农业产业化，促进农业结构调整，建构起适应时代需要，相互贯通、协调融合、城乡一体的新型农业产业化发展模式，以高效农业产业化生产体系提升农业竞争力是乡村振兴的根本 [ 2 ] 。 时势造英雄，时势更呼唤英雄。“公司 + 基地 + 农户”模式作为最常用的农业产业化经营模式能否响应时代召唤、担起“产业兴旺”大任，正面临大考。乡村振兴战略背景下“公司 + 基地 + 农户”农业产业化经营组织模式的核心竞争力到底是什么，怎么建构，可选择路径在哪里，事关“产业兴旺”全局，更事关“乡村振兴”这个党和国家业已确定的战略大局。无论是学界还是实践层面，都必须给予正面的、肯定的、明确的回答。
system: [SYSTEM]

user: 2. 我国“公司 + 基地 + 农户”农业产业化经营组织模式发展历程与现状分析 (一) 发展历程 改革40年来，中国农村农业经营方式先后经历了四次制度创新：包括从合作化运动产生的人民公社集体经营、从“包产到户”改革形成的家庭经营、在市场经济大潮催生的合作经营以及在“工商资本下乡”过程中孕育的企业规模化经营等。“公司 + 基地 + 农户”模式正是在工商资本下乡过程中由农业产业企业作为龙头牵引发起，以基地组织为“中介”，让农户为生产主体而设立的新型农业产业化组织模式。它集企业经营、合作经营、集体经营、家庭经营于一体，促进农业经济从 “一元化”逐步转移到“多元化”，已成为构建我国“四位一体”现代农业产业化经营体系的有效途径 [ 3 ] 。 在国内，“公司 + 基地 + 农户”模式曾经历过两个阶段的实践。早在1904~1914年间，英美烟公司(British American Tobacco Company)在河南许昌、安徽凤阳、山东潍坊这三个地区利用买办建立了原料基地，创立了“公司 + 基地 + 农户”的生产模式雏形，但这种殖民经济模式通过剥削、压榨农户使公司获得了超额的利润，对当时的经济社会造成了严重的负面影响 [ 4 ] 。到20世纪80年代，该模式获得重启与发展，泰国正大集团在深圳建立了第一家合资饲料厂，在国内设有100多家下属企业；瑞士的雀巢公司与黑龙江省政府达成生产奶制品协议，建立了双城雀巢公司；此后，国内农业企业相继模仿，涌现了像四川希望集团、河南双汇集团、广东温氏集团等一大批“公司 + 基地 + 农户”模式的农业产业化龙头企业。 (二) 不足之处 一方面，实践过程中，考察我国“公司 + 基地 + 农户”模式几经波折的历史进程不难看出：第一阶段外国公司以牺牲农民的廉价劳动为代价的买办方式是注定要失败的；第二阶段是以农业垄断型企业有选择性区域投资开发设立生产基地，必定是具有得天独厚优势的少数农村被选入其中，大多数农村仍在等待和观望，亟待开发。 另一方面，从理论界对农业产业化组织模式和“公司 + 基地 + 农户”模式已有研究的梳理发现：要么是偏重于产业化组织模式某一方面的功能、意义和作用，或是农业产业化发展过程中出现的具体问等，而对“公司 + 基地 + 农户”模式作为一个整体性概念嵌入社会结构、建构核心竞争力、产生良性互动社会行为和发挥经济作用的研究明显不足；要么是侧重于法律、金融、产权制度、产业组织等角度的分析，经常疏忽研究的“盲区”：农村社会、社区、人与人关系、户与户关系、农民与农村基层管理组织关系等等缺少调查与分析，从经济社会学视角对模式应涵盖、包容、吸纳农村社会关系、社会文化、社会变化等因素，反应模式的完整性、系统性、稳定性和可持续性的分析存在严重不足；要么是从欧美、日本发达农业国家的农业产业化发展模式中得出有益的启示并提出相关建议，而忽略了国外农业发展起步早起点高、农民土地私有化制度、农业垄断资本实力雄厚的特点，其实践和理论对我们来说最多只能参考，不能照搬照抄，这方面没有引起足够的警惕和重视。 乡村振兴战略是个庞大而系统的工程，也是一项全新的伟大事业，需要高举习近平新时代中国特色社会主义思想旗帜，走有中国特色的农业产业化之路，需要创新实践发展，更要理论创新指导 [ 5 ] 。 (三) 新时代赋予“公司 + 基地 + 农户”模式的使命和任务 进入新时代，乡村振兴战略五项“总要求”的关键环节是“产业兴旺”，这不仅是在有选择性的农村实现，而是要全面实现。习总书记多次强调，“全面建成小康社会”，“最艰巨、最繁重的任务在农村，特别是在贫困地区”，“小康不小康，关键看老乡”，“小康路上一个都不能少”。只有困难群众都完全实现了脱贫的目标，这才是真正的小康。总书记还强调指出“坚决打赢脱贫攻坚战，脱贫路上一个都不能少”。乡村振兴战略也如是，一个村庄也不能落下 [ 6 ] 。 乡村振兴战略涉及中国农村人口6亿多，覆盖面广，经济社会发展基础不平衡，区域特点、自然资源、社会资源等多维领域千差万别，在此基础上实施农业产业化、实现“产业兴旺”就不单单是产业或经济问题，更多意义上应该是农村社会全面发展问题，这就必然要求统筹各种生产资源、生产关系，并兼顾各相关利益方，发掘农村农业发展新动能，在发展生产的基础上培育新产业、新业态和完善产业体系，使农村经济更加繁荣，从而改善民生让农民“生活富裕” [ 7 ] 。 为此，“公司 + 基地 + 农户”模式无论从理论还是实践，应积极响应时代召唤，焕发出新的生机和活力，担纲起“农业现代化”与“产业兴旺”重任，大胆的创新发展，砥砺前行，不辱使命。
system: [SYSTEM]

user: 3. “公司 + 基地 + 农户”模式的本质内涵、行动纲领及行为特征 (一) “公司 + 基地 + 农户”模式的概念 “公司 + 基地 + 农户”的基本内容，包括公司、农产品生产基地和农户家庭经营三个层次的经济联合体。以一个技术先进、资金雄厚的公司(主要是农产品加工和流通企业)为龙头，围绕大型龙头企业建立农产品基地，形成龙头带基地、基地连农户，产加供销一条龙，农工贸技一体化经营的模式 [ 8 ] 。 (二) 模式的本质、内涵和行动要领 模式本质是以国内外市场为导向，以提高经济效益为中心，立足当地自然资源、人文环境和农村社会现状，实行区域化布局、专业化生产、一体化经营、社会化服务、企业化管理，把产供销、贸工农、经科教紧密结合起来，形成一条龙的经营体制 [ 9 ] 。内涵是从而改造传统的自给半自给的农业和农村经济，使之和市场接轨，在家庭经营的基础上逐步实现农业生产的专业化、商品化和社会化。行动纲领是：确定主导产业，凸显区域优势，依靠龙头带动，发展农业规模经营，实行市场带动龙头，龙头带动基地，基地连接农户，发挥整体生产经营效率。 (三) 模式的社会行为特征 “公司 + 基地 + 农户”生产组织模式的社会行为特征是：通过类似“虚拟的联盟”的方式实现特色农业产业发展中“公司”、“基地”和“农户”三者间“权责明确、风险共担、利益共享”。首当其中的主要是将三者的优势充分发挥同时对三者存在的不足相互弥补，共同追求利益最大化。主要的形式“农户”“基地”和“公司”各自签订关于农产品的产销合同，通过合同等法律形式确定各自权责，“农户”专业生产将市场风险部分转移给“基地”和“公司”，“基地”和“农户”签订收购合同规定价格、数量和质量等扩大生产、提高技术形成规模化和专业化、“公司”按照订立的合同收购农产品加工农产品，保证“农民”利益不受损失最终实现农产品的市场价值。
system: [SYSTEM]

user: 5. “公司 + 基地 + 农户”模式核心竞争力建构路径的案例分析 (一) 安徽霍山白芨中药材产业基地概况介绍 南京厚芽生物医药科技有限公司是一家高科技研发型公司。公司志在解决中药资源的源头生产技术问题，践行与推广道地药材的质量理念，以“厚德载物，恪守品质”为宗旨，以发展、复壮道地药材产业为己任。 2017年开始，南京厚芽生物医药科技有限公司在安徽霍山全县乡镇建设20，000亩白芨规范化种植推广基地，由农户自己经营，公司给予提供良种、技术指导培训、市场销售的支持。 白芨属植物具备药用和观赏双重属性，具有极其广阔的市场空间。白芨产业的发展具有良好的经济社会效益和生态效益。对霍山区域农业经济结构调整，推动全县农业及旅游业发展具有积极意义。 (二) 模式核心竞争力建构的创新实践 通过实施中药材规范化、标准化、程序化、科学化的栽培与种植，逐步建立健全中药材现代化和产业化生产经营体系，在组织形式、经营模式、人才结构、高新技术、投资来源和市场体系等方面不断创新实践，为中药材产业的兴起和可持续发展创立了“公司 + 基地 + 农户”的示范化模型。 1) 基地实行“四统一”管理，即统一选择种子、统一栽培技术、统一防治病虫害、统一采挖加工等，以稳定和提高药材产量和质量。 2) 以霍山县政府牵头，以南京厚芽公司为技术依托，以霍山农业局技术服务推广管理站为主体，开设全县白芨中药材科技服务中心一处，基地技术服务中心4个。按照规范要求，开展药材规范化种植技术研究，从根本上保证药材质量的稳定、可控。 3) 建立白芨仿生种植示范点、示范基地约4000亩，该用地从农户手中通过“流转合同”取得，制订一个合同标准、一种流转价格、合约30年，投资主体是南京厚芽公司，自主经营管理，起龙头带动，引领产业发展潮流，作为示范效应而精心打造，也作为农业生态观光旅游景点布局。 4) 推广白芨大面积种植20，000亩以上，由全县农户自愿加盟、跟随学习、模仿种植，示范基地提供种苗种源、技术指导、咨询服务，土地经营权仍然归农户自己，公司提供技术咨询与服务，订单式承诺最低价收购，挨家挨户签订“承诺收购合同”明确质量标准、价格、验收规则、付款时间、交货地点。 5) 南京厚芽公司负责打通外部市场，锁定国内外大型中药生产企业(如葵花药业、卫材(中国)制药、云南白药等)的采购合同；并负责承担与南京农业大学、陕西师范大学科研开发责任，建立校企联盟合作，独家享有并使用各项高科技专利技术。使公司在白芨产业方面拥有核心知识产权、保持核心竞争力。 6) 霍山县政府专门成立“白芨中药材产业管理办公室”，统筹协调各个政府职能部门，具体负责管理全县白芨产业规划与发展。专门出台相关支持政策与管理办法，采取切实可行的具体措施激励农户种植。鼓励龙头企业-南京厚芽公司打造白芨产业链延伸的深加工产业园、集商贸、物流、电商交易平台于一体供应链物流园。 (三) 经营绩效评价与社会综合效益分析 1) 转变了生产经营模式，加快了当地白芨产业化经营步伐，提高了当地农民的收人和生活水平。该项目的实施有助于促进药典白芨等相关产业链的形成，带动了农户收入大幅度增长。预计实现基地乡镇农民人均增收5000元，及区域经济可持续发展。 2) 转变了农民传统的经营方式和种植习惯，降低了农民的生产风险。通过项目的实施，带动霍山县中草药种植规模超过20，000亩，解决5万人次的务工就业问题。促进农业生产结构调整和农民增收，加快新农村建设步伐。 3) 提高了农民的素质，培育了新型农民。由于白芨种植的特殊性，农户作为生产的主体起着不可替代的作用，农户是最基层、最基本的生产单位，农户的职责是按照公司对农产品产量、质量、规格的要求进行生产，在整个农业生产过程中起着决定性的作用。 4) 产生了良好的经济、社会和生态效益，构建起了多赢格局。南京厚芽公司作为“龙头企业”，联合县农业局及其技术管理站精心打造设立“公司 + 基地 + 农户”白芨种植经营组织模式，各方利益不但得到保障、还能稳定持续地增长，既保护了当地青山绿水的生态环境，又打造了新型产业、旅游观光业、生态农业和区域特色产品品牌。
system: [SYSTEM]

user: 6. 结论与思考 在当前我国大力推进乡村振兴战略的过程中，我们将“公司 + 基地 + 农户”农业产业化模式作为一个有机的整体进行研究，通过对于不同农村社区农民生活、观念、农村社会制度、社会行为、社会结构、社会管理、社会服务，以及农村自然资源、社会资源、社会资本、网络结构等等的调查，探索我国农村经济社会的历史、内容、地位、作用、动力、道路、方式、特点、成效和走向，分析由传统农业向现代化农业的道路、模式、规律等，归纳总结模式的内涵、基本要素以及新时代的作用与意义，致力于建构模式的核心竞争力模型。 从安徽霍山白芨产业采取“公司 + 基地 + 农户”的合作模式，迅速形成当地主导产业，形成专业化、社会化的生产体系，一举多赢的经营合作方式来看，该模式的核心竞争力建构是成功的。对于公司来说，确保了有足够数量和质量的中药材成品供应，以满足客户需要，收入和利润大幅提升。对于农户来说，有了稳定可靠的产品出路，以保证和提高自己的收入水平。对于地方政府来说，使得地方特色产业得到极大的发展，形成产业链，实现区域化经济规模化发展，带动并繁荣当地经济社会发展。 基于此，也本着现代农业产业“投资项目前决策科学化、模式创建系统化、合约商洽细致化、产品创新科技化、合作盈利持续化、企业管理社会化”的“六化”标准要求，从龙头企业、基地组织、农户、产品与服务、品牌与特色效应、模式的社会行为良性发展、科学管理等途径进行剖析，发掘农村农业经济发展新动能，可以大胆提出模式的核心竞争力的解决路径如下：龙头企业投资意愿诚恳而真实、资本实力雄厚，以及科技、技术、信息、市场、管理先进；龙头企业牵引带动下，模式的资源整合能力、社会关系组织融合能力得到有效发挥；公司、基地、农户利益链条链接紧密、分配合理、契约制度健全；模式下生产的各种产品、提出的各类服务能独具特色、具备高附加值、不可复制性，开创市场能力强劲，特色品牌效应显现，市场占有率比较高、利润率比较好，持续生存发展的市场竞争力强，等等。从某种意义上来讲，“公司 + 基地 + 农户”模式起到孕育一种更符合我国现阶段农村生产力发展模式的孵化器作用，成为推进我国农业产业化进程、促进现代化农业经济快速发展的不二选择。通过模式的核心竞争力模型建构，可使“公司 + 基地 + 农户”模式日趋完善、日臻成熟，在乡村振兴战略实践中可操作、可控制、可防范、可自救。
system: [SYSTEM]

user: 摘 要 “一带一路”是我国新时期应对国际贸易挑战的重大创新战略，贸易便利化则成为切实改善贸易环境、实现经济新增长的重要路径。通过系统构建贸易便利化评价体系，客观测算了“一带一路”沿线主要国家的贸易便利化水平，并利用拓展的引力模型实证探讨了贸易便利化程度对我国出口贸易的影响。研究表明，“一带一路”沿线各国的贸易便利化现状差异显著，但普遍存在一定的“短板”，提升空间较大。而随着区域经济一体化的深入发展，传统的关税措施对国际贸易的边际效应已逐渐降低，提升“一带一路”沿线各国贸易便利化水平将对我国的出口贸易产生显著的促进作用。 关键词 :“一带一路”，贸易便利化，引力模型 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 2017年2月22日，世界贸易组织(WTO)宣布，作为该组织成立20年来达成的首个多边贸易协定，在经过长达两年多时间的交付核准之后，《贸易便利化协定》 (以下简称《协定》)终于正式生效。根据国际机构测算，有效实施《协定》将使发达国家贸易成本降低10%，发展中国家成本降低13%~15.5%。《协定》实施最高可使发展中国家出口每年增长9.9% (约5690亿美元)，发达国家增长4.5% (约4750亿美元) 1 。而贸易便利化无疑是实现这一美好前景最为关键的要素。实际上，随着近年来经济全球化和区域经济一体化的深入发展，国家和地区之间的经贸关系不断开放，世界平均关税水平大幅下降，非关税贸易壁垒也在不断削减，传统的关税和非关税措施对国际贸易的边际影响效应逐渐降低。贸易便利化则成为各国改善贸易环境、降低贸易成本、提高贸易效率、促进经济增长的新路径，其对全球贸易的重要作用得到了较为一致的认可，因此也成为全球经贸领域新的谈判焦点。 与此同时，作为最大的发展中国家，我国自2001年加入WTO以来，经过多年的发展，已成功跻身世界前列的贸易大国。但近年来以美国为主导的TPP和TTIP谈判、欧盟的债务危机等对我国的贸易环境造成了巨大的外部压力。面对新时期全球经济增长格局和经贸环境的转折性变化，习近平总书记先后提出了共建“丝绸之路经济带”和“21世纪海上丝绸之路” (以下简称“一带一路”)的重大战略倡议。通过构建全方位对外开放的新格局和国际合作新架构，实施更为主动的开放性政策，以应对国际贸易挑战。作为核心目标之一，“一带一路”战略致力于亚欧非大陆及附近海洋的互联互通，促进货物与要素流动的速度与效率，最终实现沿线国家的互利共赢。 鉴于此，本文通过构建一套完整的贸易便利化测度体系，对“一带一路”沿线主要国家的贸易便利化水平进行客观度量，掌握其分布特征，并进一步结合扩展的引力模型实证探讨贸易便利化对我国出口贸易的影响路径，以期为我国“一带一路”战略的顺利推进提供一定的理论支撑和决策参考。
system: [SYSTEM]

user: 2. 文献综述 到目前为止，包括WTO、亚太经合组织(APEC)、经济合作与发展组织(OECD)等在内的国际组织和国内外学者都围绕贸易便利化进行了较为丰富的研究，其中，贸易便利化水平的测度和影响效应是主要的研究视角之一。 在Wilson等(2003)之前，关于贸易便利化的实证研究成果非常有限，Wilson等从两个方面弥补了已有文献的不足。首先选取了包括港口效率、海关环境、规制环境和电子商务等在内的4个指标，用于测算亚太地区国家的贸易便利化水平；其次，基于扩展的引力模型估计了上述四个要素对各国贸易流量的影响 [ 1 ]。2005年，Wilson等再次对全球范围内的75个国家的贸易流量和贸易便利化之间的关系进行了实证探讨，从港口效率、海关环境、规制环境和服务业基础设施等四个方面8个指标测度了各国的贸易便利化水平 [ 2 ]。Shepherd和Wilson (2008)回顾了东南亚国家的贸易便利化演变历程，并同样基于引力模型验证了东南亚地区的贸易流量受交通基础设施和通信技术的影响最为显著 [ 3 ]。Felipe和Kumar (2010)则采用世界银行的物流绩效指标LPI来测度贸易便利化水平，进而结合引力模型检验了双边贸易流量与贸易便利化之间的关系，并对中亚国家贸易便利化水平的提升可能带来的贸易增长进行了估计 [ 4 ]。Moise和Sorescu (2013)基于OECD的研究构建了一个包含16个指标的评价体系(TFI)，用于对发展中国家的贸易便利化水平进行测度，并探讨了TFI指标对发展中国家进出口贸易的影响 [ 5 ]。Rudjanakanoknad等(2014)从货运代理人的角度构建了贸易便利化测度指标，主要包含港口基础设施、港口管理水平、海关程序和海关执法透明度四个方面，并对泰国的7个国际港口的贸易便利化水平进行了评估 [ 6 ]。 近年来，国内关于贸易便利化的研究也得到了快速发展。曾铮和周茜(2008)较早地构建了贸易便利化评价指标体系，并基于层次分析法对全球范围内的49个代表性国家的贸易便利化水平进行了测算，通过引力模型验证了贸易便利化水平是影响我国与这些贸易伙伴国之间双边贸易额的重要因素 [ 7 ]。谢娟娟和岳静(2011)、孙林和徐旭霏(2011)、方晓丽和朱明侠(2013)等都测算了中国及东盟国家贸易便利化水平，并利用引力模型验证了贸易便利化程度对我国出口的影响 [ 8 ] [ 9 ] [ 10 ]。“一带一路”倡议提出之后，王彦芳(2015)、崔日明和黄英琬(2016)尝试构建了一套评价指标体系，用于分析“一带一路”沿线不同国家的贸易便利化现状 [ 11 ] [ 12 ]。孔庆峰和董虹蔚(2015)、张亚斌等(2016)则在贸易便利化水平测度的基础上，进一步探讨了贸易便利化水平的改善对我国贸易潜力的影响 [ 13 ] [ 14 ]。此外，谭晶荣和潘华曦(2016)则重点研究了“一带一路”沿线国家的贸易便利化程度对我国农产品出口的影响 [ 15 ]。 通过对已有文献的梳理我们不难发现，在构建贸易便利化综合评价指标体系时，目前国内外主流的做法均是以Wilson等(2003， 2005)所提出的框架为基础，根据研究目的的不同，进行评价指标的调整和选择。关于贸易便利化水平对地区贸易流量等的影响，也都较为一致地采用拓展的贸易引力模型进行实证探讨。但由于国内专门针对“一带一路”区域贸易便利化的研究集中于最近两年出现，成果积累相对较少，在评价体系的构建、研究对象的范畴、研究方法的选择等方面尚存在进一步探讨的空间。
system: [SYSTEM]

user: 4. 我国出口贸易影响因素的实证分析 (一) 模型设定 Tinbergen (1962)和Poyhonen (1963)最早将引力模型应用到国际贸易领域，之后诸多学者对该模型中的因素进行了补充和拓展。如前文所述，大量研究成果表明，采用贸易引力模型分析一国贸易流量影响因素具有较好的优势和适用性，能够在一定程度上对现实问题进行解释和模拟，自20世纪60年代以来在国际贸易研究中得到了广泛应用。 标准的引力模型通常包括三类解释变量：衡量经济规模的解释变量，如一国GDP和人口总量；衡量地理位置的解释变量，如国家之间的距离；贸易政策虚拟变量，如贸易双方是否同为某国际组织成员国。在此基础上，本文进一步在模型中引入关税水平和贸易便利化水平两个变量，用于探讨贸易便利化程度对一国贸易流量可能产生的影响，并与传统的关税措施效应进行比较。由于本文着重探讨我国出口贸易的影响因素，实证分析选取的数据为2016年的截面数据，我国自身的GDP和人口总量为常数，对回归结果没有影响，因此不作为变量纳入模型。模型最终的基本表达式设定为 3 ： ln E X P i j = a 0 + a 1 ln G D P j + a 2 ln P O P j + a 3 ln D I S i j + a 4 A P E C i j + a 5 ln T a r i f f j + a 6 T F I j + u i j (2) 其中，因变量EXP ij 代表我国(第i国)对其他各国(第j国)的出口额。GDP j 和POP j 分别代表第j国的国内生产总值和人口总量；DISij是我国与第j国之间的距离；APEC ij 是虚拟变量，代表我国与第j国是否同属于亚太经合组织，如果是，取值为1，否则为0；Tariff j 是第j国的关税水平；TFI j 是前文测算得到的第j国的贸易便利化水平。上述变量的数据来源如表7所示。 Table 7 Variables in the model and data sourc 变量 含义 数据来源 EXP ij 我国对j国的出口额(单位：美元) 联合国COMTRADE数据库2016 GDP j j国的国内生产总值(单位：美元) 世界银行数据库2016 POP j j国的总人口数(单位：人) 世界银行数据库2016 DISij 我国与j国之间的距离(单位：千米) CEPII数据库 APEC ij 我国与j国是否同属于亚太经合组织 WTO官方网站 Tariff j j国的关税水平 全球竞争力报告2016-2017 TFI j j国的贸易便利化水平 前文测算结果 表7. 模型变量说明及数据来源 (二) 实证分析结果 本文收集整理了2016年各国的相关数据，由于个别国家贸易数据缺失，最终采用39个国家的截面数据进行回归分析。SPSS 19.0输出结果显示，本文建立的模型(2)拟合优度较好，调整后R 2 达到了0.91，整个回归方程通过了显著性检验，D-W值接近于2，残差分析结果也显示其基本符合正态分布假定(限于篇幅，此处未给出所有输出结果)，模型整体回归效果良好。 表8给出了主要的模型估计结果，可以看到，除了变量APEC和Tariff，其他自变量的回归系数均较为显著。 1) 国内生产总值反映了一国的经济规模和潜在市场需求，其对双边贸易的重要引力作用在已有文献中得到了较为普遍的认可。本文的最新数据分析再次验证了，作为我国贸易合作伙伴，如果“一带一路”沿线国家的国内生产总值增加1%，我国对其出口贸易额将增加0.406%。 Table 8 Estimation results of mode 自变量 回归系数 t统计量 (常量) 7.871 *** 2.776 LnGDP 0.406 ** 2.116 LnPOP 0.498 ** 2.365 LnDIS −0.596 * −1.979 APEC 0.106 0.286 LnTariff −0.119 −0.965 TFI 2.141 * 1.885 判定系数R 2 0.961 调整后R 2 0.910 D-W值 1.864 表8. 模型估计结果 注：***、**和*分别表示在1%、5%和10%的显著性水平下通过检验。 2) 关于人口对国际贸易的影响效应，理论上有着不同的解释，不同学者的实证分析结果也不尽相同。在本文的样本范围内，进口国的人口总量对于我国出口额具有显著的正向影响。在其他因素不变的情况下，如果进口国的人口数增加1%，我国对其出口贸易额将增加0.498%。这可能与“一带一路”沿线的国家构成有较为紧密的联系，当人口增加时，其市场需求规模和多样化都在增加，从而有利于促进国际贸易活动的开展。 3) 地理距离对我国的出口贸易具有较为显著的阻碍作用，显然，两国之间的距离越远，运输成本越高。实证结果显示，与“一带一路”沿线国家的地理距离增加1%，我国的出口贸易额将降低0.596%。这也进一步验证了作为贸易便利化测评的构成指标，进口国的口岸与物流效率是影响我国出口贸易的重要因素之一。可以预见，随着“一带一路”战略的推进，亚欧非大陆的设施建设日益完善，地理距离的这一约束作用将逐渐减弱。 4) APEC变量的回归系数为正，说明如果进口国与我国同属于亚太经合组织，将对我国出口贸易产生一定的积极影响，但该影响并不显著。这主要是由于“一带一路”沿线区域覆盖面广，区域内的国家分别隶属于包括东盟、欧盟等在内的不同经济合作组织，而本文是以我国的出口贸易为研究视角，因此并未将其他一体化组织的影响效应纳入考察范围。此外，区域经济一体化组织主要是通过取消区域内关税壁垒等政策来促进区域贸易发展，而随着世界范围内平均关税水平的大幅下降，关税政策对国际贸易的边际影响效应已逐渐降低(后文的分析验证了该结论)，因此该变量的正向影响也趋于弱化。 5) 关税水平的回归系数为−0.119，表示进口国的关税水平每降低1%，我国对其的出口贸易额将增加0.119%。显然，相对其他变量，传统的关税政策对于出口贸易的影响较弱，且并不显著。毫无疑问，随着“一带一路”区域经济合作的日益深化，各国之间的经贸关系不断开放，关税壁垒的阻碍作用将逐渐消除。 6) 贸易便利化作为近年来全球经贸领域新的关注焦点，其对我国出口贸易的积极影响再次得到了验证。数据分析结果显示，如果“一带一路”沿线国家的贸易便利化水平每提高0.1分，我国对其的出口贸易额将增加0.214%，且影响显著，贸易便利化已成为双边贸易的重要引力来源之一。
system: [SYSTEM]

user: 5. 结语 作为我国新时期应对国际贸易挑战的全方位扩大开放布局的重要组成部分，“一带一路”战略以创新的视野致力于培育全球经济新增长点，实现沿线国家的互利共赢。本文通过系统构建贸易便利化评价指标体系，客观测度了“一带一路”沿线各国的贸易便利化现状。由于地域覆盖面广，各国经济发展程度参差不齐，其贸易便利化水平差距显著。但普遍来说，各国均存在不同的“短板”，提升空间较大。而实证分析表明，“一带一路”沿线国家贸易便利化程度的改善对我国的出口贸易有显著的促进作用，尤其是在关税下降空间不断压缩的全球背景下，贸易便利化建设的积极效应必将日益凸显。 在“一带一路”战略的推行过程中，我国一方面应以初期阶段各类口岸基础设施的互联互通为突破口，发挥我国在口岸与物流效率方面的相对优势，加大对“一带一路”沿线发展中国家和落后地区基础设施建设的援助，弥补海关与边境管理缺陷；另一方面，通过区域合作政策的制定，督促包括我国自身在内的沿线各国加快实现政府规制环境、金融与电子商务等长期建设目标，为贸易便利化水平的综合提升，进而为我国出口贸易潜力的充分释放创造良好的“软实力”支撑。
system: [SYSTEM]

user: 摘 要 在脑机接口技术研究中，为有效提高脑电信号分类正确率，提出了一种样本容量扩增方法。该方法将一组脑电数据的所有电极数据值增加一倍，然后与原始数据一起作为样本数据进行训练，这种方法不但增加了样本容量，而且进一步提高了训练的收敛速度及测试结果的可靠性。用贝叶斯线性判别法对脑机接口行列闪范式下公共数据扩增前后的样本进行训练和测试，在两种类别的数据中将数据在直线上进行投影。经研究发现投影后每一种样本扩增数据的投影点比原始数据的投影点更接近，两种类别的数据的类别中心点也更大，分类效果更好，平均正确率相比原始公共数据平均正确率有了显著性的提高。研究表明这种用于脑机接口的训练样本集扩增方法能显著提高脑电信号分类正确率。通过实例验证了所提出的样本扩增方法的有效性，为样本容量不足提供了可行的解决方法。 关键词 :脑机接口，样本容量，贝叶斯线性判别法 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 科技进步给人们的日常生活带来了巨大的变化，不但方便了健康人群，同时也提高了残疾病患的生活质量。例如一些头脑健全但由于多种疾病引起的肢体瘫痪患者，他们无法与外界正常交流 [ 1 ] ，于是在过去几十年中，国内外众多研究者相继开始了脑机接口(brain-computer interface， BCI)技术的研究。与任何通信与控制系统一样BCI系统是由输入，输出及前者转换为后者的翻译算法几部分构成的一种全新人机交互方法 [ 2 ] ，而通用的BCI系统模型通常由数据采集和数据预处理(输入部分)、特征提取和分类识别(翻译算法)、控制接口及外围设备(输出部分)构成 [ 3 ] 。BCI系统通过实时提取反映大脑不同状态的脑电信号来实现控制，本文用到的目标脑电信号是性能较好的P300事件相关电位。BCI通过采集患者大脑皮层的脑电信号(electroencephalogram， EEG)输入计算机经过一系列预处理、特征提取和分类识别将其转化为可被外界设备识别的操作指令来确定患者的意图 [ 4 ] 。在这个过程中，对患者目标脑电信号P300的分类识别至关重要，但是通常采集到的患者脑电数据有限使得样本容量较小同时造成了P300分类正确率较低的问题 [ 5 ] 。针对此问题，本文列出了一种样本容量扩增方法，通过增加脑电信号中各电极的数据值形成新的一组数据与原始数据一起作为样本用Hoffmann等人 [ 6 ] 提出的贝叶斯线性判别方法(Bayesian Linear Discriminant Analysis， BLDA)进行训练和测试以达到提高分类正确率的目的。BLDA是一种有监督的线性分类器 [ 7 ] ，此方法的计算原理简单、稳健性较好并且具有良好的泛化能力，并且它还能解决噪声引起的过拟合问题 [ 8 ] ，因此本文我们用它来解决脑电信号的分类问题。脑电信号分类正确率的提高，不但对于实时在线脑电信号分类具有推动意义，而且使BCI系统对于肢体瘫痪患者对外界设备控制意图有一个更精确的识别，使得患者更有效的对外围设备进行操作和控制，进一步提高了残疾病患的生活质量方便了他们的日常生活。
system: [SYSTEM]

user: 2. 方法 2.1. 实验一 2.1.1. 行列闪范式 行列闪范式也被称为P300范式，它是由Farwell和Donchin于1988年提出的包含6 × 6字符矩阵的刺激范式 [ 9 ] ，范式中包括10个数字字符(0~9)和26个英文字母字符(A~Z)，范式背景为黑色，如图1所示。实验中所有字符按行和列随机高亮由灰色变为白色，当包含目标字符的行列高亮时，称之为靶刺激且其出现概率为1/6；其他不包含目标字符的行列高亮时称之为非靶刺激且其出现概率为5/6。P300电位信号的识别能够使得目标字符得以正常输出。实验过程中，被试注视目标字符并默记目标字符高亮的次数，小概率的靶刺激会诱发出P300电位但是非靶刺激不会，实验通过采集被试的脑电信号进行分类识别，检测出诱发P300电位的行和列确定目标字符的坐标即可唯一确定目标字符。 图1. 行列闪范式 2.1.2. 样本容量扩增 本文中我们使用BCI Competition III中在行列闪范式下用BCI 2000采集的Data set II数据，其行列闪频率为5.7 Hz且每个目标字符重复15次，首先对采集的波形信号采样，在64个电极上使用240 Hz的采样频率，通过0.1 Hz至60 Hz的带通滤波器 [ 10 ] ，滤出最为相关的脑电信号。Data set II中包括IIIA和IIIB两个被试的公共数据，将IIIA和IIIB的64个电极上matlab数据值都增加一倍，其中Pz电极上的P300和非P300数据在前八个采样点的数据变化如表1所示。 所有数据值增加一倍后形成新的一组实验数据，与原始数据合在一起作为训练数据进行训练，扩增数目变化如表2所示。样本数据发生变化后，其生成的64个电极上的脑电信号波形也产生了对应的变化，其中Pz电极上样本数据增加前后P300和非P300的波形对比如图2所示。 Table 1 The comparison of data changes in the first ten sampling points on the Pz electrod P300 数据增加前 14.07 14.84 14.84 16.37 15.86 10.76 4.90 5.41 数据增加后 28.15 29.67 29.67 32.73 31.71 21.52 9.81 10.83 非P300 数据增加前 3.63 0.06 −4.27 −5.29 −2.74 −1.21 −2.74 −1.72 数据增加后 7.26 0.13 −8.53 −10.57 −5.48 −2.42 −5.48 −3.44 表1. 在Pz电极上前八个采样点的数据变化对比 Table 2 The comparison of the number of samples before and after amplificatio BCI Competition III 字符数目 原始样本数目 扩增后样本数目 训练数据 100 18，000 36，000 测试数据 85 15，300 15，300 表2. 每个被试样本扩增前后数目变化对比 图2. Pz电极上靶刺激与非靶刺激数据值变化前后波形对比图。(a) 靶刺激与非靶刺激数据值变化前波形图；(b) 靶刺激与非靶刺激数据值变化后波形图 2.2. 实验二 2.2.1. 单闪实验范式 实验二作为验证实验，实验中我们不再使用行列闪范式而改为单闪范式，每个字符单独的随机进行闪烁，这样靶刺激出现的概率就会大幅度降低且目标字符闪烁时间间隔也会大幅度增加，这样就避免了传统行列闪范式中存在的“双闪”(同一个目标字符的行或列在短时间间隔内两次闪烁)问题，实验所用单闪范式如图3所示。 2.2.2. 样本容量扩增 我们在实验中采集了9名被试的脑电信号数据，经过预处理后保存为matlab格式。用本文提出的样本扩增方法进行样本扩增后样本数目变化如表3所示。样本数据发生变化后，其生成的各个电极上的脑电信号波形也产生了对应的变化，其中Pz电极上样本数据增加前后P300和非P300的波形对比如图4所示。 图3. 单闪范式 Table 3 The comparison of the number of samples before and after amplificatio 实验数据 字符数目 原始样本数目 扩增后样本数目 训练数据 15 2700 5400 测试数据 5 900 900 表3. 每个被试样本扩增前后数目变化对比 图4. Pz电极上靶刺激与非靶刺激数据值变化前后波形对比图。(a) 靶刺激与非靶刺激数据值变化前波形图；(b) 靶刺激与非靶刺激数据值变化后波形图
system: [SYSTEM]

user: 4. 讨论 根据两个实验的分类结果来看，样本扩增后的训练比原始数据的训练效果更好。从样本数据方面来说，增加原始数据的数据值一倍后的新数据与原始数据一起作为样本训练的时候数据平均值增加了一半，是原始数据平均值的1.5倍，同时训练得到的数据最大值和最小值也发生了变化，均是原始数据最值的1.5倍，如图7所示。其中a和c是用原始数据训练得到的靶刺激和非靶刺激数据最值在测试数据上对比情况，b和d是用扩增样本训练得到的靶刺激和非靶刺激数据最值在测试数据上对比情况。由图7可以看到，同一组测试数据在最值范围内的数据情况有很大的不同，使用扩增样本数据训练得到的最值间的采样点数据明显的多于用原始数据训练得到的最值间采样点数据，同时最值范围较大，落在范围内的样本点数据变多，使得靶刺激和非靶刺激的数据特征更加明显，有助于后续分类算法对靶刺激和非靶刺激的分类识别。 图7. Pz电极上样本扩增前后训练得到的最值在测试数据上的表示对比图。(a) 训练原始样本得到的最值在非靶刺激测试数据上的表示；(b) 训练扩增样本得到的最值在非靶刺激测试数据上的表示；(c) 训练原始样本得到的最值在靶刺激数据上的表示；(d) 训练扩增样本得到的最值在靶刺激测试数据上的表示 BLDA是一种有监督的学习方法，也就是说它的数据集每个样本都有类别的输出。在两种类别的数据中BLDA将数据在直线上进行投影，投影后每一种数据的投影点尽可能的接近，而两种类别的数据的类别中心点尽可能的大才会使得分类效果较好。如图8所示，我们把目标样本和非目标样本用灰色和黑色表示，用BLDA将两种类别数据投影到一条直线上。从直观上可以看出，右图比左图的投影效果要好，因为右图的黑色和灰色数据较为集中，且两种类别之间的距离也更明显。从BLDA原理来看，假设我们的数据集 D = { x 1 ， x 2 ， ⋯ ， x m } ，其中任意样本x i 为第i个数据。我们定义 X j ( j = 0 ， 1 ) 是第j类样本的集合， N j ( j = 0 ， 1 ) 是第j类样本的个数，而 ∑ j ( j = 0 ， 1 ) 是第j类样本的协方差矩阵， μ j ( j = 0 ， 1 ) 是第j类样本的均值，J(w)是BLDA投影到直线w后的损失函数。 μ j 的表达式为： μ j = 1 N j ∑ x ∈ X j x ( j = 0 ， 1 ) (16) ∑ j 的表达式为： ∑ j = ∑ x ∈ X j ( x − μ j ) T ( j = 0 ， 1 ) (17) J(w)的表达式为： J ( w ) = | μ ¯ 0 − μ ¯ 1 | 2 ∑ ¯ 0 2 + ∑ ¯ 1 2 (18) 我们实验中所用到的是两类数据，一类是目标数据，另一类是非目标数据，因此我们需要将两类数据投影到一条直线即可。假设投影的直线是向量w，则对每一个样本x i 在直线w上的投影为w T x i ，我们把目标样本和非目标样本两种类别的中心点设为μ 0 和μ 1 ，它们在直线w上的投影分别为w T μ 0 和w T μ 1 。从图8中可以看出，样本扩增后样本容量变大，数据平均值变大，由公式(16)可知样本均值向量μ j 值变大，两个类别数据的中心点μ 0 和μ 1 间距同时变大，这样使得两种类别的中心间距变大。样本数据均值变大，由公式(17)可知，协方差矩阵中的每个协方差数据值变小，说明同一类别数据的投影点更加接近。此时公式(18)中分母上每一个类别的方差之和变小，分子上两个类别各自的中心点间距的平方变大，相当于最大化了J(w)，因此即可求出比原始数据更优的w。这种效果体现在分类正确率上就是BLDA对扩增后的样本训练和测试所得平均正确率均高于对原始数据的训练和测试正确率。 图8. BLDA对扩增前后样本点在直线上的投影结果对比图。(a) BLDA对原始样本点在直线上的投影结果图；(b) BLDA对扩增样本点在直线上的投影结果图
system: [SYSTEM]

user: 5. 结束语 本文所提出的基于脑机接口技术的样本容量扩增方法，利用增加样本数据值后的新数据与原始数据一起作为训练数据，增加了样本数据的最值、均值，降低了样本数据点的离散程度，使得损失函数达到最大化，得出更优的BLDA投影直线，对测试数据分类时提高了平均分类正确率。把这种方法应用在实时分类上会使得BCI系统对于肢体瘫痪患者对外界设备控制意图有一个更精确的识别，同时让患者能够更有效的对外围设备进行操作和控制。在脑机接口技术方面，现有的分类算法准确率还是较低，不能达到实时准确分类的效果，所以在投入到实际应用中还是有一定的困难。因此未来可以从分类算法方面着手，寻求训练模型更加稳定、分类性能更好、准确率更高的脑电信号分类算法。
system: [SYSTEM]

user: 摘 要 基于机器视觉的工件瑕疵检测技术已经成为了重要的发展趋势，然而在实际的检测环境中，光照不均问题对检测结果带来了非常大的影响。因此，光照校正处理成为了机器视觉检测的关键步骤，但是现有的光照校正算法往往存在图像细节丢失的问题，直接影响了检测效果。针对这一问题，提出了一种改进型的Retinex光照校正算法，采用中值滤波代替了高斯环绕函数，从而实现了光照分量的估计，并使用信息熵函数完成了对原图光照均匀程度的判定。同时，基于所提算法搭建了一套轴承滚子瑕疵检测系统，并对所提算法的性能进行验证。实验结果表明，与现有的Retinex系列光照校正算法相比，该算法能够更好地保留工件的图像细节，从而实现更优的检测准确率。 关键词 :光照校正，中值滤波，机器视觉，缺陷检测，Retinex原理 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 轴承是现代机械设备中的重要零部件，也是我国重点发展的战略性基础产业之一 [ 1 ] [ 2 ] [ 3 ] 。“十二五”计划明确指出，轴承是加快振兴装备制造业的16个关键领域的核心基石 [ 4 ] ，它在多个关键领域有着广泛的应用。作为轴承核心零件之一的轴承滚子，其质量直接决定了轴承的质量。 目前，在轴承滚子质量检测方法中，基于机器视觉的检测方法因为速度快，准确率高等优势，已经成为了轴承滚子质量检测的主流方案 [ 5 ] 。基于机器视觉的轴承滚子检测系统通过工业相机来完成待检物表面图像的采集，随后将图像传送到图像处理单元，并完成图像中关键信息的提取和判别，进而实现对轴承表面瑕疵的检测 [ 6 ] 。 在实际场景下工作的轴承滚子检测系统中，光照不均匀问题一直是影响其检测性能的难题。由于轴承滚子的材质特性，往往会导致与滚子本身纹理相似的瑕疵难以被检测出来。另外，光照不均现象会对瑕疵的判别带来影响，在轴承滚子图像中，光照不均区域大小、灰度变化的随机性会降低瑕疵判别标准的准确性，从而降低系统的检测准确率。因此，需要在传统的基于机器视觉的轴承滚子瑕疵检测系统中加入光照校正算法，以此来消除光照影响，提高系统检测准确率。 目前应用较多的光照校正方法包括直方图均衡法 [ 7 ] 、基于照明–反射模型的同态滤波法 [ 8 ] 、基于Retinex原理的光照校正算法 [ 9 ] 以及图像边缘增强法 [ 10 ] 等。其中，基于Retinex原理的光照校正算法以处理灰度图见长 [ 11 ] 。基于Retinex原理的光照校正方法通过使用高斯环绕函数完成了光照分量的估计，而高斯环绕函数的物理意义与均值滤波类似，目标是进行图像平滑，其平滑效果相对于非线性滤波来说更加缓和。因此基于Retinex原理的光照校正算法会将图像中的细节估计为光照分量，从而容易造成光照校正后图像细节的丢失。尤其对于轴承滚子瑕疵检测这种对图像精细度要求很高的应用场景来说，影响更大。 针对当前基于Retinex原理的光照校正算法在轴承瑕疵检测应用场景中的缺陷，本文提出了一种改进型Retinex光照校正算法，采用中值滤波代替了高斯环绕函数，从而实现了光照分量的估计，并使用信息熵函数完成了对原图光照均匀程度的判定。同时，本文设计了一套基于所提算法的轴承滚子表面瑕疵检测系统，并对本文所提算法的性能进行验证。实验结果表明，本文算法能够有效滤除光照并且较好地保留图像细节，同时显著提升轴承滚子瑕疵检测系统的准确率。
system: [SYSTEM]

user: 2. 轴承滚子瑕疵检测系统 如图1所示，轴承滚子存在五个不同的表面，分别是大/小端面、大/小倒角及侧滚动面。本节将以轴承滚子大端面为例，介绍轴承滚子瑕疵检测系统。 图1. 轴承滚子各表面ROI图像 轴承滚子大端面瑕疵检测系统处理流程如图2所示，包括目标定位与分割、图像预处理、疑似缺陷定位和特征提取与缺陷确认4个主要的步骤。首先，对图像进行二值化处理，进行大端面区域的定位与分割，使用选择后的大端面区域与原图进行掩模操作得到待处理的大端面图像；随后，对定位分割后的大端面图像进行纹理过滤工作，降低纹理对检测的影响；然后，使用阈值分割对图像中出现灰度突变区域分割，定位得到疑似缺陷；最后，提取疑似缺陷特征并进行判断，得到确认缺陷。接下来，将对这4个主要步骤进行介绍。 图2. 基于机器视觉的轴承滚子检测系统处理流程 2.1. 目标定位与图像分割 本系统使用全局阈值分割方法对大端面区域的进行定位与分割。全局阈值分割分为固定阈值分割与动态阈值分割。固定阈值分割使用固定阈值对图像进行二值化操作，效果精确但是需根据灰度进行人为调整。动态阈值分割针对不同图像得到不同的阈值，效果一般但是参数不需要人为调整。大津法(OTSU) [ 12 ] 是具有代表性的动态阈值分割算法，其通过计算前景和背景的灰度方差的最大化，从而得到阈值。 如图3所示，OTSU算法虽然不需要人为调整参数来进行二值化，但是效果与固定阈值相比较差，速度相对较慢。同时，考虑到在实际的检测环境中，整体灰度变化不大，因此本文选择固定阈值进行大端面区域的定位与分割。 图3. OTSU分割与固定阈值分割对比 2.2. 图像预处理 系统中图像预处理步骤的目标是在保留或者不影响图像突变区域特征的前提下尽量滤除图像背景的纹理信息 [ 13 ] [ 14 ] 。因此，本文选取了各向异性滤波来完成图像的预处理工作，这种滤波方式能够在很好的保留图像边缘的同时，滤除图像纹理。 各向异性滤波通过各向异性扩散方程(P-M方程)完成图像滤波。P-M方程是图像去噪方法中的一个重要模型，其是传统热扩散方程的推广 [ 15 ] ，P-M方程定义如下： { ∂ u ∂ t = d i v ( g ( | ∇ u | 2 ， c ) ⋅ ∇ u ) u 0 = u ( x ， y ， 0 ) (1) 式中， g ( | ∇ u | 2 ， c ) 代表扩散系数，它根据图像的梯度来变化。 ∇ u 代表图像梯度， d i v 是散度算子， u 0 是原始图像， u ( x ， y ， 0 ) 是 t 时刻 u 0 的扩散图像。 2.3. 疑似缺陷定位 疑似缺陷是图像中的灰度突变区域。本文使用阈值分割方法来进行疑似缺陷的定位。上文中提到，鉴于固定阈值分割能够满足场景要求，因此采用固定阈值分割方法完成疑似缺陷定位。 在实际场景中，可将疑似缺陷分为疑似亮缺陷和疑似暗缺陷两类，具体定义如下： d ( x ， y ) = { flaw light ( x ， y ) ， f ( x ， y ) ≥ T l flaw dark ( x ， y ) ， f ( x ， y ) ≤ T d 0 ≤ T d < T l ≤ 255 (2) 式中， d ( x ， y ) 是已定位疑似缺陷中的像素， flaw light ( x ， y ) 是疑似亮缺陷中像素， flaw dark ( x ， y ) 是疑似暗缺陷中像素， f ( x ， y ) 是大端面中像素灰度值， T l ， T d 是进行疑似亮缺陷、疑似暗缺陷分割的阈值。本系统将大端面图像中灰度值高于 T l 的像素点定义为疑似亮缺陷 flaw light ( x ， y ) ，将大端面图像中灰度值低于 T d 的像素点定义为疑似暗缺陷 flaw dark ( x ， y ) 。 2.4. 缺陷确认 轴承滚子大端面中存在的主要缺陷如图4所示，主要有锈斑、划伤、欠磨等。锈斑具有灰度深且面积不定的特性；划伤具有灰度深、区域狭长、具有一定宽度的特性；欠磨具有其面积与非欠磨区域面积相差巨大的特性。针对各种缺陷特性，本文提取了灰度特征中的灰度均值特征、形状特征中的面积特征以及纹理特征中的对比度特征来描述疑似缺陷。 图4. 缺陷类型 面积特征通过统计疑似区域中像素数得出 [ 16 ] ；纹理特征中的对比度特征通过灰度共生矩阵GLCM (Gray-level co-occurrence matrix) [ 17 ] 计算得出，GLCM描述了存在于不同方向和距离的像素对的出现频率。根据GLCM可得对比度的定义如下： Contrast = ∑ i = 1 K ∑ j = 1 K ( i − j ) 2 p i j (3) 式中， K 是图像中的最大灰度级， i ， j 分别代表对应的像素的灰度， p i j 是在角度 θ 方向上相距 s 的两个像素灰度为 i ， j 的情况的出现频率。