研究者業績

劉 ウェン

リュウ ウェン  (Wen Liu)

基本情報

所属
千葉大学 大学院工学研究院 准教授
学位
博士(工)(千葉大学)

研究者番号
60733128
J-GLOBAL ID
201801019722087128
researchmap会員ID
B000345889

外部リンク

論文

 69
  • 山崎文雄, 劉ウェン
    日本地震工学会論文集 24(5) 309-322 2024年11月  査読有り
  • Kazuki Karimai, Wen Liu, Yoshihisa Maruyama
    Applied Sciences (Switzerland) 14(7) 2024年4月  
    Liquefaction is a significant challenge in the fields of earthquake risk assessment and soil dynamics, as it has the potential to cause extensive damage to buildings and infrastructure through ground failure. During the 2011 Great East Japan Earthquake, Urayasu City in the Chiba Prefecture experienced severe soil liquefaction, leading to evacuation losses due to the effect of the liquefaction on roads. Therefore, developing quantitative predictions of ground subsidence caused by liquefaction and understanding its contributing factors are imperative in preparing for potential future mega-earthquakes. This research is novel because previous research primarily focused on developing predictive models for determining the presence or absence of liquefaction, and there are few examples available of quantitative liquefaction magnitude after liquefaction has occurred. This research study extracts features from existing datasets and builds a predictive model, supplemented by factor analysis. Using the Cabinet Office of Japan’s Nankai Trough Megathrust Earthquake model, liquefaction-induced ground subsidence was designated as the dependent variable. A gradient-boosted decision-tree (GDBT) prediction model was then developed. Additionally, the Shapley additive explanations (SHAP) method was employed to analyze the contribution of each feature to the prediction results. The study found that the XGBoost model outperformed the LightGBM model in terms of predictive accuracy, with the predicted values closely aligned with the actual measurements, thereby proving its effectiveness in predicting ground subsidence due to liquefaction. Furthermore, it was demonstrated that liquefaction assessments, which were previously challenging, can now be interpreted using SHAP factors. This enables accountable wide-area prediction of liquefaction-induced ground subsidence.
  • 籠嶋 彩音, 劉 ウェン, 丸山 喜久, 堀江 啓
    土木学会論文集 79(13) n/a 2023年  
    2016年4月熊本地震では,熊本県熊本地方を震源とするMw6.2の地震が発生し,その約16時間後に同地域を震源とするMw7.0の地震が発生した.本研究では,地震による建物の被害状況を効率的にかつ安全に把握する方法として,航空レーザ測量データを深層学習することによって建物被害検出モデルの構築を試みた.本震前後に収集した航空レーザ測量データに対し,深層学習のアルゴリズムの一つである畳み込みニューラルネットワーク(CNN)を適用し,ネットワーク構成を変えながら最良のモデルの検討を行った.その結果,正答率が90%を超えるモデルを構築することができた.
  • 安江 崇志, 劉 ウェン, 丸山 喜久
    AI・データサイエンス論文集 4(3) 245-253 2023年  
    現在,日本の水道では年間2万件を超える漏水・破損事故が発生している.上水道管の漏水は,地上に流れ出す地上漏水と,地上には流れ出さず地下で流れている地下漏水の2種類に大別できる.地上漏水は人目に触れることから発見しやすいものの,地下漏水は漏水の状況を直接目視で確認できないため,早期発見のための技術開発が求められている.そこで本研究では,現在普及が進んでいるスマートメータを活用した水道管路のモニタリングを想定し,管網端部の水圧情報を使用した漏水位置予測に関する検討を行った.漏水シナリオや機械学習手法の異なる6つのモデルを構築し,その予測精度を比較した.水圧変化率,水圧変化量,管種情報を説明変数とし,LightGBMに基づき構築した漏水予測モデルが最も良好な結果を示した.

MISC

 78
  • Wen Liu, Fumio Yamazaki
    Joint Urban Remote Sensing Event 2013, JURSE 2013 33-36 2013年  
    Due to the remarkable improvements that have been made in radar sensors, high-resolution SAR images are now available, thus providing detailed ground surface information. In this study, a new method is developed to detect building heights automatically using 2D GIS data and a single high-resolution TerraSAR-X (TSX) intensity image. A building in a SAR image shows a layover from the actual position in the direction of the sensor, due to the side-looking nature of SAR. Since the length of the layover is proportional to the height of the building, it can be used to estimate the building height. To do this, we shift the building shape obtained from 2D GIS data in the direction of the sensor. The proposed method was tested on a TSX image of San Francisco, U.S.A. in HighSpot mode, using a resolution of about 1 m. Comparing the result with that obtained using Lidar data, the RMS error was found to be less than 3 m, which is about the height of one building story. © 2013 IEEE.
  • Wen Liu, Masashi Matsuoka, Fumio Yamazaki, Takashi Nonaka, Tadashi Sasagawa
    34th Asian Conference on Remote Sensing 2013, ACRS 2013 4 3744-3751 2013年  
    The Mw 6.1 earthquake affected Christchurch, New Zealand (NZ) on February 21, 2011. It caused widespread damage across the city, especially in the central area. Significant liquefaction occurred widely, which caused ground movement and destroyed lifelines and structures. In this study, the pre- And post-event ALOS/PALSAR and TerraSAR-X (TSX) data are used to detect the ground movements and liquefied areas. Firstly, the differential interferometric analysis (DInSAR) was applied to both the PALSAR and TSX data. The crustal movement caused by the earthquake was estimated by combining the two DInSAR results. Then the coherence value was used to detect the liquefied areas. Intensity images were also introduced to modify the detected result. Finally, the detected liquefied areas were verified by a liquefaction map which was made by field surveys. Copyright © (2013) by the Asian Association on Remote Sensing.
  • Kentaro Suzuki, Wen Liu, Miguel Estrada, Fumio Yamazaki
    34th Asian Conference on Remote Sensing 2013, ACRS 2013 2 1649-1656 2013年  
    In conducting damage assessment for scenario earthquakes in high seismic risk regions, building inventory data are required as well as building fragility functions and strong-motion distributions. But inventory data with the locations and characteristics of buildings are not so easy to construct, especially for developing countries. Hence in this study, an approach to construct building inventory data is sought as an alternative of cadastral data and field surveys. Using a high-resolution optical satellite image acquired by WorldView-2, this paper tries to develop building inventory data for earthquake damage assessment in Tacna, Peru. First, Pixel-based classification was carried out to examine basic land-cover and land-use of the urban area. Object-based building extraction was then conducted for three selected areas as an attempt to develop building inventory data.
  • Konomi Hara, Wen Liu, Fumio Yamazaki, Kentaro Suzuki, Yoshihisa Maruyama
    34th Asian Conference on Remote Sensing 2013, ACRS 2013 2 994-1001 2013年  
    Hyperspectral remote sensing makes it possible to obtain detailed spectral information of surface objects. Using airborne hyperspectral (HS) data acquired over Houston, Texas, USA, provided by the 2013 IEEE data fusion contest, the spectral reflectance characteristics of surface materials were investigated. A multispectral (MS) image acquired by WorldView-2 satellite was also introduced and it was compared with the HS image. A field measurement using a handheld spectroradiometer (EKO MS-720) was also carried out by the present authors. The irradiances of surface materials obtained by the measurement were also compared with the digital numbers of the 144 HS bands. Finally supervised classification was conducted for the HS and MS data and their results were discussed.
  • F. Yamazaki, W. Liu, E. Mas, S. Koshimura
    Safety, Reliability, Risk and Life-Cycle Performance of Structures and Infrastructures - Proceedings of the 11th International Conference on Structural Safety and Reliability, ICOSSAR 2013 5493-5498 2013年  
    Extraction of man-made structures is one of essential issues in image processing of remotely sensed data. Due to remarkable improvements in radar sensors, high-resolution Synthetic Aperture Radar (SAR) images are available, providing detail ground surface information. In this study, a new method is developed to detect building heights automatically from a 2D GIS data and a single high-resolution TerraSAR-X (TSX) intensity image.A building in a TSX image shows layover from the actual position to the direction of the sensor due to the side-looking nature of SAR. Since the length of layover on a ground-range SAR image is proportional to the building height, it can be used to estimate the building height. To do this, we shift the building footprint obtained from 2D GIS data in the direction of the sensor. The proposed method is tested on aTSX image of Lima, Peru in the High Spot mode with a resolution of about 1 m. The result showed a reasonable level of accuracy. © 2013 Taylor & Francis Group, London.
  • Wen Liu, Fumio Yamazaki, Masashi Matsuoka, Takashi Nonaka, Tadashi Sasagawa
    Conference Proceedings of 2013 Asia-Pacific Conference on Synthetic Aperture Radar, APSAR 2013 544-547 2013年  
    The Fukushima Earthquake with Mw7.1 occurred on April 11, 2011 was one of the most damaging induced events of the 2011 Tohoku earthquake with Mw9.0. It caused numerous fault scarps with a maximum displacement of 2.3 m. In this study, two methods were used to detect crustal movements from two different types of SAR images. Firstly, a differential interferometric analysis (DInSAR) was applied to pre- and post-event ALOS/PALSAR data. From the result of DInSAR, the trends of crustal movements in different areas could be grasped. Then an improved pixel-offset method was applied to pre- and post-event TerraSAR-X images. The two-dimensional movements were detected from the displacements of no-changed buildings. Finally, the detected results were compared with the observation data of GPS ground control stations. © 2013 IEICE.
  • Wen Liu, Fumio Yamazaki, Takashi Nonaka, Tadashi Sasagawa
    International Geoscience and Remote Sensing Symposium (IGARSS) 680-683 2013年  
    A method for capturing the two-dimensional (2D) surface movements from two temporal TerraSAR-X (TSX) intensity images has been proposed by the authors in previous research. However, it is impossible to detect the three-dimensional (3D) actual displacement from one pair of TSX images. Hence, three pairs of TSX images taken in ascending and descending paths were used to estimate 3D crustal movements in this study. First, the 2D crustal movements due to the 2011 Tohoku earthquake were detected from the three sets respectively. The relationship between the 3D actual displacement and 2D converted movement in SAR images was derived according to the observation model and shooting condition of the SAR sensor. Then the absolute 3D movements were estimated by the combination of the detected 2D movements that occurred within a short time interval. The results were verified by the GEONET observation records. © 2013 IEEE.
  • Fumio Yamazaki, Yoji Iwasaki, Wen Liu, Takashi Nonaka, Tadashi Sasagawa
    IMAGE AND SIGNAL PROCESSING FOR REMOTE SENSING XIX 8892 2013年  
    Building damage such as to side-walls or mid-story collapse is often overlooked in vertical optical images. Hence, in order to observe such building damage modes, high-resolution SAR images are introduced considering the side-looking nature of SAR. In the 2011 Tohoku, Japan, earthquake, a large number of buildings were collapsed or severely damaged due to repeated tsunamis. One of the important tsunami effects on buildings is that the damage is concentrated to their side-walls and lower stories. Thus this paper proposes the method to detect this kind damage from the change in layover areas in SAR intensity images. Multi-temporal TerraSAR-X images covering the Sendai-Shiogama Port were employed to detect building damage due to the tsunamis caused by the earthquake. The backscattering coefficients in layover areas of individual buildings were extracted and then, the average value in each layover area was calculated. The average value was seen to decrease in the post-event image due to the reduced backscatter from building side-walls. This example demonstrated the usefulness of high-resolution SAR intensity images to detect severe damage to building side-walls based on the changes of the backscattering coefficient in the layover areas.
  • Fumio Yamazaki, Jun Shimakage, Wen Liu, Takashi Nonaka, Tadashi Sasagawa
    2013 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM (IGARSS) 707-710 2013年  
    In this study, the flooded areas following the 2011 central Thailand flood were extracted using VNIR and TIR images of ASTER and ScanSAR-mode images of TerraSAR-X. The existence of water body was easily recognized for open spaces without trees and buildings from the NDVI value. The surface temperature was also found to be effective in detecting floods in a wide open space although it is limited by its coarse spatial resolution. The SAR intensity images were the most effective because water surfaces showed weak backscatter and they can be acquired at nighttime and under cloud-cover conditions. The extracted results were validated by a high-resolution optical satellite image.
  • 岩崎 洋志, 山崎 文雄, リュウ ウェン, 野中 崇志, 笹川 正
    日本地震工学会論文集 13(5) 5_18-5_32 2013年  
    衛星光学センサ画像を用いた建物の被害把握はこれまでに数多く行われてきたが、直下視の画像からは建物上面しか観察できないため、建物側面への被害や層崩壊などの被災パターンを判読することは困難であった。そこで本研究では、 高解像度の衛星搭載合成開口レーダ(SAR)で得られた地震前後の強度画像を用い、画像上の倒れ込み領域とレーダー影から個別建物の形状変化を把握し、被害抽出を試みた。仙台塩釜港周辺の建物と福島第一原子力発電所を対象として、多時期のTerraSAR-X画像を使用して、災害前後の倒れ込み領域とレーダー影内の後方散乱の変化を観察した。その結果、津波によって側面が損壊した建物や、爆発によって上部が吹き飛んだ原子炉建屋の被害を捉えることができた。
  • リュウ ウェン
    日本リモートセンシング学会誌 33(2) 158-158 2013年  
  • Wen Liu, Fumio Yamazaki
    IEEE GEOSCIENCE AND REMOTE SENSING LETTERS 10(1) 199-203 2013年1月  査読有り
    Significant crustal movements were caused by the 2011 Tohoku, Japan earthquake. A method for capturing the surface movements from pre- and postevent TerraSAR-X (TSX) intensity images is proposed in this letter. Because the shifts of unchanged buildings were considered as crustal movements in the two synthetic aperture radar images, we first extracted buildings from the pre- and postevent images using a segmentation approach. Then, the unchanged buildings were detected by matching the buildings in the pre- and postevent images at similar locations. Finally, the shifts were calculated by area-based matching. The method was tested on the TSX images covering the Sendai area. Compared with GPS observation records, the proposed method was found to be able to detect crustal movement at a subpixel level.
  • Wen Liu, Fumio Yamazaki, Takashi Nonaka, Tadashi Sasagawa
    33rd Asian Conference on Remote Sensing 2012, ACRS 2012 1 87-94 2012年  査読有り
    After the March 11, 2011 Tohoku, Japan earthquake (Mw9.0), numerous aftershocks and induced events occurred in the source region of the main shock and its outside areas. One of the most damaging induced earthquakes with Mw7.1 occurred in Fukushima prefecture on April 11, 2011. This event caused the numerous fault scarps, with a maximum displacement of 2.3 m. In this study, two temporal TerraSAR-X images were used to detect the crustal movements in Fukushima region due to these events by two different methods. Firstly, the interferometric (InSAR) analysis was carried out to detect the crustal movement to the radar sensor direction. Due to the temporal decorrelation, InSAR fringes could be obtained in a small part of the area. Secondly, the building-based pixel-offset method, proposed by the present authors, was applied to the SAR intensity images. 2D movements were detected for a half of the study area. Finally, the results were compared with the recorded displacements from GPS ground stations. The difference between the movements detected by the second method and the records was less than 0.3 m, showing high accuracy of our method.
  • W. Liu, F. Yamazaki, T. Nonaka, T. Sasagawa
    LAND SURFACE REMOTE SENSING 8524 2012年  査読有り
    The Tohoku earthquake on March 11, 2011 caused widespread devastation and significant crustal movements. According to the GPS Earth Observation Network System (GEONET) operated by Geospatial System Institution (GSI) of Japan, crustal movements with a maximum of 5.3 m to the horizontal direction (southeast) and a maximum of 1.2 m to the vertical direction (down) were observed over wide areas in the Tohoku (north-western) region of Japan. A method for capturing the two-dimensional (2D) surface movements from pre- and post-event TerraSAR-X (TSX) intensity images has been proposed by the present authors in our previous research. However, it is impossible to detect the three-dimensional (3D) actual displacement from one pair of TSX images. Hence, two pairs of pre- and post-event TSX images taken in ascending and descending paths respectively were used to detect 3D crustal movements in this study. First, two sets of 2D movements were detected by the authors' method. The relationship between the 3D actual displacement and 2D converted movement in SAR images was derived according to the observation model of the TSX sensor. Then the 3D movements were calculated from two sets of detected movements in a short time interval. The method was tested on the TSX images covering the Sendai area. Comparing with the GEONET observation records, the proposed method was found to be able to detect the 3D crustal movement at a sub-pixel level.
  • Liu, W, Yamazaki, F, Gokon, H, Koshimura, S
    International Geoscience and Remote Sensing Symposium (IGARSS) 6350523-4041 2012年  査読有り
  • Wen Liu, Fumio Yamazaki
    2011 Joint Urban Remote Sensing Event, JURSE 2011 - Proceedings 277-280 2011年  査読有り
    Monitoring urban growth and change is an important issue for urban planning and disaster management. In this study, three temporal TerraSAR-X images are used to monitor the urban changes. The difference and correlation coefficient between two images are calculated with a sliding window. Then a new factor that composites the difference and correlation coefficient is proposed to detect the changes. The Tokyo international airport expansion construction is used as a case study. The urban changes due to the progress of the construction are detected by the proposed approach. Then an aerial photograph taken during the construction and the visually detected results are used to verify the accuracy of the results. © 2011 IEEE.
  • Wen Liu, Fumio Yamazaki
    32nd Asian Conference on Remote Sensing 2011, ACRS 2011 1 585-590 2011年  査読有り
    The Mw 9.0 Tohoku earthquake on March 11, 2011, which occurred off the Pacific coast of the northeastern (Tohoku) Japan, caused gigantic tsunamis and brought vast devastation and a huge number of casualties. Various high-resolution satellites quickly captured the details of affected areas, and were used for emergency response. In this study, extraction of flooded areas was carried out using pre- and post-event synthetic aperture radar (SAR) images, which can observe the ground surface regardless of weather and sunlight conditions. Since the water surface generally shows very low backscattering intensity, the flooded areas could be extracted by the difference of SAR intensity between the pre- and post-event images from TerraSAR-X (X-band) and ALOS/PALSAR (L-band). Then, the characterizations of flooded areas were investigated, comparing the results from the X-band and L-band images. A pre-event ASTER DEM with 15m resolution was also employed to detect affected areas by tsunami. Finally, the results were compared with ground truth data to examine the accuracy of the proposed method.
  • Wen Liu, Fumio Yamazaki
    2011 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM (IGARSS) 2133-2136 2011年  査読有り
    Urban areas grow and change rapidly all over the world. Hence, regular and up-to-date information on urban changes is required for urban planning and disaster management. In this study, two temporal TerraSAR-X images are used to monitor urban changes. The study area is focused on a part of central Tokyo, Japan. Firstly, the changes between two images are checked by color composition. Then the difference and the correlation coefficient between the two images are calculated with a sliding window. A new factor that combines the difference and the correlation coefficient is proposed to detect changed areas. Finally, two high resolution optical images are introduced to verify the accuracy of the detection results.
  • リュウ ウェン, 山崎文雄
    日本写真測量学会学術講演会発表論文集 2010 139-140 2010年10月14日  
  • Fumio Yamazaki, Hisamitsu Inoue, Wen Liu
    2010 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM 4660-4663 2010年  査読有り
    SAR images obtained before and after a natural disaster are considered to be useful for emergency response due to its all-weather and sunlight-independent characteristics. Recently, the spatial resolutions of SAR systems have been improved significantly. In this paper, SAR intensity images acquired before and after the 2008 Iwate-Miyagi, Japan, earthquake from ALOS/PALSAR (L-band) and TerraSAR-X (X-band) are employed to investigate the radar backscattering characteristics for various acquisition and surface conditions. The spatial resolution, radar frequency, flight path, and incidence angle were shown to affect SAR backscattering echo, depending on surface materials and roughness. It is also observed that the difference of the backscattering coefficients at the pre- and post-event times gets large and their correlation coefficient becomes small at the locations of landslides and slope failures.
  • Wen Liu, Fumio Yamazaki
    2010 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM 2206-2209 2010年  査読有り
    Shadows in remote sensing images often result in problems for many applications such as land-cover classification, change detection, and damage detection in disasters. Due to these reasons, it is very useful if the radiance of shadowed areas is corrected to the same radiance as shadow-free areas. In this study, a shadow detection and correction method is proposed. Shadowed areas are detected by object-based classification, using brightness values and a neighbor relationship. Then the detected shadowed areas are corrected by a liner function to produce a shadow-free image. The shadowed areas with different darkness are corrected with different ratios to improve the accuracy of the result. The spectral characteristics of sunlit and shadowed areas in several QuickBird images were studied and then the shadow-free radiance was obtained.
  • リュウ ウェン, 山崎文雄
    日本写真測量学会学術講演会発表論文集 2009 69-70 2009年6月18日  
  • Fumio Yamazaki, Wen Liu, Makiko Takasaki
    2009 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM, VOLS 1-5 4 2806-2809 2009年  査読有り
    The effects of shadow in remote sensing imagery are investigated The measurement of radiance in sunlit and shadowed areas was carried out to investigate the spectral characteristics of sunlight Based on this observation, It is found that the radiance ratio (shadow/sunlit) increases as the sunlight gets weaker and the ratio is dependent on the wavelength of sunlight. The darkness of shadow is also found to vary depending on the surrounding condition Thus the condition to restore a shadow-free image depends on the spectral bands and the location even in one Image. A QuickBird image is then introduced and the spectral characteristics of sunlit and shadowed areas are investigated. Based on these observations, a method to detect shadowed areas and restore the shadow-free radiance for the multi-spectral bands is proposed. The effectiveness of the shadow correction method is demonstrated for the Quick Bird image.
  • Wen Liu, Fumio Yamazaki
    2009 JOINT URBAN REMOTE SENSING EVENT, VOLS 1-3 1159-1164 2009年  査読有り
    A new method is developed to extract moving vehicles and subsequently detect their speeds from a pair of QuickBird (QB) panchromatic (PAN) and multi-spectral (MS) images automatically. Since PAN and MS sensors of QB have a slight time lag, about 0.2 seconds, the speed of moving vehicles can be detected by the movement between PAN and MS images in the time lag. From a PAN image with 0.6m resolution, vehicles can be extracted by an object-based approach. But it is difficult to extract the accurate position of vehicles from a MS image with 2.4m resolution. Thus an area correlation method is proposed to estimate the location of vehicles from MS images in a sub-pixel level. Using the results of the vehicle extraction, the speed of moving vehicles can be detected. The approach is tested on several parts of the QB image covering the central Tokyo, Japan, and the accuracy of the result is demonstrated.
  • リュウ ウェン, 山崎文雄
    日本地震工学会大会梗概集 6th 348-349 2008年11月  
  • Fumio Yamazaki, Wen Liu, T. Thuy Vu
    International Geoscience and Remote Sensing Symposium (IGARSS) 3(1) 1334-1337 2008年  査読有り
    A new object-based method is developed to extract the moving vehicles and subsequently detect their speeds from two consecutive digital aerial images automatically. Several parameters of gray values and sizes are examined to classify the objects in the image. The vehicles and their associated shadows can be discriminated by removing big objects such as roads. To detect the speed, firstly the vehicles and shadows are extracted from the two images. The corresponding vehicles from these images are linked based on the order, size, and their distance within a threshold. Finally, using the distance between the corresponding vehicles and the time lag between the two images, the moving speed can be detected. Our test shows a promising result of detecting the moving vehicles' speeds. Further development will employ the proposed method for a pair of QuickBird panchromatic and multi-spectral images, which are at a coarser spatial resolution. © 2008 IEEE.
  • Wen Liu, Fumio Yamazaki
    29th Asian Conference on Remote Sensing 2008, ACRS 2008 2 693-698 2008年  査読有り
    Damage detection is carried out from ALOS optical images that captured the stricken areas by the Sichuan, China earthquake, which occurred on May 12, 2008. Since landslides occurred extensively in this earthquake, we try to extract the landslide areas comparing the pre- and post- event images of ALOS/AVNIR-2. First, a level-slice method is used to extract vegetated areas and water areas. Then, the bare ground areas are detected from the pre- and post- event images by a pixel-based classification, after masking vegetation and water. Then the difference of the bare ground between two images is considered as landslides. The rise of water-level in the river is also detected in the areas where landslides dammed off the river flow. The digital elevation model from SRTM is also employed in this study to investigate the relationship between the slope angle and the occurrence of landslide.
  • Liu W, Yamazaki F, Vu T.T, Maruyama Y
    28th Asian Conference on Remote Sensing 2007, ACRS 2007 3 1577-1582 2007年  査読有り

書籍等出版物

 2

講演・口頭発表等

 57

担当経験のある科目(授業)

 5

共同研究・競争的資金等の研究課題

 12