@@ -34,8 +34,8 @@ namespace {
3434 // const auto epsg_offset = geo_ref1.getOffset(options.epsg_code);
3535 // geo_ref1.setReferencePoint(epsg_offset);
3636
37- plateau::geometry::GeoReference geo_ref2 (options.coordinate_zone_id , options.reference_point , options.unit_scale , options.mesh_axes );
38- const auto epsg_offset = geo_ref2 .getOffset (options.epsg_code );
37+ plateau::geometry::GeoReference geo_ref (options.coordinate_zone_id , options.reference_point , options.unit_scale , options.mesh_axes );
38+ const auto epsg_offset = geo_ref .getOffset (options.epsg_code );
3939
4040 // const auto diff = options.reference_point - epsg_offset;
4141 // plateau::geometry::GeoReference geo_ref3(options.coordinate_zone_id, options.reference_point - diff , options.unit_scale, options.mesh_axes);
@@ -46,40 +46,45 @@ namespace {
4646 plateau::geometry::GeoReference geo_ref3 (options.coordinate_zone_id , TVec3d (), options.unit_scale , options.mesh_axes );
4747 const auto original_ref = geo_ref3.unproject (options.reference_point );
4848
49- plateau::geometry::GeoReference geo1 (options.coordinate_zone_id );
50- plateau::geometry::GeoReference geo2 (8 ); // 8固定
49+
50+ // int zone_id = options.coordinate_zone_id;
51+ int zone_id = 8 ; // 8固定
52+ plateau::geometry::GeoReference geo (zone_id);
5153
5254 GeoCoordinate ref_point_epsg = ReferencePointFactory::GetReferencePoint (options.epsg_code ); // EPSGの基準点
53- const auto prj_epsg = geo1 .project (ref_point_epsg);
54- plateau::geometry::GeoReference geo_ref_epsg (options. coordinate_zone_id , prj_epsg, options.unit_scale );
55+ const auto prj_epsg = geo .project (ref_point_epsg);
56+ plateau::geometry::GeoReference geo_ref_epsg (zone_id , prj_epsg, options.unit_scale );
5557
5658 GeoCoordinate ref_point_zukaku (37.4258 , 138.7378 , 0 ); // 08EE751の中心
57- const auto prj2 = geo1 .project (ref_point_zukaku);
58- plateau::geometry::GeoReference geo_ref_zukaku (options. coordinate_zone_id , prj2, options.unit_scale );
59+ const auto prj2 = geo .project (ref_point_zukaku);
60+ plateau::geometry::GeoReference geo_ref_zukaku (zone_id , prj2, options.unit_scale );
5961
60- GeoCoordinate ref_point_epsg_static (36 , 138.5 , 0 ); // EPSG:10169 の基準点
61- const auto prj3 = geo1 .project (ref_point_epsg_static);
62- plateau::geometry::GeoReference geo_ref_epsg_static (options. coordinate_zone_id , prj3, options.unit_scale );
62+ // GeoCoordinate ref_point_epsg_static(36, 138.5, 0); //EPSG:10169 の基準点
63+ // const auto prj3 = geo .project(ref_point_epsg_static);
64+ // plateau::geometry::GeoReference geo_ref_epsg_static(zone_id , prj3, options.unit_scale);
6365
6466 try {
6567 const auto pos = PolygonMeshUtils::cityObjPos (city_obj);
6668
67- const auto unprojected = geo_ref2.unproject (pos + epsg_offset);
68- // const auto unprojected2 = geo_ref2.unproject(epsg_offset - pos);
69- // const auto unprojected2 = geo_ref2.unproject(pos);
70- // const auto unprojected3 = geo_ref4.unproject(pos + epsg_offset);
71- // const auto unprojected5 = geo_ref5.unproject(epsg_offset- pos);
69+ // const auto unprojected = geo_ref.unproject(pos + epsg_offset);
70+
71+ // const auto unprojected_epsg = geo_ref_epsg.unproject(pos);
72+ // const auto unprojected_zukaku = geo_ref_zukaku.unproject(pos);
73+
74+ // const auto unprojected_epsg_offset = geo_ref_epsg.unproject(pos + epsg_offset);
75+ // const auto unprojected_zukaku_offset = geo_ref_zukaku.unproject(pos + epsg_offset);
76+
77+ // converted
78+ // const auto pos_converted = GeoReference::convertAxisToENU(geometry::CoordinateSystem::WUN, pos);
79+ const TVec3d pos_converted = { pos.y , pos.x , pos.z }; // x,y反転
7280
73- const auto unprojected_epsg = geo_ref_epsg.unproject (pos);
74- const auto unprojected_zukaku = geo_ref_zukaku.unproject (pos);
75- const auto unprojected_epsg_static = geo_ref_epsg_static.unproject (pos);
81+ const auto unprojected1 = geo_ref.unproject (pos_converted);
7682
77- const auto unprojected_epsg_offset = geo_ref_epsg.unproject (pos + epsg_offset);
78- const auto unprojected_zukaku_offset = geo_ref_zukaku.unproject (pos + epsg_offset);
79- const auto unprojected_epsg_static_offset = geo_ref_epsg_static.unproject (pos + epsg_offset);
83+ const auto unprojected_epsg2 = geo_ref_epsg.unproject (pos_converted);
84+ const auto unprojected_zukaku2 = geo_ref_zukaku.unproject (pos_converted);
8085
8186 for (const auto & extent : extents) {
82- if (extent.contains (unprojected ))
87+ if (extent.contains (unprojected_epsg2 ))
8388 return false ;
8489 }
8590
0 commit comments