<HTML dir=ltr><HEAD><TITLE>Re: R: [GRASS-dev] WinGRASS: r.in.gdal projection issue</TITLE>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.6000.16640" name=GENERATOR><BASE href=http://mailstore.rossoalice.alice.it/exchange/Alice000000001264700/Bozze/R:%20R:%20[GRASS-dev]%20WinGRASS:%20r.in.gdal%20projection%20issue.EML/1_text.htm></HEAD>
<BODY>
<DIV id=idOWAReplyText54394 dir=ltr>
<DIV dir=ltr><FONT face=Arial color=#000000 size=2>Hi Paul,</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>>How big is the file? Would it be possible to put it somewhere for me to<BR>>have a look at and try and reproduce the behaviour?</FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr><FONT size=2>the file is available here:</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><A href="http://armadillo.geog.kcl.ac.uk/portal/srtm3/srtm_data_arcascii/srtm_39_03.zip"><FONT size=2>http://armadillo.geog.kcl.ac.uk/portal/srtm3/srtm_data_arcascii/srtm_39_03.zip</FONT></A></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>(I selected the server nearest to you)</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>>gdalinfo filename</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>gdalinfo {C:/Documents and Settings/Marco/Documenti/GIS Data Repository/SRTM v3.0 DEM 90m/Z_39_3.ASC}</FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr><FONT size=2>output:</FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>Driver: AAIGrid/Arc/Info ASCII Grid<BR>Files: C:/Documents and Settings/Marco/Documenti/GIS Data Repository/SRTM v3.0 DEM 90m/Z_39_3.ASC<BR>Size is 6000, 6000<BR>Coordinate System is `'<BR>Origin = (10.000000000000000,49.999999999999979)<BR>Pixel Size = (0.000833333333333,-0.000833333333333)<BR>Corner Coordinates:<BR>Upper Left ( 10.0000000, 50.0000000) <BR>Lower Left ( 10.0000000, 45.0000000) <BR>Upper Right ( 15.0000000, 50.0000000) <BR>Lower Right ( 15.0000000, 45.0000000) <BR>Center ( 12.5000000, 47.5000000) <BR>Band 1 Block=6000x1 Type=Int16, ColorInterp=Undefined<BR> NoData Value=-9999</FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr><FONT size=2>>g.proj -p georef=filename</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>g.proj -p {georef=C:/Documents and Settings/Marco/Documenti/GIS Data Repository/SRTM v3.0 DEM 90m/Z_39_3.ASC}</FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr><FONT size=2>output:</FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>XY location (unprojected)</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>Trying to open with OGR...</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>Trying to open with GDAL...</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>...succeeded.</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2>Read of file C:/Documents and Settings/Marco/Documenti/GIS Data Repository/SRTM v3.0 DEM 90m/Z_39_3.ASC was successful, but it did not contain projection information. 'XY (unprojected)' will be used</FONT></DIV>
<DIV dir=ltr><FONT face="Courier New" size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>It seems that it cannot read the coordinate system (or that it is not present... but I beleive in the first option).</FONT></DIV>
<DIV dir=ltr><FONT size=2>Probably we should need to search the problem in the GDAL build?</FONT></DIV>
<DIV dir=ltr><FONT size=2>Here the current GDAL (1.5.1) configuration:</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><PRE><CODE>GDAL is now configured for i686-pc-mingw32
Installation directory: /usr/local
C compiler: gcc -g -O2
C++ compiler: g++ -g -O2
LIBTOOL support: yes
LIBZ support: internal
GRASS support: no
CFITSIO support: no
PCRaster support: internal
NetCDF support: no
LIBPNG support: internal
LIBTIFF support: internal (BigTIFF=yes)
LIBGEOTIFF support: internal
LIBJPEG support: internal
LIBGIF support: internal
OGDI support: no
HDF4 support: no
HDF5 support: no
Kakadu support: no
JasPer support: no
ECW support: no
MrSID support: no
GRIB support: no
cURL support (wms/wcs/...): no
PostgreSQL support: yes
MySQL support: no
Xerces-C support: no
Expat support: yes
ODBC support: no
PGeo support: no
OCI support: no
SDE support: no
DODS support: no
SQLite support: yes
DWGdirect support no
PANORAMA GIS support: no
INFORMIX DataBlade support: no
GEOS support: yes
Old-gen python no
SWIG Bindings: no
Statically link PROJ.4: no
enable OGR building: yes
enable pthread support: no
hide internal symbols: no</CODE></PRE><PRE><CODE><FONT face=Arial size=2>Regards,</FONT></CODE></PRE><PRE><CODE><FONT face=Arial size=2>Marco</FONT></CODE></PRE></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr> </DIV></DIV>
<DIV dir=ltr><BR>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>Da:</B> Paul Kelly [mailto:paul-grass@stjohnspoint.co.uk]<BR><B>Inviato:</B> sab 10/05/2008 12.11<BR><B>A:</B> marco.pasetti@alice.it<BR><B>Cc:</B> grass-dev@lists.osgeo.org<BR><B>Oggetto:</B> Re: R: [GRASS-dev] WinGRASS: r.in.gdal projection issue<BR></FONT><BR></DIV>
<DIV>
<P><FONT size=2>On Fri, 9 May 2008 marco.pasetti@alice.it wrote:<BR><BR>> Hi all,<BR>><BR>> I also tried to launch GRASS with GISBase without spaces (such as C:\GISBASE), with the demolocation, and placing the input file (Z_39_3.ASC) in a path without spaces too (C:\), but I got the same result.<BR><BR>Hello Marco,<BR>How big is the file? Would it be possible to put it somewhere for me to<BR>have a look at and try and reproduce the behaviour?<BR><BR>FWIW, other things to try would be to see if other commands can read the<BR>projection info from the file. E.g.<BR>gdalinfo filename<BR>g.proj -p georef=filename<BR><BR>Paul<BR></FONT></P></DIV></BODY></HTML>