extract bathymetry

Facts, news, and guidance about ROMS software

Moderators: arango, robertson

Post Reply
Message
Author
mistral
Posts: 23
Joined: Wed Jan 16, 2008 11:17 pm
Location: Privat

extract bathymetry

#1 Post by mistral » Sat Nov 07, 2009 1:33 pm

I can not extract a bathymetry.
I decided to follow this path: https: / / www.myroms.org / wiki / index.php / easygrid

following this path I extract the coastline: https: / / www.myroms.org / wiki / index.php / seagrid

the coastline is ok!


I have to extract the bathymetry
or used 3 methods
1) read_srtm30plus.m
2) Smith, Sandwell [http://topex.ucsd.edu/WWW_html/mar_topo.html] [ftp://topex.ucsd.edu/pub/global_topo_1min/] with extract.me topo_12.1.img
3) etopo5 and matlab script

in all three methods have problems

1)>> [lon, lat, z] = read_srtm30plus([10 20], [40 46], 30.1);

url =

http://onearth.jpl.nasa.gov/wms.cgi?req ... us&bbox=10, 40,20,46 & SRS = EPSG: 4326 & TRANSPARENT = FALSE & width = 1200 & height = 720 & format = image / GeoTIFF & version = 1.1.1

?? Error using ==> imread at 293
Can not read URL
"http://onearth.jpl.nasa.gov/wms.cgi?req ... us&bbox=10, 40,20,46 & SRS = EPSG: 4326 & TRANSPARENT = FALSE & width = 1200 & height = 720 & format = image / GeoTIFF & version = 1.1.1 .

Error in ==> read_srtm30plus at 60
pix = imread (url, 'tif');


2)>> [image_data, vlat, vlon] = mygrid_sand ([38 42 12 17], 1)

iblat =

1934


itlat =

1778


iwlon =

360


ielon =

510

Reading in bathymetry data
?? Error using ==> fseek
Invalid file identifier. Use fopen to generate a valid file identifier.

Error in ==> mygrid_sand at 91
status = fseek (fid, offset, 'bof');


3) [lon, lat, h] = x_etopo5 ([12,17,38,42], etopo5);
Warning: Function ETOPO5 is obsolete and will be removed in the next release.
Use eTopo instead.
> In etopo5 at 62
?? Error using ==> eTopo> getFileFromDir at 252
Directory "/ home / joe / Bati / matlab / bathymetry" does not contain eTopo date.

Error in ==> eTopo> parseInputs at 95
file = getFileFromDir (datadir);

Error in ==> eTopo at 67
[file, samplefactor, latlim, lonlim, readFcn] = parseInputs (varargin {:});

Error in ==> etopo5 at 66
[Z, refvec] = eTopo (varargin {:});

>




Use MATLAB 7.6.0 R2008a
I hope in your help

slvester
Posts: 32
Joined: Mon Jun 01, 2009 12:59 pm
Location: second institution of oceanography,state oceanic administration China

Re: extract bathymetry

#2 Post by slvester » Wed Nov 25, 2009 1:19 am

[quote][1)>> [lon, lat, z] = read_srtm30plus([10 20], [40 46], 30.1);
Error in ==> read_srtm30plus at 60
pix = imread (url, 'tif');
/quote]
hello mistral,I use the first method to extract the bathymetry ,and I have the proble the same with you.
Have you resolved it?
help!

Post Reply