clc;
clear;
close all;
Chprovince_0 = shaperead('省界111.shp', 'UseGeoCoords', true);%读入31省界数据,建立结构体
Chprovince_1 = shaperead('台湾和香港.shp', 'UseGeoCoords', true);
Chprovince_2 = shaperead('九段线.shp', 'UseGeoCoords', true)
Chprovince_3 = shaperead('南海诸岛及其它岛屿.shp', 'UseGeoCoords', true);
T = struct2table(Chprovince_0);
T = struct2table(Chprovince_2);
temis = xlsread('PROVINCE_MAPPING_second.xlsx',1,'E36:E78');
for i=1:31
Chprovince(i).Geometry = Chprovince_0(i).Geometry;
Chprovince(i).BoundingBox = Chprovince_0(i).BoundingBox;
Chprovince(i).Lon = Chprovince_0(i).Lon;
Chprovince(i).Lat = Chprovince_0(i).Lat;
Chprovince(i).name = Chprovince_0(i).NAME;
Chprovince(i).EmisInten2017 = tem