storing, retrieving and analysing marine ecosystem data of space. and Jan Erik Stiansen

Transcription

1 A framework for storing, retrieving and analysing marine ecosystem data of different origin with variable scale and distribution in time and space. Trond Westgård Geir Odd Johansen Cecilie Kvamme Bjørn Ådlandsvik Trond Westgård, Geir Odd Johansen, Cecilie Kvamme, Bjørn Ådlandsvik and Jan Erik Stiansen

2 Background Concepts Technical solution Example analysis Future direction

3 Background The project FishExchange h investigates t the possible effects of climate change on the spatial distributions ib ti of commercial fish stocks, fishery and related socio-economical consequences in the Barents Sea. A software system and a framework was needed to facilitate such analysis.

4 Background Concepts Technical solution Example analysis Future direction

5 Main purpose Force parameters of different origin into the same GRIDDED STORAGE horizontally (latitude, longitude), vertically and in time in order to facilitate analysis of dependencies between parameters having heterogeneous resolution in time and space.

6 Gridding in time and space Horizontally Each cell/strata is a WGS84 multipolygon. The cells may not overlap. Vertically Each layer may be free, whole watercolumn, Bottom layer with thickness or a fixed layer Time A time step may be free, day, month, quarter or year.

7 Principles A parameter may have ONE value in a cell for a given set of time steps and depth layers. A cell / strata t belongs to ONE grid / strata t system. A cell is a MULTIPOLYGON. A parameter may be used in several grids It s the responsibility of the owner of the parameter to interpolate parameter values into this gridded strait jacket. Owners of parameters who wish to analyze dependencies between their parameters must agree upon a suitable common resolution in time and space.

8 Background Concepts Technical solution Example analysis Future direction

9 Data model

10 Tools used PostgreSQL with PostGIS NetBeans 6.1 with Java and Swing R-statistical software Python with numpy ArcGIS

11 Functionality developed d User and user roles Define parameters with metadata (EdMed) Define equalarea grids (Lambert azimutal) Import/export format for parameters(.csv file) Export of grid to shapefile Calculation of area-weighted mean values

12 Users logs on with name/role

13 EDMED parameter definition

14 Definition of equal area grid

15 Import of data

16 Export of data

17 Import/export format GridName,Lat,Lon,GridCellName,Value,ParameterName,TimeFrom,TimeTo,TimePeriods,Depthfrom,DepthTo,DepthLayers FishExchange,73.5,35.8,,6.9,temperature, ,,,,, 00:00, :59,(F;M200101),0,50,(W;L0:50), ; ) FishExchange,73.5,35.8,,6.3,temperature, :00, ,(F;M200102),0,50,(W;L0:50)

18 Export of equal area grid

19 FishExchange grid (25x25) km with GEBCO depths

20 ODBC and R to get cell id library (RODBC) channel<-odbcconnect("grid", uid="postgres", pwd="******", case="postgresql") sql<-paste("select name as cellname from hcell where within(geomfromewkt('srid=4326;point(",position,")'),geoshape) (S ( )) ) and id_grid = 9 and (maxlat < ",maxlat,") and (minlat > ",minlat, ") and (maxlon < ",maxlon,") and (minlon > ",minlon,");",sep="") sqlut<- sqlquery(channel, sql) testfil$gridcellname[i]<-as.character(sqlut$cellname)

21 Background Concepts Technical solution Example analysis Future direction

22 Data into gridded format Data into gridded format Conmap (spline/laplace) = Temp atlas

23 Data into gridded format mean value per cell R.vessel

24 Data into gridded format proportion of area - fisheries

25 55 0'0"E Eastward displacement in warm years 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 76 0'0"N 1987 Q1 B50 temperature 1992 Q1 B50 temperature 76 0'0"N 67 0'0"N 67 0'0"N 0 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 76 0'0"N 76 0'0"N 1987 Q1 cod >45 cm bt rv survey 1992 Q1 cod >45 cm bt rv survey [ [ 67 0'0"N 67 0'0"N

26 Fisheries pattern do not directly correspond with temperature 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 76 0'0"N 1987 Q1 B50 temperature 1992 Q1 B50 temperature 76 0'0"N 67 0'0"N 67 0'0"N 0 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 76 0'0"N 1987 Q1 cod bt fisheries 1992 Q1 cod bt fisheries 76 0'0"N 0 0 [ [ 67 0'0"N 67 0'0"N

27 Temp. difference at bottom 50 meters (1992 minus 1987) 5 0'0"E 10 0'0"E 50 0'0"E 55 0'0"E 76 0'0"N < '0"N

28 Background Concepts Technical solution Example analysis Future direction

29 Open source project??? OpenSwing, ibatis, PostgreSQL/PostGIS,Geotools, R