2011年3月11日，日本東部海上發生九級地震引發海嘯，造成東北沿岸地區農業上巨大的損失，共23,600公頃的農地受到損害。單在宮城縣和福島縣的海嘯受災區，稻田就占受災區60%以上，稻作為主要作物，因此，了解日本災後稻田受損和復育的情況對於農政單位將會是一項重要的議題。遙測影像提供多時序、大範圍的資料，而稻作在光譜上的反應比其他土地利用更為明顯，故本研究使用2010年(海嘯前一年)、2011年(海嘯發生年)、2012年(海嘯後一年)、2013年(海嘯後二年)、2014年(海嘯後三年)的MODIS影像，探討受災區稻田於海嘯發生前、後的分布情形以及災後的復育情況。本研究主要分為四個步驟：第一步驟，使用MODIS影像建立常態化差分植生指標(Normalized Difference Index，NDVI)的時間序列資料，並以小波轉換法(wavelet transform)濾除時序資料的雜訊，保留時序資料上有關稻作的資訊。第二步驟，使用支援向量機(Support Vector Machine, SVM)進行影像分類。第三步驟為精度評估。第四步驟，分析2010年至2014年間稻田受損和復育的情況。研究成果顯示，2010年至2014年的總體精度(overall accuracy)分別為88.0%、91.7%、89.5%、92.2%和90.8%，Kappa值分別為0.76、0.83、0.78、0.84和0.81。2011年海嘯過後，受災區有84.9%的稻田受損而無法耕種，2012年日本農民在沿岸地區開始陸續種植稻作，整體復育率為27.4%，2013年整體復育率為38.5%，2014年整體復育率為33.1%。將2010年至2014年的分類成果和統計資料相比，RMSE皆僅占總面積3%以下，顯示本研究能應用於災後稻田受損和復育之相關研究。

The 2011 Great East Japan Tsunami,happened on 11 March, 2011, destroyed 23,600 hectares of cultivation areas in northeastern coastal area of Japan, especially in Fukushima Prefecture and Miyagi Prefecture. More than 60% of inundation area is the area of rice fields. Therefore, monitoring the damaged rice fields and restoration of rice area after the tsunami is critical to provide agronomic planners with valuable information for the effective crop management strategies. Satellite imagery can provide multi-temporal and wide region data, and the spectral profile of rice is different from other land use types. This study used Moderate Resolution Imaging Spectroradiometer (MODIS) data from 2010 to 2014 to investigate the damaged rice fields and restoration of rice areas after the tsunami. The procedure of data processing consists of four steps: (1) data pre-processing to produce the time-series Normalized Difference Vegetation Index (NDVI) data and data filtering of the NDVI time-series data by wavelet transform; (2) paddy rice mapping using support vector machine (SVM); (3) accuracy assessment, and (4) area estimation of damaged rice fields and their restoration. The mapping results indicated the overall accuracy is 88.0%, 91.7%, 89.5%, 92.2% and 90.8% and the Kappa coefficient is 0.76, 0.83, 0.78, 0.84 and 0.81 from 2010 to 2014, respectively. After the tsunami (2011), 84.9 % of rice fields damaged by the tsunami is estimated. In 2012, the restoration rate is 27.4%. In 2013, the restoration rate is 38.5%. In 2014, the restoration rate is 33.1%. The comparison of RMSE divided by the total area in classification results and statistic data from 2010 to 2014 are all lower than 3%. It indicated this study can apply to the analysis of damaged rice fields and their restoration.

林世峻、莊智瑋、何世華、林昭遠，2008，"植生指標對影像分類準確度影響之研究"，水土保持學報，40(2)：181-193。
陳青、陳繼藩，2010，"應用經驗模態分解法分析多時期SPOT衛星影像上一年兩期之稻作"，國立中央大學土木工程學系研究所碩士論文。
黃思維、陳繼藩，2010，"應用經驗模態分解法分析多時期SPOT衛星影像上一年兩期之稻作"，國立中央大學土木工程學系研究所碩士論文。
Boschetti, M., Stroppiana, D., Brivio, P. A., et al. (2009). Multi-year monitoring of rice crop phenology through time series analysis of MODIS images. International Journal of Remote Sensing,30(18), 4643-4662.
Chagué-Goff, C., Przemyslaw, N., Henri, K. Y. W., et al. (2012). Environmental impact assessment of the 2011 Tohoku-oki tsunami on theSendai Plain. Sedimentary Geology, 282, 175-187.
Chen, C. F., Huang, S. W., Son, N. T., et al. (2011). Mapping double-cropped irrigated rice fields in Taiwanusing time-series Satellite PourI’Observation de la Terre data. Journal of Applied Remote Sensing, 5, 1-18.
Dewi, K. S., Ishak, H. I., Widyo, N. S., et al. (2010). Detecting Rice Phenology in Paddy Fields with Complex Cropping Pattern Using Time Series MODIS Data:A Case study of Northern Part of West Java-Indonesia.ITB Journal of Science, 42A(2), 91-106.
Foerster, S., Kaden, K. ,Foerster, M., et al.(2012). Crop type mapping using spectral-temporal profiles and phenological information.Computers and Electronics in Agriculture, 89 30-40.
Galford, G. L., Mustard, J. F., Melillo, J., et al. (2006). Wavelet analysis of MODIS time series to detect expansion and intensification of row-crop agriculture in Brazil. Remote Sensing of Environment, 112, 576-587.
Goto, K., Hashimoto, K., Sugawara, D., et al. (2013). Spatial thickness variability of the 2011 Tohoku-oki tsunami deposits along the coastline of Sendai Bay. Marine Geology, 358, 38-48
Gu, Y., Brown, J. F., Verdin, J. P., et al. (2007). A five-year analysis of MODIS NDVI and NDWI for grassland drought assessment over the central Great Plains of the United States. Geophysical Research Letters, 34(6).
Gusso, A., Formaggio, A. R., Rizzi, R., et al. (2012). Soybean crop area estimation by Modis/Evi data. Pesquisa Agropecuária Brasileira, 47(3).
Han, T., Li, Y., Han, H., et al. (2005). Land cover classification in Western Loess Plateau of China with MODIS imagery. Geoscience and Remote Seneing, 1, 202-205.
Huang, C., Davis, L. S., Townshend, J. R. G. (2002). An assessment of support vectormachines for land cover classification. International Journal of Remote Sensing, 23(4), 725-749.
Karahara, I. (2012). Earthquake and plant root researches. Plant Root, 6, 1-4.
Koltunov, A., Ustin, S. L., Asner, G., et al. (2009). Selective logging changes forestphenology in the Brazilian Amazon: Evidence from MODIS image time series analysis. RemoteSening of Environment, 113(11), 2431-2440.
Lv, T., & Liu C. (2010). Study on extraction of crop information using time-series MODIS data in the Chao Phraya Basin of Thailand. Advances in Space Research, 45(6), 775-784.
MAFF. (2011). Annual report on food, agriculture and rural areas in Japan Ministry of Agriculture, Forestry and Fisheries, Japan.
MAFF. (2011). Restoration and reconstruction of agriculture, forestry and fisheries of the Great East Japan Earthquake (in Japanese).
MAFF. (2013). FY2012 Annual Reporton Food, Agriculture and Rural Areas in Japan
Mori, N., Takahashi, T., Yasuda,T., et al. (2011).Survey of 2011 Tohoku earthquake tsunami inundation and run-up. Geophysical Research Letters, 38(7).
Murakami, T., Ogawa, S., Ishitsuka, K., et al.(2001).Crop discrimination withmulti-temporal SPOT/HRV data in the Saga Plains, Japan. International Journal of Remote Sensing, 22(7),1335-1348.
Odenweller, J. B., & Johnson, K. (1984). Crop identification using Landsat temporal spectral profiles. Remote Sensing of Environment, 14, 39-54.
Plaza, A., Martinez, P., Perez, R., et al. (2004). A new approach to mixed pixel classification of hyperspectral imagery based on extended morphological profiles. Pattern Recognition, 37, 1097-1116.
Richmond, B., Szczuciński, W., Chagué-Goff, C., et al. (2012). Erosion, deposition and landscape change on the Sendai coastal plain, Japan, resultingfrom the March 11, 2011 Tohoku-oki tsunami. Sedimentary Geology, 282, 27-39.
Sakamoto, T., Nguyen, N.V., Ohno, H., et al. (2005). Spatio-Temporal Distribution of Rice Phenology and Cropping Systems in The Mekong Delta With Special Reference to The Seasonal Water Flow of The Mekong and Bassac Rivers,Remote Sensing of Environment, 100, 1-16.
Sakamoto, T., Yokozawa, M., Toritani, H., et al. (2005). A Crop Phenology Detection Method using Time-Series MODIS data.Remote Sensing of Environment, 96, 366-374.
Son, N. T., Chen, C. F., Chen, C. R., et al. (2014). A Phenology-Based Classification of Time-Series MODIS Data for Rice Crop Monitoring in Mekong Delta, Vietnam.Remote Sensing, 6(1), 135-156.
Tomoko, M. N., Keitaro, T. (2013). Agricultural Implications of the Fukushima Nuclear Accident.Graduate School of Agricultural and Life Sciences.
Van Niel, T., & McVicar, T. (2004). Determining temporal windows for crop discrimination with remote sensing: a case study in south-eastern Australia. Computers and Electronics in Agriculture, 45, 91-108.
Vermote, E. F., Kotchenova, S. Y. (2008). MOD09 (Surface Reflectance) User’s Guide, MODIS Land Surface Reflectance Science Computing Facility.
Vinciková, H., Hais,M., Brom,J., et al. (2010). Use of remote sensing methods instudying agricultural landscapes-a review.Journal of LandscapeStudies, 3, 53-63.
Wardlow, B. D., & Egbert, S. L. (2008). Large-area crop mapping using time-series MODIS 250 m NDVI data: An assessment for the U.S. Central Great Plains. Remote Sensing of Environment, 112, 1096-1116.
Wardlow, B. D., & Egbert,S. L. (2003). A State-Level Comparative Analysis of the GAP and NLCD Land-Cover Data Sets. Photogrammetric Engineering and Remote Sensing, 69(12), 1387-1397.
Yang, C., Everitt, J. H., Murden, D. (2011). Evaluating high resolution SPOT 5 satellite imagery for crop identification.Computers and Electronics in Agriculture, 75, 347-354.
Yoshikawa, N., Obara, H., Ogasa, M., et al. (2014).Cs-137 in irrigation water and its effect on paddy fields in Japan after the Fukushima nuclear accident. Science of The Total Environment, 3(16).