def test_osr_GetCRSInfoListFromDatabase():
l = osr.GetCRSInfoListFromDatabase('EPSG')
found = False
for record in l:
if record.auth_name == 'EPSG' and record.code == '2154':
assert record.name == 'RGF93 / Lambert-93'
assert record.type == osr.OSR_CRS_TYPE_PROJECTED
assert not record.deprecated
assert record.bbox_valid
assert record.west_lon_degree == -9.86
assert record.south_lat_degree == 41.15
assert record.east_lon_degree == 10.38
assert record.north_lat_degree == 51.56
assert record.area_name == 'France'
assert record.projection_method == 'Lambert Conic Conformal (2SP)'
found = True
assert found
elif output_format == '-postgis':
authorities = ['EPSG', 'ESRI']
else:
authorities = ['EPSG', 'ESRI', 'IGNF']
set_srid = set()
for authority in authorities:
if authority in ('EPSG', 'ESRI'):
set_codes_geographic = set()
set_codes_geographic_3d = set()
set_codes_projected = set()
set_codes_geocentric = set()
set_codes_compound = set()
set_deprecated = set()
for crs_info in osr.GetCRSInfoListFromDatabase(authority):
code = int(crs_info.code)
if crs_info.type == osr.OSR_CRS_TYPE_COMPOUND:
set_codes_compound.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOGRAPHIC_3D:
set_codes_geographic_3d.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOGRAPHIC_2D:
set_codes_geographic.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_PROJECTED:
set_codes_projected.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOCENTRIC:
set_codes_geocentric.add(code)
if crs_info.deprecated:
set_deprecated.add(code)
def main(argv):
output_format = '-pretty_wkt'
authority = None
argv = gdal.GeneralCmdLineProcessor(argv)
if argv is None:
sys.exit(0)
# Parse command line arguments.
i = 1
while i < len(argv):
arg = argv[i]
if arg == '-wkt' or arg == '-pretty_wkt' or arg == '-proj4' \
or arg == '-postgis' or arg == '-xml' or arg == '-copy':
output_format = arg
elif arg == '-authority':
i = i + 1
authority = argv[i]
elif arg[0] == '-':
Usage()
else:
Usage()
i = i + 1
# Output BEGIN transaction for PostGIS
if output_format == '-postgis':
print('BEGIN;')
# loop over all codes to generate output
if authority:
authorities = [ authority ]
elif output_format == '-postgis' :
authorities = [ 'EPSG', 'ESRI' ]
else:
authorities = [ 'EPSG', 'ESRI', 'IGNF' ]
set_srid = set()
for authority in authorities:
if authority in ('EPSG', 'ESRI'):
set_codes_geographic = set()
set_codes_geographic_3d = set()
set_codes_projected = set()
set_codes_geocentric = set()
set_codes_compound = set()
set_deprecated = set()
for crs_info in osr.GetCRSInfoListFromDatabase(authority):
code = int(crs_info.code)
if crs_info.type == osr.OSR_CRS_TYPE_COMPOUND:
set_codes_compound.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOGRAPHIC_3D:
set_codes_geographic_3d.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOGRAPHIC_2D:
set_codes_geographic.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_PROJECTED:
set_codes_projected.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOCENTRIC:
set_codes_geocentric.add(code)
if crs_info.deprecated:
set_deprecated.add(code)
set_codes_geographic = sorted(set_codes_geographic)
set_codes_geographic_3d = sorted(set_codes_geographic_3d)
set_codes_projected = sorted(set_codes_projected)
set_codes_geocentric = sorted(set_codes_geocentric)
set_codes_compound = sorted(set_codes_compound)
for typestr, set_codes in (('Geographic 2D CRS', set_codes_geographic),
('Projected CRS', set_codes_projected),
('Geocentric CRS', set_codes_geocentric),
('Compound CRS', set_codes_compound),
('Geographic 3D CRS', set_codes_geographic_3d)):
if set_codes and output_format == '-postgis':
print('-' * 80)
print('--- ' + authority + ' ' + typestr)
print('-' * 80)
for code in set_codes:
srs = osr.SpatialReference()
srs.SetFromUserInput(authority + ':' + str(code))
deprecated = False
if code in set_deprecated:
deprecated = True
trHandleCode(set_srid, srs, authority, str(code), deprecated, output_format)
else:
for crs_info in osr.GetCRSInfoListFromDatabase(authority):
srs = osr.SpatialReference()
srs.SetFromUserInput(authority + ':' + crs_info.code)
trHandleCode(set_srid, srs, authority, crs_info.code, crs_info.deprecated, output_format)
# Output COMMIT transaction for PostGIS
if output_format == '-postgis':
print('COMMIT;')
print('VACUUM ANALYZE spatial_ref_sys;')
def epsg_tr(output_format: str = '-pretty_wkt',
authority: Optional[str] = None):
# Output BEGIN transaction for PostGIS
if output_format == '-postgis':
print('BEGIN;')
# loop over all codes to generate output
if authority:
authorities = [authority]
elif output_format == '-postgis':
authorities = ['EPSG', 'ESRI']
else:
authorities = ['EPSG', 'ESRI', 'IGNF']
set_srid = set()
for authority in authorities:
if authority in ('EPSG', 'ESRI'):
set_codes_geographic = set()
set_codes_geographic_3d = set()
set_codes_projected = set()
set_codes_geocentric = set()
set_codes_compound = set()
set_deprecated = set()
for crs_info in osr.GetCRSInfoListFromDatabase(authority):
code = int(crs_info.code)
if crs_info.type == osr.OSR_CRS_TYPE_COMPOUND:
set_codes_compound.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOGRAPHIC_3D:
set_codes_geographic_3d.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOGRAPHIC_2D:
set_codes_geographic.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_PROJECTED:
set_codes_projected.add(code)
elif crs_info.type == osr.OSR_CRS_TYPE_GEOCENTRIC:
set_codes_geocentric.add(code)
if crs_info.deprecated:
set_deprecated.add(code)
set_codes_geographic = sorted(set_codes_geographic)
set_codes_geographic_3d = sorted(set_codes_geographic_3d)
set_codes_projected = sorted(set_codes_projected)
set_codes_geocentric = sorted(set_codes_geocentric)
set_codes_compound = sorted(set_codes_compound)
for typestr, set_codes in (('Geographic 2D CRS',
set_codes_geographic),
('Projected CRS', set_codes_projected),
('Geocentric CRS',
set_codes_geocentric),
('Compound CRS', set_codes_compound),
('Geographic 3D CRS',
set_codes_geographic_3d)):
if set_codes and output_format == '-postgis':
print('-' * 80)
print('--- ' + authority + ' ' + typestr)
print('-' * 80)
for code in set_codes:
srs = osr.SpatialReference()
srs.SetFromUserInput(authority + ':' + str(code))
deprecated = False
if code in set_deprecated:
deprecated = True
trHandleCode(set_srid, srs, authority, str(code),
deprecated, output_format)
else:
for crs_info in osr.GetCRSInfoListFromDatabase(authority):
srs = osr.SpatialReference()
srs.SetFromUserInput(authority + ':' + crs_info.code)
trHandleCode(set_srid, srs, authority, crs_info.code,
crs_info.deprecated, output_format)
# Output COMMIT transaction for PostGIS
if output_format == '-postgis':
print('COMMIT;')
print('VACUUM ANALYZE spatial_ref_sys;')
return 0