Final Project: EPS 209Carolina Munoz Saez

IntroductionLast March 13th, 2011 Japan was affected for a M 8.9 offshore earthquake, which was followed for a tsunami

Tsunami waves were estimated to have reached 37.9 meters (124 ft) on a coastal city in Iwate Prefecture.The tsunami inundated a total area of approximately 500 km2 (181.5 sq mi) in Japan.

Earth movement seen from

spaceNASA satellites have snapped pictures of the devastating floods and fires in Japan from the deadly earthquake and tsunami

http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=Japan

Band 7-2-1 Combination: Vegetation, burned areas & water

Feb. 23, 2011

March 13, 2011

Pixel Size: 250m

Objective: Calculate area affected by tsunami

Inland (Band 721)Oceanic floor (True Color)

March 13, 2011 True-color imagery uses MODIS Bands 1, 4, and 3, corresponding to the red, green, and blue range of the light spectrum, respectively

Methodology:

Retrieve satellite image from a WMSCrop interest areaApply filtersGrain Analysis (regionprops)Mahalanobis distance

Results Inland

I1=imread ('Japan.2011054.aqua.721.250m.tif'); %cut image. M1 = [4278.5 3100.5 636 1716]G1 = imcrop(I1,M1) %Convert it to black & white:level1 = graythresh(G1); bw1 = im2bw(G1,level1);

Results Inland

graindata1=regionprops(cc1,'all');graindata2=regionprops(cc2,'all');grain_area1=[graindata1.Area];grain_area2=[graindata2.Area];Area_inland=grain_area2-grain_area1

Results Offshore: Mask

Results Offshore

% FOURTH STEP: Single grain dilatationSE = strel('disk', 150, 8);dsg1= imdilate(1-sbw1,SE);dsg2= imdilate(1-sbw2,SE);figure, imshow(dsg1)figure, imshow(dsg2)D=dsg1-(1-sbw1); %dilatation areafigure, imshow(D) %FIFTH: Ocean Floor: extract information for ocean floor on dilated grain MaskD2 = logical(D);D3 = cat(3, D2, D2, D2);A3 = immultiply(G3,D3);A4 = immultiply(G4,D3);figure;imshow(A3);figure;imshow(A4);

function [roiPix myStd] = getPixStd(rgbImg, filtSize)% read and check image, get size [nRows nCols nLyrs] = size(rgbImg);if (nLyrs ~= 3)error('Image is not RGB!');endnPix = nRows * nCols;% display original imagefigure;imshow(rgbImg);title('RGB image'); % smooth imagefprintf('Smoothing image ...\n');myFilt = ones(filtSize) ./ filtSize^2;smoothImg = imfilter(rgbImg, myFilt); …………………

% SIXTH: Ocean Floor: Aply Mahalanobis distance on the mask in both imagesfiltSize = 25;if (exist('myPatches.mat')) load('myPatches.mat');else % call read patch coderoiPixA = getPixStd(A3, filtSize);roiPixB = getPixStd(A4, filtSize); roiPix = cat(1, roiPixA, roiPixB);myStd = max(std(roiPix)); save('myPatches', 'roiPix', 'myStd');end

Final ResultsArea_inland= 9285 pixelArea_offshore = 22267 pixel

Image Pixel Size: 0.25 KmArea_inland= 580 Km2Area_offshore = 1392 Km2