CN105259602A - Novel aspheric panoramic reflector design method - Google Patents

Novel aspheric panoramic reflector design method Download PDF

Info

Publication number
CN105259602A
CN105259602A CN201510843620.6A CN201510843620A CN105259602A CN 105259602 A CN105259602 A CN 105259602A CN 201510843620 A CN201510843620 A CN 201510843620A CN 105259602 A CN105259602 A CN 105259602A
Authority
CN
China
Prior art keywords
formula
equation
curve
sampled point
thickness
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
CN201510843620.6A
Other languages
Chinese (zh)
Other versions
CN105259602B (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.)
Heilongjiang University of Science and Technology
Original Assignee
Heilongjiang University of Science and Technology
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 Heilongjiang University of Science and Technology filed Critical Heilongjiang University of Science and Technology
Priority to CN201510843620.6A priority Critical patent/CN105259602B/en
Publication of CN105259602A publication Critical patent/CN105259602A/en
Application granted granted Critical
Publication of CN105259602B publication Critical patent/CN105259602B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B5/00Optical elements other than lenses
    • G02B5/08Mirrors
    • G02B5/10Mirrors with curved faces
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Lenses (AREA)
  • Optical Elements Other Than Lenses (AREA)
  • Closed-Circuit Television Systems (AREA)

Abstract

The invention discloses a novel aspheric reflector design method for an omnidirectional vision system. The method is characterized by comprising steps: the diameter, the thickness and the maximum incidence angle of the reflector are determined, a proper sampling point number is selected according to the maximum incidence angle and required processing precision, the diameter and the thickness are subdivided according to the sampling point number, as for a cluster of coaxial hyperbolic curves with the same focus and with real axis lengths increased linearly, curve segments are intercepted horizontally for finite element splicing to obtain a surface-type curve meeting certain requirements, an iterative formula is substituted, coordinates of the sampling points are obtained, the sampling points generated by an iterative equation are used for fitting an aspheric contour curve equation, and the equation is converted into ZMAX recognition optical equation.

Description

A kind of new aspheric panorama mirror design method
Technical field
The invention belongs to omni-directional visual process field, particularly relate to a kind of new aspheric panorama mirror design method.
Background technology
Constantly perfect along with overall view visual system, panoramic vision technology has had preliminary application in vision and optical field.In indoor and outdoor surroundings, panoramic vision is just progressively applied to the field such as monitoring and long-range reality.Due to its angular field of view of 360 degree, overall view visual system in outdoor environment the most directly application be monitoring field application.Video camera based on omnibearing vision sensor has very outstanding advantage, its looking away (360 degree), Information Compression in a hemisphere visual field can be become piece image, the quantity of information of piece image is larger, any point in panoramic picture all corresponds to a certain point determined monitored in space, therefore space position calibration algorithm is simple, without run-home during monitoring, making the algorithm of the moving object in the acquisition of visual information, detection and tracking monitoring range more simple, is a kind of quick, reliable panoramic vision information acquisition means.But current panorama system, carrying out the very large pattern distortion of panoramic picture expansion existence, propose a lot of method for reducing and eliminating this distortion, improving curved surface of reflector profile is one relatively directly effective method, greatly can improve the linearity that mirror image launches.
Summary of the invention
The nonlinear distortion that method object in the present invention exists in being to provide a kind of direct effective method abatement panoramic picture to launch, realizes the accurate reproduction of image information.
The technical solution adopted in the present invention is: the new method that a kind of non-spherical reflector for fully-directional visual system designs, and determines the diameter of catoptron, thickness and maximum incident angle three parameters.According to the machining precision of maximum incident angle and requirement, choose suitable sampling number, diameter and thickness determine the size of the useful area of minute surface, according to sampling number, diameter and thickness are segmented, substitute into iterative formula, ask for the coordinate of sampled point, carry out curve fitting according to the coordinate of all sampled points afterwards, and be converted into the Optical equations of ZMAX identification.The coaxial hyperbolic curve that real axis length linear increases progressively to cluster focal length is identical, level intercepts segment of curve and carries out finite element splicing, expands the method for the good imaging area of the linearity under the condition meeting single view constraint.
The invention has the beneficial effects as follows: improve image quality, the distortion of horizontal direction and vertical direction in removal of images unfolding calculation by the curvature of curve improving mirror surface band portions corresponding.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is curved design conditional curve schematic diagram.
Embodiment
Below in conjunction with embodiment, the present invention is further described.
Embodiment 1
The new method that a kind of non-spherical reflector for fully-directional visual system designs, determine the diameter of catoptron, thickness and maximum incident angle three parameters, according to the machining precision of maximum incident angle and requirement, choose suitable sampling number, according to sampling number, diameter and thickness are segmented, the coaxial hyperbolic curve that real axis length linear increases progressively to cluster focal length is identical, level intercepting segment of curve carries out finite element splicing and is met necessarily required type curve, substitute into iterative formula, ask for the coordinate of sampled point, the sampled point generated by iterative equation simulates aspheric contour curve equation, and be converted into the Optical equations of ZMAX identification.
Adjacent two incident ray space angles of regulation curved surface are △ θ=(β-α)/m, i.e. unlimited segmentation of catoptron actual observation angle, and m is sampling number; Separately establish △ x=r/m, for sampled point is in the step-length of minute surface horizontal direction displacement, with minute surface bottom centre for system coordinates initial point, n=1 ~ m, iteration asks for sampled point successively, and formula (3) (4) calculate the coordinate of current sampling point, and formula (7) calculates θ nafter, n adds up automatically, n=n+1, gets back to the coordinate that formula (3) calculates next sampled point, and so forth,
x 0=y 0=0(1)
θ 0=0.5×(90°-α)(2)
x n=n×△x(3)
y n=y n-1+tanθ n-1×△x(4)
θ sn=(90°-α)+n×△θ+λ n(6)
θ n=0.5×θ snn(7)
λ nfor the angle of reflection ray and main shaft, θ snbe the angle of n-th incident ray and corresponding reflection ray, the shape obtained through curve is as the equation of formula (8) form
x 2=a 1z+a 2z 2+…+a kz k+…(8)
A 1=2r 0, r 0for the radius-of-curvature on summit, the secondary surface of revolution is had:
x 2=2r 0z-(1+e)z 2(9)
Formula (9) is made into as shown in the formula
z=c 1x 2+c 2z 2(10)
Wherein: c 1=1/2r 0, c 2=(1+e)/2r 0, make x 1=x 2, x 2=z 2, turn to set { x 1i, x 2i, z imultilinear fitting, obtain following formula
z=c 1x 1+c 2x 2(11)
Derived (14) by formula (12) (13)
r 0=1/c 1(12)
e=2r 0c 2-1(13)
z = cx 2 1 + 1 - ( 1 + e ) c 2 x 2 + H ( x 2 ) - - - ( 14 )
Formula (14) is aspheric contour curve equation.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included in protection of the present invention.

Claims (2)

1. the new method that designs of the non-spherical reflector for fully-directional visual system, it is characterized in that, determine the diameter of catoptron, thickness and maximum incident angle three parameters, according to the machining precision of maximum incident angle and requirement, choose suitable sampling number, according to sampling number, diameter and thickness are segmented, the coaxial hyperbolic curve that real axis length linear increases progressively to cluster focal length is identical, level intercepting segment of curve carries out finite element splicing and is met necessarily required type curve, substitute into iterative formula, ask for the coordinate of sampled point, the sampled point generated by iterative equation simulates aspheric contour curve equation, and be converted into the Optical equations of ZMAX identification.
2. the new method of the non-spherical reflector design of fully-directional visual system according to claim 1, is characterized in that: the concrete steps of described aspheric contour curve equation are:
Adjacent two incident ray space angles of regulation curved surface are △ θ=(β-α)/m, i.e. unlimited segmentation of catoptron actual observation angle, and m is sampling number; Separately establish △ x=r/m, for sampled point is in the step-length of minute surface horizontal direction displacement, with minute surface bottom centre for system coordinates initial point, n=1 ~ m, iteration asks for sampled point successively, and formula (3) (4) calculate the coordinate of current sampling point, and formula (7) calculates θ nafter, n adds up automatically, n=n+1, gets back to the coordinate that formula (3) calculates next sampled point, and so forth,
x 0=y 0=0(1)
θ 0=0.5×(90°-α)(2)
x n=n×△x(3)
y n=y n-1+tanθ n-1×△x(4)
θ sn=(90°-α)+n×△θ+λ n(6)
θ n=0.5×θ snn(7)
λ nfor the angle of reflection ray and main shaft, θ snbe the angle of n-th incident ray and corresponding reflection ray, the shape obtained through curve is as the equation of formula (8) form
x 2=a 1z+a 2z 2+…+a kz k+…(8)
A 1=2r 0, r 0for the radius-of-curvature on summit, the secondary surface of revolution is had:
x 2=2r 0z-(1+e)z 2(9)
Formula (9) is made into as shown in the formula
z=c 1x 2+c 2z 2(10)
Wherein: c 1=1/2r 0, c 2=(1+e)/2r 0, make x 1=x 2, x 2=z 2, turn to set { x 1i, x 2i, z imultilinear fitting, obtain following formula
z=c 1x 1+c 2x 2(11)
Derived (14) by formula (12) (13)
r 0=1/c 1(12)
e=2r 0c 2-1(13)
z = cx 2 1 + 1 - ( 1 + e ) c 2 x 2 + H ( x 2 ) - - - ( 14 )
Formula (14) is aspheric contour curve equation.
CN201510843620.6A 2015-11-26 2015-11-26 A kind of aspherical panorama mirror design method Expired - Fee Related CN105259602B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510843620.6A CN105259602B (en) 2015-11-26 2015-11-26 A kind of aspherical panorama mirror design method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510843620.6A CN105259602B (en) 2015-11-26 2015-11-26 A kind of aspherical panorama mirror design method

Publications (2)

Publication Number Publication Date
CN105259602A true CN105259602A (en) 2016-01-20
CN105259602B CN105259602B (en) 2018-06-22

Family

ID=55099362

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510843620.6A Expired - Fee Related CN105259602B (en) 2015-11-26 2015-11-26 A kind of aspherical panorama mirror design method

Country Status (1)

Country Link
CN (1) CN105259602B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017214774A1 (en) * 2016-06-12 2017-12-21 合肥工业大学 Deflection film design method and liquid crystal display device
CN113552717A (en) * 2021-07-19 2021-10-26 森思泰克河北科技有限公司 Lens design method, lens and laser radar

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS56126803A (en) * 1980-03-11 1981-10-05 Kunihiko Hosaka Reflector having multifaceted reflection element
US20070002444A1 (en) * 2004-04-05 2007-01-04 Patricia Piers Ophthalmic lenses capable of reducing chromatic aberration
CN102620683A (en) * 2012-03-31 2012-08-01 中国科学院长春光学精密机械与物理研究所 Aspheric surface adjustment error compensation method for sub-aperture stitching detection

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS56126803A (en) * 1980-03-11 1981-10-05 Kunihiko Hosaka Reflector having multifaceted reflection element
US20070002444A1 (en) * 2004-04-05 2007-01-04 Patricia Piers Ophthalmic lenses capable of reducing chromatic aberration
CN102620683A (en) * 2012-03-31 2012-08-01 中国科学院长春光学精密机械与物理研究所 Aspheric surface adjustment error compensation method for sub-aperture stitching detection

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
沙定国,全书学,朱秋东,苏大图: "一种光学非球面检测技术中的数据处理问题", 《光子学报》 *
沙定国,全书学,朱秋东,苏大图: "光学非球面面形拟合方法研究", 《光子学报》 *
沙定国,郭剑锋: "光学非球面检测系统的精度仿真分析", 《宇航计测技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017214774A1 (en) * 2016-06-12 2017-12-21 合肥工业大学 Deflection film design method and liquid crystal display device
CN113552717A (en) * 2021-07-19 2021-10-26 森思泰克河北科技有限公司 Lens design method, lens and laser radar

Also Published As

Publication number Publication date
CN105259602B (en) 2018-06-22

Similar Documents

Publication Publication Date Title
CN109003311B (en) Calibration method of fisheye lens
CN109903227B (en) Panoramic image splicing method based on camera geometric position relation
CN106871787B (en) Large space line scanning imagery method for three-dimensional measurement
KR100599423B1 (en) An omnidirectional imaging apparatus
CN102509261B (en) Distortion correction method for fisheye lens
CN109859272B (en) Automatic focusing binocular camera calibration method and device
CN106679676B (en) A kind of monoscopic multifunctional optical sensor and implementation method
CN102927982B (en) Double-spectrum autonomous navigation sensor and design method of double-spectrum autonomous navigation sensor
CN106971408A (en) A kind of camera marking method based on space-time conversion thought
US20180274819A1 (en) Calibration method for heliostats
CN102080969A (en) Apparatus for rapidly correcting installation position of panoramic vision measuring system
CN105259602A (en) Novel aspheric panoramic reflector design method
CN105093523A (en) Multi-scale multi-aperture optical imaging system
CN102780834A (en) Semi-cylindrical panorama upwrapping method of pano-eye image
CN109471252B (en) Complex optical scene space motion simulation device
CN104483750A (en) Large-image-plane interscan time-sharing imaging optical system
CN112964238B (en) Device and method for improving pointing accuracy of optical telescope auxiliary radio telescope
CN107498156B (en) A kind of multi-angle electric arc light intensity acquisition method based on CCD camera
CN113920201A (en) Polar line geometric constraint fisheye camera calibration method
Xu Comparison and research of fisheye image correction algorithms in coal mine survey
CN113989105B (en) Single-camera spherical mirror reflection imaging projection device
CN105511060B (en) The global big visual field moon edge optical image-forming objective lens of face ring shape
CN209877942U (en) Image distance measuring instrument
US10511794B2 (en) Wide field of view optical module for linear sensor
CN111327796A (en) Optical imaging system based on concentric sphere secondary imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180622

Termination date: 20211126