Lower Guadiana basin, Spain
This is my work area on geodata collection for Guadiana, Spain. In support of visiting scholar Irene Blanco
####Soils
###Sources
1. European Soils Database
Link:http://eusoils.jrc.ec.europa.eu/data.html
Resolution:
Format: vector; also raster files of 70+ properties in ESRI GRID format, Lambert Azimuthal coordinates;
Issues:
Status: pending;
####Rivers
###Sources
1. Hydro1K dataset for Europe
Link:http://edc.usgs.gov/products/elevation/gtopo30/hydro/eu_streams.html
Resolution: NA
Format: shapefile
Issues: This is based on flow accumulation analysis on GTOPO elevation data
Status: archive downloaded. unprocessed.
2. Digital Chart of the World
Link:http://www.maproom.psu.edu/dcw/
Resolution: NA
Format: arc info format (.eoo);Latlong Clarke1866
Issues: must be reprojected; dataset is old
Status: selected archive for Spain downloaded – drainage, also landcover, roads, populated places. unprocessed
Scripts:
#converting the drainage into shapefiles and also reprojecting to wgs84
ogr2ogr -f “ESRI Shapefile” -t_SRS EPSG:4326 DCW_drainage dnnet.e00
<ERROR 6: Can’t create fields of type IntegerList on shapefile layers.>
what is this?? try importing directly in Manifold/Arcview/ArcGIS
3.
http://dataservice.eea.europa.eu/dataservice/metadetails.asp?id=1053
4. European Commission Joint Research Center
Link:http://desert.jrc.ec.europa.eu/action/php/index.php?action=view&id=24
Resolution:NA
Format: Latlongwgs84; ESRI ArcGIS geodatabase; cannot use in GRASS, must switch to ArcGIS
Issues: archive downloaded for “SouthWest” tile which includes Spain. unprocessed; handover to Irene to continue in ArcGIS.
####DEM
###Sources
1. GTOPO30
Link:http://eros.usgs.gov/products/elevation/gtopo30/gtopo30.html
Resolution: 30 arc seconds (~1km)
Format: Binary, 16bit integer, Motorola MSB
Issues: Failed to import using several software; 2 tiles needed to cover lower Guadiana basin
Status: 2 tiles raw data downloaded, unprocessed
2. SRTM
Link:http://eros.usgs.gov/products/elevation/gtopo30/gtopo30.html
Resolution: 3 arc seconds (~90m)
Format: SRTM DEM format
Issues: Will need 20-30 tiles downloaded, imported and mosaiced
Status: 2 tiles raw data downloaded by Irene; unprocessed; abandoned for – no need for such high resolution for WEAP work
3. GLOBEDEM
Link:http://www.ngdc.noaa.gov/mgg/topo/
Resolution: 30 arc seconds (~1km)
Format: Binary, 16bit integer, Motorola MSB
Issues: Boundary coordinates fed to download rectangular region for most of Spain
Status: Successfully processed in GRASS; and also output to ESRI Ascii grid format which can be used in ArcGIS and Arcview.
Scripts:
###GRASS GIS
##March 31 Guadiana basin DEM work
#import GlobeDEM using r.in.bin (http://grass.itc.it/gdp/html_grass63/r.in.bin.html)
#r.in.gdal did not work
r.in.bin -sb input=Vishal_Guadiana.bin output=GlobeDEM north=44.0 south=35.9 east=2.0 west=-9.3 r=972 c=1356 bytes=2 anull=-500
# creating watershed at 1000km2 threshold
r.watershed GlobeDEM thresh=1000 accum=FlowAccum1000 drain=draindir1000 basin=basin_1000 stream=rivers1000
#creating and displaying shaded relief
r.shaded.relief GlobeDEM shadedmap=GlobeShaded
d.his h=elevation i=GlobeShaded
# clean up border bad watersheds
r.mapcalc basins=basin_1000
r.null map=basins setnull=98,100,102,104,106,108,110,112,114,116,118,120,662,660,798,1082,800,802,1078,804,936,944,1060,1068,1072,1100,1106,1110,1114,1118,1122,1126,1130,1166,1168,1170
r.null map=basins setnull=4,6,8,10,12,14,16,20,22,24,26,28,30,416,414,640,644,648,652,650,656,654
# convert basins to vector
r.to.vect -s basins out=basins feature=area
#display vector watersheds map
d.vect -c basins
#export to shapefile; -ep flags generate .prj and convert lines to polygons respectively.
v.out.ogr -ep input=basins type=area dsn=Watersheds.shp format=ESRI_Shapefile
v.info basins
# convert GlobeDEM into Esri Ascii grid; also the partially cleaned up ‘basins’ watersheds raster layer. the basins will need to be converted to vector and cleaned up later in Manifold/Arcview or ArcGIS
r.out.gdal format=AAIGrid input=GlobeDEM output=GlobeDEM.asc nodata=-9999
r.out.gdal format=AAIGrid input=basins output=basins.asc nodata=-9999
#
Leave a comment