Global datasets: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
(→‎ETOPO1: instructs)
Line 14: Line 14:


  # create a new simple XY location and mapset
  # create a new simple XY location and mapset
 
  # Import grid registered binary float
  # Import grid registered binary float
  r.in.bin -f in=etopo1_bed_g.flt out=etopo1_bed_g \
  r.in.bin -f in=etopo1_bed_g.flt out=etopo1_bed_g \
Line 20: Line 20:
     w=-180.00833333334 rows=10801 cols=21601 anull=-9999
     w=-180.00833333334 rows=10801 cols=21601 anull=-9999
  r.colors out=etopo1_bed_g color=etopo2
  r.colors out=etopo1_bed_g color=etopo2
 
  # reduce region by 1 cell
  # reduce region by 1 cell
  g.region rast=etopo1_bed_g
  g.region rast=etopo1_bed_g
  eval `g.region -g`
  eval `g.region -g`
  g.region n=n-$nsres s=s+$nsres e=e-$ewres -p
  g.region n=n-$nsres s=s+$nsres e=e-$ewres -p
 
  # save smaller raster and remove original
  # save smaller raster and remove original
  r.mapcalc "etopo1_bed_g.crop = etopo1_bed_g"
  r.mapcalc "etopo1_bed_g.crop = etopo1_bed_g"
  g.remove etopo1_bed_g
  g.remove etopo1_bed_g
 
  # change the location to lat/lon by restarting GRASS
  # change the location to lat/lon by restarting GRASS
  #  in the PERMANENT mapset and running {{cmd|g.setproj}}
  #  in the PERMANENT mapset and running {{cmd|g.setproj}}

Revision as of 07:30, 11 October 2008

Blue Marble imagery

ETOPO

ETOPO1

Cell registered version can be loaded directly into a lat/lon region.

Grid registered version must be loaded into a simple XY location and have the 90deg N,S rows cropped away.

# create a new simple XY location and mapset

# Import grid registered binary float
r.in.bin -f in=etopo1_bed_g.flt out=etopo1_bed_g \
   n=90.008333333335 s=-90.008333333335 e=180.00833333334 \
   w=-180.00833333334 rows=10801 cols=21601 anull=-9999
r.colors out=etopo1_bed_g color=etopo2

# reduce region by 1 cell
g.region rast=etopo1_bed_g
eval `g.region -g`
g.region n=n-$nsres s=s+$nsres e=e-$ewres -p

# save smaller raster and remove original
r.mapcalc "etopo1_bed_g.crop = etopo1_bed_g"
g.remove etopo1_bed_g

# change the location to lat/lon by restarting GRASS
#  in the PERMANENT mapset and running g.setproj

ETOPO2

ETOPO2v2 data download

Miscellaneous

LANDSAT

MODIS

Pathfinder AVHRR SST

merge that info here?

QuickBird

  • ?

SeaWiFS

merge that info here?

Smith and Sandwell

SRTM