CN106202506A - Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing - Google Patents

Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing Download PDF

Info

Publication number
CN106202506A
CN106202506A CN201610575917.3A CN201610575917A CN106202506A CN 106202506 A CN106202506 A CN 106202506A CN 201610575917 A CN201610575917 A CN 201610575917A CN 106202506 A CN106202506 A CN 106202506A
Authority
CN
China
Prior art keywords
noise
item
map
attentuation
dimensional traffic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610575917.3A
Other languages
Chinese (zh)
Inventor
蔡铭
张智伟
马侠霖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN201610575917.3A priority Critical patent/CN106202506A/en
Publication of CN106202506A publication Critical patent/CN106202506A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46Multiprogramming arrangements
    • G06F9/50Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5083Techniques for rebalancing the load in a distributed system
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/28Indexing scheme for image data processing or generation, in general involving image processing hardware

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • Databases & Information Systems (AREA)
  • Architecture (AREA)
  • Computer Hardware Design (AREA)
  • Data Mining & Analysis (AREA)
  • Traffic Control Systems (AREA)

Abstract

The present invention proposes the three-dimensional traffic Noise map update method of a kind of combination offline storage and instant computing, and the calculating of noise is divided into the noise emission item relevant to traffic flow and noise attentuation item two parts unrelated with traffic flow;Only relevant to geographic element noise attentuation item is calculated and by its offline storage in advance in data base, when three-dimensional traffic Noise map needs to be updated, change according to urban road traffic flow modes, only calculate the noise emission item of change, noise attentuation item in conjunction with offline storage is calculated final noise figure, reduces amount of calculation.For accelerating the renewal speed of three-dimensional traffic Noise map, this method uses the mode of computer cluster parallel computation to be accelerated noise calculation and the renewal of Noise map.The features such as it is good that the distributed computing system that this method is used in combination has load equilibrium, and motility is high and autgmentability is strong, coordinate the Noise map update algorithm after optimizing can greatly reduce the calculating time, are obviously improved renewal efficiency.

Description

Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing
Technical field
The present invention relates to road traffic three-dimensional noise and calculate renewal and computer cluster distributed computing technology field, more In particular it relates to a kind of combination offline storage and the three-dimensional traffic Noise map Rapid Updating of instant computing.
Technical background
Three-dimensional traffic Noise map is possible not only to show the distribution situation of traffic noise on ground, simultaneously can also be to towering The noise profile situation of numerous building surfaces on the ground is shown.In view of the noise source affecting people's daily life Mainly road traffic noise, and the state of road traffic is to have regular hour characteristic.Therefore, according to road traffic shape Noise map is quickly updated by state, and applies three-dimensional traffic Noise map to open up the noise profile situation of survey region Show, analyze and assess, the spatial and temporal distributions of the road traffic noise of research specific region is had important effect.
Research and application that current three-dimensional traffic Noise map quickly updates are substantially by using supercomputer parallel Calculate or Distributed Calculation is to reduce the time that big noise region calculates.Such as Noise map business software SoundPlan, Cadna/ A, Lima etc. provide the function of Distributed Calculation: multinode disposes software for calculation, subtask distribution, each node subtask meter Calculate, subtask result collects and shows.But these business softwares also exist, and deployment is complicated, system flexibility is poor, calculate task management The shortcoming such as dumb, cannot play the computing capability of each computer node during Distributed Calculation completely.
As can be seen here, the renewal of the three-dimensional traffic Noise map in the biggest region still suffers from updating computationally intensive, consumption Take the problems such as calculating resource is many.
Summary of the invention
It is an object of the invention to the existing issue overcoming the three-dimensional traffic Noise map in big region quickly to update, it is provided that one Planting and have ease-to-operate, the combination offline storage of the extensibility of load equilibrium and height is handed over the three-dimensional of instant computing Logical Noise map Rapid Updating.
In order to realize above-mentioned goal of the invention, the technical solution used in the present invention is as follows:
A kind of combination offline storage and the three-dimensional traffic Noise map Rapid Updating of instant computing, including following step Rapid:
A. the calculating of noise is divided into calculating noise emission item and two parts of noise attentuation item.The calculating of noise emission item Only relevant to the traffic flow modes of urban road, noise attentuation item geography only and between traffic noise sound source and noise Rx point Factor is correlated with;
Further, the bicycle noise emission of large, medium and small three kinds of different automobile types is considered among noise emission item;Noise declines Deduction item includes under the noise range attenuation under three dimensions, the buildings noise attentuation under three dimensions and three dimensions Sound barrier noise attentuation;
B. computer cluster Distributed Calculation noise is utilized to decline item.The distributed computing system that method uses is based on tile ground Noise calculation region is cut into tiling area, using tiling area as computing unit, and distributed by figure service (TMS) thought Use the method calculating node active request task to replace traditional method passively received an assignment during calculating, make task divide Join and there is harmony, and by noise attentuation item offline storage to data base.
Further, method uses region segmentation mechanism to split basic data with tiling area for unit, and point Do not store.18th level of the size correspondence tile map service standard of every map tile, and by each noise Rx point Noise Rx radius is set to 500m.
Further, during Distributed Calculation, each calculating node uses Data cache technology, reduce redundancy Data transmission in a network.
C. calculate noise emission item, and combine noise attentuation item and be calculated the noise figure of noise Rx point.Because of traffic The change of stream and when carrying out the renewal of Noise map, noise emission item the most again, and combine offline storage by energy supposition method Noise attentuation item be calculated new noise profile result, update three-dimensional traffic Noise map with this.
Compared with prior art, the invention have the advantages that
The present invention proposes the three-dimensional traffic Noise map Rapid Updating of a kind of combination offline storage and instant computing. The method is by combining instant computing with offline storage noise attentuation item, it is possible to be used in combination distributed computing system etc. The three-dimensional traffic Noise map in big region is quickly updated by means.The Noise map update algorithm that method is used is counted in advance Calculate each noise attentuation item received between point and noise source, and by its offline storage in data base.Additionally, the method is also The calculating speed of noise can be further speeded up in conjunction with distributed computing system.This distributed computing system uses tile map clothes Noise calculation region is cut into tiling area, using tiling area as computing unit, and in Distributed Calculation by business (TMS) thought During use the method calculating node active request task to replace the traditional method passively received an assignment, have easy and simple to handle Property, load equilibrium and the extensibility of height.This Rapid Updating can be applied not only to have mainframe computer collection In the middle of the laboratory of group or supercomputer, it is also possible to be applied to common individual's electricity of the most several phase the Internets On brain, by the method using instant computing and offline storage to combine, while significantly speeding up noise renewal, have very Strong adaptability and autgmentability.
The present invention is the three-dimensional traffic Noise map Rapid Updating of a kind of combination offline storage and instant computing, solves The renewal that the urban road three-dimensional traffic Noise map in big region faces when updating is computationally intensive, expends that to calculate resource many etc. Problem.
Accompanying drawing explanation
Fig. 1 is method overall structure schematic diagram in the present invention.
Fig. 2 is world map map tile numbering result schematic diagram when level is 2 in the present invention.
Fig. 3 is that in the present invention, regional road traffic noise profile formula calculates design flow diagram.
Detailed description of the invention
Accompanying drawing being merely cited for property explanation, it is impossible to be interpreted as the restriction to this patent;In order to the present embodiment is more preferably described, attached Scheme some parts to have omission, zoom in or out, do not represent the size of actual product;
To those skilled in the art, in accompanying drawing, some known features and explanation thereof may be omitted is to be appreciated that 's.
As it is shown in figure 1, the three-dimensional traffic Noise map of the combination offline storage of present invention proposition and instant computing is the most more New method, mainly contains 3 modules: data acquisition process module, noise profile formula computing module and Noise map update Module.Whole Noise map calculates and more new technological process is as follows:
Step 1: data acquisition process module, by GIS platform, uses Arcgis Engine API coding, from GIS Data base reads urban geography element information.
Step 2: utilize the gps data of Floating Car in city, in conjunction with by Green's primary speed-density logarithmic relationship model The Speed-flow Relationship model derived is to obtain the vehicle flowrate data of each road section in city.
Step 3: use region segmentation mechanism for unit, basic data to be split with tiling area, tile is compiled Number, and set up latitude and longitude coordinates to the corresponding relation between tile number.Such as, tile number coordinate origin is scheduled on the world The upper left corner of map, numbering X-axis overlap with north latitude 85.05 ° and direction to the left, Y-axis and east longitude 180 ° (also be 180 ° of west longitude) are heavily Close and direction is downward.In conjunction with Mercator projection formula, the latitude and longitude coordinates of its correspondenceTo tile number (Xtile, Ytile) conversion formula is as follows:
Wherein, n=2level, level is tile level;λ andIt is respectively longitude and latitude (unit radian).Work as level When=2, numbering X-direction and Y-direction are respectively arranged with tile quantity n=4, and total tile data is 16, and its numbering result is as shown in Figure 2.Will Various geographic elements store in corresponding tile, and are stored in data storage server by all of tile data.
Step 4: use Distributed Computing Platform to carry out the attenuation term of noise calculating and offline storage: to use multiple calculating Node, and on each calculating node, run distributed-computation program client so that it is participate in the middle of Distributed Calculation, such as figure Shown in 3.Wherein, the process of the participation Distributed Calculation of each calculating node is as follows:
(1) before calculating starts, calculate and control server according to zoning scope and level formation zone, tiling area The tile number set inside comprised, is formed and calculates task list, the corresponding tile number of each task;
(2) when calculating node A participation and entering calculating, its active calculates task from calculating control server request;
(3) calculate control server from be not fully complete task list choose a tile number T and distribute to calculate node A;
(4) calculate the tile number T that node A arrives according to request, from data server requests, the noise in this tile is connect All tile datas (including this tile itself) in the influential regional extent of sink;
(5) after receiving all relevant basic datas, the node A setting according to reception point generating mode is calculated, raw Become the noise Rx point in region corresponding to tile number T, and calculate each noise figure receiving point;
(6) after completing to calculate, noise calculation result being returned to data server, data server is equally with tiling area The noise result of this tile is stored for unit;
(7) after confirming that result of calculation returns successfully, send message informing and calculate control server by tile number T never Complete task list is rejected;
Now, calculate node A be in can duty, jump to second step, control the new tile of server request from calculating Number and circulate said process, until calculating the task of controlling not to be not fully complete in server.
Step 5: calculate noise emission item, and combine noise attentuation item and be calculated the noise figure of noise Rx point.Below The flow process quickly updated for Noise map:
(1) when regional road traffic stream changes and needs to update noise result, every road section is calculated it and hands over Through-flow item;
(2) for a certain noise Rx point, its section tag set stored is utilized to index road influential on it Final noise figure, to obtain corresponding traffic flow item, in conjunction with the fixterm of offline storage, is tried to achieve in section, updates three-dimensional friendship with this Logical Noise map.
Obviously, the above embodiment of the present invention is only for clearly demonstrating example of the present invention, and is not right The restriction of embodiments of the present invention.For those of ordinary skill in the field, the most also may be used To make other changes in different forms.Here without also cannot all of embodiment be given exhaustive.All at this Any amendment, equivalent and the improvement etc. made within the spirit of invention and principle, should be included in the claims in the present invention Protection domain within.

Claims (8)

1. the three-dimensional traffic Noise map update method combining offline storage and instant computing, it is characterised in that described side Method comprises the following steps:
A, the calculating of noise is divided into calculating noise emission item and two parts of noise attentuation item;
B, computer cluster distributed speed-up computation noise is utilized to decline item, and by noise attentuation item offline storage to data base;
C, calculate noise emission item, and the noise attentuation item combining offline storage is calculated the noise figure of noise Rx point, with This updates three-dimensional traffic Noise map.
Three-dimensional traffic Noise map update method the most according to claim 1, it is characterised in that: in described step a, noise Discharge item calculating only relevant to the traffic flow modes of urban road, noise attentuation item only with traffic noise sound source and noise Rx Geographic factor between point is correlated with, and the calculating of noise emission item and noise attentuation item is non-interference.
Three-dimensional traffic Noise map update method the most according to claim 2, it is characterised in that: examine among noise emission item Consider the bicycle noise emission of large, medium and small three kinds of different automobile types.
Three-dimensional traffic Noise map update method the most according to claim 2, it is characterised in that: noise attentuation item includes three Sound barrier noise under noise range attenuation under dimension space, the buildings noise attentuation under three dimensions and three dimensions Decay.
Three-dimensional traffic Noise map update method the most according to claim 1, it is characterised in that: in described step b, utilize Computer cluster distributed speed-up computation noise declines item, is that the distributed computing system used is based on tile map service (TMS) Noise calculation region is cut into tiling area by thought, using tiling area as computing unit, and during Distributed Calculation Use the method calculating node active request task, make task distribution have harmony.
Three-dimensional traffic Noise map update method the most according to claim 5, it is characterised in that: use region segmentation mechanism For unit, basic data is split with tiling area, and store respectively, the size correspondence tile map of every map tile 18th level of service standard, and the noise Rx radius of each noise Rx point is set to 500m.
Three-dimensional traffic Noise map update method the most according to claim 5, it is characterised in that: in Distributed Calculation process In, each calculating node uses Data cache technology, reduces redundant data transmission in a network.
Three-dimensional traffic Noise map update method the most according to claim 1, it is characterised in that: in described step c, because of The change of traffic flow and when carrying out the renewal of Noise map, noise emission item the most again, and combine off-line by energy supposition method The noise attentuation item of storage is calculated new noise profile result, updates three-dimensional traffic Noise map with this.
CN201610575917.3A 2016-07-20 2016-07-20 Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing Pending CN106202506A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610575917.3A CN106202506A (en) 2016-07-20 2016-07-20 Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610575917.3A CN106202506A (en) 2016-07-20 2016-07-20 Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing

Publications (1)

Publication Number Publication Date
CN106202506A true CN106202506A (en) 2016-12-07

Family

ID=57491833

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610575917.3A Pending CN106202506A (en) 2016-07-20 2016-07-20 Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing

Country Status (1)

Country Link
CN (1) CN106202506A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106611438A (en) * 2016-12-27 2017-05-03 广州都市圈网络科技有限公司 Method and device for updating and slicing local area of three-dimensional simulation map
CN106960097A (en) * 2017-03-27 2017-07-18 北京图安世纪科技股份有限公司 One kind builds three-dimensional city level road network and microscopic traffic flow simulation method and system
CN107564098A (en) * 2017-08-17 2018-01-09 中山大学 A kind of Fast rendering method of big Local Area Network three-dimensional noise map
CN108334726A (en) * 2017-08-17 2018-07-27 中山大学 A kind of urban area road noise distribution forecasting method
CN108801449A (en) * 2018-07-09 2018-11-13 滨州职业学院 A kind of urban traffic noise inquiry system of wireless monitor
CN109063288A (en) * 2018-07-19 2018-12-21 中山大学 A kind of city area-traffic noise profile parallel calculating method
CN111492323A (en) * 2017-12-20 2020-08-04 Wing航空有限责任公司 Mitigating noise exposure of unmanned aerial vehicles
CN112287062A (en) * 2020-12-23 2021-01-29 智道网联科技(北京)有限公司 Real-time traffic information synchronization method and device and electronic equipment
CN112463898A (en) * 2020-10-21 2021-03-09 中山大学 Noise map updating method combining speed and noise monitoring data

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070031237A1 (en) * 2005-07-29 2007-02-08 General Electric Company Method and apparatus for producing wind energy with reduced wind turbine noise
CN103440411A (en) * 2013-08-15 2013-12-11 中山大学 Traffic noise pollution model based on exposed crowd/area/acoustic environment functional area

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070031237A1 (en) * 2005-07-29 2007-02-08 General Electric Company Method and apparatus for producing wind energy with reduced wind turbine noise
CN103440411A (en) * 2013-08-15 2013-12-11 中山大学 Traffic noise pollution model based on exposed crowd/area/acoustic environment functional area

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CNTCHEN: "国内主要地图瓦片坐标系定义及计算原理", 《HTTP://CNTCHEN.GITHUB.IO/2016/05/09/%E5%9B%BD%E5%86%85%E4%B8%BB%E8%A6%81%E5%9C%B0%E5%9B%BE%E7%93%A6%E7%89%87%E5%9D%90%E6%A0%87%E7%B3%BB%E5%AE%9A%E4%B9%89%E5%8F%8A%E8%AE%A1%E7%AE%97%E5%8E%9F%E7%90%86/》 *
CSDN: "OpenStreetMap/Google/百度/Bing瓦片地图服务(TMS)", 《HTTPS://BLOG.CSDN.NET/YOUNGKINGYJ/ARTICLE/DETAILS/23365849》 *
李楠等: "基于面向服务对象体系结构的交通噪声地图分布式计算方法", 《计算机应用》 *
蔡铭等: "广州市道路交通噪声地图的绘制", 《环境工程》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106611438B (en) * 2016-12-27 2019-12-24 广州都市圈网络科技有限公司 Local area updating and map cutting method and device of three-dimensional simulation map
CN106611438A (en) * 2016-12-27 2017-05-03 广州都市圈网络科技有限公司 Method and device for updating and slicing local area of three-dimensional simulation map
CN106960097A (en) * 2017-03-27 2017-07-18 北京图安世纪科技股份有限公司 One kind builds three-dimensional city level road network and microscopic traffic flow simulation method and system
CN106960097B (en) * 2017-03-27 2020-06-16 北京图安世纪科技股份有限公司 Simulation method and system for constructing three-dimensional urban road network and microcosmic traffic flow
CN108334726A (en) * 2017-08-17 2018-07-27 中山大学 A kind of urban area road noise distribution forecasting method
CN107564098A (en) * 2017-08-17 2018-01-09 中山大学 A kind of Fast rendering method of big Local Area Network three-dimensional noise map
CN111492323A (en) * 2017-12-20 2020-08-04 Wing航空有限责任公司 Mitigating noise exposure of unmanned aerial vehicles
CN111492323B (en) * 2017-12-20 2023-08-01 Wing航空有限责任公司 Mitigating noise exposure of unmanned aerial vehicles
CN108801449A (en) * 2018-07-09 2018-11-13 滨州职业学院 A kind of urban traffic noise inquiry system of wireless monitor
CN109063288A (en) * 2018-07-19 2018-12-21 中山大学 A kind of city area-traffic noise profile parallel calculating method
CN112463898A (en) * 2020-10-21 2021-03-09 中山大学 Noise map updating method combining speed and noise monitoring data
CN112287062A (en) * 2020-12-23 2021-01-29 智道网联科技(北京)有限公司 Real-time traffic information synchronization method and device and electronic equipment
CN112287062B (en) * 2020-12-23 2021-04-02 智道网联科技(北京)有限公司 Real-time traffic information synchronization method and device and electronic equipment

Similar Documents

Publication Publication Date Title
CN106202506A (en) Three-dimensional traffic Noise map update method in conjunction with offline storage Yu instant computing
Yang et al. A spatiotemporal and graph-based analysis of dockless bike sharing patterns to understand urban flows over the last mile
Yin et al. Group-wise itinerary planning in temporary mobile social network
CN103871234A (en) Grid mapping growth-based traffic network division method and configuration server
CN105160707B (en) Threedimensional model quick visualization method based on viewpoint index
Kim Identifying the structure of cities by clustering using a new similarity measure based on smart card data
CN103337161A (en) Optimization method of intersection dynamic comprehensive evaluation and signal control system based on real-time simulation model
WO2013020611A1 (en) Generating an od matrix
CN104750895A (en) Real-time city emergency evacuation simulating method and system based on mobile phone data
Chen et al. Discovering urban traffic congestion propagation patterns with taxi trajectory data
Higuchi et al. How to keep a vehicular micro cloud intact
CN103968854A (en) Route cluster planning method
CN105427210A (en) Tourism resource information collection and evaluation network system
CN101719148A (en) Three-dimensional spatial information saving method, device, system and dispatching system
CN104573140A (en) Layered dynamic path planning method applied to virtual simulation
Wang et al. User recruitment system for efficient photo collection in mobile crowdsensing
CN105574194A (en) Coordinate point processing method and device for electronic map interface
CN109274645A (en) A kind of hierarchical layered access implementation method of smart city space-time cloud platform
CN106130110B (en) The electric taxi charging station constant volume method on trip ground is selected based on stratified probability
Boutsis et al. Efficient event detection by exploiting crowds
Chen et al. QarSUMO: a parallel, congestion-optimized traffic simulator
Jung et al. Simulation framework for modeling large-scale flexible transit systems
Xie et al. An energy-efficient high definition map data distribution mechanism for autonomous driving
Tao et al. Vehicular data cloud platform with 5G support: Architecture, services, and challenges
CN107588779B (en) Intelligent vehicle navigation method based on travel time between any two nodes

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161207

RJ01 Rejection of invention patent application after publication