переход плоск-сферич Формула 25 действующего ГОСТа

[ivkomn]

неизвестная СК - это СК параметры которой не известны...

СК DST UTM(или ЛЮБАЯ ПРЯМОУГОЛЬНАЯ СК БЛИЗКАЯ К НИЖЕГОРОДСКОЙ ОБЛ)
СК SRC МЕЖДУНАРОДНАЯ СК с GPS координатами(что б эти координаты в яндекс подставить и он точку отобразил)

[ivkomn]

Можно чётко написать, какие конкретно параметры нужно передать?

//МЕЖДУНАРОДНАЯ СК с GPS координатами(что б эти координаты в яндекс подставить и он точку отобразил)
CoordinateReferenceSystem src = new CoordinateReferenceSystem(
"что это за строка",
"что это за массив строк",
Datum.WGS84,
что за Projection);

//UTM(или ЛЮБАЯ ПРЯМОУГОЛЬНАЯ СК БЛИЗКАЯ К НИЖЕГОРОДСКОЙ ОБЛ)
CoordinateReferenceSystem dst = new CoordinateReferenceSystem(
"что это за строка",
"что это за массив строк",
Datum.WGS84,
что за Projection);
CoordinateTransform ct = new CoordinateTransformFactory().createTransform(src,dst);

[freeExec]

На пальцах можно как-то объяснить:
1. выбираешь в своей системе координат нулевую точку и получаешь её GPS координаты в градусах
2. берешь класс такой-то библиотеки https://github.com/locationtech/proj4j
2. вставляешь в него GPS координаты в градусах(из п1), с такими-то параметрами
3. получаешь координаты в МЕЖДУНАРОДНОЙ прямоугольной системе координат
4. высчитываешь разницу по Х и У между своей СК и международной
5. добавляешь эту разницу ко всем своим координатам
6. переводишь свои координаты в GPS в градусы

Вы собрались условно "к милям прибавить километры" и радоваться получившемуся числу.

[trir]

CoordinateReferenceSystem dst = new CoordinateReferenceSystem("не важно", неизвестно, неизвестно, Projection.tmerc);
CoordinateReferenceSystem src = new CoordinateReferenceSystem("EPSG:4326", null, Datum.WGS84, Projection.longlat);

[ivkomn]

Projection.tmerc

Projection.longlat

что это такое? Как получить эти объекты?

[trir]

tmerc
longlat

[ivkomn]

Так?

Код: Выделить всё

	Projection tmerc = new TransverseMercatorProjection();
        Projection longlat = new LongLatProjection();
        longlat.setEllipsoid(Ellipsoid.KRASSOVSKY);
        longlat.setLonCDegrees(54);
        longlat.setProjectionLatitude(43);        
        CoordinateReferenceSystem src = new CoordinateReferenceSystem(
                "EPSG:4326",
                null,
                Datum.WGS84,
                longlat);
        CoordinateReferenceSystem dst = new CoordinateReferenceSystem(
                "UTM",
                null,
                Datum,
                tmerc);
        CoordinateTransform ct = new CoordinateTransformFactory().createTransform(src,dst);

Если хочу получить два числа double, отражающие плоские координаты в ск UTM, мне нужно что-то подставить в Datum и правильно инициализировать tmerc ? С какими параметрами?

[trir]

facepalm

Код: Выделить всё

        longlat.setEllipsoid(Ellipsoid.KRASSOVSKY);
        longlat.setLonCDegrees(54);
        longlat.setProjectionLatitude(43); 

зачем?
в UTM - Datum.WGS84

[ivkomn]

Код: Выделить всё

        Projection tmerc = new TransverseMercatorProjection();
        Projection longlat = new LongLatProjection();
        CoordinateReferenceSystem src = new CoordinateReferenceSystem(
                "EPSG:4326",
                null,
                Datum.WGS84,
                longlat);
        CoordinateReferenceSystem dst = new CoordinateReferenceSystem(
                "UTM",
                null,
                Datum.WGS84,
                tmerc);
        CoordinateTransform ct = new CoordinateTransformFactory().createTransform(src,dst);

??
Но когда же указывать gps координаты?

[trir]

метод transform

[ivkomn]

Код: Выделить всё

Projection tmerc = new TransverseMercatorProjection();
        Projection longlat = new LongLatProjection();
        CoordinateReferenceSystem src = new CoordinateReferenceSystem(
                "EPSG:4326",
                null,
                Datum.WGS84,
                longlat);
        CoordinateReferenceSystem dst = new CoordinateReferenceSystem(
                "UTM",
                null,
                Datum.WGS84,
                tmerc);
        CoordinateTransform ct = new CoordinateTransformFactory().createTransform(src,dst);
        ProjCoordinate wgs84Coord = new ProjCoordinate(54,43);
        ProjCoordinate utcCoord = new ProjCoordinate(0,0);
        ct.transform(wgs84Coord,utcCoord);
        System.out.println("X->"+utcCoord.x+":Y->"+utcCoord.y);

почему тогда:
X->NaN:Y->NaN

?

[trir]

попробуй сделать:
dst.setUTMZone(38);

[ivkomn]

Код: Выделить всё

TransverseMercatorProjection tmerc = new TransverseMercatorProjection();
        tmerc.setUTMZone(5);
        LongLatProjection longlat = new LongLatProjection();
        CoordinateReferenceSystem src = new CoordinateReferenceSystem(
                "EPSG:4326",
                null,
                Datum.WGS84,
                longlat);
        CoordinateReferenceSystem dst = new CoordinateReferenceSystem(
                "UTM",
                null,
                Datum.WGS84,
                tmerc);

        CoordinateTransform ct = new CoordinateTransformFactory().createTransform(src,dst);
        ProjCoordinate wgs84Coord = new ProjCoordinate(54,43);
        ProjCoordinate utcCoord = new ProjCoordinate(0,0);

        utcCoord = ct.transform(wgs84Coord,utcCoord);
        System.out.println("X->"+utcCoord.x+":Y->"+utcCoord.y);

тот же результат

[trir]

38

[ivkomn]

Код: Выделить всё

CRSFactory factory = new CRSFactory();

        CoordinateReferenceSystem srcCrs = factory.createFromName("EPSG:4326");
        CoordinateReferenceSystem dstCrs = factory.createFromName("EPSG:32638");

        ProjCoordinate srcCoord = new ProjCoordinate(54.9850521593413, 43.2725334459904);
        ProjCoordinate dstCoord = new ProjCoordinate();

        CoordinateTransform ct_my = new CoordinateTransformFactory().createTransform(srcCrs,dstCrs);
        ct_my.transform(srcCoord,dstCoord);

        System.out.println("X:"+dstCoord.x+":Y:"+dstCoord.y);

не мытьём, так катаньем ))
Получил результат

Код: Выделить всё

X:1310495.622518064:Y:4839745.462222159

Вопрос: это координаты на плоскости, выраженные в метрах относительно какой-то точки 0,0?