geofiles 1.2.0a2

Geofiles is a project for importing and exporting geo-referenced 3D geometry files.

Geo-referenced Geometry File Formats

Classic geometry file formats as .obj, .off, .ply, .stl or .dae do not support the utilization of coordinate systems besides from a local system, that can not be defined more precisely. This feature is a major requirement for global applications, exchanging geo-referenced models e.g. in the context of outdoor augmented reality applications.

For this reason, the present project evaluates different possibilities of geo-referenced geometry files.

Next to wide-spread standards as GeoVRML, X3D, CityJSON, GeoJSON, GML or KML, we introduce four geo-referenced extensions called .geoobj, .geooff, .geoply and .geostl. While, the named standard formats allow defining objects with multiple additional features, they also come with an overhead according to the file size and an increased structural complexity with disadvantages to the read performance. This overhead has to be minimized in many use cases and for this reason, some file formats are more suitable than others. For completeness: the overhead can be further reduced using binary representations (like binary obj or stl) instead of textual ones, but since not all named formats support a binary mode, this is ignored in favor of human-readability within this project and comparison.

Geo-referenced Extensions

In this chapter you can find multiple extensions of classic geometry file formats. All of these extensions support two major features:

1. Defining the coordinate reference system (CRS) of the used vertices' coordinates
2. Defining an optional origin, which represents an absolute geo-referenced location and an implicit local Cartesian coordinate frame for the models with a metric unit.

Additional features

Some proposed file formats support additional features such as:

• Geographical extent
• Transformation information (local per object or global for all vertices)
  • Scaling factor
  • Translation in meters
  • Rotation in degrees

GeoOBJ

The .geoobj file format extends the classic .obj file format with two line-types:

1. The crs line type is used to define the coordinate system which is used within the file
2. The or line type is used to define the optional origin of the vertices

Example:

crs urn:ogc:def:crs:EPSG::4326
or 48.3028533074941 14.2842865755919 279.307006835938

Addons

Next to the geo-referencing features, the .geoobj extension also supports exchanging scale, rotation, as well as translation information using the following line prefixes:

1. sc for adding scale information (s is already defined in .obj for smoothing groups)
2. t for translation information
3. r for rotation angular information

Example:

sc 1.5 2 5
t 10 -5 4
r 90 45 10

Note that: if either a scaling, a translation or a rotation information is stated after an object definition with the prefix o, the transformation is not interpreted globally, but locally for the specific object.

Next to that it also supports the optional geographical extent meta information, containing the minimal (first three values) and maximal (remaining three values) coordinate value expressions using the line prefix e. This information can be useful for filtering geo-referenced files without any need to iterate all vertices.

e -0.5 -0.5 -0.5 0.5 0.5 0.5

GeoOFF

The .geooff file format extends the classic .off file. For this we introduce a new file header using GeoOFF instead of OFF. The next non-empty line after this header is used to define the crs.

Example:

GeoOFF
urn:ogc:def:crs:EPSG::4326

GeoOFF supports alternative headers, using different postfix values. The pattern of the header is based on the OFF header definition and is defined like [ST][C][N][4][n]GeoOFF[o][e][s][t][r]. Like this GeoOFF is able to support:

• An absolute origin using the o header postfix symbol
• Extent information using the e header postfix symbol
• Global scaling information using the s header postfix symbol
• Global translation information using the t header postfix symbol
• Global rotation information using the r header postfix symbol

Example:

GeoOFFostr
urn:ogc:def:crs:OGC:2:84
14.2842798233032 48.30284881591775 279.807006835938
2 2 2
10 50 100
90 0 0

Note: The .off prefixes are currently not supported in the implementations.

GeoPLY

The .geoply file format extends the classic .ply file format with two header-line-types:

1. The crs line type is used to define the coordinate system which is used within the file
2. The origin line type is used to define the optional origin of the vertices

In addition to those types it also changes the header from ply to geoply.

geoply
crs urn:ogc:def:crs:EPSG::4326
o 48.3028533074941 14.2842865755919 279.307006835938
...
end_header

Addons

Next to the geo-referencing features, the .geply extension also supports exchanging global scale, rotation, as well as translation information using the following line prefixes:

1. scale for adding scale information
2. translate for translation information
3. rotate for rotation angular information

Example:

geoply
...
scale 1.5 2 5
translate 10 -5 4
rotate 90 45 10
...
end_header

Next to that it also supports the optional geographical extent meta information, containing the minimal (first three values) and maximal (remaining three values) coordinate value expressions using the extent header. This information can be useful for filtering geo-referenced files without any need to iterate all vertices.

geoply
...
extent -0.5 -0.5 -0.5 0.5 0.5 0.5
...
end_header

GeoSTL

The .geostl file format extends the classic .stl file using the geosolid root element. Followed by the geosolid prefix of the file format a meta-data tuple is introduced. This tuple consists of the crs at the first position, followed by the optional origin coordinates and finally the optional stl name.

geosolid urn:ogc:def:crs:OGC:2:84 48.3028533074941 14.2842865755919 279.307006835938 fileName

File format comparison

The different named file formats come with a variable amount of features according to e.g. the supported CRS, an origin for specifying a local coordinate system, the meta information of the minimal and maximal coordinates (geographical extent), multi-object support or also differ in the representation of vertices. To take up the last point, some formats use e.g. central vertex lists with referencing indices in the face definition and others re-define the vertices within every indiviudal face. Next to that there are many different other features (e.g. smoothing groups in .obj, geographicalExtent in CityJSON or exact property definitions in .ply), which vary between the file formats and lead to a diverse semantic expressiveness.

Format	Base format	Encoding	Coordinate Reference System	Multiple Objects	Vertex References	Origin Support	Transformation Information	Geographical Extent	Semantic Expressiveness
CityJSON	JSON	Text	Any	Yes	Yes	No	No	Yes	++
GeoJSON	JSON	Text	WGS 84	Yes	No	No	No	No	+
GeoObj	OBJ	Text *	Any	Yes	Yes	Yes	Yes	Yes	~
GeoOFF	OFF	Text *	Any	No	Yes	Yes	Yes	Yes	-
GeoPLY	PLY	Text *	Any	No	Yes	Yes	Yes	Yes	~
GeoSTL	STL	Text *	Any	No	No	Yes	No	No	--
GML	XML	Text	Any	Yes	No	No	No	No	++
KML	XML	Text or compressed (KMZ)	WGS 84	Yes	No	No	No	No	++
GeoVRML	VRML	Text or Binary	WGS84 / EPSG4326 / UTM	Yes	Yes	Yes	Yes	Yes	+
X3D	XML, VRML	Text or Binary	Multiple supported CRS	Yes	Yes	Yes	Yes	Yes	++

* Currently, only text support, but the base format would support binary, so a binary extension would be possible.

Size comparison

The following file comparison uses the minimal required sub-set of the specific file formats to represent geo-referenced 3D models. Additional features (e.g. GeographicalExtent in CityJSON) are not considered as far as possible, to avoid a negative bias of the comparison. Next to that also the used reader does not support all features of the input files (e.g. .mtb material information) and may use different line-ending symbols (Windows vs Linux). This would lead to the situation that in some cases the input file size is greater than the geo-referenced version. For this reason we have decided to do a normalization first by reading the input files with our custom reader and exporting the normalized representation using our custom obj writer. The exported and normalized files are used as basis of comparison.

You can find the used obj files here. Also note that, X3D is currently not explicitly supported by the present framework since it can be encoded using XML or VRML. Next to that it is the successor of VRML and supports many features of GeoVRML within the geospatial extension with only minor differences. For this reason it would result in a comparable size when using the VRML encoding and is left out in the following comparison.

file	original	cityjson	geojson	geoobj	geoobj-origin	geooff	geooff-origin	geoply	geoply-origin	geostl	geostl-origin	geovrml	geovrml-origin	gml	kml
alligator.obj	184.07	308.53	1402.25	278.69	184.16	272.42	177.89	272.57	178.04	1396.25	877.53	1298.15	779.87	2202.64	2237.52
armadillo.obj	3608.54	5279.43	23622.39	4791.05	3608.62	4693.41	3510.98	4693.56	3511.14	23476.92	16381.74	22142.54	15047.79	36998.27	37583.91
beast.obj	3332.12	2850.39	9550.49	4087.92	3332.20	2566.04	1810.32	2566.19	1810.48	9686.46	6660.98	9516.06	6491.00	13723.15	13944.12
beetle-alt.obj	1364.18	2037.85	9120.80	1847.80	1364.26	1808.95	1325.41	1809.10	1325.57	9067.75	6249.00	8517.67	5699.35	14292.74	14519.08
beetle.obj	125.25	105.10	485.01	155.55	125.33	92.38	62.16	92.53	62.32	599.52	437.23	444.46	282.61	759.87	771.74
cheburashka.obj	420.93	667.26	3148.23	601.94	421.01	588.91	407.98	589.06	408.14	3129.51	2044.29	2914.98	1830.20	4932.36	5010.33
cow.obj	182.93	286.76	1369.35	258.21	183.01	252.53	177.33	252.68	177.49	1361.40	909.95	1264.59	813.57	2146.05	2179.90
fandisk.obj	391.15	646.64	3047.66	583.21	391.23	570.56	378.58	570.71	378.73	3031.65	1879.89	2823.32	1671.99	4779.88	4855.58
happy.obj	3443.42	5200.30	23271.05	4718.75	3443.50	4622.54	3347.30	4622.69	3347.45	23134.59	15475.17	21818.18	14159.19	36462.98	37040.56
homer.obj	378.70	600.33	2832.50	541.53	378.78	529.80	367.05	529.95	367.21	2815.80	1839.43	2622.71	1646.77	4438.16	4508.31
horse.obj	4847.43	5161.59	23251.47	4687.92	4847.51	4593.21	4752.80	4593.36	4752.96	23025.46	23982.47	21729.74	22687.18	36224.65	36792.65
igea.obj	9716.50	14490.75	63413.22	13178.60	9716.58	12916.19	9454.18	12916.34	9454.33	63040.76	42268.45	59765.41	38993.53	99360.66	100934.83
lucy.obj	3728.41	5310.72	23867.91	4822.38	3728.49	4724.74	3630.85	4724.89	3631.01	23660.91	17097.44	22327.21	15764.17	37242.99	37828.60
max-planck.obj	3628.48	5296.99	23729.97	4808.40	3628.57	4710.58	3530.75	4710.73	3530.91	23558.65	16499.04	22223.06	15163.88	37107.86	37693.59
nefertiti.obj	3621.60	5288.27	23695.48	4800.08	3621.68	4702.47	3524.07	4702.62	3524.23	23529.73	16459.26	22196.35	15126.31	37066.29	37651.70
ogre.obj	3820.09	6569.39	29260.92	5963.30	3820.17	5841.82	3698.69	5841.97	3698.85	29090.92	16271.91	27450.40	14631.82	45852.03	46578.48
rocker-arm.obj	666.93	1008.57	4742.64	910.28	667.02	890.65	647.39	890.80	647.55	4714.26	3254.51	4394.58	2935.25	7430.39	7547.93
spot.obj	334.84	289.43	1380.09	409.11	334.92	254.89	180.70	255.04	180.86	1372.51	927.39	1275.06	830.38	2163.75	2197.91
stanford-bunny.obj	2442.78	3719.00	16390.43	3377.29	2442.86	3307.08	2372.64	3307.23	2372.80	16294.39	10877.71	15350.55	9934.31	25682.41	26089.19
suzanne.obj	48.62	41.65	149.19	62.89	48.70	37.04	22.85	37.19	23.01	175.03	119.80	144.37	89.58	213.99	217.22
teapot.obj	202.44	341.20	1489.62	309.19	202.52	302.07	195.40	302.22	195.55	1481.45	937.85	1377.36	834.20	2335.36	2372.23
woody.obj	36.33	63.44	297.24	56.92	36.42	55.56	35.06	55.71	35.22	295.92	186.12	272.52	163.16	466.95	474.21
xyzrgb_dragon.obj	9283.49	13455.47	59107.69	12234.45	9283.57	11990.17	9039.29	11990.32	9039.44	58735.04	41041.19	55643.17	37949.76	92555.01	94019.69

Note: Sizes are in KiloByte (KB)

An interactive visualization of the size comparison can be found here.

Getting started

To set up and use the project have a look at the detailed description here

You can install the framework via pip:

pip install geofiles

Importing files

The present project supports multiple reader implementations for importing (geo-referenced) geometry files (.obj, .geoobj, .ply, .geoply, .off, .geooff, .stl, .geostl). Using one of these readers is the entrypoint to the framework and allows to create an in-memory geometry model using the GeoObjectFile class. Note that only a subset of the features of the named files are currently supported. So reading files with non-supported features may result in a loss of information (e.g. smoothing groups in .obj, or exact property definitions of .ply)

reader = GeoObjReader()
path = "mygreatfile.geoobj"
with open(path) as file:
    geoObjFile: GeoObjectFile = reader.read(file)

Converting

The present framework supports different conversion methodologies as converting from one to another coordinate reference system. Next to the CRS-conversion the framework also supports to transform between origin and non-origin based representations, as well as between geo-referenced and local representations. Finally, there is also basic support for transforming (scale, rotation, translation) a model's vertices.

# 1. Change coordinate system
converter = CrsConverter()
converted = converter.convert(geoObjFile, "urn:ogc:def:crs:OGC:2:84")

# 2. Change to origin-based representation
originconverter = OriginConverter()
origin_based = originconverter.to_origin(converted)

# 3. Apply transformation
transformer = Transformer()
origin_based.scaling = [2, 2, 2]
origin_based.translation = [5, 10, -5]
origin_based.rotation = [45, 30, 90]
transformed = transformer.transform(origin_based)

Determining Geographical Extent

The geographical extent of a file can be determined in two ways:

1. Only searching for the min and max coordinates of the values (not considering transformation or origin information)
2. Determining the geographical extent considering all available meta information

Depending on the use case one variant is more suitable than the other. If you are going to ignore transformation information in your application, the additional overhead of the second method is not required, otherwise if you want to know the extents considering this meta information you have to use the ExtentCalculator class.

# 1. Classic geographical extent 
geoObjFile: GeoObjectFile = ...
geoObjFile.update_extent()

# 2. Using the ExtentCalculator
objFileWithExtents = ExtentCalculator.update_extent(geoObjFile, True, True)

Exporting files

Finally, the in-memory model representations can be re-written to your hard drive using one of the writer implementations. Note:

1. That you can find more writers than readers in this project for evaluation purposes. The additional writers only use a minimal subset of the specific file format's features as base for comparison and are not complete.
2. That some file formats presuppose a specific CRS (e.g. .kml requires vertices in Wgs84 representation)
3. Most file formats do not support transformation (scale, rotation, translation) meta-information. A model's vertices have to be transformed first before exported to such a file format.
4. The writers will automatically append the specific file type (unless you set append_file_type to False)

writer = GeoPlyWriter()
writer.write("mygreatfile.geoply", transformed, append_file_type=False)

FAQ

• Why yet another 3D geometry file format like .geoobj?
  • During our research in the context of outdoor augmented reality applications, we were looking for a possibility for exchanging geo-referenced geometry models. In this context, the other named file formats come with a too high overhead (e.g. XML tags or not required meta information as object types like in CityJson) in our opinion and are for this reason not ideal.
• You describe multiple geo-referenced file formats. Which one should I use for geo-referenced 3D models?
  • This depends on the use case. If you have to exchange the models with as little overhead as possible we recommend using the proposed .geoply, .geoobj or .geooff format extensions. If you require semantic expressiveness, you should prefer other formats like CityJson or GML.
• How are vertices defined, if I use the origin-based approach of .geoobj, .geoply, .geooff or .geostl?
  • In the origin-based version, vertices are represented within a local Cartesian coordinate system with the defined origin as coordinate system origin (0, 0, 0).
  • The local coordinate system is intended as a left handed system.
  • The local coordinate system uses the x-axis as abscissa axis (width information), y-axis as ordinate axis (length information) and z-axis as applicate axis (height information).
  • The units used in this type of coordinate system are assumed to be in meters.
• How is the transformation information defined?
  • The proposed transformation information is separated into tuples (one value per axis) for translation, rotation and scale. For the translation, meter based offsets are intended to be used, the rotation is based on degrees and the scale tuple is represented using numeric factors.

Contributing

First make sure to read our general contribution guidelines.

Licence

Copyright (c) 2021 the original author or authors. DO NOT ALTER OR REMOVE COPYRIGHT NOTICES.

This Source Code Form is subject to the terms of the Mozilla Public License, v. 2.0. If a copy of the MPL was not distributed with this file, You can obtain one at https://mozilla.org/MPL/2.0/.

Research

If you are going to use this project as part of a research paper, we would ask you to reference this project by citing it.