Global datasets: Difference between revisions

From GRASS-Wiki
Jump to navigation Jump to search
(+ import from WCS)
Line 165: Line 165:
* {{cmd|r.in.gdal}} will read .BIL format, see the {{cmd|r.in.srtm}} shell script for hints
* {{cmd|r.in.gdal}} will read .BIL format, see the {{cmd|r.in.srtm}} shell script for hints
* binary format is 2 byte integer. Multiply by 10 using {{cmd|r.mapcalc}} to convert units. See http://www.worldclim.org/format.htm for more information and the [[MODIS]] help page for example of converting raw to data units.
* binary format is 2 byte integer. Multiply by 10 using {{cmd|r.mapcalc}} to convert units. See http://www.worldclim.org/format.htm for more information and the [[MODIS]] help page for example of converting raw to data units.
=== OGC WCS - Albedo example ===
GRASS imports OGC Web Coverage Service data. Example server (please suggest a better one!)
  <WCS_GDAL>
  <ServiceURL>http://laits.gmu.edu/cgi-bin/NWGISS/NWGISS?</ServiceURL>
  <CoverageName>AUTUMN.hdf</CoverageName>
  <Timeout>90</Timeout>
  <Resample>nearest</Resample>
  </WCS_GDAL>
Save this as albedo.xml. Import into a LatLong WGS84 location:
  r.in.gdal albedo.xml out=albedo
Unfortunately this server sends out the map shifted by 0.5 pixel. This requires a fix to the map boundary coordinates:
  r.region albedo n=90 s=-90 w=-180 e=180
Now apply color table and look at the map:
  r.colors albedo color=byr
  d.mon x0
  d.rast albedo


== Vector data ==
== Vector data ==

Revision as of 16:03, 16 May 2009

Raster data

AVHRR

Blue Marble imagery

NASA's Blue Marble is a 500m-8 degree per-cell world wide visual image of the Earth from space, with the clouds removed.

Natural Earth imagery

True Marble imagery

  • True Marble: 250m world wide visual image of the Earth from space, with the clouds removed. GeoTIFF (use the r.in.gdal module)

ETOPO

The ETOPO datasets provide global topography and bathymetry at 1', 2', and 5' per-cell resolutions.

ETOPO1

The cell registered version can be loaded directly into a lat/lon region. GRASS raster data is cell registered (see the GRASS raster semantics page)


The grid registered version must be loaded into a simple XY location using the parameters found in the .hdr file, and needs to 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
#  or move the entire mapset into a lat/lon location
#  and manually edit the $MAPSET/cellhd/ files (dirty!)

ETOPO2

ETOPO2v2 data download

CleanTOPO2

GEBCO

  • The General Bathymetric Chart of the Oceans (original 1' release 2003, new 1' and 30" releases 2008)
http://www.gebco.net/data_and_products/gridded_bathymetry_data/
http://www.bodc.ac.uk/data/online_delivery/gebco/

r.in.gdal can be used to import the GMT netCDF files directly, or if that doesn't work you can use GMT tools to convert to an old-style native GMT format and import that with r.in.bin.

example: (GEBCO 2003 1' data)
# convert to an old style GMT binary .grd using grdreformat
$ grdreformat 3n24s47w14w.grd 3n24s47w14w_Native.grd=bs

# then import into GRASS,
GRASS> r.in.bin -h -s bytes=2 in=3n24s47w14w_Native.grd out=3n24s47w14w

# and set some nice colors
GRASS> r.colors 3n24s47w14w rules=- << EOF
nv magenta
0% black
-7740 0:0:168
0 84:176:248
0 40:124:0
522 68:148:24
1407 148:228:108
1929 232:228:108
2028 232:228:92
2550 228:160:32
2724 216:116:8
2730 grey
2754 grey
2760 252:252:252
2874 252:252:252
2883 192:192:192
2913 192:192:192
100% 252:252:252
EOF

Miscellaneous

Data sources


Modules

  • The r.in.gdal modules may be used to import data of many formats, including GMT netCDF
  • The r.in.bin module may be used to import raw binary files

LANDSAT

Since October 1, 2008 all Landsat 7 ETM+ scenes held in the USGS EROS archive are available for download at no charge.

  • Download via the Glovis online search tool (req. Java)
  • Download via the USGS's EarthExplorer interface

Modules

  • r.in.gdal - Main import tool for complete multiband scenes
  • r.in.wms - Download data covering current map region via WMS server
  • r.in.onearth - WMS frontend for NASA's OnEarth Global Landsat Mosaic
  • i.landsat.rgb - Color balancing/enhancement tool

See also

  • Processing tips can be found on the LANDSAT wiki page

MODIS

Pathfinder AVHRR SST

merge that info here?

QuickBird

SeaWiFS

Smith and Sandwell

SRTM

Space Shuttle Radar Topography Mission

WorldClim

WorldClim is a set of global climate layers (climate grids) with a spatial resolution of a square kilometer.


OGC WCS - Albedo example

GRASS imports OGC Web Coverage Service data. Example server (please suggest a better one!)

 <WCS_GDAL>
 <ServiceURL>http://laits.gmu.edu/cgi-bin/NWGISS/NWGISS?</ServiceURL>
 <CoverageName>AUTUMN.hdf</CoverageName>
 <Timeout>90</Timeout>
 <Resample>nearest</Resample>
 </WCS_GDAL>

Save this as albedo.xml. Import into a LatLong WGS84 location:

 r.in.gdal albedo.xml out=albedo

Unfortunately this server sends out the map shifted by 0.5 pixel. This requires a fix to the map boundary coordinates:

 r.region albedo n=90 s=-90 w=-180 e=180

Now apply color table and look at the map:

 r.colors albedo color=byr
 d.mon x0
 d.rast albedo

Vector data

CDC Geographic Boundary and Public Health Maps

Global Administrative Areas

  • GADM is a database of the location of the world's administrative areas (boundaries) available in shapefiles.
http://biogeo.berkeley.edu/gadm/ (extracted by country here)
  • World Borders Dataset including ISO 3166-1 Country codes available in shapefiles.
http://thematicmapping.org/downloads/world_borders.php
  • Free GIS data from Mapping Hacks
http://mappinghacks.com/data/

GSHHS World Coastline

GSHHS is a high resolution shoreline dataset. It is derived from data in the public domain and licensed as GPL. The shorelines are constructed entirely from hierarchically arranged closed polygons. It is closely linked to the GMT project.

For GRASS 6 you can download 1:250,000 shoreline data from NOAA's site in Mapgen format, which can be imported with the v.in.mapgen module.

SALB

Second Administrative Level Boundaries: "The SALB dataset is a global digital dataset consisting of digital maps and codes that can be downloaded on a country by country basis."

VMap0

1:1 million vector data. Formerly known as Digital Chart of the World

Misc