ZEMĚMĚŘICKÝ ÚŘADZEMĚMĚŘICKÝ ÚŘAD

Ing. Pavel Šidlichovský

INSPIRE conference 2016

Ortoimagery and Elevation coverages forINSPIRE

www.cuzk.cz

• Description of source data (OI, „elevation dataset“ )• Licenses, prices• View service OI• View service EL• Alternative to view service EL• Predefined coverages (OI, EL)• Difficulties in setting up WCS• Conclusions

2

Content

www.cuzk.cz 3

Airborne laserscanning of CZ territory

•Project in 2008•Scanning 22.6.2010‐11.11.2013•Postprocessing:  until july 2016

Figure:  © RIEGL, GmbH,  Austria

BASIC PARAMETRS 

• level of flight (AGL) maximum 1500 m (1200)

• distance of flight lines 830 m

• width of strips about 1400 m 

• overlap of the strips 35 ‐ 50%

• minimum density 1ps/m2

• average density 1,5 ps/m2  

• RMSE of vertical measurement < 0,10 m

www.cuzk.cz 4

ELEVATION PRODUCT LINE

DTM 4G DTM 5G DSM 1G

4

1. Digital Terrain Model 4G (grid 5 x 5 m) with RMSE = 0.30 m in uncovered areas and1 m in wooded and shrubed areasTerm: finished in 2013

2. Digital Terrain Model 5G (TIN form) with RMSE = 0,18 m in uncovered areas and0.30 m in wooded and shrubed areasTerm: 30.6.2016

3. Digital Surface Model 1G (TIN form) with RMSE = 0.4 m for well shaped feature and 0.7 m for other featuresTerm: finished just recently

www.cuzk.cz

• 1998: B&W  first tests• 1999‐2001: B&W production run• 2002: color tests• 2003‐2011:  33% of CZ area per 

year (three “coverages”)• 2012‐2016: 50% per year, 

transition to full digital

• 2016: acquisition of ADS80 Leicaline scanner – new challenges

Digital Orthoimagery 1998‐2016

www.cuzk.cz 6

CURRENT ORTOPHOTO DATASET

PARAMETERS OF ORTHOPHOTOS

• size of pixel 0,25 m• high color quality• RMSE = 0,3 m in horizontal plane• two years periodicity • national JTSK & UTM/ETRS89 CRS• metadata according ISO 19115 are 

provided through the Geoportal ČÚZK

PARAMETERS OF DIGITAL IMAGES

• using digital format cameras• georeferencing using GPS, IMU and

GCP   • scale of images 1: 34700• frame of images 103,9 x 67,8 mm • covered area 3608 x 2354 m• size of pixel Ø 0,2 m in ground• RGB and NIR in TIFF format

www.cuzk.cz

Data access policy for OI, EL datasets

(In fact its similar for most datasets of LSU)

• View services are FREE, non registered• Data and services are available to goverment and municipal authorities

FREE of charge , based on licence agreement• Both download services and access to predefined datasets are provided

for a fee

DMR 4G or INSPIRE EL:        cca. 1,0 EUR/km2

DMR 5G  cca. 1,2 EUR/km2

DMP 1G  cca. 1,5EUR/km2

OROIMAGERY cca. 1,1EUR/km2

www.cuzk.cz

OI view service

www.cuzk.cz

Native CRS :national : S‐JTSK         INSPIRE: ETRS‐89/TM33N 

Usage: (eg. August 2016)national:  32 milion requests INSPIRE:  13 225 requests

9

Main differences in OI view services

www.cuzk.cz 10

Ortophoto mosaic ‐ seamlines

Mosaic seamlines are available inside production workflow but not pulled to publication

www.cuzk.cz 11

Traditional 2D view service for elevation

www.cuzk.cz 12

INSPIRE view service

www.cuzk.cz 13

„Elevation Analysis“  application

Initial capability:                      http://ags.cuzk.cz/dmr/#    

Based on „Image“ ArcGIS services + WebAppBuilder

www.cuzk.cz 14

Future development of DTM Analysis:LoS

www.cuzk.cz 15

Future development of DTM Analysis: Visibility

www.cuzk.cz 16

INSPIRE  ELEVATION‐ COVERAGE as predefined dataset

GML „envelope“ + TIFF LZW compressed

Provided in ETRS‐89/TM33N  planar CS

Tiles according to INSPIRE GRID ETRS89‐GRS80(2 sizes ) 

Irregular image area

www.cuzk.cz 17

INSPIRE Ortoimagery as predefined dataset

www.cuzk.cz

IR RequirementAnnex I, Section 6.2.1.2 Rectified Grid Coverage (RectifiedGridCoverage) • (…) • Grid points of a RectifiedGridCoverage shall coincide with the centres of cells of the 

geographical grids defined in Section 2.2 of Annex II of the IR at any resolution level.

Annex III, Section 3.5.2 Requirements for Orthoimage Coverages(1) By way of derogation from the requirement in Section 2.2 of Annex II, any grid 

compatible with one of the following coordinate reference systems may be used for making gridded Orthoimagery data available: 

• ….• ….• – plane coordinates using the ETRS89 Transverse Mercator coordinate reference 

system.

18

Is it legal ?

www.cuzk.cz 19

Grid_ETRS89‐GRS80

www.cuzk.cz

• Describe Spatial Object Type (IR 1088/2010) – should be supported for direct download servicesOI:  ( Coverage, Single Mosaic Element, Agregated Mosaic Element)EL:  (Elevation Grid Coverage)‐ How to map on WCS operations?

• Mosaic Element subtypes:  set up separate WFS service?

• Aditional Coverage attributes:Store in /gml‐cov:metadata or /gmlcov:Extension ?

20

IR and TG  requirement for download service

www.cuzk.cz

/el‐cov:beginLifespanVersion Constant value

/el‐cov:domainExtent/gml:EX_Extent/../gml:posList

reuse bbox ‐ extent of subset

/el‐cov:endLifespanVersion No meaning unless archived versions are provided

/base:localId Spatial Object = Coverage,  localId= ID of coverage?

/base:namespace Constant value

/el‐cov:propertyType Constant value

/el‐cov:surfaceType Constant value

21

INSPIRE attribute extension to EL coverages:

All could be done by XSLT transformation of output GML of WCS service.

www.cuzk.cz

/oi:beginLifespanVersion Constant value

/oi:domainExtent/gml:EX_Extent/../gml:posList reuse bbox ‐ extent of subset

/oi:footprint Has to be computed  from  request extent and mosaic elements / coverage footprint

/oi:endLifespanVersion No meaning unless archive versions are provided / still constant value

/base:localId Spatial Object = Coverage,  localId= ID of coverage?

/base:namespace Constant value

/oi:interpolationType Constant value

/oi:phenomenonTime/gml:begin/gml:TimeInstant/gml:timePosition

Given difficulties with “Mosaic Elements”It`s tempting to provide constant range for coverage/oi:phenomenonTime/gml:end/gml:TimeInstant

/gml:timePosition

22

INSPIRE attribute extension to OI coverages:

www.cuzk.cz

• What is “coverage” ?  • Size of LZW compressed TIFF OI coverage is ~ 3TB

• Can we use ECW then?   

• Service cannot be free:• Autentization• Spatial and temporal restriction

• Limit on GB per single request, similar to WFS restrictions• Cannot offer transformation 

23

Aditional :national: requirements for WCS 

= full country

( while national dataset source data in ECW format ~ 600 GB)No, only JPEG2000

www.cuzk.cz

• We cannot do WCS as INSPIRE download  service without substantial development (middleware between WCS and user)

• Not clear “added value” of INSPIRE extensions to justify the development

• Wait for TG for WCS

• Low usage of INSPIRE services, which do not encourage investments

24

Conclusions

Anybody has “out of the box” solution for WCS ?