Marine Science: Difference between revisions
Jump to navigation
Jump to search
(→Bathymetric data: still trying to get proj right) |
|||
Line 47: | Line 47: | ||
To load it into GRASS lat/lon location (spherical): | To load it into GRASS lat/lon location (spherical): | ||
: Location setup: | : Location setup: | ||
: http://thread.gmane.org/gmane.comp.gis.gmt.user/918 | |||
: http://article.gmane.org/gmane.comp.gis.proj-4.devel/192/ | |||
# ( | Process with GMT's img2grd to convert from spherical Mercator projection to geographic coordinates, then import into GRASS | ||
img2grd topo_9.1b.img -T1 -S1 -V -R0/360/-80.738/80.738 -m1 -D -Gtopo_all.grd | |||
# (out of memory, needs 1.4gb) | |||
# try just for NZ (W/E/S/N bounds) | |||
REGION=160/180/-50/-30 | |||
img2grd topo_9.1b.img -T1 -S1 -V -R"$REGION" -m1 -D -Gtopo_NZ.grd | |||
grd2xyz topo_NZ.grd -S > topo_NZ.xyz | |||
# check that rows and columns match (some NS res problem there) | |||
g.region n=30S s=50S w=160E e=180 ewres=0:01 nsres=0.01260 -p | |||
r.in.xyz in=topo_NZ_180.xyz out=topo_NZ_1min fs=tab | |||
Is it possible to load directly into GRASS? | |||
How to set up Merc/Sphere proj? | |||
# the following does not work correctly, just a trial | |||
# offset n,s,e,w by 1/2 a grid cell? | # offset n,s,e,w by 1/2 a grid cell? | ||
r.in.bin input=topo_9.1b.img output=topo_9.1b \ | r.in.bin input=topo_9.1b.img output=topo_9.1b \ |
Revision as of 12:47, 24 April 2008
this page is a work in progress
Tools for marine scientists
Bathymetry processing
- Basic DEM creation: v.in.ascii + v.surf.rst
- r.in.xyz for processing multibeam sonar swaths
- r.surf.nnbathy add-on script for creating bathymetic DEMs from input x,y,z data
- Multibeam sonar swath transect planning script
Bathymetric data
- See the ETOPO2 article by M.H. Bowman in the GRASS Newsletter, 1:8-11, August 2004.
- Smith and Sandwell 1-minute global elevation, v9.1, August 21, 2007
global_topo_1min/README_V9.1.txt file:
Version 9.1 has a very different FORMAT than V8.2 The main differences are that the grid spacing in longitude is now 1 minute rather than 2 minutes. In addition, the latitude range is increased to +/- 80.738. Like the old versions, the elevation(+) and depth(-) are stored as 2-byte integers to the nearest meter. An odd depth of say -2001m signifies that this pixel was constrained by a real depth sounding while an even depth of say -2000m is a predicted depth. Here are the parameters for the old and new versions: param V8.2 V9.2 ___________________________ nlon 10800 21600 nlat 12672 17280 rlt0 -72.006 -80.738 rltf 72.006 80.738 ___________________________ The binary format of the integers is bigendian so the bytes need to be swapped if you are running on an Intel processor. Here is a typical command for swapping bytes: dd if=topo_9.1.img of=topo_9.1.img.swab bs=21600 conv=swab.
To load it into GRASS lat/lon location (spherical):
- Location setup:
- http://thread.gmane.org/gmane.comp.gis.gmt.user/918
- http://article.gmane.org/gmane.comp.gis.proj-4.devel/192/
Process with GMT's img2grd to convert from spherical Mercator projection to geographic coordinates, then import into GRASS
img2grd topo_9.1b.img -T1 -S1 -V -R0/360/-80.738/80.738 -m1 -D -Gtopo_all.grd # (out of memory, needs 1.4gb) # try just for NZ (W/E/S/N bounds) REGION=160/180/-50/-30 img2grd topo_9.1b.img -T1 -S1 -V -R"$REGION" -m1 -D -Gtopo_NZ.grd grd2xyz topo_NZ.grd -S > topo_NZ.xyz # check that rows and columns match (some NS res problem there) g.region n=30S s=50S w=160E e=180 ewres=0:01 nsres=0.01260 -p r.in.xyz in=topo_NZ_180.xyz out=topo_NZ_1min fs=tab
Is it possible to load directly into GRASS?
How to set up Merc/Sphere proj?
# the following does not work correctly, just a trial # offset n,s,e,w by 1/2 a grid cell? r.in.bin input=topo_9.1b.img output=topo_9.1b \ title="1' worldwide relief (1.852 km-sq)" \ -b -s bytes=2 rows=17280 cols=21600 \ n=80.738 s=-80.738 w=0 e=360 # offset n,s,e,w by 1/2 a grid cell?
r.colors output=topo_9.1b color=etopo2
Sidescan sonar processing
- i.gdalwarp script for georectifying and mosaicking scanned paper rolls into a GeoTIFF
- v.swathwidth script for planning swath bathymetry surveys
Wave exposure
- Using GRASS to prepare and process data for the SWAN Wave Model
- preparing input DEM
- r.in.mat and r.out.mat
Circulation models
- Preparing input grids
- r.in.mat and r.out.mat
Tutorials
Remote Sensing
- Importing MODIS Aqua SST and Chlorophyll-a data