JP2009038439A - Imaging method performing spatial filtering and its imaging apparatus - Google Patents

Imaging method performing spatial filtering and its imaging apparatus Download PDF

Info

Publication number
JP2009038439A
JP2009038439A JP2007198603A JP2007198603A JP2009038439A JP 2009038439 A JP2009038439 A JP 2009038439A JP 2007198603 A JP2007198603 A JP 2007198603A JP 2007198603 A JP2007198603 A JP 2007198603A JP 2009038439 A JP2009038439 A JP 2009038439A
Authority
JP
Japan
Prior art keywords
image data
image
imaging
optical system
numerical aperture
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2007198603A
Other languages
Japanese (ja)
Inventor
Kenji Amaya
賢治 天谷
Masashi Uejima
正史 上島
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.)
Tokyo Institute of Technology NUC
Original Assignee
Tokyo Institute of Technology NUC
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 Tokyo Institute of Technology NUC filed Critical Tokyo Institute of Technology NUC
Priority to JP2007198603A priority Critical patent/JP2009038439A/en
Publication of JP2009038439A publication Critical patent/JP2009038439A/en
Pending legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide an imaging apparatus in which the cut-off frequency when the low frequency components of a captured image are removed is set such that the high frequency noise, e.g. an quantization error, can be reduced effectively and surely. <P>SOLUTION: An image processing means for picking up an image by receiving the light from an optical system and processing the imaging output, i.e. image data, performs spatial filtering for removing the components of higher frequency than a cut-off frequency dependent on the numerical aperture of an optical system from the image data. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

本発明は顕微鏡等の高倍率の画像撮影に好適な撮像方法及びその撮像装置に関し、特に光学系の開口数NAから決定されるカットオフ空間周波数より高い周波数成分を除去する空間フィルタ処理を行うようにした撮像方法及びその撮像装置に関する。   The present invention relates to an imaging method suitable for high-magnification imaging such as a microscope and an imaging apparatus thereof, and in particular, to perform spatial filter processing for removing frequency components higher than a cutoff spatial frequency determined from the numerical aperture NA of an optical system. The present invention relates to an imaging method and an imaging apparatus thereof.

従来、撮像装置として、撮像光学系を通じて入射される光をCCD(charge−Coupled Device)で光電変換し出力するものが知られている。例えば特開2004−289786号公報(特許文献1)には、図10に示すように、撮像部101及び画像処理部110を具備した撮像装置100が開示されている。撮像部101は、撮像対象からの入射光を受光して光電変換する撮像手段として機能するものであり、例えばCCDカメラが用いられる。撮像部101にはCCDが撮像素子として用いられ、そのCCDの前方には撮像光学系が設置されている。   2. Description of the Related Art Conventionally, as an imaging device, a device that photoelectrically converts light incident through an imaging optical system using a CCD (charge-coupled device) and outputs the same is known. For example, Japanese Unexamined Patent Application Publication No. 2004-289786 (Patent Document 1) discloses an imaging apparatus 100 including an imaging unit 101 and an image processing unit 110 as shown in FIG. The imaging unit 101 functions as an imaging unit that receives incident light from an imaging target and performs photoelectric conversion, and for example, a CCD camera is used. The imaging unit 101 uses a CCD as an imaging element, and an imaging optical system is installed in front of the CCD.

画像処理部110は、撮像部101から出力される撮影画像を画像処理する画像処理手段としての機能を有しており、空間フィルタ111及びD/A変換部112で構成されている。空間フィルタ111は撮影信号の空間周波数における所定の低周波成分を除去し、撮影映像における高輝度部の周辺のハレーションを低減する空間フィルタである。D/A変換部112は、空間フィルタ111の出力信号をデジタル信号からアナログ信号に変換し、映像信号として出力する。   The image processing unit 110 has a function as an image processing unit that performs image processing on a captured image output from the imaging unit 101, and includes a spatial filter 111 and a D / A conversion unit 112. The spatial filter 111 is a spatial filter that removes a predetermined low-frequency component in the spatial frequency of the photographic signal and reduces halation around the high-luminance portion in the photographic image. The D / A conversion unit 112 converts the output signal of the spatial filter 111 from a digital signal to an analog signal and outputs it as a video signal.

空間フィルタ111としては、離散フーリエ変換処理、低周波除去処理及び逆離散フーリエ変換処理を行うもの、或いは1次元デジタルフィルタ又は2次元デジタルフィルタを用いて撮影画像の空間周波数における低周波成分を除去するものが用いられる。   The spatial filter 111 performs a discrete Fourier transform process, a low-frequency removal process, and an inverse discrete Fourier transform process, or removes a low-frequency component in the spatial frequency of a captured image using a one-dimensional digital filter or a two-dimensional digital filter. Things are used.

上述の特許文献1に記載の撮像装置は、撮影映像における高輝度部の周辺のハレーションを低減するために、撮影信号の空間周波数における所定の低周波成分を除去するものであるが、撮影画像をデジタル変換する際に混入する量子化誤差を低減するために、撮影信号の高周波成分を除去する撮像装置も知られている。
特開2004−289786号公報
The imaging device described in Patent Document 1 described above removes a predetermined low-frequency component in the spatial frequency of a photographic signal in order to reduce halation around a high-luminance portion in a photographic image. An imaging apparatus that removes a high-frequency component of a photographic signal is also known in order to reduce a quantization error mixed during digital conversion.
JP 2004-289786 A

上述の特許文献1に記載の撮影画像の低周波成分を除去するようにした撮像装置や、量子化誤差を低減するために撮影画像の高周波成分を除去するようにした撮像装置では、除去する際のカットオフ空間周波数の決定は、経験や勘に基づいて行われることが多い。このようにカットオフ空間周波数を経験や勘に基づいて決定していたのでは、量子化誤差を効果的に低減することができず、適切な撮影ができないという問題がある。   In the image pickup apparatus that removes the low-frequency component of the photographed image described in Patent Document 1 described above or the image pickup apparatus that removes the high-frequency component of the photographed image in order to reduce the quantization error, The cut-off spatial frequency is often determined based on experience and intuition. Thus, if the cut-off spatial frequency is determined based on experience and intuition, the quantization error cannot be effectively reduced, and there is a problem that appropriate photographing cannot be performed.

本発明は上述のような事情によりなされたものであり、本発明の目的は、撮影画像の低周波成分を除去する際のカットオフ空間周波数を、量子化誤差等の高周波ノイズを効果的かつ確実に低減できるように設定する空間フィルタ処理を行う撮像方法及びその撮像装置を提供することにある。   The present invention has been made for the above-described circumstances, and an object of the present invention is to effectively and reliably reduce a cutoff spatial frequency when removing a low-frequency component of a captured image and a high-frequency noise such as a quantization error. It is an object of the present invention to provide an imaging method and an imaging apparatus for performing a spatial filter process that is set so as to be reduced.

本発明は、光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法に関し、本発明の上記目的は、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去して空間フィルタ処理を行うことによって達成される。   The present invention relates to an imaging method for receiving and imaging light from an optical system and performing image processing on image data as an imaging output. The object of the present invention is a cutoff spatial frequency determined by the numerical aperture NA of the optical system. This is achieved by removing higher frequency components from the image data and performing spatial filtering.

本発明は、光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法に関し、本発明の上記目的は、前記画像データを離散フーリエ変換して第1画像データを生成し、前記第1画像データから前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を除去して第2画像データを生成し、前記第2画像データを逆離散フーリエ変換して第3画像データを生成して出力することによって空間フィルタ処理を行うことにより、或いは前記画像データに1次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去することによって空間フィルタ処理を行い、或いは前記画像データに2次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去することによって空間フィルタ処理を行い、或いは前記撮像を複数の撮像素子により行い、前記複数の撮像素子は互いの距離が、前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように配置され、前記画像データから前記カットオフ空間周波数より高い周波数成分を除去して空間フィルタ処理を行うことにより、或いは撮像素子を走査することにより前記画像データを生成し、前記撮像素子が受光する点間の距離が前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように、前記撮像素子の走査速度及び前記撮像素子の受光の時間間隔を設定し、前記画像データが前記カットオフ空間周波数より高い周波数成分を含まないように空間フィルタ処理を行うことによって達成される。   The present invention relates to an imaging method for receiving and imaging light from an optical system and performing image processing on image data which is an imaging output. The object of the present invention is to perform first Fourier transform on the image data by performing a discrete Fourier transform. The second image data is generated by removing a frequency component higher than the cutoff spatial frequency determined by the numerical aperture NA of the optical system from the first image data, and the second image data is subjected to inverse discrete Fourier transform. Generating a third image data and outputting it, or performing a spatial filtering process, or performing a one-dimensional digital filtering process on the image data, and a frequency higher than a cutoff spatial frequency determined by the numerical aperture NA of the optical system Spatial filtering is performed by removing components from the image data, or two-dimensional digital filtering is applied to the image data. And performing spatial filtering by removing a frequency component higher than a cutoff spatial frequency determined by the numerical aperture NA of the optical system from the image data, or performing the imaging with a plurality of imaging elements, and performing the plurality of imaging The elements are arranged such that the mutual distance is twice the cutoff spatial frequency determined by the numerical aperture NA of the optical system, and spatial filtering is performed by removing frequency components higher than the cutoff spatial frequency from the image data. Or by scanning the image sensor so that the distance between the points received by the image sensor is twice the cut-off spatial frequency determined by the numerical aperture NA of the optical system. In addition, a scanning speed of the image sensor and a time interval of light reception of the image sensor are set, and the image data is It is achieved by performing spatial filtering to be free of frequency components higher than-off spatial frequency.

また、本発明は、光学系からの光を受光して撮像し、撮像出力である第1画像データを画像処理する撮像方法に関し、本発明の上記目的は、撮像素子を走査することにより前記第1画像データを生成し、前記光学系の開口数NAに基づいてカットオフ空間周波数を決定し、前記撮像素子が受光する点間の距離が前記カットオフ空間周波数の2倍となるデータを前記第1画像データから間引いて第2画像データを生成し、前記第2画像データが前記カットオフ周波数より高い周波数成分を含まないようにして空間フィルタ処理を行うことによって達成され、更に前記開口数NAを外部から入力することによって、或いは前記画像データ若しくは前記第1画像データの光強度がピーク値から略半分になる位置までをスポット像とし、前記スポット像から前記スポット像の直径dを算出し、前記開口数NAを
NA=cλ/πd
(ただし、cは(2J(x))=0.5(J(x)は第1種ベッセル関 数の最小の解であり、約1.61634の定数、λは光の波長である。)
によって求めることによって、或いは前記画像データの若しくは前記第1画像データの前記スポット像の第1暗環をエッジ抽出法等で抽出し、前記第1暗環の直径dを求め、前記開口数NAを
NA=cλ/πd
(ただし、cは(2J(x))=0(J(x)は第1種ベッセル関数) の最小の解であり、約3.832の定数、λは光の波長である。)
によって求めることによって、より効果的に達成される。
The present invention also relates to an imaging method for receiving and imaging light from an optical system and performing image processing on first image data that is an imaging output. The object of the present invention is to perform the first by scanning an imaging device. One image data is generated, a cutoff spatial frequency is determined based on the numerical aperture NA of the optical system, and data in which a distance between points received by the image sensor is twice the cutoff spatial frequency is the first data. The second image data is generated by decimating from one image data, and the second image data does not contain a frequency component higher than the cutoff frequency, and is subjected to a spatial filter process, and the numerical aperture NA is further reduced. The spot image is formed by inputting from the outside, or a position where the light intensity of the image data or the first image data is approximately half of the peak value, and the spot image Calculating the diameter d of the spot image from the numerical aperture NA
NA = cλ / πd
(Where c is (2J 1 (x)) 2 = 0.5 (J 1 (x) is the smallest solution of the first type Bessel function, a constant of about 1.61634, and λ is the wavelength of light. is there.)
Or the first dark ring of the spot image of the image data or the first image data is extracted by an edge extraction method or the like, the diameter d of the first dark ring is obtained, and the numerical aperture NA is determined.
NA = cλ / πd
Where c is the minimum solution of (2J 1 (x)) 2 = 0 (J 1 (x) is the first kind Bessel function), a constant of about 3.832, and λ is the wavelength of light. )
Is achieved more effectively.

更に、本発明は、光学系からの光を受光して撮像する撮像手段と、前記撮像手段から出力される画像データを画像処理する画像処理手段とを具備した撮像装置に関し、本発明の上記目的は、前記画像処理手段が、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去する空間フィルタを具備することによって、或いは前記画像処理手段が、前記画像データを離散フーリエ変換して第1画像データを生成するフーリエ変換部と、前記光学系の開口数NAに基づいて定まるカットオフ空間周波数より高い周波数成分を前記第1画像データから除去して第2画像データを生成する空間フィルタと、前記第2画像データを逆離散フーリエ変換して第3画像データを生成して出力する逆フーリエ変換部とを具備することによって、或いは前記画像処理手段が前記画像データに1次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去する機能を具備することによって、或いは前記画像処理手段が前記画像データに2次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去する機能を具備することによって、或いは前記撮像手段は複数の撮像素子から成り、前記複数の撮像素子は互いの距離が前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように配置され、前記カットオフ空間周波数より高い周波数成分が前記画像データから除去されることによって、或いは前記撮像手段は撮像素子で成り、前記撮像素子を走査することにより前記画像データを生成し、前記撮像素子が受光する点間の距離が、前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように、前記撮像素子の走査速度及び前記撮像素子の受光の時間間隔を設定して前記走査を行い、前記画像データが前記カットオフ空間周波数より高い周波数成分を含まないようにすることによって、或いは前記撮像手段が、撮像素子を走査することにより前記第1画像データを生成して出力する手段で成り、前記画像処理手段は、前記光学系の開口数NAに基づいてカットオフ空間周波数を決定する機能と、前記撮像素子が受光を行う点間の距離が前記カットオフ空間周波数の2倍となるデータを前記第1画像データから間引いて第2画像データを生成する機能とを具備し、前記第2画像データが前記カットオフ空間周波数より高い周波数成分を含まないようにすることによって達成される。   Furthermore, the present invention relates to an imaging apparatus including an imaging unit that receives and captures light from an optical system, and an image processing unit that performs image processing on image data output from the imaging unit. The image processing means comprises a spatial filter that removes from the image data a frequency component higher than a cutoff spatial frequency determined by the numerical aperture NA of the optical system, or the image processing means comprises the image data A Fourier transform unit for generating a first image data by performing a discrete Fourier transform on the second image by removing a frequency component higher than a cutoff spatial frequency determined based on a numerical aperture NA of the optical system from the first image data. A spatial filter for generating data, and inverse Fourier transform for generating and outputting third image data by performing inverse discrete Fourier transform on the second image data. Or the image processing means performs one-dimensional digital filter processing on the image data, and removes a frequency component higher than a cutoff spatial frequency determined by the numerical aperture NA of the optical system from the image data. A function of removing a frequency component higher than a cutoff spatial frequency determined by a numerical aperture NA of the optical system from the image data by providing a function or by the image processing means performing two-dimensional digital filter processing on the image data Or the imaging means comprises a plurality of imaging devices, and the plurality of imaging devices are arranged such that the distance between them is twice the cutoff spatial frequency determined by the numerical aperture NA of the optical system. A frequency component higher than the cut-off spatial frequency is removed from the image data. Or the imaging means comprises an imaging device, and the image data is generated by scanning the imaging device, and the distance between points received by the imaging device is determined by the numerical aperture NA of the optical system. The scanning is performed by setting the scanning speed of the image sensor and the light reception time interval of the image sensor so as to be twice the cut-off spatial frequency, and the image data has a frequency component higher than the cut-off spatial frequency. The image processing means includes means for generating and outputting the first image data by scanning the image sensor, or the image processing means is configured to adjust the numerical aperture NA of the optical system. A function for determining a cut-off spatial frequency based on the data, and data in which a distance between points at which the image sensor receives light is twice the cut-off spatial frequency. And the second image data is generated by thinning out from one image data, and the second image data does not include a frequency component higher than the cut-off spatial frequency.

本発明によれば、撮影画像から高周波の成分を除去する際のカットオフ空間周波数を光学系の開口数NAに基づいて客観的な数式で決定しているので、光学系を通過する光の情報は維持し、下流で混入する種々の高周波ノイズを効果的かつ確実に低減することができ、適切な撮影画像を生成することができる。   According to the present invention, the cut-off spatial frequency for removing high-frequency components from the photographed image is determined by an objective mathematical expression based on the numerical aperture NA of the optical system, so that information on the light passing through the optical system is obtained. Can be maintained, various high-frequency noises mixed in downstream can be effectively and reliably reduced, and an appropriate photographed image can be generated.

本発明は、CCD等の撮像素子に結像した光には光学系の開口数NAから決められるカットオフ空間周波数より高い空間周波数は含まれていないので、カットオフ空間周波数より高い空間周波数成分の信号はノイズとみなすことができることに着眼している。
本発明では、量子化誤差等の高周波ノイズを効果的かつ確実に低減できるように、光学系からの光を受光して撮像した画像データの画像処理において、光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を画像データから除去して空間フィルタ処理を行うようにしている。空間フィルタ処理は1次元又は2次元デジタルフィルタ処理、撮像素子の走査速度や受光の時間間隔の設定で行っている。
In the present invention, since the light imaged on the image sensor such as a CCD does not include a spatial frequency higher than the cutoff spatial frequency determined from the numerical aperture NA of the optical system, a spatial frequency component higher than the cutoff spatial frequency is included. We focus on the fact that the signal can be regarded as noise.
In the present invention, in order to effectively and reliably reduce high-frequency noise such as quantization error, in image processing of image data picked up by receiving light from the optical system, a cutoff determined by the numerical aperture NA of the optical system Spatial filtering is performed by removing frequency components higher than the spatial frequency from the image data. Spatial filter processing is performed by setting one-dimensional or two-dimensional digital filter processing, scanning speed of the image sensor, and time interval of light reception.

以下、本発明の実施形態について図面を参照して説明する。   Embodiments of the present invention will be described below with reference to the drawings.

本発明に係る撮像装置は、図1に示されるように光学系1からの入射光を受光する撮像部2と、撮像部2から出力される画像データIM1を処理する画像処理部3とから構成されており、画像処理部3から出力される画像データIM2は、光学系1の開口数NAによって定まるカットオフ空間周波数より高い周波数成分が除去されている。光学系1は、例えば顕微鏡等の高倍率の像を形成するための光学系であり、1枚又は複数枚のレンズの組み合わせから成り、撮像対象及び光学系1を透過した光が撮像部2に入射する。撮像部2には、光学系からの入射光を受光し光電変換するための撮像手段としての撮像素子が備えられており、撮像素子としては例えば1次元又は2次元のCCDカメラが用いられる。   As shown in FIG. 1, the imaging apparatus according to the present invention includes an imaging unit 2 that receives incident light from the optical system 1 and an image processing unit 3 that processes image data IM1 output from the imaging unit 2. Thus, the image data IM2 output from the image processing unit 3 has frequency components higher than the cutoff spatial frequency determined by the numerical aperture NA of the optical system 1 removed. The optical system 1 is an optical system for forming a high-magnification image such as a microscope, for example. Incident. The imaging unit 2 is provided with an imaging device as an imaging means for receiving and photoelectrically converting incident light from the optical system. As the imaging device, for example, a one-dimensional or two-dimensional CCD camera is used.

画像処理部3には空間フィルタ部31が備えられ、撮像部2から出力された画像データIM1を入力し、画像処理されカットオフ空間周波数より高い周波数成分が除去された画像データIM2を出力するようになっている。空間フィルタ部31は、入力された画像データIM1のカットオフ空間周波数より高い周波数成分を除去するもので、量子化誤差等の高周波ノイズを効果的に低減するために、本発明ではカットオフ空間周波数は光学系1の開口数NAに基づいて定まる最大空間周波数を用いる。即ち、光の波長をλとすると、最大空間周波数fは下記数1によって与えられる。
(数1)
=2NA/λ
図2に示されるように、空間フィルタ部31からの出力信号IMS、つまり画像データIM1からカットオフ空間周波数(f)より高い周波数成分を除去されたデジタル信号をアナログ信号に変換するために、画像処理部3にD/A変換部32を具備するようにしても良い。なお、図2に示される実施形態では、デジタル信号の画像データIM1を空間フィルタ部31で処理してからD/A変換部32でアナログ信号に変換するようにしているが、画像データIM1がアナログ信号であり、空間フィルタ部31での空間フィルタ処理をデジタル信号によって行う場合、A/D変換部を空間フィルタ部31の前に挿入し画像データIM1をデジタル信号に変換してから空間フィルタ部31で画像データを処理するようにしても良い。また、D/A変換部及びA/D変換部は本発明に必須のものではなく、撮像部2からの画像データIM1の形態或いは画像処理部3の出力である画像データIM2の利用形態に応じて適宜設ければ良い。
The image processing unit 3 includes a spatial filter unit 31 that receives the image data IM1 output from the imaging unit 2, and outputs the image data IM2 that has been subjected to image processing and from which frequency components higher than the cutoff spatial frequency have been removed. It has become. The spatial filter unit 31 removes a frequency component higher than the cutoff spatial frequency of the input image data IM1, and in order to effectively reduce high frequency noise such as a quantization error, the spatial filter unit 31 according to the present invention has a cutoff spatial frequency. Uses the maximum spatial frequency determined based on the numerical aperture NA of the optical system 1. That is, assuming that the wavelength of light is λ, the maximum spatial frequency f m is given by the following equation (1).
(Equation 1)
f m = 2NA / λ
As shown in FIG. 2, in order to convert an output signal IMS from the spatial filter unit 31, that is, a digital signal in which a frequency component higher than the cutoff spatial frequency (f m ) is removed from the image data IM1, into an analog signal, The image processing unit 3 may include a D / A conversion unit 32. In the embodiment shown in FIG. 2, the digital signal image data IM1 is processed by the spatial filter unit 31 and then converted to an analog signal by the D / A conversion unit 32. However, the image data IM1 is analog. When the spatial filter processing in the spatial filter unit 31 is performed with a digital signal, the A / D conversion unit is inserted before the spatial filter unit 31 to convert the image data IM1 into a digital signal, and then the spatial filter unit 31 In this case, the image data may be processed. Further, the D / A conversion unit and the A / D conversion unit are not essential to the present invention, depending on the form of the image data IM1 from the imaging unit 2 or the usage form of the image data IM2 that is the output of the image processing unit 3. May be provided as appropriate.

本発明の第1実施形態では、空間フィルタ部31は図3に示すように、撮像部2からの画像データIM1を離散フーリエ変換するフーリエ変換部311と、離散フーリエ変換された画像データIMS1からカットオフ空間周波数より高い周波数成分の除去を行う空間フィルタ312と、空間フィルタ312で高周波数成分を除去された画像データIMS2を逆離散フーリエ変換して画像データIMSを出力する逆フーリエ変換部313とで構成されている。   In the first embodiment of the present invention, as shown in FIG. 3, the spatial filter unit 31 includes a Fourier transform unit 311 that performs discrete Fourier transform on the image data IM1 from the imaging unit 2, and a cut from the image data IMS1 that has been subjected to discrete Fourier transform. A spatial filter 312 that removes frequency components higher than the off-spatial frequency, and an inverse Fourier transform unit 313 that performs inverse discrete Fourier transform on the image data IMS2 from which the high frequency component has been removed by the spatial filter 312 and outputs the image data IMS. It is configured.

図4(A)〜(D)はこのような空間フィルタ処理の過程を示しており、各過程における画像が示されている。先ず、図4(A)に示す撮像部2からの原画像10がフーリエ変換部331で離散フーリエ変換され、図4(B)に示す画像11となり、原画像10におけるスポット10Aは、離散フーリエ変換された画像11では4隅に分離され、図4(B)に示す符号部11Aになる。そして、空間フィルタ312が画像11に対して高周波成分の除去を行うと、図4(C)に示すような画像12が得られ、図4(B)の符号11Bに現れているノイズが低減され、図4(C)の符号12Bに示されるようになる。この画像12をフーリエ変換部313で逆離散フーリエ変換すると、高周波成分のノイズが除去された 図4(D)に示すような画像13が得られる。ここで、除去される高周波成分のカットオフ周波数は、前記数1で与えられる最大空間周波数fである。 4A to 4D show the process of such spatial filter processing, and images in each process are shown. First, the original image 10 from the imaging unit 2 shown in FIG. 4A is subjected to discrete Fourier transform by the Fourier transform unit 331 to become an image 11 shown in FIG. 4B, and the spot 10A in the original image 10 is a discrete Fourier transform. The image 11 is separated into four corners and becomes the encoding unit 11A shown in FIG. Then, when the spatial filter 312 removes high frequency components from the image 11, an image 12 as shown in FIG. 4C is obtained, and noise appearing at reference numeral 11B in FIG. 4B is reduced. , As indicated by reference numeral 12B in FIG. When this image 12 is subjected to inverse discrete Fourier transform by the Fourier transform unit 313, an image 13 as shown in FIG. 4D from which high-frequency component noise has been removed is obtained. Here, the cutoff frequency of the high frequency components removed is the maximum spatial frequency f m given by the number 1.

一方、空間フィルタ部31において、1次元デジタルフィルタ処理又は2次元デジタルフィルタ処理を行うようにしても良い(第2実施形態)。   On the other hand, the spatial filter unit 31 may perform one-dimensional digital filter processing or two-dimensional digital filter processing (second embodiment).

1次元デジタルフィルタ処理では例えば図5に示されるような3画素のオペレータ(マスク)を用い、2次元デジタルフィルタ処理では図6に示されるような9画素のオペレータ(マスク)を用いる。共に画像データの平滑化を行い、前記数1で与えられる最大空間周波数fより高い周波数成分を除去することができる。図5及び図6のフィルタのタップ数は“3”であるが、このフィルタのタップ数はできるだけ大きくすることが好ましい。 In the one-dimensional digital filter processing, for example, a three-pixel operator (mask) as shown in FIG. 5 is used, and in the two-dimensional digital filter processing, a nine-pixel operator (mask) as shown in FIG. 6 is used. Both perform smoothing of the image data, the frequency components higher than the maximum spatial frequency f m given by the number 1 may be removed. The number of taps of the filter of FIGS. 5 and 6 is “3”, but it is preferable to increase the number of taps of this filter as much as possible.

図7に示される第3実施形態では、撮像部2を構成する多数の撮像素子2Aを複数箇所に配設し、撮像素子2A間の距離を前記数1で与えられる最大空間周波数fの2倍程度にすることにより、撮像素子2Aで生成される画像データIM3に対して最大空間周波数fより高い周波数成分が除去されるようにしている。撮像素子2Aは2次元のエリア素子であっても、線状のリニア素子であっても良く、リニア素子の場合には空間的若しくは時間的に走査する必要がある。また、撮像素子2AはCCDであっても良い。 In the third embodiment shown in FIG. 7, a large number of imaging elements 2A constituting the imaging unit 2 are arranged at a plurality of locations, and the distance between the imaging elements 2A is 2 at the maximum spatial frequency f m given by the above equation 1. by approximately doubled, frequency components higher than the maximum spatial frequency f m is to be removed with respect to the image data IM3 generated by an imaging element 2A. The imaging element 2A may be a two-dimensional area element or a linear linear element. In the case of a linear element, it is necessary to scan spatially or temporally. The image sensor 2A may be a CCD.

撮像素子2Aは、光学系1からの入射光を受光するために2次元平面内に配設され、互いの距離が最大空間周波数fの2倍になるように配設される。このようにすることによって、画像処理部3に送られる画像データに、最大空間周波数fより高い周波数成分が含まれないようにすることができる。 Imaging element 2A is disposed in a two-dimensional plane for receiving incident light from the optical system 1 is arranged such that the distance therebetween is 2 times the maximum spatial frequency f m. By doing so, the image data sent to the image processing unit 3, can be prevented contain frequency components higher than the maximum spatial frequency f m.

また、図8に示される第4実施形態では、単一の撮像素子2Aが、光学系1からの入射光を受光するために2次元平面内を走査することにより、画像処理部3に画像データIM4を送るようにしている。2次元平面内の走査は、線状の走査を少しずつ空間的若しくは時間的にずらして繰り返し行うことによって実現できる。   In the fourth embodiment shown in FIG. 8, the single image sensor 2A scans the two-dimensional plane in order to receive the incident light from the optical system 1, whereby the image processing unit 3 receives image data. IM4 is sent. Scanning in a two-dimensional plane can be realized by repeatedly performing linear scanning with a slight shift in space or time.

本実施形態では、撮像素子2Aが光学系1からの入射光を受光する点間の距離が前記数1で与えられる最大空間周波数fの2倍になるように、撮像素子2Aの走査速度と撮像素子2Aが受光を行う時間間隔(サンプリング間隔)とを設定する。ナイキスト・シャノンのサンプリング定理にあるように、このようにすることによって撮像素子2Aによって生成される画像データIM4が、最大空間周波数fより高い周波数成分を含まないようにすることができる。 In this embodiment, as the image pickup device 2A is 2 times the maximum spatial frequency f m of distance between a point for receiving incident light is given by the number 1 from the optical system 1, and the scanning speed of the image pickup device 2A A time interval (sampling interval) at which the image sensor 2A receives light is set. As in the Nyquist Shannon sampling theorem, the image data IM4 generated by the imaging element 2A by doing so it can be made to not contain frequency components higher than the maximum spatial frequency f m.

また、撮像素子2Aの走査速度と撮像素子2Aが受光を行う時間間隔を任意に設定した場合でも、画像処理部3において、受光された点間の距離が最大空間周波数fの2倍になるデータのみを間引くようにして、間引かれたデータのみから画像データIM2を生成することによって、画像データIM2が最大空間周波数fより高い周波数成分を含まないようにしても良い。 Further, even when the scanning speed and the imaging device 2A of the imaging element 2A was arbitrarily set the time interval of light receiving, in the image processing unit 3, the distance between the light receiving been points is 2 times the maximum spatial frequency f m The image data IM2 may not include a frequency component higher than the maximum spatial frequency f m by generating only the data and generating the image data IM2 from only the thinned data.

光学系1の開口数NAは撮像部2に予め設定しておいても良いが、開口数NAは光学系1の絞り等の条件によって変化するため、計算して求めた開口数NAを外部からその都度入力するようにしても良い。また、光学系の開口数NAが完全に既知でない場合には、想定できる開口数NAのうち、最大の開口数NAを用いてカットオフ空間周波数を決めても良い。例えば観察対象の光の開口数NA1が未知であり、顕微鏡などの測定装置の光学系の結像部における最大実効開口数NA2が既知である場合には、開口数NA1が十分に大きい場合について、結像部における実効開口数NA0を算出し、実効開口数NA0を用いてカットオフ空間周波数を決定する。   The numerical aperture NA of the optical system 1 may be set in the imaging unit 2 in advance. However, since the numerical aperture NA varies depending on conditions such as the diaphragm of the optical system 1, the numerical aperture NA obtained by calculation is externally applied. You may make it input each time. When the numerical aperture NA of the optical system is not completely known, the cutoff spatial frequency may be determined using the maximum numerical aperture NA among the possible numerical apertures NA. For example, when the numerical aperture NA1 of the light to be observed is unknown and the maximum effective numerical aperture NA2 in the imaging unit of the optical system of the measuring device such as a microscope is known, the numerical aperture NA1 is sufficiently large. An effective numerical aperture NA0 in the imaging unit is calculated, and a cutoff spatial frequency is determined using the effective numerical aperture NA0.

更に、撮像部2にスポット像解析部を具備するようにして、スポット像から開口数NAを計算するようにすることによって、光学系1の所望の条件における開口数NAを撮像部2に設定することができる。   Furthermore, the numerical aperture NA under the desired conditions of the optical system 1 is set in the imaging unit 2 by calculating the numerical aperture NA from the spot image by providing the imaging unit 2 with a spot image analysis unit. be able to.

光学系1の開口数NAは、スポット像から求めることができる。即ち、図9に示されるように、光学系1から像面6に結像するスポット像に対して、入射光の入射角をθ、媒質の屈折率をnとすると、開口数NAは下記数2によって定義される。
(数2)
NA=n×sinθ

スポット像解析部は、例えばスポット画像の光強度がピーク値から略半分になる位置までをスポット像とし、スポット像の面積sを求める。そして、スポット像の面積sから、スポット像の直径dを下記数3により算出する。
(数3)
d=√(4s/π)

上記数3によってスポット像の直径dが求まると、下記数4により開口数NAを求めることができる。
(数4)
NA=cλ/πd
ただし、cは(2J(x))=0.5(J(x)は第1種ベッセル関数 )の最小の解であり、約1.61634の定数、λは光の波長である。
なお、スポット像の直径dは、スポット像の面積sから求めるのではなく、スポット像の直径から直接求めるようにしても良い。
The numerical aperture NA of the optical system 1 can be obtained from the spot image. That is, as shown in FIG. 9, when the incident angle of incident light is θ and the refractive index of the medium is n with respect to the spot image formed on the image plane 6 from the optical system 1, the numerical aperture NA is the following number. Defined by 2.
(Equation 2)
NA = n × sin θ

For example, the spot image analysis unit obtains the area s of the spot image by setting the spot image to a position where the light intensity of the spot image is approximately half from the peak value. Then, the diameter d of the spot image is calculated from the area s of the spot image by the following formula 3.
(Equation 3)
d = √ (4 s / π)

When the diameter d of the spot image is obtained by the above equation 3, the numerical aperture NA can be obtained by the following equation 4.
(Equation 4)
NA = cλ / πd
Where c is the minimum solution of (2J 1 (x)) 2 = 0.5 (J 1 (x) is the first kind Bessel function), a constant of about 1.61634, and λ is the wavelength of light .
The diameter d of the spot image may be obtained directly from the diameter of the spot image rather than from the area s of the spot image.

また、スポット像解析部において、スポット画像の第1暗環をエッジ抽出法等で抽出し、抽出された第1暗環の直径dを求める方法もある。この場合、直径dを求めた後、下記数5によって開口数NAを求めることができる。
(数5)
NA=cλ/πd
ただし、cは(2J(x))=0(J(x)は第1種ベッセル関数)の 最小の解であり、約3.832の定数である。

なお、開口数NAに誤差が含まれている場合を想定する場合には、安全係数βを乗算してカットオフ空間周波数fを下記数6で決定する。
(数6)
=2β・NA/λ

例えば開口数NAに20%の誤差が想定される場合には、安全係数βを“1.2”とする。
Also, there is a method in which the spot image analysis unit extracts the first dark ring of the spot image by an edge extraction method or the like and obtains the diameter d of the extracted first dark ring. In this case, after obtaining the diameter d, the numerical aperture NA can be obtained by the following equation (5).
(Equation 5)
NA = cλ / πd
However, c is the minimum solution of (2J 1 (x)) 2 = 0 (J 1 (x) is the first kind Bessel function), and is a constant of about 3.832.

In the case where it is assumed that an error is included in the numerical aperture NA determines the cut-off spatial frequency f m by the following Equation 6 is multiplied by a safety factor beta.
(Equation 6)
f m = 2β · NA / λ

For example, when an error of 20% is assumed in the numerical aperture NA, the safety coefficient β is set to “1.2”.

以上、本発明の実施形態について具体的に説明したが、本発明はこれら実施形態に限定されるものではなく、その趣旨を逸脱しない範囲で適宜変更可能である。特に光学系の瞳については、本発明の実施形態は、任意の瞳形状に適用することができる。例えば、瞳の開口部で“1”、それ以外の領域では“0”を持つ瞳のマスク関数f(x,y)の自己相関関数の台(関数が値を持つ領域)の空間周波数信号を通過させ、台以外の帯域の信号をカットするローパスフィルタを構成すれば、任意の瞳形状に対して本発明を適用することができる。   As mentioned above, although embodiment of this invention was described concretely, this invention is not limited to these embodiment, In the range which does not deviate from the meaning, it can change suitably. Particularly for the pupil of the optical system, the embodiments of the present invention can be applied to any pupil shape. For example, the spatial frequency signal of the autocorrelation function platform (region in which the function has a value) of the pupil mask function f (x, y) having “1” in the pupil opening and “0” in the other regions is obtained. The present invention can be applied to an arbitrary pupil shape by configuring a low-pass filter that passes and cuts signals in bands other than the table.

本発明に係る撮像装置の構成例を示すブロック図である。It is a block diagram which shows the structural example of the imaging device which concerns on this invention. 本発明に係る撮像装置の画像処理部にD/A変換部を具備するようにした場合の構成例を示すブロック図である。It is a block diagram which shows the structural example at the time of making it comprise a D / A conversion part in the image processing part of the imaging device which concerns on this invention. 本発明の空間フィルタ部の構成例を示すブロック図である。It is a block diagram which shows the structural example of the spatial filter part of this invention. 本発明の画像処理を説明するための図である(第1実施形態)。It is a figure for demonstrating the image processing of this invention (1st Embodiment). 本発明で実行する1次元デジタルフィルタ処理のオペレータの例を示す図である(第2実施形態)。It is a figure which shows the example of the operator of the one-dimensional digital filter process performed by this invention (2nd Embodiment). 本発明で実行する2次元デジタルフィルタ処理のオペレータの例を示す図である(第2実施形態)。It is a figure which shows the example of the operator of the two-dimensional digital filter process performed by this invention (2nd Embodiment). 複数の撮像素子を最大空間周波数の2倍の間隔で配設した撮像装置(第3実施形態)の構成例を示すブロック図である。It is a block diagram which shows the structural example of the imaging device (3rd Embodiment) which has arrange | positioned several image pick-up element by the space | interval twice as large as the maximum spatial frequency. 撮像素子を走査させることによって画像データを生成する撮像装置(第4実施形態)の構成例を示すブロック図である。It is a block diagram which shows the structural example of the imaging device (4th Embodiment) which produces | generates image data by scanning an imaging device. 光学系の開口数のスポット像による定義を説明するための図である。It is a figure for demonstrating the definition by the spot image of the numerical aperture of an optical system. 従来の撮像装置の構成例を示すブロック図である。It is a block diagram which shows the structural example of the conventional imaging device.

符号の説明Explanation of symbols

1 光学系
2 撮像部
2A 撮像素子
3 画像処理部
31 空間フィルタ部
32 D/A変換部
311 フーリエ変換部
312 空間フィルタ
313 逆フーリエ変換部
6 像面
DESCRIPTION OF SYMBOLS 1 Optical system 2 Image pick-up part 2A Image pick-up element 3 Image processing part 31 Spatial filter part 32 D / A conversion part 311 Fourier transform part 312 Spatial filter 313 Inverse Fourier transform part 6 Image surface

Claims (20)

光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法において、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去することを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged, and image data as an imaging output is processed, a frequency component higher than a cutoff spatial frequency determined by the numerical aperture NA of the optical system is removed from the image data. An imaging method for performing spatial filter processing characterized by the above. 光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法において、前記画像データを離散フーリエ変換して第1画像データを生成し、前記第1画像データから前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を除去して第2画像データを生成し、前記第2画像データを逆離散フーリエ変換して第3画像データを生成して出力することを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged, and image data that is an imaging output is image-processed, first image data is generated by performing discrete Fourier transform on the image data, and the first image data is used to generate the first image data. The second image data is generated by removing frequency components higher than the cutoff spatial frequency determined by the numerical aperture NA of the optical system, and the second image data is subjected to inverse discrete Fourier transform to generate and output third image data. An imaging method for performing spatial filter processing characterized by the above. 光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法において、前記画像データに1次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去することを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged, and image data that is an imaging output is image-processed, a one-dimensional digital filter process is performed on the image data, and a cutoff space determined by a numerical aperture NA of the optical system An imaging method for performing spatial filtering, wherein a frequency component higher than a frequency is removed from the image data. 光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法において、前記画像データに2次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去することを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged, and image data that is an imaging output is image-processed, a two-dimensional digital filter process is performed on the image data, and a cut-off space determined by a numerical aperture NA of the optical system An imaging method for performing spatial filtering, wherein a frequency component higher than a frequency is removed from the image data. 光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法において、前記撮像を複数の撮像素子により行い、前記複数の撮像素子は互いの距離が、前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように配置され、前記画像データから前記カットオフ空間周波数より高い周波数成分を除去することを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged, and image data that is an imaging output is image-processed, the imaging is performed by a plurality of imaging elements, and the plurality of imaging elements are spaced apart from each other by the optical system. An imaging method for performing spatial filter processing, characterized in that a frequency component higher than the cut-off spatial frequency is removed from the image data, which is arranged so as to be twice the cut-off spatial frequency determined by the numerical aperture NA. 光学系からの光を受光して撮像し、撮像出力である画像データを画像処理する撮像方法において、撮像素子を走査することにより前記画像データを生成し、前記撮像素子が受光する点間の距離が前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように、前記撮像素子の走査速度及び前記撮像素子の受光の時間間隔を設定し、前記画像データが前記カットオフ空間周波数より高い周波数成分を含まないようにすることを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged, and image data that is an imaging output is image-processed, the image data is generated by scanning the image sensor, and the distance between points received by the image sensor Is set to be twice the cutoff spatial frequency determined by the numerical aperture NA of the optical system, the scanning speed of the imaging element and the time interval of light reception of the imaging element are set, and the image data is the cutoff spatial frequency An imaging method for performing spatial filtering, characterized by not including a higher frequency component. 光学系からの光を受光して撮像し、撮像出力である第1画像データを画像処理する撮像方法において、撮像素子を走査することにより前記第1画像データを生成し、前記光学系の開口数NAに基づいてカットオフ空間周波数を決定し、前記撮像素子が受光する点間の距離が前記カットオフ空間周波数の2倍となるデータを前記第1画像データから間引いて第2画像データを生成し、前記第2画像データが前記カットオフ周波数より高い周波数成分を含まないようにしたことを特徴とする空間フィルタ処理を行う撮像方法。 In an imaging method in which light from an optical system is received and imaged and image processing is performed on first image data that is an imaging output, the first image data is generated by scanning an image sensor, and the numerical aperture of the optical system Cut-off spatial frequency is determined based on NA, and second image data is generated by thinning out data from which the distance between the points received by the image sensor is twice the cut-off spatial frequency from the first image data. An imaging method for performing spatial filtering, wherein the second image data does not include a frequency component higher than the cutoff frequency. 前記開口数NAを外部から入力するようになっている請求項1乃至7のいずれかに記載の空間フィルタ処理を行う撮像方法。 The imaging method for performing spatial filter processing according to claim 1, wherein the numerical aperture NA is input from the outside. 前記画像データ若しくは前記第1画像データの光強度がピーク値から略半分になる位置までをスポット像とし、前記スポット像から前記スポット像の直径dを算出し、前記開口数NAを
NA=cλ/πd
(ただし、cは(2J(x))=0.5(J(x)は第1種ベッセル関 数)の最小の解であり、約1.61634の定数、λは光の波長である。)
として求める請求項1乃至7のいずれかに記載の空間フィルタ処理を行う撮像方法。
A spot image is defined as a spot image where the light intensity of the image data or the first image data is approximately half of the peak value, the diameter d of the spot image is calculated from the spot image, and the numerical aperture NA is calculated.
NA = cλ / πd
Where c is the minimum solution of (2J 1 (x)) 2 = 0.5 (J 1 (x) is the first type Bessel function), a constant of about 1.61634, and λ is the wavelength of light .)
The imaging method which performs the spatial filter process in any one of Claims 1 thru | or 7 calculated | required as.
前記画像データ若しくは前記第1画像データのスポット像の第1暗環をエッジ抽出法等で抽出し、前記第1暗環の直径dを求め、前記開口数NAを
NA=cλ/πd
(ただし、cは(2J(x))=0(J(x)は第1種ベッセル関数) の最小の解であり、約3.832の定数、λは光の波長である。)
として求める請求項1乃至7のいずれかに記載の空間フィルタ処理を行う撮像方法。
A first dark ring of the image data or the spot image of the first image data is extracted by an edge extraction method or the like to obtain a diameter d of the first dark ring, and the numerical aperture NA is calculated.
NA = cλ / πd
(However, c is the smallest solution of (2J 1 (x)) 2 = 0 (J 1 (x) is the first kind Bessel function), about 3.832 constants, lambda is the wavelength of light. )
The imaging method which performs the spatial filter process in any one of Claims 1 thru | or 7 calculated | required as.
光学系からの光を受光して撮像する撮像手段と、前記撮像手段から出力さ
れる画像データを画像処理する画像処理手段とを具備した撮像装置において、前記画像処
理手段が、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数
成分を前記画像データから除去する空間フィルタを具備していることを特徴とする空間フ
ィルタ処理を行う撮像装置。
An imaging apparatus comprising: an imaging unit that receives and captures light from an optical system; and an image processing unit that performs image processing on image data output from the imaging unit. The image processing unit includes an opening of the optical system. An imaging apparatus for performing spatial filter processing, comprising: a spatial filter that removes a frequency component higher than a cutoff spatial frequency determined by a number NA from the image data.
光学系からの光を受光して撮像する撮像手段と、前記撮像手段から出力される画像データを画像処理する画像処理手段とを具備した撮像装置において、前記画像処理手段が、前記画像データを離散フーリエ変換して第1画像データを生成するフーリエ変換部と、前記光学系の開口数NAに基づいて定まるカットオフ空間周波数より高い周波数成分を前記第1画像データから除去して第2画像データを生成する空間フィルタと、前記第2画像データを逆離散フーリエ変換して第3画像データを生成して出力する逆フーリエ変換部とを具備していることを特徴とする空間フィルタ処理を行う撮像装置。 An imaging apparatus comprising: an imaging unit that receives light from an optical system to capture an image; and an image processing unit that performs image processing on image data output from the imaging unit. The image processing unit discretely processes the image data. A Fourier transform unit that generates first image data by performing Fourier transform, and removes a frequency component higher than a cutoff spatial frequency determined based on the numerical aperture NA of the optical system from the first image data to obtain second image data. An imaging apparatus for performing spatial filter processing, comprising: a spatial filter to be generated; and an inverse Fourier transform unit configured to generate and output third image data by performing inverse discrete Fourier transform on the second image data . 光学系からの光を受光して撮像する撮像手段と、前記撮像手段から出力される画像データを画像処理する画像処理手段とを具備した撮像装置において、前記画像処理手段が前記画像データに1次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去する機能を具備していることを特徴とする空間フィルタ処理を行う撮像装置。 An image pickup apparatus including an image pickup unit that receives light from an optical system and picks up an image and an image processing unit that performs image processing on image data output from the image pickup unit. An imaging apparatus for performing spatial filter processing, comprising a function of performing digital filter processing and removing a frequency component higher than a cutoff spatial frequency determined by a numerical aperture NA of the optical system from the image data. 光学系からの光を受光して撮像する撮像手段と、前記撮像手段から出力される画像データを画像処理する画像処理手段とを具備した撮像装置において、前記画像処理手段が前記画像データに2次元デジタルフィルタ処理を行い、前記光学系の開口数NAによって定まるカットオフ空間周波数より高い周波数成分を前記画像データから除去する機能を具備していることを特徴とする空間フィルタ処理を行う撮像装置。 An imaging apparatus comprising: an imaging unit that receives light from an optical system to capture an image; and an image processing unit that performs image processing on image data output from the imaging unit. An imaging apparatus for performing spatial filter processing, comprising a function of performing digital filter processing and removing a frequency component higher than a cutoff spatial frequency determined by a numerical aperture NA of the optical system from the image data. 光学系からの光を受光して撮像する撮像手段と、前記撮像手段による画像データを画像処理する画像処理手段とを具備した撮像装置において、前記撮像手段は複数の撮像素子から成り、前記複数の撮像素子は互いの距離が前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように配置され、前記カットオフ空間周波数より高い周波数成分が前記画像データから除去されるようになっていることを特徴とする空間フィルタ処理を行う撮像装置。 An imaging apparatus comprising: an imaging unit that receives and captures light from an optical system; and an image processing unit that performs image processing on image data obtained by the imaging unit. The imaging unit includes a plurality of imaging elements, and The imaging elements are arranged such that the mutual distance is twice the cutoff spatial frequency determined by the numerical aperture NA of the optical system, and frequency components higher than the cutoff spatial frequency are removed from the image data. An imaging apparatus that performs spatial filter processing. 光学系からの光を受光して撮像する撮像手段と、前記撮像手段による画像データを画像処理する画像処理手段とを具備した撮像装置において、前記撮像手段は撮像素子で成り、前記撮像素子を走査することにより前記画像データを生成し、前記撮像素子が受光する点間の距離が、前記光学系の開口数NAによって定まるカットオフ空間周波数の2倍になるように、前記撮像素子の走査速度及び前記撮像素子の受光の時間間隔を設定して前記走査を行い、前記画像データが前記カットオフ空間周波数より高い周波数成分を含まないようにすることを特徴とする空間フィルタ処理を行う撮像装置。 An imaging apparatus comprising: an imaging unit that receives light from an optical system to capture an image; and an image processing unit that performs image processing on image data obtained by the imaging unit. The imaging unit includes an imaging element, and scans the imaging element. The image data is generated, and the scanning speed of the image sensor and the distance between the points received by the image sensor are twice the cutoff spatial frequency determined by the numerical aperture NA of the optical system. An imaging apparatus for performing spatial filtering, wherein the scanning is performed by setting a time interval of light reception of the imaging element so that the image data does not include a frequency component higher than the cutoff spatial frequency. 光学系からの光を受光して撮像する撮像手段と、前記撮像手段による第1画像データを画像処理する画像処理手段とを具備した撮像装置において、前記撮像手段が、撮像素子を走査することにより前記第1画像データを生成して出力する手段で成り、前記画像処理手段は、前記光学系の開口数NAに基づいてカットオフ空間周波数を決定する機能と、前記撮像素子が受光を行う点間の距離が前記カットオフ空間周波数の2倍となるデータを前記第1画像データから間引いて第2画像データを生成する機能とを具備し、前記第2画像データが前記カットオフ空間周波数より高い周波数成分を含まないようにしたことを特徴とする空間フィルタ処理を行う撮像装置。 In an imaging apparatus including an imaging unit that receives light from an optical system and captures an image, and an image processing unit that performs image processing on the first image data by the imaging unit, the imaging unit scans the imaging element. The first image data is generated and output. The image processing unit has a function of determining a cutoff spatial frequency based on a numerical aperture NA of the optical system, and a point where the image sensor receives light. The second image data is generated by thinning out data from which the distance is twice the cut-off spatial frequency from the first image data, and the second image data has a frequency higher than the cut-off spatial frequency. An imaging apparatus that performs spatial filter processing characterized by not including a component. 前記開口数NAを外部から入力するようになっている請求項11乃至17のいずれかに記載の空間フィルタ処理を行う撮像装置。 The imaging device for performing spatial filter processing according to claim 11, wherein the numerical aperture NA is input from the outside. 前記画像データ若しくは前記第1画像データの光強度がピーク値から略半分になる位置までをスポット像とし、前記スポット像から前記スポット像の直径dを算出し、前記開口数NAを
NA=cλ/πd
(ただし、cは(2J(x))=0.5の最小の解であり、約1.616 34の定数、λは光の波長である。)
として求めるようになっている請求項11乃至17のいずれかに記載の空間フィルタ処理を行う撮像装置。
A spot image is defined as a spot image where the light intensity of the image data or the first image data is approximately half of the peak value, the diameter d of the spot image is calculated from the spot image, and the numerical aperture NA is calculated.
NA = cλ / πd
(Where c is the minimum solution of (2J 1 (x)) 2 = 0.5, a constant of about 1.61634, and λ is the wavelength of light.)
An imaging apparatus for performing spatial filter processing according to claim 11, wherein the imaging apparatus performs spatial filter processing according to claim 11.
前記画像データのスポット像の第1暗環をエッジ抽出法等で抽出し、前記第1暗環の直径dを求め、前記開口数NAを
NA=cλ/πd
(ただし、定数cは(2J(x))=0の最小の解であり、約3.832 の定数、λは光の波長である。)
として求めるようになっている請求項11乃至17のいずれかに記載の撮像装置。
The first dark ring of the spot image of the image data is extracted by an edge extraction method or the like, the diameter d of the first dark ring is obtained, and the numerical aperture NA is calculated.
NA = cλ / πd
(However, the constant c is the minimum solution of (2J 1 (x)) 2 = 0, a constant of about 3.832, and λ is the wavelength of light.)
The imaging device according to claim 11, wherein the imaging device is obtained as follows.
JP2007198603A 2007-07-31 2007-07-31 Imaging method performing spatial filtering and its imaging apparatus Pending JP2009038439A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007198603A JP2009038439A (en) 2007-07-31 2007-07-31 Imaging method performing spatial filtering and its imaging apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007198603A JP2009038439A (en) 2007-07-31 2007-07-31 Imaging method performing spatial filtering and its imaging apparatus

Publications (1)

Publication Number Publication Date
JP2009038439A true JP2009038439A (en) 2009-02-19

Family

ID=40440014

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007198603A Pending JP2009038439A (en) 2007-07-31 2007-07-31 Imaging method performing spatial filtering and its imaging apparatus

Country Status (1)

Country Link
JP (1) JP2009038439A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014096640A (en) * 2012-11-07 2014-05-22 Canon Inc Image processing method, program, image processing apparatus, imaging device, and network equipment
WO2016098348A1 (en) * 2014-12-19 2016-06-23 Canon Kabushiki Kaisha Image processing apparatus, image processing method, image capturing apparatus and image processing program

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014096640A (en) * 2012-11-07 2014-05-22 Canon Inc Image processing method, program, image processing apparatus, imaging device, and network equipment
WO2016098348A1 (en) * 2014-12-19 2016-06-23 Canon Kabushiki Kaisha Image processing apparatus, image processing method, image capturing apparatus and image processing program
JP2016119542A (en) * 2014-12-19 2016-06-30 キヤノン株式会社 Image processing method, image processing program, image processing unit, and imaging apparatus
US10326951B2 (en) 2014-12-19 2019-06-18 Canon Kabushiki Kaisha Image processing apparatus, image processing method, image capturing apparatus and image processing program

Similar Documents

Publication Publication Date Title
EP2300983B1 (en) System and method for producing an optically sectioned image using both structured and uniform illumination
JP4914303B2 (en) Image processing apparatus and imaging apparatus, image processing method and imaging method, and image processing program
US9049356B2 (en) Image processing method, image processing apparatus and image processing program
JP2011211669A (en) Image processing apparatus, image pickup apparatus and image processing program
JP2012005056A (en) Image processing device, image processing method and program
JP2011166588A (en) Imaging apparatus and imaging method, and image processing method for the imaging apparatus
JP2011095027A (en) Imaging device
JP2021007247A5 (en)
US7733411B2 (en) Image capturing apparatus performing filtering process with variable cut-off frequency
JP4544336B2 (en) Image processing apparatus, imaging apparatus, image processing method, and program
US8648936B2 (en) Image processing apparatus and image processing method
JP2009038439A (en) Imaging method performing spatial filtering and its imaging apparatus
JP6682184B2 (en) Image processing method, image processing program, image processing device, and imaging device
JP5682443B2 (en) Image processing apparatus, image processing method, and image processing program
JP5115297B2 (en) Image processing apparatus, imaging apparatus, image processing method, and program
JP2007094236A (en) Imaging apparatus and method
JP2006229702A (en) Imaging apparatus and focus control method
JP2009133642A (en) Apparatus and program for image processing and imaging apparatus
JP2019029834A (en) Imaging apparatus
CN113419307B (en) Optical fiber imaging system
JP5914092B2 (en) Image processing system and microscope system including the same
JP6468751B2 (en) Image processing apparatus, image processing method, and program
JP5868758B2 (en) Image processing system and microscope system including the same
JP2017053710A (en) Image processing apparatus and image processing method
JP2017123532A (en) Image processing device, imaging device, image processing method, image processing program and recording medium