def render(self, params): map_ = self.create_map() use_name = (params.version == Version(1, 0)) for coverage in params.coverages: # ReferenceableDatasets are not supported in WCS < 2.0 if issubclass(coverage.real_type, models.ReferenceableDataset): raise NoSuchCoverageException((coverage.identifier,)) data_items = self.data_items_for_coverage(coverage) native_format = self.get_native_format(coverage, data_items) layer = self.layer_for_coverage( coverage, native_format, params.version ) map_.insertLayer(layer) for outputformat in self.get_all_outputformats(not use_name): map_.appendOutputFormat(outputformat) request = ms.create_request(params) raw_result = ms.dispatch(map_, request) result = result_set_from_raw_data(raw_result) return result
def render(self, params): map_ = self.create_map() use_name = (params.version == Version(1, 0)) for coverage in params.coverages: # ReferenceableDatasets are not supported in WCS < 2.0 if coverage.grid.is_referenceable: raise NoSuchCoverageException((coverage.identifier,)) data_locations = self.arraydata_locations_for_coverage(coverage) native_format = self.get_native_format(coverage, data_locations) layer = self.layer_for_coverage( coverage, native_format, params.version ) map_.insertLayer(layer) for outputformat in self.get_all_outputformats(not use_name): map_.appendOutputFormat(outputformat) request = ms.create_request(params) raw_result = ms.dispatch(map_, request) result = result_set_from_raw_data(raw_result) # load XML using lxml # find and exclude <metadataLink> nodes if present # re-encode xml_result = etree.fromstring(result[0].data) for elem in xml_result.xpath('//*[local-name() = "metadataLink"]'): elem.getparent().remove(elem) xml_result_data = etree.tostring(xml_result, pretty_print=True, encoding='UTF-8', xml_declaration=True) result[0] = ResultBuffer(xml_result_data, result[0].content_type) return result
def render(self, params): map_ = self.create_map() use_name = (params.version == Version(1, 0)) for coverage in params.coverages: # ReferenceableDatasets are not supported in WCS < 2.0 if issubclass(coverage.real_type, models.ReferenceableDataset): raise NoSuchCoverageException((coverage.identifier, )) data_items = self.data_items_for_coverage(coverage) native_format = self.get_native_format(coverage, data_items) layer = self.layer_for_coverage(coverage, native_format, params.version) map_.insertLayer(layer) for outputformat in self.get_all_outputformats(not use_name): map_.appendOutputFormat(outputformat) request = ms.create_request(params) raw_result = ms.dispatch(map_, request) result = result_set_from_raw_data(raw_result) return result
def render(self, params): # get coverage related stuff coverage = params.coverage # ReferenceableDataset are not supported in WCS < 2.0 if issubclass(coverage.real_type, models.ReferenceableDataset): raise NoSuchCoverageException((coverage.identifier,)) data_items = self.data_items_for_coverage(coverage) range_type = coverage.range_type bands = list(range_type) subsets = params.subsets if subsets: srid = subsets.srid if srid is not None: if not crss.validateEPSGCode(srid): raise RenderException( "Failed to extract an EPSG code from the CRS URI " "'%s'." % srid, "subset" ) # create and configure map object map_ = self.create_map() # configure outputformat native_format = self.get_native_format(coverage, data_items) if get_format_by_mime(native_format) is None: native_format = "image/tiff" frmt = params.format or native_format if frmt is None: raise RenderException("Format could not be determined", "format") mime_type, frmt = split_format(frmt) imagemode = ms.gdalconst_to_imagemode(bands[0].data_type) time_stamp = datetime.now().strftime("%Y%m%d%H%M%S") basename = "%s_%s" % (coverage.identifier, time_stamp) of = create_outputformat( mime_type, frmt, imagemode, basename, getattr(params, "encoding_params", {}) ) map_.appendOutputFormat(of) map_.setOutputFormat(of) # TODO: use layer factory here layer = self.layer_for_coverage(coverage, native_format, params.version) map_.insertLayer(layer) for connector in self.connectors: if connector.supports(data_items): break else: raise OperationNotSupportedException( "Could not find applicable layer connector.", "coverage" ) try: connector.connect(coverage, data_items, layer) # create request object and dispatch it against the map request = ms.create_request( self.translate_params(params, range_type) ) request.setParameter("format", mime_type) raw_result = ms.dispatch(map_, request) finally: # perform any required layer related cleanup connector.disconnect(coverage, data_items, layer) result_set = result_set_from_raw_data(raw_result) if params.version == Version(2, 0): if getattr(params, "mediatype", None) in ("multipart/mixed", "multipart/related"): encoder = WCS20EOXMLEncoder() is_mosaic = issubclass( coverage.real_type, models.RectifiedStitchedMosaic ) if not is_mosaic: tree = encoder.alter_rectified_dataset( coverage, getattr(params, "http_request", None), etree.parse(result_set[0].data_file).getroot(), subsets.bounding_polygon(coverage) if subsets else None ) else: tree = encoder.alter_rectified_stitched_mosaic( coverage.cast(), getattr(params, "http_request", None), etree.parse(result_set[0].data_file).getroot(), subsets.bounding_polygon(coverage) if subsets else None ) result_set[0] = ResultBuffer( encoder.serialize(tree), encoder.content_type ) # "default" response return result_set
def render(self, layers, bbox, crs, size, frmt, time, elevation, styles): if not time: raise RenderException("Missing mandatory 'time' parameter.") try: time = time.value except AttributeError: raise RenderException( "Parameter 'time' must be a slice and not a range." ) llbbox = self.get_llbbox(bbox, crs) mapfile_path = get_eoxserver_config().get("wmm", "mapfile") map_ = ms.mapObj(mapfile_path) #TODO: path to map map_.setMetaData("ows_enable_request", "*") map_.setProjection("EPSG:4326") map_.imagecolor.setRGB(0, 0, 0) # set supported CRSs decoder = crss.CRSsConfigReader(get_eoxserver_config()) crss_string = " ".join( map(lambda crs: "EPSG:%d" % crs, decoder.supported_crss_wms) ) map_.setMetaData("ows_srs", crss_string) map_.setMetaData("wms_srs", crss_string) datasources = [] datasets = [] for layer_name in layers: layer = map_.getLayerByName(layer_name) if not layer: continue product = layer.metadata.get("wmm_product") filename = self.generate_filename("tif") ds = self.create_dataset( llbbox, time, elevation, size, product, filename ) datasets.append(ds) if layer.type == ms.MS_LAYER_LINE: flavor = layer.metadata.get("wmm_flavor") contour_steps = self.get_contour_intervals( flavor, llbbox, size ) filename = self.generate_filename("shp") self.generate_contours( ds, contour_steps, filename ) layer.connectiontype = ms.MS_OGR layer.connection = filename layer.data, _ = splitext(basename(filename)) datasources.append(filename) ms_request = ms.create_request(( ("service", "WMS"), ("version", "1.3.0"), ("request", "GetMap"), ("layers", ",".join(layers)), ("bbox", "%f,%f,%f,%f" % (bbox[1], bbox[0], bbox[3], bbox[2])), ("crs", crs), ("width", str(size[0])), ("height", str(size[1])), ("styles", ",".join(styles)), ("format", frmt) )) raw_result = ms.dispatch(map_, ms_request) result = result_set_from_raw_data(raw_result) shp_drv = ogr.GetDriverByName("ESRI Shapefile") # cleanup datasources and datasets for filename in datasources: shp_drv.DeleteDataSource(filename) for ds in datasets: driver = ds.GetDriver() for filename in ds.GetFileList(): os.remove(filename) return result, get_content_type(result)
def render(self, params): # get coverage related stuff coverage = params.coverage # ReferenceableDataset are not supported in WCS < 2.0 if issubclass(coverage.real_type, models.ReferenceableDataset): raise NoSuchCoverageException((coverage.identifier, )) data_items = self.data_items_for_coverage(coverage) range_type = coverage.range_type bands = list(range_type) subsets = params.subsets if subsets: subsets.srid # this automatically checks the validity # create and configure map object map_ = self.create_map() # configure outputformat native_format = self.get_native_format(coverage, data_items) if get_format_by_mime(native_format) is None: native_format = "image/tiff" frmt = params.format or native_format if frmt is None: raise RenderException("Format could not be determined", "format") mime_type, frmt = split_format(frmt) imagemode = ms.gdalconst_to_imagemode(bands[0].data_type) time_stamp = datetime.now().strftime("%Y%m%d%H%M%S") basename = "%s_%s" % (coverage.identifier, time_stamp) of = create_outputformat(mime_type, frmt, imagemode, basename, getattr(params, "encoding_params", {})) map_.appendOutputFormat(of) map_.setOutputFormat(of) # TODO: use layer factory here layer = self.layer_for_coverage(coverage, native_format, params.version) map_.insertLayer(layer) for connector in self.connectors: if connector.supports(data_items): break else: raise OperationNotSupportedException( "Could not find applicable layer connector.", "coverage") try: connector.connect(coverage, data_items, layer, {}) # create request object and dispatch it against the map request = ms.create_request( self.translate_params(params, range_type)) request.setParameter("format", mime_type) raw_result = ms.dispatch(map_, request) finally: # perform any required layer related cleanup connector.disconnect(coverage, data_items, layer, {}) result_set = result_set_from_raw_data(raw_result) if params.version == Version(2, 0): if getattr(params, "mediatype", None) in ("multipart/mixed", "multipart/related"): encoder = WCS20EOXMLEncoder() is_mosaic = issubclass(coverage.real_type, models.RectifiedStitchedMosaic) if not is_mosaic: tree = encoder.alter_rectified_dataset( coverage, getattr(params, "http_request", None), etree.parse(result_set[0].data_file).getroot(), subsets.bounding_polygon(coverage) if subsets else None) else: tree = encoder.alter_rectified_stitched_mosaic( coverage.cast(), getattr(params, "http_request", None), etree.parse(result_set[0].data_file).getroot(), subsets.bounding_polygon(coverage) if subsets else None) result_set[0] = ResultBuffer(encoder.serialize(tree), encoder.content_type) # "default" response return result_set
def render(self, params): # get coverage related stuff coverage = params.coverage # ReferenceableDataset are not supported in WCS < 2.0 if params.coverage.grid.is_referenceable and params.version: raise NoSuchCoverageException((coverage.identifier, )) data_locations = self.arraydata_locations_for_coverage(coverage) range_type = coverage.range_type bands = list(range_type) subsets = params.subsets if subsets: subsets.srid # this automatically checks the validity # create and configure map object map_ = self.create_map() env = {} for data_location in data_locations: env.update(data_location.env) gdal.set_env(env, False) # configure outputformat native_format = self.get_native_format(coverage, data_locations) if native_format and get_format_by_mime(native_format) is None: native_format = "image/tiff" frmt = params.format or native_format if frmt is None: raise RenderException("Format could not be determined", "format") mime_type, frmt = split_format(frmt) imagemode = ms.gdalconst_to_imagemode(bands[0].data_type) time_stamp = datetime.now().strftime("%Y%m%d%H%M%S") basename = "%s_%s" % (coverage.identifier, time_stamp) of = create_outputformat(mime_type, frmt, imagemode, basename, getattr(params, "encoding_params", {})) map_.appendOutputFormat(of) map_.setOutputFormat(of) # TODO: use layer factory here layer = self.layer_for_coverage(coverage, native_format, params.version) map_.insertLayer(layer) connector = get_connector_by_test(coverage, data_locations) if not connector: raise OperationNotSupportedException( "Could not find applicable layer connector.", "coverage") try: connector.connect(coverage, data_locations, layer, {}) # create request object and dispatch it against the map request = ms.create_request( self.translate_params(params, range_type)) request.setParameter("format", mime_type) raw_result = ms.dispatch(map_, request) finally: # perform any required layer related cleanup connector.disconnect(coverage, data_locations, layer, {}) result_set = result_set_from_raw_data(raw_result) if params.version == Version(2, 0): mediatype = getattr(params, "mediatype", None) if mediatype in ("multipart/mixed", "multipart/related"): with vsi.TemporaryVSIFile.from_buffer(result_set[1].data) as f: ds = gdal.Open(f.name) grid = objects.Grid.from_gdal_dataset(ds) # get the output CRS definition crs = params.outputcrs or subsets.crs or 'imageCRS' if crs == 'imageCRS': crs = coverage.grid.coordinate_reference_system grid._coordinate_reference_system = crs origin = objects.Origin.from_gdal_dataset(ds) size = [ds.RasterXSize, ds.RasterYSize] range_type = coverage.range_type if params.rangesubset: range_type = range_type.subset(params.rangesubset) coverage._grid = grid coverage._origin = origin coverage._size = size coverage._range_type = range_type if isinstance(result_set[1].filename, binary_type): file_name = result_set[1].filename.decode() else: file_name = result_set[1].filename reference = 'cid:coverage/%s' % file_name encoder = WCS20EOXMLEncoder() if not isinstance(coverage, objects.Mosaic): tree = encoder.encode_rectified_dataset( coverage, getattr(params, "http_request", None), reference, mime_type, subsets.bounding_polygon(coverage) if subsets else None) else: tree = encoder.encode_rectified_stitched_mosaic( coverage, getattr(params, "http_request", None), reference, mime_type, subsets.bounding_polygon(coverage) if subsets else None) result_set[0] = ResultBuffer(encoder.serialize(tree), encoder.content_type) # "default" response return result_set
def render(self, layers, bbox, crs, size, frmt, time, elevation, styles): if not time: raise RenderException("Missing mandatory 'time' parameter.") try: time = time.value except AttributeError: raise RenderException( "Parameter 'time' must be a slice and not a range.") llbbox = self.get_llbbox(bbox, crs) mapfile_path = get_eoxserver_config().get("wmm", "mapfile") map_ = ms.mapObj(mapfile_path) #TODO: path to map map_.setMetaData("ows_enable_request", "*") map_.setProjection("EPSG:4326") map_.imagecolor.setRGB(0, 0, 0) # set supported CRSs decoder = crss.CRSsConfigReader(get_eoxserver_config()) crss_string = " ".join( map(lambda crs: "EPSG:%d" % crs, decoder.supported_crss_wms)) map_.setMetaData("ows_srs", crss_string) map_.setMetaData("wms_srs", crss_string) datasources = [] datasets = [] for layer_name in layers: layer = map_.getLayerByName(layer_name) if not layer: continue product = layer.metadata.get("wmm_product") filename = self.generate_filename("tif") ds = self.create_dataset(llbbox, time, elevation, size, product, filename) datasets.append(ds) if layer.type == ms.MS_LAYER_LINE: flavor = layer.metadata.get("wmm_flavor") contour_steps = self.get_contour_intervals( flavor, llbbox, size) filename = self.generate_filename("shp") self.generate_contours(ds, contour_steps, filename) layer.connectiontype = ms.MS_OGR layer.connection = filename layer.data, _ = splitext(basename(filename)) datasources.append(filename) ms_request = ms.create_request( (("service", "WMS"), ("version", "1.3.0"), ("request", "GetMap"), ("layers", ",".join(layers)), ("bbox", "%f,%f,%f,%f" % (bbox[1], bbox[0], bbox[3], bbox[2])), ("crs", crs), ("width", str(size[0])), ("height", str(size[1])), ("styles", ",".join(styles)), ("format", frmt))) raw_result = ms.dispatch(map_, ms_request) result = result_set_from_raw_data(raw_result) shp_drv = ogr.GetDriverByName("ESRI Shapefile") # cleanup datasources and datasets for filename in datasources: shp_drv.DeleteDataSource(filename) for ds in datasets: driver = ds.GetDriver() for filename in ds.GetFileList(): os.remove(filename) return result, get_content_type(result)