Simulation of spatial sensor characteristics in the context of the EnMAP Hyperspectral mission
(2010)
The simulation of remote sensing images is a valuable tool for defining future Earth observation systems, optimizing instrument parameters, and developing and validating data-processing algorithms. A scene simulator for optical Earth observation data has been developed within the Environmental Mapping and Analysis Program (EnMAP) hyperspectral mission. It produces EnMAP-like data following a sequential processing approach consisting of five independent modules referred to as reflectance, atmospheric, spatial, spectral, and radiometric modules. From a modeling viewpoint, the spatial module is the most complex. The spatial simulation process considers the satellite-target geometry, which is adapted to the EnMAP orbit and operating characteristics, the instrument spatial response, and the sources of spatial nonuniformity (keystone, telescope distortion and smile, and detector coregistration). The spatial module of the EnMAP scene simulator is presented in this paper. The EnMAP spatial and geometric characteristics will be described, the simulation methodology will be presented in detail, and the capability of the EnMAP simulator will be shown by illustrative examples.
Abbildende Spektrometrie
(2006)
Hyperspectral REE (Rare Earth Element) Mapping of Outcrops-Applications for Neodymium Detection
(2015)
In this study, an in situ application for identifying neodymium (Nd) enriched surface materials that uses multitemporal hyperspectral images is presented (HySpex sensor). Because of the narrow shape and shallow absorption depth of the neodymium absorption feature, a method was developed for enhancing and extracting the necessary information for neodymium from image spectra, even under illumination conditions that are not optimal. For this purpose, the two following approaches were developed: (1) reducing noise and analyzing changing illumination conditions by averaging multitemporal image scenes and (2) enhancing the depth of the desired absorption band by deconvolving every image spectrum with a Gaussian curve while the rest of the spectrum remains unchanged (Richardson-Lucy deconvolution). To evaluate these findings, nine field samples from the Fen complex in Norway were analyzed using handheld X-ray fluorescence devices and by conducting detailed laboratory-based geochemical rare earth element determinations. The result is a qualitative outcrop map that highlights zones that are enriched in neodymium. To reduce the influences of non-optimal illumination, particularly at the studied site, a minimum of seven single acquisitions is required. Sharpening the neodymium absorption band allows for robust mapping, even at the outer zones of enrichment. From the geochemical investigations, we found that iron oxides decrease the applicability of the method. However, iron-related absorption bands can be used as secondary indicators for sulfidic ore zones that are mainly enriched with rare earth elements. In summary, we found that hyperspectral spectroscopy is a noninvasive, fast and cost-saving method for determining neodymium at outcrop surfaces.
Knowledge of the origin of suspended sediment is important for improving our understanding of sediment dynamics and thereupon support of sustainable watershed management. An direct approach to trace the origin of sediments is the fingerprinting technique. It is based on the assumption that potential sediment sources can be discriminated and that the contribution of these sources to the sediment can be determined on the basis of distinctive characteristics (fingerprints). Recent studies indicate that visible-near-infrared (VNIR) and shortwave-infrared (SWIR) reflectance characteristics of soil may be a rapid, inexpensive alternative to traditional fingerprint properties (e.g. geochemistry or mineral magnetism).
To further explore the applicability of VNIR-SWIR spectral data for sediment tracing purposes, source samples were collected in the Isabena watershed, a 445 km(2) dryland catchment in the central Spanish Pyrenees. Grab samples of the upper soil layer were collected from the main potential sediment source types along with in situ reflectance spectra. Samples were dried and sieved, and artificial mixtures of known proportions were produced for algorithm validation. Then, spectral readings of potential source and artificial mixture samples were taken in the laboratory. Colour coefficients and physically based parameters were calculated from in situ and laboratory-measured spectra. All parameters passing a number of prerequisite tests were subsequently applied in discriminant function analysis for source discrimination and mixing model analyses for source contribution assessment.
The three source types (i.e. badlands, forest/grassland and an aggregation of other sources, including agricultural land, shrubland, unpaved roads and open slopes) could be reliably identified based on spectral parameters. Laboratory-measured spectral fingerprints permitted the quantification of source contribution to artificial mixtures, and introduction of source heterogeneity into the mixing model decreased accuracies for some source types. Aggregation of source types that could not be discriminated did not improve mixing model results. Despite providing similar discrimination accuracies as laboratory source parameters, in situ derived source information was found to be insufficient for contribution modelling.
The laboratory mixture experiment provides valuable insights into the capabilities and limitations of spectral fingerprint properties. From this study, we conclude that combinations of spectral properties can be used for mixing model analyses of a restricted number of source groups, whereas more straightforward in situ measured source parameters do not seem suitable. However, modelling results based on laboratory parameters also need to be interpreted with care and should not rely on the estimates of mean values only but should consider uncertainty intervals as well.
Remote sensing analysis is a crucial tool for monitoring the extent of mine waste surfaces and their mineralogy in countries with a long mining history, such as South Africa, where gold and platinum have been produced for over 90 years. These mine waste sites have the potential to contain problematic trace element species (e. g., U, Pb, Cr). In our research, we aim to combine the mapping and monitoring capacities of multispectral and hyperspectral spaceborne sensors. This is done to assess the potential of existing multispectral and hyperspectral spaceborne sensors (OLI and Hyperion) and future missions, such as Sentinel-2 and EnMAP (Environmental Mapping and Analysis Program), for mapping the spatial extent of these mine waste surfaces. For this task we propose a new index, termed the iron feature depth (IFD), derived from Landsat-8 OLI data to map the 900-nm absorption feature as a potential proxy for monitoring the spatial extent of mine waste. OLI was chosen, because it represents the most suitable sensor to map the IFD over large areas in a multi-temporal manner due to its spectral band layout; its (183 km x 170 km) scene size and its revisiting time of 16 days. The IFD is in good agreement with primary and secondary iron-bearing minerals mapped by the Material Identification and Characterization Algorithm (MICA) from EO-1 Hyperion data and illustrates that a combination of hyperspectral data (EnMAP) for mineral identification with multispectral data (Sentinel-2) for repetitive area-wide mapping and monitoring of the IFD as mine waste proxy is a promising application for future spaceborne sensors. A maximum, absolute model error is used to assess the ability of existing and future multispectral sensors to characterize mine waste via its 900-nm iron absorption feature. The following sensor-signal similarity ranking can be established for spectra from gold mining material: EnMAP 100% similarity to the reference, ALI 97.5%, Sentinel-2 97%, OLI and ASTER 95% and ETM+ 91% similarity.
Modern imaging spectrometers produce an ever-growing amount of data, which increases the need for automated analysis techniques. The algorithms employed, such as the United States Geological Survey (USGS) Tetracorder and the Mineral Identification and Characterization Algorithm (MICA), use a standardized spectral library and expert knowledge for the detection of surface cover types. Correct absorption feature definition and isolation are key to successful material identification using these algorithms. Here, a new continuum removal and feature isolation technique is presented, named the 'Geometric Hull Technique'. It is compared to the well-established, knowledge-based Tetracorder feature database together with the adapted state of the art techniques scale-space filtering, alpha shapes and convex hull.
The results show that the geometric hull technique yields the smallest deviations from the feature definitions of the MICA Group 2 library with a median difference of only 8nm for the position of the features and a median difference of only 15% for the feature shapes. The modified scale-space filtering hull technique performs second best with a median feature position difference of 16nm and a median difference of 25% for the feature shapes. The scale-space alpha hull technique shows a 23nm median position difference and a median deviation of 77% for the feature shapes. The geometric hull technique proposed here performs best amongst the four feature isolation techniques and may be an important building block for next generation automatic mapping algorithms.
In this study, an in situ application for identifying neodymium (Nd) enriched surface materials that uses multitemporal hyperspectral images is presented (HySpex sensor). Because of the narrow shape and shallow absorption depth of the neodymium absorption feature, a method was developed for enhancing and extracting the necessary information for neodymium from image spectra, even under illumination conditions that are not optimal. For this purpose, the two following approaches were developed: (1) reducing noise and analyzing changing illumination conditions by averaging multitemporal image scenes and (2) enhancing the depth of the desired absorption band by deconvolving every image spectrum with a Gaussian curve while the rest of the spectrum remains unchanged (Richardson-Lucy deconvolution). To evaluate these findings, nine field samples from the Fen complex in Norway were analyzed using handheld X-ray fluorescence devices and by conducting detailed laboratory-based geochemical rare earth element determinations. The result is a qualitative outcrop map that highlights zones that are enriched in neodymium. To reduce the influences of non-optimal illumination, particularly at the studied site, a minimum of seven single acquisitions is required. Sharpening the neodymium absorption band allows for robust mapping, even at the outer zones of enrichment. From the geochemical investigations, we found that iron oxides decrease the applicability of the method. However, iron-related absorption bands can be used as secondary indicators for sulfidic ore zones that are mainly enriched with rare earth elements. In summary, we found that hyperspectral spectroscopy is a noninvasive, fast and cost-saving method for determining neodymium at outcrop surfaces
Discovery of high and very high-energy emission from the BL Lacertae object SHBL J001355.9-185406
(2013)
The detection of the high-frequency peaked BL Lac object (HBL) SHBL J001355.9-185406 (z = 0.095) at high (HE; 100 MeV < E < 300 GeV) and very high-energy (VHE; E > 100 GeV) with the Fermi Large Area Telescope (LAT) and the High Energy Stereoscopic System (H.E.S.S.) is reported. Dedicated observations were performed with the H. E. S. S. telescopes, leading to a detection at the 5.5 sigma significance level. The measured flux above 310 GeV is (8.3 +/- 1.7(stat) +/- 1.7(sys)) x 10(-13) photons cm(-2) s(-1) (about 0.6% of that of the Crab Nebula), and the power-law spectrum has a photon index of Gamma = 3.4 +/- 0.5(stat) +/- 0.2(sys). Using 3.5 years of publicly available Fermi-LAT data, a faint counterpart has been detected in the LAT data at the 5.5 sigma significance level, with an integrated flux above 300 MeV of (9.3 +/- 3.4(stat) +/- 0.8(sys)) x 10(-10) photons cm(-2) s(-1) and a photon index of Gamma = 1.96 +/- 0.20(stat) +/- 0.08(sys). X-ray observations with Swift-XRT allow the synchrotron peak energy in vF(v) representation to be located at similar to 1.0 keV. The broadband spectral energy distribution is modelled with a one-zone synchrotron self-Compton (SSC) model and the optical data by a black-body emission describing the thermal emission of the host galaxy. The derived parameters are typical of HBLs detected at VHE, with a particle-dominated jet.
A deep observation campaign carried out by the High Energy Stereoscopic System (HESS) on Centaurus A enabled the discovery of gamma-rays from the blazar 1ES 1312-423, 2 degrees away from the radio galaxy. With a differential flux at 1 TeV of phi(1 TeV) = (1.9 +/- 0.6(stat) +/- 0.4(sys)) x 10(-13) cm(-2) s(-1) TeV-1 corresponding to 0.5 per cent of the Crab nebula differential flux and a spectral index Gamma = 2.9 +/- 0.5(stat) +/- 0.2(sys), 1ES 1312-423 is one of the faintest sources ever detected in the very high energy (E > 100 GeV) extragalactic sky. A careful analysis using three and a half years of Fermi Large Area Telescope (Fermi-LAT) data allows the discovery at high energies (E > 100 MeV) of a hard spectrum (Gamma = 1.4 +/- 0.4(stat) +/- 0.2(sys)) source coincident with 1ES 1312-423. Radio, optical, UV and X-ray observations complete the spectral energy distribution of this blazar, now covering 16 decades in energy. The emission is successfully fitted with a synchrotron self-Compton model for the non-thermal component, combined with a blackbody spectrum for the optical emission from the host galaxy.
Context. Globular clusters (GCs) are established emitters of high-energy (HE, 100 MeV < E < 100 GeV) gamma-ray radiation which could originate from the cumulative emission of the numerous millisecond pulsars (msPSRs) in the clusters’ cores or from inverse Compton (IC) scattering of relativistic leptons accelerated in the GC environment. These stellar clusters could also constitute a new class of sources in the very-high-energy (VHE, E > 100 GeV) gamma-ray regime, judging from the recent detection of a signal from the direction of Terzan 5 with the H.E.S.S. telescope array. Aims. To search for VHE gamma-ray sources associated with other GCs, and to put constraints on leptonic emission models, we systematically analyzed the observations towards 15 GCs taken with the H. E. S. S. array of imaging atmospheric Cherenkov telescopes. Methods. We searched for point-like and extended VHE gamma-ray emission from each GC in our sample and also performed a stacking analysis combining the data from all GCs to investigate the hypothesis of a population of faint emitters. Assuming IC emission as the origin of the VHE gamma-ray signal from the direction of Terzan 5, we calculated the expected gamma-ray flux from each of the 15 GCs, based on their number of millisecond pulsars, their optical brightness and the energy density of background photon fields. Results. We did not detect significant VHE gamma-ray emission from any of the 15 GCs in either of the two analyses. Given the uncertainties related to the parameter determinations, the obtained flux upper limits allow to rule out the simple IC/msPSR scaling model for NGC6388 and NGC7078. The upper limits derived from the stacking analyses are factors between 2 and 50 below the flux predicted by the simple leptonic scaling model, depending on the assumed source extent and the dominant target photon fields. Therefore, Terzan 5 still remains exceptional among all GCs, as the VHE gamma-ray emission either arises from extra-ordinarily efficient leptonic processes, or from a recent catastrophic event, or is even unrelated to the GC itself.
Axionlike particles (ALPs) are hypothetical light (sub-eV) bosons predicted in some extensions of the Standard Model of particle physics. In astrophysical environments comprising high-energy gamma rays and turbulent magnetic fields, the existence of ALPs can modify the energy spectrum of the gamma rays for a sufficiently large coupling between ALPs and photons. This modification would take the form of an irregular behavior of the energy spectrum in a limited energy range. Data from the H. E. S. S. observations of the distant BL Lac object PKS 2155 - 304 (z = 0.116) are used to derive upper limits at the 95% C. L. on the strength of the ALP coupling to photons, g(gamma a) < 2.1 x 10(-11) GeV-1 for an ALP mass between 15 and 60 neV. The results depend on assumptions on the magnetic field around the source, which are chosen conservatively. The derived constraints apply to both light pseudoscalar and scalar bosons that couple to the electromagnetic field.
The gamma-ray spectrum of the low-frequency-peaked BL Lac (LBL) object AP Librae is studied, following the discovery of very-high-energy (VHE; E > 100 GeV) gamma-ray emission up to the TeV range by the H.E.S.S. experiment. Thismakes AP Librae one of the few VHE emitters of the LBL type. The measured spectrum yields a flux of (8.8 +/- 1.5(stat) +/- 1.8(sys)) x 10(-12) cm(-2) s(-1) above 130 GeV and a spectral index of Gamma = 2.65 +/- 0.19(stat) +/- 0.20(sys). This study also makes use of Fermi-LAT observations in the high energy (HE, E > 100 MeV) range, providing the longest continuous light curve (5 years) ever published on this source. The source underwent a flaring event between MJD 56 306-56 376 in the HE range, with a flux increase of a factor of 3.5 in the 14 day bin light curve and no significant variation in spectral shape with respect to the low-flux state. While the H.E.S.S. and (low state) Fermi-LAT fluxes are in good agreement where they overlap, a spectral curvature between the steep VHE spectrum and the Fermi-LAT spectrum is observed. The maximum of the gamma-ray emission in the spectral energy distribution is located below the GeV energy range.
Aims. Previous observations with the High Energy Stereoscopic System (H.E.S.S.) have revealed an extended very-high-energy (VHE; E > 100 GeV) gamma-ray source, HESS J1834-087, coincident with the supernova remnant (SNR) W41. The origin of the gamma-ray emission was investigated in more detail with the H.E.S.S. array and the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope.
Methods. The gamma-ray data provided by 61 h of observations with H.E.S.S., and four years with the Fermi LAT were analyzed, covering over five decades in energy from 1.8 GeV up to 30 TeV. The morphology and spectrum of the TeV and GeV sources were studied and multiwavelength data were used to investigate the origin of the gamma-ray emission toward W41.
Results. The TeV source can be modeled with a sum of two components: one point-like and one significantly extended (sigma(TeV) = 0.17 degrees +/- 0.01 degrees), both centered on SNR W41 and exhibiting spectra described by a power law with index Gamma(TeV) similar or equal to 2.6. The GeV source detected with Fermi LAT is extended (sigma(GeV) = 0.15 degrees +/- 0.03 degrees) and morphologically matches the VHE emission. Its spectrum can be described by a power-law model with an index Gamma(GeV) = 2.15 +/- 0.12 and smoothly joins the spectrum of the whole TeV source. A break appears in the gamma-ray spectra around 100 GeV. No pulsations were found in the GeV range.
Conclusions. Two main scenarios are proposed to explain the observed emission: a pulsar wind nebula (PWN) or the interaction of SNR W41 with an associated molecular cloud. X-ray observations suggest the presence of a point-like source (a pulsar candidate) near the center of the remnant and nonthermal X-ray diffuse emission that could arise from the possibly associated PWN. The PWN scenario is supported by the compatible positions of the TeV and GeV sources with the putative pulsar. However, the spectral energy distribution from radio to gamma-rays is reproduced by a one-zone leptonic model only if an excess of low-energy electrons is injected following a Maxwellian distribution by a pulsar with a high spin-down power (> 10(37) erg s(-1)). This additional low-energy component is not needed if we consider that the point-like TeV source is unrelated to the extended GeV and TeV sources. The interacting SNR scenario is supported by the spatial coincidence between the gamma-ray sources, the detection of OH (1720 MHz) maser lines, and the hadronic modeling.
Context. Puppis A is an interesting similar to 4 kyr-old supernova remnant (SNR) that shows strong evidence of interaction between the forward shock and a molecular cloud. It has been studied in detail from radio frequencies to high-energy (HE, 0.1-100 GeV) gamma-rays. An analysis of the Fermi-LAT data has shown extended HE gamma-ray emission with a 0.2-100 GeV spectrum exhibiting no significant deviation from a power law, unlike most of the GeV-emitting SNRs known to be interacting with molecular clouds. This makes it a promising target for imaging atmospheric Cherenkov telescopes (IACTs) to probe the gamma-ray emission above 100 GeV.
Aims. Very-high-energy (VHE, E >= 0.1 TeV) gamma-ray emission from Puppis A has been, for the first time, searched for with the High Energy Stereoscopic System (HESS.).
Methods. Stereoscopic imaging of Cherenkov radiation from extensive air showers is used to reconstruct the direction and energy of the incident gamma-rays in order to produce sky images and source spectra. The profile likelihood method is applied to find constraints on the existence of a potential break or cutoff in the photon spectrum.
Results. The analysis of the HESS. data does not reveal any significant emission towards Puppis A. The derived upper limits on the differential photon flux imply that its broadband gamma-ray spectrum must exhibit a spectral break or cutoff. By combining Fermi-LAT and HESS. measurements, the 99% confidence-level upper limits on such a cutoff are found to be 450 and 280 GeV, assuming a power law with a simple exponential and a sub-exponential cutoff, respectively. It is concluded that none of the standard limitations (age, size, radiative losses) on the particle acceleration mechanism, assumed to be continuing at present, can explain the lack of VHE signal. The scenario in which particle acceleration has ceased some time ago is considered as an alternative explanation. The HE/VHE spectrum of Puppis A could then exhibit a break of non-radiative origin (as observed in several other interacting SNRs, albeit at somewhat higher energies), owing to the interaction with dense and neutral material, in particular towards the NE region.
G349.7+0.2 is a young Galactic supernova remnant (SNR) located at the distance of 11.5 kpc and observed across the entire electromagnetic spectrum from radio to high energy (HE; 0.1 GeV < E < 100 GeV) gamma-rays. Radio and infrared observations indicate that the remnant is interacting with a molecular cloud. In this paper, the detection of very high energy (VHE, E > 100 GeV) gamma-ray emission coincident with this SNR with the High Energy Stereoscopic System (HESS.) is reported. This makes it one of the farthest Galactic SNR ever detected in this domain. An integral flux F(E > 400 GeV) = (6.5 +/- 1.1(stat) +/- 1.3(syst)) x 10-11 ph cm(-2) s(-1) corresponding to similar to 0.7% of that of the Crab Nebula and to a luminosity of similar to 10(34) erg s(-1) above the same energy threshold, and a steep photon index Gamma(VHE) = 2.8 +/- 0.27(stat) +/- 0.20(syst) are measured. The analysis of more than 5 yr of Fermi-LAT data towards this source shows a power-law like spectrum with a best-fit photon index Gamma(HE) = 2.2 +/- 0.04.2(stat-0.31sys)(+0.13), The combined gamma-ray spectrum of 0349.7+0.2 can be described by either a broken power law (I3PL) or a power law with exponential (or sub exponential) cutoff (PLC). In the former case, the photon break energy is found at E-br,E-gamma = 551(-30)(+70) GeV, slightly higher than what is usually observed in the HE/VHE gamma-ray emitting middle-aged SNRs known to be interacting with molecular clouds. In the latter case. the exponential (respectively sub-exponential) cutoff energy is measured at E-cat,E-gamma = 1.4(-0.55)(+1.6) (respectively 0.35(-0.21)(+0.75)) TeV. A pion decay process resulting from the interaction of the accelerated protons and nuclei with the dense surrounding medium is clearly the preferred scenario to explain the gamma-ray emission. The BPL with a spectral steepening of 0.5-1 and the PLC provide equally good fits to the data. The product or the average gas density and the total energy content of accelerated protons and nuclei amounts to nu W-p similar to 5 x 10(51) erg cm(-3)
Very high energy (VHE, E > 100 GeV)gamma-ray flaring activity of the high-frequency peaked BL Lac object PG 1553 + 113 has been detected by the H.E.S.S. telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with a hint of intra-night variability. No counterpart of this event has been detected in the Fermi-Large Area Telescope data. This pattern is consistent with VHE gamma(-)ray flaring being caused by the injection of ultrarelativistic particles, emitting.-rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at z = 0.49 +/- 0.04 using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz invariance violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are E-QG,E- 1 > 4.10 x 10(17) GeV and E-QG,E- 2 > 2.10 x 10(10) GeV for linear and quadratic LIV effects, respectively.
An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated on-off observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of similar to 9 h of on-off observations. Upper limits on the velocity averaged cross section, <sigma upsilon >, for the annihilation of dark matter particles with masses in the range of similar to 300 GeV to similar to 10 TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of <sigma upsilon > that are larger than 3 x 10(-24) cm(3)/s are excluded for dark matter particles with masses between similar to 1 and similar to 4 TeV at 95% C.L. if the radius of the central dark matter density core does not exceed 500 pc. This is the strongest constraint that is derived on <sigma upsilon > for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.
The Large Magellanic Cloud, a satellite galaxy of the Milky Way, has been observed with the High Energy Stereoscopic System (H.E.S.S.) above an energy of 100 billion electron volts for a deep exposure of 210 hours. Three sources of different types were detected: the pulsar wind nebula of the most energetic pulsar known, N 157B; the radio-loud supernova remnant N 132D; and the largest nonthermal x-ray shell, the superbubble 30 Dor C. The unique object SN 1987A is, unexpectedly, not detected, which constrains the theoretical framework of particle acceleration in very young supernova remnants. These detections reveal the most energetic tip of a g-ray source population in an external galaxy and provide via 30 Dor C the unambiguous detection of g-ray emission from a superbubble.
This Letter reports the discovery of a remarkably hard spectrum source, HESS J1641-463, by the High Energy Stereoscopic System (H.E.S.S.) in the very high energy (VHE) domain. HESS J1641-463 remained unnoticed by the usual analysis techniques due to confusion with the bright nearby source HESS J1640-465. It emerged at a significance level of 8.5 standard deviations after restricting the analysis to events with energies above 4 TeV. It shows a moderate flux level of phi(E > 1TeV) = (3.64 +/- 0.44(stat)+/- 0.73(sys)) x 10(-13) cm(-2) s(-1), corresponding to 1.8% of the Crab Nebula flux above the same energy, and a hard spectrum with a photon index of Gamma = 2.07 +/- 0.11(stat)+/- 0.20(sys). It is a point-like source, although an extension up to a Gaussian width of sigma = 3 arcmin cannot be discounted due to uncertainties in the H.E.S.S. point-spread function. The VHE gamma-ray flux of HESS J1641-463 is found to be constant over the observed period when checking time binnings from the year-by-year to the 28 minute exposure timescales. HESS J1641-463 is positionally coincident with the radio supernova remnant SNR G338.5+0.1. No X-ray candidate stands out as a clear association; however, Chandra and XMM-Newton data reveal some potential weak counterparts. Various VHE gamma-ray production scenarios are discussed. If the emission from HESS J1641-463 is produced by cosmic ray protons colliding with the ambient gas, then their spectrum must extend close to 1 PeV. This object may represent a source population contributing significantly to the galactic cosmic ray flux around the knee.