从表数据创建等值线图的地图

读取数据到表

一个逗号分隔的文本文件读入一个表。文件包含农田统计数据对于我们每一个州和哥伦比亚特区,基于国家农业统计服务(NASS)农场和土地2019年农场总结。有关统计数据的更多信息,请参阅[1]。显示表的第一行。

农场= readtable (“farmland.txt”,“TextType”,“字符串”);:农场(1)

ans =1×4表StateName NumberOfFarms AcresInFarms AverageFarmSize _____ _________________ _______售予“阿拉巴马”38800 8300 214

从表中创建一个地图,你也需要空间信息。阅读包含多边形形状的shapefile美国每一个州和哥伦比亚特区的地理空间表。一个地理空间表有一个形状为每一行包含形状信息的变量。显示地理空间表的第一行。

州= readgeotable (“usastatelo.shp”);州(1:)

ans =表1×5形状名称LabelLat LabelLon PopDens2000 _______ _____说___________ geopolyshape 32.283 -86.921 87.6“阿拉巴马”

的形状表包含的变量geopolyshape对象,多边形的地理坐标。显示的属性geopolyshape对象在第一行。

states.Shape (1)

ans = geopolyshape属性:NumRegions: 1 NumHoles: 0几何:“多边形”CoordinateSystemType:“地理”GeographicCRS: [1×1 geocrs]

连接表

表的农田统计确定了哥伦比亚特区华盛顿特区,它不匹配的相应行表的多边形形状。改变农田的名字表中的数据表中,这样它的名称匹配的多边形形状。

dcRow =农场。StateName = =“华盛顿特区”;farms.StateName (dcRow) =“特区”;

加入表,使用国家名称变量的名字和StateName的关键变量。加入表的第一行显示。

statesFarms = outerjoin(州、农场、“LeftKey”,“名称”,“RightKey”,“StateName”);:statesFarms (1)

ans =1×9表形状名称LabelLat LabelLon PopDens2000 StateName NumberOfFarms AcresInFarms AverageFarmSize _______ _____说___________ _____ _________________ _______售予geopolyshape“阿拉巴马”32.283 -86.921 87.6 38800 8300 214“阿拉巴马”

关于如何结合矢量数据使用的更多信息表连接,看到的结合矢量数据表连接。

地图创建等值线图

创建一个等值线图的地图,显示了每个州的农场的平均尺寸相接的我们。

图行= statesFarms。名字~ =“夏威夷”& statesFarms。名字~ =“阿拉斯加”;:statesFarmsConus = statesFarms(行);geoplot (statesFarmsConus ColorVariable =“AverageFarmSize”);

改变技术。添加一个标题和标签colorbar。

geobasemapdarkwater标题(“平均每2019年美国国家农场大小”)cb = colorbar;cb.Label。字符串=“英亩大小”;

改变colormap这样,随着平均农场规模的增加,美国的颜色从黄色转变为绿色。

提出=翻转(夏季);colormap城市规划机构(cmap)

出口地图

导出映射到一个PNG文件。

甘氨胆酸gx =;exportgraphics (gx“AverageFarmSizeForUS.png”)

引用

[1]国家农业统计服务。“农场,农场的土地,平均农场规模-州和美国:2018 - 2019。”2019年农场农场和土地的总结6。美国农业部,2020年2月。https://www.nass.usda.gov/Publications/Todays_Reports/reports/fnlo0220.pdf。