Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Wed Jul 18, 2018 6:40 pm




Post new topic Reply to topic  [ 2 posts ] 

All times are UTC

Author Message
 Post subject: extract bathymetry
PostPosted: Sat Nov 07, 2009 1:33 pm 
Offline

Joined: Wed Jan 16, 2008 11:17 pm
Posts: 23
Location: Privat
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?request=GetMap&styles=short_int&layers=srtmplus&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


Top
 Profile  
Reply with quote  
 Post subject: Re: extract bathymetry
PostPosted: Wed Nov 25, 2009 1:19 am 
Offline

Joined: Mon Jun 01, 2009 12:59 pm
Posts: 32
Location: second institution of oceanography,state oceanic administration China
[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!


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group