<!doctype html public "-//W3C//DTD W3 HTML//EN">
<html><head><style type="text/css"><!--
blockquote, dl, ul, ol, li { padding-top: 0 ; padding-bottom: 0 }
--></style><title>ASTER HDF L1_B problem in Mac OS
X</title></head><body>
<div>Hi</div>
<div><br></div>
<div>When I use gdalinfo for viewing subdataset in HDF file (utm zone
33 with coordinate in lat-lon I have an error to the beginning:</div>
<div><br></div>
<blockquote><font face="Courier">[lorenzo:/usr/local/bin] ben%
./gdalinfo
HDF4_SDS:ASTER_L1B:"/Volumes/ESTERNO/ASTER_FO_gen02/AST_L1B_003_<span
></span>01102002101254_05222002150244.hdf":2</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">ERROR 6: Unable to load PROJ.4
library (libproj.so), creation of</font></blockquote>
<blockquote><font face="Courier">OGRCoordinateTransformation
failed.</font></blockquote>
<blockquote><font face="Courier">ERROR 6: Unable to load PROJ.4
library (libproj.so), creation of</font></blockquote>
<blockquote><font face="Courier">OGRCoordinateTransformation
failed.</font></blockquote>
<blockquote><font face="Courier">ERROR 6: Unable to load PROJ.4
library (libproj.so), creation of</font></blockquote>
<blockquote><font face="Courier">OGRCoordinateTransformation
failed.</font></blockquote>
<blockquote><font face="Courier">ERROR 6: Unable to load PROJ.4
library (libproj.so), creation of</font></blockquote>
<blockquote><font face="Courier">OGRCoordinateTransformation
failed.</font></blockquote>
<blockquote><font face="Courier">Driver: HDF4Image/HDF4
Dataset</font></blockquote>
<blockquote><font face="Courier">Size is 4980,
4200</font></blockquote>
<blockquote><font face="Courier">Coordinate System
is:</font></blockquote>
<blockquote><font face="Courier">PROJCS["UTM Zone 33, Northern
Hemisphere",</font></blockquote>
<blockquote><font face="Courier"> GEOGCS["WGS
84",</font></blockquote>
<blockquote><font face="Courier">
DATUM["WGS_1984",</font></blockquote>
<blockquote><font
face="Courier"
>
SPHEROID["WGS 84",6378137,298.257223563,</font></blockquote>
<blockquote><font
face="Courier"
> <span
></span>
AUTHORITY["EPSG","7030"]],</font></blockquote>
<blockquote><font
face="Courier"
>
TOWGS84[0,0,0,0,0,0,0],</font></blockquote>
<blockquote><font
face="Courier"
>
AUTHORITY["EPSG","6326"]],</font></blockquote>
<blockquote><font face="Courier">
PRIMEM["Greenwich",0,</font></blockquote>
<blockquote><font
face="Courier"
>
AUTHORITY["EPSG","8901"]],</font></blockquote>
<blockquote><font face="Courier">
UNIT["degree",0.0174532925199433,</font></blockquote>
<blockquote><font
face="Courier"
>
AUTHORITY["EPSG","9108"]],</font></blockquote>
<blockquote><font face="Courier">
AXIS["Lat",NORTH],</font></blockquote>
<blockquote><font face="Courier">
AXIS["Long",EAST],</font></blockquote>
<blockquote><font face="Courier">
AUTHORITY["EPSG","4326"]],</font></blockquote>
<blockquote><font face="Courier">
PROJECTION["Transverse_Mercator"],</font></blockquote>
<blockquote><font face="Courier">
PARAMETER["latitude_of_origin",0],</font></blockquote>
<blockquote><font face="Courier">
PARAMETER["central_meridian",15],</font></blockquote>
<blockquote><font face="Courier">
PARAMETER["scale_factor",0.9996],</font></blockquote>
<blockquote><font face="Courier">
PARAMETER["false_easting",500000],</font></blockquote>
<blockquote><font face="Courier">
PARAMETER["false_northing",0],</font></blockquote>
<blockquote><font face="Courier">
UNIT["Meter",1]]</font></blockquote>
<blockquote><font face="Courier">GeoTransform =</font></blockquote>
<blockquote><font face="Courier"> 12.19647665061962,
0.0001816592369477914, -3.696047619047622e-05</font></blockquote>
<blockquote><font face="Courier"> 43.49892337703098,
-2.646787148594408e-05, -0.0001322861904761889</font></blockquote>
<blockquote><font face="Courier">GCP Projection = GEOGCS["WGS
84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030<span
></span
>"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","632<span
></span
>6"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG",<span
></span
>"8901"]],UNIT["degree",0.0174532925199433,AUTHOR<span
></span
>ITY["EPSG","9108"]],AXIS["Lat",NORTH],<span
></span
>AXIS["Long",EAST],AUTHORITY["EPSG","4326<span
></span>"]]</font></blockquote>
<blockquote><font face="Courier">GCP[ 0]: Id=,
Info=</font></blockquote>
<blockquote><font
face="Courier">
(0.5,0.5) -> (12.1965,43.4988,0)</font></blockquote>
<blockquote><font face="Courier">GCP[ 1]: Id=,
Info=</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">...</font></blockquote>
<blockquote><font face="Courier">snip</font></blockquote>
<blockquote><font face="Courier">...</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">
SOUTHBOUNDINGCOORDINATE=42.811432<br>
TIMEOFDAY=101254443000Z<br>
CALENDARDATE=20020110<br>
FUTUREREVIEWDATE=20020312<br>
SCIENCEREVIEWDATE=20011201<br>
QAPERCENTMISSINGDATA=0.000006<br>
QAPERCENTOUTOFBOUNDSDATA=0.000006<br>
QAPERCENTINTERPOLATEDDATA=0.000000<br>
REPROCESSINGACTUAL=not reprocessed<br>
PGEVERSION=04.00R06<br>
PROCESSINGLEVELID=1B<br>
MAPPROJECTIONNAME=Universal Transverse Mercator<br>
Corner Coordinates:<br>
Upper Left ( 12.1964767, 43.4989234)<br>
Lower Left ( 12.0412427, 42.9433214)<br>
Upper Right ( 13.1011397, 43.3671134)<br>
Lower Right ( 12.9459057, 42.8115114)<br>
Center ( 12.5711912,
43.1552174)</font></blockquote>
<blockquote><font face="Courier">Band 1 Block=4980x1 Type=Byte,
ColorInterp=Gray</font></blockquote>
<blockquote><br></blockquote>
<div>When I use gdal_translate I have the same error.</div>
<div>When I use Grass with r.in.gdal module I have the same
error.</div>
<div><br></div>
<div>I see the code and in the file</div>
<div>../ogr/ogrct.cpp</div>
<div><br></div>
<div>there is</div>
<div><br></div>
<blockquote><font face="Courier">CPL_CVSID("$Id: ogrct.cpp,v 1.23
2004/01/24 09:35:00 warmerda Exp $");</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">/*
====================================================================
*/</font></blockquote>
<blockquote><font face="Courier">/*
PROJ.4 interface
stuff. <span
></span
> <span
></span
> <span
></span>
*/</font></blockquote>
<blockquote><font face="Courier">/*
====================================================================
*/</font></blockquote>
<blockquote><font face="Courier">#ifndef
PROJ_STATIC</font></blockquote>
<blockquote><font face="Courier">typedef struct { double u, v; }
projUV;</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">#define projPJ void
*</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">#define
RAD_TO_DEG
57.29577951308232</font></blockquote>
<blockquote><font face="Courier">#define
DEG_TO_RAD
.0174532925199432958</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">#endif</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">static
projPJ (*pfn_pj_init_plus)(const
char *) = NULL;</font></blockquote>
<blockquote><font face="Courier">static
projPJ (*pfn_pj_init)(int, char**)
= NULL;</font></blockquote>
<blockquote><font face="Courier">static
projUV (*pfn_pj_fwd)(projUV,
projPJ) = NULL;</font></blockquote>
<blockquote><font face="Courier">static
projUV (*pfn_pj_inv)(projUV,
projPJ) = NULL;</font></blockquote>
<blockquote><font face="Courier">static void
(*pfn_pj_free)(projPJ) = NULL;</font></blockquote>
<blockquote><font face="Courier">static int
(*pfn_pj_transform)(projPJ, projPJ, long, int,</font></blockquote>
<blockquote><font
face="Courier"
> <span
></span
> <span
></span
> <span
></span> double *, double *, double * ) =
NULL;</font></blockquote>
<blockquote><font face="Courier">static
int
*(*pfn_pj_get_errno_ref)(void) = NULL;</font></blockquote>
<blockquote><font face="Courier">static
char
*(*pfn_pj_strerrno)(int) = NULL;</font></blockquote>
<blockquote><font face="Courier">static
char
*(*pfn_pj_get_def)(projPJ,int) = NULL;</font></blockquote>
<blockquote><font face="Courier">static
void
(*pfn_pj_dalloc)(void *) = NULL;</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">#ifdef WIN32</font></blockquote>
<blockquote><font face="Courier"># define
LIBNAME
"proj.dll"</font></blockquote>
<blockquote><font face="Courier">#else</font></blockquote>
<blockquote><font face="Courier"># define
LIBNAME
"libproj.so"</font></blockquote>
<blockquote><font face="Courier">#endif</font></blockquote>
<blockquote><font face="Courier">....</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">snip</font></blockquote>
<blockquote><font face="Courier"><br></font></blockquote>
<blockquote><font face="Courier">....</font></blockquote>
<blockquote><font
face="Courier"
>/*******************************************************************<span
></span>*****/<br>
/* <span
></span>
OGRCreateCoordinateTransformation() <span
></span
>
*/<br>
/********************************************************************<span
></span>****/<br>
<br>
/**<br>
* Create transformation object.<br>
*<br>
* This is the same as the C function
OCTNewCoordinateTransformation().<br>
*<br>
* The delete operator, or OCTDestroyCoordinateTransformation()
should<br>
* be used to destroy transformation objects.<br>
*<br>
* @param poSource source spatial reference system.<br>
* @param poTarget target spatial reference system.<br>
* @return NULL on failure or a ready to use transformation
object.<br>
*/<br>
<br>
OGRCoordinateTransformation*<br>
OGRCreateCoordinateTransformation( OGRSpatialReference *poSource,<br>
<span
></span
> <span
></span
>
OGRSpatialReference *poTarget )<br>
<br>
{<br>
OGRProj4CT *poCT;<br>
<br>
if( !LoadProjLibrary() )<br>
{<br>
CPLError( CE_Failure,
CPLE_NotSupported,<br>
<span
></span> "Unable to load
PROJ.4 library (%s), creation of\n"<br>
<span
></span>
"OGRCoordinateTransformation failed.",<br>
<span
></span> LIBNAME );<br>
return NULL;<br>
}<br>
<br>
poCT = new OGRProj4CT();<br>
<br>
if( !poCT->Initialize( poSource, poTarget )
)<br>
{<br>
delete poCT;<br>
return NULL;<br>
}<br>
else<br>
{<br>
return poCT;<br>
}</font></blockquote>
<blockquote><font face="Courier">}<br>
</font></blockquote>
<blockquote>.....</blockquote>
<blockquote><br></blockquote>
<blockquote><br></blockquote>
<blockquote><br></blockquote>
<div>In gdal config there is'n --with-proj=DIR and the error is from
this file.</div>
<div>I have PROJ4.4.7 installed as shared library (libproj.dylib in
Mac OS X in /usr/local/lib)</div>
<div>I have changed the name in source (libproj.so changed
libproj.dylib) but the error remains ...</div>
<div><br></div>
<div>Is this file important for coordinate recognition ?</div>
<div><br></div>
<div>Can you help me ?</div>
<div> </div>
<div><br></div>
<div>Thanx</div>
<x-sigsep><pre>--
</pre></x-sigsep>
<div><font face="Verdana"
color="#000000"
>____________________________________________________________________<span
></span>____</font></div>
<div><font face="Courier New"
color="#000000">|| Lorenzo
Moretti e-mail:
lorenzo.moretti@bologna.enea.it </font></div
>
<div><font face="Courier New" color="#000000">||/|/| ENEA
prot <span
></span> Web:
http://wwwamb.bologna.enea.it/ </font></div>
<div><font face="Courier New" color="#000000">|| |
via Don Fiammelli, 2 FTP:
ftp://ftpamb.bologna.enea.it/ (ris.)</font></div>
<div><font face="Courier New" color="#000000">~~~~~~ 40128
BOLOGNA - ITALY Ph: +39-0516098086 Fax:
+39-0516098131</font></div>
<div><font face="Courier New"
color="#000000"
>____________________________________________________________________<span
></span>____</font></div>
</body>
</html>