0

ここに画像の説明を入力してください lzw 圧縮されたラスター (tiff) をテキスト形式に変換しようとしています (つまり、各ピクセル重心とバンド 1 値を抽出します..変換中にラスターの一部が失われていることがわかります..しかし、ミスにパターンがあります。データ (データは不均等な縞模様で入ってきます) スナップショットをご覧ください..黄色の点は変換された重心を示しますが、黄色の点は変換中に失われたデータを示しています..

タイル サイズ 128 と 256 を試しましたが、どちらの場合も出力が異なりました。

List<String> tiffExtensions = new ArrayList<>();
tiffExtensions.add(".tif");
tiffExtensions.add(".TIF");
tiffExtensions.add(".tiff");
tiffExtensions.add(".TIFF");

final scala.Option<CRS> crsNone = scala.Option.apply(null);
final scala.Option<Object> objectNone = scala.Option.apply(null);
final scala.Option<Object> numPartitionsObject = scala.Option.apply(new Integer(10));
final scala.Option<Object> tileSize = scala.Option.apply(Integer.parseInt("256"));
final scala.Option<Object> partitionBytes = scala.Option.apply(128l * 1024 * 1024);
HadoopGeoTiffRDD.Options options = new HadoopGeoTiffRDD.Options$().apply(
        JavaConverters.asScalaIteratorConverter(tiffExtensions.iterator()).asScala().toSeq(), crsNone,
        HadoopGeoTiffRDD.GEOTIFF_TIME_TAG_DEFAULT(), HadoopGeoTiffRDD.GEOTIFF_TIME_FORMAT_DEFAULT(), tileSize,
        numPartitionsObject, partitionBytes, objectNone);
RDD<Tuple2<ProjectedExtent, Tile>> rasterImageRdd = HadoopGeoTiffRDD.spatial(new Path(rastedImageDir), options,
        javaSparkContext.sc());
JavaRDD<Tuple2<ProjectedExtent, Tile>> rasterImageJavaRdd = rasterImageRdd.toJavaRDD();
JavaRDD<String> pixelRdd = rasterImageJavaRdd
        .flatMap(new FlatMapFunction<Tuple2<ProjectedExtent, Tile>, String>() {
            private static final long serialVersionUID = -6395159549445351446L;

            public Iterator<String> call(Tuple2<ProjectedExtent, Tile> v1) throws Exception {
                ArrayList<String> list = new ArrayList<String>();
                Tile t = v1._2;
                ProjectedExtent projectedExtent = v1._1;
                ProjectedRaster<CellGrid> r = new ProjectedRaster<CellGrid>(
                        new Raster<CellGrid>(t, projectedExtent.extent()), projectedExtent.crs());
                GeometryFactory geometryFactory = new GeometryFactory(new PrecisionModel(), 4283);
                WKTWriter wktWriter = new WKTWriter();
                for (int i = 0; i < t.rows(); i++) {
                    for (int j = 0; j < t.cols(); j++) {
                        StringBuilder sb = new StringBuilder();
                        if (!Double.isNaN(t.getDouble(j, i))) {
                            Double elevation = t.getDouble(j, i);
                            Tuple2<Object, Object> longLatTupel = r.raster().rasterExtent().gridToMap(j, i);
                            if (longLatTupel._2() != null && longLatTupel._1() != null) {
                                Double latitude = Double.parseDouble(longLatTupel._2() + "");
                                Double longitude = Double.parseDouble(longLatTupel._1() + "");
                                Point point = geometryFactory.createPoint(new Coordinate(longitude, latitude));
                                sb.append(elevation).append(TAB_SEP);
                                sb.append(point);
                                list.add(sb.toString());
                            }
                        }
                    }
                }
                return list.iterator();
            }
        });

_

4

0 に答える 0