[add]tif裁剪

张洋洋

2025-02-12 75a3851294865b31d228872ea677f80cd0cda028

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
package com.se.simu.utils;
 
import com.alibaba.fastjson.JSONArray;
import org.osgeo.proj4j.BasicCoordinateTransform;
import org.osgeo.proj4j.CRSFactory;
import org.osgeo.proj4j.CoordinateReferenceSystem;
import org.osgeo.proj4j.ProjCoordinate;
 
public class ProjectionToGeographicUtil {
    public static JSONArray getPoint(Double x,Double y) {
        // 创建 CRSFactory 对象
        CRSFactory crsFactory = new CRSFactory();
        // 定义投影坐标系统，这里以 UTM 投影为例，zone 33N
        CoordinateReferenceSystem sourceCRS = crsFactory.createFromName("EPSG:4548");
        // 定义地理坐标系统，这里使用 WGS84
        CoordinateReferenceSystem targetCRS = crsFactory.createFromName("EPSG:4326");
        // 创建坐标转换对象
        BasicCoordinateTransform transform = new BasicCoordinateTransform(sourceCRS, targetCRS);
        //lon=116.683795  lat=39.90849042
        // 输入投影坐标，以 UTM 投影坐标为例，单位是米
        ProjCoordinate sourceCoord = new ProjCoordinate(x, y);
        ProjCoordinate targetCoord = new ProjCoordinate();
 
        // 执行坐标转换
        transform.transform(sourceCoord, targetCoord);
        JSONArray vertice = new JSONArray();
        vertice.add(targetCoord.x);
        vertice.add(targetCoord.y);
        vertice.add(0.0);
        return vertice;
    }
 
    public static JSONArray get4548Point(Double x,Double y) {
        // 创建 CRSFactory 对象
        CRSFactory crsFactory = new CRSFactory();
        // 定义投影坐标系统，这里以 UTM 投影为例，zone 33N
        CoordinateReferenceSystem sourceCRS = crsFactory.createFromName("EPSG:4326");
        // 定义地理坐标系统，这里使用 WGS84
        CoordinateReferenceSystem targetCRS = crsFactory.createFromName("EPSG:4548");
        // 创建坐标转换对象
        BasicCoordinateTransform transform = new BasicCoordinateTransform(sourceCRS, targetCRS);
        // 输入投影坐标，以 UTM 投影坐标为例，单位是米
        ProjCoordinate sourceCoord = new ProjCoordinate(x, y);
        ProjCoordinate targetCoord = new ProjCoordinate();
        // 执行坐标转换
        transform.transform(sourceCoord, targetCoord);
        JSONArray vertice = new JSONArray();
        vertice.add(targetCoord.x);
        vertice.add(targetCoord.y);
        return vertice;
    }
 
    public static void main(String[] args) {
        System.out.println(get4548Point(116.66,39.8962));
    }
}