java 读取las文件

要读取LAS文件，需要使用LAS Java工具包。这个工具包提供了一个LASReader类，可以读取LAS文件中的数据。以下是一个简单的示例代码：

import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;
import org.geolatte.geom.G2D;
import org.geolatte.geom.Point;
import org.geolatte.geom.builder.DSL;
import org.geolatte.geom.builder.PointBuilder;

import java.io.File;
import java.io.IOException;
import java.util.List;

import org.geolatte.geom.builder.DSL;
import org.geolatte.geom.builder.GeometryBuilder;
import org.geolatte.geom.crs.CoordinateReferenceSystems;
import org.geolatte.geom.crs.CrsId;
import org.geolatte.geom.jts.JTS;

import com.google.common.collect.Lists;
import com.google.common.primitives.Floats;

import io.vertx.core.json.JsonArray;
import io.vertx.core.json.JsonObject;
import las.LASHeader;
import las.LASReader;

public class ReadLasFile {

  public static void main(String[] args) throws IOException {

    // 读取LAS文件
    LASReader reader = new LASReader(new File("example.las"));
    LASHeader header = reader.getHeader();

    // 获取点云数据
    List<double[]> pointCloud = reader.readPointCloud();

    // 将点云数据转换为Point对象
    List<Point<G2D>> points = Lists.newArrayList();
    for (double[] point : pointCloud) {
      PointBuilder<G2D> builder = DSL.point();
      builder.coordinate(DSL.g(2).p(point[0], point[1]));
      Point<G2D> p = builder.build();
      points.add(p);
    }

    // 创建多边形对象
    GeometryBuilder<G2D> builder = DSL.geometry(CoordinateReferenceSystems.PROJECTED_2D_METER);
    builder.polygon(DSL.linearRing(points.toArray(new Point[points.size()])));
    Geometry<G2D> polygon = builder.build();

    // 将多边形对象转换为JSON字符串
    JsonObject json = new JsonObject();
    json.put("type", "Feature");
    json.put("geometry", JTS.toGeoJson(polygon));
    System.out.println(json.toString());
  }
}

这个示例代码读取了一个名为"example.las"的LAS文件，并将其转换为多边形对象。该代码使用了LAS Java工具包、Apache Commons Math、Geolatte Geometry和Guava库。