• 제목/요약/키워드: correction parameter

검색결과 299건 처리시간 0.028초

Application of X-band polarimetric radar observation for flood forecasting in Japan

  • Kim, Sun-Min;Yorozu, Kazuaki;Tachikawa, Yasuto;Shiiba, Michiharu
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2011년도 학술발표회
    • /
    • pp.15-15
    • /
    • 2011
  • The radar observation system in Japan is operated by two governmental groups: Japan Meteorological Agency (JMA) and the Ministry of Land, Infrastructure, Transport and Tourism (MLIT) of Japan. The JMA radar observation network is comprised of 20 C-band radars (with a wavelength of 5.6 cm), which cover most of the Japan Islands and observe rainfall intensity and distribution. And the MLIT's radar observation system is composed of 26 C-band radars throughout Japan. The observed radar echo from each radar unit is first modified, and then sent to the National Bureau of Synthesis Process within the MLIT. Through several steps for homogenizing observation accuracy, including distance and elevation correction, synthesized rainfall intensity maps for the entire nation of Japan are generated every 5 minutes. The MLIT has recently launched a new radar observation network system designed for flash flood observation and forecasting in small river basins within urban areas. It is called the X-band multi parameter radar network, and is distinguished by its dual polarimetric wave pulses of short length (3cm). Attenuation problems resulting from the short wave length of radar echo are strengthened by polarimetric wavelengths and very dense radar networks. Currently, the network is established within four areas. Each area is observed using 3-4 X-band radars with very fine resolution in spatial (250 m) and temporal (1 minute intervals). This study provides a series of utilization procedures for the new input data into a real-time forecasting system. First of all, the accuracy of the X-band radar observation was determined by comparing its results with the rainfall intensities as observed by ground gauge stations. It was also compared with conventional C-band radar observation. The rainfall information from the new radar network was then provided to a distributed hydrologic model to simulate river discharges. The simulated river discharges were evaluated again using the observed river discharge to estimate the applicability of the new observation network in the context of operations regarding flood forecasting. It was able to determine that the newly equipped X-band polarimetric radar network shows somewhat improved observation accuracy compared to conventional C-band radar observation. However, it has a tendency to underestimate the rainfall, and the accuracy is not always superior to that of the C-band radar. The accuracy evaluation of the X-band radar observation in this study was conducted using only limited rainfall events, and more cases should be examined for developing a broader understanding of the general behavior of the X-band radar and for improving observation accuracy.

  • PDF

고해상도 강수정보 생산을 위한 레이더 반사도-강수량 관계식 매개변수 보정 및 불확실성 평가 (Estimation of reflectivity-rainfall relationship parameters and uncertainty assessment for high resolution rainfall information)

  • 김태정;김장경;김진국;권현한
    • 한국수자원학회논문집
    • /
    • 제54권5호
    • /
    • pp.321-334
    • /
    • 2021
  • 일반적으로 레이더 강수량 추정에 활용되는 Marshall-Palmer 관계식은 강수현상의 계절적 변동성을 고려하지 않고 선별된 강수사상에 대하여 시공간적으로 고정된 매개변수를 적용하여 레이더 강수량을 추정하므로 실제 강수량과 추정된 레이더 강수량은 정량적인 오차가 발생할 수 있다. 이러한 제약성을 극복하고자 본 연구는 장기간 레이더 반사도 인자를 가용하여 레이더 반사도-강수량 관계식 매개변수를 Bayesian 추론기법으로 보정하고 불확실성을 정량화하여 레이더 강수량의 편의 보정을 수행하였다. Bayesian 추론기법 기반으로 추정된 레이더 반사도-강수량 관계식의 보정 매개변수는 계절성이 규명되었으며 지역적 특성이 존재하였다. Bayesian 추론기법을 통하여 산정된 레이더 강수량은 Marshall-Palmer 관계식의 과소추정 문제를 극복하고 지상 강수특성을 정량적으로 현실성 있게 재현하였다. 본 연구결과는 집중호우 발생 시 능동적인 유역단위 수자원 해석 시스템을 구축하여 국가적 레이더 자원의 가치를 향상할 수 있을 것으로 판단된다.

탄성파 탐사를 활용한 산지사면 토심 추정 및 예측모델 보정 (Soil Depth Estimation and Prediction Model Correction for Mountain Slopes Using a Seismic Survey)

  • 봉태호;임상준;서정일;김동엽;허준
    • 한국산림과학회지
    • /
    • 제112권3호
    • /
    • pp.340-351
    • /
    • 2023
  • 산사태는 매년 막대한 재산 피해와 인명 피해를 유발하는 주요 자연재해 중 하나이며, 기후변화의 영향으로 산사태 취약성은 더욱 증대되고 있다. 사면안정해석을 위한 다양한 매개변수 중 토심은 산사태 및 토석류 분석 시 주요 매개변수이며, 사면 안정성에 영향을 미치는 유역의 수문학적 과정을 평가하는 데 중요한 역할을 한다. 토심을 추정하는 정확한 방법은 현장에서 직접 지층을 조사하는 것이다. 하지만 이를 위해서는 많은 시간과 비용이 요구되므로 다양한 토심 예측 모델들이 제안되었으나 실용성 및 정확성 측면에서 다양한 한계가 존재한다. 본 연구에서는 산지사면에 대한 토심을 추정 하기 위하여 국내 산지를 대상으로 수행된 71개의 탄성파 탐사 결과를 수집하였으며 탄성파 속도 700 m/s를 기준으로 토심을 추정하였다. 이에 따라 사면의 경사, 고도, 토심 자료를 구축하고 토심에 대한 통계적 특성을 파악하였으며 경사와 토심 및 고도와 토심 간 상관관계를 규명하였다. 또한, 사면경사를 기반으로 한 다양한 토심 예측 모델을 조사하고 이에 대한 비교·분석을 수행하였으며 경사를 활용한 보정된 토심 예측식을 제안하였다.

레이저 출력에 따른 난류 모사 위상판 측정 (Measurement of a Phase Plate Simulates Atmospheric Turbulence Depending on Laser Power)

  • 오한결;강필성;이재현;이혁교;김영식
    • 한국광학회지
    • /
    • 제34권3호
    • /
    • pp.99-105
    • /
    • 2023
  • 천체망원경의 성능은 여러가지 요소에 의해 결정된다. 대기 난류도 그 중 하나인데, 대기 난류는 망원경으로 수집한 빛을 왜곡시켜 이미지의 선명도와 해상도를 저하시킨다. 때문에 대기 난류를 보정하기 위한 기술이 연구되어 왔다. 보정 기술을 연구하기 위해서는 대기 난류를 실험실에서 모사해야 하며, 그 중 가장 실용적인 방법으로 위상판을 이용한 방법이 있다. 심한 난기류를 모사한 위상판을 측정할 때에는 주로 샥하트만 파면 센서로 측정하게 된다. 이 때, 레이저 광원은 위상판을 거쳐 샥-하트만 파면 센서로 들어가게 되는데 위상판을 거치면서 레이저의 세기가 줄어들고, 이로 인해 샥-하트만 파면 센서가 위상판을 측정하지 못하는 경우가 발생한다. 본 논문에서는 난기류를 모사한 위상판 측정 시 레이저 출력 조절의 필요성과 레이저 출력이 측정된 파면에 어떤 영향을 미치는지를 알아본다. 프라이드 파라미터 r0이 상대적으로 낮은 위상판의 경우 레이저 출력으로 인해 10% 이상 r0이 변화하였다. r0이 상대적으로 높은 위상판의 경우 레이저 출력으로 인한 변화가 5% 미만으로 r0이 거의 변하지 않음을 보였다. 따라서 난기류가 심한 대기 상태를 모사한 위상판일수록 레이저 출력의 영향이 미미함을 알 수 있었다. 또한, 본 논문의 시스템을 기준으로 레이저 출력 5 mW 이상에서 난기류를 모사한 위상판을 측정할 수 있었다.

수질 매개변수 추정에 있어서 항공 초분광영상의 가용성 고찰 (Airborne Hyperspectral Imagery availability to estimate inland water quality parameter)

  • 김태우;신한섭;서용철
    • 대한원격탐사학회지
    • /
    • 제30권1호
    • /
    • pp.61-73
    • /
    • 2014
  • 본 연구는 항공 초분광영상을 사용한 수질추정 활용을 검토하고 한강일부분에 대해 가용한 측정자료를 이용하여 초분광영상 기반의 수질추정을 테스트하였다. 원격탐사에 의한 수질추정은 수체에 대한 downwelling과 수체 내에서의 산란과 반사에 대한 관측정보를 이용하는 방법과 원격탐사 센서에 도달하는 upwelling과 수질측정정보와의 선형적 회귀분석을 구하는 방법이 선호된다. 두 방법 모두 유의미한 결과를 도출하지만 수질정보나 산란정보 등 추정에 필요한 보조자료에 의한 영향이 더 클 것으로 판단되었다. 수질 추정 테스트는 팔당댐 하류에 위치한 한강의 일부분에 대해서 적용되었다. AISA eagle 초분광센서로 취득된 자료와 수질관측정보를 선형적 회귀분석을 통한 방법을 적용하였다. 기존 문헌에서 제시된 밴드조합에 대해서 회귀분석한 결과 유의미한 밴드조합으로 $-24.847+0.013L_{560}$의 회귀식을 얻었다 ($L_{560}$은 560 nm 파장에서의 radiance로 $R^2$=0.985). 다중분광영상을 이용했을 경우의 결과와 비교하기 위해서 spectral resampling을 통해 Landsat TM 영상을 생성하여 -55.932 + 33.881(TM, TM3)의 회귀식을 얻을 수 있었다(TM, TM3는 radiance로, $R^2$=0.968). 부유물질 농도는 수질측정지점에서 약 3.75 mg/l 이고, 초분광영상으로 추정된 농도는 약 3.65 mg/l, 시뮬레이션된 TM은 약 5.85 mg/l 로 다중분광영상을 이용했을 경우 과대 추정하는 경향을 보였다. 항공 초분광영상의 활용가치를 높이고 보다 정밀한 값을 추정하기 위해서 영상 전반에 걸친 sun glint 와 같은 영향을 최소화하기 위해 태양고도각을 고려하여 정교한 비행계획을 구성하고 체계적 전처리와 검 보정 체계를 갖출 필요가 있다고 사료된다. 일반적으로 적용된 방법에 따른 테스트로, 대기보정의 정밀성과 부족한 수질측정 샘플자료, 분광밴드의 검색, 적합한 선형회귀모델의 선택, 그리고 정량적 검증방법과 같은 몇 가지 문제점과 제약사항들을 발견할 수 있었다.

고 에너지 광자선의 표준측정법에 대한 선량 교정 프로그램 개발 (Development of a Dose Calibration Program for Various Dosimetry Protocols in High Energy Photon Beams)

  • 신동오;박성용;지영훈;이창건;서태석;권수일;안희경;강진오;홍성언
    • Radiation Oncology Journal
    • /
    • 제20권4호
    • /
    • pp.381-390
    • /
    • 2002
  • 목적 : 고 에너지 광자선에 대한 기준점에서의 물 흡수선량 계산을 절차상 또는 계산상의 오류를 피하기 위해 공기커마(혹은 조사선량) 교정정수에 토대를 두고 있는 IAEA TRS-277과 AAPM TG-21 및 최근 발표된 새로운 개념의 물 흡수선량 교정정수에 토대를 두고 있는 IAEA TRS-398과 AAPM TG-51 표준측정법에 기초한 고 에너지 광자선의 선량 교정 프로그램을 개발하고자 한다. 대상 및 방법 : 현재 국내외에서 널리 사용되고 있는 고 에너지 광자선에 대한 흡수선량 표준측정법은 IAEA TRS-277과 AAPM TG-21로서 공기커마(혹은 조사선량) 교정정수에 토대를 두고 있어 수식 체계가 복잡하고, 사용된 물리량에 대한 불확정도가 커서 선량측정의 정확성을 향상시키는데 한계가 있다. 최근 국제원자력기구와 미국의학물리학회에서는 새로운 개념의 물 흡수선량 교정정수에 토대를 두고 있는 IAEA TRS-398과 AAPM TG-51을 발표하였다. 개발된 네 종류의 선량 교정 프로그램은 이들 표준측정법에서 사용되고 있는 수식체계와 물리적인 매개변수를 엄격하게 적용하였고, 선량계에 대한 정보 및 물리적인 값에 대한 표와 그래프 값은 수치화하여 데이터베이스화하였다. 이들 프로그램은 윈도우 환경에서 사용이 용이하도록 비쥬얼 $C^{++}$ 언어를 사용하여 각각의 표준측정법에서 권고하고 있는 방법 및 절차에 따라 사용자의 편의성을 고려하여 개발하였다. 결과 : 네 종류의 표준측정법에 대하여 개발된 고 에너지 광자선에 대한 선량 교정 프로그램은 사용자가 병원에서 사용하고 있는 표준측정법을 선택하여 선량측정 절차에 따라 선량계, 선질 특성 및 측정 조건에 관한 정보와 측정 결과를 입력하고, 순차적으로 수행하도록 되어 있어 절차상 혹은 선량 계산에 있어서 사용자간의 오차 및 실수를 최소화할 수 있었다 또한 서로 다른 개념의 네 종류의 표준측정법에 대한 기준점에서의 선량값을 상호 비교할 수 있었다. 결론 : 이 프로그램은 이온함에 대한 정보와 물리적인 자료에 대한 표와 그래프 값들을 수식화하여 데이터베이스함으로써 수작업으로 각 프로토콜의 수행 절차상 혹은 사용자간의 발생할 수 있는 개인적인 실수 및 오차를 줄일 수 있었다. 또한 이 프로그램은 사용자 편의성을 고려하였고, 모든 보정계수와 물흡수선량을 정확하게 계산할 수 있기 때문에 각 표준측정법에 대한 주요한 차이점을 비교 분석할 수 있어 사용자가 적당한 표준측정법을 선택하여 수행하므로써 고 에너지 광자선 선량 교정에 이용시 매우 유익할 것으로 사료된다.

선천성 만곡족의 장기 수술적 치료 결과 (Long-Term Results of Surgical Treatment for the Idiopathic Clubfoot)

  • 김휘택;김인희;조윤재;안태영
    • 대한정형외과학회지
    • /
    • 제54권6호
    • /
    • pp.547-556
    • /
    • 2019
  • 목적: 도수조작 및 석고 교정을 통해 보존적 치료를 시행하였으나 잔여 혹은 재발변형으로 수술적 치료를 받았던 특발성 만곡족 환아의 장기 치료 결과를 분석하였다. 대상 및 방법: 방사선상 골 성숙이 완료된 32명의 환자 총 51예를 대상으로 후향적인 연구를 시행하였다. 최종 추시 시 평균 연령은 18.7세였다. 수술은 선택적 혹은 광범위 연부조직 유리술과 힘줄 연장, 힘줄 이전, 다양한 절골술 등을 시행하였다. 방사선적 분석은 최종 촬영한 체중부하 전후면 사진에서 거종골 각(talocalcaneal angle), 거골-제1 중족골 각(talo-first metatarsal angle), 측면 사진에서 측정한 거종골 각, 종골 경사각(calcaneal pitch)을 측정하여 성인정상 측정값과 비교하였다. 임상적 평가는 American Orthopaedic Foot and Ankle Society (AOFAS)의 ankle-hindfoot score와 midfoot score를 이용하여 우수(>85), 양호(71-85), 보통(56-70), 불량(<56)으로 평가하였다. 결과: 최종 방사선 계측치는 전후면 거종골 각에서 각각 41.2%, 전후면 거골과 제1 중족골 간의 각에서 90.2%, 측면 거종골 각에서 84.3%, 측면 종골 경사각에서 각각 61%가 정상범위에 포함되었다. AOFAS 평균은 ankle-hindfoot score에서 88.1±10.7점 midfoot score에서 86.7±11.5점이었다. 결론: 수술적 치료를 받았던 특발성 만곡족 환아의 장기 치료 결과에서 정상범위에 속하는 방사선 지표 비율은 43%-90%로 확인되었으며 임상적으로 AOFAS 평균값은 우수에 속하였다. 따라서 비수술적 교정 방법에 실패한 경우라도 변형의 요소를 보다 더 정확히 분석한 후 선택적 수술 방법을 통할 경우 임상적으로 만족스러운 결과를 가져올 수 있다.

반복적 대화식 통합 탄성파 속도분석 (An Iterative, Interactive and Unified Seismic Velocity Analysis)

  • 서상용;정부흥;장성형
    • 지구물리와물리탐사
    • /
    • 제2권1호
    • /
    • pp.26-32
    • /
    • 1999
  • 탄성파 속도분석법은 일괄식 속도분석법과 대화식 속도 분석등 두 가지가 있다. 일괄식 속도분석법에서는 각 속도 분석점마다 셈블런스 컨투어, 슈퍼게더 및 중합 패널등을 일괄 작성하여 도면화 시킨 후 분석자가 그 도면을 보고 속도 함수를 결정하는 방법이다. 과거 유행한 전산처리 소프트웨어들이 이 방법을 사용하고 있다. 그러나 이 방법은 도면 분석시 아주 많은 수작업이 필요하고 속도분석 결과도 정밀치 못하다는 단점이 있다. 최근에는 워크스테이션의 고속 그래픽 기능을 이용한 대화식 속도분석 기술이 개발되었다. 그런데 이들 프로그램은 기존 일괄식 속도분석법과 대동소이한 내용을 그래픽 화면으로 처리할 수 있도록 함으로써 종이 절약 외에는 특별히 나아진 것이 없다. 프로그램의 주 기능은 속도 스펙트럼에서 속도점 노드를 선택하는 것이며, 입력자료에 있을 수 있는 잡음을 제거하여 다시 속도 스펙트럼을 수정하는 기능은 없다. 잡음의 제거없이 계산한 부정확한 속도 스펙트럼을 이용해서 속도 함수를 선정한다면 정밀 속도분석은 불가능할 것이다. 방대한 탄성파 탐사자료에 대한 속도분석을 신속 정확하게 수행하기 위해서는 속도 분석과 밀접한 관련이 있는 전산처리 공정들 즉, 슈퍼게더 조립, 셈블런스 계산, 동보정, 뮤트, 중합등을 동시에 지원하는 통합된 반복적 대화식 속도분석 프로그램이 필요하다. 분석 구간의 속도와 뮤트함수를 변화시켰을 때 그로부터 얻어지는 셈블런스와 동보정 및 중합을 검토하고 이러한 수정과 검토를 신속히 반복할 수 있도록 함으로써 정확한 속도분석이 가능하기 때문이다. 여기에서는 속도분석을 신속 정확하게 수행하기 위해 속도 분석과 밀접한 관련이 있는 전산처리 공정들 즉, 슈퍼게더 조립, 셈블런스 계산, 동보정, 뮤트, 중합등을 동시에 지원하는 대화식 속도분석 프로그램 xva를 작성하였다. 대화식 속도분석에서는 분석 구간의 트레이스들을 고속으로 참조해야 하는데 이를 위해 간단한 트레이스 인덱스 파일을 설계하여 사용하였다. 직접파와 굴절파등 천부 잡음을 제거하기 위한 효과적인 수단인 뮤트 함수 영역 변환법을 새로 고안하였으며, 본 프로그램은 이 기법을 이용하고 있다. 본 영 역 변환법은 기존 알려진 역동보정법과 같이 정밀 전산처리가 가능할 뿐만 아니라 동보정과 역동보정시 발생하는 자료의 내삽 오차가 없으며 계산 시간이 크게 단축되기 때문에 정밀 대화식 속도 분석에 사용 가능하다. 프로그램 xva는 28개의 소스 파일로 구성된 패키지인데 줄 수는 12,029, 단어 수는 34,990, 글자 수는 304,073이다. 프로그램 xva는 X-Window와 Motif 환경하에서 작동한다. 프로그램 메뉴는 Motif 표준 스타일에 따라 작성하였는 바 그 사용법을 간략히 기술하였다. 본 프로그램이 완성됨으로 인하여 정밀 탄성파 속도 분석이 가능하게 되었고 그 결과 가스층의 존재 여부를 직접 확인할 수 있는 AVO(Amplitude Versus Offset)단면도등의 제작에 활용할 수 있었다.

  • PDF

N-13 암모니아와 양전자방출단층촬영 동적영상을 이용하여 심근혈류량을 정량화하는 새로운 방법 개발에 관한 연구 (A Refined Method for Quantification of Myocardial Blood Flow using N-13 Ammonia and Dynamic PET)

  • 김준영;이경한;김상은;최연성;주희경;김영진;김병태;최용
    • 대한핵의학회지
    • /
    • 제31권1호
    • /
    • pp.73-82
    • /
    • 1997
  • N-13 암모니아 PET 동적영상을 이용하여 심근혈류량을 측정할 때, 부분용적효과와 스필오버현상을 보정하는 새로운 방법을 고안하고, 이 방법을 이용하여 측정한 심근혈류량을 2구획모델만을 이용한 종래의 방법과 비교, 분석하여 새로운 방법에 대한 효율성과 정확성을 고찰한다. 9명의 관상동맥 환자에서 20mCi의 N-13 암모니아를 안정상태와 부하상태에서 주사한 후, PET 동적영상을 얻어 심근의 3부분(중격, 전면벽, 측면벽)과 좌심실방에 관심영역을 그려 시간-방사능곡선을 얻었다. 심근의 3부분에 대한 심근혈류량과 스필오버분획을 기하학적 관심영역 모델 개념을 2구획모델에 포함시킨 새로운 방법과 종래의 2구획모델 방법으로 각각 구하였다. 이때 관심영역의 위치에 따른 심근혈류량을 알아보기 위해서 심근 밖에서 좌심실 중심쪽으로 20개의 관심영역을 그리고, 위의 두 방법을 이용하여 심근혈류량을 구하였다. 종래의 N-13 암모니아 2구획모델과 기하학적 모델 개념을 2구획모델에 포함시킨 새로운 방법으로 모델 합치 곡선을 구하였다. 관심영역 위치 변화에 대한 심근혈류량의 값을 종래의 방법과 새로운 방법으로 구하고, 부분용적효과를 정확하게 보정하기 위하여 관심영역을 심내막쪽에 설정하였다. 총 108개의 심근관심영역에서 두 가지 방법을 이용하여 구한 심근혈류량 사이의 회귀곡선 기울기는 1.57, 상관계수는 0.88이었다. 그리고 같은 방법으로 얻어진 스필오버분획들도 선형적 상관관계(r=1.00, 기울기=0.98)가 있었다. 결론적으로 N-13 암모니아 PET 동적영상과 기하학적 모델 개념을 2구획모델에 포함시킨 새로운 방법을 이용하여 구한 심근혈류량은 종래의 방법보다 더 효율적이며 정확하게 정량화되어짐을 알 수 있었으며, 앞으로 임상환경에서 심근혈류량 정량분석연구에 유용할 것이라고 생각한다. 환자 군에 비해 유의하게 낮았다(p<0.01). 결론적으로, 흉통을 호소한 환자에서 심근관류 신티그라피상 정상인 경우에는 심장사건의 발생율이 낮음을 알 수 있었고, 특히 관동맥조영술에 관동맥병변이있는 환자와 없는 환자간에 1차 심장사건의 발생율에 유의한 차이가 없는 것을 고려하면, 심근관류 신티그라피가 정상소견을 보이는 흉통환자는 정상 관동맥조영술 소견을 보이는 환자에 준하여 치료하여도 좋을 것으로 사료되었다.리고 0.19, 0.40, 0.53, 0.61이었다. 실제 속도상수의 비 $k_3/k_4$에 대한 (BG-OCC)/OCC와 $R_A,\;R_v$간의 상관계수는 각각 0.983, 0.984, 0.999이었으며 그때의 기울기는 각각 1.76, 0.47, 1.25이었다. 결 론 : IPT 약역학은 시간이 흐름에 따라 혈류량의 변동에 비해 도파민 운반체량의 변동에 더욱 민감한 경향을 보였으며 $k_3/k_4$에 대한 (BG-OCC)/OCC, $R_A,\;R_v$의 결과간에 좋은 상관관계를 가졌다. 따라서 이러한 약역학 컴퓨터시뮬레이션이 SPECT 영상을 이용한 도파민 운반체 또는 수용체 정량분석을 최적화하는데 매우 유용할 것으로 생각된다.TEX>-CIT SPECT는 파킨슨병의 조기진단 및 진행 추적에 임상적으로 유용할 것으로 판단된다., SCC 4예, AC 1예)였으며, 11예 중 9예(81.8%)에서 방사선학적 검사결과와 Tc-99m MIBI섭취율의 변화가 일치하였다. 결론적으로, Tc-99m MIBI SPECT는 폐암병소의 국소화 및 방사선치료 효과의 판정에 어느정도 유용하리라 사료되었다.냈고 4명에서는 low CBD obstruction을 나타내었으며 후에 CBD stone, CBD carcinoma,

  • PDF