Index: /issm/trunk-jpl/src/m/coordsystems/epsg2proj.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/epsg2proj.m	(revision 24777)
+++ /issm/trunk-jpl/src/m/coordsystems/epsg2proj.m	(revision 24777)
@@ -0,0 +1,17 @@
+function string = epsg2proj(epsg)
+%use gdalsrsinfo to provide Proj.4 compatible string from EPSG code
+%
+%   Usage:
+%      proj4string = epsg2proj(4326);
+%
+%   Example: 
+%      proj4string = epsg2proj(4326); 
+%      return proj4string='+proj=longlat +datum=wgs84 +no_defs '
+%
+
+	[s,r]=system(['gdalsrsinfo  epsg:' num2str(epsg)  ' | grep PROJ.4']);
+
+	if s~=0, 
+		error(r);
+	end
+	string=r(11:end-2);
Index: /issm/trunk-jpl/src/m/coordsystems/gdaltransform.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/gdaltransform.m	(revision 24776)
+++ /issm/trunk-jpl/src/m/coordsystems/gdaltransform.m	(revision 24777)
@@ -30,5 +30,5 @@
 	fclose(fid);
 
-	[s,r]=system(['gdaltransform -s_srs ',proj_in,' -t_srs ',proj_out,'  < ' filename_in ' > ' filename_out]);
+	[s,r]=system(['gdaltransform -s_srs "',proj_in,'" -t_srs "',proj_out,'"  < ' filename_in ' > ' filename_out]);
 	if s~=0 | ~isempty(deblank(r)),
 		error(r);
Index: /issm/trunk-jpl/src/m/coordsystems/laea.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/laea.m	(revision 24777)
+++ /issm/trunk-jpl/src/m/coordsystems/laea.m	(revision 24777)
@@ -0,0 +1,11 @@
+function string = laea(lat,long)
+%Lambert Azimuthal Equal Area projection at lat,long projection center. 
+%
+%   Usage:
+%      string = laea(45,-90);
+%
+%   Example: 
+%      string = laea(45,-90); 
+%      return string='+proj=laea +lat_0=45 +lon_0=-90 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs'
+
+	string=sprintf('+proj=laea +lat_0=%i +lon_0=%i +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs ',lat,long);