抓取地震台网近一年的地震数据，并自定义筛选保留(java)

public void getdata(Integer siteId) throws Exception {
		DecimalFormat df = new DecimalFormat("0.00");
		String str = "http://www.ceic.ac.cn/speedsearch?time=5&&page=";//地震台数据地址
		//添加字符输出缓冲流
		BufferedWriter bw = new BufferedWriter(new FileWriter("F:\\123.txt"));//临时存储到本地
		//定义统一资源标识符
		URL url = null;
		//定义字符缓冲类
		StringBuffer sb = new StringBuffer();
		//定义字节输入流
		InputStream is = null;
		//定义字符输入流，可修改编码表
		InputStreamReader isr =null;
		//定义字符输入缓冲流
		BufferedReader bf =null;
		//定义文件类
		Document document = null;
		Elements element =null;
		url = new URL(str+""+1);
		is = url.openStream();
		isr = new InputStreamReader(is,"UTF-8");//使用"UTF-8"编码
		bf = new BufferedReader(isr);
		String line =null;
		while((line = bf.readLine())!= null){
			//System.out.println(line);
			sb.append(line+"\n");  //数据写入缓冲流中
		}
		//使用Jsoup进行html的解析
		document = Jsoup.parse(sb.toString()) ;
		element = document.select("tr");//网页节点
		for(int i=1;i<element.size();i++){
			bw.write(element.get(i).text()); //将数据写入文件中
			bw.newLine(); //写换行
			bw.flush();   //刷新
			//System.out.println();
		}
		Element ele = document.getElementById("pagenum");//网页节点
		//目的为了读取翻页网页的数量，即该数据在很多相似的网页上
		String string = ele.text().toString();
		String[] arraystr = string.split(":");//分割字符串
		Integer isign = new Integer(arraystr[2]);//字符变数组
		int intsign = isign.intValue();
		bf.close();
		sb.delete(0, sb.length());
		//利用上面读到的翻页数量，进行网页的循环读取
		for(int sign=2; sign<=intsign; sign++){
			url = new URL(str+""+sign);
			is = url.openStream();
			isr = new InputStreamReader(is,"UTF-8");
			bf = new BufferedReader(isr);
			line =null;
			while((line = bf.readLine())!= null){
				sb.append(line+"\n");
			}
			document = Jsoup.parse(sb.toString()) ;
			element =document.select("tr");
			for(int i=1;i<element.size();i++){
				bw.write(element.get(i).text());
				bw.newLine();
				bw.flush();
				String text = element.get(i).text().trim();
				String[] split = text.split(" ");

				//如果地震源范围300公里里包含监测点
				double lon = 0.0; //需要判断的位置经度
				double lat = 0.0; //需要判断的位置监测点纬度
				
					//计算震源地与需要判断的位置之间的距离(m)
				double v = gpsDistanceUtil.GetDistance(lon, lat, Double.parseDouble(split[4]), Double.parseDouble(split[3]));//返回m 

				if(v<=(300*1000)){//范围定义300公里范围内
					//如果数据库存在相同记录，不再记录
					
					//*******************以下代码可根据实际业务进行调整*******************
					QueryWrapper queryWrapper=new QueryWrapper();
					queryWrapper.eq("time",String.valueOf(DateUtil.StringToTimestamp(split[1]+" "+split[2])));
					queryWrapper.eq("address",split[6]);
					List<TEarthquakeReal> list = tEarthquakeRealServicep.list(queryWrapper);
					if(list.size()==0){//数据库不存在
						
						//在震源范围内进行入库
						TEarthquakeReal tEarthquakeReal=new TEarthquakeReal();
						tEarthquakeReal.setGrade(Double.parseDouble(split[0]));
						tEarthquakeReal.setTime(String.valueOf(DateUtil.StringToTimestamp(split[1]+" "+split[2])));
						tEarthquakeReal.setLat(Double.valueOf(split[3]));
						tEarthquakeReal.setLon(Double.valueOf(split[4]));
						tEarthquakeReal.setDepth(Double.valueOf(split[5]));
						tEarthquakeReal.setAddress(split[6]);
						tEarthquakeReal.setSiteId(siteId);
						tEarthquakeReal.setDistance(Double.valueOf(df.format(v/1000)));
						tEarthquakeRealServicep.save(tEarthquakeReal);
						//*************************************************************
					}
				}

			}
			bf.close();
			sb.delete(0, sb.length());

		}
		bw.close();
		//System.out.println("已全部加载完:震级(M)\t发震时刻(UTC+8)\t纬度(°)\t经度(°)\t深度(千米)\t参考位置");

	}

gpsDistanceUtil类：计算两个经纬度之间的距离

public class GPSDistanceUtil {

	/**
	 * 默认地球半径
	 */
	private static double EARTH_RADIUS = 6371000;//赤道半径(单位m)

	/**
	 * 转化为弧度(rad)
	 * */
	private static double rad(double d)
	{
		return d * Math.PI / 180.0;
	}
	/** 方式一
	 * @param lon1 第一点的经度
	 * @param lat1 第一点的纬度
	 * @param lon2 第二点的经度
	 * @param lat2 第二点的纬度
	 * @return 返回的距离，单位m
	 * */
	public double GetDistance(double lon1,double lat1,double lon2, double lat2) {
		double radLat1 = rad(lat1);
		double radLat2 = rad(lat2);
		double a = radLat1 - radLat2;
		double b = rad(lon1) - rad(lon2);
		double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) + Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b / 2), 2)));
		s = s * EARTH_RADIUS;
		s = Math.round(s * 10000) / 10000;
		return s;
	}

	/**
	 * 计算中心经纬度与目标经纬度的距离（米）  方式二
	 *
	 * @param centerLon
	 *            中心精度
	 * @param centerLat
	 *            中心纬度
	 * @param targetLon
	 *            需要计算的经度
	 * @param targetLat
	 *            需要计算的纬度
	 * @return 米
	 */
	private static double distance(double centerLon, double centerLat, double targetLon, double targetLat) {

		double jl_jd = 102834.74258026089786013677476285;// 每经度单位米;
		double jl_wd = 111712.69150641055729984301412873;// 每纬度单位米;
		double b = Math.abs((centerLat - targetLat) * jl_jd);
		double a = Math.abs((centerLon - targetLon) * jl_wd);
		return Math.sqrt((a * a + b * b));
	}

//	public static void main(String []args){
//		double lon = 116.510958; //当前经度
//		double lat = 39.90786; //当前纬度
//		double lon1=116.510842;
//		double lat1=39.90777;
//		double dist;
//		double dist2;
//		dist=GetDistance(lon, lat, lon1, lat1);
//		dist2=distance(lon, lat, lon1, lat1);
//		System.out.println("1.两点相距：" + dist + " 米");
//		System.out.println("2.两点相距：" + dist2 + " 米");
//	}
}