最近地震地图
我们想读的位置最近的地震地质勘探局网站和情节在一个交互式地图。
内容
从美国地质调查局读取位置
加载数据的实时数据处理。阅读所有的4.5级地震在过去一个月,大。在这里找到的信息数据:http://earthquake.usgs.gov/earthquakes/feed/v1.0/glossary.php
找到过去30天,4.5或更高版本。http://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/4.5_month.geojson
可以寻找上周如果我们优先。http://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/4.5_week.geojson
选择= weboptions (“超时”10);quakeDataJSON = webread (“http://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/4.5_month.geojson”、选择);
提取每个地震信息
quakeDataInfo = [quakeDataJSON.features.properties];quakeDataLocation = [quakeDataJSON.features.geometry];
转换表结构
开放quakeTable看看它的样子
quakeTable = struct2table (quakeDataInfo);:quakeTable (1:5)
ans =杂志的地方时间更新tz url细节感到cdi mmi警戒状态海啸团体净代码ids数据源类型望远镜dmin rms差距magType类型标题____ ______________________________________ __________ __________ ____ _____________________________________________________________ ___________________________________________________________________________ ___ _____ ___ _____ __________ _________ ___ _____ _____ * * * ___________________________________________________________________________ ___累积___ _____ _______ ______________________________________________ 5.2 Tairua东北偏东约122公里,新西兰的12 1.4247 1.4247 e + e + 12 720“http://earthquake.usgs.gov/earthquakes/eventpage/usc000tsfu”“http://earthquake.usgs.gov/earthquakes/feed/v1.0/detail/usc000tsfu.geo..”。[0][1][][]“综述”0 416“我们”“c000tsfu”“usc000tsfu,”“,我们,”帽,戴菲现实,geoserve,附近城市,起源、相位数据,我们……”[][2.622]1.17 [124]mb‘地震’‘M 5.2 - Tairua东北偏东约122公里,新西兰“4.9”Taron之西南偏南约134公里,巴布亚新几内亚的1.4247 e + 12 1.4247 e + 12 600 http://earthquake.usgs.gov/earthquakes/eventpage/usc000tsfg“http://earthquake.usgs.gov/earthquakes/feed/v1.0/detail/usc000tsfg.geo..。”[0][1][][]”了“0 369”我们“c000tsfg ' ', usc000tsfg, ' ',, ' ',帽,戴菲,现实,geoserve,附近城市,起源、相位数据,我们……”[][1.559]0.66 [90]mb‘地震’‘M 4.9 - Taron之西南偏南约134公里,巴布亚新几内亚“4.6”Oami之西南偏南约3公里,日本“1.4247 e + 12 1.4247 e + 12 540 http://earthquake.usgs.gov/earthquakes/eventpage/usc000tsf4“http://earthquake.usgs.gov/earthquakes/feed/v1.0/detail/usc000tsf4.geo..。”[26][3.1][][]”了“0 334”我们“c000tsf4 ' ', usc000tsf4, ' ',, ',帽,戴菲,现实,geoserve,附近城市,起源、相位数据,我们……”[][0.748]1.1 [115]mb‘地震’‘M Oami之西南偏南约4.6 - 3公里,日本“4.5”布莱尔港西南119公里,印度的1.4247 e + 12 1.4247 e + 12 360 http://earthquake.usgs.gov/earthquakes/eventpage/usc000tses“http://earthquake.usgs.gov/earthquakes/feed/v1.0/detail/usc000tses.geo..。”[0][1][][]”了“0 312”我们“c000tses ' ', usc000tses, ' ',, ' ',帽,戴菲,现实,geoserve,附近城市,起源、相位数据,我们……”[][1.057]1.05 [120]mb‘地震’‘M 4.5 -布莱尔港西南119公里,印度“4.5”41公里ESE的磐,日本“1.4247 e + 12 1.4247 e + 12 540 http://earthquake.usgs.gov/earthquakes/eventpage/usc000tsdx“http://earthquake.usgs.gov/earthquakes/feed/v1.0/detail/usc000tsdx.geo..。“[][][][]”了“0 312”我们“c000tsdx ' ', usc000tsdx, ' ',, ',帽,现实,geoserve,附近城市,起源、相位数据,tectonic-su……“[][1.047]1.28[205]的4.5 mb‘地震’‘M - 41公里ESE的磐,日本的
表中是什么?
quakeTable.Properties.VariableNames
ans =列1到7“杂志”“地方”“时间”“更新”“tz”“url”“细节”列8到14”感到“cdi的mmi的“警报”“状态”“海啸”“团体”列15至21“净”“代码”“id”“源”“类型”“望远镜”“dmin”列22到26 rms的“差距”“magType”“类型”“标题”
添加位置信息表
首先让我们得到地震坐标。然后将它们附加到信息收集表。
eqcoordinates = [quakeDataLocation.coordinates] ';quakeTable。朗= eqcoordinates (: 1);quakeTable。Lat = eqcoordinates (:, 2);quakeTable.depth = eqcoordinates (: 3);quakeTable.Properties.VariableNames
ans =列1到7“杂志”“地方”“时间”“更新”“tz”“url”“细节”列8到14”感到“cdi的mmi的“警报”“状态”“海啸”“团体”列15至21“净”“代码”“id”“源”“类型”“望远镜”“dmin”列22到28 rms的“差距”“magType”“类型”“标题”“朗”“纬度”列29“深度”
按震级表排序
首先看看无序表的一部分
quakeTable.mag (1:5)
ans = 5.2 4.9 4.6 4.5 - 4.5
它是更难保持一起像数据结构,而不是表。
quakeTable = sortrows (quakeTable,“杂志”,“提升”);quakeTable.mag (end-4:结束)
ans = 6.3 6.4 6.7 6.7 - 7.1
找到最大的地震
quakeTable(最终,{“地方”,“纬度”,“朗”,“深度”})
ans Lat =地方经度深度_____________________________ ______和_____的大西洋中脊的北部52.635 -31.884 16.4
首先找到范围的地震情况
我们可以选择销maxmag我们最大的期待和minmag 4.5我们可以把每月的地图与颜色协调。
minmag = min (quakeTable.mag);maxmag = max (quakeTable.mag);
情节上webmap地震位置
发现纬度和经度范围图标颜色大小
厘米= parula (10);iconColor =厘米(装天花板(1 + 9 * (quakeTable.mag-minmag) / (maxmag-minmag)),:);
把名称的位置放在hoverover标签
名称= quakeTable.title;
转换quakeTable geopint向量
quakePoints = geopoint (table2struct (quakeTable));
地图与海洋技术符合一些断裂区
webmap (“海洋技术”)wmmarker (quakePoints“OverlayName”,“地震点”,…“FeatureName”、名称、“颜色”iconColor,“最适合的”、假);wmzoom (1) snapnow
负载在板块边界
我这里可以说明文件导入。板块边界的数据参考http://geoscience.wisc.edu/查克/ MORVEL / PltBoundaries.html引用:Argus, d F。、戈登·r·G。,and DeMets, C., 2011. Geologically current motion of 56 plates relative to the no-net-rotation reference frame, Geochemistry, Geophysics, Geosystems, accepted for publication, September, 2011.
(纬度、经度)= importPlates (“All_boundaries.txt”);海岸=负载(“海岸”);图worldmap世界setm (gca),“mlabelparallel”,-90,“mlabellocation”,90)plotm (coast.lat coast.long,“颜色”,“k”)plotm(纬度、经度、“线宽”,2)
找到第一个板
印第安纳州=找到(isnan (lat), 1,“第一”)plotm(纬度(1:印第安纳州)、经度(1:印第安纳州),“颜色”,“红色”,“线宽”3)
印第安纳州= 829
使一系列geopoints板块边界
边界= geopoint(纬度、经度);
画地图板块边界
地图中心的最大地震的经度。
wmcenter (0, quakeTable.Lon(结束))wmline(范围,“FeatureName”,“板块边界”,“颜色”,“米”,“最适合的”假)wmzoom (1) snapnow
