KR20040083368A - 스테레오-ｐｉｖ 방법의 실행에 있어 자동 캘리브레이션을위한 공액 거리 방정식 결정 방법 - Google Patents

스테레오-ｐｉｖ 방법의 실행에 있어 자동 캘리브레이션을위한 공액 거리 방정식 결정 방법 Download PDF

Info

Publication number
KR20040083368A
KR20040083368A KR1020040018601A KR20040018601A KR20040083368A KR 20040083368 A KR20040083368 A KR 20040083368A KR 1020040018601 A KR1020040018601 A KR 1020040018601A KR 20040018601 A KR20040018601 A KR 20040018601A KR 20040083368 A KR20040083368 A KR 20040083368A
Authority
KR
South Korea
Prior art keywords
cameras
camera
light
determined
piv
Prior art date
Application number
KR1020040018601A
Other languages
English (en)
Inventor
베르나르트 비네케
Original Assignee
라비숑 게엠바하
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 라비숑 게엠바하 filed Critical 라비숑 게엠바하
Publication of KR20040083368A publication Critical patent/KR20040083368A/ko

Links

Classifications

    • A—HUMAN NECESSITIES
    • A23—FOODS OR FOODSTUFFS; TREATMENT THEREOF, NOT COVERED BY OTHER CLASSES
    • A23L—FOODS, FOODSTUFFS, OR NON-ALCOHOLIC BEVERAGES, NOT COVERED BY SUBCLASSES A21D OR A23B-A23J; THEIR PREPARATION OR TREATMENT, e.g. COOKING, MODIFICATION OF NUTRITIVE QUALITIES, PHYSICAL TREATMENT; PRESERVATION OF FOODS OR FOODSTUFFS, IN GENERAL
    • A23L11/00—Pulses, i.e. fruits of leguminous plants, for production of food; Products from legumes; Preparation or treatment thereof
    • A23L11/05—Mashed or comminuted pulses or legumes; Products made therefrom
    • A23L11/07—Soya beans, e.g. oil-extracted soya bean flakes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/18—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
    • G01P5/20—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance using particles entrained by a fluid stream
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/001—Full-field flow measurement, e.g. determining flow velocity and direction in a whole region at the same time, flow visualisation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/18—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
    • G01P5/22—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
    • A—HUMAN NECESSITIES
    • A23—FOODS OR FOODSTUFFS; TREATMENT THEREOF, NOT COVERED BY OTHER CLASSES
    • A23V—INDEXING SCHEME RELATING TO FOODS, FOODSTUFFS OR NON-ALCOHOLIC BEVERAGES AND LACTIC OR PROPIONIC ACID BACTERIA USED IN FOODSTUFFS OR FOOD PREPARATION
    • A23V2002/00—Food compositions, function of food ingredients or processes for food or foodstuffs
    • A—HUMAN NECESSITIES
    • A23—FOODS OR FOODSTUFFS; TREATMENT THEREOF, NOT COVERED BY OTHER CLASSES
    • A23V—INDEXING SCHEME RELATING TO FOODS, FOODSTUFFS OR NON-ALCOHOLIC BEVERAGES AND LACTIC OR PROPIONIC ACID BACTERIA USED IN FOODSTUFFS OR FOOD PREPARATION
    • A23V2200/00—Function of food ingredients
    • A23V2200/20—Ingredients acting on or related to the structure
    • A23V2200/218—Coagulant

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Botany (AREA)
  • Agronomy & Crop Science (AREA)
  • Health & Medical Sciences (AREA)
  • Nutrition Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Polymers & Plastics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Measurement Of Optical Distance (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

본 발명은 적어도 두 개의 카메라와 하나의 영상영역을 포함하고, 카메라들은 광 슬릿의 거의 같은 영역을 다른 방향에서 바라보며, 두 카메라 사이의 점대응이 광학적인 상호관계를 사용하여 카메라 영상들에서 각 인터로게이션 필드의 이동을 측정하여 결정되고, 알려진 내부와 외부 카메라 파라미터를 사용하여 근사법으로 공액 거리 방정식이 결정되는 것을 특징으로 하는 시각화된 유동에서 스테레오-PIV 방법의 실행에 있어 자동 캘리브레이션을 위한 공액 거리 방정식 결정 방법에 관한 것이다.

Description

스테레오-ＰＩＶ 방법의 실행에 있어 자동 캘리브레이션을 위한 공액 거리 방정식 결정 방법{METHOD OF DETERMINATION OF CONJUGATE DISTANCE EQUATION FOR THE SELF-CALIBRATION REGARDING THE CARRYING OUT STEREO-PIV-METHOD}
본 발명은 스테레오-PIV 방법을 실행하는데 있어서 자동 캘리브레이션을 위한 공액 거리 방정식(conjugate distance equation)을 결정하기 위한 방법에 관한것이다.
먼저 PIV 방법에 대해 설명한다. PIV는 입자영상유속계(Particle Image Velocimetry)의 약자이다. PIV를 이용하여 공간에서 기체나 유체의 유동 상황을 나타낼 수 있다(예를 들면, DE 199 28 698 A1). 이러한 PIV 방법을 실행하기 위하여 먼저 레이저 또는 다른 광원이 필요하다. 광원은 기체나 액체와 같은 매개물의 유동에서 이른바 광 슬릿(light slit)이 생기게 하며, 광 슬릿은 적어도 하나의 카메라에 의해 관찰된다. 광 슬릿을 그에 수직으로 배열된 하나의 카메라만으로 관찰하면 빛 평면의 양 속도 성분만 결정될 수 있으나, 서로 다른 각도에서 관찰하는 적어도 두개의 카메라(스테레오-PIV)를 사용하면 세 개의 모든 성분을 결정할 수 있다. 언급한 바와 같이 PIV 방법은 2차원과 3차원의 속도장을 측정하는데 의미가 있다. 즉, 상기 매개물의 속도를 공간에서 시각화하기 위해 액체나 기체에 유동에 직접적으로 따라오는 작은 입자가 첨가된다. 이러한 스테레오-PIV 방법을 실행하기 위하여 먼저 캘리브레이션, 즉 광 절단면에 대한 카메라의 위치 결정이 필요하다. 이것은 결국 공액 거리 방정식
에 의해 결정된다. 여기서, xi, yi는 카메라 1과 2의 영상에서(도 1 참조) 복셀(x,y,z)의 영상좌표를 나타낸다. 일반적으로 좌표계는 광 절단면이 상수 z(예를 들면, z=0)에 일치되도록 정해진다. 영상함수 M 으로는 종종 핀홀카메라(Pinhole camera)모델이 쓰인다. 핀홀카메라 모델에서는 영상이 외부 파라미터 - 카메라들의 서로에 대한 방향과 위치 그리고 광 슬릿에 대한 방향과 위치 - 와 내부 파라미터 - 무엇보다도, 이미지너리 조리개(영상거리)에 대한 카메라 칩의 간격과 카메라 칩에서 광 주축의 기초점(main point) - 에 의해 결정된다. 왜곡 파라미터가 추가 감소는 영상을 0,1화소 보다 좋은 정확도로 결정하는 것이 가능하게 한다. 캘리브레이션은 선행기술에 의해 이른바 캘리브레이션 플레이트에 의해 실행된다. 캘리브레이션 플레이트는 양 카메라에 의해 공간의 한 위치 또는 여러 위치로 삽입된다. 이때, 그 위치는 광 절단면에 정확히 일치해야만 한다.
공액 거리 방정식은 캘리브레이션 플레이트와 공간에서 양 카메라의 절대 위치에 의하여 얻을 수 있다. 또는, 캘리브레이션 플레이트, 캘리브레이션 플레이트에 대한 카메라들의 각도와 방향, 캘리브레이션 플레이트에 대한 카메라들의 간격에 의하여 얻을 수 있다. 또는, 둘 또는 여러 z-위치들에서 카메라들에 의해 파악되는 캘리브레이션 플레이트에 의하여 얻을 수 있다.
또한, 공액 거리 방정식을 결정하기 위해 이른바 3차원의 캘리브레이션 플레이트를 사용하는 것이 알려져 있다. 3차원의 캘리브레이션 플레이트는 예를 들면 두 개의 평면을 가지며, 각 평면에는 일정한 간격으로 10 x 10 마킹의 스크린이 제공된다. 종래의 이와 같은 캘리브레이션 방법은 여러 단점들을 갖는다. 즉, 캘리브레이션 플레이트는 동일 위치에 배치되어야 하고, 빛과 정확히 평행하게 위치해야만 한다. 이것을 실현하려면 비용이 매우 많이 든다. 0.6°의 작은 편차조차도 두 개의 영상영역에서 벡터결정시 픽쳐 보더(picture border)에서 10 화소의 위치 오류를 초래한다. 두 영상영역은 높은 속도 구배에서 큰 퍼센트 범위의 오류를 일으킨다. 캘리브레이션을 위한 비용이 많이 든다. 큰 관찰필드에서는 이에 맞는 큰 캘리브레이션 플레이트가 제작되어야 하고, 경우에 따라 Z-방향으로 정확한 양만큼 일정하게 이동되어야 한다. 또는, 각도결정 및/또는 간격결정이 필요하다. 이것도 역시 비용이 많이 들며 오류가 발생한다. 예를 들면, 간격 결정시 캘리브레이션 플레이트 위의 제로 포인트와 이미지너리 카메라 구멍위치 사이의 거리를 결정하는 것이 어렵다. 하지만 이것은 다중 렌즈를 가진 일반적인 대물렌즈에서는 대물렌즈의 특정 위치에 있다. 관 내부와 같은 폐쇄된 공간에서 캘리브레이션을 하거나 PIV 방법을 실행할 때는 캘리브레이션 플레이트를 배치하기 위하여 관 내부로 들어갈 수 있어야 한다. 동시에, 캘리브레이션이 측정과 유사한 광학적 조건에서 이루어져야 한다. 즉, 경우에 따라서는 캘리브레이션은 후에 측정이 이루어지는 관 내부와 같은 액체의 조건하에서 이루어져야 한다.
마이크로 채널에서와 같은 많은 대물렌즈에서는 소위 인-시츄-캘리브레이션 (in-situ-calibration)을 전혀 실행할 수 없거나 실행을 위하여 높은 비용을 필요로 한다. 캘리브레이션 플레이트를 배치하는 것이 매우 어렵기 때문이다.
같은 또는 유사한 이유에서, 컴퓨터 영상 및 사진측량분야에서는 캘리브레이션의 도움없이도 충분히 정확한 공액 거리 방정식에 이르는 여러 방법들이 개발되었다. 자동 캘리브레이션이라 불리우는 이와 같은 방법은, 같은 점들, 이른바 점대응이 두 개의 카메라 영상들에서 찾아진다는 것에 기초를 두고 있다. 상기 점들은 같은 복셀(voxel)에 속한다. 충분히 많은 점대응들을 알 수 있고 개별적인 내부 카메라 파라미터들을 전부 또는 부분적으로 알수 있으며 정확한 척도를 알수 있으면, 상기의 공액 거리 방정식, 즉 나머지 내부 카메라 파라미터와 카메라들 서로의 방향과 간격을 결정하는 것이 가능하다. 하지만 이러한 방법을 곧바로 스테레오-PIV 방법에 적용할 수 없다. 왜냐하면 카메라들 간의 점대응들을 결정하는 것이 부분적으로 어렵기 때문이다. 또한, 단단한 구조를 가진 단단한 표면이 관찰되는 것이 아니라, 광 슬릿에 의해 주어진 체적에서 움직이는 입자들이 관찰되기 때문이다.
또한, DE 198 01 615 A1에는 광절단법(light cut method)을 위한 캘리브레이션 방법이 개시되어 있다. 이 방법에서는 평가단위의 캘리브레이션이, 영상축척으로 카메라에 의해 유동물체에서 찍힌 영상을 또 다른 영상축척으로 유동물체 외부에서 찍힌 영상과 양적으로 비교하여 실행된다. 이와 같은 방법에서는 카메라들이 아주 빨리 움직여져야 한다는 단점이 있다.
따라서, 본 발명의 목적은 상기 언급한 단점들을 제거하여 캘리브레이션에 대한 비용을 낮추고 폐쇄 공간 또는 마이크로 채널에서도 실행할 수 있는 스테레오-PIV 방법의 캘리브레이션 가능성을 제공하는 데에 있다.
도 1은 일반적인 스테레오-PIV 구조의 개략도,
도 2는 카메라 1과 2의 상관관계로부터 상관관계 필드의 생성을 도시한 도면,
도 3은 첫 번째 레이저(왼쪽)와 두 번째 레이저(오른쪽)에서의 도 2의 상관관계 필드를 도시한 도면,
도 4는 가장 높은 상관관계 정점의 위치로부터 계산된 이동벡터를 도시한 도면(보다 나은 시각화를 위해 특정 요소만큼 확대됨)이다.
상기 목적은 본 발명에 따른 스테레오-PIV 방법의 실행에 있어 자동 캘리브레이션을 위한 공액 거리 방정식을 결정하기 위한 방법에 의하여 달성될 수 있다.이 방법은 시각화된 유동에서 적어도 두 개의 카메라와 하나의 광 슬릿을 포함한다. 카메라들은 광 슬릿의 거의 같은 영역을 다른 방향에서 관찰하며, 적어도 두개의 카메라들 사이의 점대응의 결정은 광학적 상호관계에 의하여 카메라 영상에서 각 인터로게이션 필드(interrogation field)의 이동 측정을 통해 이루어지며, 잘 알려져 있는 내부와 외부 카메라 파라미터에 의하여 근사법(approximaiton method)을 적용하여 공액 거리 방정식이 결정된다. 본 발명에 의한 방법에서 핵심이 되는 것은 상기 언급한 적어도 두 카메라 사이의 점대응을 결정하는 것이다. 상기한 바와 같이 점대응의 결정은 이른바 광학적 상호관계에 의해 이루어진다. 카메라 영상은 광학적 상호관계를 가지고 특정 시점 t에서 카메라에 의해 찍힌다. 이때, 동일한 카메라 영상이 두 번째 카메라에 의해 같은 시점 t에서 다른 방향에서 찍힌다. 즉, 카메라 영상들은 각 경우 같은 섹션(section)을 나타내나, 관찰하는 카메라들의 렌즈때문에 영상들은 상대적으로 서로 이동되거나 비틀어지거나 일그러뜨려져 보인다. 카메라 영상의 이동 정도를 측정하기 위하여 각 카메라 영상은 개별 부분, 소위 인터로게이션 필드로 분할된다. 즉, 하나의 카메라 영상은 예를 들면 20 x 20 인터로게이션 필드로 구성된다. 하나의 인터로게이션 필드가 첫 번째 카메라 영상에서 결정되고, 이에 대응하는 인터로게이션 필드가 두 번째 카메라 영상에서 결정된다. 두 번째 영상 영역의 인터로게이션 필드에 대한 첫 번째 카메라 영상의 인터로게이션 필드의 거리는, 카메라 렌즈에 의해 관찰된 카메라 영상의 이동의 크기를 형성한다. 이러한 거리의 크기는 2차원의 상관함수(dx, dy)에서 가장 높은 상관관계 정점을 형성하고, 상관관계 필드에서의 이러한 정점의 위치는 각 카메라의 위치(x1, y1);(x2, y2)를 나타낸다. 이로써, 각 인터로게이션 필드에 대한 점 상관관계 (x1, y1<-> x2, y2)가 생긴다.
그 후, 하나 또는 복수의 내부 카메라 파라미터와 점대응, 절대적 길이 축척을 이용하여 나머지 내부 및 외부 파라미터가 결정될 수 있으며, 완전한 공액 거리 방정식은 레벤베르그-마쿼트-알고리즘(Levenberg-Marquardt-Algorithm)과 같은 근사법에 의하여 결정될 수 있다.
캘리브레이션이 이미 일어나 이미지의 내부와 외부 파라미터가 이미 인식되었으나 공간에서 광 슬릿의 위치가 아직 알려지지 않은 두 번째 전형적 경우에는, 양 레이저의 하나 또는 양쪽의 광 슬릿의 위치는 점대응에 의하여 공간에서 일반적인 삼각법(triangular method)에 의해 결정될 수 있다.
점대응이 잘못 결정될 위험이 항상 있다. 즉, 한 카메라 영상의 각 인터로게이션 필드가 다른 카메라 영상의 각 인터로게이션 필드와 맞지 않을 수 있다. 그러나, 알려져 있는 RANSAC(RAndom SAmpling Consensus) 기법에 의하여 이러한 잘못된 점대응은 제거될 수 있다. 즉, RANSAC법이 알고리즘으로서 본래의 근사법에 적용됨으로써 잘못된 점대응이 제거될 수 있다.
단단한 표면이 아니라 광 슬릿에 의해 한정되는 체적에서 분리된 입자들이 관찰되기 때문에, 입자들은 경우에 따라 한 카메라에서 볼때와 다른 카메라에 의해 관찰될 때 매우 다르게 보인다. 입자들은 공간-광 슬릿의 깊이-에 배치되기 때문이다. 따라서, 양 카메라 영상들 사이의 상호관계는 오류를 갖기가 매우 쉽다. 정확한 상관관계 정점이 매우 선명치 않게 되며 종종 소음 정점보다 낮아서 인식되지 않기 때문이다. 이러한 오류 가능성은 바람직하게는 t0에서 tn까지 연속 시간 동안 적어도 두 개의 카메라에서 매번 두 개 또는 그 이상의 카메라 영상들이 찍힘으로써 제거된다. 이때, 각 시점 t0에서 tn까지에서 이 영상들에 의하여 C0(dx, dy)부터 Cn(dx, dy)까지 2차원 상관함수가 광학적 상호관계를 이용하여 결정된다. 이때, 상관함수 C0부터 Cn까지가 축적되며, 가장 높은 상관관계 정점의 결정 후에 각 인터로게이션 필드의 이동 dx, dy와 점대응이 결정된다.
이하, 본 발명을 첨부된 도면을 참조한 실시예를 통하여 상세히 설명한다.
실시예 1:
두 개의 카메라를 가진 일반적인 스테레오-PIV 구조(도 1과 도2)가 기초가 된다. 카메라들은 x축을 따라 위치되며, 마주보는 면에서 30 내지 50°의 각도 하에서 광 절단면으로 향해 있으며, 광 절단면은 z=0인 x-y평면으로 정의된다. 따라서 양 카메라는 z=-Zcam에 놓여 있다. 카메라들의 광 주축은 동일평면에 있는 공통 x-y평면에 놓여 있다. 두 개의 펄스 레이저(3)가 광 절단 렌즈(4)를 써서 같은 위치에서 짧게 잇달아서 광 슬릿(5)을 생성한다. 이때, 양 카메라는 두 개의 영상(6)을 짧게 연속하여 찍는다(각 영상에 하나의 레이저 펄스).
본 실시예는, 예를 들면 3차원 캘리브레이션 플레이트가 양 카메라에 의해동시에 취해지는 것과 같이, 체적 캘리브레이션이 본래의 광 슬릿에 독립하여 이미 일어난 것으로 가정한다. 따라서, 모든 내부와 외부 영상 파라미터는 좌표계에 관하여 캘리브레이션 플레이트의 위치를 기초로 한다.
실제 측정시 동시에 찍힌 입자 영상들 간의 광학적 상호관계를 사용하여, 서로 다른 시간에 찍힌 상관관계 필드에 대한 전달시 각 인터로게이션 필드(도 3, 1번)에 대해 증가된 상관관계 필드가 결정되며, 가장 높은 상관관계 정점(도 3, 2번-영상 3에 있는 화살표에 상응함)의 위치로 부터 카메라 1과 2사이의 점대응들(영상 4)이 생긴다. 화살표의 시작점은 카메라의 1의 영상에서 인터로게이션 윈도우의 위치를 나타내고, 끝점은 카메라 2의 영상에서 대응하는 점을 나타낸다. 시작점과 끝점은 함께 점대응을 형성한다.
그 후 점대응으로부터 알려져 있는 공액 거리 방정식을 사용하여 삼각법에 의하여 광 슬릿의 절대적 위치가 공간에서 결정되며, 적절한 좌표 변환에 의해 광 절단면이 z=0으로 정의된다. 따라서 광 절단면을 위한 영상이 결정되고, 실제 스테레오-PIV 평가를 위해 사용될 수 있다. 이 방법의 장점은 캘리브레이션 플레이트가 광 절단면에 정확하게 위치될 필요없이 공간 어디엔가에 위치되면 되고, 그럼에도 불구하고 광 절단면에서 아주 정확한 캘리브레이션이 계산될 수 있다는 데에 있다.
또한, 상관관계 정점의 너비(도 3, 3번)와, 쉽게 계산될 수 있는 기하학적 요소로 부터 직접 광 슬릿의 두께가 나온다. 도 3 왼쪽은 레이저 1의 상관관계 필드를 나타내며, 도 3 오른쪽은 레이저 2의 상관관계 필드를 나타낸다. 따라서, 공간에서의 양 광 슬릿의 상대적인 위치와 두께로 부터, 양 광 슬릿이 어느정도로 서로 겹치며 PIV 측정에 적합한지가 나타난다.
실시예 2:
실시예 1과 같은 구조를 전제로 한다. 이와 동시에, 샤임플러그 조건을 충족시키기 위해 카메라 대물렌즈가 카메라 평면에 대한 앵글(angel)에 있다고 가정하며, 따라서 모든 입자들이 광 절단면 초점에 있다. 하지만 본 실시예에서는 사전의 캘리브레이션이 없고, 점대응 자체로 부터 공액 거리 방정식이 결정되어야 한다. 이것은 주어지지 않은 영상 파라미터들을 조절하여 알아내는 직접적인 근사법에 의해 실행된다. 자유 파라미터들이 너무 많기 때문에 해답에 수렴하기 위하여 특정한 가정이 필요하다. 또한, 자유 파라미터의 수를 줄이기 위하여 알려진 상황들에 의한 다른 가능성들이 있다.
단지 한번 실행될 수 있는 샤임플러그 어댑터의 사전 캘리브레이션으로 부터, 기초점이 각도의 함수로써 어떻게 이동되는지 알려진다고 가정한다. 또는,샤임플러그 조건이 직접 기하 형상으로부터 계산된다고 가정한다. 따라서, 기초점은 조절하여 알아낼 필요가 없고, 조절하여 측정된 외부 관찰 각도의 함수이다.
영상거리도 마찬가지 방식으로 한다. 영상거리는 렌즈방정식 1/B + 1/G = 1/f로부터 계산된다(B=영상거리, G=물체거리, f=카메라 대물렌즈의 알려진 초점거리). 근사법 중 초점거리는 물체거리의 함수로 계산된다. 이때, G은 자유 외부 파라미터로써 서서히 알아져야 한다. 렌즈방정식을 사용하지 않고, 사전에 각 카메라에 대해 분리하여 물체거리에 대한 영상거리의 종속성을 경험적으로 캘리브레이션할 수도 있다.
또한, 이 경우 광주축이 동일평면인 것을 사용하여 자유 파라미터의 수를 더 줄일 수 있는 가능성이 있다.
이 방법의 장점은, 인시투 캘리브레이션 없이도 완전히 실행될 수 있으며, 공액 거리 방정식이 전적으로 몇 가지 알려진 상황을 가정하여 계산된 점대응으로부터 계산될 수 있다는 데에 있다. 따라서, PIV 측정과정이 매우 쉬워진다.
본 발명에 따른 방법은 캘리브레이션 비용을 줄일 수 있고, 폐쇄 공간 또는 마이크로 채널에서도 실행할 수 있는 장점이 있다.

Claims (14)

  1. 적어도 두 개의 카메라와 하나의 영상영역을 포함하고, 상기 카메라들은 광 슬릿의 거의 같은 영역을 다른 방향에서 바라보며, 두 카메라 사이의 점대응이 광학적인 상호관계를 사용하여 카메라 영상들에서 각 인터로게이션 필드의 이동을 측정함으로써 결정되고, 알려진 내부와 외부 카메라 파라미터를 사용하여 근사법으로 공액 거리 방정식이 결정되는 것을 특징으로 하는 시각화된 유동에서 스테레오-PIV 방법의 실행에 있어 자동 캘리브레이션을 위한 공액 거리 방정식 결정 방법.
  2. 제 1 항에 있어서,
    상기 내부 파라미터는 초점거리, 광축들의 위치(x0, y0), 카메라 렌즈의 왜곡 파라미터를 포함하는 것을 특징으로 하는 방법.
  3. 제 1 항에 있어서,
    상기 외부 파라미터는 카메라들의 서로에 대한 위치와 방향을 포함하는 것을 특징으로 하는 방법.
  4. 제 1 항 내지 제 3 항 중 어느 한 항에 있어서,
    알려진 공액 거리 방정식의 좌표계에 관한 광 슬릿의 알려지지 않은 위치에있어서, 광 슬릿의 위치는 점대응을 사용하여 결정되는 것을 특징으로 하는 방법.
  5. 제 1 항 내지 제 4 항 중 어느 한 항에 있어서,
    공액 거리 방정식을 결정하기 위하여, 점대응을 사용한 내부 파라미터의 하나 또는 여러 개의 정보와 절대적 길이 축척의 정보를 가지고, 나머지 내부와 외부 카메라 파라미터들을 결정하는 것을 특징으로 하는 방법.
  6. 제 1 항 내지 제 5 항 중 어느 한 항에 있어서,
    t0에서tn까지의 연속 시간 중 적어도 두 개의 카메라에서 매번 두 개 또는 그 이상의 카메라 영상들이 촬영되고, 이때 t0에서 tn까지의 각 시점에서 이 영상들에 의하여 C0(dx, dy)부터 Cn(dx, dy)까지 2차원 상관함수가 광학적인 상호관계에 의하여 결정되며, 이때 C0부터 Cn까지의 상관함수가 축적되고, 가장 높은 상관관계 정점의 결정 후에 각 인터로게이션 필드의 이동 dx, dy와 점대응이 결정되는 것을 특징으로 하는 방법.
  7. 제 1 항 내지 제 6 항 중 어느 한 항에 있어서,
    상기 근사법은 레벤베르그-마쿼트-알고리즘(Levenberg-Marquardt-Algorithm)에 기초하는 것을 특징으로 하는 방법.
  8. 제 7 항에 있어서,
    상기 레벤베르그-마쿼트-알고리즘에 RANSAC-알로리즘이 적용되는 것을 특징으로 하는 방법.
  9. 제 1 항에 있어서,
    상기 각 카메라에 의하여 두 개의 영상이 짧게 잇달아서 촬영되며, 영상들 간의 상호관계를 사용하여 시점 t와 t+dt에서 추가적인 점대응들이 결정되는 것을 특징으로 하는 방법.
  10. 제 1 항에 있어서,
    상기 적어도 두 개의 카메라의 광축들은 동일평면에 배치되는 것을 특징으로 하는 방법.
  11. 제 6 항에 있어서,
    상관관계 정점의 너비와 기하 요소에 의해 양 광 슬릿의 광 절단 두께가 결정되며, 따라서 공간에서의 광 슬릿의 위치와 함께 광 슬릿들이 서로 어느 정도로 겹치며 PIV 측정에 적합한지가 정확히 결정될 수 있는 것을 특징으로 하는 방법.
  12. 제 5 항에 있어서,
    상기 근사법에서 광 슬릿에서 입자들에 초점을 맞추는 것을 가정한 경우, 영상거리는 대물렌즈의 초점거리와 카메라에 대한 광 슬릿의 간격의 함수로 계산되며 더 이상 조절하여 알아낼 필요가 없는 것을 특징으로 하는 방법.
  13. 제 5 항에 있어서,
    상기 근사법에서 샤임플러그어댑터의 최적 조절을 가정하여 샤임플러그어댑터를 사용할때, 카메라 칩과 주축 간의 각도와 카메라 칩에서 기초점의 위치가 외부 영상 파라미터에 의해 계산되며 더 이상 조절하여 알아낼 필요가 없는 것을 특징으로 하는 방법.
  14. 제 6 항에 있어서,
    상관관계 정점의 너비와 영상 기하 형상에 의해 양 광 슬릿의 광 절단 두께가 결정되며, 따라서 공간에서의 광 슬릿들의 위치와 함께 광 슬릿들이 서로 어느 정도로 겹치며 스테레오-PIV 방식에 적합한지 정확히 결정될 수 있는 것을 특징으로 하는 방법.
KR1020040018601A 2003-03-21 2004-03-18 스테레오-ｐｉｖ 방법의 실행에 있어 자동 캘리브레이션을위한 공액 거리 방정식 결정 방법 KR20040083368A (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10312696.1 2003-03-21
DE10312696A DE10312696B3 (de) 2003-03-21 2003-03-21 Verfahren zur Bestimmung der Abbildungsgleichung für die Selbstkalibrierung in Bezug auf die Durchführung von Stereo-PIV-Verfahren

Publications (1)

Publication Number Publication Date
KR20040083368A true KR20040083368A (ko) 2004-10-01

Family

ID=32798025

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020040018601A KR20040083368A (ko) 2003-03-21 2004-03-18 스테레오-ｐｉｖ 방법의 실행에 있어 자동 캘리브레이션을위한 공액 거리 방정식 결정 방법

Country Status (6)

Country Link
US (1) US20040183909A1 (ko)
EP (1) EP1460433B1 (ko)
JP (1) JP2004286733A (ko)
KR (1) KR20040083368A (ko)
CN (1) CN1536366A (ko)
DE (1) DE10312696B3 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150124069A (ko) * 2014-04-25 2015-11-05 주식회사 고영테크놀러지 카메라 파라미터 산출 방법

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10343160B3 (de) 2003-09-18 2005-05-25 Lavision Gmbh Verfahren zur Bestimmung eines dreidimensionalen Geschwindigkeitsfeldes in einem Volumen
CN100443899C (zh) * 2005-08-19 2008-12-17 北京航空航天大学 一种叶轮机械内部流场测量方法
JP4545666B2 (ja) * 2005-09-20 2010-09-15 株式会社フォトロン 流体計測装置および流体計測方法
DE102006002794A1 (de) 2006-01-20 2007-07-26 Wash Tec Holding Gmbh Verfahren und Vorrichtung zur Steuerung einer Fahrzeugwaschanlage
DE102006055746A1 (de) 2006-11-25 2008-05-29 Lavision Gmbh Verfahren zur Korrektur einer Volumenabbildungsgleichung zur Bestimmung eines Geschwindigkeitsfeldes von Teilchen in einem Volumen
DE102007013221B3 (de) * 2007-03-15 2008-08-14 Ksa Kugelstrahlzentrum Aachen Gmbh Verfahren zur Bestimmung des Abstandes zwischen zwei Erfassungsstellen einer Geschwindigkeitsmesseinrichtung für Partikel
CN100458372C (zh) * 2007-03-22 2009-02-04 同济大学 一种对建筑、城市空间进行精确测量的粒子图像测速方法
US20080306708A1 (en) * 2007-06-05 2008-12-11 Raydon Corporation System and method for orientation and location calibration for image sensors
JP5312236B2 (ja) * 2009-07-08 2013-10-09 本田技研工業株式会社 ３次元空間の粒子画像流速測定装置
DE112010002857T8 (de) 2009-07-08 2013-02-07 Honda Motor Co., Ltd. Teilchenbildgeschwindigkeitsmessungsverfahren, Teilchenbildgeschwindigkeitsmessungsverfahren für einen dreidimensionalen Raum, Teilchenbildgeschwindigkeitsmessungssystem und Tracerteilchen-Erzeugungsvorrichtung in einem Teilchenbildgeschwindigkeitsmessungssystem
WO2011058899A1 (ja) * 2009-11-10 2011-05-19 本田技研工業株式会社 ３次元空間の音源分布測定装置
TWI404894B (zh) * 2009-12-22 2013-08-11 Ind Tech Res Inst 照明系統
DE102009060580B4 (de) * 2009-12-23 2013-02-21 Lavision Gmbh Verfahren zur Bestimmung eines Satzes von optischen Abbildungsfunktionen für die 3D-Strömungsmessung
CN102331510B (zh) * 2011-06-09 2013-02-13 华南理工大学 纸浆两相流piv测量的图像处理方法
CN102291530A (zh) * 2011-06-17 2011-12-21 河海大学 自动调节piv摄像机位置的方法及其装置
CN102331511B (zh) * 2011-06-17 2014-05-07 河海大学 Piv图像高频采集方法
TW201500735A (zh) * 2013-06-18 2015-01-01 Nat Applied Res Laboratories 行動式影像流速辨識之方法及其裝置
DE202016100728U1 (de) 2015-05-06 2016-03-31 Lavision Gmbh Scheimpflugadapter und Verwendung
CN105785066B (zh) * 2016-03-28 2019-03-05 上海理工大学 凸曲面容器内流场微粒成像测速技术的径向畸变校正方法
CN106127724B (zh) * 2016-05-06 2018-09-21 北京信息科技大学 用于场相关畸变模型的标定场设计及标定方法
CN106895232B (zh) * 2017-03-07 2018-08-10 清华大学 一种用于tpiv测量的4台相机协同调节的安装平台
DE102017002235A1 (de) * 2017-03-08 2018-09-13 Blickfeld GmbH LIDAR-System mit flexiblen Scanparametern
US10186051B2 (en) 2017-05-11 2019-01-22 Dantec Dynamics A/S Method and system for calibrating a velocimetry system
DE102017110533A1 (de) 2017-05-15 2018-11-15 Lavision Gmbh Verfahren zum Kalibrieren eines optischen Messaufbaus
DE102018131059A1 (de) * 2018-12-05 2020-06-10 SIKA Dr. Siebert & Kühn GmbH & Co. KG Strömungsmessverfahren und Strömungsmessvorrichtung zur optischen Strömungsmessung
DE102019103441A1 (de) 2019-02-12 2020-08-13 Voith Patent Gmbh Verfahren zur Kalibrierung einer PIV Messanordnung
CN109946478A (zh) * 2019-03-24 2019-06-28 北京工业大学 一种针对空气静压主轴内部气体流速的检测系统
CN115932321B (zh) * 2022-12-22 2023-10-10 武汉大学 基于粒子图像测速的微观溶蚀可视化装置及方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4980690A (en) * 1989-10-24 1990-12-25 Hughes Aircraft Company Bistatic radar seeker with range gating
DE4408072C2 (de) * 1994-02-01 1997-11-20 Deutsche Forsch Luft Raumfahrt Verwendung einer elektronischen Hochgeschwindigkeitskamera bei einem Verfahren zur Bestimmung von Strömungsgeschwindigkeiten in einer Strömung
US5699444A (en) * 1995-03-31 1997-12-16 Synthonics Incorporated Methods and apparatus for using image data to determine camera location and orientation
EP0829726B1 (de) * 1996-09-05 2003-04-16 Robert Bosch Gmbh Verfahren und Vorrichtung zum Erfassen 3-dimensionaler Strömungsstrukturen
US5905568A (en) * 1997-12-15 1999-05-18 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Stereo imaging velocimetry
DE19801615A1 (de) * 1998-01-17 1999-07-29 Bosch Gmbh Robert Kalibrierverfahren für Laserlichtschnittverfahren
US6278460B1 (en) * 1998-12-15 2001-08-21 Point Cloud, Inc. Creating a three-dimensional model from two-dimensional images
DE19928698A1 (de) * 1999-06-23 2000-09-21 Deutsch Zentr Luft & Raumfahrt Vorrichtung zur Durchführung von PIV-Messungen
US6990228B1 (en) * 1999-12-17 2006-01-24 Canon Kabushiki Kaisha Image processing apparatus
US6789039B1 (en) * 2000-04-05 2004-09-07 Microsoft Corporation Relative range camera calibration
WO2002093916A2 (en) * 2001-05-14 2002-11-21 Elder James H Attentive panoramic visual sensor
US6542226B1 (en) * 2001-06-04 2003-04-01 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Planar particle imaging and doppler velocimetry system and method
US7257237B1 (en) * 2003-03-07 2007-08-14 Sandia Corporation Real time markerless motion tracking using linked kinematic chains
US8064684B2 (en) * 2003-04-16 2011-11-22 Massachusetts Institute Of Technology Methods and apparatus for visualizing volumetric data using deformable physical object
DE10343160B3 (de) * 2003-09-18 2005-05-25 Lavision Gmbh Verfahren zur Bestimmung eines dreidimensionalen Geschwindigkeitsfeldes in einem Volumen

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150124069A (ko) * 2014-04-25 2015-11-05 주식회사 고영테크놀러지 카메라 파라미터 산출 방법

Also Published As

Publication number Publication date
EP1460433B1 (de) 2012-02-08
CN1536366A (zh) 2004-10-13
US20040183909A1 (en) 2004-09-23
DE10312696B3 (de) 2004-12-23
JP2004286733A (ja) 2004-10-14
EP1460433A2 (de) 2004-09-22
EP1460433A3 (de) 2007-01-24

Similar Documents

Publication Publication Date Title
KR20040083368A (ko) 스테레오-ｐｉｖ 방법의 실행에 있어 자동 캘리브레이션을위한 공액 거리 방정식 결정 방법
Prasad Stereoscopic particle image velocimetry
KR100927096B1 (ko) 기준 좌표상의 시각적 이미지를 이용한 객체 위치 측정방법
US20030035100A1 (en) Automated lens calibration
US20020149691A1 (en) Aperture coded camera for three dimensional imaging
JP6503221B2 (ja) ３次元情報取得装置、及び、３次元情報取得方法
CN101539422B (zh) 一种单目视觉实时测距方法
CN111192235B (zh) 一种基于单目视觉模型和透视变换的图像测量方法
US20080123939A1 (en) Method of correcting a volume imaging equation for more accurate determination of a velocity field of particles in a volume
CN111325801A (zh) 一种激光雷达和相机的联合标定方法
CN102810205A (zh) 一种摄像或照相装置的标定方法
Pedersini et al. Accurate and simple geometric calibration of multi-camera systems
JP2008267843A (ja) トンネル切羽面の測量システム
ES2894935T3 (es) Aparato de medición de distancias tridimensionales y procedimiento para el mismo
Legarda et al. A new method for Scheimpflug camera calibration
Zhao et al. Volumetric calibration for scheimpflug light-field PIV
Diebold et al. Light-field camera design for high-accuracy depth estimation
US8629978B1 (en) Curved laser sheet for conformal optical diagnostics
JP2003346130A (ja) ３次元情報処理装置および３次元情報処理方法
US6616347B1 (en) Camera with rotating optical displacement unit
WO2020075213A1 (ja) 計測装置、計測方法および顕微鏡システム
CN115049784A (zh) 一种基于双目粒子图像的三维速度场重建方法
Iovenitti et al. Three-dimensional measurement using a single image
Lai et al. Volumetric three-component velocimetry: a new tool for 3d flow measurement
JPS6256814A (ja) ３次元位置計測カメラ較正方式

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E601 Decision to refuse application