Changeset 17004

Timestamp:

12/04/13 09:11:36 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: moved print/export_geotiff.m

File:

• issm/trunk-jpl/src/m/print/export_geotiff.m (moved) (moved from issm/trunk-jpl/src/m/interp/export_geotif.m ) (1 diff)

issm/trunk-jpl/src/m/print/export_geotiff.m

@@ -1 @@
-function export_geotif(filename,ref)
+function export_geotiff(filename,ref)
 %EXPORT_GEOTIF - export geotiff 
 %
+%   Usage:
+%      export_geotif(filename,ref);
 %
-% Usage: export_geotif('Articles/figures/figure','EPSG:3411');
+%      This function must be called after plotmodel
+%      filname: no extension
+%      ref:     UPS Greenland  EPSG:3413 (http://www.spatialreference.org/ref/epsg/3413/)
+%               UPS Antarctica EPSG:3031 (http://www.spatialreference.org/ref/epsg/3031/)
 %
+%   Example:
+%      export_geotif('Greenland','EPSG:3413');
 
+%Get axis limits and convert to strings
+XLIM = xlim(); x0 = num2str(XLIM(1)); x1 = num2str(XLIM(2));
+YLIM = ylim(); y0 = num2str(YLIM(1)); y1 = num2str(YLIM(2));
 
-        %first export the figure
-        export_fig([filename '.jpg']);
+%first export the figure
+export_fig([filename '.jpg']);
 
-        %call gdal on this: 
-        system(sprintf('gdal_translate -a_srs %s -of GTiff -co "INTERLEAVE=PIXEL" -a_ullr %g %g %g %g %s.jpg %s.tif', ref, min(xlim), max(ylim),max(xlim),min(ylim),filename,filename));
-        system(['rm -rf ' filename '.jpg']);
-
-end
+%call gdal on this: 
+system(['gdal_translate -a_srs ' ref ' -of GTiff -co "INTERLEAVE=PIXEL" -a_ullr ' x0 ' ' y1 ' ' x1 ' ' y0 ' ' filename '.jpg ' filename '.tif']);
+system(['rm -rf ' filename '.jpg']);