clc
clear all
close all
%读取excel中数据
data = xlsread('D:/metero.xlsx', 'sheet3', 'A2:G579');
%获取精度纬度
lon= data(:, 6);
lat = data(:, 5);
%获取降水天数的斜率
ppt = 10.*data(:, 2);
%先创建图像窗口句柄
figure1=figure;
%创建axes句柄
%展示中国边界
geoshow('D:/LST/china_line.shp', 'color', 'black');
%设置坐标轴
%设置坐标轴标签
xlabel('Longitude(°E)', 'fontsize', 12, 'fontweight', 'b');
ylabel('Latitude(°N)', 'fontsize', 12, 'fontweight', 'b');
%显示四周边框
box on;
hold on;
%将流域边界加进去使用geoshow函数
geoshow('D:/LST/river_border.shp', 'color', 'blue');
%MarkerEdgeColor表示标记边缘颜色,将其设为白色,即[1.0 1.0 1.0]
scatter(lon, lat, 50, ppt, 'filled&