Hi,<br><br>some month ago I had to make some maps with the HWSD database. I used the mdbtool (<a href="http://mdbtools.sourceforge.net/">http://mdbtools.sourceforge.net/</a>) to extract all tables from the .mdb file to some sql files. Then I just imported all the table in a postgres database, and it is quite easy to access them in GRASS ...<br>
<br><br>Regards,<br>Sylvain<br><br><br><div class="gmail_quote">2011/3/24 Luigi Ponti <span dir="ltr">&lt;<a href="mailto:lponti@inbox.com">lponti@inbox.com</a>&gt;</span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Dear listers,<br>
<br>
I seem to have successfully imported the Harmonized World Soil Database (HWSD) in BIL format to GRASS 6.4.svn via r.in.gdal. The data is available at<br>
<a href="http://www.iiasa.ac.at/Research/LUC/luc07/External-World-soil-database/HTML/index.html" target="_blank">http://www.iiasa.ac.at/Research/LUC/luc07/External-World-soil-database/HTML/index.html</a><br>
<br>
The HWSD is composed of a raster image file and a linked attribute database stored in Microsoft Access 2003 format. According to the documentation (page 28)<br>
<a href="http://www.iiasa.ac.at/Research/LUC/External-World-soil-database/HWSD_Documentation.pdf" target="_blank">http://www.iiasa.ac.at/Research/LUC/External-World-soil-database/HWSD_Documentation.pdf</a><br>
<br>
using the HWSD database in a GIS should be straightforward. The MU_GLOBAL attribute column in the database is also stored in each raster cell, and may therefore be considered as the link to the raster.<br>
<br>
My goal would be to obtain more raster maps that show other (soil) attributes from the database or combination of different attributes based on computations. The documentation suggests the following work flow:<br>
<br>
- if necessary, realize the appropriate calculations (ex: after exporting from Access to Excel);<br>
- convert final attributes table to a compatible GIS format;<br>
- join the MU_GLOBAL attribute and the GRID value (dbf or txt formats);<br>
- convert the attribute to a new GRID (in the case it is needed).<br>
<br>
What comes to mind is value replacement using r.recode but maybe there is a better course of action.<br>
<br>
This is seemingly the best available soil GIS database with global coverage, so I thought my question may be useful to others.<br>
<br>
Kind regards and thanks in advance for any hint,<br>
<br>
Luigi<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote></div><br>