I installed grass 6.2.1 from the ubuntu packages provide by les-ejk
and I'm trying to create a location from a MODIS geotiff with no
success. I select the tif file, click on create location, a black
screen quickly flashes and Grass tells me it needs to shut down. But
the location is not created.
The geotiff was created using the Modis Reprojection Tool. gdalinfo on
the geotiff file gives me:
daniel@mancalinha:~/modis$ gdalinfo importada/lai.2007.01.01.Lai_1km.tif
Driver: GTiff/GeoTIFF
Size is 1307, 1196
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235629972,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-60.921414118540142,-0.004182589339916)
Pixel Size = (0.008365178679831,-0.008365178679831)
Metadata:
AREA_OR_POINT=Point
TIFFTAG_SOFTWARE=MODIS Reprojection Tool v3.3a January 2006
Corner Coordinates:
Upper Left ( -60.9214141, -0.0041826) ( 60d55'17.09"W, 0d 0'15.06"S)
Lower Left ( -60.9214141, -10.0089363) ( 60d55'17.09"W, 10d 0'32.17"S)
Upper Right ( -49.9881256, -0.0041826) ( 49d59'17.25"W, 0d 0'15.06"S)
Lower Right ( -49.9881256, -10.0089363) ( 49d59'17.25"W, 10d 0'32.17"S)
Center ( -55.4547699, -5.0065594) ( 55d27'17.17"W, 5d 0'23.61"S)
Band 1 Block=1307x1 Type=Byte, ColorInterp=Gray
daniel@mancalinha:~/modis$
Any clues?
If I try the same with the original hdf file the same problem happens
I just noticed that grass is looking for the epsg file in
/usr/local/share/proj and it's installed in /usr/share/proj. Could
that be related? Also, the geotiff has several epsg codes in it and
the only one valid is 4326. If I create a location with this EPSG code
all goes fine. If I use the others, a tmp.<somenumber> directory is
created.
Thanks
Daniel
On 6/16/07, Daniel Victoria <daniel.victoria@gmail.com> wrote:
Hi all,
I installed grass 6.2.1 from the ubuntu packages provide by les-ejk
and I'm trying to create a location from a MODIS geotiff with no
success. I select the tif file, click on create location, a black
screen quickly flashes and Grass tells me it needs to shut down. But
the location is not created.
The geotiff was created using the Modis Reprojection Tool. gdalinfo on
the geotiff file gives me:
daniel@mancalinha:~/modis$ gdalinfo importada/lai.2007.01.01.Lai_1km.tif
Driver: GTiff/GeoTIFF
Size is 1307, 1196
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235629972,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-60.921414118540142,-0.004182589339916)
Pixel Size = (0.008365178679831,-0.008365178679831)
Metadata:
AREA_OR_POINT=Point
TIFFTAG_SOFTWARE=MODIS Reprojection Tool v3.3a January 2006
Corner Coordinates:
Upper Left ( -60.9214141, -0.0041826) ( 60d55'17.09"W, 0d 0'15.06"S)
Lower Left ( -60.9214141, -10.0089363) ( 60d55'17.09"W, 10d 0'32.17"S)
Upper Right ( -49.9881256, -0.0041826) ( 49d59'17.25"W, 0d 0'15.06"S)
Lower Right ( -49.9881256, -10.0089363) ( 49d59'17.25"W, 10d 0'32.17"S)
Center ( -55.4547699, -5.0065594) ( 55d27'17.17"W, 5d 0'23.61"S)
Band 1 Block=1307x1 Type=Byte, ColorInterp=Gray
daniel@mancalinha:~/modis$
Any clues?
If I try the same with the original hdf file the same problem happens
I just noticed that grass is looking for the epsg file in
/usr/local/share/proj and it's installed in /usr/share/proj. Could
that be related?
You should put the path to where your "epsg" file resides into the
"Path to the EPSG codes file" dialog. The file comes with PROJ.4.
Also, the geotiff has several epsg codes in it and
the only one valid is 4326.
Your gdalinfo output seems corrupted. Compare with the output of
"epsg_tr.py 4326" to see the difference. It is missing EPSG codes for
prime meridian and unit. That's propably why GRASS fails to creata a
location using it.
Also please note that only the last EPSG code designates the coordinate
system. Those coming before it are for the CS components, like
spheroid, datum etc.
If I create a location with this EPSG code
all goes fine. If I use the others, a tmp.<somenumber> directory is
created.
That's because you are passing wrong EPSG coordinate system code, which
GRASS 6.2 can't handle and creates a corrupted location. As expected,
only "4326" works OK. In GRASS 6.3, using an invalid coordinate system
EPSG code yields an error message like on the attached screendump.
Please fill in a wish in the tracker for backporting this usefull
feature to GRASS 6.2. Maybe it will get included in 6.2.2 release.
I just noticed that grass is looking for the epsg file in
/usr/local/share/proj and it's installed in /usr/share/proj. Could
that be related? Also, the geotiff has several epsg codes in it and
the only one valid is 4326. If I create a location with this EPSG code
all goes fine. If I use the others, a tmp.<somenumber> directory is
created.
the default dir should be set by the packager before compile with:
./configure --with-proj-share=
Or by the user at run time with "Define new location with... -> Path to
the EPSG-codes file"
But note that's only for the startup GUI EPSG search tool (AFAIK), not
the rest of GRASS. The rest of GRASS (ie "g.proj -c" used by
make_location_epsg.sh.in) depends on PROJ.4 being installed properly.
(so e.g. GDAL can find it too)
