geodata for guadiana, spain

April 3, 2009 at 3:29 pm · Filed under geodata, GRASS GIS, Linux, Uncategorized ·Tagged , , , ,

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