Earth at night - Revision history

⚠️AnneGhisla: /* Import */ fixed map name according to link fix

2011-09-15T12:51:36Z

Import: fixed map name according to link fix

← Older revision		Revision as of 12:51, 15 September 2011
Line 18:		Line 18:
	# we do know exactly what it should be).		# we do know exactly what it should be).

	r.in.gdal -o -l in=~~land_ocean_ice_lights_8192~~.tif out=land_ocean_ice_lights		r.in.gdal -o -l in=land_ocean_ice_lights_2048.tif out=land_ocean_ice_lights

	for BAND in red green blue ; do		for BAND in red green blue ; do

⚠️AnneGhisla: /* Get data */ fixed broken link

2011-09-15T12:50:06Z

Get data: fixed broken link

← Older revision		Revision as of 12:50, 15 September 2011
Line 6:		Line 6:

	Download		Download
	http://~~veimages~~.gsfc.nasa.gov//1438/~~land_ocean_ice_lights_8192~~.tif		http://eoimages.gsfc.nasa.gov/ve/1438/land_ocean_ice_lights_2048.tif
	from NASA's		from NASA's
	[http://visibleearth.nasa.gov/view_set.php?categoryID=2363 Visible Earth] website.		[http://visibleearth.nasa.gov/view_set.php?categoryID=2363 Visible Earth] website.

⚠️NikosA at 00:30, 5 July 2009

2009-07-05T00:30:51Z

← Older revision		Revision as of 00:30, 5 July 2009
Line 1:		Line 1:
	== Import NASA's Night on Earth and Reproject in Polar ~~Sterographic~~ ==		== Import NASA's Night on Earth and Reproject in Polar Stereographic ==

	''This micro-tutorial requires GRASS 6.5''		''This micro-tutorial requires GRASS 6.5''

⚠️HamishBowman at 14:38, 25 May 2009

2009-05-25T14:38:06Z

← Older revision		Revision as of 14:38, 25 May 2009
Line 39:		Line 39:
	#### NORTHERN HEMISPHERE #################		#### NORTHERN HEMISPHERE #################
	# pull in grids		# pull in grids
	v.proj in=grid_15 location=~~xy_to_ll~~		v.proj in=grid_15 location=world_ll
	v.proj in=grid_30 location=~~xy_to_ll~~		v.proj in=grid_30 location=world_ll

	g.region n=10000000 s=-10000000 w=-10000000 e=10000000 res=5000 -p		g.region n=10000000 s=-10000000 w=-10000000 e=10000000 res=5000 -p
Line 48:		Line 48:
	for BAND in red green blue ; do		for BAND in red green blue ; do
	echo "Processing [$BAND] ..."		echo "Processing [$BAND] ..."
	r.proj -n in=$MAP.$BAND location=~~xy_to_ll~~ \		r.proj -n in=$MAP.$BAND location=world_ll \
	method=cubic_f memory=300		method=cubic_f memory=300
	r.mapcalc "${MAP}_$BAND = if($MAP.$BAND > 255, 255, if($MAP.$BAND < 0, 0, int($MAP.$BAND)))"		r.mapcalc "${MAP}_$BAND = if($MAP.$BAND > 255, 255, if($MAP.$BAND < 0, 0, int($MAP.$BAND)))"
Line 79:		Line 79:
	#### SOUTHERN HEMISPHERE #################		#### SOUTHERN HEMISPHERE #################
	# pull in grids		# pull in grids
	v.proj in=grid_15 location=~~xy_to_ll~~		v.proj in=grid_15 location=world_ll
	v.proj in=grid_30 location=~~xy_to_ll~~		v.proj in=grid_30 location=world_ll

	g.region n=10000000 s=-10000000 w=-10000000 e=10000000 res=5000 -p		g.region n=10000000 s=-10000000 w=-10000000 e=10000000 res=5000 -p
Line 88:		Line 88:
	for BAND in red green blue ; do		for BAND in red green blue ; do
	echo "Processing [$BAND] ..."		echo "Processing [$BAND] ..."
	r.proj -n in=$MAP.$BAND location=~~xy_to_ll~~ \		r.proj -n in=$MAP.$BAND location=world_ll \
	method=cubic_f memory=600		method=cubic_f memory=600
	r.mapcalc "${MAP}_$BAND = if($MAP.$BAND > 255, 255, if($MAP.$BAND < 0, 0, int($MAP.$BAND)))"		r.mapcalc "${MAP}_$BAND = if($MAP.$BAND > 255, 255, if($MAP.$BAND < 0, 0, int($MAP.$BAND)))"

⚠️HamishBowman: /* Reproject South side */

2009-05-25T14:24:57Z

Reproject South side

← Older revision		Revision as of 14:24, 25 May 2009
Line 83:		Line 83:

	g.region n=10000000 s=-10000000 w=-10000000 e=10000000 res=5000 -p		g.region n=10000000 s=-10000000 w=-10000000 e=10000000 res=5000 -p

	MAP=land_ocean_ice_lights		MAP=land_ocean_ice_lights

⚠️HamishBowman: /* Reproject South side */

2009-05-25T14:24:05Z

Reproject South side

← Older revision		Revision as of 14:24, 25 May 2009
Line 70:		Line 70:

	* Create new location based on the inverse of EPSG #3413		* Create new location based on the inverse of EPSG #3413
	#proj : stere		+proj : stere
	#datum : wgs84		+datum : wgs84
	#lat_0 : -90		+lat_0 : -90
	#lat_ts : -70		+lat_ts : -70
	#lon_0 : +135		+lon_0 : +135
	#k : 1		+k : 1

	#### SOUTHERN HEMISPHERE #################		#### SOUTHERN HEMISPHERE #################

⚠️HamishBowman: /* Reproject South side */

2009-05-25T14:23:27Z

Reproject South side

← Older revision		Revision as of 14:23, 25 May 2009
Line 72:		Line 72:
	#proj : stere		#proj : stere
	#datum : wgs84		#datum : wgs84
	#lat_0 : 90		#lat_0 : -90
	#lat_ts : 70		#lat_ts : -70
	#lon_0 : +135		#lon_0 : +135
	#k : 1		#k : 1

⚠️HamishBowman: /* Reproject South side */

2009-05-25T14:23:08Z

Reproject South side

← Older revision		Revision as of 14:23, 25 May 2009
Line 69:		Line 69:
	=== Reproject South side ===		=== Reproject South side ===

	~~#### SOUTHERN HEMISPHERE #################~~		* Create new location based on the inverse of EPSG #3413
	~~#create~~ new ~~loc'n~~ based on inverse of ~~epsg~~ #3413
	#proj : stere		#proj : stere
	#datum : wgs84		#datum : wgs84
Line 77:		Line 76:
	#lon_0 : +135		#lon_0 : +135
	#k : 1		#k : 1

			#### SOUTHERN HEMISPHERE #################
	# pull in grids		# pull in grids
	v.proj in=grid_15 location=xy_to_ll		v.proj in=grid_15 location=xy_to_ll

⚠️HamishBowman: /* Reproject North side */

2009-05-25T14:22:18Z

Reproject North side

← Older revision		Revision as of 14:22, 25 May 2009
Line 31:		Line 31:

	=== Reproject North side ===		=== Reproject North side ===

			* Create new location based on EPSG #3413
			# WGS 84 / NSIDC Sea Ice Polar Stereographic North
			<3413> +proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 \
			+y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs <>

	#### NORTHERN HEMISPHERE #################		#### NORTHERN HEMISPHERE #################
	~~#create new loc'n based on epsg #3413~~
	~~# WGS 84 / NSIDC Sea Ice Polar Stereographic North~~
	~~# <3413> +proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs <>~~

	# pull in grids		# pull in grids
	v.proj in=grid_15 location=xy_to_ll		v.proj in=grid_15 location=xy_to_ll

⚠️HamishBowman: /* Import */ r.in.gdal -l

2009-05-25T14:20:58Z

Import: r.in.gdal -l

← Older revision		Revision as of 14:20, 25 May 2009
Line 11:		Line 11:

	=== Import ===		=== Import ===
	# ~~create simple XY location~~		* Create a WGS84 Lat/Lon location (EPSG # 4326)
	r.in.gdal in=land_ocean_ice_lights_8192.tif out=land_ocean_ice_lights
			# import into lat/lon, overriding projection check and
			# setting dummy region bounds which fit into geographic
			# space (the file is missing georeferencing metadata but
			# we do know exactly what it should be).

			r.in.gdal -o -l in=land_ocean_ice_lights_8192.tif out=land_ocean_ice_lights

	for BAND in red green blue ; do		for BAND in red green blue ; do
Line 18:		Line 24:
	done		done

	~~#FIXME~~		g.region n=90 s=45 w=-180 e=180
	~~#hack PROJ WIND and cellhd to be 3 not 0~~
	~~# then run~~ g.~~setproj to set datum~~ == ~~wgs84 and create PROJ_INFO files~~
	~~# now this location is lat/lon~~

	~~g.region n=90 s=45 w=-180 e=180~~		# create 15° and 30° grid line vectors
	v.mkgrid map=grid_15 box=15,15 grid=12,24 pos=coor coor=-180,-90 breaks=30		v.mkgrid map=grid_15 box=15,15 grid=12,24 pos=coor coor=-180,-90 breaks=30
	v.mkgrid map=grid_30 box=30,30 grid=6,12 pos=coor coor=-180,-90 breaks=30		v.mkgrid map=grid_30 box=30,30 grid=6,12 pos=coor coor=-180,-90 breaks=30