"
]
@@ -2216,11 +2200,15 @@
"emodnet_open_maplibrary__Modiolus_habitatsuitabilitymodel_uk_50\n",
"emodnet_open_maplibrary__Sabellaria_habitatsuitabilitymodel_uk_50\n",
"emodnet_open_maplibrary__Zostera_habitatsuitabilitymodel_uk_50\n",
+ "emodnet_open_maplibrary__Hierarchical_class_Level1\n",
+ "emodnet_open_maplibrary__Hierarchical_class_Level2\n",
+ "emodnet_open_maplibrary__Hierarchical_class_Level3\n",
"emodnet_open_maplibrary__LAMpred_Nordic\n",
"emodnet_open_maplibrary__SACpred_Nordic\n",
"emodnet_open_maplibrary__mediseh_cora\n",
"emodnet_open_maplibrary__mediseh_maerl\n",
"emodnet_open_maplibrary__mediseh_posidonia\n",
+ "emodnet_open_maplibrary__Non_hierarchical_class\n",
"emodnet_open_maplibrary__DK004008_EFH_Cod_Spawning grounds\n",
"emodnet_open_maplibrary__DK004009_EFH_Cod_Feeding grounds\n",
"emodnet_open_maplibrary__DK004010_EFH_Cod_Nursery grounds\n",
@@ -2280,7 +2268,9 @@
"emodnet_open_maplibrary__PT005003_EFH_Sardine_SpawningGrounds\n",
"emodnet_open_maplibrary__PT005004_EFH_Sardine_SpawningGrounds\n",
"emodnet_open_maplibrary__PT005005_EFH_Sardine_SpawningGrounds\n",
- "emodnet_open_maplibrary__PT005006_EFH_Sardine_SpawningGrounds\n"
+ "emodnet_open_maplibrary__PT005006_EFH_Sardine_SpawningGrounds\n",
+ "emodnet_open_maplibrary__PT005007\n",
+ "emodnet_open_maplibrary__PT005008\n"
]
}
],
@@ -2464,8 +2454,9 @@
"outputs": [],
"source": [
"# get parameters for grid\n",
- "bbox2 = Study_area.bounds.values[0].tolist() # get bbox for the Western Mediterranean Sea study area \n",
- "crs2 = Study_area.crs.srs # get crs of the Study Area file"
+ "Study_area_3857 = Study_area.to_crs('EPSG:3857')\n",
+ "bbox2 = Study_area_3857.bounds.values[0].tolist() # get bbox for the Western Mediterranean Sea study area\n",
+ "crs2 = Study_area_3857.crs.srs"
]
},
{
@@ -2477,15 +2468,19 @@
"name": "stdout",
"output_type": "stream",
"text": [
- "[-7.6, 33.07066846, 18.22312212, 46.42720294]\n"
+ "[('X', -845389.1484423319, 2027949.6934754928), ('Y', 3951488.8573080758, 5753792.596194154)]\n"
]
}
],
"source": [
- "#Adding 2 decimal degrees to each side of the bbox\n",
+ "# Adding 2 decimal degrees to each side of the bbox\n",
"\n",
- "bbox2 = [bbox2[0]-2,bbox2[1]-2,bbox2[2]+2,bbox2[3]+2]\n",
- "print(bbox2)"
+ "buffer = 111_000 * 2\n",
+ "subsets = [\n",
+ " ('X', bbox2[0] - buffer, bbox2[2] + buffer),\n",
+ " ('Y', bbox2[1] - buffer, bbox2[3] + buffer)\n",
+ "]\n",
+ "print(subsets)"
]
},
{
@@ -2497,22 +2492,20 @@
"name": "stdout",
"output_type": "stream",
"text": [
- "https://ows.emodnet-seabedhabitats.eu/geoserver/emodnet_open_maplibrary/wcs?version=2.0.1&request=GetCoverage&service=WCS&CoverageID=emodnet_open_maplibrary__mediseh_cora&crs=EPSG%3A4326&format=image%2Ftiff\n"
+ "https://ows.emodnet-seabedhabitats.eu/geoserver/emodnet_open_maplibrary/wcs?version=2.0.1&request=GetCoverage&service=WCS&CoverageID=emodnet_open_maplibrary__mediseh_cora&format=image%2Ftiff&subset=X%28-845389.1484423319%2C2027949.6934754928%29&subset=Y%283951488.8573080758%2C5753792.596194154%29\n"
]
}
],
"source": [
"# GetCoverage request\n",
+ "\n",
"response = wcs.getCoverage(identifier=[dataset],\n",
- " bbox=bbox2,\n",
" format='image/tiff',\n",
- " crs=crs2,\n",
- " resx=0.00833333,resy=0.00833333,\n",
- " timeout = 120\n",
- " )\n",
+ " subsets=subsets,\n",
+ " timeout=120)\n",
"print(response.geturl())\n",
"with open(f'data/{dataset.replace(\":\",\"_\")}.tif','wb') as outfile:\n",
- " outfile.write(response.read())\n"
+ " outfile.write(response.read())"
]
},
{
@@ -2539,10 +2532,9 @@
"output_type": "stream",
"text": [
"Raster CRS: EPSG:3857\n",
- "Study Area CRS: EPSG:4326\n",
+ "Study Area CRS: EPSG:3857\n",
"Raster Bbox: (-44.0000000000998, 30.2670495513959, 68.5000001347984, 90.0000000001001)\n",
- "Study Area Bbox: [-7.6, 33.07066846, 18.22312212, 46.42720294]\n",
- "Need to reproject raster to match raster CRS.\n"
+ "Study Area Bbox: [-623389.1484423319, 4173488.8573080758, 1805949.6934754928, 5531792.596194154]\n"
]
}
],
@@ -2569,7 +2561,7 @@
},
{
"cell_type": "code",
- "execution_count": 39,
+ "execution_count": 35,
"metadata": {},
"outputs": [
{
@@ -2631,15 +2623,15 @@
},
{
"cell_type": "code",
- "execution_count": 40,
+ "execution_count": 36,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
- "{'driver': 'GTiff', 'dtype': 'float32', 'nodata': -3.4028234663852886e+38, 'width': 14294, 'height': 5155, 'count': 1, 'crs': CRS.from_epsg(4326), 'transform': Affine(0.0034645994105386834, 0.0, -6.557315834740004,\n",
- " 0.0, -0.0034645994105386834, 47.300772718499964)}\n"
+ "{'driver': 'GTiff', 'dtype': 'float32', 'nodata': -3.4028234663852886e+38, 'width': 7364, 'height': 3690, 'count': 1, 'crs': CRS.from_wkt('GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AXIS[\"Latitude\",NORTH],AXIS[\"Longitude\",EAST],AUTHORITY[\"EPSG\",\"4326\"]]'), 'transform': Affine(0.0033642872069247045, 0.0, -6.557315834740004,\n",
+ " 0.0, -0.0033642872069247045, 45.83349537089768)}\n"
]
}
],
@@ -2697,27 +2689,76 @@
" */\n",
"\n",
":root {\n",
- " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
- " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
- " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
- " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
- " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
- " --xr-background-color: var(--jp-layout-color0, white);\n",
- " --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
- " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
+ " --xr-font-color0: var(\n",
+ " --jp-content-font-color0,\n",
+ " var(--pst-color-text-base rgba(0, 0, 0, 1))\n",
+ " );\n",
+ " --xr-font-color2: var(\n",
+ " --jp-content-font-color2,\n",
+ " var(--pst-color-text-base, rgba(0, 0, 0, 0.54))\n",
+ " );\n",
+ " --xr-font-color3: var(\n",
+ " --jp-content-font-color3,\n",
+ " var(--pst-color-text-base, rgba(0, 0, 0, 0.38))\n",
+ " );\n",
+ " --xr-border-color: var(\n",
+ " --jp-border-color2,\n",
+ " hsl(from var(--pst-color-on-background, white) h s calc(l - 10))\n",
+ " );\n",
+ " --xr-disabled-color: var(\n",
+ " --jp-layout-color3,\n",
+ " hsl(from var(--pst-color-on-background, white) h s calc(l - 40))\n",
+ " );\n",
+ " --xr-background-color: var(\n",
+ " --jp-layout-color0,\n",
+ " var(--pst-color-on-background, white)\n",
+ " );\n",
+ " --xr-background-color-row-even: var(\n",
+ " --jp-layout-color1,\n",
+ " hsl(from var(--pst-color-on-background, white) h s calc(l - 5))\n",
+ " );\n",
+ " --xr-background-color-row-odd: var(\n",
+ " --jp-layout-color2,\n",
+ " hsl(from var(--pst-color-on-background, white) h s calc(l - 15))\n",
+ " );\n",
"}\n",
"\n",
- "html[theme=dark],\n",
- "body[data-theme=dark],\n",
+ "html[theme=\"dark\"],\n",
+ "html[data-theme=\"dark\"],\n",
+ "body[data-theme=\"dark\"],\n",
"body.vscode-dark {\n",
- " --xr-font-color0: rgba(255, 255, 255, 1);\n",
- " --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
- " --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
- " --xr-border-color: #1F1F1F;\n",
- " --xr-disabled-color: #515151;\n",
- " --xr-background-color: #111111;\n",
- " --xr-background-color-row-even: #111111;\n",
- " --xr-background-color-row-odd: #313131;\n",
+ " --xr-font-color0: var(\n",
+ " --jp-content-font-color0,\n",
+ " var(--pst-color-text-base, rgba(255, 255, 255, 1))\n",
+ " );\n",
+ " --xr-font-color2: var(\n",
+ " --jp-content-font-color2,\n",
+ " var(--pst-color-text-base, rgba(255, 255, 255, 0.54))\n",
+ " );\n",
+ " --xr-font-color3: var(\n",
+ " --jp-content-font-color3,\n",
+ " var(--pst-color-text-base, rgba(255, 255, 255, 0.38))\n",
+ " );\n",
+ " --xr-border-color: var(\n",
+ " --jp-border-color2,\n",
+ " hsl(from var(--pst-color-on-background, #111111) h s calc(l + 10))\n",
+ " );\n",
+ " --xr-disabled-color: var(\n",
+ " --jp-layout-color3,\n",
+ " hsl(from var(--pst-color-on-background, #111111) h s calc(l + 40))\n",
+ " );\n",
+ " --xr-background-color: var(\n",
+ " --jp-layout-color0,\n",
+ " var(--pst-color-on-background, #111111)\n",
+ " );\n",
+ " --xr-background-color-row-even: var(\n",
+ " --jp-layout-color1,\n",
+ " hsl(from var(--pst-color-on-background, #111111) h s calc(l + 5))\n",
+ " );\n",
+ " --xr-background-color-row-odd: var(\n",
+ " --jp-layout-color2,\n",
+ " hsl(from var(--pst-color-on-background, #111111) h s calc(l + 15))\n",
+ " );\n",
"}\n",
"\n",
".xr-wrap {\n",
@@ -2758,7 +2799,7 @@
".xr-sections {\n",
" padding-left: 0 !important;\n",
" display: grid;\n",
- " grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
+ " grid-template-columns: 150px auto auto 1fr 0 20px 0 20px;\n",
"}\n",
"\n",
".xr-section-item {\n",
@@ -2766,11 +2807,14 @@
"}\n",
"\n",
".xr-section-item input {\n",
- " display: none;\n",
+ " display: inline-block;\n",
+ " opacity: 0;\n",
+ " height: 0;\n",
"}\n",
"\n",
".xr-section-item input + label {\n",
" color: var(--xr-disabled-color);\n",
+ " border: 2px solid transparent !important;\n",
"}\n",
"\n",
".xr-section-item input:enabled + label {\n",
@@ -2778,6 +2822,10 @@
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
+ ".xr-section-item input:focus + label {\n",
+ " border: 2px solid var(--xr-font-color0) !important;\n",
+ "}\n",
+ "\n",
".xr-section-item input:enabled + label:hover {\n",
" color: var(--xr-font-color0);\n",
"}\n",
@@ -2799,7 +2847,7 @@
"\n",
".xr-section-summary-in + label:before {\n",
" display: inline-block;\n",
- " content: 'â–º';\n",
+ " content: \"â–º\";\n",
" font-size: 11px;\n",
" width: 15px;\n",
" text-align: center;\n",
@@ -2810,7 +2858,7 @@
"}\n",
"\n",
".xr-section-summary-in:checked + label:before {\n",
- " content: 'â–¼';\n",
+ " content: \"â–¼\";\n",
"}\n",
"\n",
".xr-section-summary-in:checked + label > span {\n",
@@ -2882,15 +2930,15 @@
"}\n",
"\n",
".xr-dim-list:before {\n",
- " content: '(';\n",
+ " content: \"(\";\n",
"}\n",
"\n",
".xr-dim-list:after {\n",
- " content: ')';\n",
+ " content: \")\";\n",
"}\n",
"\n",
".xr-dim-list li:not(:last-child):after {\n",
- " content: ',';\n",
+ " content: \",\";\n",
" padding-right: 5px;\n",
"}\n",
"\n",
@@ -2907,7 +2955,9 @@
".xr-var-item label,\n",
".xr-var-item > .xr-var-name span {\n",
" background-color: var(--xr-background-color-row-even);\n",
+ " border-color: var(--xr-background-color-row-odd);\n",
" margin-bottom: 0;\n",
+ " padding-top: 2px;\n",
"}\n",
"\n",
".xr-var-item > .xr-var-name:hover span {\n",
@@ -2918,6 +2968,7 @@
".xr-var-list > li:nth-child(odd) > label,\n",
".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
" background-color: var(--xr-background-color-row-odd);\n",
+ " border-color: var(--xr-background-color-row-even);\n",
"}\n",
"\n",
".xr-var-name {\n",
@@ -2967,8 +3018,15 @@
".xr-var-data,\n",
".xr-index-data {\n",
" display: none;\n",
- " background-color: var(--xr-background-color) !important;\n",
- " padding-bottom: 5px !important;\n",
+ " border-top: 2px dotted var(--xr-background-color);\n",
+ " padding-bottom: 20px !important;\n",
+ " padding-top: 10px !important;\n",
+ "}\n",
+ "\n",
+ ".xr-var-attrs-in + label,\n",
+ ".xr-var-data-in + label,\n",
+ ".xr-index-data-in + label {\n",
+ " padding: 0 1px;\n",
"}\n",
"\n",
".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
@@ -2981,6 +3039,12 @@
" float: right;\n",
"}\n",
"\n",
+ ".xr-var-data > pre,\n",
+ ".xr-index-data > pre,\n",
+ ".xr-var-data > table > tbody > tr {\n",
+ " background-color: transparent !important;\n",
+ "}\n",
+ "\n",
".xr-var-name span,\n",
".xr-var-data,\n",
".xr-index-name div,\n",
@@ -3040,45 +3104,55 @@
" stroke: currentColor;\n",
" fill: currentColor;\n",
"}\n",
- "<xarray.DataArray (band: 1, y: 2702, x: 6300)>\n",
- "[17022600 values with dtype=float32]\n",
+ "\n",
+ ".xr-var-attrs-in:checked + label > .xr-icon-file-text2,\n",
+ ".xr-var-data-in:checked + label > .xr-icon-database,\n",
+ ".xr-index-data-in:checked + label > .xr-icon-database {\n",
+ " color: var(--xr-font-color0);\n",
+ " filter: drop-shadow(1px 1px 5px var(--xr-font-color2));\n",
+ " stroke-width: 0.8px;\n",
+ "}\n",
+ "<xarray.DataArray (band: 1, y: 2782, x: 6488)> Size: 72MB\n",
+ "[18049616 values with dtype=float32]\n",
"Coordinates:\n",
- " * band (band) int32 1\n",
- " * x (x) float64 -5.599 -5.596 -5.592 -5.589 ... 16.22 16.22 16.22\n",
- " * y (y) float64 44.43 44.42 44.42 44.42 ... 35.08 35.08 35.07 35.07\n",
- " spatial_ref int32 0\n",
+ " * band (band) int64 8B 1\n",
+ " * x (x) float64 52kB -5.6 -5.597 -5.593 -5.59 ... 16.22 16.22 16.22\n",
+ " * y (y) float64 22kB 44.43 44.42 44.42 44.42 ... 35.08 35.07 35.07\n",
+ " spatial_ref int64 8B 0\n",
"Attributes:\n",
" AREA_OR_POINT: Area\n",
" _FillValue: -3.4028235e+38\n",
" scale_factor: 1.0\n",
- " add_offset: 0.0PandasIndex
PandasIndex(Index([1], dtype='int32', name='band'))
PandasIndex
PandasIndex(Index([ -5.599354097726058, -5.595889498315519, -5.592424898904981,\n",
- " -5.588960299494442, -5.585495700083904, -5.5820311006733645,\n",
- " -5.578566501262826, -5.575101901852287, -5.571637302441749,\n",
- " -5.56817270303121,\n",
+ " add_offset: 0.0
PandasIndex
PandasIndex(Index([1], dtype='int64', name='band'))
PandasIndex
PandasIndex(Index([ -5.600176124369926, -5.596811837163002, -5.5934475499560765,\n",
+ " -5.590083262749152, -5.586718975542228, -5.5833546883353025,\n",
+ " -5.579990401128378, -5.576626113921454, -5.573261826714528,\n",
+ " -5.569897539507604,\n",
" ...\n",
- " 16.192976194562263, 16.196440793972798, 16.19990539338334,\n",
- " 16.203369992793874, 16.206834592204416, 16.21029919161495,\n",
- " 16.213763791025492, 16.217228390436034, 16.22069298984657,\n",
- " 16.22415758925711],\n",
- " dtype='float64', name='x', length=6300))PandasIndex
PandasIndex(Index([44.426887507458126, 44.423422908047584, 44.41995830863705,\n",
- " 44.41649370922651, 44.41302910981597, 44.40956451040543,\n",
- " 44.4060999109949, 44.402635311584355, 44.39917071217381,\n",
- " 44.39570611276328,\n",
+ " 16.193676402088307, 16.197040689295232, 16.200404976502156,\n",
+ " 16.203769263709084, 16.20713355091601, 16.210497838122933,\n",
+ " 16.213862125329857, 16.21722641253678, 16.220590699743706,\n",
+ " 16.22395498695063],\n",
+ " dtype='float64', name='x', length=6488))
PandasIndex
PandasIndex(Index([ 44.42554117479969, 44.42217688759276, 44.41881260038584,\n",
+ " 44.41544831317891, 44.41208402597199, 44.408719738765065,\n",
+ " 44.405355451558144, 44.401991164351216, 44.398626877144295,\n",
+ " 44.39526258993737,\n",
" ...\n",
- " 35.10018589428799, 35.096721294877455, 35.09325669546691,\n",
- " 35.08979209605637, 35.086327496645836, 35.082862897235294,\n",
- " 35.07939829782476, 35.07593369841422, 35.07246909900368,\n",
- " 35.06900449959314],\n",
- " dtype='float64', name='y', length=2702))
- AREA_OR_POINT :
- Area
- _FillValue :
- -3.4028235e+38
- scale_factor :
- 1.0
- add_offset :
- 0.0
- AREA_OR_POINT :
- Area
- _FillValue :
- -3.4028235e+38
- scale_factor :
- 1.0
- add_offset :
- 0.0