Can't read URL when using "read_srtm30plus.m" function

[kee]

Hi all,
I tried to get the bathy data via read_srtm30plus.m, and it returned wrong mesages like below.
can anyone tell me why?
and how to fix it?
thanks a lot.
I feel bad these days since always having problems...

Code: Select all

>> [lon,lat,z]=read_srtm30plus([-5 30],[95 135],60,1);

url =

http://onearth.jpl.nasa.gov/wms.cgi?request=GetMap&styles=short_int&layers=srtmplus&bbox=-5,95,30,135&SRS=EPSG:4326&TRANSPARENT=FALSE&width=2100&height=2400&format=image/geotiff&version=1.1.1

??? Error using ==> imread at 353
Can't read URL
"http://onearth.jpl.nasa.gov/wms.cgi?request=GetMap&styles=short_int&layers=srtmplus&bbox=-5,95,30,135&SRS=EPSG:4326&TRANSPARENT=FALSE&width=2100&height=2400&format=image/geotiff&version=1.1.1".

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

[wilkin]

If you go to the top level of that URL http://onearth.jpl.nasa.gov/ you will see

STATUS:

The OnEarth server for MODIS is currently being phased out and replaced by the Global Imagery Browse Services (GIBS) provided by the Earth Science Data and Information System (ESDIS) Project

There are other places to get SRTM bathy data, starting from:
http://topex.ucsd.edu/WWW_html/srtm30_plus.html
for the most recent version. The ftp site gives the data in regional subsets in netcdf format.

Older versions of SRTM30+ are available via THREDDS at http://oceanwatch.pfeg.noaa.gov/thredds/catalog.html

[kee]

Hi Wilkin
Thanks a lot!

[pilipalapipa]

Dear friends,

I tried the method below, but I got the wrong message like this,
-----------------------------------------
>> [x,y,z]=read_srtm30plus(ax(1:2),ax(3:4),30);
Warning: The netcdf-java cdm contains no coordinate information
associated with the variable. Returning ncvariable instead of
ncgeovariable object. (Methods that rely on coordinate information
like 'grid' or 'geosubset' are not available.
> In ncgeodataset>ncgeodataset.geovariable at 215
In ncgeodataset>ncgeodataset.subsref at 378
In read_srtm30plus2 at 32
??? Error using ==> subsref
No appropriate method, property, or field geosubset for class
ncvariable.

Error in ==> ncvariable>ncvariable.subsref at 295
sref = builtin('subsref',obj,s);

Error in ==> read_srtm30plus2 at 36
g=topovar.geosubset(s);

------------------------------
I will be pretty appreciated if anyone gives some suggestions,

Tony

Hi all,
I tried to get the bathy data via read_srtm30plus.m, and it returned wrong mesages like below.
can anyone tell me why?
and how to fix it?
thanks a lot.
I feel bad these days since always having problems...

Code: Select all

>> [lon,lat,z]=read_srtm30plus([-5 30],[95 135],60,1);

url =

http://onearth.jpl.nasa.gov/wms.cgi?request=GetMap&styles=short_int&layers=srtmplus&bbox=-5,95,30,135&SRS=EPSG:4326&TRANSPARENT=FALSE&width=2100&height=2400&format=image/geotiff&version=1.1.1

??? Error using ==> imread at 353
Can't read URL
"http://onearth.jpl.nasa.gov/wms.cgi?request=GetMap&styles=short_int&layers=srtmplus&bbox=-5,95,30,135&SRS=EPSG:4326&TRANSPARENT=FALSE&width=2100&height=2400&format=image/geotiff&version=1.1.1".

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