CN102353519B - Resolving power evaluation method for three-generation dim light image intensifier - Google Patents
Resolving power evaluation method for three-generation dim light image intensifier Download PDFInfo
- Publication number
- CN102353519B CN102353519B CN 201110150196 CN201110150196A CN102353519B CN 102353519 B CN102353519 B CN 102353519B CN 201110150196 CN201110150196 CN 201110150196 CN 201110150196 A CN201110150196 A CN 201110150196A CN 102353519 B CN102353519 B CN 102353519B
- Authority
- CN
- China
- Prior art keywords
- resolving power
- target line
- image
- value
- line unit
- 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.)
- Active
Links
Images
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention discloses a resolving power measuring device and a resolving power evaluation method for a three-generation dim light image intensifier, and belongs to the field of optical measurement and metering. The resolving power measuring device is characterized by consisting of a light source component, a resolving power target, a collimator tube, an imaging objective lens, a test camera obscura, a charge coupled device (CCD) camera and a computer. The resolving power evaluation method comprises the following steps of: imaging the resolving power target irradiated by a standard light source to a fluorescent screen by using a measured image intensifier; converting into a frame image of a target line by using the CCD camera, and transmitting into the computer; successively processing a single frame image of the target line by using a normalized cross correlation model and an optical modulation degree model through internal image processing software of the computer to acquire a single frame processing result; and analyzing multi-frame processing results and performing corresponding supplement operation to acquire a final resolving power evaluation result. By the device and the method, the problem of objective evaluation during resolving power measurement of the three-generation dim light image intensifier is solved; and the device and the method can be popularized to other measurement fields such as an intensified charge coupled device (ICCD) measurement field and the like where resolving power is required to be objectively evaluated, and have wide application prospects.
Description
Technical field
The invention belongs to field of optical measuring technologies, relate generally to a kind of three generations's gleam image intensifier resolution measurement instrument, relate in particular to the evaluation method of cover three generations's gleam image intensifier resolution measurement device and three generations's gleam image intensifier resolving power.
Background technology
In recent years, along with scientific-technical progress, more and more higher to the requirement of low-light technology in fields such as space astronomy detection, fluorescence detection, military night vision reconnaissances, three generations's gleam image intensifier becomes increasingly conspicuous as this its effect of field core devices.The accurate measurement of three generations's gleam image intensifier parameter has significant meaning to improving its performance in research and production, also be an urgent demand of national defense industry numerous areas.Resolving power is the spatial frequency of gleam image intensifier MTF curve 2%~3% degree of modulation correspondence, is one of important parameter of reflection gleam image intensifier performance, is determining the low-light system 10
-31x~10
-1Operating distance and image definition during the above illumination of 1x.The development that is measured as gleam image intensifier, production and the application of gleam image intensifier resolving power provides and measured the calibration means accurately.
Traditional two generation gleam image intensifier resolution measurement method adopt the visualization method usually, Institutes Of Technology Of Nanjing, Shijiazhuang College Of Arms and China Arms Industry the 205 research institute once developed two generation gleam image intensifier resolution test equipment, all be to adopt the visualization method.Two generation gleam image intensifier resolution test equipment comprise standard sources, integrating sphere, neutral attenuator, power supply, resolving power target, parallel light tube, image-forming objective lens, gleam image intensifier to be measured, eyepiece, by resolving power target image on the eye-observation image intensifier video screen, judge the image intensifier resolving power at last.This method has measures advantage simply and intuitively, but influenced by people's subjective factor, such as to same target, different observer's observed result differences, same observer is because continue the difference of observation time length, age, health, and observed result is difference to some extent all; So visualization method accuracy of measurement is not high, repeatability is relatively poor.Therefore, the low-light level night vision device objective evaluation is very important.
At " optical technology " in the May issue, 2000, the 21st the 2nd phase of volume, among the P451-453, people such as Liu Zhengyun use gray scale-gradient co-occurrence matrix model to two generation low-light level night vision device resolution measurement image texture features analyze, with auxiliary artificial select and turning point is judged as the objective criterion of resolving power of five streak feature parameters such as gray scale entropy, gradient entropy, gray scale mean square deviation.Its algorithm computational accuracy is not high, can't provide quantitative conclusion voluntarily by algorithm, must auxiliary subjective judgement, and when in the visual field, having bright or dark spot or stripe pattern itself small-sized, can cause all that occurrence law is poor as a result.
Three generations's gleam image intensifier belongs to the photoelectronic imaging device, all the video screen imagings of its image and other that becomes are the same, have problems such as screen flicker, dark background noise, ion burn, therefore, performance and index that three generations's gleam image intensifier is new are had higher requirement to resolution measurement, press for and adopt resolving power objective evaluation algorithm development gleam image intensifier resolution measurement system, repeatability and accuracy when improving three generations's gleam image intensifier resolution measurement.
Not seeing both at home and abroad at present has the open report that adopts method for objectively evaluating to measure three generations's gleam image intensifier resolving power.
Summary of the invention
The technical problem to be solved in the present invention is at the deficiencies in the prior art, and the device of a cover measurement three generations gleam image intensifier resolving power and the method for objective evaluation three generations gleam image intensifier resolving power are provided.
Three generations's gleam image intensifier resolution measurement device provided by the invention, comprise light source assembly, the resolving power target, parallel light tube, the image-forming objective lens that has zoom function, the test camera bellows that has input and outgoing window, base platform, the ccd video camera that has function of temperature control, computing machine, described light source assembly, parallel light tube and test camera bellows are fixed on the described base platform by corresponding bracing frame, described resolving power target is connected on the parallel light tube, its target surface is positioned on the object space focal plane of parallel light tube and target surface is centered close on the optical axis of measuring light path, described image-forming objective lens and described ccd video camera are installed in described base platform by two-dimension translational mechanism and D translation mechanism respectively, the target surface center of the optical axis of image-forming objective lens and ccd video camera all is positioned on the optical axis of measuring light path, and is placed on the video screen of the image intensifier to be measured in the described test camera bellows and the object space focal plane and focus that the center correspondence is positioned at described ccd video camera; The light beam that described light source assembly sends illuminates described resolving power target, and the time of the target line pattern of resolving power target through focused on described three generations's gleam image intensifier to be measured behind the described parallel light tube collimation by described image-forming objective lens is on the pole-face; Three generations's gleam image intensifier to be measured forms brighter target line image to collimation target line pattern multiplication back at its video screen; Described ccd video camera is gathered the target line image on three generations's gleam image intensifier video screen to be measured and is converted thereof into electric signal and transfers in the described computing machine; The described ccd video camera of described computer control is finished the setting of dependence test parameter; Gather the target line image of cold-scarce scape image, hot background image and the resolving power target of ccd video camera output; The corresponding signal of gathering is carried out a series of image handle, finally obtain the resolving power evaluation result of described three generations's gleam image intensifier to be measured.
The three generations's gleam image intensifier resolving power evaluation method that adopts three generations's gleam image intensifier resolution measurement device of the present invention to realize may further comprise the steps:
The first step, the running parameter of the described ccd video camera of initialization;
In second step, according to keyboard instruction, gather the cold-scarce scape image A of described ccd video camera output in succession
LWith hot background image A
RAnd all deposit in the storer,
The 3rd step, according to keyboard instruction, gather a frame resolving power target image F of described ccd video camera output and deposit in the storer, show this image at display screen simultaneously;
In the 4th step, according to keyboard instruction, from storer, call cold-scarce scape image A
L, hot background image A
RWith resolving power target image F, resolving power target image F is deducted cold-scarce scape image A according to pursuing grey scale pixel value
LWith hot background image A
RThe gray-scale value of corresponding pixel points obtains testing image G and corresponding two dimensional gray matrix thereof, stores the two dimensional gray matrix and shows testing image G at display screen;
In the 5th step, according to keyboard commands, it is the ROI district that the image processing region is set, the zone that the ROI district limits for the highest distinguishable or inferior high distinguishable resolving power target line group of tested image intensifier;
In the 6th step, generate the standard form that has minimum resolving power target line group and have all target line unit correspondences in the inferior low resolution target line group in the ROI zone;
The 7th step, travel through all pixels in ROI zone successively in the mode of pursuing pixel line by line with each standard form that generates, simultaneously, adopt the normalized crosscorrelation formula to calculate the cross-correlation coefficient NC value of standard form and subgraph in each pixel position, find out the corresponding subgraph that is the coupling target line unit that have maximum cross correlation coefficient NC value with each standard form one by one, record these maximum cross correlation coefficient NC values and corresponding coupling target line unit thereof sequence number (i, j); Cross-correlation coefficient NC value and resolution threshold value NC with these coupling target line unit
DifferentiateCompare, in cross-correlation coefficient NC value greater than differentiating threshold value NC
DifferentiateCoupling target line unit in, will have best result and distinguish that the coupling target line unit of power is as optimum matching target line unit;
In the 8th step, in optimum matching target line unit, get the gray-scale value I of a row pixel of horizontal target line group switching centre position
1I
2... I
kAnd these gray-scale values are carried out size ordering, rejects acquisition selection gray value sequence I after k/Q the maximal value
ChoosingAnd Q is positive integer, if k/Q is not integer, then units rounds downwards, and calculates the degree of modulation of the horizontal target line of optimum matching group with following formula:
In the formula, I
AvgFor selecting the average of each gray-scale value in the gray value sequence, I
MaxFor selecting in the gray value sequence mean value greater than each gray-scale value of average, I
MinFor selecting in the gray value sequence mean value less than each gray-scale value of average; In like manner, calculate the degree of modulation of the perpendicular target line group of optimum matching; Ask the average of optimum matching target line unit horizontal target line group and perpendicular target line group degree of modulation
That is the optical modulation degree of optimum matching target line unit;
The 9th step is with the optical modulation degree of optimum matching target line unit
With threshold modulation M
cCompare: if
Continue to judge the optical modulation degree of optimum matching target line unit
Whether apparently higher than threshold value M
c, if
Then the resolving power value of optimum matching target line unit correspondence is that single frames is handled net result; If
Then calculate the optical modulation degree that resolving power is higher than the adjacent cells of optimum matching target line unit
If the optical modulation degree of high resolution adjacent target line unit
Satisfy
Then the resolving power value of high resolution adjacent target line unit correspondence is the net result that single frames is handled, otherwise the resolving power value of keeping optimum matching target line unit correspondence is that single frames is handled net result; If
Calculate the optical modulation degree that resolving power is lower than the adjacent cells of optimum matching target line unit
If satisfy
And
Then the resolving power value of low resolution adjacent target line unit correspondence is that single frames is handled net result; If satisfy
Then utilize the cross-correlation coefficient NC between normalized crosscorrelation formula calculating low resolution adjacent target line unit and the corresponding standard form thereof again
AdjacentIf, NC
Adjacent〉=0.97NC, then the resolving power value of optimum matching target line unit correspondence is that single frames is handled net result, if NC
Adjacent<0.97NC, then the resolving power value of low resolution adjacent target line unit correspondence is that single frames is handled net result and shown at display screen;
The tenth step, repeatedly repeating the 3rd step to the 9th step handles follow-up multiframe resolving power target image respectively, thereby obtain the single frames resolving power result of respective numbers, if have half and above single frames resolving power result identical, with the net result of this result as three generations's gleam image intensifier resolving power evaluation method; If the single frames resolving power result of multiple image is in a discrete distribution, then the optical modulation degree is asked by the 8th step respectively in the corresponding target line of each single frames resolving power result unit
And compare the optical modulation degree with degree of modulation criterion Mc one by one
Near M
cThe corresponding resolving power value in target line unit be considered as the net result of described three generations's gleam image intensifier resolving power evaluation method and show at display screen.
Overall technology effect of the present invention is embodied in the following aspects.
(1) the present invention has set up a cover three generations gleam image intensifier resolution measurement device by light source assembly, resolving power target, parallel light tube, image-forming objective lens, ccd video camera and computing machine, wherein, light source assembly has adopted the weak illumination light source of even diffusion, solved the problem of prior art light source lack of homogeneity under the low-light-level measurement background, measured the resolving power of three generations's gleam image intensifier for objective evaluation and lay a good foundation; In addition, used ccd video camera has refrigeration unit and temperature controller, has eliminated the electronic noise of CCD itself and thermonoise etc. greatly, for software analysis and the noise handled as three generations's low-light booster provide hardware guarantee.
(2) in the present invention, be equipped with image processing software in the computing machine, this image processing software carries out objective evaluation to the resolving power target image that becomes with ccd video camera through three generations's low-light booster in succession, its evaluation algorithms adopts the bimodel based on different principle, be normalized crosscorrelation model and optical modulation degree model, two models are independent mutually to complement each other again, obtains final conclusion by dual independent criterion.Thereby the present invention can accurately measure three generations's gleam image intensifier resolving power, the accuracy of measurement height, and good reproducibility is for development, production and the application of three generations's gleam image intensifier provides reliable technical support.
(3) in image processing process, with the matching factor in the normalized crosscorrelation model as independent criterion, not only can qualitative measurement target surface groove group whether clear, but also can search the corresponding target line of a certain particular resolution group position, be subsequent calculations lock image processing region.Thus, not only reduced the operand that image is handled, and the error result of having avoided unnecessary calculating to cause.
(4) the present invention introduces optical modulation degree model in the resolving power objective evaluation originally, because under the low contrast condition of low-light field, the image resolvability that the judgement of employing optical modulation degree criterion draws and the result of eye-observation have the height consistance, therefore, with the optical modulation degree as another independent criterion of the present invention, the resolving power data be can quantitatively obtain, thereby objectivity, accuracy and the accuracy of evaluation method of the present invention improved.
Description of drawings
Fig. 1 is the formation synoptic diagram of three generations's gleam image intensifier resolution measurement device of the present invention.
Fig. 2 is the measurement index path of three generations's gleam image intensifier resolution measurement device.
Fig. 3 is the workflow diagram of three generations's gleam image intensifier resolving power evaluation method of the present invention.
Embodiment
The present invention is described in further detail below in conjunction with accompanying drawing and preferred embodiment.
As shown in Figure 1, the preferred embodiment of three generations's gleam image intensifier resolution measurement device mainly comprises light source assembly 1, resolving power target 2, parallel light tube 3, image-forming objective lens 4, test camera bellows 5, ccd video camera 6, base platform 7, computing machine 8 and power supply 9.
Resolving power target 2 is the USAF1951 resolving power target that the low-light field tests is often used, and is of a size of 101.6mm * 82.6mm * 1.5mm, has 10 groups of target line unit on it, i.e. width such as grade and the equally spaced bright dark fringe that forms at target surface.Every group of target line unit constitutes by isometric three horizontal target lines and three vertical target lines, and the length of target line is five times of target line width, and target line width and adjacent target line equate at interval.The spacing of horizontal target line and vertical target line is the twice of target line width, and from maximum one group of target line unit, the target line unit that per two packet sizes approach is all pressed
Ratio dwindle.
According to shown in Figure 2, the light beam that standard sources 1-1 sends is in succession through forming uniform diffusion light beam and illuminating resolving power target 2 behind neutral colour filter and iris 1-2 and the integrating sphere 1-3; Whole target lines unit on the resolving power target 2 is focused on the time of the three generations's gleam image intensifier 5-1 to be measured that tests in the camera bellows 5 on the pole-face by image-forming objective lens 4 through parallel light tube 3 collimation back; Three generations's gleam image intensifier 5-1 to be measured forms brighter target line image to collimation target line pattern multiplication back at its video screen; Ccd video camera 6 is gathered the target line image on the video screen and is converted thereof into electric signal and transfers in the computing machine 8.
Ccd video camera 6 is selected the COOLSNAP K4 video camera of U.S. Photometric company production for use and is subjected to the control of computing machine 8, and this video camera has refrigeration unit and temperature controller, the highest sample frequency 20MHz; Pixel Dimensions 7.4 μ m * 7.4 μ m, pixel count is 2048 * 2048, chip area 15.16m * 15.16mm, frame rate 3fps.In running order following at ccd video camera 6, the image that collects when its shutter close is cold-scarce scape image, and this cold-scarce scape image is real to be ccd video camera 6 intrinsic heat noise own and electronic noises; Open and standard sources 1-1 power supply when not opening when its shutter, its image that collects is hot background image, and this image is produced by experimental enviroment; Open and standard sources 1-1 when opening when its shutter, its image that collects is the target line image.Cold-scarce scape image, hot background image and target line image are bitmap (bitmap) form.
Resolving power evaluation method of the present invention realizes that by computing machine 8 after three generations's gleam image intensifier 5-1 to be measured was placed into the resolution measurement device, image processing software was according to the following operation steps of flow performing shown in Figure 3.
The first step, initialization
After start powers on, load the driver of ccd video camera, capture card, mouse, keyboard, initialization apparatus.In conjunction with the input data setting ccd video camera 6 of keyboard and the duty of capture card, as the on-keyboard input, then call ccd video camera and capture card initialization status data in the storer, set the duty of the two; This acquiescence duty is focus pattern and video monitoring mode, and lightness and contrast parameter are 50, and working temperature is-25 ℃, and frequency acquisition was 3 frame/seconds, and drainage pattern is the single frames collection, and transfer rate is 20M/s;
In second step, gather cold-scarce scape image and hot background image
After the acquisition of receiving the keyboard input, at first send cold-scarce scape acquisition to ccd video camera 6, gather the cold-scarce scape image A of ccd video camera 6 outputs by pci bus
LAnd deposit in the storer; Then, send hot background acquisition to ccd video camera 6, under the situation that light source 1-1 does not open, gather the hot background image A of ccd video camera 6 outputs by pci bus
RAnd deposit in the storer.
In the 3rd step, gather resolving power target image
After the acquisition of receiving the keyboard input, send resolving power target acquisition to ccd video camera 6, under the situation that light source 1-1 opens, gather a frame resolving power target image F of ccd video camera 6 outputs and deposit in the storer by pci bus, show this resolving power target image F at display screen simultaneously.
The 4th step, the image pre-service
In the present invention, the image pre-service refers to resolving power target image F according to deduct cold-scarce scape image A by grey scale pixel value
LWith hot background image A
RA kind of processing mode of corresponding pixel points gray-scale value, its specific algorithm formula is as follows:
G (s, t)=n (s, t)-n
L(s, t)-n
R(s, t) in the formula, n (s, t), n
L(s, t), n
R(s t) is followed successively by resolving power target image F, cold-scarce scape image A
L, hot background image A
RPixel (s, gray-scale value t), g (s, t) for the testing image G that obtains after the pre-service at pixel (s, gray-scale value t).After whole two field picture pre-service finishes, just obtained the two dimensional gray matrix of testing image G and with this matrix stores in storer, simultaneously show testing image G at display screen.
In the 5th step, set the ROI zone
Behind the ROI zone setting command of receiving the keyboard input, calculate the size in ROI zone according to corresponding formula.The ROI zone is the perform region, refers to the image-region that need handle, and this zone is the square area that the group of the highest distinguishable resolving power target of tested image intensifier limits.Since present embodiment at the resolving power target be the USAF1951 target, this target comprises-2 ,-1,0, totally 10 groups of 1......6,7, the length of side L in the ROI zone of every group of correspondence
iDetermined by following formula:
In the formula, i represents the group of resolving power target, i.e. i=-2 ,-1,0,1......6,7.At present, the three-generation image enhancer best result is distinguished power generally at second group, and therefore, in the present embodiment, correspondence is got i=2, obtains L by calculating
2=3.56.Calculated amount is reduced in the target line zone that definition ROI can avoid handling low resolution in the zone.
In the 6th step, generate target line unit standard form
Standard form refers to target line number and the shape and size generated ideal target line image according to each discrete target line unit on the used resolving power target.
For the present embodiment, the width of standard form is determined by following formula:
In the formula, j represents each target line unit number in every group of resolution target line.In the present embodiment, minimum resolving power target line group in the ROI zone is second group, i.e. i=2, and the span of j is 1,2,3......6, the target line unit that has maximum line width in the corresponding resolving power group of j=1, the target line unit that has minimum feature in the corresponding resolving power group of j=6.For USAF1951 resolving power target, its standard form height is 5/12 of width, and three horizontal three perpendicular rectangular target lines are arranged in the template, six roots of sensation target line equal and opposite in direction, and target line length breadth ratio is 5: 1, target line width E
I, j(mm of unit) determined by following formula:
E
I, jAlso be simultaneously the spacing between three horizontal target lines being equally spaced or the three perpendicular target lines, the distance between horizontal target line group and the perpendicular target line group is 2E
I, j
According to above-mentioned computing formula and corresponding target line parameters relationship, generate the standard form that has minimum resolving power target line group and have all target line unit correspondences in the inferior low resolution target line group in the ROI zone.For present embodiment, generate i=2 and j=1~6 successively, the standard form of each unit correspondence of i=3 and j=1~6.
In the 7th step, calculate the normalized crosscorrelation coefficient
The standard form that generates is placed the upper left corner, ROI zone of testing image G, and according to moving this standard form by the mode of pixel line by line, every movement once, utilize the normalized crosscorrelation formula to calculate normalized crosscorrelation coefficient (abbreviation cross-correlation coefficient) NC of current standard form and subgraph, all pixels until current standard form traversal ROI zone.Subgraph is testing image G by the part that standard form covered.Find out the subgraph of corresponding maximum cross correlation coefficient NC value then, record maximum cross correlation coefficient NC value and corresponding subgraph position thereof.Specific practice is, according to resolving power standard form order from low to high, travel through all pixels in ROI zone one by one with standard form, finally find out the corresponding subgraph that has maximum cross correlation coefficient NC value with each standard form, that is coupling target line unit, record these maximum cross correlation coefficient NC values and the corresponding target line of corresponding subgraph unit number thereof (i, j).
Cross-correlation coefficient NC value is more big, illustrate that subgraph is more similar to standard form, when cross-correlation coefficient NC value is 1, it is the perfect match position, and in actual conditions, cross-correlation coefficient NC can not be 1, therefore only need find at testing image G to have the corresponding subgraph of maximum cross correlation coefficient NC value position, can think that then the target line unit of this position is exactly the coupling target line unit of current standard form.Like this, the corresponding coupling target line unit of each standard form.Cross-correlation coefficient NC value and resolution threshold value NC with these coupling target line unit
DifferentiateCompare, in cross-correlation coefficient NC value greater than differentiating threshold value NC
DifferentiateCoupling target line unit in, find out and have that best result is distinguished the target line unit of power and as optimum matching target line unit.
In these coupling target line unit corresponding with whole standard forms, its cross-correlation coefficient NC value too hour, this target line unit is actually and can't differentiates, and cross-correlation coefficient NC value is when too big, the resolving power that can not truly reflect tested image intensifier, therefore a suitable cross-correlation coefficient NC value need be set namely differentiates threshold value NC
Differentiate, to find out not only distinguishable and resolving power but also high target line unit.Usually, need at first obtain not only distinguishable and resolving power but also high target line unit by a large amount of experiments, then, with this unit corresponding cross-correlation coefficient NC value as differentiating threshold value NC
DifferentiateNC
DifferentiateSpan is generally 0.5~0.7, NC in the present embodiment
DifferentiateGet 0.68.
Because tested gleam image intensifier is the video screen imaging, there are the bright phaeism flicker of intrinsic screen and screen shake scintillation.For this reason, the present invention utilize normalized crosscorrelation model accuracy height, adaptability good, to the insensitive characteristics of the linear transformation of gradation of image value, subgraph to the testing image G of standard form and covering thereof carries out matching operation, thereby can reduce the influence that the bright phaeism flicker of tested image intensifier screen is estimated resolution measurement; In addition, because the result of calculation of normalized crosscorrelation model is to find the solution optimum matching target line unit, therefore, can solve the influence that screen shake flicker is estimated resolution measurement well.
The 8th step, the calculating optical degree of modulation
The size of optical modulation degree has reflected the resolvability of image to a great extent, can judge the readability of image on the whole.The sharpness of different target lines and the optical modulation degree of its grey scale curve are proportional.Simultaneously, the clear picture degree of the size of optical modulation degree and human eye observation has good consistance, and the optical modulation degree is more high, and image is more clear, and vice versa.Among the present invention, the optical modulation degree calculates according to the following steps:
8.1 in optimum matching target line unit, get the gray-scale value I of a row pixel of horizontal target line group switching centre position
1I
2... I
kAnd these gray-scale values are carried out size ordering, rejects acquisition selection gray value sequence I after k/Q the maximal value
ChoosingAnd Q is positive integer.If k/Q is not integer, then units rounds downwards.The Q value depends on the fluoroscopic ion burn quantity of tested image intensifier grade, and its size should guarantee effectively to reject ion burn, guarantees as far as possible that again leaving more pixel count is used for subsequent treatment.In the present embodiment, get Q=12.
To selecting gray value sequence I
ChoosingAsk the arithmetic mean value to obtain I
AvgTo select gray value sequence I
ChoosingIn greater than average I
AvgGray-scale value ask arithmetic mean to obtain I
Max, and to less than I
AvgGray-scale value ask the arithmetic mean value to obtain I
MinAfterwards, calculate the degree of modulation of horizontal target line group with following formula:
8.2 get the gray-scale value of perpendicular target line group switching centre position, optimum matching target line unit one-row pixels, the method with 8.1 is calculated the degree of modulation of perpendicular target line group.
8.3 ask the average of optimum matching target line unit horizontal target line group and perpendicular target line group degree of modulation
It is the optical modulation degree of optimum matching target line unit.Optical modulation degree with this target line unit
With the threshold modulation M in the storer
cCompare, if
Show that this optimum matching target line unit is distinguishable under the degree of modulation criterion, then enter 9.1 steps; If
Show that this optimum matching target line unit can not differentiate under the degree of modulation criterion, then enter 9.2 steps.
Threshold modulation M
cIt is through test of many times and the resolution critical point of determining after verifying repeatedly.Generally speaking, degree of modulation is lower than 0.1, and resolving power target image will be difficult to differentiate, thereby, threshold modulation M
cValue generally at 0.1≤M
c≤ 0.2 scope.
Because three generations's gleam image intensifier is the video screen imaging, has intrinsic ion burn.For this reason, the present invention by removing k/Q maximum gradation value, can effectively reject the big unusual gray-scale value that ion burn causes, thereby remove the interference that ion burn is estimated resolution measurement in the calculating of optical modulation degree.
In the 9th step, obtain single-frame images resolving power result calculated
9.1 work as
Continue to judge the optical modulation degree of optimum matching target line unit
Whether apparently higher than threshold value M
c, if
Then the sequence number of optimum matching target line unit is the single frames result; If
Then calculate the optical modulation degree that resolving power is higher than the adjacent cells of optimum matching target line unit
If the optical modulation degree of high resolution adjacent target line unit
Satisfy
Then the sequence number of this high resolution target line unit is the single frames result, otherwise the sequence number of keeping optimum matching target line unit is the single frames result.
9.2 work as
Calculate the optical modulation degree that resolving power is lower than the adjacent cells of optimum matching target line unit
If satisfy
And
Then the sequence number of this low resolution adjacent target line unit is the single frames result; If satisfy
Then further utilize the cross-correlation coefficient NC between normalized crosscorrelation formula calculating low resolution adjacent target line unit and the corresponding standard form thereof
AdjacentIf, cross-correlation coefficient NC
AdjacentThere are not bigger variation, i.e. NC with the cross-correlation coefficient NC of optimum matching target line unit more yet
Adjacent〉=0.97NC, then the sequence number of optimum matching target line unit is the single frames result, if NC
Adjacent<0.97NC, then the sequence number of low resolution adjacent target line unit is the single frames result.
The sequence number of above-mentioned resolving power target line unit shows as i and organizes the j unit, also need to search the corresponding resolving power value of this sequence number by (table 1) mode of tabling look-up, and with this resolving power value as the final process result of single frames and be presented on the screen.
In the tenth step, obtain multiple image resolving power result calculated
For making result of calculation more accurate, one width of cloth resolving power target image is carried out the multiframe collection, usually select for use 10 frames to get final product, then repeatedly the 3rd step of (10 times) repetition to the 9th step handles every frame resolving power target image respectively, thereby obtains the single frames resolving power result of respective numbers.Like this, will there be following two kinds of situations:
1). there are half and above single frames resolving power final process result identical, can be with the net result of this result as three generations's gleam image intensifier resolving power evaluation method.
2). the single frames resolving power final process result of multiple image is in a discrete distribution, and does not namely have a kind of single frames resolving power final process result ratio to surpass half, at this moment, the optical modulation degree is asked by the 8th step respectively in the corresponding target line of each single frames resolving power result unit
And one by one with degree of modulation criterion M
cRelatively, optical modulation degree
Near M
cThe resolving power value of target line unit correspondence be considered as the net result of three generations's gleam image intensifier resolving power evaluation method.
Table 1 resolving power target line unit and resolving power value corresponding tables
Unit: line is right
Claims (2)
1. three generations's gleam image intensifier resolving power evaluation method, this method realizes at three generations's gleam image intensifier resolution measurement device, described three generations's gleam image intensifier resolution measurement device comprises light source assembly (1), resolving power target (2), parallel light tube (3), the image-forming objective lens (4) that has zoom function, the test camera bellows (5) that has input and outgoing window, base platform (7), the ccd video camera (6) that has function of temperature control, computing machine (8), described light source assembly (1), parallel light tube (3) and test camera bellows (5) are by corresponding bracing frame (7-1,7-2,7-4) be fixed on the described base platform (7), described resolving power target (2) is connected on the parallel light tube (3), its target surface is positioned on the object space focal plane of parallel light tube (3) and target surface is centered close on the optical axis of measuring light path, described image-forming objective lens (4) and described ccd video camera (6) are installed in described base platform (7) by two-dimension translational mechanism (7-3) and D translation mechanism (7-5) respectively, the target surface center of the optical axis of image-forming objective lens (4) and ccd video camera (6) all is positioned on the optical axis of measuring light path, and is placed on the video screen of the image intensifier to be measured (5-1) in the described test camera bellows (5) and the object space focal plane and focus that the center correspondence is positioned at described ccd video camera (6); The light beam that described light source assembly (1) sends illuminates described resolving power target (2), and the time of the target line pattern of resolving power target (2) through focused on described three generations's gleam image intensifier to be measured (5-1) behind described parallel light tube (3) collimation by described image-forming objective lens (4) is on the pole-face; Three generations's gleam image intensifier to be measured (5-1) forms brighter target line image to collimation target line pattern multiplication back at its video screen; Described ccd video camera (6) is gathered the target line image on three generations's gleam image intensifier to be measured (5-1) video screen and is converted thereof into electric signal;
It is characterized in that: described computing machine (8) is estimated the resolving power of three generations's gleam image intensifier by following operation steps:
The first step, the running parameter of the described ccd video camera of initialization (6);
In second step, according to keyboard instruction, gather the cold-scarce scape image A of described ccd video camera (6) output in succession
LWith hot background image A
RAnd all deposit in the storer;
The 3rd step, according to keyboard instruction, gather a frame resolving power target image F of described ccd video camera (6) output and deposit in the storer, show this image at display screen simultaneously;
In the 4th step, according to keyboard instruction, from storer, call cold-scarce scape image A
L, hot background image A
RWith resolving power target image F, resolving power target image F is deducted cold-scarce scape image A according to pursuing grey scale pixel value
LWith hot background image A
RThe gray-scale value of corresponding pixel points obtains testing image G and corresponding two dimensional gray matrix thereof, stores the two dimensional gray matrix and shows testing image G at display screen;
In the 5th step, according to keyboard commands, it is the ROI district that the image processing region is set, the zone that the ROI district limits for the highest distinguishable or inferior high distinguishable resolving power target line group of tested image intensifier;
In the 6th step, generate the standard form that has minimum resolving power target line group and have all target line unit correspondences in the inferior low resolution target line group in the ROI zone;
The 7th step, travel through all pixels in ROI zone successively in the mode of pursuing pixel line by line with each standard form that generates, simultaneously, adopt the normalized crosscorrelation formula to calculate the cross-correlation coefficient NC value of standard form and subgraph in each pixel position, find out the corresponding subgraph that is the coupling target line unit that have maximum cross correlation coefficient NC value with each standard form one by one, record these maximum cross correlation coefficient NC values and corresponding coupling target line unit thereof sequence number (i, j); Cross-correlation coefficient NC value and resolution threshold value NC with these coupling target line unit
DifferentiateCompare, in cross-correlation coefficient NC value greater than differentiating threshold value NC
DifferentiateCoupling target line unit in, will have best result and distinguish that the coupling target line unit of power is as optimum matching target line unit;
In the 8th step, in optimum matching target line unit, get the gray-scale value I of a row pixel of horizontal target line group switching centre position
1I
2... I
kAnd these gray-scale values are carried out size ordering, rejects acquisition selection gray value sequence I after k/Q the maximal value
ChoosingAnd Q is positive integer, if k/Q is not integer, then units rounds downwards, and calculates the degree of modulation of the horizontal target line of optimum matching group with following formula:
In the formula, I
MaxFor selecting in the gray value sequence mean value greater than each gray-scale value of selecting the gray value sequence average, I
MinFor selecting in the gray value sequence mean value less than each gray-scale value of selecting the gray value sequence average; In like manner, calculate the degree of modulation of the perpendicular target line group of optimum matching; Ask the average of optimum matching target line unit horizontal target line group and perpendicular target line group degree of modulation
That is the optical modulation degree of optimum matching target line unit;
The 9th step is with the optical modulation degree of optimum matching target line unit
With threshold modulation M
cCompare: if
Continue to judge the optical modulation degree of optimum matching target line unit
Whether apparently higher than threshold value M
c, if
Then the resolving power value of optimum matching target line unit correspondence is that single frames is handled net result; If
Then calculate the optical modulation degree that resolving power is higher than the adjacent cells of optimum matching target line unit
If the optical modulation degree of high resolution adjacent target line unit
Satisfy
Then the resolving power value of high resolution adjacent target line unit correspondence is the net result that single frames is handled, otherwise the resolving power value of keeping optimum matching target line unit correspondence is that single frames is handled net result; If
Calculate the optical modulation degree that resolving power is lower than the adjacent cells of optimum matching target line unit
If satisfy
And
Then the resolving power value of low resolution adjacent target line unit correspondence is that single frames is handled net result; If satisfy
Then utilize the cross-correlation coefficient NC between normalized crosscorrelation formula calculating low resolution adjacent target line unit and the corresponding standard form thereof again
AdjacentIf, NC
Adjacent〉=0.97NC, then the resolving power value of optimum matching target line unit correspondence is that single frames is handled net result, if NC
Adjacent<0.97NC, then the resolving power value of low resolution adjacent target line unit correspondence is that single frames is handled net result and shown at display screen;
The tenth step, repeatedly repeating the 3rd step to the 9th step handles follow-up multiframe resolving power target image respectively, thereby obtain the single frames resolving power result of respective numbers, if have half and above single frames resolving power result identical, with the net result of this result as three generations's gleam image intensifier resolving power evaluation method; If the single frames resolving power result of multiple image is in a discrete distribution, then the optical modulation degree is asked by the 8th step respectively in the corresponding target line of each single frames resolving power result unit
And one by one with degree of modulation criterion M
cRelatively, optical modulation degree
Near M
cThe corresponding resolving power value in target line unit be considered as the net result of described three generations's gleam image intensifier resolving power evaluation method and show at display screen.
2. three generations's gleam image intensifier resolving power evaluation method according to claim 1 is characterized in that: described resolution threshold value NC
DifferentiateValue in 0.5~0.7 scope; Described threshold modulation M
cValue at 0.1≤M
c≤ 0.2 scope.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110150196 CN102353519B (en) | 2011-06-03 | 2011-06-03 | Resolving power evaluation method for three-generation dim light image intensifier |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110150196 CN102353519B (en) | 2011-06-03 | 2011-06-03 | Resolving power evaluation method for three-generation dim light image intensifier |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102353519A CN102353519A (en) | 2012-02-15 |
CN102353519B true CN102353519B (en) | 2013-07-10 |
Family
ID=45577130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110150196 Active CN102353519B (en) | 2011-06-03 | 2011-06-03 | Resolving power evaluation method for three-generation dim light image intensifier |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102353519B (en) |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798514A (en) * | 2012-08-08 | 2012-11-28 | 中国科学院长春光学精密机械与物理研究所 | Method for detecting imaging resolution of aerial camera in high-altitude environment condition |
CN103048605B (en) * | 2012-12-22 | 2015-10-14 | 三峡电力职业学院 | The detection screening technique that a kind of LED is aging |
CN103105286B (en) * | 2013-01-24 | 2015-09-23 | 中国兵器工业第二0五研究所 | Imaging electric system spectral response nonuniform measurement method |
CN103268020B (en) * | 2013-05-10 | 2015-11-18 | 中国人民解放军济南军区72465部队 | Hologram diffraction formula low-light level night vision device resolving power inspection method and device |
CN103826118B (en) * | 2014-03-18 | 2016-02-03 | 西安工业大学 | A kind of color TV imaging system spatial discrimination force checking device |
CN104065957B (en) * | 2014-06-13 | 2016-08-24 | 中国科学院西安光学精密机械研究所 | Camera resolving power testing device and method |
CN104296968B (en) * | 2014-10-10 | 2016-12-07 | 中国科学院长春光学精密机械与物理研究所 | The modulation transfer function test method of multichannel CCD |
CN104614154B (en) * | 2015-02-03 | 2017-06-27 | 南京理工大学 | A kind of radiance gain test device of ultraviolet imaging enhancer |
CN105203458A (en) * | 2015-10-23 | 2015-12-30 | 无锡溥汇机械科技有限公司 | Image processing device of blood analyzer |
CN105319052A (en) * | 2015-11-30 | 2016-02-10 | 武汉光驰科技有限公司 | Low-light-level image intensifier characteristic test system |
CN105391998B (en) * | 2015-12-24 | 2017-05-24 | 无锡市星迪仪器有限公司 | Automatic detection method and apparatus for resolution of low-light night vision device |
CN105973570B (en) * | 2016-04-27 | 2018-07-31 | 西安应用光学研究所 | Low-light ICCD resolution measurements device and measurement method |
CN105929351B (en) * | 2016-06-03 | 2019-06-04 | 中国电力科学研究院 | Test device and the test method of ultraviolet imager sensitivity |
CN106053031B (en) * | 2016-06-29 | 2019-01-04 | 南京理工大学 | A kind of gleam image intensifier single tube resolving power testing device |
CN107560828A (en) * | 2016-07-01 | 2018-01-09 | 北京振兴计量测试研究所 | A kind of gleam image intensifier signal to noise ratio calibrating installation |
CN106226037A (en) * | 2016-08-20 | 2016-12-14 | 南京理工大学 | A kind of large aperture gleam image intensifier imagewise uniform system safety testing device |
CN108303236B (en) * | 2017-01-11 | 2019-12-17 | 福建省莆田市中涵机动力有限公司 | Laser scanning head installation calibration inspection system |
CN107421722A (en) * | 2017-07-01 | 2017-12-01 | 南京理工大学 | Based on CCD camera image intensifier modulation transfer function test method and system |
CN107449587B (en) * | 2017-08-16 | 2023-05-05 | 中国科学院国家天文台 | Device and method for testing non-uniformity in detector pixel |
CN107884160A (en) * | 2017-09-25 | 2018-04-06 | 杭州浙大三色仪器有限公司 | Virtual image photoelectric measuring instrument |
CN107884159A (en) * | 2017-09-25 | 2018-04-06 | 杭州浙大三色仪器有限公司 | virtual image display device photoelectric measuring device |
CN108801594B (en) * | 2018-05-24 | 2024-04-30 | 中国人民解放军陆军沈阳军事代表局驻长春地区军事代表室 | Quantitative evaluation device for image quality of boundary temperature condition photoelectric imaging system |
CN109151451B (en) * | 2018-05-28 | 2020-09-25 | 上海久航电子有限公司 | Binocular night vision device resolution detection method |
CN108982990B (en) * | 2018-05-30 | 2020-12-04 | 中国人民解放军陆军工程大学 | Image intensifier tester |
CN109141829B (en) * | 2018-07-23 | 2020-10-09 | 北京大恒图像视觉有限公司 | Target line definition calculation and industrial camera rear intercept detection method |
CN109002823B (en) * | 2018-08-09 | 2020-11-10 | 歌尔科技有限公司 | Region-of-interest determining method, device, equipment and readable storage medium |
TWI669962B (en) * | 2018-12-07 | 2019-08-21 | 致伸科技股份有限公司 | Method for detecting camera module |
CN110346120B (en) * | 2019-08-05 | 2021-03-09 | 北方夜视技术股份有限公司 | System and method for testing strong light resolution and dynamic range of automatic gate control image intensifier |
CN110657959A (en) * | 2019-10-22 | 2020-01-07 | 长春军晟科技有限公司 | Universal device for detecting performance of photoelectric night vision instrument |
CN111610002B (en) * | 2020-05-27 | 2021-11-05 | 北方夜视技术股份有限公司 | Method for measuring cathode close-proximity distance of image intensifier |
CN113834635B (en) * | 2020-06-24 | 2024-08-16 | 浙江宇视科技有限公司 | Image acquisition virtual focus test method, device, equipment and storage medium |
CN112763189A (en) * | 2020-12-24 | 2021-05-07 | 松山湖材料实验室 | Measuring device for EBCMOS resolution parameter |
CN113375908B (en) * | 2021-05-26 | 2022-09-23 | 北方夜视技术股份有限公司 | Dark box system for automatic test of low-light-level image intensifier |
CN113432839B (en) * | 2021-06-09 | 2022-09-16 | 北方夜视技术股份有限公司 | System and method for comprehensively testing image quality of low-light-level image intensifier |
CN113432833B (en) * | 2021-06-15 | 2022-09-16 | 北方夜视技术股份有限公司 | Device and method for testing stability of photo-cathode of image intensifier tube after illumination |
CN113686551A (en) * | 2021-08-20 | 2021-11-23 | 北方夜视技术股份有限公司 | Portable image intensifier resolution measuring device and measuring method |
CN113781414B (en) * | 2021-08-27 | 2024-07-26 | 歌尔光学科技有限公司 | Lens resolution testing method and device and electronic equipment |
CN113873231B (en) * | 2021-09-26 | 2024-05-03 | 江西盛泰精密光学有限公司 | Monitoring system and method for baking camera module |
CN114935446A (en) * | 2022-05-05 | 2022-08-23 | 西安应用光学研究所 | Miniaturized low-light-level night vision device resolution and field-of-view test system |
CN114770067B (en) * | 2022-06-21 | 2022-09-02 | 昆明昆科测控技术有限公司 | Equipment suitable for picture tube and sleeve automatic assembly |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1384345A (en) * | 2001-04-28 | 2002-12-11 | 慧生科技股份有限公司 | Automatic measurement system and method for optical characteristics of optical unit |
CN1648630A (en) * | 2005-02-04 | 2005-08-03 | 沈阳敏像科技有限公司 | Optical detecting method and device for mini photographic module |
CN101308059A (en) * | 2008-07-08 | 2008-11-19 | 中国科学院长春光学精密机械与物理研究所 | Ultraviolet optics instrument resolution tester |
JP4204342B2 (en) * | 2003-02-07 | 2009-01-07 | 株式会社Adeka | Flame retardant synthetic resin composition |
CN101923000A (en) * | 2010-07-13 | 2010-12-22 | 中国兵器工业第二〇五研究所 | Optical measuring device with high reflectivity and high transmissivity |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02146571A (en) * | 1988-11-29 | 1990-06-05 | Fuji Xerox Co Ltd | Modulation measuring instrument for optical image |
JPH04204342A (en) * | 1990-11-30 | 1992-07-24 | Meisho Koki:Kk | Device for measuring resolution of binoculars |
-
2011
- 2011-06-03 CN CN 201110150196 patent/CN102353519B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1384345A (en) * | 2001-04-28 | 2002-12-11 | 慧生科技股份有限公司 | Automatic measurement system and method for optical characteristics of optical unit |
JP4204342B2 (en) * | 2003-02-07 | 2009-01-07 | 株式会社Adeka | Flame retardant synthetic resin composition |
CN1648630A (en) * | 2005-02-04 | 2005-08-03 | 沈阳敏像科技有限公司 | Optical detecting method and device for mini photographic module |
CN101308059A (en) * | 2008-07-08 | 2008-11-19 | 中国科学院长春光学精密机械与物理研究所 | Ultraviolet optics instrument resolution tester |
CN101923000A (en) * | 2010-07-13 | 2010-12-22 | 中国兵器工业第二〇五研究所 | Optical measuring device with high reflectivity and high transmissivity |
Non-Patent Citations (4)
Title |
---|
一种新的微光像增强器综合测试系统;刘宇等;《应用光学》;19970625;第18卷(第3期);表1及图1 * |
刘宇等.一种新的微光像增强器综合测试系统.《应用光学》.1997,第18卷(第3期), |
宽光谱像增强器分辨力测试技术研究;王守为;《信息科技辑》;20081115;第4页及图1.2.3和图1.2.4 * |
王守为.宽光谱像增强器分辨力测试技术研究.《信息科技辑》.2008, |
Also Published As
Publication number | Publication date |
---|---|
CN102353519A (en) | 2012-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102353519B (en) | Resolving power evaluation method for three-generation dim light image intensifier | |
CN102175701B (en) | System and method for online flaw detection of industrial X-ray machine | |
CN105547342B (en) | Industrial lens test device and method based on liquid crystal panel | |
CN106596073A (en) | Method and system for detecting image quality of optical system, and testing target plate | |
CN104634449B (en) | Low-light ICCD signal-to-noise test system and method for testing | |
CN110324611B (en) | Camera module detection system and detection method | |
CN107403177A (en) | Brightness measurement method based on industrial camera | |
CN109767425B (en) | Machine vision light source uniformity evaluation device and method | |
CN107884414B (en) | System and method for detecting surface defects of mirror surface object by eliminating influence of dust | |
CN110958450B (en) | Imaging system space testing device, contrast and frequency testing method | |
CN107025891B (en) | A kind of display module defect fast repairing method and system | |
Koshti | Assessing visual and system flaw detectability in nondestructive evaluation | |
CN109632269A (en) | Method based on image grayscale infomation detection diffractive-optical element performance | |
CN106248634A (en) | Fruit exocuticle glossiness measurement apparatus and method | |
CN109709460A (en) | A kind of blind type ultraviolet imagery system Photoncounting methods | |
CN101915734A (en) | Line scanning collecting method and device for magnetic trace images in magnetic particle inspection | |
CN102998026A (en) | Automatic detection device for liquid glass thermometers | |
CN106500577A (en) | A kind of clinac vane grating method for detecting position | |
CN103383353A (en) | Optical vortices-based dynamic speckle test method | |
CN210867988U (en) | Imaging system space testing device | |
CN105894482B (en) | Balance the method and system of luminance difference between each visual point image | |
CN111479097B (en) | Scattering lens imaging system based on deep learning | |
CN204044085U (en) | A kind of foreign matter detection system based on light spectrum image-forming | |
CN206906245U (en) | A kind of Regional Rainfall uniformity measurement apparatus based on laser refraction principle | |
Wang et al. | Objective evaluation of the resolution of low-light-level image intensifiers based on fast Fourier transform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |