Hyperion - Revision history

Neteler: use generic manual URL; add categories

2023-12-03T13:28:41Z

use generic manual URL; add categories

← Older revision		Revision as of 13:28, 3 December 2023
Line 1:		Line 1:
	== Hyperion Hyperspectral Data ==		== Hyperion Hyperspectral Data ==
	[[File:eo1_satellite.jpg]][[File:HyperionR50G35B100.png\|\|100px]]		[[File:eo1_satellite.jpg]][[File:HyperionR50G35B100.png\|\|100px]]

	=== EO-1 platform ===		=== EO-1 platform ===
	The NASA EO-1 satellite was launched on November 21, 2000 as part of a one-year technology validation/demonstration mission ([https://www.usgs.gov/centers/eros/science/earth-observing-1-eo-1 EO-1 Mission Page]).		The NASA EO-1 satellite was launched on November 21, 2000 as part of a one-year technology validation/demonstration mission ([https://www.usgs.gov/centers/eros/science/earth-observing-1-eo-1 EO-1 Mission Page]).
Line 53:		Line 54:
	Of which ''EO1H1260582015073110T1_MTL_L1GST.TXT'' is a Landsat type metadata header file (''MTL type'') well known to Landsat users.		Of which ''EO1H1260582015073110T1_MTL_L1GST.TXT'' is a Landsat type metadata header file (''MTL type'') well known to Landsat users.

	=== Data ~~perculiarity~~ ===		=== Data peculiarity ===
	It seems that bands 1-8 are filled with NODATA, the same for bands 58 to 76, and again for 225 to 242.		It seems that bands 1-8 are filled with NODATA, the same for bands 58 to 76, and again for 225 to 242.

Line 126:		Line 127:

	# CONFIGURATION OF 6S		# CONFIGURATION OF 6S
	# https://grass.osgeo.org/~~grass84~~/manuals/i.atcorr.html		# https://grass.osgeo.org/grass-stable/manuals/i.atcorr.html

	# Atmospheric mode		# Atmospheric mode
Line 442:		Line 443:

	</source>		</source>

			[[Category: Documentation‏‎]]
			[[Category: Image processing]]

⚠️Ychemin: /* Signatures preparation and launching i.spec.sam */

2023-12-02T12:15:37Z

Signatures preparation and launching i.spec.sam

← Older revision		Revision as of 12:15, 2 December 2023
Line 431:		Line 431:
	n=$(echo "$n + 1" \| bc)		n=$(echo "$n + 1" \| bc)
	done		done
			# i.spec.sam does not like a space at the end of a line
			sed -i 's/\s*$//' $sigfname

	# Define the imagery working group		# Define the imagery working group
Line 437:		Line 439:

	# Run a first set of Spectral Angle convergence		# Run a first set of Spectral Angle convergence
	i.spec.sam group=temp input=$sigfname result=s		i.spec.sam group=temp input=$sigfname result=s --verbose

	</source>		</source>

⚠️Ychemin: /* Signatures preparation and launching i.spec.sam */

2023-12-02T11:50:24Z

Signatures preparation and launching i.spec.sam

← Older revision		Revision as of 11:50, 2 December 2023
Line 368:		Line 368:

	<source lang="shell">		<source lang="shell">

	# Define the Hyperion files available for processing		# Define the Hyperion files available for processing
	listf="ref.8 ref.9 ref.10 ref.11 ref.12 ref.13 ref.14 ref.15 ref.16 ref.17 ref.18 ref.19 ref.20 ref.21 ref.22 ref.23 ref.24 ref.25 ref.26 ref.27 ref.28 ref.29 ref.30 ref.31 ref.32 ref.33 ref.34 ref.35 ref.36 ref.37 ref.38 ref.39 ref.40 ref.41 ref.42 ref.43 ref.44 ref.45 ref.46 ref.47 ref.48 ref.49 ref.50 ref.51 ref.52 ref.53 ref.54 ref.55 ref.56 ref.57 ref.77 ref.78 ref.79 ref.80 ref.81 ref.82 ref.83 ref.84 ref.85 ref.86 ref.87 ref.88 ref.89 ref.90 ref.91 ref.92 ref.93 ref.94 ref.95 ref.96 ref.97 ref.98 ref.99 ref.100 ref.101 ref.102 ref.103 ref.104 ref.105 ref.106 ref.107 ref.108 ref.109 ref.110 ref.111 ref.112 ref.113 ref.114 ref.115 ref.116 ref.117 ref.118 ref.119 ref.120 ref.121 ref.122 ref.123 ref.124 ref.125 ref.126 ref.127 ref.128 ref.129 ref.130 ref.131 ref.132 ref.133 ref.134 ref.135 ref.136 ref.137 ref.138 ref.139 ref.140 ref.141 ref.142 ref.143 ref.144 ref.145 ref.146 ref.147 ref.148 ref.149 ref.150 ref.151 ref.152 ref.153 ref.154 ref.155 ref.156 ref.157"		listf="ref.8 ref.9 ref.10 ref.11 ref.12 ref.13 ref.14 ref.15 ref.16 ref.17 ref.18 ref.19 ref.20 ref.21 ref.22 ref.23 ref.24 ref.25 ref.26 ref.27 ref.28 ref.29 ref.30 ref.31 ref.32 ref.33 ref.34 ref.35 ref.36 ref.37 ref.38 ref.39 ref.40 ref.41 ref.42 ref.43 ref.44 ref.45 ref.46 ref.47 ref.48 ref.49 ref.50 ref.51 ref.52 ref.53 ref.54 ref.55 ref.56 ref.57 ref.77 ref.78 ref.79 ref.80 ref.81 ref.82 ref.83 ref.84 ref.85 ref.86 ref.87 ref.88 ref.89 ref.90 ref.91 ref.92 ref.93 ref.94 ref.95 ref.96 ref.97 ref.98 ref.99 ref.100 ref.101 ref.102 ref.103 ref.104 ref.105 ref.106 ref.107 ref.108 ref.109 ref.110 ref.111 ref.112 ref.113 ref.114 ref.115 ref.116 ref.117 ref.118 ref.119 ref.120 ref.121 ref.122 ref.123 ref.124 ref.125 ref.126 ref.127 ref.128 ref.129 ref.130 ref.131 ref.132 ref.133 ref.134 ref.135 ref.136 ref.137 ref.138 ref.139 ref.140 ref.141 ref.142 ref.143 ref.144 ref.145 ref.146 ref.147 ref.148 ref.149 ref.150 ref.151 ref.152 ref.153 ref.154 ref.155 ref.156 ref.157"
Line 393:		Line 392:
	lno=$(cat temp \| wc -l)		lno=$(cat temp \| wc -l)
	sed -i 's/$.$,,$.$/\2/' temp		sed -i 's/$.$,,$.$/\2/' temp
	tr '\n' ' ' < ~~temp~~ > $fname		rm -f temp1
			for number in $(cat temp)
			do
			printf "%.8f " $number >> temp1
			done
			tr '\n' ' ' < temp1 > $fname
	echo "# $classname" >> $sigfname		echo "# $classname" >> $sigfname

Line 407:		Line 411:
	done		done
	sed -i 's/$.$,,$.$/\2/' temp		sed -i 's/$.$,,$.$/\2/' temp
	tr '\n' ' ' < ~~temp~~ > $fname		rm -f temp1
			for number in $(cat temp)
			do
			printf "%.8f " $number >> temp1
			done
			tr '\n' ' ' < temp1 > $fname
	echo "# $classname" >> $sigfname		echo "# $classname" >> $sigfname

Line 415:		Line 424:
	# Complete the sig file		# Complete the sig file
	echo "# " >> $sigfname		echo "# " >> $sigfname
	echo "Matrix $signo by $lno" >> $sigfname		echo "Matrix: $signo by $lno" >> $sigfname
	n=0		n=0
	for file in *.sign		for file in *.sign
Line 425:		Line 434:
	# Define the imagery working group		# Define the imagery working group
	listmaps=$(echo $listf \| sed 's/\ /,/g')		listmaps=$(echo $listf \| sed 's/\ /,/g')
	i.group group=temp input=$listmaps		i.group group=temp input=$listmaps --quiet

	# Run a first set of Spectral Angle convergence		# Run a first set of Spectral Angle convergence
	i.spec.sam group=temp input=$sigfname result=s		i.spec.sam group=temp input=$sigfname result=s
	</source>		</source>

⚠️Ychemin: /* Signatures preparation and launching i.spec.sam */

2023-12-02T11:07:57Z

Signatures preparation and launching i.spec.sam

← Older revision		Revision as of 11:07, 2 December 2023
Line 424:		Line 424:

	# Define the imagery working group		# Define the imagery working group
	~~# If already existing and loaded, then empty it~~		listmaps=$(echo $listf \| sed 's/\ /,/g')
	~~i.group -l group=temp > temp~~		i.group group=temp input=$listmaps
	~~val=$(cat temp \| wc -l)~~
	~~if [ $val -ne "" ]~~
	~~then i.group -r group=temp file=temp~~
	fi
	~~# Load it with the new set from $listf~~
	~~i.group group=temp input~~=$(echo $listf \| sed '/\ /,/')

	# Run a first set of Spectral Angle convergence		# Run a first set of Spectral Angle convergence
	i.spec.sam group=temp input=$sigfname result=s		i.spec.sam group=temp input=$sigfname result=s

	</source>		</source>

⚠️Ychemin: /* Hyperion Hyperspectral Data */

2023-12-02T11:01:01Z

Hyperion Hyperspectral Data

← Older revision		Revision as of 11:01, 2 December 2023
Line 1:		Line 1:
	== Hyperion Hyperspectral Data ==		== Hyperion Hyperspectral Data ==
	[[File:eo1_satellite.jpg]][[File:HyperionR50G35B100.png]]		[[File:eo1_satellite.jpg]][[File:HyperionR50G35B100.png\|\|100px]]
	=== EO-1 platform ===		=== EO-1 platform ===
	The NASA EO-1 satellite was launched on November 21, 2000 as part of a one-year technology validation/demonstration mission ([https://www.usgs.gov/centers/eros/science/earth-observing-1-eo-1 EO-1 Mission Page]).		The NASA EO-1 satellite was launched on November 21, 2000 as part of a one-year technology validation/demonstration mission ([https://www.usgs.gov/centers/eros/science/earth-observing-1-eo-1 EO-1 Mission Page]).

⚠️Ychemin at 11:00, 2 December 2023

2023-12-02T11:00:31Z

← Older revision		Revision as of 11:00, 2 December 2023
Line 1:		Line 1:
	== Hyperion Hyperspectral Data ==		== Hyperion Hyperspectral Data ==
	[[File:eo1_satellite.jpg]]		[[File:eo1_satellite.jpg]][[File:HyperionR50G35B100.png]]
	=== EO-1 platform ===		=== EO-1 platform ===
	The NASA EO-1 satellite was launched on November 21, 2000 as part of a one-year technology validation/demonstration mission ([https://www.usgs.gov/centers/eros/science/earth-observing-1-eo-1 EO-1 Mission Page]).		The NASA EO-1 satellite was launched on November 21, 2000 as part of a one-year technology validation/demonstration mission ([https://www.usgs.gov/centers/eros/science/earth-observing-1-eo-1 EO-1 Mission Page]).

⚠️Ychemin: Adding signatures to i.spec.sam

2023-12-02T10:56:52Z

Adding signatures to i.spec.sam

@@ Line 363: / Line 363: @@
          routgdal(input=b_out, output=HYPERIONDir+"_"+b_out+".tif", overwrite=OVR, quiet=QUIET)
          satbandno = satbandno + 1
 </source>

⚠️Ychemin: /* Import in GRASS, DN2Rad & Rad2Ref */

2023-12-02T08:51:49Z

Import in GRASS, DN2Rad & Rad2Ref

← Older revision		Revision as of 08:51, 2 December 2023
Line 76:		Line 76:
	# PURPOSE		# PURPOSE
	# This script processes HYPERION images		# This script processes HYPERION images
	# 1 - ~~untar~~ *.zip files		# 1 - unzip *.zip files
	# 2 - Connect files in GRASS GIS Location of your choice		# 2 - Connect files in GRASS GIS Location of your choice
	# 3 - DN to Radiance (Apply Gain)		# 3 - DN to Radiance (Apply Gain)
Line 331:		Line 331:
	rinfo_out = parse_key_val(rinfo.outputs.stdout)		rinfo_out = parse_key_val(rinfo.outputs.stdout)
	# Extract mean values from Dict		# Extract mean values from Dict
	minimum = rinfo_out["min"]		minimum = float(rinfo_out["min"])
	maximum = rinfo_out["max"]		maximum = float(rinfo_out["max"])
	if (minimum < 0.0):		if (minimum < 0.0):
	minimum = 0.0		minimum = 0.0

⚠️Ychemin: /* Import in GRASS, DN2Rad & Rad2Ref */

2023-12-02T08:49:27Z

Import in GRASS, DN2Rad & Rad2Ref

← Older revision		Revision as of 08:49, 2 December 2023
Line 76:		Line 76:
	# PURPOSE		# PURPOSE
	# This script processes HYPERION images		# This script processes HYPERION images
	# 1 - ~~unzip~~ *.zip files		# 1 - untar *.zip files
	# 2 - Connect files in GRASS GIS Location of your choice		# 2 - Connect files in GRASS GIS Location of your choice
	# 3 - DN to Radiance (Apply Gain)		# 3 - DN to Radiance (Apply Gain)
Line 162:		Line 162:
	for HYPERIONDir in HYPERIONDirs:		for HYPERIONDir in HYPERIONDirs:
	# Unzip all of your HYPERION images in all your directories		# Unzip all of your HYPERION images in all your directories
	print("Unzip HYPERION files in ", HYPERIONDir)		#print("Unzip HYPERION files in ", HYPERIONDir)
	p = os.system("unzip "+HYPERIONDir+".zip")		#p = os.system("unzip "+HYPERIONDir+".zip")
	# Go into the new created directory		# Go into the new created directory
	path = rsdatapath+'/'+HYPERIONDir+'/'+HYPERIONDir.split('_')[0]		path = rsdatapath+'/'+HYPERIONDir+'/'+HYPERIONDir.split('_')[0]
Line 188:		Line 188:
	rmapcalc.wait()		rmapcalc.wait()
	# Import HYPERION in GRASS		# Import HYPERION in GRASS
	count = 0		count = 1
	# To remove nodata from DN maps		# To remove nodata from DN maps
	rnull = Module('r.null')		rnull = Module('r.null')
Line 241:		Line 241:
	# HYPERION DN->Rad->Ref		# HYPERION DN->Rad->Ref
	print("Convert DN to Rad")		print("Convert DN to Rad")
	# Apply Gain ~~from file~~		# Apply Gain for VNIR
	for j in range(~~len~~(~~gain~~)):		for j in range(8, 57+1):
	exp = str("rad."+str(j+1)+"=dn."+str(j+1)+"*10.0/"+str(~~gain[j]~~))		exp = str("rad."+str(j)+"=dn."+str(j)+"*10.0/"+str(gain_vnir))
			#print(exp)
			rmapcalc(expression=exp, overwrite=OVR, quiet=QUIET)
			rmapcalc.wait()
			# Apply Gain for VNIR
			for j in range(77, 224+1):
			exp = str("rad."+str(j)+"=dn."+str(j)+"*10.0/"+str(gain_swir))
	#print(exp)		#print(exp)
	rmapcalc(expression=exp, overwrite=OVR, quiet=QUIET)		rmapcalc(expression=exp, overwrite=OVR, quiet=QUIET)
Line 267:		Line 273:
	# Create the handle for r.out.gdal		# Create the handle for r.out.gdal
	routgdal = Module('r.out.gdal')		routgdal = Module('r.out.gdal')

	# Process VNIR		# Process VNIR
	# satellite band number (HYPERION VNIR [433-481])		# satellite band number (HYPERION VNIR [433-481])
Line 280:		Line 285:
	rinfo_out = parse_key_val(rinfo.outputs.stdout)		rinfo_out = parse_key_val(rinfo.outputs.stdout)
	# Extract mean values from Dict		# Extract mean values from Dict
	minimum = rinfo_out["min"]		minimum = float(rinfo_out["min"])
	maximum = rinfo_out["max"]		maximum = float(rinfo_out["max"])
	if (minimum < 0.0):		if (minimum < 0.0):
	minimum = 0.0		minimum = 0.0
Line 358:		Line 363:
	routgdal(input=b_out, output=HYPERIONDir+"_"+b_out+".tif", overwrite=OVR, quiet=QUIET)		routgdal(input=b_out, output=HYPERIONDir+"_"+b_out+".tif", overwrite=OVR, quiet=QUIET)
	satbandno = satbandno + 1		satbandno = satbandno + 1

	</source>		</source>

⚠️Ychemin: /* Data Download */

2023-12-02T07:01:48Z

Data Download

← Older revision		Revision as of 07:01, 2 December 2023
Line 17:		Line 17:
	The EO-1 Extended Mission was chartered to collect and distribute ALI multispectral and Hyperion hyperspectral products in response to Data Acquisition Requests (DARs). Under the Extended Mission provisions, image data acquired by EO-1 were archived and distributed by the USGS Center for Earth Resources Observation and Science (EROS) and placed in the public domain.		The EO-1 Extended Mission was chartered to collect and distribute ALI multispectral and Hyperion hyperspectral products in response to Data Acquisition Requests (DARs). Under the Extended Mission provisions, image data acquired by EO-1 were archived and distributed by the USGS Center for Earth Resources Observation and Science (EROS) and placed in the public domain.

	== Data Download ==		== Hyperion Data Download ==
	Level 1 GST is terrain corrected and provided in 16-bit radiance values. The data are available in Geographic Tagged Image-File Format (GeoTIFF) and are distributed via download at no charge through either [https://earthexplorer.usgs.gov/ EarthExplorer] or [https://glovis.usgs.gov/ USGS Global Visualization Viewer (GloVis)].		Level 1 GST is terrain corrected and provided in 16-bit radiance values. The data are available in Geographic Tagged Image-File Format (GeoTIFF) and are distributed via download at no charge through either [https://earthexplorer.usgs.gov/ EarthExplorer] or [https://glovis.usgs.gov/ USGS Global Visualization Viewer (GloVis)].