Line | |
---|
1 | function string = epsg2proj(epsg)
|
---|
2 | %EPSG2PROJ - uses gdalsrsinfo to provide PROJ.4 compatible string from
|
---|
3 | %EPSG code
|
---|
4 | %
|
---|
5 | % Usage:
|
---|
6 | % proj4string = epsg2proj(4326);
|
---|
7 | %
|
---|
8 | % Example:
|
---|
9 | % proj4string = epsg2proj(4326);
|
---|
10 | % return proj4string='+proj=longlat +datum=wgs84 +no_defs'
|
---|
11 | %
|
---|
12 |
|
---|
13 | %Call PROJ library
|
---|
14 | [status, string]=system(['projinfo -o PROJ -q epsg:' num2str(epsg)]);
|
---|
15 |
|
---|
16 | %Check status
|
---|
17 | if status~=0; error(string); end
|
---|
18 |
|
---|
19 | %remove trailing blanks
|
---|
20 | string = deblank(string);
|
---|
Note:
See
TracBrowser
for help on using the repository browser.