搭建GDAL JAVA环境；DXF转KML；坐标转换；PROJ: proj_create_from_database - ToB企服应用市场:ToB评测及商务社交产业平台

import org.gdal.gdal.gdal;
import org.gdal.ogr.ogr;
//测试案例
class GdalDemo {
//获取所有驱动的名称并打印
public void printDriver() {
gdal.AllRegister();
int count = ogr.GetDriverCount();
for (int i = 0; i < count; i++) {
String driverName = ogr.GetDriver(i).getName();
System.out.print(driverName + "\t");
}
gdal.GDALDestroyDriverManager();
}
public static void main(String[] args) {
GdalDemo demo = new GdalDemo();
demo.printDriver();
}
}

复制代码

/**
* 使用gdal 解析dxf，并将数据读取转为KML 或 GEOJSON
* 包含坐标转换
*/
public class DXFToGeoJSONConverter {
public static void main(String[] args) {
Boolean converter = converter("E://test//1.kml", "E://test//白马2000-3-108.dxf");
if (converter) {
System.out.println("成功");
} else {
System.out.println("失败");
}
}
/**
*
* @param outPath 输出文件全路径
* @param inPath 读取文件全路径
* @return
*/
public static Boolean converter(String outPath, String inPath) {
// 注册所有的驱动
//gdal.AllRegister();
ogr.RegisterAll();
gdal.SetConfigOption("DWG_ENCODING", "UTF-8");//DWG文件编码
gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");//支持中文路径
gdal.SetConfigOption("SHAPE_ENCODING", "CP936");//属性表字段支持中文
DataSource dataSource = ogr.Open(inPath, 0);
if (dataSource == null) {
System.out.println("打开文件失败！" + gdal.GetLastErrorMsg());
return false;
}
System.out.println("打开文件成功！");
//构造坐标系转换对象
SpatialReference src = new SpatialReference();
// src.ImportFromProj4("+proj=tmerc +lat_0=0 +lon_0=108 +k=1 +x_0=36500000 +y_0=0 +ellps=GRS80 +units=m +no_defs");
src.ImportFromEPSG(4545); // CGCS2000 中央经线 108 3分带 36带
System.out.println("原坐标系：" + src.ExportToProj4());
SpatialReference dest = new SpatialReference();
dest.ImportFromEPSG(4326); // WGS84
System.out.println("目标坐标系：" + dest.ExportToProj4());
// SpatialReference dest = dest1.CloneGeogCS();
CoordinateTransformation transformation = new CoordinateTransformation(src, dest);//源,目标
//1.创建驱动程序
Driver dv = ogr.GetDriverByName("KML"); //直接转KML
// Driver dv = ogr.GetDriverByName("GeoJSON"); //直接转GEOJSON
int layerCount = dataSource.GetLayerCount();
System.out.println("总图层数：" + layerCount);
DataSource ds = dv.CreateDataSource(outPath);
for (int i = 0; i < layerCount; i++) {
// 2.循环读取每个图层
Layer inLayer = dataSource.GetLayerByIndex(i);
if (inLayer == null) {
continue;
}
String layerName = inLayer.GetName();
//3.每个DataSource创建一个输出的layer
Layer outLayer = ds.CreateLayer(layerName, dest, ogrConstants.wkbUnknown, null);
long featureCount = inLayer.GetFeatureCount();
System.out.println("总图形数：" + featureCount);
for (long j = 0; j < featureCount; j++) {
Feature feature = inLayer.GetFeature(j);
feature.GetStyleString();
//4.获取输入的几何对象，调用其转换方法，转换后存入新的layer中
Geometry geometry = feature.GetGeometryRef();
geometry.Transform(transformation);
//注意，转换为后时经度纬度的顺序，如果需要纬度在前，需要交换xy
// geometry.SwapXY();
geometry.CloseRings();
outLayer.CreateFeature(feature);
}
}
System.out.println("转换成功!!!");
return true;
}
}

复制代码

/**
* 用gdal进行坐标转换
*/
public class CoordinateTransformationExample {
public static void main(String[] args) {
// 注册GDAL驱动
gdal.AllRegister();
// 输入的坐标
double inputX = 36587998.113885;
double inputY = 3796047.384520;
// 创建源坐标系统 (CGCS2000 3度带，中央经线108)
SpatialReference srcSRS = new SpatialReference();
srcSRS.ImportFromProj4("+proj=tmerc +lat_0=0 +lon_0=108 +k=1 +x_0=36500000 +y_0=0 +ellps=GRS80 +units=m +no_defs");
// 创建目标坐标系统 (WGS84)
SpatialReference tgtSRS = new SpatialReference();
tgtSRS.ImportFromEPSG(4326);
// 创建坐标转换对象
CoordinateTransformation transformation = new CoordinateTransformation(srcSRS, tgtSRS);
// 执行坐标转换
double[] transformedCoords = transformation.TransformPoint(inputX, inputY);
// 输出转换后的经纬度（转换后经纬度位置交换了，注意不要取错参数）
System.out.println("Transformed Coordinates (WGS84):");
System.out.println("Longitude: " + transformedCoords[1]);
System.out.println("Latitude: " + transformedCoords[0]);
}
}

复制代码

ERROR 1: PROJ: proj_create_from_database: D:\ideawork\pkgs\proj-6.2.1-h9f7ef89_0\Library\share\proj\proj.db lacks DATABASE.LAYOUT.VERSION.MAJOR / DATABASE.LAYOUT.VERSION.MINOR metadata. It comes from another PROJ installation.

复制代码