| | |
| | | import java.awt.*; |
| | | import java.awt.image.BufferedImage; |
| | | import java.io.*; |
| | | import java.math.BigDecimal; |
| | | import java.math.RoundingMode; |
| | | import java.text.ParseException; |
| | | import java.text.SimpleDateFormat; |
| | | import java.util.*; |
| | | import java.util.List; |
| | | |
| | |
| | | @Resource |
| | | PropertiesConfig config; |
| | | |
| | | public final static double MAX_X_OFFSET = 0.002; |
| | | public final static double MAX_X_OFFSET = 0; |
| | | |
| | | public void process(DataPo data) throws Exception { |
| | | String basePath = config.getInPath() + File.separator + data.getInPath() + File.separator; |
| | |
| | | copeWater(dto, layer); |
| | | copeFlow(dto, layer); |
| | | copeLayerJson(dto, layer); |
| | | copeRainFallJson(dto, layer); |
| | | } finally { |
| | | File dir = new File(dto.getTemp()); |
| | | if (dir.exists()) { |
| | |
| | | Dataset ds = null; |
| | | try { |
| | | ds = gdal.Open(dto.getTerrainFile(), gdalconst.GA_ReadOnly); |
| | | if (null == ds || 0 == ds.getRasterCount() || null == ds.GetSpatialRef()) return; |
| | | if (null == ds || 0 == ds.getRasterCount()) return; |
| | | |
| | | setTerrainInfo(ds, layer); |
| | | setWaterInfo(dto, layer); |
| | |
| | | } |
| | | } |
| | | |
| | | /** |
| | | * 重采样 |
| | | */ |
| | | private static void Resample(Dataset ds, String dest, int width, int height, LayerDto layer) { |
| | | Vector<String> vector = new Vector<>(); |
| | | vector.add("-t_srs"); |
| | | vector.add("EPSG:" + 4326); |
| | | //vector.add("-s_srs"); |
| | | //vector.add("EPSG:" + 4548); |
| | | //vector.add("-t_srs"); |
| | | //vector.add("EPSG:" + 4326); |
| | | vector.add("-ts"); |
| | | vector.add("" + width); |
| | | vector.add("" + height); |
| | | vector.add("-te"); |
| | | vector.add("" + layer.getExtension().getMinx()); |
| | | vector.add("" + layer.getExtension().getMiny()); |
| | | vector.add("" + layer.getExtension().getMaxx()); |
| | | vector.add("" + layer.getExtension().getMaxy()); |
| | | vector.add("-te_srs"); |
| | | vector.add("EPSG:" + 4326); |
| | | //vector.add("-te"); |
| | | //vector.add("" + layer.getExtension().getMinx()); |
| | | //vector.add("" + layer.getExtension().getMiny()); |
| | | //vector.add("" + layer.getExtension().getMaxx()); |
| | | //vector.add("" + layer.getExtension().getMaxy()); |
| | | //vector.add("-te_srs"); |
| | | //vector.add("EPSG:" + 4326); |
| | | vector.add("-r"); |
| | | vector.add("bilinear"); |
| | | vector.add("-of"); |
| | |
| | | Dataset ds = null; |
| | | try { |
| | | ds = gdal.Open(file, gdalconst.GA_ReadOnly); |
| | | if (null == ds || 0 == ds.getRasterCount() || null == ds.GetSpatialRef()) return; |
| | | if (null == ds || 0 == ds.getRasterCount()) return; |
| | | |
| | | double[] mm = new double[2]; |
| | | ds.GetRasterBand(1).ComputeRasterMinMax(mm, 0); |
| | |
| | | Dataset ds = null; |
| | | try { |
| | | ds = gdal.Open(files.get(i), gdalconst.GA_ReadOnly); |
| | | if (null == ds || 0 == ds.getRasterCount() || null == ds.GetSpatialRef()) return; |
| | | if (null == ds || 0 == ds.getRasterCount()) return; |
| | | |
| | | createWaterPng(dto, ds, layer, layer.getWaters().getData().get(i)); |
| | | } finally { |
| | |
| | | |
| | | double X = transform[0] + x * transform[1] + y * transform[2]; |
| | | double Y = transform[3] + x * transform[4] + y * transform[5]; |
| | | BuildingDto building = intersects(dto, X, Y); |
| | | if (null != building) continue; |
| | | //BuildingDto building = intersects(dto, X, Y); |
| | | //if (null != building) continue; |
| | | |
| | | int r = 0, g, b; |
| | | if (buffer[offset] - layer.getExtension().getMaxHeight() > 0) { |
| | |
| | | Dataset ds = null; |
| | | try { |
| | | ds = gdal.Open(files.get(i), gdalconst.GA_ReadOnly); |
| | | if (null == ds || 0 == ds.getRasterCount() || null == ds.GetSpatialRef()) return; |
| | | if (null == ds || 0 == ds.getRasterCount()) return; |
| | | |
| | | createFlowPng(dto, ds, layer, layer.getWaters().getData().get(i)); |
| | | } finally { |
| | |
| | | */ |
| | | private void copeLayerJson(ResultDto dto, LayerDto layer) throws IOException { |
| | | layer.getWaters().setFiles(null); |
| | | //String json = JSONUtil.toJsonPrettyStr(layer); |
| | | layer.getTerrain().setEpsg(null); |
| | | |
| | | String json = JSON.toJSONString(layer); |
| | | //String json = JSONUtil.toJsonPrettyStr(layer); |
| | | String filePath = dto.getOutPath() + File.separator + "layer.json"; |
| | | |
| | | FileWriter fw = new FileWriter(filePath); |
| | |
| | | bw.close(); |
| | | fw.close(); |
| | | } |
| | | |
| | | /** |
| | | * 处理降水曲线文件曲线图 |
| | | */ |
| | | public void copeRainFallJson(ResultDto dto, LayerDto layer) throws IOException, ParseException { |
| | | String rainGageFilePath = config.getInPath() + File.separator + dto.getServiceName() + File.separator + "RainGage.dat"; |
| | | String filePath = dto.getOutPath() + File.separator + "rainfall.json"; |
| | | |
| | | String line; |
| | | Map<String, Double> rainFallJsons = new LinkedHashMap<>(); |
| | | SimpleDateFormat sdf = new SimpleDateFormat("yyyy-MM-dd HH:mm"); |
| | | |
| | | BufferedReader br = new BufferedReader(new FileReader(rainGageFilePath)); |
| | | // 处理第一行数据 |
| | | if ((line = br.readLine()) != null) { |
| | | while ((line = br.readLine()) != null) { |
| | | // 处理每一行数据 |
| | | String[] rainFall = line.split(" "); |
| | | |
| | | if (rainFall.length < 7) continue; |
| | | |
| | | String sdt = rainFall[1] + "-" + rainFall[2] + "-" + rainFall[3] + " " + rainFall[4] + ":" + rainFall[5]; |
| | | BigDecimal num = new BigDecimal(rainFall[6]); |
| | | rainFallJsons.put("" + sdf.parse(sdt).getTime(), num.setScale(2, RoundingMode.HALF_UP).doubleValue()); |
| | | } |
| | | } |
| | | |
| | | FileWriter fw = new FileWriter(filePath); |
| | | BufferedWriter bw = new BufferedWriter(fw); |
| | | bw.write(JSON.toJSONString(rainFallJsons)); |
| | | bw.close(); |
| | | fw.close(); |
| | | } |
| | | } |