CN110472179A - Round bullet impact and the case suspect of the perpendicular window in ground check method - Google Patents

Round bullet impact and the case suspect of the perpendicular window in ground check method Download PDF

Info

Publication number
CN110472179A
CN110472179A CN201910572332.XA CN201910572332A CN110472179A CN 110472179 A CN110472179 A CN 110472179A CN 201910572332 A CN201910572332 A CN 201910572332A CN 110472179 A CN110472179 A CN 110472179A
Authority
CN
China
Prior art keywords
image
glass
trajectory
angle
search
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.)
Granted
Application number
CN201910572332.XA
Other languages
Chinese (zh)
Other versions
CN110472179B (en
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.)
Nanjing Huya Mdt Infotech Ltd
Gulou Branch Of Nanjing Public Security Bureau
Original Assignee
Nanjing Huya Mdt Infotech Ltd
Gulou Branch Of Nanjing Public Security Bureau
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 Nanjing Huya Mdt Infotech Ltd, Gulou Branch Of Nanjing Public Security Bureau filed Critical Nanjing Huya Mdt Infotech Ltd
Priority to CN201910572332.XA priority Critical patent/CN110472179B/en
Publication of CN110472179A publication Critical patent/CN110472179A/en
Application granted granted Critical
Publication of CN110472179B publication Critical patent/CN110472179B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Optimization (AREA)
  • Remote Sensing (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Biology (AREA)
  • Computational Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a kind of round bullet impacts and the case suspect of the perpendicular window in ground to check method, and ballistic trajectory parameter calculates, and obtains the numerical value of the angle γ of instantaneous velocity V and glass planar;Measure the numerical value of shell hole image long axis and glass planar y' axis angle α;Calculate the angle beta of trajectory plane and glass plane normal;Trajectory parabolic path is calculated, calculates ballistic trajectory towards module;With establishing emphasis personnel's household register and residence latitude and longitude coordinates database;Formulate three-level search range;The suspect of different regions of search is checked using means of different.Operation of the present invention is easy, it is easy to accomplish, maximum precise restoration ballistic trajectory parabola, whole-process automatic search investigation uses manpower and material resources sparingly;The positional relationship that can show suspect's Location Yu specific case projectile's path track makes search investigation have more specific aim;According to the suspicion grade of search grade classification suspect, keep search investigation more scientific.

Description

Round bullet impact and the case suspect of the perpendicular window in ground check method
Technical field
The present invention relates to a kind of round bullet impacts and the case suspect of the perpendicular window in ground to check method, belongs to Criminal investigation technical field.
Background technique
In recent years, with the extensive use of social economy's fast development and internet, illegal activities people can be increasingly lighter It changes places and obtains gun-simulation, air gun and this kind of contraband of bow and crossbow, leading to steel ball class illegal activities in recent years is in continue high-incidence situation.
In such case, often only a glass hole can be left in aggrieved glass in other.To such case of investigating and prosecuting Part needs to solve the problems, such as three: 1, reduction transmitting steel ball trajectory parabola judges to emit position;2, to periphery emphasis suspicion crowd Seeking investigation is carried out, the doubtful personnel of emphasis are found;3, research is compared to periphery history case information, finds out case in region Feature occurs.
Tradition is handled a case, and there are following disadvantages for means: 1, trajectory parabola and county of emission judgement lack effectively quantitative calculate Method relies primarily on personal experience;2, the sequence of suspicion degree can not be carried out to multiple suspects according to on-site actual situations;3, for Periphery history case information is only capable of making the record on when and where, can not further be compared and be divided according to actual transmission trajectory Analysis;4, three above problem needs to solve respectively by the not homologous ray of different departments, mutually isolates and lacks to each other effectively Communication and coordination mechanism.
How to pass through analysis steel ball and beat the trace after glass, to determine the ballistic trajectory of transmitting steel ball, to be further It determines lawbreaker's location-based service, is those skilled in the art's problem in urgent need to solve.
Calculate that the parabolical calculation method of trajectory is primarily present following difficult point by shell hole trace:
1, it can not be measured for the vector component velocity of trajectory parabola calculation and the angular relationship of projectile flight speed;
2, the orientation angles of the trajectory plane of projectile flight can not measure;
3, the different angle relationship between glass planar and ground can not measure the influence that ballistic trajectory calculates.
It is glass planar and the perpendicular window in ground that the bullet to be studied of the present invention, which beats object, as shown in Figs. 1-2, including extremely Few one of which: sliding window and casement window.
Summary of the invention
Purpose: in order to overcome the deficiencies in the prior art, the present invention provides a kind of round bullet impact with ground and is in The case suspect of vertical window checks method, and such bullet shape is fixed, and flight is by air in calm or gentle breeze environment Drag effects are smaller, and for velocity attenuation less than 10%, trajectory is similar to parabolic motion, it then follows parabolic motion rule, it is public Formula isV in formulaLAnd VHThe respectively component velocity of bullet in the horizontal and vertical directions, the two with There are geometry teaching relationships by speed V in projectile flight, and are influenced with the angular values of V by following situations: 1, projectile impact The instantaneous velocity and glass planar angle of glass;2, inclined degree of the shell hole trace in glass planar;3, glass planar and ground The tilt angle in face.
Technical solution: in order to solve the above technical problems, the technical solution adopted by the present invention are as follows:
A kind of case suspect investigation method of the impact of circle bullet and the perpendicular window in ground, includes the following steps:
Ballistic trajectory parameter calculates:
Obtain the numerical value of the angle γ of instantaneous velocity V and glass planar;
Measure the numerical value of shell hole image long axis and glass planar y' axis angle α;
Calculate the angle beta of trajectory plane and glass plane normal;
Calculate trajectory parabolic path:
The instantaneous stroke speed V of bullet is set to limit variable, then trajectory parabola indicates are as follows:
T in formula is time variable, and g is acceleration of gravity, shot orbit horizontal axis flying distance L, the shot orbit longitudinal axis Flying height H, h are shell hole o point apart from ground level;
Set the horizontal maximum distance L of trajectory parabolamaxTo limit variable, then trajectory parabola indicates are as follows:
T in formula is time variable, and g is acceleration of gravity, shot orbit horizontal axis flying distance L, the shot orbit longitudinal axis Flying height H, h are shell hole o point apart from ground level;
Set ballistic trajectory apogee distance shell hole o point maximum normal distance h1To limit variable, then trajectory parabola table It is shown as:
T in formula is time variable, and g is acceleration of gravity, shot orbit horizontal axis flying distance L, the shot orbit longitudinal axis Flying height H, h are shell hole o point apart from ground level;
Calculate ballistic trajectory direction:
Plane is on the right side of plane where glass planar normal where trajectory parabola, then the parabolical azimuth of trajectory For δ+β, plane is in the left side of plane where glass planar normal, the then parabolical azimuth of trajectory where trajectory parabola δ-β, δ are the azimuth of the normal of glass planar where shell hole;
With establishing emphasis personnel's household register and residence latitude and longitude coordinates database;
Formulate three-level search range;
By the latitude and longitude coordinates mark of the household register address of illegal previous conviction personnel or residence location or household register address and residence location It infuses in the region of search on map, wherein the suspect of different regions of search is checked using means of different.
Preferably, described
Preferably, steps are as follows for the numerical value acquisition of the angle γ of instantaneous velocity V and glass planar:
Step 1: hole on glass image recognition: passing through image algorithm and the convolutional neural networks bullet based on TensorFlow The angle γ of instantaneous velocity and glass when hitting glass is identified;
1.1 sample image pre-processings:
Angle γ is set as 0-90 degree, average every image increases by 1 degree, and totally 91 hole on glass images are input to machine In, binary conversion treatment is carried out, is removed variegated;
Using the outermost pixel on shell hole image four sides up and down as basic point, expands a certain distance respectively and intercept, obtain To master sample image;
1.2 sample image feature extractions:
Master sample image is input in the convolutional neural networks of TensorFlow and carries out the operation of n times convolution pondization, it will Sample image gradually decomposes, and calculates the characteristic value of sample image;
1.3 sample images substitute into training:
The hole on glass image pattern of multiple same angle γ angles is input in machine, 1.1 and 1.2 behaviour is repeated Make, machine is allowed to carry out a large amount of repetition training;
Step 2: the identification of live image:
2.1 live image pre-processings:
The hole on glass image of live shooting is input in machine, image binaryzation processing and image interception is carried out, obtains Standard to live shell hole image identifies format;
2.2 live image feature extractions:
The live shell hole image of reference format is input in the convolutional neural networks of TensorFlow, n times convolution is carried out Chi Hua calculates the characteristic value of live shell hole image;
2.3 aspect ratios pair obtain angle γ:
The feature of different angle sample image in the characteristic value and sample database of the live shell hole image that system-computed is gone out Value is compared, and carries out result output for the highest angle of degree of fitting as the angle of live shell hole image.
Preferably, described includes: to public affairs with residence latitude and longitude coordinates database with establishing emphasis personnel's household register Pacify internal database carry out information retrieval, find out because gun-related case handled or be listed in accordance with the law major tasks object it is illegal before Section personnel, household register address and residence location to illegal previous conviction personnel carry out inverse geocoding, obtain its household register address and now live The latitude and longitude coordinates information of address, and database table is established, table field includes: name, gender, citizenship number, case Part number, the detailed location of household register, household register X-coordinate, household register Y-coordinate, the detailed location in personnel address, address X-coordinate, is now lived at brief merit Location Y-coordinate.
Preferably, the formulation three-level search range includes:
2-1 determines glass planar X-axis both ends azimuth:
2-2 determines three-level search mechanisms:
The latitude and longitude coordinates of shell hole are labeled on map by 2-2-1, are set as O;
For 2-2-2 using shell hole O as the center of circle, certain distance is radius in map upper drawing circle;
2-2-3 draws straight line using O as starting point, to glass planar normal towards the azimuth direction δ, meets at point F with circle;It is with O Starting point draws straight line to the both direction vertical with azimuth δ respectively, meets at point A, B respectively with circle;Using O as starting point, thrown to trajectory Object line draws straight line towards the azimuth direction δ+β or δ-β, meets at point D with circle;Radius OD is rotated into certain angle to two sides respectively Degree meets at point C, E with circle respectively;
Region of search is arranged in 2-3:
The fan-shaped region 2-3-1 COE is set as Primary search region;
The fan-shaped region 2-3-2 BOC, fan-shaped region AOE are set as secondary search region;
The fan-shaped region 2-3-3 BOA is set as three-level region of search.
Preferably, glass planar X-axis both ends azimuth preparation method is as follows: two sides of glass planar X-axis To azimuth be+90 °, δ -90 ° of δ.
Preferably, the certain angle is set as a value in section (15 ° -20 °) range.
Preferably, the certain distance is set as a value in section (300M-400M) range.
Preferably, the suspect of the different regions of search includes: using means of different investigation
Suspect in 3-1 Primary search region, being defined as highest-order emission suspicion is that emphasis checks object;
Suspect in 3-2 secondary search region is defined as intermediate transmitting suspicion;
Suspect in 3-3 three-level region of search, is defined as lower-order emission suspicion.
The utility model has the advantages that the case suspect investigation side of circle bullet impact and the perpendicular window in ground provided by the invention Method is α by estimation, shell hole image long axis and the angle of glass planar y' axis to instantaneous velocity and the angle γ of glass planar Measurement, calculate the angle beta of trajectory plane Yu glass plane normal, and combine live actual environment, set following three item numbers One in:
1, speed when projectile impact glass;
2, the parabolical maximum height in place;
3, the parabolical longest distance in place.
Parabola is calculated, to calculate the direction δ ± β and distance in transmitting bullet place.
Latitude and longitude coordinates information in the household register address and residence location for grasping the illegal emphasis personnel of Related Cases previous conviction, with And on the multi-pass decoding manufacturing basis made towards azimuth using projectile's path track, potential suspicion personnel are flown with bullet Row trajectory makees the association analysis on geographical location.
There are following advantages for this method:
1, related fields blank has been filled up, a kind of calculating side for extrapolating entire ballistic trajectory from shell hole trace has been invented Method;
2, this method only need to make DATA REASONING for the shell hole after projectile impact, without measuring entire ballistic flight rail Mark, it is easy to operate, it is easy to accomplish;
3, this method can be according to live different situations setup parameter, and all algorithms meet physics by stringent operation And mathematical operation rule, it can maximum precise restoration ballistic trajectory parabola;
4, whole-process automatic search investigation, uses manpower and material resources sparingly;
5, the positional relationship that can show suspect's Location Yu specific case projectile's path track makes search investigation more Has specific aim;
6, according to the suspicion grade of search grade classification suspect, keep search investigation more scientific.
Detailed description of the invention
Fig. 1 is the structural schematic diagram of sliding window;
Fig. 2 is the structural schematic diagram of casement window;
Fig. 3 is sample image pre-processing schematic diagram;
Fig. 4 is master sample image procossing schematic diagram;
Fig. 5 is sample image feature extraction processing schematic;
Fig. 6 is the schematic diagram of the angle γ of instantaneous velocity and glass planar;
Fig. 7 is the schematic diagram that the angle of shell hole image long axis and glass planar y' axis is α;
Fig. 8 is the schematic diagram of ballistic trajectory plane and glass plane normal y'o'z' angle β;
Fig. 9 is the schematic diagram that instantaneous velocity V is η with the y-axis angle in three-dimensional system of coordinate;
Figure 10 is the triangle relation schematic diagram of angle α, γ, β, η;
Figure 11 is the schematic diagram of shot orbit;
Figure 12 is the relation schematic diagram of ballistic trajectory direction;
Figure 13 is region of search schematic diagram;
Figure 14 is three-level region of search schematic diagram.
Specific embodiment
The present invention will be further explained with reference to the accompanying drawing.
A kind of case suspect investigation method of the impact of circle bullet and the perpendicular window in ground, includes the following steps:
Step 1: hole on glass image recognition: passing through image algorithm and the convolutional neural networks bullet based on TensorFlow The angle γ of instantaneous velocity and glass when hitting glass is identified.
1.1 sample image pre-processings:
As shown in figure 3, setting angle γ as 0-90 degree, average every image increases by 1 degree, and totally 91 hole on glass images are defeated Enter into machine, carry out binary conversion treatment, removes variegated.
As shown in figure 4, expanding a certain distance respectively using the outermost pixel on shell hole image four sides up and down as basic point And intercept, obtain master sample image.
1.2 sample image feature extractions:
As shown in figure 5, master sample image is input to progress n times convolution pond in the convolutional neural networks of TensorFlow Change operation, sample image is gradually decomposed, calculates the characteristic value of sample image by such as n=5.
1.3 sample images substitute into training:
The hole on glass image pattern of multiple same angle γ angles is input in machine, 1.1 and 1.2 behaviour is repeated Make, allows machine to carry out a large amount of repetition training, keep characteristic value more accurate.
Step 2: the identification of live image:
2.1 live image pre-processings:
The hole on glass image of live shooting is input in machine, image binaryzation processing and image interception is carried out, obtains Standard to live shell hole image identifies format.
2.2 live image feature extractions:
The live shell hole image of reference format is input in the convolutional neural networks of TensorFlow, n times convolution is carried out Chi Hua calculates the characteristic value of live shell hole image.
2.3 aspect ratios pair obtain angle γ:
The feature of different angle sample image in the characteristic value and sample database of the live shell hole image that system-computed is gone out Value is compared, and carries out result output for the highest angle of degree of fitting as the angle of live shell hole image.
Step 3: ballistic trajectory parameter calculates:
3.1 obtain the angle γ of instantaneous velocity and glass planar, as shown in fig. 6, o ' is shell hole, V is that steel ball hits glass Instantaneous velocity when plane, x ' are the trunnion axis of glass planar, and y ' is the vertical axis of glass planar, and o ' a is V in glass planar Projection, angle of the angle γ between V and o ' a.
When 3.2 bullet tilt angles hit glass, shell hole image is an ellipse, and oval long axis is that shell hole image is long Axis, if the angle of shell hole image long axis and glass planar y' axis is α, as shown in fig. 7, wherein shell hole image long axis direction is wink Projecting direction of the Shi Sudu V in glass planar.
3.3 calculate the angle beta of trajectory plane and glass plane normal.
Instantaneous velocity is V when 3.3.1 setting steel ball shock glass;
3.3.2 setting when hoodle hits glass has glass planar x' axle speed Vx', glass planar y' axle speed Vy' and perpendicular to The z' axle speed Vz' of glass planar;
3.3.3 as shown in figure 8, setting ballistic trajectory plane and plane x'o'z' intersecting straight lines o'b, then line o'b and line o'z' Angle is ballistic trajectory plane and glass plane normal y'o'z' angle β.
3.3.4 obtaining following relationship:
Formula 3 will be substituted by formula 1,2 above.
3.4 calculate y-axis angle η in trajectory parabola and three-dimensional system of coordinate.
Since glass planar is perpendicular to ground, x' axis, y' axis and perpendicular to the z' axis of glass planar and three in glass planar X-axis, y-axis and z-axis in dimension coordinate system coincide, so:
As shown in figure 9, set trajectory parabola most end point, in the instantaneous velocity V and three-dimensional system of coordinate when projectile impact glass Y-axis angle be η, ob is projection of the instantaneous velocity V on plane xoz, and plane yob is plane where trajectory parabola;
As shown in Figure 10, following relational expression can be obtained:
∴vy=vsin γ tan β cot α
∵ vcos η=vy
∴ cos η=sin γ tan β cot α
Wherein, willSubstitute into above formula.
∴ cos η=cos γ cos α
∴ η=arccos (cos γ cos α)
It can be seen that it is unrelated with the numerical value of instantaneous velocity V to calculate angle η.
Step 4: calculating trajectory parabolic path.
It is shot orbit flying distance L that trajectory parabolic path plane, which is set as horizontal axis, and the longitudinal axis is that shot orbit flight is high The two-dimensional plane coordinate system of H is spent, as shown in figure 11.
It is hole on glass position that o point, which is projectile impact point, in figure, and c point is coordinate origin, and the angle η is projectile impact The angle of instantaneous velocity V and y-axis both H axis when glass, a point are ballistic trajectory apart from ground highest point, and b point is the minimum transmitting of bullet Both ground launch position, h are shell hole o point apart from ground level, h for position1It hangs down for ballistic trajectory apogee distance shell hole o point maximum Straight distance, LmaxThe maximum distance of farthest horizontal distance for ballistic trajectory apart from shell hole o point both pellet injection.
When projectile impact glass, using rum point shell hole o point as track starting point, trajectory parabolic path is calculated.Due to nothing Spot speed V is hit in method measurement, needs to introduce restriction variable, various respectively to hit spot speed V, ballistic trajectory apart from shell hole below The farthest horizontal distance L of o pointmax, ballistic trajectory vertex distance shell hole o point maximum normal distance h1To limit variable, both basis Live actual environment sets each numberical range for limiting variable.
4.1: step 1.1 make shell hole sample image during, measurement record each shell hole sample angle be γ when The instantaneous stroke speed V of bullet, instantaneous velocity V under same glass material and thickness when projectile impact glass and shell hole angle The rule of γ is stablized, so when live window-glass is approximate with sample glass in material and thickness, it can be according to angle γ's Recognition result directlys adopt the speed V recorded in laboratory.Both the speed of V was drawn according to live glass material and projected angle of impact It is fixed.Then trajectory parabola is represented by following formula:
Wherein, the t in formula is time variable, indicates the value of the L and H under different time.
It willThe substitution of η value can obtain following formula:
As H=0, indicate that both ballistic trajectory has extended to ground to ballistic trajectory apart from shell hole o point vertical range for-h, this When be set as trail termination.
4.2: specific by some such as rivers, city wall such as discovery scene because particular surroundings limits according to on-site actual situations L can not then be set as limiting variable by landform barrier in place transmitting bullet of the distance greater than L, otherwise horizontal with trajectory parabola Maximum distance L is to limit variable, then trajectory parabola is represented by following formula:
T in above formulacbTrack time required for from L axis c point to b point in indicates coordinate system, in this period rail The horizontal distance that mark is passed through is Lmax
Eliminate tcb, then:
By η value, cos η=cos γ cos α, sin2η=1-cos2η,
It substitutes intoFollowing relationship can be obtained:
As H=0, indicate that both ballistic trajectory has extended to ground to ballistic trajectory apart from shell hole o point vertical range for-h, this When be set as trail termination.
4.3: in conjunction with on-site actual situations, such as finding live periphery but the special launch site of high building, plateau, with trajectory Maximum vertical height of the parabola apart from shell hole o point is to limit variable, both sets track apex height depending on the site environment, then bullet Road parabola is represented by following relationship:
T in above formulacaIndicate that track is the time required to from L axis c point to a point in a coordinate system, a point is track vertex Coordinate, in this coordinate, the speed of track vertical direction is 0, the wink according to free-fall velocity formula, in shell hole position when o point The vertical component of Shi Sudu is vcos η=gtca
Result above substitution coordinate formula can be obtained into following relationship:
η value is substituted into above formula, following relationship can be obtained:
As H=0, indicate that both ballistic trajectory has extended to ground to ballistic trajectory apart from shell hole o point vertical range for-h, this When be set as trail termination.
Step 5, ballistic trajectory direction is calculated.
As shown in figure 12, δ is the azimuth of the normal of glass planar where shell hole, plane and glass where β trajectory parabola The angle of plane where glass plane normal, plane is on the right side of plane where glass planar normal where trajectory parabola, then The parabolical azimuth of trajectory is δ+β, and plane is in the left side of plane where glass planar normal where trajectory parabola, then bullet The parabolical azimuth in road is δ-β, the as direction of trajectory parabolic path.
Step 6: with establishing emphasis personnel's household register and residence latitude and longitude coordinates database;
Information retrieval is carried out to public security internal database, is found out because gun-related case is handled or be listed in accordance with the law major tasks The illegal previous conviction personnel of object, household register address and residence location to illegal previous conviction personnel carry out inverse geocoding, obtain its family The latitude and longitude coordinates information of nationality address and residence location, and database table is established, table field includes: name, gender, citizen Identification number, docket, brief merit, the detailed location of household register, household register X-coordinate, household register Y-coordinate, the detailed location in personnel address, address X-coordinate, address Y-coordinate, as shown in table 1.
1 database table of table
Step 7: formulate three-level search range:
7-1 determines glass planar X-axis both ends azimuth:
Since glass planar and glass plane normal are to be mutually perpendicular to relationship, so the side of glass planar X-axis both direction Parallactic angle is+90 °, δ -90 ° of δ.
7-2 determines three-level search mechanisms:
The latitude and longitude coordinates of shell hole are labeled on map by 7-2-1, are set as O;
For 7-2-2 using shell hole O as the center of circle, taking a value in section (300M-400M) range is radius in map upper drawing circle.
As shown in figure 13,7-2-3 draws straight line to glass planar normal towards the azimuth direction δ, hands over circle using O as starting point In point F;Using O as starting point, straight line is drawn to the both direction vertical with azimuth δ respectively, meets at point A, B respectively with circle;It is with O Point draws straight line towards the azimuth direction δ+β or δ-β to trajectory is parabolical, meets at point D with circle;By radius OD respectively to two sides It rotates by a certain angle, meets at point C, E respectively with circle.
As shown in figure 14, region of search is arranged in 7-3:
The fan-shaped region 7-3-1 COE is set as Primary search region;
The fan-shaped region 7-3-2 BOC, fan-shaped region AOE are set as secondary search region;
The fan-shaped region 7-3-3 BOA is set as three-level region of search.
Step 8: by the longitude and latitude of the household register address of illegal previous conviction personnel or residence location or household register address and residence location In region of search on Labeling Coordinate to map, wherein the suspect of different regions of search is checked using means of different:
8-1 Primary search region is pellet injection region, and the suspect that this area marking goes out can emit bullet from place dwelling Ball hits live glass and bullet institute judges unanimously, to be defined as highest-order emission suspicion as emphasis and check object through trajectory and system.
8-2 secondary search region is shell hole glass planar towards region, and the suspect that this area marking goes out is from place dwelling Live glass can be observed, emitted by the mobile i.e. reachable emitting area of small range, be defined as intermediate transmitting suspicion.
8-3 three-level region of search is glass back towards region, and the suspect that this area marking goes out can not from place dwelling Shooting can not also observe live glass, but can observe the periphery situation in live glass location, using to scene Solution can be moved to transmitting domain and be emitted, and be defined as lower-order emission suspicion.
The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (9)

1. a kind of circle bullet impact and the case suspect of the perpendicular window in ground check method, it is characterised in that: including Following steps:
Ballistic trajectory parameter calculates:
Obtain the numerical value of the angle γ of instantaneous velocity V and glass planar;
Measure the numerical value of shell hole image long axis and glass planar y' axis angle α;
Calculate the angle beta of trajectory plane and glass plane normal;
Calculate trajectory parabolic path:
The instantaneous stroke speed V of bullet is set to limit variable, then trajectory parabola indicates are as follows:
T in formula is time variable, and g is acceleration of gravity, shot orbit horizontal axis flying distance L, the flight of the shot orbit longitudinal axis Height H, h are shell hole o point apart from ground level;
Set the horizontal maximum distance L of trajectory parabolamaxTo limit variable, then trajectory parabola indicates are as follows:
T in formula is time variable, and g is acceleration of gravity, shot orbit horizontal axis flying distance L, the flight of the shot orbit longitudinal axis Height H, h are shell hole o point apart from ground level;
Set ballistic trajectory apogee distance shell hole o point maximum normal distance h1To limit variable, then trajectory parabola indicates are as follows:
T in formula is time variable, and g is acceleration of gravity, shot orbit horizontal axis flying distance L, the flight of the shot orbit longitudinal axis Height H, h are shell hole o point apart from ground level;
Calculate ballistic trajectory direction:
Plane is on the right side of plane where glass planar normal where trajectory parabola, then the parabolical azimuth of trajectory be δ+ β, plane is in the left side of plane where glass planar normal where trajectory parabola, then the parabolical azimuth of trajectory is δ-β, δ is the azimuth of the normal of glass planar where shell hole;
With establishing emphasis personnel's household register and residence latitude and longitude coordinates database;
Formulate three-level search range;
The latitude and longitude coordinates of the household register address of illegal previous conviction personnel or residence location or household register address and residence location mark is arrived In region of search on map, wherein the suspect of different regions of search is checked using means of different.
2. circle bullet impact according to claim 1 and the case suspect of the perpendicular window in ground check method, It is characterized by: described
3. circle bullet impact according to claim 1 and the case suspect of the perpendicular window in ground check method, It is characterized by: the numerical value of the angle γ of instantaneous velocity V and glass planar obtains, steps are as follows:
Step 1: hole on glass image recognition: passing through image algorithm and the convolutional neural networks projectile impact based on TensorFlow The angle γ of instantaneous velocity and glass when glass is identified;
1.1 sample image pre-processings:
Angle γ is set as 0-90 degree, average every image increases by 1 degree, and totally 91 hole on glass images are input in machine, into Row binary conversion treatment removes variegated;
Using the outermost pixel on shell hole image four sides up and down as basic point, expands a certain distance respectively and intercept, marked Quasi- sample image;
1.2 sample image feature extractions:
Master sample image is input in the convolutional neural networks of TensorFlow and carries out the operation of n times convolution pondization, by sample Image gradually decomposes, and calculates the characteristic value of sample image;
1.3 sample images substitute into training:
The hole on glass image pattern of multiple same angle γ angles is input in machine, 1.1 and 1.2 operation is repeated, allows Machine carries out a large amount of repetition training;
Step 2: the identification of live image:
2.1 live image pre-processings:
The hole on glass image of live shooting is input in machine, image binaryzation processing and image interception is carried out, is showed The standard of field shell hole image identifies format;
2.2 live image feature extractions:
The live shell hole image of reference format is input in the convolutional neural networks of TensorFlow, n times convolution pond is carried out, Calculate the characteristic value of live shell hole image;
2.3 aspect ratios pair obtain angle γ:
By system-computed go out live shell hole image characteristic value and sample database in different angle sample image characteristic value into Row compares, and carries out result output for the highest angle of degree of fitting as the angle of live shell hole image.
4. circle bullet impact according to claim 1 and the case suspect of the perpendicular window in ground check method, It is characterized by: described includes: to public security internal data with residence latitude and longitude coordinates database with establishing emphasis personnel's household register Library carries out information retrieval, finds out because gun-related case is handled or be listed in accordance with the law the illegal previous conviction personnel of major tasks object, right The household register address of illegal previous conviction personnel and residence location carry out inverse geocoding, obtain the longitude and latitude of its household register address and residence location Coordinate information is spent, and establishes database table, table field includes: name, gender, citizenship number, docket, briefly The detailed location of merit, household register, household register X-coordinate, household register Y-coordinate, the detailed location in personnel address, address X-coordinate, address Y-coordinate.
5. circle bullet impact according to claim 1 and the case suspect of the perpendicular window in ground check method, It is characterized by: the formulation three-level search range includes:
2-1 determines glass planar X-axis both ends azimuth:
2-2 determines three-level search mechanisms:
The latitude and longitude coordinates of shell hole are labeled on map by 2-2-1, are set as O;
For 2-2-2 using shell hole O as the center of circle, certain distance is radius in map upper drawing circle;
2-2-3 draws straight line using O as starting point, to glass planar normal towards the azimuth direction δ, meets at point F with circle;Using O as starting point, Straight line is drawn to the both direction vertical with azimuth δ respectively, meets at point A, B respectively with circle;Using O as starting point, to trajectory parabola Towards the azimuth direction δ+β or δ-β draw straight line, with circle meet at point D;Radius OD is rotated by a certain angle to two sides respectively, with Circle meets at point C, E respectively;
Region of search is arranged in 2-3:
The fan-shaped region 2-3-1 COE is set as Primary search region;
The fan-shaped region 2-3-2 BOC, fan-shaped region AOE are set as secondary search region;
The fan-shaped region 2-3-3 BOA is set as three-level region of search.
6. circle bullet impact according to claim 5 and the case suspect of the perpendicular window in ground check method, It is characterized by: glass planar X-axis both ends azimuth preparation method is as follows: the azimuth of glass planar X-axis both direction It is+90 °, δ -90 ° of δ.
7. circle bullet impact according to claim 5 and the case suspect of the perpendicular window in ground check method, It is characterized by: the certain angle is set as a value in section (15 ° -20 °) range.
8. circle bullet impact according to claim 5 and the case suspect of the perpendicular window in ground check method, It is characterized by: the certain distance is set as a value in section (300M-400M) range.
9. circle bullet impact according to claim 5 and the case suspect of the perpendicular window in ground check method, It is characterized by: the suspect of the difference region of search includes: using means of different investigation
Suspect in 3-1 Primary search region, being defined as highest-order emission suspicion is that emphasis checks object;
Suspect in 3-2 secondary search region is defined as intermediate transmitting suspicion;
Suspect in 3-3 three-level region of search, is defined as lower-order emission suspicion.
CN201910572332.XA 2019-06-27 2019-06-27 Case suspicion checking method for round shot hitting window perpendicular to ground Active CN110472179B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910572332.XA CN110472179B (en) 2019-06-27 2019-06-27 Case suspicion checking method for round shot hitting window perpendicular to ground

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910572332.XA CN110472179B (en) 2019-06-27 2019-06-27 Case suspicion checking method for round shot hitting window perpendicular to ground

Publications (2)

Publication Number Publication Date
CN110472179A true CN110472179A (en) 2019-11-19
CN110472179B CN110472179B (en) 2023-07-18

Family

ID=68507297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910572332.XA Active CN110472179B (en) 2019-06-27 2019-06-27 Case suspicion checking method for round shot hitting window perpendicular to ground

Country Status (1)

Country Link
CN (1) CN110472179B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6125308A (en) * 1997-06-11 2000-09-26 The United States Of America As Represented By The Secretary Of The Army Method of passive determination of projectile miss distance
CA2484674A1 (en) * 2003-10-15 2005-04-15 Dimitri Petrov Method and apparatus for locating the trajectory of an object in motion
CN103914532A (en) * 2014-03-31 2014-07-09 安徽新华博信息技术股份有限公司 Analyzing method for automatically troubleshooting suspects
CN108036776A (en) * 2017-12-15 2018-05-15 沈阳盛之威靶业工程有限公司 It is a kind of to utilize radiographic measurement velocity of shot and the method for posture

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6125308A (en) * 1997-06-11 2000-09-26 The United States Of America As Represented By The Secretary Of The Army Method of passive determination of projectile miss distance
CA2484674A1 (en) * 2003-10-15 2005-04-15 Dimitri Petrov Method and apparatus for locating the trajectory of an object in motion
CN103914532A (en) * 2014-03-31 2014-07-09 安徽新华博信息技术股份有限公司 Analyzing method for automatically troubleshooting suspects
CN108036776A (en) * 2017-12-15 2018-05-15 沈阳盛之威靶业工程有限公司 It is a kind of to utilize radiographic measurement velocity of shot and the method for posture

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
于遨洋: ""平板玻璃上弹孔形态与射击方向的相关性检验"", 《中国刑警学院学报》 *
朱敏洁等: "实时仿真系统中外弹道仿真的一种方法", 《系统仿真学报》 *
杜周绩等: "钢珠弹丸侵彻玻璃的痕迹研究", 《科技视界》 *
汪信言: ""弹道曲线族中运动特征与几何性质"", 《数学通报》 *
田会等: "抛物线弹道弹丸飞行参数测量模型与精度分析", 《仪器仪表学报》 *
顾建农等: "旋转弹丸入水侵彻规律", 《爆炸与冲击》 *

Also Published As

Publication number Publication date
CN110472179B (en) 2023-07-18

Similar Documents

Publication Publication Date Title
CN108958292B (en) Aircraft penetration trajectory planning method based on RRT (rapid return) algorithm
CN109597423A (en) A kind of design method of the multiple constraint Terminal Guidance Laws based on reference sight angle signal
CN103984936A (en) Multi-sensor multi-feature fusion recognition method for three-dimensional dynamic target recognition
CN103996027B (en) Space-based space target recognizing method
CN105716524B (en) A kind of laser multiple target standing target measurement device and measurement method
CN106155081A (en) A kind of rotor wing unmanned aerial vehicle target monitoring on a large scale and accurate positioning method
CN106651956A (en) Method for calibrating parabolic refraction and reflection camera via single ball and parallel circle properties
CN109815810A (en) A kind of biopsy method based on single camera
WO2015184768A1 (en) Method and device for generating video abstract
CN107063267A (en) A kind of quickly localization method based on sun shadow information
CN107436434A (en) Track initiation method based on two-way Doppler estimation
CN110472180A (en) The ballistic trajectory calculation method of round bullet impact and the perpendicular window in ground
Wu [Retracted] Monitoring System of Key Technical Features of Male Tennis Players Based on Internet of Things Security Technology
CN110472185A (en) The ballistic trajectory calculation method of round bullet impact and the window at the inclined angle in ground
Qiu et al. Radio-assisted human detection
George et al. Shooter Localization using a Wireless Sensor Network of Soldier-Worn Gunfire Detection Systems.
CN110472179A (en) Round bullet impact and the case suspect of the perpendicular window in ground check method
Hiebert " CHINA'S RELATIONS WITH BURMA, MALAYSIA, AND VIETNAM".
CN110472181A (en) Round bullet impact and the case suspect of the window at the inclined angle in ground check method
Fan et al. Global contextual attention for pure regression object detection
CN105069197B (en) A kind of emulation mode of aerial platform radar detection dynamic process
CN105911524A (en) Ultra-sparse radar data association coupling method
CN110457409A (en) A kind of method that bullet beats glass case suspect search investigation
Weiss et al. Technology evaluations and performance metrics for soldier-worn sensors for assist
Wang et al. Multi Target Detection of UAV Video Based on Deep Convolution Neural Network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant