问题:有大范围地区的雨量站数据,如何批量提取某个指定地区的雨量站点数据?
这是收集到的逐小时雨量数据格式,分别包含了站码、经纬度、高度、降雨量。
现在需要提取出shapefile范围内的站点。
1、首先,使用shaperead函数读取shapefile文件,识别出shp文件的经纬度信息。
%读取shape经纬度,需要注意的是shp文件坐标为WGS1984,这是用的是广州的shp文件
GZ=shaperead('GZ1984.shp');
下图的X、Y代表了shp的经度、纬度
2、然后提取逐小时数据站点的经纬度。
%读取站点经纬度,经纬度分别存储在第4行以下,第2、3列范围内。
Sx=shuju(4:length(shuju),2); %经度
Sy=shuju(4:length(shuju),3); %纬度
3、使用inpolygon函数判断站点经纬度是否在shp范围内。
%判断是否在shape范围内,其中GZ.X和GZ.Y为shp经度、纬度
GZs=inpolygon(Sx,Sy,GZ.X,GZ.Y);
此