tag:blogger.com,1999:blog-43007200093185607262024-02-08T17:17:31.611+03:00Land Use Planning - Arumeru District -TanzaniaUnknownnoreply@blogger.comBlogger10125tag:blogger.com,1999:blog-4300720009318560726.post-85810101587374903842009-06-18T18:59:00.001+03:002009-06-18T18:59:44.408+03:00WORLDCLIM Data<div xmlns='http://www.w3.org/1999/xhtml'><span style='font-family: sans-serif;'>All WORLDCLIM data for Tanzania (Zone 36 and 37) downloaded and imported in GRASS.<br/><br/>WorldClim is a set of global climate layers (climate grids) with a spatial resolution of a square kilometer. The current version is Version 1.4 (release 3).<br/><br/>Some notes: tiles are in Generic grid (raster) format. EPSG:4326. Spatial resolutions: 30 seconds (0.93 x 0.93 = 0.86 km<sup>2 </sup>at the equator). Data files (.BIL) are sequential binary files in which values are stored line by line from the upper-left to the lower-right corner. Each cell (pixel) is an integer value (2 bytes; -32768..32767). If you simply open the data as images or raster data sets, there is a problem. ESRI software assumes that BIL files do not have negative values. These values (x) are replaced by (65535 - x); E.g., -10 becomes 65525. The nodata value of -9999 is not recognized... a BIG mess results... I got the same problem with GRASS. For a workaround:<br/>Use <a href='http://www.diva-gis.org/'>DIVA-GIS</a> (free download)<br/>Data\Import to gridfile\Multiple Files (BIL/BIP/BSQ) and then<br/>Data\Export to gridfile\Multiple Files (ESRI Ascii)<br/>Now GRASS works, with module r.in.arc. Or Qgis GRASS plugin, r.in.arc module.<br/><br/>Informations are:<br/><b>tmean</b> = average monthly mean temperature (°C * 10)<br/><b>tmin</b> = average monthly minimum temperature (°C * 10)<br/><b>tmax</b> = average monthly maximum temperature (°C * 10)<br/><b>prec</b> = average monthly precipitation (mm)<br/><b>bio</b> = <a href='http://www.worldclim.org/bioclim.htm'>bioclimatic variables</a> derived from the tmean, tmin, tmax and prec<br/><b>alt</b> = altitude (elevation above sea level) (m) (from <a href='http://www2.jpl.nasa.gov/srtm/'>SRTM</a>)<br/></span><span style='font-family: sans-serif;'><br/>Note that temperature is expressed as degree C * 10.<br/><br/><b>BIOCLIM</b><br/>Bioclimatic variables are derived from the monthly temperature and rainfall values in order to generate more biologically meaningful variables. These are often used in ecological niche modeling (e.g., BIOCLIM, GARP). The bioclimatic variables represent annual trends (e.g., mean annual temperature, annual precipitation) seasonality (e.g., annual range in temperature and precipitation) and extreme or limiting environmental factors (e.g., temperature of the coldest and warmest month, and precipitation of the wet and dry quarters). A quarter is a period of three months (1/4 of the year).<br/>They are coded as follows:<br/>BIO1 = Annual Mean Temperature<br/>BIO2 = Mean Diurnal Range (Mean of monthly (max temp - min temp))<br/>BIO3 = Isothermality (P2/P7) (* 100)<br/>BIO4 = Temperature Seasonality (standard deviation *100)<br/>BIO5 = Max Temperature of Warmest Month<br/>BIO6 = Min Temperature of Coldest Month<br/>BIO7 = Temperature Annual Range (P5-P6)<br/>BIO8 = Mean Temperature of Wettest Quarter <br/>BIO9 = Mean Temperature of Driest Quarter<br/>BIO10 = Mean Temperature of Warmest Quarter<br/>BIO11 = Mean Temperature of Coldest Quarter<br/>BIO12 = Annual Precipitation<br/>BIO13 = Precipitation of Wettest Month<br/>BIO14 = Precipitation of Driest Month<br/>BIO15 = Precipitation Seasonality (Coefficient of Variation)<br/>BIO16 = Precipitation of Wettest Quarter<br/>BIO17 = Precipitation of Driest Quarter<br/>BIO18 = Precipitation of Warmest Quarter<br/>BIO19 = Precipitation of Coldest Quarter<br/><br/>This scheme follows that of ANUCLIM, except that for temperature seasonality the standard deviation was used because a coefficient of variation does not make sense with temperatures between -1 and 1).<br/><br/><span style='font-weight: bold;'>Links:</span><br/>WORLDCLIM Docs:<br/><a href='http://www.worldclim.org/methods.htm'>http://www.worldclim.org/methods.htm</a><br/><a href='http://www.worldclim.org/format.htm'>http://www.worldclim.org/format.htm</a><br/><a href='http://www.worldclim.org/bioclim.htm'>http://www.worldclim.org/bioclim.htm</a><br/><a href='http://www.diva-gis.org/'>DIVA-GIS</a><br/>WORLDCLIM Data:<br/><a href='http://www.worldclim.org/tiles.php'>http://www.worldclim.org/tiles.php</a></span><br/><br/><br /><br /><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-89945040766627746502009-06-18T13:25:00.001+03:002009-06-18T15:35:09.264+03:00AFRICOVER Data<div xmlns='http://www.w3.org/1999/xhtml'>All AFRICOVER data for Tanzania downloaded and imported in GRASS:<br/><br/> • tz_boundaries<br/> ⁃ tz_politicalboundary<br/> ⁃ tz_adm<br/><br/> • tz_towns<br/> ⁃ tz_majortowns<br/> ⁃ tz_othertowns<br/><br/> • tz_roads<br/> • tz_rv<br/> • tz_landform_agg<br/> • tz_lithology_agg<br/> • tz_spatial_agg <br/> • tz_cultiv_agg <br/> • tz_grass_agg<br/> • tz_woody_agg<br/><span style='font-weight: bold;'><br/></span>All vector in shape format, EPSG:4326. No problem with GRASS Qgis plugin, Module v.in.ogr.qgis.<br/>Database uses sqlite driver, set up by GRASS Qgis plugin, Module db.connect.<br/>AFRICOVER.sqlite created BEFORE importing data, with sqlite-manager extension for firefox.<br/><span style='font-weight: bold;'><br/>Links:</span><br/>AFRICOVER: <a href='http://www.africover.org/'>http://www.africover.org/</a><br/>Sqlite-Manager: <a href='http://code.google.com/p/sqlite-manager/'>http://code.google.com/p/sqlite-manager/</a><br/>Mozilla Firefox: <a href='http://www.mozilla.com/en-US/'>http://www.mozilla.com/en-US/</a><br/><br/><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p><br/><br/></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-57922465308653975692009-06-18T10:29:00.003+03:002009-06-18T10:55:21.023+03:00SRTM DEM<div xmlns='http://www.w3.org/1999/xhtml'>All SRTM tiles for Tanzania downloaded and imported in GRASS (EPSG:4326):<br/><br/><span style='text-decoration: underline;'>42_13 </span>|<span style='text-decoration: underline;'> 43_13 </span>|<span style='text-decoration: underline;'> 44_13</span><br/><span style='text-decoration: underline;'>42_14</span> | <span style='text-decoration: underline;'>43_14 </span>|<span style='text-decoration: underline;'> 44_14</span><br/>42_15 | 43_15 | 44_15 |45_15<br/><br/>All tiles patched by:<br/>r.patch input="srtm_42_13,srtm_42_14,srtm_42_15,srtm_43_13,srtm_43_14,srtm_43_15,srtm_44_13,srtm_44_14,srtm_44_15,srtm_45_15" output="srtm_full"<br/><br/><span style='font-weight: bold;'>Links:</span><br/>SRTM (Original dataset Docs): <a href='http://www2.jpl.nasa.gov/srtm/index.html' target='_blank'>http://www2.jpl.nasa.gov/srtm/index.html</a><br/>SRTM (V4): <a href='http://srtm.csi.cgiar.org/' target='_blank'>http://srtm.csi.cgiar.org/</a><br/><a href='ftp://xftp.jrc.it/pub/srtmV4/tiff/srtm_44_13.zip'><br/></a><p class='scribefire-powered'><a href='ftp://xftp.jrc.it/pub/srtmV4/tiff/srtm_44_13.zip'>Powered by </a><a href='http://www.scribefire.com/'>ScribeFire</a>.</p><br/><br/><br/></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-21940874135035908972009-06-18T09:44:00.001+03:002009-06-18T09:44:33.460+03:00Updates<div xmlns='http://www.w3.org/1999/xhtml'>Switched to:<br/><ul><li>GRASS 6.4 RC5</li><li>Qgis 1.1.0</li></ul>QGIS/GRASS and Frameworks Mac binaries: <a href='http://www.kyngchaos.com/software:unixport' target='_blank'>http://www.kyngchaos.com/software:unixport</a><br/><br/><br/><br /><br /><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-84885153508291584092009-03-31T13:48:00.001+03:002009-03-31T13:48:25.746+03:00Reprojecting data<div xmlns='http://www.w3.org/1999/xhtml'>Another map has been were georeferenced in WGS84LL using Georeferencing Qgis Plugin:<br/><ul><li>"<a href='http://www.agriculture.go.tz/Research-Training/soil-agro-map.htm#agroeco'>Agroecological zones</a> of Tanzania",<span style='font-family: monospace;'> </span>Ministry Of Agriculture, Food and Cooperative</li></ul>Some reprojection has been done on raster maps using gdal tools:<br/><br/><span style='font-weight: bold;'>No Datum change:</span><br style='font-weight: bold;'/><br/><blockquote>gdalwarp -t_srs EPSG:32737 4326_HRImage.tif 32737_HRImage.tif<br/></blockquote><br/>move the source WGS84LL source GeoTiff file to the WGS84UTM37 huge, uncompressed target file, so you need to reduce it with <br/><br/><blockquote>gdal_translate -co "COMPRESS=JPEG" -co "PHOTOMETRIC=YCBCR" -co "JPEG_QUALITY=100" 32737_HRImage.tif 32737_HRImage2.tif<br/></blockquote>and add pyramid overviews (load very fast in Qgis) with <br/><blockquote>gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL 32737_HRImage2.tif 2 4 8 16<br/></blockquote><br/><span style='font-weight: bold;'>Datum change:</span><br/><br/><blockquote>gdalwarp -t_srs EPSG:32737 -s_srs "1285_TZ.prf" 21037_Topo.tif 32737_Topo.tif<br/></blockquote>move the source ARC60UTM37 source GeoTiff file to the WGS84UTM37 file, using the right transformation parameters written in:<br/><blockquote><span style='font-weight: bold;'>1285_TZ.prf</span><br/>+proj=utm +zone=37 +south +a=6378249.145 +b=6356514.96582849 +units=m +towgs84=-175,-23,-303<br/></blockquote>The process works backward: from standard EPSG_code --> to custom_CRS<br/><br/>In this case you can add right info to the header of the result with:<br/><blockquote>gdal_translate -a_srs "21037_WK.txt" <span style='font-style: italic;'>source</span>.tif 21037_<span style='font-style: italic;'>target</span>.tif</blockquote>with <span style='font-weight: bold;'>21037_WK.txt</span>:<br/><blockquote>PROJCS[" Projection Name = Arc_1960_UTM_Zone_37S Units = meters GeoTIFF Units = meters",<br/> GEOGCS["Arc 1960",<br/> DATUM["Arc_1960",<br/> SPHEROID["Clarke 1880 (RGS)",6378249.145,293.4649999999983,<br/> AUTHORITY["EPSG","7012"]],<br/> AUTHORITY["EPSG","6210"]],<br/> PRIMEM["Greenwich",0],<br/> UNIT["degree",0.0174532925199433],<br/> AUTHORITY["EPSG","4210"]],<br/> PROJECTION["Transverse_Mercator"],<br/> PARAMETER["latitude_of_origin",0],<br/> PARAMETER["central_meridian",39],<br/> PARAMETER["scale_factor",0.9996],<br/> PARAMETER["false_easting",500000],<br/> PARAMETER["false_northing",10000000],<br/> UNIT["metre",1,<br/> AUTHORITY["EPSG","9001"]]]<br/></blockquote>that can be extract with gdalinfo>info.txt on the original toposheets (see <a href='http://luparumeru.blogspot.com/2009/03/defining-area-of-intervention.html' target='_blank'>Defining the Area of Intervention</a>)<br/><br/>In this way we obtained:<br/><br/>./Raster/21037:<br/>21037_Agro_Eco_Zones.tif<br/>21037_Arusha_soil.tif<br/>21037_EcoClimatic_Zones.tif<br/>21037_HRImage.tif<br/>21037_Topo.tif<br/><br/>The same structure is repeated for 32737 and 4326. <br/><br/>With<br/><blockquote>gdaltindex index.shp *.tif</blockquote>we obtained the coverage of all the rasters in every folder.<br/><br/>The same transformation procedures work with vector data:<br/><blockquote>ogr2ogr 4326_Index.shp 21037_Index.shp -t_srs EPSG:4326 -s_srs "1285.prf"<br/></blockquote>or<br/><blockquote>ogr2ogr mosaic2.shp mosaic.shp -t_srs EPSG:4326 -s_srs "+init=EPSG:21037 +wgs84=-175,-23,-303"<br/></blockquote><span style='font-weight: bold;'>Links:</span><br/><a href='http://www.gdal.org/gdalwarp.html' target='_blank'>http://www.gdal.org/gdalwarp.html</a><br/><br/><div class='zemanta-pixie'><img src='http://img.zemanta.com/pixy.gif?x-id=6b6ea6b8-3a16-8afb-a3fb-7b3428d26bbf' class='zemanta-pixie-img'/></div><br /><br /><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-62688961981870336582009-03-27T16:02:00.001+03:002009-03-27T16:02:17.432+03:00Georeferencing raster data<div xmlns='http://www.w3.org/1999/xhtml'>Two maps:<br/><ul><li>"Soils of Arusha and Monduli Area 1:250.000" (Soil Research Institute, Ottawa, Canada, 1974)</li><li>"Tanzania Eco-Climatic Zones" (Institute of Resource Assestment). <br/></li></ul>were georeferenced in WGS84LL using Georeferencing Qgis Plugin.<br/><br/>Other sources can be georeferenced using gdal utilities.<br/>For testing purposes is possible to download high resolution imagery covering part of Area of Intervention from Google dataset using GoogleMV. Level 17 has enough resolution.<br/>The result is a big *.jpg file and a *.map file (Ozi Explorer format, see <a href='http://www.rus-roads.ru/gps/help_ozi/map_file_format.html' target='_blank'>here</a> info on the format).<br/>The *.map is simple text, and coordinates of the corners of the image in latlongwgs84 can be found near the end of the file:<br/><blockquote><span style='color: rgb(255, 0, 0);'>MMPLL,1, 36.754167, -3.076833</span><br/>MMPLL,2, 36.903752, -3.076833<br/><span style='color: rgb(255, 0, 0);'>MMPLL,3, 36.903752, -3.233316</span><br/>MMPLL,4, 36.754167, -3.233316<br/></blockquote>(in red are marked the Upper Left and Lower Right corners coordinates).<br/><br/><a target='_blank' href='http://www.gdal.org/gdal_translate.html'>gdal_translate</a> can be used now:<br/><br/>gdal_translate -a_srs EPSG:4326 -a_ullr 36.754167 -3.076833 36.903752 -3.233316 -co "COMPRESS=JPEG" -co "PHOTOMETRIC=YCBCR" -co "JPEG_QUALITY=100" <span style='font-style: italic;'>source</span>.jpg <span style='font-style: italic;'>target</span>.tif<br/><br/><ul><li>the -a_srs parameter assign an EPSG code to the image</li><li>the -a_ullr parameter define upper left and lower right corners coordinates (taken from the *.map file)</li><li>the -co parameters define the compression algoritm used</li></ul>The result is a GeoTIFF with this header (read it with <a target='_blank' href='http://www.gdal.org/gdalinfo.html'>gdalinfo</a>):<br/><br/><blockquote>Driver: GTiff/GeoTIFF<br/>Files: target.tif<br/>Size is 13943, 14608<br/>Coordinate System is:<br/>GEOGCS["WGS 84",<br/> DATUM["WGS_1984",<br/> SPHEROID["WGS 84",6378137,298.2572235630016,<br/> AUTHORITY["EPSG","7030"]],<br/> AUTHORITY["EPSG","6326"]],<br/> PRIMEM["Greenwich",0],<br/> UNIT["degree",0.0174532925199433],<br/> AUTHORITY["EPSG","4326"]]<br/>Origin = (36.754167000000002,-3.076833000000000)<br/>Pixel Size = (0.000010728322456,-0.000010712144031)<br/>Metadata:<br/> AREA_OR_POINT=Area<br/>Image Structure Metadata:<br/> SOURCE_COLOR_SPACE=YCbCr<br/> COMPRESSION=JPEG<br/> INTERLEAVE=PIXEL<br/>Corner Coordinates:<br/>Upper Left ( 36.7541670, -3.0768330) ( 36d45'15.00"E, 3d 4'36.60"S)<br/>Lower Left ( 36.7541670, -3.2333160) ( 36d45'15.00"E, 3d13'59.94"S)<br/>Upper Right ( 36.9037520, -3.0768330) ( 36d54'13.51"E, 3d 4'36.60"S)<br/>Lower Right ( 36.9037520, -3.2333160) ( 36d54'13.51"E, 3d13'59.94"S)<br/>Center ( 36.8289595, -3.1550745) ( 36d49'44.25"E, 3d 9'18.27"S)<br/>Band 1 Block=13943x16 Type=Byte, ColorInterp=Red<br/>Band 2 Block=13943x16 Type=Byte, ColorInterp=Green<br/>Band 3 Block=13943x16 Type=Byte, ColorInterp=Blue<br/></blockquote><br/>The file is huge, so pyramid layers can be added to fasten the opening in GIS with <a target='_blank' href='http://www.gdal.org/gdaladdo.html'>gdaladdo</a>:<br/><br/><blockquote>gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL <span style='font-style: italic;'>target</span>.tif 2 4 8 16<br/></blockquote><br/>In this way the last lines of the header of GeoTIFF are substituted with:<br/><br/><blockquote>Band 1 Block=13943x16 Type=Byte, ColorInterp=Red<br/> Overviews: 6972x7304, 3486x3652, 1743x1826, 872x913<br/>Band 2 Block=13943x16 Type=Byte, ColorInterp=Green<br/> Overviews: 6972x7304, 3486x3652, 1743x1826, 872x913<br/>Band 3 Block=13943x16 Type=Byte, ColorInterp=Blue<br/> Overviews: 6972x7304, 3486x3652, 1743x1826, 872x913<br/></blockquote>You can obtain a shapefile with the coverage of the GeoTiff with:<br/><ul><li>gdaltindex <span style='font-style: italic;'>target</span>.shp <span style='font-style: italic;'>target</span>.tif</li></ul><span style='font-weight: bold;'>Links:</span><br/><br/>GoogleMV:<br/><ul><li><a target='_blank' href='http://www.silber2004.narod.ru/trevel.htm'>http://www.silber2004.narod.ru/trevel.htm</a> (this link is in Russian, you can download the program <a target='_blank' href='http://www.silber2004.narod.ru/GoogleMV1.zip'>here</a>)</li></ul>GDAL Utilities:<br/><ul><li><a target='_blank' href='http://www.gdal.org/gdal_utilities.html'>http://www.gdal.org/gdal_utilities.html</a></li></ul>OZI *.map format:<br/><ul><li><a href='http://www.rus-roads.ru/gps/help_ozi/map_file_format.html' target='_blank'>http://www.rus-roads.ru/gps/help_ozi/map_file_format.html</a></li></ul><br/><br/><div class='zemanta-pixie'><img src='http://img.zemanta.com/pixy.gif?x-id=8ba64bee-d6e2-8b9e-bcdb-da794d00b61d' class='zemanta-pixie-img'/></div><br /><br /><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-53350994906433687972009-03-25T15:01:00.001+03:002009-03-25T15:01:47.045+03:00Defining the Area of Intervention<div xmlns='http://www.w3.org/1999/xhtml'>The Project AI (Area of Intervention) is defined by four of the Standard Topo Sheets "East Africa 1:50.000 (United Republic of Tanzania)":<br/><ul><li>55/1 Oldonyo Sambu</li><li>55/2 Engare Nanyuki</li><li>55/3 Arusha</li><li>55/4 Usa River.</li></ul>These maps (geotiff) are georeferenced in Standard Arc 1960 UTM (EPSG:21037).<br/><br/>To obtain a single raster mosaic covering all the area:<br/><ol><li>ls -1 *.tif > tiff_list.txt</li><li>gdal_merge.py -v -o mosaic.tif --optfile tiff_list.txt</li></ol>Then with:<br/><ul><li>gdalinfo mosaic.tif</li></ul>You can obtain:<br/><br/>Driver: GTiff/GeoTIFF<br/>Files: mosaic.tif<br/>Size is 13109, 13042<br/>Coordinate System is:<br/>PROJCS[" Projection Name = Arc_1960_UTM_Zone_37S Units = meters GeoTIFF Units = meters",<br/> GEOGCS["Arc 1960",<br/> DATUM["Arc_1960",<br/> SPHEROID["Clarke 1880 (RGS)",6378249.145,293.4649999999983,<br/> AUTHORITY["EPSG","7012"]],<br/> AUTHORITY["EPSG","6210"]],<br/> PRIMEM["Greenwich",0],<br/> UNIT["degree",0.0174532925199433],<br/> AUTHORITY["EPSG","4210"]],<br/> PROJECTION["Transverse_Mercator"],<br/> PARAMETER["latitude_of_origin",0],<br/> PARAMETER["central_meridian",39],<br/> PARAMETER["scale_factor",0.9996],<br/> PARAMETER["false_easting",500000],<br/> PARAMETER["false_northing",10000000],<br/> UNIT["metre",1,<br/> AUTHORITY["EPSG","9001"]]]<br/>Origin = (222098.049999999988358,9668234.250000000000000)<br/>Pixel Size = (4.250000000000000,-4.250000000000000)<br/>Metadata:<br/> AREA_OR_POINT=Area<br/>Image Structure Metadata:<br/> INTERLEAVE=PIXEL<br/>Corner Coordinates:<br/>Upper Left ( 222098.050, 9668234.250) ( 36d30'0.01"E, 2d59'56.28"S)<br/>Lower Left ( 222098.050, 9612805.750) ( 36d29'55.57"E, 3d29'59.98"S)<br/>Upper Right ( 277811.300, 9668234.250) ( 37d 0'3.48"E, 3d 0'0.01"S)<br/>Lower Right ( 277811.300, 9612805.750) ( 36d59'59.92"E, 3d30'4.33"S)<br/>Center ( 249954.675, 9640520.000) ( 36d44'59.74"E, 3d15'0.26"S)<br/>Band 1 Block=13109x1 Type=Byte, ColorInterp=Red<br/>Band 2 Block=13109x1 Type=Byte, ColorInterp=Green<br/>Band 3 Block=13109x1 Type=Byte, ColorInterp=Blue<br/><br/>These are the limits (and other useful information) of the AI in Arc60, UTM and LatLong.<br/><br/>You can obtain a shapefile covering AI (EPSG:21037)with:<br/><ul><li>gdaltindex mosaic.shp mosaic.tif</li></ul>Now you can add a *.prj file to the shape with:<br/><ul><li>ogr2ogr mosaic2.shp mosaic.shp -a_srs EPSG:21037</li></ul>And read parameters with:<br/><ul><li>ogrinfo -al -so mosaic2.shp<br/></li></ul><span style='font-weight: bold;'>Links:</span><br style='font-weight: bold;'/><br/>GDAL:<br/><ul><li><a href='http://www.gdal.org/' target='_blank'>http://www.gdal.org/</a></li><li><a href='http://www.gdal.org/gdal_utilities.html' target='_blank'>http://www.gdal.org/gdal_utilities.html</a></li><li><a href='http://trac.osgeo.org/gdal/wiki/UserDocs' target='_blank'>http://trac.osgeo.org/gdal/wiki/UserDocs</a></li></ul>OGR:<br/><ul><li><a href='http://www.gdal.org/ogr/index.html' target='_blank'>http://www.gdal.org/ogr/index.html</a></li><li><a href='http://www.gdal.org/ogr/index.html' target='_blank'>http://www.gdal.org/ogr_utilities.html</a><br/></li></ul><br/><br/><div class='zemanta-pixie'><img src='http://img.zemanta.com/pixy.gif?x-id=b9eb359e-e193-45cd-b1bc-5945421d9eee' class='zemanta-pixie-img'/></div><br /><br /><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-48915614867247110452009-03-24T12:35:00.001+03:002009-03-24T12:35:21.755+03:00EPSG Codes and GRASS LOCATIONS<div xmlns='http://www.w3.org/1999/xhtml'><b>Software used:<br/></b><br/>Mac OS X 10.5.6 on Intel Macbook (late 2007)<br/>Quantum GIS (1.0.1-Kore)<br/>GRASS GIS (6.4.0RC3)<br/><br/><br/><b>CRS used:</b><ul><li>EPSG:4326 - WGS84LL</li><li>EPSG:32737 - WGS84UTM37S</li><li>EPSG:21037 - ARC60UTM37S</li><li>EPSG:4210 - ARC60LL</li><li>Custom ARC60UTM37S (Kenya_Tanzania, based on EPSG:1122): +proj=utm +zone=37 +south +ellps=clrk80 +units=m +towgs84=-160,-6,-302</li><li>Custom ARC60UTM37S (Tanzania Onshore and Offshore, based on EPSG:1285): +proj=utm +zone=37 +south +a=6378249.145 +b=6356514.96582849 +units=m +towgs84=-175,-23,-303</li></ul><span style='font-weight: bold;'>Projection Parameters</span> (as written in the Standard Tanzanian Topo Sheets):<br/><br/><ul><li>Grid: UTM Zone 37</li><li>Projection: Transverse Mercator</li><li>Spheroid: Clarke 1880 (Modified)</li><li>Unit of Measurement: metre</li><li>Meridian of Origin: 39°00' East of Greenwich</li><li>Latitude of Origin: Equator</li><li>Scale Factor at Origin: 0.9996</li><li>False Coords of Origin: 500.000 Easting; 10.000.000 Northing</li><li>Datum: New (1960) Arc<br/></li></ul><b>GRASS Locations</b><br/><br/>Standard EPSG 4326, 21037, 32737<br/><br/><ol><li>In QGIS/GRASS Plugin Create new mapset</li><li>Choose the database (path_to_GIS_data)</li><li>Create New Location: WGS84LL ("Projected": EPSG:4326)*<br/></li><li>Default Wind: Tanzania <br/></li><li>Choose a name for the Mapset (ARUMERU)</li></ol>* Repeat for WGS84UTM37S (EPSG:32737) and for ARC60UTM37S (EPSG:21037)<br/><br/><b>Links:</b><br/>QGIS: <a target='_blank' href='http://www.qgis.org/'>http://www.qgis.org/</a><br/>GRASS GIS: <a target='_blank' href='http://grass.itc.it/'>http://grass.itc.it/</a><br/>QGIS/GRASS and Frameworks Mac binaries: <a target='_blank' href='http://www.kyngchaos.com/software:unixport'>http://www.kyngchaos.com/software:unixport</a><br/>EPSG Codes: <a target='_blank' href='http://www.epsg-registry.org/'>http://www.epsg-registry.org/</a><br/><br/><p class='scribefire-powered'>Powered by <a href='http://www.scribefire.com/'>ScribeFire</a>.</p><br/><br/><div class='zemanta-pixie'><img src='http://img.zemanta.com/pixy.gif?x-id=924b6ed5-5edf-4758-98d4-4ca95e029713' class='zemanta-pixie-img'/></div></div>Unknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-52661953342545543202009-03-24T10:00:00.009+03:002009-03-24T12:37:02.215+03:00Village Land DemarcationVillage land demarcation activity was done to define the village administrative boundaries. This is important step before starting the village land use planning exercise. The challenges of village Land <span class="blsp-spelling-corrected" id="SPELLING_ERROR_0">demarcation</span> activities were <span class="blsp-spelling-corrected" id="SPELLING_ERROR_1">immense</span>. The land issue in the <span class="blsp-spelling-error" id="SPELLING_ERROR_2"><span class="blsp-spelling-error" id="SPELLING_ERROR_0">Meru</span></span> and <span class="blsp-spelling-error" id="SPELLING_ERROR_3"><span class="blsp-spelling-error" id="SPELLING_ERROR_1">Maasai</span></span> communities is very sensitive due to complexity and historical issues related to land tenure. The villagers decided the boundaries of the Villages to follow the colonial settler farms that were in the old map projections (<span class="blsp-spelling-error" id="SPELLING_ERROR_4"><span class="blsp-spelling-error" id="SPELLING_ERROR_2">TTM</span></span>), so we had to <span class="blsp-spelling-corrected" id="SPELLING_ERROR_5">re-project</span> the geographical data using free <span class="blsp-spelling-error" id="SPELLING_ERROR_6"><span class="blsp-spelling-error" id="SPELLING_ERROR_3">GIS</span></span> software to <span class="blsp-spelling-error" id="SPELLING_ERROR_7"><span class="blsp-spelling-error" id="SPELLING_ERROR_4">UTM</span></span> ARC 60. We had to seek for data related to the boundaries of the <span class="blsp-spelling-corrected" id="SPELLING_ERROR_8">settler</span> farms, West Kilimanjaro Ranch, <span class="blsp-spelling-error" id="SPELLING_ERROR_9"><span class="blsp-spelling-error" id="SPELLING_ERROR_5">Arusha</span></span> National Park and the <span class="blsp-spelling-error" id="SPELLING_ERROR_10"><span class="blsp-spelling-error" id="SPELLING_ERROR_6">Meru</span></span> District Council boundaries as stipulated in the Government Notice (<span class="blsp-spelling-error" id="SPELLING_ERROR_11"><span class="blsp-spelling-error" id="SPELLING_ERROR_7">GN</span></span>). Getting all this old data was a big <span class="blsp-spelling-corrected" id="SPELLING_ERROR_12">challenge</span> and it takes a lot of time and other resources. After all the challenges and with right projections of the coordinates, the beacons were set on the village boundary points. The village leaders, village land adjudication committees, the district surveyor and representatives from <span class="blsp-spelling-corrected" id="SPELLING_ERROR_13">institutions</span> such as West Kilimanjaro Ranch were involved. out of this exercise, Village maps were produced.<br /><br /><span style="color: rgb(0, 0, 0);">Legal and Institutional issues related to Village Land Demarcation.</span><br /><span style="color: rgb(51, 0, 0);">- Village Councils</span><br />- Village land adjudication committees<br />- Village Assembly<br />- District Council<br />- Ministry of Land and human Settlements Development<br /><span class="blsp-spelling-corrected" id="SPELLING_ERROR_8">Important</span> documents<br />- National Land Policy, 1995<br />- Village Land Act, 1999<br />- Village Land Act Regulations, 2002<br />- Village Land Act Guidelines developed by ministry of land.<br />Important trainings<br />- Stakeholders workshop - target (all stakeholders)<br />- Village land Act, 1999 - Target (district staff, ward, villages)<br />- Village boundary agreements procedures and processes -Target ( Village councils, Village Land Committees)<br />- Basic <span class="blsp-spelling-error" id="SPELLING_ERROR_9">GIS</span> training - Target (district council staff)<br />- Compass use and basic Map reading - Target (village game scouts [VGS]), village leadersUnknownnoreply@blogger.com0tag:blogger.com,1999:blog-4300720009318560726.post-60674694954899928942009-03-23T11:32:00.001+03:002009-03-23T11:32:51.973+03:00Helo world<div xmlns='http://www.w3.org/1999/xhtml'>Here I (and the rest of the MAE Team, I hope) will publish progresses, constraints, challenges on project implementation.<br/><br/>At least, I'll try to do it...<br/><br/>Rock<br/>and<br/>Roll<br/><br/><div class='zemanta-pixie'><img src='http://img.zemanta.com/pixy.gif?x-id=88ecfaae-1347-41b5-910f-e8ce8a985008' class='zemanta-pixie-img'/></div></div>Unknownnoreply@blogger.com0